Recent advances in physical understanding and quantitative prediction of impinging-jet dynamics and atomization

Impinging-jet injectors are widely used in liquid propulsion applications, since their simple configuration provides reliable and efficient atomization. The flowfield involves a series of complicated spatio-temporal evolutions. Much effort has been directed toward understanding the underlying physics and developing quantitative predictions of impinging-jet atomization. This paper summarizes the recent advances in this direction, including state-of-the-art theoretical, experimental, and numerical studies, along with representative results. Finally, concluding remarks address remaining challenges and highlight modeling capabilities of high-fidelity simulations.     

Combining MHD and kinetic modelling of solar flares

Solar flares are explosive events in the solar corona, representing fast conversion of magnetic energy into thermal and kinetic energy, and hence radiation, due to magnetic reconnection. Modelling is essential for understanding and predicting these events. However, self-consistent modelling is extremely difficult due to the vast spatial and temporal scale separation between processes involving thermal plasma (normally considered using magnetohydrodynamic (MHD) approach) and non-thermal plasma (requiring a kinetic approach). In this mini-review we consider different approaches aimed at bridging the gap between fluid and kinetic modelling of solar flares. Two types of approaches are discussed: combined MHD/test-particle (MHDTP) models, which can be used for modelling the flaring corona with relatively small numbers of energetic particles, and hybrid fluid-kinetic methods, which can be used for modelling stronger events with higher numbers of energetic particles. Two specific examples are discussed in more detail: MHDTP models of magnetic reconnection and particle acceleration in kink-unstable twisted coronal loops, and a novel reduced-kinetic model of particle transport in converging magnetic fields.     

CFD predictions of LBO limits for aero-engine combustors using fuel iterative approximation

Lean blow-out (LBO) is critical to operational performance of combustion systems in propulsion and power generation. Current predictive tools for LBO limits are based on decadesold empirical correlations that have limited applicability for modern combustor designs. According to the Lefebvre’s model for LBO and classical perfect stirred reactor (PSR) concept, a load parameter (LP) is proposed for LBO analysis of aero-engine combustors in this paper. The parameters contained in load parameter are all estimated from the non-reacting flow field of a combustor that is obtained by numerical simulation. Additionally, based on the load parameter, a method of fuel iterative approximation (FIA) is proposed to predict the LBO limit of the combustor. Compared with experimental data for 19 combustors, it is found that load parameter can represent the actual combustion load of the combustor near LBO and have good relativity with LBO fuel/air ratio (FAR). The LBO FAR obtained by FIA shows good agreement with experimental data, the maximum prediction uncertainty of FIA is about ±17.5%. Because only the non-reacting flow is simulated, the time cost of the LBO limit prediction using FIA is relatively low (about 6 h for one combustor with computer equipment of CPU 2.66 GHz · 4 and 4 GB memory), showing that FIA is reliable and efficient to be used for practical applications.     

Numerical investigation on jet interaction with a compression ramp

A numerical investigation on jet interaction in supersonic laminar flow with a compres- sion ramp is performed utilizing the AUSMDV scheme and a parallel solver. Several parameters dominating the interference flowfield are studied after defining the relative increment of normal force and the jet amplification factor as the evaluation criterion of jet control performance. The computational results show that most features of the interaction flowfield between the transverse jet and the ramp are similar to those between a jet and a flat plate, except that the flow structures are more complicated and the low-pressure region behind the jet is less extensive. The relative force increment and the jet amplification factor both increase with the distance between the jet and the ramp shortening till quintuple jet diameters. Inconspicuous difference is observed between the jet-before-ramp and jet-on-ramp cases. The variation of the injection angle changes the extent of the separation region, the plateau pressure, and the peak pressure near the jet. In the present computational conditions, 120 is indicated relatively optimal among all the injection angles studied. For cold gas simulations, although little influence of the jet temperature on the pressure distribution near the jet is observed under the computation model and the flow parameters studied, reducing jet temperature somehow benefits the improvement of the normal force and the jet efficiency. When the pressure ratio of jet to freestream is fixed, the relative force increment varies little when increasing the freestream Mach number, while the jet amplification factor increases.     

气动噪声风洞试验技术发展概述

在大量查阅国内外相关技术文献的基础上,对风洞声学试验技术的发展历程进行了分析与总结,并对目前试验的两大关键技术领域——声学风洞与麦克风相阵列作了比较详细的介绍与讨论。通过对国内外气动噪声试验技术状况的调查研究,充分阐释了目前中国发展该项技术的必要性和迫切性,并提出了相应的建议。     

Delaying stall of morphing wing by periodic trailing-edge deflection

Morphing wings can improve aircraft performance during different flight phases. Recently research has focused on steady aerodynamic characteristics of the morphing wing with a flexible trailing-edge, and the unsteady aerodynamic and stall characteristics in the deflection process of the morphing wing are worthy further investigation. The effects of the angle of attack and deflection rate on aerodynamic characteristics were examined, and based on the aerodynamic characteristics of the morphing wing, a method was developed to delay stall by using the flexible periodic trailing-edge deflection. The numerical results show that the lift coefficients in the deflection process are smaller than those in the static situation at small angles of attack, and that the higher the deflection rate is, the smaller the lift coefficients will be. On the contrary, at large angles of attack, the lift coefficients are higher than those in the static case, and they become larger with the increase of the deflection rate. Further, the periodic deflection of the flexible trailing-edge with a small deflection amplitude and high deflection rate can increase lift coefficients at the critical stall angle.     

From geometry to CFD grids—An automated approach for conceptual design

The CEASIOM software developed in the EU-funded collaborative research project SimSAC generates stability and control data for preliminary aircraft design using different methods of varying fidelity. In order to obtain the aerodynamic derivatives by CFD, the aircraft geometry must be defined, computational meshes generated, and numerical parameters set for the flow solvers. An approach to automation of the process is discussed, involving geometry generation and mesh generation for inviscid as well as RANS flow models.     

Transonic aeroelastic simulation for instability searches and uncertainty analysis

In this paper the use of eigenvalue stability analysis of very large dimension aeroelastic numerical models arising from the exploitation of computational fluid dynamics is reviewed. A formulation based on a block reduction of the system Jacobian proves powerful to allow various numerical algorithms to be exploited, including frequency domain solvers, reconstruction of a term describing the fluid-structure interaction from the sparse data which incurs the main computational cost, and sampling to place the expensive samples where they are most needed. The stability formulation also allows non-deterministic analysis to be carried out very efficiently through the use of an approximate Newton solver. Finally, the system eigenvectors are exploited to produce nonlinear and parameterised reduced order models for computing limit cycle responses. The performance of the methods is illustrated with results from a number of academic and large dimension aircraft test cases.     

Numerical investigation of unsteady vortex breakdown past 80°/65° double-delta wing

An improved delayed detached eddy simulation （IDDES） method based on the k-x-SST （shear stress transport） turbulence model was applied to predict the unsteady vortex breakdown past an 80o/65o double-delta wing （DDW）, where the angles of attack （AOAs） range from 30° to 40°. Firstly, the IDDES model and the relative numerical methods were validated by simulating the massively separated flow around an NACA0021 straight wing at the AOA of 60°. The fluctuation properties of the lift and pressure coefficients were analyzed and compared with the available measurements. For the DDW case, the computations were compared with such mea-surements as the mean lift, drag, pitching moment, pressure coefficients and breakdown locations. Furthermore, the unsteady properties were investigated in detail, such as the frequencies of force and moments, pressure fluctuation on the upper surface, typical vortex breakdown patterns at three moments, and the distributions of kinetic turbulence energy at a stream wise section. Two dominated modes are observed, in which their Strouhal numbers are 1.0 at the AOAs of 30°, 32° and 34° and 0.7 at the AOAs of 36o, 38° and 40°. The breakdown vortex always moves upstream and downstream and its types change alternatively. Furthermore, the vortex can be identified as breakdown or not through the mean pressure, root mean square of pressure, or even through correlation analysis.     

基于浸入式边界方法的运动物体声辐射数值模拟

为了在声学计算中简化网格的处理,从而发展出适用于运动复杂边界问题模拟的计算格式,基于浸入式边界方法,提出了一种适用于无黏流动、无穿透边界条件的处理方法,并对一些简单的运动物体声辐射问题进行了数值模拟.数值结果与解析解得到了令人满意的对比,初步验证了所提出的研究方案的可行性.