l1-based calibration of POD-Galerkin models of two-dimensional unsteady flows

This paper discusses a physics-informed methodology aimed at reconstructing efficiently the fluid state of a system. Herein, the generation of an accurate reduced order model of two-dimensional unsteady flows from data leverages on sparsity-promoting statistical learning techniques. The cornerstone of the approach is l₁ regularised regression, resulting in sparsely-connected models where only the important quadratic interactions between modes are retained. The original dynamical behaviour is reproduced at low computational costs, as few quadratic interactions need to be evaluated. The approach has two key features. First, interactions are selected systematically as a solution of a convex optimisation problem and no a priori assumptions on the physics of the flow are required. Second, the presence of a regularisation term improves the predictive performance of the original model, generally affected by noise and poor data quality. Test cases are for two-dimensional lid-driven cavity flows, at three values of the Reynolds number for which the motion is chaotic and energy interactions are scattered across the spectrum. It is found that: (A) the sparsification generates models maintaining the original accuracy level but with a lower number of active coefficients; this becomes more pronounced for increasing Reynolds numbers suggesting that extension of these techniques to real-life flow configurations is possible; (B) sparse models maintain a good temporal stability for predictions. The methodology is ready for more complex applications without modifications of the underlying theory, and the integration into a cyber-physical model is feasible.     

Recent experience with different methods of drag prediction

During the recent past the process of engineering design has changed entirely because of rapid developments in computational simulations. In aerodynamic design, computational fluid dynamics (CFD) methods are slowly superseding empirical methods and design engineers are spending more and more time applying CFD tools to analyze and predict the aerodynamic characteristics of vehicles. The purpose of this paper is to review recent experiences in CFD-based drag prediction with an emphasis on flow solutions governed by the Euler and Reynolds-averaged Navier–Stokes equations. For these types of flow solutions, various drag analysis methodologies are outlined and applied to determine the drag of components of as well as whole-body airplanes, helicopters, and ground-based vehicles at subsonic and transonic flow conditions. The review demonstrates that although significant progress has been made, CFD-based drag prediction still faces a number of hurdles that must be dealt with before it will become more widely accepted.     

太阳光球磁场特征物理量在质子事件短期预报中的应用

利用SOHO/MDI全日面纵向磁图, 计算了三个描述太阳活动区磁场复杂性和非势性的特征物理量, 即纵向磁场最大水平梯度Bz, 强梯度中性线长度L, 孤立奇点数目η. 为检验太阳光球磁场特征在质子事件短期预报中是否有效, 采用BP神经网络方法, 建立了基于这三个磁场特征物理量简单的太阳质子事件短期(24h)预报模型. 模 型在对2002年和2003年连续两年的样本检测中, 有很高的准确率(2002年和2003年 分别为90 %, 87.54 %)和较高的 质子事件报准率(2002年和2003年分别为60 %, 75 %),从而为光球磁场特征物理量作为质子事件预报的有效因子提供了依据.     

Influence of non-equilibrium reactions on the optimization of aerothrust aeroassisted maneuver with orbital change

The spaceplane is perspective vehicle due to wide maneuverability in comparison with a space capsule. Its maneuverability is expressed by the larger flight range and also by a possibility to rotate orbital inclination in the atmosphere by the aerodynamic and thrust forces. Orbital plane atmospheric rotation maneuvers can significantly reduce fuel costs compared to rocket-dynamic non-coplanar maneuver. However, this maneuver occurs at Mach numbers about 25, and such velocities lead to non-equilibrium chemical reactions in the shock wave. Such reactions change a physicochemical air property, and it affects aerodynamic coefficients. This paper investigates the influence of non-equilibrium reactions on the aerothrust aeroassisted maneuver with orbital change. The approach is to solve an optimization problem using the differential evolution algorithm with a temperature limitation. The spaceplane aerodynamic coefficients are determined by the numerical solution of the Reynolds-averaged Navier-Stokes equations. The aerodynamic calculations are conducted for the cases of perfect and non-equilibrium gases. A comparison of optimal trajectories, control laws, and fuel costs is made between models of perfect and non-equilibrium gases. The effect of a chemically reacting gas on the finite parameters is also evaluated using control laws obtained for a perfect gas.     

Aeroelastic simulation of the first 1.5-stage aeroengine fan at rotating stall

A massive parallel aeroelastic simulation platform has been built to investigate the first 1.5-stage fan of an aeroengine at rotating stall. The Computational Fluid Dynamics (CFD) solver and Computational Structural Dynamics (CSD) solver are coupled directly by non-matching mesh interfaces. The unsteady rotor/stator interaction is solved by the Sliding Mesh Interface method. The original rotor blades are shrouded by the midspan shrouds. An unshrouded fan is also created to investigate the effects of the midspan shrouds. Both the shrouded fan and unshrouded fan have stable aeroelasticity at the designed state. At rotating stall, the stalled region rotates at 30% of the rotor speed on the absolute reference frame. The energy spectrum of the rotating stalled flow is measured quantitatively. It shows that the first two order excitations are much stronger than the higher order excitations. In the flow of rotating stall, the fifth backward travelling wave mode of shrouded fan is resonated by the fifth excitation of the rotational stalled flow because the rotational speed of the stalled region coincides with the modal rotational speed, while for the unshrouded fan, the first bending mode is resonated by the second excitation of the rotational stalled flow, forming two waves in the circumference of the annulus blades. At rotating stall, the vibration of the shrouded blades is still under control but the vibration of the unshrouded blades is diverged and out of control. A novel tool, i.e., resonance map, is proposed to predict the resonance. It provides a perspective to explain the effects of midspan shrouds at a theoretical level, and it would also be helpful in the structural design of blades.     

Nacelle-airframe integration design method for blended-wing-body transport with podded engines

Strong shock waves and flow separation often occur during the integration of nacelle and airframe for blended-wing-bodies with podded engines. To address this problem, this paper presents an integration method with numerical simulations. The philosophy of channeling flow and avoiding the throat effect on the nacelle and airframe is established based on the analysis of flow interference in the initial configuration. A parametric integration design method is proposed from twodimensional plane to three-dimensional space with control mechanisms and selection principles of the key parameters determined by their influences. Results show that strong shock waves and flow separation can be successfully eliminated under the influence of both the reshaped channel and decelerated inflow below the nacelle. Supersonic regions around the nacelle are effectively reduced, concentrating mainly on the lip position. Thus, a significant cruise drag reduction(8.7%) is achieved though the pressure drag of the nacelle increases.     

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.     

受限控制直接分配新算法

对于三维目标的受限控制量分配问题,给出了直接控制分配方案的一种新的算法--相邻面搜索算法.直接控制分配方法的关键是找到期望目标向量与目标可达集合外表面的交点.该算法不需要确定目标可达集的所有表面,而是由目标可达集表面一个面出发,逐步确定相邻的面,直到找到与期望向量相交的面.算法放宽了原算法中控制效率矩阵任意三列线性无关的条件,并且能够在各种情况下得到最优分配结果.经验证,算法满足精确性与实时性要求.该算法的有效性在某多操纵面飞机的飞行控制系统仿真中得到了验证.     

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.