圆柱容器液体晃动问题的数值计算

本文将二维VOF方法推广到三维情形,研究了圆柱容器遭受水平强迫振荡时的液体晃动问题。在数值模型中,采用了改进的施主与受主单元体方法计算流体在单元体之间的交换。同时对处理任意容器边界的部分单元体方法进行了改进,使之满足边界面上的速度边界条件。应用改进后的VOF方法得到的结果与试验结果符合较好,能够成功地数模出由于动力不稳定性引起的圆柱容器遭受纵向水平振荡时液体表面的旋转运动响应。     

Regularities between kinematic and aerodynamic characteristics of flexible membrane wing

The kinematic characteristics of flexible membrane wing have vital influences on its aerodynamic characteristics. To deeply explore the regularities between them, time-resolved aerodynamic forces and deformations at different aeroelastic parameters and angles of attack(α) were measured synchronously by wind tunnel experiments. The membrane motion can be mainly divided into two states at α > 0° with various lift-enhancement regularities: Deformed-Steady State(DSS)at pre-stall, and Dynamic Bala...     

圆柱贮箱内液体非线性稳态晃动实验及动力学特性分析

针对圆柱贮箱内液体非线性稳态晃动力和力矩动态信号的采集与分析问题,设计了一种能同时实现充液贮箱整体固定安装和液体动态晃动力、力矩精确测量的实验装置,给出了相应的实验数据采集、分析和处理方法。通过开展大量实验,将实验数据与CFD仿真结果进行细致对比后发现:在旋转晃动状态下,CFD仿真结果存在明显的计算偏差,不便被独立的用于完整揭示贮箱内液体晃动的真实特性。最后,归纳了贮箱内液体的非线性稳态晃动行为随横向外激励频率的变化规律,并对比分析了驻波晃动、拍振晃动、多类复杂的旋转晃动等多种不同形式的非线性稳态晃动行为对运载工具控制系统的任务执行效率和稳定性的影响。     

充液航天器质量特性计算方法研究

充液航天器的质量特性对飞行动力学建模和控制设计有着重要影响。质量特性测量装置通常只能测出干重状态下航天器的质量特性，而充液状态下的航天器质量特性只能靠计算给出，其中如何计及各贮箱内液体的惯性张量是此类工程计算的难点。首先采用势流理论推导了任意形状贮箱内液体质量特性的一般表达式，包括全充液情况和部分充液情况，其中对部分充液情况考虑了自由液面小幅晃动的影响；然后采用三维有限元方法建立了液体质量特性的数值计算方法，并通过一个验证算例确认了该计算方法的有效性；最后应用该方法计算了一个充液航天器的质量特性，并研究了不同的液体处理方法对计算结果的影响。研究表明，工程上通常采用的“固化液体”处理方法所计算出来的转动惯量明显偏大，而所提出的计算方法具有更好的理论正确性和工程适用性。     

多充液贮腔航天器耦合动力学与姿态控制

基于已有的质心面液体大幅晃动等效力学模型，采用拉格朗日方法，系统地建立了任意复杂激励下四贮腔航天器刚-液耦合动力学精确模型。由本文方法所得到的等效晃动力及晃动力矩与通过Flow3D计算出的相应结果吻合良好，验证了所建模型的可靠性。根据所建立的大幅晃动类充液航天器耦合系统动力学模型进一步研究了该类航天器在进行三轴稳定姿态机动时的主要动力学耦合特性，并设计了一类抑制姿态振荡的补偿控 制器。     

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.     

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.     

Effect of pitch angle on initial stage of a transport airplane ditching

Airworthiness regulations require that the transport airplane should be proved to ensure the survivability of the ditching for the passengers. The planned ditching of a transport airplane on the calm water is numerically simulated. The effect of pitch angle on the impact characteristics is especially investigated by a subscaled model. The Reynolds-averaged Navier-Stokes (RANS) equations of unsteady compressible flow are solved and the realizable j-e equations are employed to model the turbulence. The transformation of the air-water interface is tracked by volume of fluid (VOF) model. The motion of the rigid body is modeled by dynamic mesh method. The initial ditching stage of the transport airplane is analyzed in detail. The numerical results show that as the pitching angle increases, the maximal normal force decreases and the pitching motion becomes much gentler. The aft fuselage would be sucked down by the water and lead to pitching up, whereas the low horizontal tail prevents this trend. Consequently, the transport aircraft with low horizontal tail should ditch on the water at an angle between 10 and 12 as a recommendation.     

飞机抛物线飞行试验中的液体行为仿真

微重力环境下液体的流动特性和晃动行为是航天工程上颇为关心的问题,也是进行推进剂管理的基础。飞机抛物线飞行是创建微重力环境地面实验方法之一,可实现对一定时间内液体行为的观测,但其成本较高且试验时长受限。基于失重飞机的一般运动规律,对单次及多次抛物线飞行的任务剖面进行了设计,并以单次抛物线飞行为例、利用CFD软件对液体行为进行了数值仿真,分析了失重飞行段初始条件对微重力液体响应的影响。仿真结果可为失重飞机飞行任务剖面的设计和优化提供参考。     

BAFFLE EFFECT OF SLOSHING INDUCEDFLOWS IN RESPONSE TO AXIAL IMPULSE IN MICROGRAVITY

Inordertocarryoutscientificexperiments,someexperimentalspacecraftusecryogeniccooIingforobservationinstrumentationandtelescope,superconductingsensorsforgyroread-outandmaintainverylowtemperaturesnearabsolutezeroformechanicalstability.Theap-proachestobothcoollngandcontrolinvolvetheuseofsuperfluidliquidheliumn.Inthisstudy,theresponseofcryogenicsystemstosloshingdynamicsassociatedwithaxialimpulseisinvestigated.cryogenicheliumn,atatemperatureofl.8K,isconsideredinthisstudy.withitssu-perfluidbehavior,t…