炸药爆轰物质点法三维数值模拟

炸药爆轰以及多点起爆所产生的爆轰波汇聚问题很难应用有限元法进行模拟分析,尤其当网格发生畸变时,导致有限元法计算效率和数值精度严重下降,甚至无法得到正确结果.为此,该文应用显式积分算法的物质点法对炸药两点起爆和按时间序列的多点起爆的爆轰过程进行数值模拟,与炸药爆轰的理论计算结果相吻合.物质点法不但可以有效地避免网格畸变问题,而且为炸药爆轰的数值模拟提供了新的思路.     

超音速尖锥边界层中扰动演化特征的数值研究

采用高精度紧致格式, 对超音速尖锥边界层中二维扰动的空间演化, 进行了直接数值模拟．结果表明,虽然尖锥边界层流动存在一定的锥面法向速度,但小扰动的幅值及相位的演化都与由平行流假设得到的线性理论结果吻合．还研究了有限幅值扰动的演化,给出了其演化规律．并在扰动幅值增长到一定值时,发现了小激波．     

隐式TVD格式在气液两相爆轰数值模拟中的应用

为研究FAE(燃料-空气炸药)爆炸参数规律,运用二维轴对称气液两相方程组模型,针对其中的气相方程组采用高分辨率的隐式TVD格式,液相方程组采用MacCormack格式,较好地对FAE气液两相爆轰产生的爆轰波的发展和传播过程进行了数值模拟,计算结果与国内外的研究结果符合良好.     

超声速边界层中小幅值T-S波的数值研究

对来流马赫数Ma_∞=45的平板边界层中,幅值A分别为来流速度的0.01,0.001,0.0001倍的扰动波传播的物理过程进行了直接数值模拟。计算采用NND格式。模拟中发现即使扰动幅值尚小时,流场中即已出现小激波。     

一种治愈强激波数值不稳定性的混合方法

HLLC（Harten-Lax-Leer-contact）格式是一种高分辨率格式,能够准确捕捉激波、接触间断和稀疏波.但是使用HLLC格式计算多维问题时,在强激波附近会出现激波不稳定现象.FORCE（first-order centred）格式在强激波附近表现出很好的稳定性,并且其数值耗散比HLL（Harten-Lax-Leer）格式小.分析了HLLC格式和FORCE格式在特定流动条件下的稳定性,构造了HLLC-FORCE混合格式并且进一步结合开关函数来消除HLLC格式的激波不稳定现象.数值试验表明新构造的混合格式不仅能够消除HLLC格式的激波不稳定现象,还最大程度地保留HLLC格式高分辨率的优点.     

耗散比拟法数值求解航天飞机超声速粘性绕流问题

文中利用特殊系数矩阵分裂法和简化了的耗散比拟法数值求解了M_∞=7,Re=148000和攻角α=5°的航天飞机超声速粘性绕流问题.简化了的耗散比拟法既提高后体侧翼附近复杂流场数值解的分辨能力,又能减少数值计算工作量.     

带激波的一维非定常两相流动的数值模拟

本文将处理带激波的单相气体非定常流动问题的两种高分辨数值方法(随机取样法和二阶GRP有限差分法)推广应用于气固悬浮体(亦称含灰气体)两相情况,计算了含灰气体激波管中两相激波特性、波后流场结构及气固两相流动参数随时间的变化.数值结果表明:这两种方法均能给出带有尖锐间断阵面的两相激波松弛结构.二阶GRP方法在计算精度和机时耗用等方面优于随机取样法.     

转捩边界层中次生流向涡演化的数值研究

采用高精度直接数值模拟方法和高效的特征无反射边界条件,进行可压缩流转捩边界层中出现的次生流向涡演化的数值研究.精细的数值模拟结果清楚地揭示了转捩边界层的复杂流场中次生流向涡的形成和演化过程,探讨它对转捩至关重要的环状涡生成的影响,发现次生流向涡和主流向涡的共同作用形成环状涡的一种新机理.     

粉尘气体中激波绕射的数值研究*

本文针对稀相气固两相体系,选取双流体耦合模型,综合运用算子分裂技术和高精度高分辨率数值方法,研究了激波在粉尘气体中沿90°拐角的绕射特性,揭示了固相颗粒及其物性改变对激波绕射特征和波后流场结构的影响.     

可压缩均匀各向同性湍流的直接数值模拟

采用8阶精度的中心差分格式及7阶精度的迎风偏斜格式对Re_l = 72~153, M_t= 0.2~0.7的均匀各向同性湍流进行了直接数值模拟, 建立了湍流数据库. 与他人的计算结果吻合十分理想, 说明方法的有效性. 数值结果表明, 采用适当的迎风型差分格式可以克服起动问题(start-up problem)对湍流Mach数的限制, 提高可计算的湍流Mach数, 是可压湍流直接数值模拟的有效方法. 分析了压缩性效应对湍流统计量的影响, 发现压缩性使得湍动能的衰减加快. 探讨了可压湍流中微激波产生的机理, 对流场进行了标度律分析. 发现在本文的Reynolds数和湍流Mach数条件下, 流场中扩展自相似性仍然成立, 同时发现压缩性对标度指数影响不大.     

Numerical Investigation of the Dynamics of Thick-Walled Cylindrical Shells with a Helical Reinforcement

The Wilkins method is modified for numerically investigating the dynamics of thick-walled cylindrical shells with a helical orthotropy. The results of calculations by this algorithm agree well with known analytical solutions. Some new physical effects, which are absent in isotropic and cylindrically orthotropic shells, are investigated.__________Translated from Mekhanika Kompozitnykh Materialov, Vol. 41, No. 2, pp. 225–236, March–April, 2005.     

A Numerical and Asymptotic Investigation of Boundary-layer Wake Evolution

The evolution of four characteristically different velocityprofiles in the boundary layer at a trailing edge is studiedboth numerically and analytically; the complete wake flow iscomputed using a finite difference procedure and this is comparedwith asymptotic expansions of the flow in both the near wakeand the far wake. Two of the profiles have direct applicationto physical situations which include the flow in the wake ofa finite flat plate and the flow in the free shear layer downstreamof the point of separation of a supersonic stream from a planewall. The other two profiles demonstrate the structural changeswhich occur in the near wake when part of the initial profileis on the point of separation.     

直线爆炸所形成的柱面波的近似分析解

由直线爆炸所形成的柱面波,文献[1]得到了数值解。本文对于绝热指数γ=C_p/C_v较大的情形,取ε=1/γ2为小参数,用伸缩坐标法得到了问题的一级近似分析解。算例(γ=3)表明,近似解与准确解(数值解)十分符合。     

Numerical Simulation of the Effect of Deformation Rates on the Dynamic Strength of GFRP Cylindrical Shells

Mechanics of Composite Materials - A technique for the numerical analysis of nonlinear dynamic deformation and progressive failure of cylindrical glass-fiber plastic shells is developed with...     

悬浮流中柱状粒子取向分布函数的数值模拟

推导并数值模拟了悬浮流中柱状粒子取向的分布函数,并将其用于楔形流场计算,得到了取向分布函数与粒子方向角的关系。说明与短时间内形成的最可几角分布相比,稳定状态的分布变化幅度不大,主要变化在于右上区域的逆时针转动,速度梯度值大的地方,转动角度也越大;粒子最可几方向接近于流线方向:在流线方向上,随着极径减小,最可几角增大,同时它与流线夹角减小。     

关于定积分数值计算的探讨

本文概括地讨论了现有的等距内插求积公式,并从定积分定义出发建立了一个更一般的等距内插求积公式。