Inverse design of low boom configurations using proper orthogonal decomposition and augmented Burgers equation

Mitigation of sonic boom to an acceptable stage is a key point for the next generation of supersonic transports. Meanwhile, designing a supersonic aircraft with an ideal ground signature is always the focus of research on sonic boom reduction. This paper presents an inverse design approach to optimize the near-field signature of an aircraft, making it close to the shaped ideal ground signature after the propagation in the atmosphere. Using the Proper Orthogonal Decomposition(POD) method, a guessed input of augmented Burgers equation is inversely achieved. By multiple POD iterations, the guessed ground signatures successively approach the target ground signature until the convergence criteria is reached. Finally, the corresponding equivalent area distribution is calculated from the optimal near-field signature through the classical Whitham F-function theory. To validate this method, an optimization example of Lockheed Martin 1021 is demonstrated. The modified configuration has a fully shaped ground signature and achieves a drop of perceived loudness by 7.94 PLdB. This improvement is achieved via shaping the original near-field signature into wiggles and damping it by atmospheric attenuation. At last, a nonphysical ground signature is set as the target to test the robustness of this inverse design method and shows that this method is robust enough for various inputs.     

Solitary and periodic waves in magneto–rotating plasmas: Effect of the Kappa–Cairns distributed electrons

The nonlinear propagation of ion–acoustic (IA) waves in a magneto–rotating plasma is studied by considering the Kappa-Cairns electron distribution. Employing the perturbation scheme, Korteweg–de Vries equation is derived. It is seen that both positive and negative potential solitons can be supported in the present plasma model. The numerical results reveal that the Kappa-Cairns distributed electrons modify features of the electrostatic waves significantly. The effects of non–thermal parameters, plasma rotation frequency, ion temperature, and the wave propagation angle on electrostatic solitary wave structures are also discussed here. It is found that the plasma parameters considerably influence the propagation of IA waves in rotating plasmas. Furthermore, using the bifurcation theory of planar dynamical systems to the K-dV equation, we have presented the existence of solitary and periodic traveling waves. Our study may be helpful to understand the behavior of ion–acoustic wave in the rotating plasma.     

基于目标加速度方向观测的微分对策制导律

针对延迟信息条件下的机动目标拦截问题,研究了一种考虑目标加速度方向观测的微分对策制导律。首先,针对引入目标方向观测信息末制导博弈问题,建立了微分对策数学模型。随后,利用状态依赖黎卡提方程求解微分对策问题,得到包含无延迟目标机动项的微分对策制导律,可以更为有效地拦截机动目标。最后,在延迟信息条件下,利用目标加速度方向观测的方法补偿延迟的目标加速度,再将补偿后的目标加速度信息取代无延迟的目标加速度,得到延迟信息条件下考虑目标加速度方向观测的微分对策制导律。仿真结果表明,新的制导律可以在延迟信息条件下有效地拦截机动目标。     

Positron nonextensivity contributions on the rational solitonic,periodic, dissipative structures for MKP equation described critical plasmas

New closed forms have been exposed that rational, trigonometric, periodical, explosive, hyperbolic and shock solutions can be usable in the ionosphere plasma of earth. To explore the nonextensive impacts on the features of nonlinear waves in this plasma model, using Riccati-Bernoulli sub-ODE process, the MKP equation has been solved. This method is very important in the study of dynamics and motion in fluids. Some of the obtained new potential solutions are prefect achievements in plasma observations and applications in ionosphere.     

Reevaluation of thermosphere heating by auroral electrons

This paper presents a new calculation of neutral gas heating by precipitating auroral electrons. It is found that the heating rate of the neutral gas is significantly lower than previous determinations below 200 km altitude. The neutral gas heating arises from the many exothermic chemical reactions that take place from the ions and excited species created by the energetic electrons. The calculations show that less than half the energy initially deposited ends up heating the neutral gases. The rest is radiated or lost in the dissociation of O₂ because the O atoms do not recombine in the thermosphere. This paper also presents a new way of calculating the heating rate per ionization that can be used for efficient determination of the overall neutral gas heating for global thermosphere models. The heating rates are relatively insensitive to the neutral atmosphere when plotted against pressure rather than altitude coordinates. At high altitudes, the heating rates are sensitive to the thermal electron density and long-lived species. The calculations were performed with the Field Line Interhemispheric Plasma (FLIP) model using a 2-stream auroral electron precipitation model. The heating rate calculations in this paper differ from previous heating rate calculations in the treatment of backscattered electrons to produce better agreement with observed flux spectra. This paper shows that more realistic model auroral electron spectra can be obtained by reflecting the up going flux back to the ionosphere at the upper boundary of the model. In this case, the neutral gas heating rates are 20%–25% higher than when the backscattered flux escapes from the ionosphere.     

Consequences of the application of the streamer fluid model to the study of the sprite inception mechanism

We present the results of a streamer-fluid model used to investigate the electrodynamical coupling between the troposphere and upper atmosphere due to the penetration of lightning electric fields into the mesosphere and the lower ionosphere, generating sprites. The model solves the continuity equation for electrons and ions coupled to Poisson equation. The dominant physical response of the atmosphere is the formation of a screening-ionization wave. The wave shields the atmosphere above it from the action of the lightning field and, together with the conductivity reduction below it due to attachment, the wave amplifies the total field below it, allowing for the penetration of intense electric fields in the mesosphere as it propagates downwards into regions of higher density that compress the wave. This is the key physical mechanism for sprite inception. We evaluated the effects of the thundercloud charge geometry, lightning current waveshape, atmospheric conductivity, via different electron density profiles, and the effect of ionization, attachment and electron mobility coefficients in the electrical breakdown process, related to halo production, and sprite streamer initiation. The results showed that electrons with higher mobility are more efficient in shielding the lightning electric field before breakdown, causing delay, and they contribute to the formation of the streamer seed after breakdown, anticipating the sprite streamer inception. Similarly, a higher effective ionization rate, produced by modifications in the attachment and ionization coefficients, anticipates sprite inception. The simulations with 6 different electron density profiles, and therefore conductivities, spanning 4 orders of magnitude, showed that the altitude of breakdown and sprite initiation, as well as their time delays from the lightning discharge are directly related to atmospheric conductivity: higher conductivities produce halo and sprite inception at lower altitudes with longer delays and may hinder sprite formation. We document that variations of 30 times in the lightning current leads to sprite initiation altitudes in the range 66.0–73.5 km, with delays between 1.550 and 34.500 ms, while variations of 4 orders of magnitude in the conductivity profile lead to initiation altitudes 61.0–70.6 km, with delays in the range 3.825–9.825 ms. Consequently, we suggest that lightning characteristics dominate over atmospheric parameters in determining sprites’ initiation altitude and delay. The simulation of a −CG, with a constant current of 30 kA, did not produce a sprite seed, confirming an asymmetry in the response of the atmosphere to positive and negative lightning. This is due to the free electron drift direction that is away from the screening ionization wave, preventing the formation of the streamer seed for the great majority of −CGs. The same does not apply to halos, which depend on the occurrence of breakdown and can be produced by discharges of both polarities.     

一种滤波SST方法在翼型深失速模拟中的应用

 为了提高原始剪切应力输运(SST)湍流模型对于分离流动的求解精度,将大涡模拟(LES)中的滤波因子和SST方程相结合构造出一种滤波SST方法,利用湍流尺度对流场求解区域进行划分,近壁面附近的稳态流动由湍流模型控制,远壁面采用LES方法进行模拟。与传统混合RANS/LES方法相比,该方法的特点是:滤波因子的选取不再依赖于网格尺度,可以有效地降低网格诱导分离现象发生的概率。采用该方法对NACA0021翼型深失速特性进行了仿真研究,对比了非定常雷诺平均Navier-Stokes(URANS)方法和SST-DES方法,从仿真结果可以看出滤波SST方法有效地激活了分离区域的脉动,充分展现了分离的三维特性;同时算例求解结果证明该方法的精度高于URANS方法,与试验结果吻合较好,显示其具有一定的工程应用价值。     

T型尾翼风洞颤振实验保护装置绕流特性分析

采用SST两方程湍流模型,通过求解非定常Navier-Stokes（N-S）方程,对T型尾翼风洞实验流场进行了模拟,分析了保护装置对T型尾翼风洞实验流场的影响,研究了保护装置几何外形和保护装置基座后移距离对流场影响。通过对平尾气动力的分析以及对非定常流场的对比,可以得出：采用NACA系列翼型对基座进行气动整流后,基座两侧局部超声速区显著减小,局部激波减弱甚至消失,流场品质得到改善。且采用NACA0010翼型对基座修形后的结果最理想。随着保护装置基座后移距离的增加,平尾气动力均方根值和波动值先是急剧减小,达到0.85倍平均气动弦长后开始有所增大,在2.45～4.05倍平均气动弦长范围基本不再变化,稳定到单独T型尾翼模型相应系数1倍左右。此结论对T型尾翼风洞颤振实验保护装置设计具有一定的指导意义。     

Satellite formation reconfiguration and station-keeping using state-dependent Riccati equation technique

The current paper presents optimal reconfigurations and formation-keeping for formation flying satellites. The state-dependent Riccati equation (SDRE) technique is utilized as a non-linear controller for both the reconfiguration problem and formation-keeping problem. For the SDRE controller, a state-dependent coefficient (SDC) form is formulated to include non-linearities in the relative dynamics and J₂ orbital perturbation. The Taylor series and a transformation matrix are used to establish the SDC form. Optimal reconfiguration trajectories that minimize energy in satellite formation flying are obtained by the SDRE controller and compared with those obtained from a linear quadratic regulator (LQR) and a linear parameter varying (LPV) control method. It is illustrated that the SDRE non-linear controller of the current study obtains relocation accuracy of less than 0.1% of formation base-line length, while the LQR controller and LPV controller yield relatively large relocation errors. The formation-keeping controller developed using the SDRE technique in the current study also provides robustness under severe orbital perturbations.     

基于节点电压增量方程的电路故障子网络诊断

基于节点电压增量方程的电路故障诊断方法通常用来诊断线性电路中的元器件级的故障情况.本文将其推广应用于诊断大规模线性模拟电路中子网络级故障以简化测前仿真,提高诊断效率.