Search results

The purpose of this paper is to use dynamic meshing to perform CFD analyses of a NACA 0012 airfoil fitted with a morphing trailing edge (TE) flap when it undergoes static and time-dependent morphing. The steady CFD predictions of the original and morphing airfoils are validated against published data. The study also investigates an airfoil with a hinged TE flap for aerodynamic performance comparison. The study further extends to an unsteady CFD analysis of a dynamically morphing TE flap for proof-of-concept and also to realise its potential for future applications.

An existing parametrization method was modified and implemented in a user-defined function (UDF) to perform dynamic meshing which is essential for morphing airfoil unsteady simulations. The results from the deformed mesh were verified to ensure the validity of the adopted mesh deformation method. ANSYS Fluent software was used to perform steady and unsteady analysis and the results were compared with computational predictions.

Steady computational results are in good agreement with those from OpenFOAM for a non-morphing airfoil and for a morphed airfoil with a maximum TE deflection equal to 5 per cent of the chord. The results obtained by ANSYS Fluent show that an average of 6.5 per cent increase in lift-to-drag ratio is achieved, compared with a hinged flap airfoil with the same TE deflection. By using dynamic meshing, unsteady transient simulations reveal that the local flow field is influenced by the morphing motion.

An airfoil parametrisation method was modified to introduce time-dependent morphing and used to drive dynamic meshing through an in-house-developed UDF. The morphed airfoil’s superior aerodynamic performance was demonstrated in comparison with traditional hinged TE flap. A methodology was developed to perform unsteady transient analysis of a morphing airfoil at high angles of attack beyond stall and to compare with published data. Unsteady predictions have shown signs of rich flow features, paving the way for further research into the effects of a dynamic flap on the flow physics.

The paper aims to develop an efficient data-driven modeling approach for the hydroelastic analysis of a semi-circular pipe conveying fluid with elastic end supports. Besides the structural displacement-dependent unsteady fluid force, the steady one related to structural initial configuration and the variable structural parameters (i.e. the variable support stiffness) are considered in the modeling.

The steady fluid force is treated as a pipe preload, and the elastically supported pipe-fluid model is dealt with as a prestressed hydroelastic system with variable parameters. To avoid repeated numerical simulations caused by parameter variation, structural and hydrodynamic reduced-order models (ROMs) instead of conventional computational structural dynamics (CSD) and computational fluid dynamics (CFD) solvers are utilized to produce data for the update of the structural, hydrodynamic and hydroelastic state-space equations. Radial basis function neural network (RBFNN), autoregressive with exogenous input (ARX) model as well as proper orthogonal decomposition (POD) algorithm are applied to modeling these two ROMs, and a hybrid framework is proposed to incorporate them.

The proposed approach is validated by comparing its predictions with theoretical solutions. When the steady fluid force is absent, the predictions agree well with the “inextensible theory”. The pipe always loses its stability via out-of-plane divergence first, regardless of the support stiffness. However, when steady fluid force is considered, the pipe remains stable throughout as flow speed increases, consistent with the “extensible theory”. These results not only verify the accuracy of the present modeling method but also indicate that the steady fluid force, rather than the extensibility of the pipe, is the leading factor for the differences between the in- and extensible theories.

The steady fluid force and the variable structural parameters are considered in the data-driven modeling of a hydroelastic system. Since there are no special restrictions on structural configuration, steady flow pattern and variable structural parameters, the proposed approach has strong portability and great potential application for other hydroelastic problems.

The purpose of this paper is to investigate and improve a new method of unstructured rotational dynamic overset grids, which can be used to simulate the unsteady flows around rotational parts of aircraft.

The computational domain is decomposed into two sub-domains, namely, the rotational sub-domain which contains the rotational boundaries, and the stationary sub-domain which contains the remainder flow field including the stationary boundaries. The artificial boundaries and restriction boundaries are used as the restriction condition to generate the entire computational grid, and then the overset grids are established according to the radius parameters of artificial boundaries set previously. The deformation of rotational boundary is treated by using the linear spring analogy method which is suitable for the dynamic unstructured grid. The unsteady Navier-Stokes/Euler equations are solved separately in the rotational sub-domain and stationary sub-domain, and data coupling is accomplished through the overlapping area. The least squares method is used to interpolate the flow variables for the artificial boundary points with a higher calculating precision. Implicit lower-upper symmetric-Gauss-Seidel (LU-SGS) time stepping scheme is implemented to accelerate the inner iteration during the unsteady simulation.

The airfoil steady flow, airfoil pitching unsteady flow, three-dimensional (3-D) rotor flow field, rotor-fuselage interaction unsteady flow field and the flutter exciting system unsteady flow field are numerically simulated, and the results have good agreements with the experimental data. It is shown that the present method is valid and efficient for the prediction of complicated unsteady problems which contain rotational dynamic boundaries.

The results are entirely based on computational fluid dynamics (CFD), and the 3D simulations are based on the Euler equations in which the viscous effect is ignored. The current work shows further applicable potential to simulate unsteady flow around rotational parts of aircraft.

The current study can be used to simulate the two-dimensional airfoil pitching, 3-D rotor flow field, rotor-fuselage interaction and the flutter exciting system unsteady flow. The work will help the aircraft designer to get the unsteady flow character around rotational parts of aircraft.

A new type of rotational dynamic overset grids is presented and validated, and the current work has a significant contribution to the development of unstructured rotational dynamic overset grids.

This paper aims to describe a research effort towards the comprehension of the unsteady phenomena due to the deployment of high-lift devices at approach/landing conditions.

The work starts from a preexisting platform based on an immersed boundary (IB) method whose capabilities are extended to study compressible and viscous flows around moving/deforming objects. A hybrid Lagrangian-Eulerian approach is designed to consider the motion of multiple bodies through a fixed Cartesian mesh. That is, the cells’ volumes do not move in space but rather they observe the solid walls crossing themselves. A dynamic discrete forcing makes use of a moving least-square procedure which has been validated by simulating well-known benchmarks available for rigid body motions. Partitioned fluid-structure interactions (FSI) strategies are explored to consider aeroelastic phenomena. A shared platform, between the aerodynamic and the structural solvers, fulfils the loads’ transfer and drives the sequence of the operating steps.

The first part of the results is devoted to a basic two-dimensional study aiming at evaluating the accuracy of the method when simple rigid motions are prescribed. Afterwards, the paper discusses the solution obtained when applying the dynamic IB method to the rigid deployment of a Krueger-flap. The final section discusses the aeroelastic behaviour of a three-element airfoil during its deployment phase. A loose FSI coupling is applied for estimating the possible loads’ downgrade.

The IB surfaces are allowed to move less than one IB-cell size at each time-step de-facto restricting the Courant-Friedrichs-Lewy (CFL) based on the wall velocity to be smaller than unity. The violation of this constraint would impair the explicit character of the method.

The proposed method improves automation in FSI numerical analysis and relaxes the human expertise/effort for meshing the computational domain around complex three-dimensional geometries. The logical consequence is an overall speed-up of the simulation process.

The value of the paper consists in demonstrating the applicability of dynamic IB techniques for studying high-lift devices. In particular, the proposed Cartesian method does not want to compete with body-conforming ones whose accuracy remains generally superior. Rather, the merit of this research is to propose a fast and automatic simulation system as a viable alternative to classic multi-block structured, chimaera or unstructured tools.

The purpose of this paper is to solve unsteady compressible Navier–Stokes equations without the commonly used dual-time loop. The authors would like to use an adaptive time-stepping (ATS)-based local error control instead of CFL-based time-stepping technique. Also, an all-speed flow algorithm is implemented with simple low dissipation AUSM convective scheme, which can be computed without preconditioning which in general destroys the time accuracy.

In transient flow computations, the time-step is generally determined from the CFL condition. In this paper, the authors demonstrate the usefulness of ATS based on local time-stepping previously used extensively in ordinary differential equations (ODE) integration. This method is implemented in an implicit framework to ensure the numerical domain of dependence always contains the physical domain of dependence.

In this paper, the authors limit their focus to capture the unsteady physics for three cases: Sod’s shock-tube problem, Stokes’ second problem and a circular cylinder. The use of ATS with local truncation error control enables the solver to use the maximum allowable time-step, for the prescribed tolerance of error. The algorithm is also capable of converging very rapidly to the steady state (if there is any) after the initial transient phase. The authors present here only the first-order time-stepping scheme. An algorithmic comparison is made between the proposed adaptive time-stepping method and the commonly used dual time-stepping approach that indicates the former will be more efficient.

The original method of ATS based on local error control is used extensively in ODE integration, whereas, this method is not so popular in the computational fluid dynamics (CFD) community. In this paper, the authors investigate its use in the unsteady CFD computations. The authors hope that it would provide CFD researchers with an algorithm based on an adaptive time-stepping approach for unsteady calculations.

This study aims to propose an aerodynamic force decomposition which, for the first time, allows for thrust/drag bookkeeping in two-dimensional viscous and unsteady flows. Lamb vector-based far-field methods are used at the scope, and the paper starts with extending recent steady compressible formulas to the unsteady regime.

Exact vortical force formulas are derived considering inertial or non-inertial frames, viscous or inviscid flows, fixed or moving bodies. Numerical applications to a NACA0012 airfoil oscillating in pure plunging motion are illustrated, considering subsonic and transonic flow regimes. The total force accuracy and sensitivity to the control volume size is first analysed, then the axial force is decomposed and results are compared to the inviscid force (thrust) and to the steady force (drag).

Two total axial force decompositions in thrust and drag contributions are proposed, providing satisfactory results. An additional force decomposition is also formulated, which is independent of the arbitrary pole appearing in vortical formulas. Numerical inaccuracies encountered in inertial reference frames are eliminated, and the extended formulation also allows obtaining an accurate force prediction in presence of shock waves.

No thrust/drag bookkeeping methodology was actually available for oscillating airfoils in viscous and compressible flows.

The purpose of this work is to investigate the similarity requirements for the application of multifidelity modeling (MFM) for the prediction of airfoil dynamic stall using computational fluid dynamics (CFD) simulations.

Dynamic stall is modeled using the unsteady Reynolds-averaged Navier–Stokes equations and Menter's shear stress transport turbulence model. Multifidelity models are created by varying the spatial and temporal discretizations. The effectiveness of the MFM method depends on the similarity between the high- (HF) and low-fidelity (LF) models. Their similarity is tested by computing the prediction error with respect to the HF model evaluations. The proposed approach is demonstrated on three airfoil shapes under deep dynamic stall at a Mach number 0.1 and Reynolds number 135,000.

The results show that varying the trust-region (TR) radius (λ) significantly affects the prediction accuracy of the MFM. The HF and LF simulation models hold similarity within small (λ ≤ 0.12) to medium (0.12 ≤ λ ≤ 0.23) TR radii producing a prediction error less than 5%, whereas for large TR radii (0.23 ≤ λ ≤ 0.41), the similarity is strongly affected by the time discretization and minimally by the spatial discretization.

The findings of this work present new knowledge for the construction of accurate MFMs for dynamic stall performance prediction using LF model spatial- and temporal discretization setup and the TR radius size. The approach used in this work is general and can be used for other unsteady applications involving CFD-based MFM and optimization.

For higher swirling flows (swirl > 0.5), flow confinement significantly impacts fluid flow, flame stability, flame length and heat transfer, especially when the confinement ratio is less than 9. Past numerical studies on helical axial swirler type systems are limited to non-reacting or reacting flows type Reynolds averaged Navier Stokes closure models, mostly are non-parametric studies. Effects of parametric studies like swirl angle and confinement on the unsteady flow field, either numerical or experimental, are very minimal. The purpose of this paper is to document modeling practices for a large eddy simulation (LES) type grid, predict the confinement effects of a single swirler lean direct injection (LDI) system and validate with literature data.

The first part of the paper discusses the approach followed for numerical modelling of LES with the minimum number of cells required across critical sections to capture the spectrum of turbulent energy with good accuracy. The numerical model includes all flow developing sections of the LDI swirler, right from the axial setting chamber to the exit of the flame tube, and its length is effectively modelled to match the experimental data. The computational model predicts unsteady features like vortex breakdown bubble, represented by a strong recirculation zone anchored downstream of the fuel nozzle. It is interesting to note that the LES is effective in predicting the secondary recirculation zones in the divergent section as well as at the corners of the tube wall.

The predictions of a single helical axial swirler with a vane tip angle of 60°, with a duct size of 2 × 2 square inches, are compared with the experimental data at several axial locations as well as with centerline data. Both mean and unsteady turbulent quantities obtained through the numerical simulations are validated with the experimental data (Cai et al., 2005). The methodology is extended to the confinements effect on mean flow characteristics. The time scale and length scale are useful parameters to get the desired results. The results show that with an increase in the confinement ratio, the recirculation length increases proportionally. A sample of three cases has been documented in this paper.

The novelty of the paper is the modelling practices (grid/unsteady models) for a parametric study of LDI are established, and the mean confinement effects are validated with experimental data. The spectrum of turbulent energies is well captured by LES, and trends are aligned with experimental data. The methodology can be extended to reacting flows also to study the effect of swirl angle, fuel injection on aerodynamics, droplet characteristics and emissions.

The purpose of this paper is to present the design, analysis and experiments of the active trailing-edge flap (Seoul National University Flap, SNUF) for vibration reduction in the helicopter rotor prior to the small-scale blades planned to test in a whirl tower.

The predictions of the hinge moment in both steady and unsteady flows were obtained through computational fluid dynamics calculations. When compared with the results originated from analytical formulations, the proposed method showed improved prediction capabilities. To validate the deflection of the flap under the centrifugal load by rotating, static analysis was conducted using both contact and rotating condition of MSC NASTRAN. The corresponding experiment also was performed using the vertical frame for simulating the effect of the centrifugal force.

The hinge moment of the flap is predicted through unsteady analysis in the actuation frequency of 3/rev. The material of the guide in the flap mechanism was selected through static analysis under both contact and rotating condition. Finally, reduction of the deflection occurred because of the load in the axial direction of the hinge like the centrifugal load.

The important aspects, such as design, analysis, and experiments for the active trailing-edge flap were shown.

This paper showed the relationship of the displacement, block force and voltage of the piezo-actuator, combined with the hinge moment predicted. The methodology and the experiment were presented for simulating the centrifugal force acting on the flap.

The purpose of this paper is to investigate airfoil self-noise generation and propagation by using a hybrid method based on the large-eddy simulation (LES) approach and Curle’s acoustic analogy as implemented in OpenFOAM.

Large-eddy simulation of near-field flow over a NACA6512-63 airfoil at zero angle of attack with a boundary layer trip at Re_c = 1.9 × 10⁵ has been carried out using the OpenFOAM® computational fluid dynamics (CFD) code. Calculated flow results are compared with published experimental data. The LES includes the wind tunnel installation effects by using appropriate inflow boundary conditions obtained from a RANS κ – ω SST model computation of the whole wind tunnel domain. Far-field noise prediction was achieved by an integral method based on Curle’s acoustic analogy. The predicted sound pressure levels are validated against the experimental data at various frequency ranges.

The numerical results presented in this paper show that the flow features around a NACA6512-63 airfoil have been correctly captured in OpenFOAM LES calculations. The mean surface pressure distributions and the local pressure peaks for the step trip setup agree very well with the experimental measurements. Aeroacoustic prediction using Curle’s analogy shows an overall agreement with the experimental data. The sound pressure level-frequency spectral analysis produces very similar data at low to medium frequency, whereas the experimentally observed levels are slightly over predicted at a higher frequency range.

This study has achieved and evaluated an alternative aeroacoustic simulation method based on the combination of LES with a simple Smagorinsky SGS model and Curle’s analogy, as implemented in the OpenFOAM CFD code. The unsteady velocity/pressure source data produced can be used for any simpler analytically based far-field noise prediction scheme.

A complete integration of the LES and Curle’s acoustic analogy for aeroacoustic simulations has been achieved in OpenFOAM. The capability and accuracy of the hybrid method are fully evaluated for high-camber airfoil self-noise predictions. Wind tunnel installation effects have been incorporated properly into the LES.