球壳结构馈通增长的瑞利-泰勒不稳定性

利用小扰动分析法 ,导出不可压缩球壳结构的馈通增长方程 ,数值模拟了高压气体驱动外表面有初始扰动的明胶球壳的瑞利 泰勒不稳定性模型。计算结果表明 :对于低波数扰动 ,外界面比较稳定 ,内表面的馈通增长较快 ,具有比较明显的三个演化阶段和波形反转现象。高波数扰动的增长恰好与低波数相反。球壳会聚结构比柱壳会聚结构的界面稳定性要好些。     

NUMERICAL STUDY ON THE EFFECTS OF SMALL-AMPLITUDE INITIAL PERTURBATIONS ON RM INSTABILITY

采用Navier-Stokes 方程对入射激波及其反射激波连续诱导小振幅扰动界面的Richtmyer-Meshkov 不稳定性增长过程进行了二维数值模拟,分析了单模和随机多模两种不同初始形态的界面上钉结构和泡结构在反射激波作用前后的发展特性. 研究结果发现：单模扰动的初始界面形态对反射激波前、后界面的扰动增长都有影响,反射激波前的界面形态信息可以通过钉和泡结构之间的反转传递到反射激波过后. 扰动界面上钉结构的发展速度控制了界面混合区总体的发展速度,反射激波前界面上发展成具有完整冠部形态的钉,在反射激波后会反转成复杂的泡结构,此泡结构对反射激波后钉的发展不利. 随机多模界面显示了与单模界面类似的发展规律,但随机多模界面上的复杂泡结构分布的不对称性使得其对钉结构增长的拖曳效应相对要弱,这导致了相似扰动波长下多模随机界面的扰动发展相对单模界面扰动发展要快.     

小扰动界面形态对RM不稳定性影响的数值分析

采用Navier-Stokes 方程对入射激波及其反射激波连续诱导小振幅扰动界面的Richtmyer-Meshkov 不稳定性增长过程进行了二维数值模拟,分析了单模和随机多模两种不同初始形态的界面上钉结构和泡结构在反射激波作用前后的发展特性. 研究结果发现:单模扰动的初始界面形态对反射激波前、后界面的扰动增长都有影响,反射激波前的界面形态信息可以通过钉和泡结构之间的反转传递到反射激波过后. 扰动界面上钉结构的发展速度控制了界面混合区总体的发展速度,反射激波前界面上发展成具有完整冠部形态的钉,在反射激波后会反转成复杂的泡结构,此泡结构对反射激波后钉的发展不利. 随机多模界面显示了与单模界面类似的发展规律,但随机多模界面上的复杂泡结构分布的不对称性使得其对钉结构增长的拖曳效应相对要弱,这导致了相似扰动波长下多模随机界面的扰动发展相对单模界面扰动发展要快.      

入口扰动频率对矩形射流影响的大涡模拟

为研究入口扰动对矩形射流的影响,对不同扰动频率下的三维矩形射流进行了大涡模拟,用分布投影法求解动量方程,对压力泊松方程采用FACR直接解法.计算结果预报的速度在横向的单峰分布和速度半宽沿流向的近似线性增长与实验结果吻合良好.计算结果表明扰动频率为射流的本征频率fe=0.22时射流的卷吸作用最强,而高频扰动则抑制卷吸作用;激励频率接近本征频率时射流提前出现速度在展向的马鞍型双峰分布,更早趋于圆射流形态.进一步的流场二维和三维瞬态涡结构图谱表明:上述现象与f=0.1和f=0.22时射流近场更频繁卷起、更丰富的大涡结构存在密切关系,正是这些大涡运动不断地在横向卷吸入外界流体,实现射流在横向的迅速扩展.     

圆筒内旋转细长管与流体界面耦合的数值方法

圆筒内旋转细长管是石油钻采工程中特有结构,细长管不仅与圆筒发生碰撞接触,还与管内流体和管外环空流体耦合,是一个复杂的非线性流固耦合系统。细长管固体域离散成梁单元,采用非线性碰撞接触动力学方程描述;管内外流体离散成六面体单元,采用计算流体动力学方程描述,在耦合界面处用任意拉格朗日欧拉法动网格来处理运动界面。根据梁单元位移...     

蜿蜒边界下平面流的线性流动稳定性

为便于数值分析,蜿蜒河流水动力和演变模型中一般隐性假设二次时均流-二次涡的关系与明渠流时均流-明渠湍流的关系相同,但由于高雷诺数下的DNS算力限制和实验尺度限制,这种隐含假设是否成立目前尚无相关湍流研究来支撑.文章试图通过分析明渠湍流和二次湍流发展初期的研究,侧面揭示其湍流结构的异同.通过对曲线正交坐标系下的平面二维NS方程使用双参数摄动的方法,建立了一种求解蜿蜒边界弱非线性层流的摄动解法,并推导得出一个适用于蜿蜒边界的EOS方程以及其特征值问题的解法.蜿蜒边界下弱非线性层流解为一系列蜿蜒谐波分量的叠加,其中线性部分使得两壁产生流速差,非线性部分随着雷诺数增大呈指数增长.水流的扰动增长率特征谱的第一模态与直道流相似,由3条曲线、4个波段合成,但其长波段和短波段的扰动流场与直道流不同,所有短波段的扰动流速近似于KH涡.蜿蜒边界对内部水流扰动有一定的选择性.偏角幅值越大扰动增长越快;蜿蜒波数的影响则为先增后减,有一个使扰动增长最快的蜿蜒波数.扰动流场由一个典型的TS波和一对波包形式的二次涡叠加而成,波包只有纵向流速分量,包络线由蜿蜒波数控制,波包内是与直道扰动波参数相同的TS波.     

临界爆轰的稳定条件和螺旋爆轰波

本文应用耗散结构热力学理论和微扰的方法,研究临界爆轰波对二维扰动的稳定性问题。在忽略粘性、热传导、浓度扩散、外力和交叉效应的情况下,作出了扰动在爆轰波结构内传播的稳定性判据。对Arrhenius反应率,当化学反应活化能E大于某一临界值E_c之后,对反应率随温度增加而增加的放热反应,临界爆轰波对二维扰动是不稳定的,扰动的振幅在反应燃烧区内随时间的增加而增长,直至扰动最后离开爆轰波结构进入波后产物区。当考虑了粘性的影响之后,扰动振幅的增长和衰减依赖于扰动本身的频率,在反应放热量超过某临界值后,频率愈小亦即扰动波长愈长的扰动振幅,随时间增长愈快,以致最长波长的扰动增长掩盖了其它波长的扰动,或者只有最长波长的扰动振幅维持不变,其它波长的扰动振幅都逐渐衰减,最后形成有规则的螺旋爆轰波。所得结果当忽略化学反应以及粘性对扰动传播相速度的微小影响之后,结论与N.Monson和J.A.Fay等所作的声波理论结果相一致,比较成功地解释了螺旋爆轰的一些实验现象。     

内爆炸加载条件下圆筒的膨胀、破裂规律研究

以经典热粘塑性本构关系 (John Cook本构关系和Powerlaw本构关系 )为基础 ,建立了更适当表征材料特征的新的、解耦的本构关系。采用SHPB(分离式霍普金森杆 )技术测定了常用弹箭材料 35CrMnSiA的不同应变率下的应力、应变关系 ,并采用拟合的方法确定了本构方程中的材料常数。与高速摄影技术测得的径向位移函数及圆筒材料的损伤演化方程相结合 ,建立了控制内爆炸加载圆筒膨胀 ,直至破裂的完备方程组 ,完成了内爆炸加载圆筒问题的数值模拟 ,计算结果与以往有关箭弹材料圆筒膨胀的实验结果符合较好。     

与温度相关的非均质圆筒的非线性耦合广义热弹性分析

基于能量守恒方程和L–S广义热弹性理论,借助状态空间技术和Newmark法求解了材料性质沿径向任意梯度分布同时又与温度相关的非均质圆筒非线性耦合广义热弹性问题。通过对材料性质与温度无关和相关功能梯度圆筒的算例分析,给出了在线性和非线性耦合下圆筒内温度和应力沿径向和随时间的变化关系,验证了本文解的正确性和有效性。数值结果表明,考虑材料性质是否与温度相关,能量守恒方程中耦合项是线性还是非线性,得到的温度与应力均存在不同程度的差异。本文解可方便地应用于不同边界条件和初始条件下圆筒的广义热弹性分析。     

圆筒内横向流体与弹性管耦合的分域求解算法

海洋丛式井组隔水管和换热器管束等在流体作用下,会诱导管束振动及碰撞,导致管束断裂失效。将弹性管离散成梁单元,采用非线性结构动力学方程描述;将圆筒流体域离散成实体单元,采用计算流体动力学方程描述。在流固耦合界面处,推导了界面位移、速度和载荷计算公式及收敛判断准则,建立了圆筒内横向流体与弹性管耦合的分域求解算法。算例表明,其分域与全域求解计算结果吻合较好,本文算法为复杂流体域内多根管束的振动及碰撞问题求解提供行之有效的计算方法。     

Lie symmetries,symmetrical perturbation and a new adiabatic invariant for disturbed nonholonomic systems

For a nonlinear nonholonomic constrained mechanical system with the action of small forces of perturbation, Lie symmetries, symmetrical perturbation, and a new type of non-Noether adiabatic invariants are presented in general infinitesimal transformation of Lie groups. Based on the invariance of the equations of motion for the system under general infinitesimal transformation of Lie groups, the Lie symmetrical determining equations, constraints restriction equations, additional restriction equations, and exact invariants of the system are given. Then, under the action of small forces of perturbation, the determining equations, constraints restriction equations, and additional restriction equations of the Lie symmetrical perturbation are obtained, and adiabatic invariants of the Lie symmetrical perturbation, the weakly Lie symmetrical perturbation, and the strongly Lie symmetrical perturbation for the disturbed nonholonomic system are obtained, respectively. Furthermore, a set of non-Noether exact invariants and adiabatic invariants are given in the special infinitesimal transformations. Finally, one example is given to illustrate the application of the method and results.     

On localization in saturated porous continua

Localization in elastic-plastic saturated porous media is investigated here using a linear perturbation approach. The adopted localization criterion corresponds to unbounded growth of perturbations. The critical conditions are compared with those obtained by a classical band analysis. While for one phase materials these conditions coincide, in the present context the linear perturbation approach leads in the limit of unbounded growth to the singularity of the undrained acoustic tensor, while the band analysis leads to the singularity of the drained acoustic tensor. Some general results clarifying the hierarchy of these two conditions are established for a quite general set of constitutive equations.     

Application of the TVD scheme to the nonlinear instability analysis of a capillary jet

In the past, when either the perturbation‐type method or direct‐simulation approach was used to analyse capillary jets, the governing equations, which are parabolic in time and elliptic in space, were simplified or linearized. In the present study, the convective derivative term and a full, nonlinear form of the capillary pressure term are retained in the governing equations to investigate nonlinear effects on the break‐up of capillary jets. In this work, the TVD (i.e. total variation diminishing) scheme with flux‐vector splitting is applied to obtain the solutions of the system of nonlinear equations in a matrix form. Numerical results show that the present nonlinear model predicts longer jet break‐up lengths and slower growth rates for capillary jets than the previous linear model does. Comparing with other measurements from past literatures, the nonlinear results are consistent with the experimental data and appear more accurate than the linear analysis. In the past, the classic perturbation‐type analyses assumed constant growth rates for the fundamental and all harmonic components. By contrast, the present model is able to capture the local features of growth rates, which are not spatially and temporally constant. Copyright © 2006 John Wiley & Sons, Ltd.     

剪切湍流大尺度相干结构的模式研究

发展了一种计算剪切湍流大尺度相干结构的新模式．该模式的基础是认为大尺度相干结构为湍流场中流体脉动能量增长最快的那部分，且包含大部分的湍流脉动能量．在此基础上。通过对湍流相干能量方程的推演。建立了描述大尺度相干结构的特征控制方程，并应用Ｃｈｅｂｙｓｈｅｖ多项式方法求得湍流相干能量的最大增长率在波数空间的分布，从而获得对应的大尺度相干结构．应用该模式研究了槽流和一自然对流中的大尺度相干结构，得到的近壁区流动结构与实验现象十分接近．     

Numerical study of three-dimensional developments of premixed flame induced by multiple shock waves

The three-dimensional interactions of a perturbed premixed flame interface with a planar incident shock wave and its reflected shock waves are numerically simulated by solving the compressible, reactive Navier–Stokes equations with the high-resolution scheme and a single-step chemical reaction. The effects of the initial incident shock wave strength (Mach number) and the initial perturbation pattern of interface on the interactions are investigated. The distinct properties of perturbation growth on the flame interface during the interactions are presented. Our results show that perturbation growth is mainly attributed to the flame stretching and propagation. The flame stretching is associated with the larger-scale vortical flow due to Richtmyer–Meshkov instability while the flame propagation is due to the chemical reaction. The mixing properties of unburned/burned gases on both sides of the flame are quantitatively analyzed by using integral and statistical diagnostics. The results show that the large-scale flow due to the vortical motion always plays a dominating role during the reactive interaction process; however, the effect of chemistry becomes more important at the later stage of the interactions, especially for higher Mach number cases. The scalar dissipation due to the molecular diffusion is always small in the present study and can be negligible.     

Shock waves in diatomic chains—II. nonlinear analysis

The effectsof nonlinear terms in the force-displacement relations of a general diatomic chain are studied by methods closely analogous to the analysis of a monatomic chain by J. Tasi (1972, 1973). Two-time parameter perturbation solutions of the equations of motion are obtained for the nonlinear growth of the shock response at the head of the wave. These solutions are similar in form to the asymptotic results obtained by M.J.P. Musgrave and J. Tasi (1976) but depend upon a dispersive wave coordinate modified by the nonlinear parameters. Partial differential equations are derived for the far-field velocity response of each typical particle in a unit cell. In the ease of a general diatomic chain, these equations are not identical and discontinuity in strain within the unit cell occurs. Graphical plots of solutions are presented.     

大间隙环流中偏心转子动特性系数的数值分析方法

基于作者用整体流动理论和Moody壁面摩擦系数方程建立的大间隙环流中转子动特性系数数值计算模型,应用摄动方法推导了大间隙环流流场非线性控制方程组的一阶摄动方程,提出了求解大间隙环流中偏心转子动力学特性系数的数值分析方法。用该方法得到的数值结果与已有的解析解和实验结果具有较好的一致性。     

On methods for continuous systems with quadratic and cubic nonlinearities

Methods for determining the response of continuous systems with quadratic and cubic nonlinearities are discussed. We show by means of a simple example that perturbation and computational methods based on first discretizing the systems may lead to erroncous results whereas perturbation methods that attack directly the nonlinear partial-differential equations and boundary conditions avoid the pitfalls associated with the analysis of the discretized systems. We describe a perturbation technique that applies either the method of multiple scales or the method of averaging to the Lagrangian of the system rather than the partial-differential equations and boundary conditions.     

Reductive Perturbation Method of Super KdV Equations

ＩｎｔｒｏｄｕｃｔｉｏｎＳｕｐｅｒＫｄＶｅｑｕａｔｉｏｎｓａｒｅｕｔ-ｂｕｕｘ 3ｈｈｘ ｕｘｘｘ =0 ,ｈｔ-ｂｕｈｘ-ａｈｕｘ ｃｈｘｘ =0 ,( 1 )ｗｈｅｒｅａ ,ｂ,ｃ(ｃ≠ 0 )ａｒｅｃｏｎｓｔａｎｔｓ.Ｔｈｅｓｉｍｉｌａｒｓｏｌｕｔｉｏｎｓｏｆ ( 1 )ｗｅｒｅｇｉｖｅｎｉｎ [1 ]ｂｙｄｉｒｅｃｔｍｅｔｈｏｄｐｒｅｓｅｎｔｅｄｂｙＣｌａｒｋｓｏｎａｎｄＫｒｕｓｋａｌ.Ｉｎｔｈｉｓｐａｐｅｒ,( 1 )ａｒｅｃｈａｎｇｅｄｉｎｔｏｏｒｄｉｎａｒｙＫｄＶｅｑｕａｔｉｏｎｓｂｙｒｅｄｕｃｔｉｖｅｐｅｒｔｕｒｂａｔｉｏｎｍ…