Search results

The purpose of this paper is the presentation of a technique to be integrated in a numerical airfoil optimization scheme, for the exact satisfaction of a strict equality cross-sectional area constraint.

An airfoil optimization framework is presented, based on Area-Preserving Free-Form Deformation (AP FFD) technique. A parallel metamodel-assisted differential evolution (DE) algorithm is used as an optimizer. In each generation of the DE algorithm, before the evaluation of the fitness function, AP FFD is applied to each candidate solution, via coupling a classic B-Spline-based FFD with an area correction step. The area correction step is achieved by solving a sub problem, which consists of computing and applying the minimum possible offset to each one of the free-to-move control points of the FFD lattice, subject to the area preservation constraint.

The proposed methodology is able to obtain better values of the objective function, compared to both a classic penalty function approach and a generic framework for handling constraints, which suggests the separation of constraints and objectives (separation-sub-swarm), without any loss of the convergence capabilities of the DE algorithm, while it also guarantees an exact area preservation. Due to the linearity of the area constraint in each axis, the extraction of an inexpensive closed-form solution to the sub problem is possible by using the method of Lagrange multipliers.

AP FFD can be easily incorporated into any 2D shape optimization/design process, as it is a time-saving and easy-to-implement repair algorithm, independent from the nature of the problem at hand.

The proposed methodology proved to be an efficient tool in facing airfoil design problems, enhancing the rigidity of the optimal airfoil by preserving its cross-sectional area to a predefined value.

This study aims to feature the application of the artificial compressibility method (ACM) for the numerical prediction of two-dimensional (2D) axisymmetric swirling flows.

The respective academic numerical solver, named IGal2D, is based on the axisymmetric Reynolds-averaged Navier–Stokes (RANS) equations, arranged in a pseudo-Cartesian form, enhanced by the addition of the circumferential momentum equation. Discretization of spatial derivative terms within the governing equations is performed via unstructured 2D grid layouts, with a node-centered finite-volume scheme. For the evaluation of inviscid fluxes, the upwind Roe’s approximate Riemann solver is applied, coupled with a higher-order accurate spatial reconstruction, whereas an element-based approach is used for the calculation of gradients required for the viscous ones. Time integration is succeeded through a second-order accurate four-stage Runge-Kutta method, adopting additionally a local time-stepping technique. Further acceleration, in terms of computational time, is achieved by using an agglomeration multigrid scheme, incorporating the full approximation scheme in a V-cycle process, within an efficient edge-based data structure.

A detailed validation of the proposed numerical methodology is performed by encountering both inviscid and viscous (laminar and turbulent) swirling flows with axial symmetry. IGal2D is compared against the commercial software ANSYS fluent – by using appropriate metrics and characteristic flow quantities – but also against experimental measurements, confirming the proposed methodology’s potential to predict such flows in terms of accuracy.

This study provides a robust methodology for the accurate prediction of swirling flows by combining the axisymmetric RANS equations with ACM. In addition, a detailed description of the convective flux Jacobian is provided, filling a respective gap in research literature.

The purpose of this paper is to design an S‐duct intake for unmanned aerial vehicles (UAVs) applications with good efficiency in a wide range of operating conditions.

A fully‐parametric 3‐D CAD model of the intake was constructed in order to produce different intake configurations, within specific geometric constraints, and to study the influence of geometry variation on efficiency. O‐type blocking methodology was adopted in order to construct the block‐structured mesh of hexahedral elements, used in the simulations. The commercial CFD code ANSYS‐CFX was used to compute the flow field inside the flow domain of each case considered. The Reynolds averaged Navier‐Stokes (RANS) equations are discretized using an implicit, vertex‐based finite volume method, combined with the shear stress transport (SST) two‐equation turbulence model and an automatic wall treatment.

By shortening the axial length the flow separation after the first turning becomes more pronounced and the losses are increasing. For very long ducts the increased internal wall area leads to increased wall friction and, consequently, to increased loss production.

The adoption of Gerlach‐shaped profiles for the design of the S‐duct resulted in a low pressure loss level for the optimal shape, although more uniform distribution of total pressure losses resulted for ducts longer than the optimal one, which should be taken into account in the design process.