微重力下液封液桥内振荡热毛细对流的数值模拟

对微重力条件下液封液桥(流体采用KF-96硅油和FC-70氟化物,桥高1.4~4毫米,内层直径2或3毫米,外层直径5毫米)和单层液桥(流体采用KF-96硅油,桥高1~1.6毫米,直径2毫米)内的热毛细对流进行了数值模拟,模拟条件与Majima等的实验条件相同。从理论上证实了大温差条件下将出现振荡热毛细对流,确定了发生振荡的临界条件并与相关实验结果进行了对比,同时还计算了振荡频率。     

复杂界面(界面层)条件下的弹性波传播问题研究综述

界面（界面层）广泛存在于各种材料和结构中,并具有千差万别的形态．各种界面条件下的波动问题有着重要的理论与实际意义．综述了完好粘接界面、弱连接界面（或界面层）以及接触界面等复杂界面模型下弹性波传播问题的研究现状,主要集中于界面模型的建立、波传播问题的研究方法及主要结论．并提出了值得进一步研究的问题．     

气相爆轰波冲击气固界面的透反射特性

为了研究气相爆轰波冲击气固界面过程中透射波和反射波的相关特性,建立爆轰波冲击气固界面的一维理论模型,对不同初始压力条件下爆轰波到达气固界面后的界面两侧的压力和界面速度变化进行分析。利用时空守恒元求解元方法对气相爆轰波冲击气固界面过程进行数值模拟,分析气体部分反射波的压力分布和速度变化规律及透射入固体中应力波的波形和波速特征,并搭建气相爆轰波冲击活塞实验装置进行进一步验证。结果表明：气体爆轰波到达气固界面后,在固体中透射指数形式的弹性波,并在界面处向气体区反射一道激波。爆轰波后的稀疏波与反射激波相交,削弱反射激波,最终形成稳定激波回传。气固界面在稀疏波和反射稀疏波的作用下,压力和速度逐渐下降,最终也形成稳定状态。在不同混气初始压力情况下,爆轰波冲击过程中产生的最高压力和爆压的比值基本保持不变。理论模型对特征点相关物理量的计算值和实验数据符合的较好。

     

爆炸焊接界面波的数值模拟

借助动力学分析软件ANSYS/LS-DYNA 11.0,运用光滑粒子流体动力学方法(SPH方法),建立以Johnson-Cook材料模型和Grüneisen状态方程为基础的热塑性流体力学模型,对同种钢板爆炸焊接界面波进行了数值模拟.结果表明,运用SPH方法可以得到清晰的爆炸焊接界面波形貌.与实验结果的比较表明:模拟误...     

单侧界面条件下弹性波传播规律的研究

界面是介质间的分界面,把界面看成是双侧的,只是真实界面模型表述的一个极端,在很多情况下,它的另一个极端,即单侧界面模型可能更近于实际．在单侧界面条件下弹性波的传播问题,与双侧界面问题相比,还是一个尚未充分研究的领域．本文意在对这一领域的研究现状及问题作一评述和展望      

多孔介质——自由流界面应力跳跃条件下流动特性解析解

对非对称多孔介质--自由流复合通道内多孔介质内部及多孔介质与自由流体界面处复杂质量、动量输运特性进行研究. 在多孔介质区采用Brinkman-extended Darcy模型并结合速度连续,剪切应力跳跃的界面条件对此复合通道内流体的传递现象进行求解,提出了考虑界面应力跳跃时非对称复合通道各区域流体运动速度及摩擦系数的解析式,分析了界面应力跳跃系数,达西数及无量纲多孔层偏心厚度对流体速度及摩擦系数的影响. 结果表明:改变界面性质可在一定条件下明显控制各区域流体速度分布;在达西数、多孔层偏心厚度一定情况下,界面应力系数的增大会使界面流速减小,而使流体摩擦系数增大,特别是界面应力系数小于0的情况下变化更明显,此时若不考虑界面应力系数则会造成较大误差. 当界面应力系数及多孔层偏心厚度均为较小负数值时,改变多孔层偏心厚度对界面速度的影响要大于改变界面应力系数的情况;而当界面应力系数及多孔层偏心厚度为较大正数值时,情况则相反. 较大达西数下,界面应力系数及多孔层偏心厚度对流体摩擦系数的影响均较大,继续减小达西数至一定程度时,界面应力系数对流体摩擦系数的影响可忽略不计而认为只与多孔层偏心厚度相关,且对较大多孔层偏心厚度更敏感.      

有限厚度流体层界面运动Rayleigh-Taylor不稳定性的数值模拟

采用二阶ＴＶＤ格式及ＬｅｖｅｌＳｅｔ方法计算了二维可压缩有限厚度流体层Ｒａｙｌｅｉｇｈ Ｔａｙｌｏｒ流体不稳定性。计算结果与Ｔａｙｌｏｒ的线性解和Ｏｔ的薄层非线性解析解符合很好。     

多流体界面不稳定性的守恒和非守恒高精度数值模拟方法

对多流体界面问题守恒和非守恒格式(M)WENO重构方法进行探讨,采用虚拟流动方法并用Level-Set函数捕捉界面的运动变化。数值模拟结果表明本文的数值方法具有较高的分辨率,并能有效地抑制界面附近的非物理振荡。     

不同节理角度和地应力条件下岩石双孔爆破的数值模拟

在采矿工程、地下交通工程和水利水电工程等岩体开挖工程中,爆破法仍是一种主要的破岩方法. 实际天然岩体中存在的裂隙、节理等不连续面和所处的地应力场环境,都会对爆破荷载传播过程和最终的破碎效果产生重要影响. 本文把岩石爆破视为爆炸应力波动态作用和爆生气体压力准静态作用两个独立的先后作用过程,同时考虑初始应力场的静态作用,建立岩石爆破的力学模型. 基于渗流方程描述爆生气体在爆生裂纹中的传播过程,进而基于流固耦合理论实现爆生气体准静态压力对裂纹尖端的力学作用. 爆炸应力波主要在孔周引起压碎区和径向微裂纹区,随后爆生气体压力楔入径向裂纹尖端,促使裂纹进一步扩展,形成径向主裂纹. 数值模型可以再现孔周压碎区、径向微裂纹区、径向主裂纹区萌生、扩展的完整演化过程. 针对不同节理角度和地应力条件下岩石双孔爆破过程的数值模拟结果表明,初始地应力场的压应力作用不利于爆生裂纹的萌生与扩展,但节理的存在对裂纹的扩展具有明显的导向和促进作用,有利于爆生裂纹沿节理面方向的扩展.      

微重条件下两底为半球面腰为圆柱面的贮箱中液体的晃动特性

本文研究了微重条件下两底为半球面腰为圆柱面的贮箱中液体的晃动特性。首先分析了静液面形状。然后用一类特征函数展开的方法求出液体自由振动频率、液体对腔壁的反作用力和力矩,在此基础上得到了等效力学模型.最后给出了数值计算结果。     

THE POROUS EFFECT ON THE BOTTOM INSTABILITY UNDER PARTIALLY STANDING WAVES

In this paper, we extended Blondeaux [1], Liu, et al. [2] and Mei and Yu‘s [3] the ories to study theporous effect on instability of rippled bed under the partially standing surface waves. Comparisons are made with both Blondeaux and Mei & Yu‘s theoretical predictions and experimental data performed by many other authors. We found that based on the linear instability analysis and the parameter regimes we explored the percolation effect can change the threshold conditions by 10-20 percent and the seepage force is negligible up to O(∈). Furthermore, the steady streamings with or without porous effect show the considerable difference in intensities at the same locations and it might change the sediment transport rate especially the suspended load.     

零重力条件下低温射流抑制大尺寸液氢储罐热分层的数值研究

液氢是一种常用的沸点低、易蒸发的空间低温推进剂.空间微重力环境中浮力对流被极大减弱,当推进剂储罐壁面存在局部漏热时,储罐内部气液两相流体系会出现环绕漏热源的热分层现象,引起局部过热沸腾,导致储罐内部压力急剧增大,危害系统结构安全.利用低温射流抑制储罐热分层现象是一种有效手段.低温流体通过设置在储罐内部的射流喷嘴与储罐内...     

微重力环境下Marangoni对流的有限元数值模拟

本文以悬浮区为背景研究液桥中气液交界面上由表面张力所驱动的流体对流。我们采用有限元方法对轴对称、定常运动方程,能量方程,自由面切向、法向应力平衡条件迭代求解,首次考虑了边界形状的影响,得到了自洽的流区位形和流场、温度场、表面压强分布。一般说来,流区自由面呈弯月形。结果表明,只要表面张力数 C_a<1,静态平衡的流区位形就是动态情况的良好近似。本文还分析了C_a数及G_r 数对流区位形的影响,得到了不同 M_a 数及散热条件下的温度场和流函数分布。     

多孔材料冲击特性的数值研究

使用物质点方法对多孔材料的冲击响应过程进行模拟研究.重点考察了孔隙度、冲击波强度和空隙平均大小对波后物质状态参量的影响,计算了局域旋度和散度随时间的演化及其涨落. 研究表明,在孔隙度较低时, 波后系统可以达到一个近似动力学平稳态,下游区域的空隙反射回的稀疏波使得物理量随时间轻微振荡.在孔隙度较大时, 波后物质不易被压实, 冲击波的衰减明显.在孔隙度一定的前提下, 更大的平均空隙导致波后系统温度的进一步提高.冲击波的衰减速率依赖于孔隙度、平均空隙大小和冲击波强度.局部湍流和体积耗散是非均匀材料中机械能转化为热能的重要机制;加载的冲击波越强, 材料的多孔特征越显著.压强、密度、温度、粒子速度趋于或远离平衡的节奏并不同步.      

A NUMERICAL STUDY OF CYLINDERS IN WAVES AND CURRENTS

The research on combination flow of planar oscillatory flow plus an in-line steady stream is of importance to the situation of structures in waves and current. The combination flow has not been studied extensively. There is still little information about the effect of current and about combined effect of current and waves on hydrodynamic loading of the structures. The present study investigates the combination flow around a circular cylinder using a vortex-based method incorporating vortex moving particles (discrete vortices) with a finite-difference scheme for the vorticity diffusion. The main attention is paid to the effects of a small current on in-line fluid forces and vortex patterns in the wake. Morison's equation and an equation with two drag terms are examined. The results show that the presence of a small current in an oscillatory flow can reduce the drag coefficient significantly. Morison's equation gives reasonably good predictions for the in-line forces for an oscillatory flow plus a small current. The current tends to bring the whole vortex wake downstream and tries to form the stable Karman asymmetrical form in the downstream wake. The present results show certain agreement with some previous experimental results.     

非周期波浪与直墙作用的非线性数值研究

基于时域高阶边界元方法,建立了完全非线性二维数值波浪水槽,对非周期波浪与直墙的相互作用问题进行了模拟和研究.自由表面满足完全非线性自由水面运动学和动力学边界条件,采用混合欧拉-拉格朗日方法追踪瞬时自由面流体质点,采用四阶Runge-Kutta法对下一时间步的波面和自由面速度势进行更新.采用加速度式法求解直墙表面速度势的时间导数,对瞬时物体湿表面上的水动力压强积分,得到作用在物体上的瞬时波浪力.首先,将全非线性与Serre-Green-Naghdi(SGN)模型的结果进行了对比分析,发现对于大幅值双入射波问题,仅满足弱色散关系的SGN模型大大低估了最大波浪爬高;其次,研究了双入射波与直墙的非线性作用问题,发现线性预报对波浪最大爬高有较大低估,而波浪的非线性成分不只导致了自由面爬高的异常增大,也引起了局部自由面的高频振荡,该物理过程中,直墙所受的波浪载荷,也展示出了与波浪爬高相似的非线性特性;最后,对波浪爬升和波浪力的时间历程进行了频谱分析,发现入射主频波的部分能量传递给了更高频的波浪成分,反映出该问题具有典型的非线性特性.     

NUMERICAL SIMULATION STUDY ON DYNAMIC RESPONSE UNDER RUDDER CONTROL

针对超声速机动飞行器舵面操纵动态响应问题,开展了相应的算法研究,建立了适用于模拟舵面与弹体相对运动及弹体气动力动态响应的非定常计算方法和局部动态网格变形技术. 以方形弹超声速绕流为例,验证了所提出的定常和非定常计算方法. 对方形弹舵面操纵的3种基本过程(脉冲、阶跃和谐波) 进行了数值模拟,获得了舵面操纵下导弹的俯仰运动对操纵的响应过程. 研究结果表明,快速的舵面操纵过程,甚至会导致导弹的俯仰运动失稳.     

电场作用下气泡上升行为特性的数值计算研究

利用电场控制气泡形态及运动,强化气液相间传热传质是电流体动力学的重要研究内容之一.然而目前多数研究集中在非电场下的气泡动力学上,对于电场下的气泡行为特性及电场的作用机制仍需开展深入研究.本研究对电场作用下单个气泡在流体中上升过程的动力学行为进行了数值模拟研究.在建立二维模型的基础上求解电场方程与Navier-Stokes方程,并采用水平集方法捕捉了上升气泡的位置及形状.模拟结果的准确性与有效性通过与前人实验和数值结果进行对比得到了验证.通过改变雷诺数、邦德数和电邦德数等不同参数研究了电场下液体黏度、表面张力和电场力对气泡运动变形的影响.计算结果表明,电场对气泡的动态特性有显著影响.非电场情况下液体黏度和表面张力较大时气泡基本维持球状,反之气泡发生变形并逐步达到稳定状态.此外,电场作用使气泡在初始上升阶段发生剧烈形变,随着不断上升,气泡形变程度不断减小,且气泡的上升速度和长径比均出现振荡.垂直电场使气泡的上升速度有较大的提高,且随着电邦德数的增大,难以达到相对稳定的状态.