In order to optimize the tip of the wing turbine blade, a two-stage approach was developed. At first, the external aerodynamic characteristics of the tip were modified by employing a shape optimization algorithm aiming to maximize the generated torque without major changes on the overall length of the blade and the implemented blended curved tip [
7]. Subsequently, the internal structure of the blade underwent computational topology optimization study towards an improved geometry aiming for maximum compliance at least material condition [
13]. During this step, the design flexibility offered by AM was taken into consideration, presenting a significant contribution towards the manufacturability of the optimized component [
14]. The overall roadmap of the aforementioned approach is presented in
Figure 1.
2.1. Shape Optimization
Shape design optimization is a vital part of the product design workflow, which approaches the design process in a systemic manner, investigating the effect of all input and output parameters of the given system and adjusting each accordingly with the scope of producing the best possible solution for a given design problem [
17,
18].
The adjoint solver methodology is exploited as a shape optimization technique, which makes use of a specialized mathematical tool that provides detailed sensitivity analysis regarding the performance of a fluid system with solid elements that is subjected to specific boundary conditions [
19]. The adjoint optimization method considers the computational flow solution
and defines the problem inputs in the form of a vector
. The quantity of interest is
, and the residuals of the Navier Stokes equations are
. By defining the Lagrangian L with the vector of Lagrange multipliers
(the adjoint solution variable), the following equation takes form:
By choosing a
such that:
the problem is reduced to linear form:
where the left-hand side indicates the adjoint sensitivities. The sensitivity equation is evaluated at each given mesh node in the CFD model. For shape sensitivity, the input vector
is considered as the (x, y, z) coordinates of every node in the model. In the same equation, the
factor is the change in J value due to changes in the (x, y, z) coordinates, and
refers to the changes induced in J due to the sensitivity of the flow solution at given nodes locations, depending on the adjoint solution itself. Both parameters are calculated using expressions deriving from the observable definitions and the discretized computational model equations. Thus, by solving the adjoint problem:
the computation of the shape derivative takes place in even large-scale optimization cycles. This solution is carried out as follows. The adjoint residuals are calculated using the formula:
The correction to the solution is determined:
The adjoint equations are updated using the corrected vector:
This iterative process repeats until the specified threshold is reached or the predefined maximum number of iterations is complete. The aim of this work is to extend the scope of a CFD study by extracting data with regard to the system’s sensitivity against geometrical changes by computing the derivative of engineering quantities, taking into consideration every system input, including flow geometry [
20]. This capability is utilized to realize intelligent design modifications for shape optimization of any geometric feature in the defined computational domain and further identify the optimal shape for given operating conditions, deriving from a baseline CFD flow calculation. As such, no model parameterization is necessary, making it a suitable candidate for intricate non-parametric geometries [
21]. At first, the study makes use of user-defined relevant observables. Various scalar quantities of interest can act as observables, indicatively forces exerted on the component under investigation, flow parameters and geometric features. Either a continuous or discrete solver is selected, depending on the case, with the first option focusing on speed and memory efficiency and the second achieving higher accuracy by increasing the computational cost. Subsequently, accordingly tuned controls for the optimization procedure are utilized, aiming to produce a stable and robust calculation. Predetermined residuals with suitable convergence values balance precision and computational cost. After the completion of the aforementioned steps, the solver performs the calculation in order to identify the system sensitivities, essentially predicting how geometry alterations would affect the performance of the part. Optimization cycles are suggested to be carried out in short cycles, with incremental sets of moderate or small improvement each time, up to 10% of the initial value, as larger increments pose risk of causing solution overshoot and potentially guiding the solver to converge on suboptimal results. The described workflow is presented schematically in
Figure 2. Limiting the optimization area locally to the features of interest, such as the curved tip of the wind turbine blade in this specific scenario (
Figure 3), can further aid in minimizing the computational cost and isolating specific areas that require improvement while preserving others. A stabilization scheme was implemented to support the concurrence of the study within the preset number of iterations by monitoring local flow rate and air velocity residual value convergence around the tip. The solver was indicatively set to a maximum of 1000 calculation iterations, and the residual convergence was set at 10
−5 for increased accuracy over the default values of 10
−3.
After the completion of an appropriate number of calculation rounds, the user is able to carry out geometry adjustments, based on the calculation results, and extract the updated, potentially optimized geometry. This geometry has to undergo further simulations to validate the actual improvement.
The simulation scenario that was carried out as part of this optimization study to produce initial and post optimization performance results assumed constant rated wind velocity and rotor speed, set at 8.56 m/s and 6.22 rpm, respectively, in an indicative scenario aiming to support the two stage optimization workflow. A single blade of a three-blade rotor with a 240 m total diameter was examined, considering the symmetry of the rotor’s layout, with a −4° precone angle. The study utilized a finite conical mesh of tetrahedral elements (
Figure 4), with a maximum element size of 24 m at the domain that adapted to smaller values near the Body of Interest (BOI), maintaining a steady size of 1.8 m around the blade and reaching down to 0.008 m near intricate features, such as the trailing edge of the blade. The element size values were selected through a mesh independence study, undergoing 4 consecutive refinement steps.
The CFD study furthermore utilized a Shear Stress Transfer (SST) k-ω turbulence model to capture turbulent flows due to expected Reynolds numbers up to 2 × 10
6, considering the geometric and operating characteristics of the wind turbine under investigation. The studies were carried out in Ansys Fluent 2021 R1, employing time-dependent Reynolds Averaged Navier Stokes (RANS) equations [
22]. The accuracy and reliability of the model was verified through comparison of the values for the tip speed and power coefficient. Whereas the tip speed remained within 1% of the nominal design value during every mesh refinement step, the difference between rated and calculated power coefficient decreased from 30.4% to 19.6% after the 4th refinement step. Further refinement of the mesh was not made possible at this point due to computational limitations.
2.2. Topology Optimization
Topology optimization (TO) is employed to achieve material reduction, essential for structural applications, where weight reduction is of the essence, and the quality and robustness of the structure may also require enhancement or preservation to a certain degree. It takes into consideration the material and operating conditions of the designed component in order to produce the most efficient geometry that meets its requirements, also meeting manufacturing constraints [
23]. TO increases the design flexibility, as well as shortens the design process, and presents the potential to decrease lead times in product development, although it is a generally time intensive, iterative computational procedure [
24].
One of the most widely employed algorithms to perform topology optimization is the Solid Isotropic Material with Penalization method (SIMP) [
25]. Initially introduced by Bendsoe and Kikuchi in 1988, this methodology predicts the best available material distribution within a given design space, for predefined load cases, boundary conditions, manufacturing constraints and performance requirements and updated versions still widely employed in various optimization applications [
26,
27]. According to SIMP, the domain is discretized as a grid of finite elements named isotropic microstructures. These elements can be assigned a binary density value, either containing material (state 1) or being void (state 0). SIMP introduces a continuous relative density distribution varying between 1 and a specified minimum value ρ
min. A penalty factor p decreases the contribution of elements with intermediate densities to the total stiffness of the component. Using an iterative process, the optimization algorithm performs a sensitivity analysis to evaluate the impact of the variable material densities on the objective function to meet the objective (e.g., maximization of part stiffness) [
28]. Mathematically, this sensitivity analysis is expressed as a derivative of the objective function with respect to the material densities:
where ρ
ε indicates the density vector value of the elements; K
e describes the stiffness matrix for the elements; and u
e is the nodal displacement vectors calculated during the iterative procedure [
29]. The tools utilized within the context of this study employ the SIMP combined with elements of an alternative technique named a level set method. Level set methods define the three-dimensional design space as an open bounded
. The final, optimized shape is described by an additional set Ω, which is a subset of D. A level set function is defined such that:
The evolution of the shape Ω depends on the Hamilton Jacobi equation:
where v is called the descent direction and is calculated by the optimizer in an attempt to morph the initial geometry towards a predefined goal [
30,
31].
The overall workflow for a structural TO study resembles the one described in the previous paragraph up to a degree and is depicted in
Figure 5. An initial geometry is required, as well as information regarding material properties. Boundary conditions are specified, such as part mechanical or thermal loading and fixtures. Within the context of this specific approach, loading data is derived from the first stage, where the CFD analysis resulted in comprehensive pressure distribution profiles along the surface of the component (
Figure 6). A design space in which shape changes can occur is also determined, allowing in a similar manner to external shape optimization to isolate regions of interest and potentially exclude sensitive areas, such as features with limited thickness, where material removal would lead to non-manufacturable parts or preserve vital segments of the initial geometry, such as bolt orifices or threads, necessary for assembly. A benchmarking Finite Element Analysis (FEA) is carried out to determine the behavior of the initial geometry. Following this step, an optimization goal must be user-defined, such as a mass reduction percentage or a compliance minimization target. The optimization algorithm follows this study, calculating the system sensitivities against shape alteration and material removal. A filtering software tool is then utilized to remove potential gray areas produced during this step, aiming for a clean solution geometry. An iterative optimization procedure is subsequently performed, incrementally calculating the shape alteration towards the predefined goals up to the point of convergence between the set and targeted values or until no further significant changes occur, signifying saturation. A validation FEA study confirms the effectiveness of the optimization, and the user receives the updated, enhanced geometry.
In order to determine the limitation point of the study, two distinct response constraints were set with regard to the amount of material to be removed and the strength of the resulting part. A 40% mass reduction limit was placed, and a safety factor of value 2 concerning the equivalent Von Mises stresses was developed on the component. The maximum number of optimization iterations was set at 20. This value is typically acceptable for a topology optimization analysis, with many instances usually converging in fewer. If no convergence is detected in the minimization of the compliance or the constraints set around mass removal and the safety factor, the process can continue for a larger user-defined number of iterations. A program-controlled solver type was selected as no specific case constraints dictated for the application of an alternative method. The convergence accuracy was retained at the default value of 0.1%. Aiming to reduce computational time and avoid implications caused by implementing thin and small features in the study, the bladelet geometry was compartmentalized into 4 distinct sub-areas: the joining base; the main body; the trailing edge; and the thin tip of the bladelet (
Figure 7). The region in which the solver will perform alterations and the areas to be excluded from the study are defined. In this case the hollowed-out base of the bladelet (
Figure 7b) is excluded from the study, as much material has already been removed to facilitate attachment of the tip to the blade during assembly. Moreover, the trailing edge of the bladelet and the tip were also excluded due to their thinness, which would potentially increase the computational cost and produce parts with non-manufacturable thin features (
Figure 7c,d). The external surface of the selected region for the optimization was also preserved so as to avoid removing material from the outer shell of the bladelet. The objective of the optimization was set in the form of minimization of the compliance of the given design throughout the material removal process.
2.3. Validation through 3D Printing
Whereas structured industrial design usually consists of well-defined parametric features, shapes produced via shape and topology optimization are mostly comprised of free-form constituents, which are often challenging to manufacture using traditional, conventional fabrication techniques, such as CNC milling or composite fabric layup. Pairing optimization results with additive manufacturing technologies due to their ability to directly realize these complex, usually organic shapes can potentially provide the adopted methodology with a significantly higher degree of flexibility in comparison to conventional production workflows. In this current study, the development of a prototype for further experimentation and physical validation of the modeling process dictated the manufacturing of a scaled-down version of the optimized bladelet through 3D printing. The aforementioned investigation is not included in this publication. Vat Photopolymerization (VPP) 3D printing utilizes a light source, a vertically translating build plate, which is immersed in a transparent basin with liquid photopolymer resin. The light source is projected through an LCD screen, which acts as a masking layer, allowing for selective emission over the photopolymer, procedurally solidifying the resin in a layer-by-layer manner, leading to a 3D object (
Figure 8). It is a 3D printing technology capable of attaining complex parts with intricate internal features and high surface quality. This technology further boasts a good degree of compatibility in terms of commercially available materials and high fabrication speed, compared to other AM technologies and compatibility with several translucent resins in order to fabricate partially transparent components, aiming to support the visualization of the internal configuration that occurred as a result of the TO process.
The demonstration of the methodology results presented in this study is carried out using a LC Magna VPP 3D printing unit, whose large build volume of 510 × 270 × 350 mm designates it as a suitable match for the fabrication of scaled-down prototypes for wind turbine components. The material selected for printing was the Daylight Magna Draft Resin by Photocentric, due to its notable printing quality for rapid prototyping. This specific material displays notable mechanical properties, presenting high strength and stiffness, comparable to rigid engineering thermoplastics, such as Polyamide 12. 3D printing of this resin also results in translucent parts where any internal features are easily visible.