二维谐振腔内声流的气-动格式模拟

在利用气-动格式模拟谐振腔内非线性声振荡行为的基础上建立了声流分析的气-动格式模型,并模拟了声谐振腔内的二维瞬时流场及平均流场(声流)。重点研究了有限振幅非线性情况下,激波对轴向速度、径向速度及流场的影响,进而分析了速度反转、涡流、不规则Rayleigh流等非线性现象形成的原因。     

基于标准格子的低马赫可压缩Boltzmann模型及应用于气体谐振的模拟

提出一种基于D2Q9标准格子的低马赫数可压缩格子Boltzmann模型,模拟了管内一维气体谐振。气体振荡由位于谐振管左端的活塞振动引起,传播至右端并反射回来。利用本文模型捕捉到了谐振频率附近不同频率下的激波现象,得到了谐振频率下不同位置处轴向速度、密度和温度随时间和沿轴向的变化规律。模拟结果和理论分析及前人文献结果相一致,验证了本文所提出的基于标准格子的低马赫数可压缩格子Boltzmann模型的有效性。     

热声起振现象的模拟研究

热声现象中广泛存在非线性效应,其中最典型的就是热声自激起振.本文建立了一个四分之三波长的Rijke管计算模型,基于并行计算,采用数值模拟对其起振现象进行了模拟.研究了不同初场对起振行为的影响,得出起振后的气体振荡频率,并对管内的流动与换热规律进行了详尽分析.模拟结果与文献结果相吻合.本文工作有助于更多地了解热声起振现象,同时也为我们开展热声问题的介观模拟提供了有力的参考.     

掺杂硅纳米梁谐振频率的理论模型及分子动力学模拟

通过理论计算与模拟,研究分析了P元素替代掺杂单晶硅纳米梁的谐振频率.计算模拟了两端固支单晶硅纳米梁的谐振频率随尺寸、掺杂浓度与温度的变化.通过对计算结果与模拟结果的分析得到:单晶硅纳米梁的谐振频率随着硅纳米梁长度尺寸的增大而减小;硅纳米梁的谐振频率随着掺杂浓度的增大而增大,但变化趋势并不明显;最后考虑了温度效应,发现掺杂硅纳米梁的谐振频率随着温度的增大而减小,但从谐振频率的数值来看,硅梁的谐振频率随温度的变化趋势并不明显,即温度对硅梁谐振频率基本无影响.由此得出结论:掺杂浓度与温度对硅纳米梁谐振频率的影响很小,影响单晶硅纳米梁谐振频率的主要因素是尺寸大小,掺杂单晶硅纳米梁的谐振频率具有尺寸效应.     

格子气的声振荡衰减

采用格子气数值方法模拟了共振腔中气体声振荡衰减过程．计算结果表明：格子气中声振荡衰减的规律和经典理论结果一致；衰减系数随频率增加而减小，并且随密度增加而减小．该结果可以用来说明一些大学统计物理教科书中的一个经典例子存在不完善之处．本中采用了一种新的扭转9-bit格子气模型和并行计算方法，这些结果对今后研究格子气中的热声问题有参考价值．     

无网格方法中的背景积分方案及单颗粒下降过程的数值模拟

在用无网格方法对有边界移动的气固两相流动进行数值模拟时发现,如果采用背景积分网格对求解域进行积分,将会产生压力场的数值振荡,使计算结果失真.详细分析了产生这种数值振荡的原因,提出了解决方案,并计算了一个颗粒在管道内的下降过程.计算结果表明,提出的解决方案能有效地降低压力场的数值振荡,从而使无网格方法能应用于气固两相流动的直接数值模拟.     

飞行器大攻角非定常气动特性神经网络建模

大攻角气动特性预测与气动建模是新型飞行器提升飞行性能的重要内容.以轴对称导弹简化模型为研究对象，首先采用计算流体力学方法，对70°大攻角状态的非定常气动特性进行数值模拟，计算方法基于RANS的N-S方程，湍流模型采用SA模型，对流场采用有限体积法离散，无黏项采用Roe通量差分分裂格式，黏性项采用中心差分，时间推进采用LU-SGS格式的双时间步法.飞行器运动模式采用强迫振荡的方式，对5种不同振荡频率进行了非定常数值计算，并记录每一内迭代周期最终的气动力和力矩数值.其次，以CFD预测结果作为气动建模的样本，采用动导数模型、多项式模型等传统方法，进行气动建模，并分析其有效性和精度.最后采用神经网络方法对大攻角非定常气动力进行建模，并和动导数模型、多项式模型进行精度对比.结果表明，基于神经网络的人工智能气动建模方法具有较高的精度和适应性.该方法为飞行器大攻角非定常非线性气动建模，大攻角飞行稳定性分析与控制提供理论参考.      

QUICK格式在湍流旋流流动数值模拟中的应用

为研究不同精度的离散格式对湍流旋流流动数值模拟结果的影响,同时应用QUICK格式和混合格式对同轴射流旋流燃烧室内的湍流流动进行了数值模拟.在采用k-ε湍流模型的条件下,QUICK格式计算得到的燃烧室内气体轴向与切向速度及轴向脉动速度均方根值分布与实验数据符合较好,而混合格式给出的数值模拟结果则与实验有一定的偏差.     

封闭圆内开缝圆自然对流的非线性特性研究

本文通过数值计算对封闭圆内开缝圆自然对流的非线性特性进行了研究。数值计算以整个流场为计算区域,采用非稳态数学模型和具有QUICK差分格式的SIMPLE算法。计算结果表明,在不同参数下流动和换热存在稳态定常解、周期性振荡解、拟周期性振荡解和非周期性振荡解。稳态定常解的相图是一个点;周期性振荡解的相图是一个极限环,对应功率谱含一个基频及其谐波;拟周期性振荡解的相图为环面,对应功率谱含两个不相关的频率及它们的线性组合频率;非周期性振荡解的功率谱为无规则的宽带连续谱。     

Rijke管自激热声振荡的数值模拟

为了对热声不稳定的发生及控制机理进行研究,对Rijke管内的自激热声振荡现象进行了数值模拟。采用具有低频散低耗散特点的计算气动声学方法,对带有非线性热源项的声波方程进行数值求解,并比较了不同的热源模型及边界条件对非线性效应的影响。结果表明,计算气动声学方法可以成功捕捉到Rijke管内压力的起振过程,而且在速度扰动达到平均流速度的1/3时,振荡会由线性增长转为非线性增长,最终达到有限幅值极限循环。相比热源项,考虑管口辐射耗散的非线性边界条件在振荡幅值和频谱方面对结果的影响都比较小。数值模拟得到的结果与实验符合较好,表明计算气动声学方法适合于热声振荡问题的研究。     

Nonlinear gas oscillations in an open tube under anharmonic excitation

Nonlinear gas oscillations excited in a tube with an open end by a piston driven by a crank mechanism are investigated. For the open end of the tube, a nonlinear boundary condition is formulated with allowance for oscillations at the subharmonic resonance frequency. Both first- and second-order approximations to the oscillations at the fundamental frequency and at half this frequency are calculated. The results of theoretical calculations are compared with experimental data.     

Mathematical simulation of a gas flow in the vicinity of the open end of a tube for harmonic oscillations of a piston at the resonance frequency at the other end of the tube

The structure of the gas flow in the vicinity of the open end of a tube for the oscillating gas flow caused by piston oscillations at the first resonance frequency at the other end of the tube has been determined by numerical integration of the Navier–Stokes equations using the ANSYS FLUENT program package. For the variant of the tube with an infinitely long flange and a sharp edge, the influence of the piston displacement amplitude on the gas flow rate in the tube is investigated, and the phases of gas inflow and outflow during the period of oscillation have been determined.     

Resonant oscillations of a gas in an open tube in the shock-free wave mode

Nonlinear gas oscillations excited in an open tube by a flat piston at one of the tube ends are studied. The sinusoidal piston oscillations in the shock-free wave mode are created by a vibration exciter near the first eigenfrequency. Expressions for gas pressure oscillations are obtained for a tube with a nonrounded end without a flange and secondary flow velocity components. The influence of the piston displacement amplitude on the pressure range and secondary flow velocity of gas is studied. The theoretical calculations of the gas pressure are compared with experimental data. An estimate for the velocity of particle motion along the tube axis is presented with calculated values of the secondary flow velocity.     

Particle drift in a resonance tube--a numerical study

The paper considers longitudinal drift of small particles in a resonance tube, caused by periodic shock waves, and its effect on particle agglomeration. It is found that depending on particle size, drift is caused by shock waves and/or gas acceleration and compression. It is also shown that the drift velocity and direction can be controlled by the frequency of the piston that causes gas oscillations in the resonance tube. The obtained numerical solutions indicate that particle drift in a resonance tube enhances aerosol agglomeration. An agglomeration kernel is derived for this case, accounting for particle drift, leading to an estimate of agglomeration time. The time predicted by present model is of the same order of magnitude as that obtained from experiments in the literature.     

Violent folding of a flame front in a flame-acoustic resonance

The first direct numerical simulations of violent flame folding because of the flame-acoustic resonance are performed. Flame propagates in a tube from an open end to a closed one. Acoustic amplitude becomes extremely large when the acoustic mode between the flame and the closed tube end comes in resonance with intrinsic flame oscillations. The acoustic oscillations produce an effective acceleration field at the flame front leading to a strong Rayleigh-Taylor instability during every second half period of the oscillations. The Rayleigh-Taylor instability makes the flame front strongly corrugated with elongated jets of heavy fuel mixture penetrating the burnt gas and even with pockets of unburned matter separated from the flame front.     

The evolution of nonlinear resonant oscillations in closed tubes: Solution of the standard mapping

The concern of this paper is the evolution of small-amplitude resonant oscillations of an inviscid gas in a closed tube. The evolution of the oscillations of the gas generated at and near the fundamental frequency and half the fundamental frequency, where experiment shows that shocks are a feature of the final periodic motion, is examined. The basis for the analysis of this nonlinear initial value, boundary value problem on a semi-infinite strip is the Dissipative Standard Map. Since the purpose is to elucidate how the initial rest state of the gas evolves to the final periodic state, the focus of the analysis is on how a prescribed initial curve is mapped, under the Dissipative Standard Map, onto an invariant curve. The method used is to approximate the Dissipative Standard Map, in the small rate limit, by the partial differential equation appropriate to the resonance in question.     

弹性微管内气泡的非线性受迫振动

将弹性管壁视为膜弹性结构, 探索在外部声场作用下弹性微管内液柱-气泡-管壁构成耦合振动系统的非线性特征. 利用逐级近似法对系统非线性共振频率、基频和三倍频振动幅值响应、 分频激励共振机理等进行了理论分析. 基频和三倍频振动的幅-频响应数值结果表明: 气泡的轴向共振和管壁共振不能同时出现; 两垂直方向的振动均表现出幅值响应多值性, 进而可能引起系统的不稳定声响应; 三倍频振动在低频区的声响应强于高频区. 关键词： 弹性微管 受迫振动 非线性振动 气泡声响应     

Hydrodynamics of resonance oscillations of columns of inelastic particles

We study oscillations of a one-dimensional (1D) column of N slightly inelastic particles, produced by a piston vibrating at one end of a closed tube. It is found that for large enough vibrational amplitudes of the piston, the column oscillates periodically with the period equal to the vibrational period. The oscillation patterns are governed by the shock waves propagating across the column. The averaged kinetic energy per particle is shown to be proportional to the square of the vibrational frequency, omega. This energy also strongly depends on the vibrational amplitude. The maximal value of this kinetic energy achievable by these external vibrations is found to be of order omega(2)L(2), where L is the total volume (length) of the tube free of particles. The above results on the column resonance oscillations are also predicted by a 3D hydrodynamic model of an inelastic granular gas.     

等截面驻波管内大振幅驻波场的实验研究

通过改进等截面驻波管实验系统,在1阶峰值共振频率激励下获得了182.1 dB大振幅驻波场,并对1~5阶峰值共振频率激励下的大振幅驻波场谐波饱和情况以及波形畸变进行了实验研究。研究发现,尽管1阶峰值共振频率激励下声压级已达到182.1 dB,但波形畸变最小,谐波并未表现出饱和现象,而3阶峰值共振频率激励下的大振幅驻波场表现出了饱和趋势。对谷值共振频率激励下获得的大振幅驻波场进行对比实验研究,发现谷值共振频率激励下,1阶谷值共振频率所获得的驻波场声压级最大,但波形畸变也最大。在相同声源驱动电压下,1阶峰值共振频率激励下获得的驻波场声压级始终大于1阶谷值共振频率激励下获得的驻波场声压级。由此可见,利用扬声器在等截面驻波管中获取大振幅驻波场,驻波管由1阶峰值共振频率激励较为合适。      

A new approach to treating pressure oscillations in combustion instability phenomena

Developing exact models of combustion instabilities is not an easy task to carry out and requires a great deal of time prior to obtaining success. The present study proposes a low-order model for pressure oscillations that does not require any knowledge of the systems, any new physical findings nor intricate details regarding its operating condition. This new approach is obtained using a Modified Van der Pol’s equation (MVDP) which is tuned by use of a Dual Extended Kalman Filter (DKEF) as a recursive estimator with perspectives in control by computer. This phenomenological model is used to predict the pressure signal from a variety of different combustors. Input data were taken from experimental cases such as a Rijke tube, a gas turbine and a liquid-fuel aero-engine combustor. Furthermore, a simulation considering high frequency oscillations to show the capability of the new approach is presented. In all cases, the results demonstrated the feasibility of applying the tractable model MVDP and DKEF running together to investigate pressure oscillations in practical cases.