三维非稳定渗流自由面边界积分项的精确数值计算

利用固定网格法分析三维非稳定渗流问题时,将要面对两项积分难题:以自由面及单元表面为边界的空间积分及以自由面为边界的曲面积分。针对常用的任意8结点6平面三维普通单元,提出采用坐标变换及等参变换技术求取空间积分项的精确数值解;至于曲面积分项,建议改用单元非饱和区部分表面作为积分边界,经过坐标变换及等参变换处理积分边界后,利用高斯数值积分可求出曲面积分项的精确数值解。通过一个普通单元及一项均质半无限边界堤坝的实例分析,表明此方法的精确性和稳定性良好。     

基于神光Ⅲ原型的激光直接驱动自由面准等熵实验研究

长脉冲激光打靶的直接驱动方式效率较高,但对激光强度变化非常敏感。利用测量自由面速度历史的方法对激光直接驱动准等熵压缩实验技术进行了分析,为改进激光直接驱动的准等熵压缩实验技术提供技术基础。介绍了国内在神光III原型装置上首次开展的激光直接驱动自由面准等熵压缩实验。对实验靶型、激光波形、典型实验结果以及实验中的关键技术进行了分析。在激光能量为1000 J的实验中,使用成像型速度干涉仪获得了较为理想的准等熵压缩条纹图。实验发现,自由面速度达到11.3 km/s时,自由面反射的探针光信号消失。通过对比Multi-1D软件的理论模拟数据与实验处理的结果,发现由于Al自由面缺少约束,这类实验会造成波系复杂化的问题,但是对主冲击的研究还是有意义的。对某些特殊材料,在很难找到阻抗匹配窗口的条件下,利用激光直接驱动方式研究材料高压响应特性提供了另一条可实施的技术途径。     

“聚龙一号”装置上磁驱动铝飞片实验的数值模拟

采用二维磁驱动数值模拟程序MDSC2,对大电流脉冲装置"聚龙一号"上的151发次磁驱动370μm厚铝飞片实验、164发次磁驱动330μm厚铝飞片实验进行了数值模拟和分析。数值模拟表明:370μm厚铝飞片和330μm厚铝飞片的磁驱动过程中,VISAR测量的速度不是飞片自由面的速度,而是飞片中邻近自由面最近的固体反射面的速度。这是由于磁驱动飞片发射过程中,飞片自由面部分被烧蚀,密度低于固体密度状态,而飞片自由面和加载面中间的飞片还保持固体密度状态。VISAR测量的激光将穿过自由面的低于固体密度状态的飞片部分,到达飞片自由面最近的固体密度位置再反射回去,获得这一位置的飞片速度。数值模拟的飞片固体反射面速度历史和VISAR测量的速度历史相吻合。     

基于移动粒子半隐式法的表面张力模拟

采用移动粒子半隐式法(MPS)模拟了受表面张力影响的自由面流动。表面张力的计算采取了一种较适合于MPS方法的表面自由能模型。方形液滴振荡和射流断裂的模拟结果分别与理论分析和试验结果一致,同时进行了三维射流注水模拟,从而验证了MPS方法结合该表面张力模型可以有效、方便地进行自由面流动中表面张力现象的模拟。     

无单元法在有自由面渗流计算中的应用

针对有自由面渗流分析中的有限元固定网各法存在的不足，利用无单元法中积分网格和结点相互独立的优点，提出了有自由面渗流的无单元法。计算结果表明，无单元法可以方便地解决迭代计算中的自由面变化问题，实现了真正意义上的网格固定。     

采用等效渗透系数法搜索渗流自由面

采用一个具体的实例对比分析了当前最为常用的几种自由面确定方法 ,阐明了各自的优缺点 ;然后在此基础上对单元传导矩阵调整法作了简化和改进 ,提出了等效渗透系数法。应用表明 ,等效渗透系数法不仅具有足够的计算精度 ,而且还具备算法简单、易推广等优点     

实践十号卫星项目——热毛细对流振荡特征的地面研究

上部开口环形液池在水平径向温度梯度作用下会出现内部温度和自由面的振荡,本文研究了二者发生的临界条件.环形液池内柱加热外壁制冷,以0.5℃/min的速率线性升温得到水平径向温差,T型热电偶测量液层内部单点温度,高精度激光位移传感器测量液层自由面某点形变.随两端温差增加,当超过某一临界温度△T_(cr)时,开始出现振荡.实验结果表明,对同一种硅油,两种振荡的临界条件随液层厚度具有相同的变化趋势.对不同普朗特(Prandtl)数(Pr)的硅油,振荡临界条件临界马兰哥尼数(Ma_(cr))随着邦德数(Bo)的增加而变大.本文工作是中国科学院科学先导专项SJ-10返回式科学实验卫星项目-热毛细对流表面波空间实验研究的地面研究结果,该工作为空间实验提供前期的基础科学研究数据和实验保障.     

关于沿内角的自发毛细流动现象回顾

文章从静力和动力学的角度简要回顾了关于沿内角的自发毛细流动研究的最近进展.作为一个通用几何形状,内角在地面微观尺度下或处于失重状态的航天飞行器系统内大尺度下为液体提供有效的输运通道.当一定的几何条件得到满足并且当毛细力远远大于体力比如重力的时候,沿着内角会发生自发毛细力驱动流动现象.从静力学的角度来说,本文讨论的自发毛细驱动流动和当特定的边界条件发生突然变化,比如重力作用突然消失时带有内角的容器内部单值有限高度的平衡自由面的非存在性有关系. Concus-Finn方法可以用来确定这样的平衡自由面在一个横截面处处一致的柱形容器内的非存在性.用这个方法可以推导出在失重状态下一个内角为2α的通常柱形容器里,当接触角小于π/2-α时,平衡曲面不存在.通常来说，沿内角的自发毛细驱动流动属于层流.利用尺度分析和摄动法，成功分析了该流动的动力学特性，并且推导出对设计有用的封闭形式的解析解.一个典型的结果是在黏性流的范畴里毛细面端点的移动和t1/2成正比      

格子Boltzmann方法模拟双液滴撞击液膜的流动过程

采用单相自由面格子Boltzmann方法中流场计算,自由面条件和表面张力处理方法,研究双液滴撞击液膜的流动过程,实现对垂直间隔的两液滴相继撞击液膜流动过程的模拟.数值结果表明,液滴相继冲击液膜时,液滴间上下距离对液膜状态有很大的影响,液滴间距不同,液膜形状会出现很大差异.数值结果与实验结论定性一致.     

三维波浪-流-串行桩柱相互作用的数值分析

为探讨波浪-流-串行桩柱的相互作用,建立了三维数值模型,采用大涡模拟(LES)方法对湍流进行模拟。研究了波浪-流-串行桩柱相互作用下的拖曳系数,对波浪-流靠近结构时的非线性相互作用引起的拖曳力和涡的变化,以及在结构前后自由面的变化进行了分析,并且与波浪-流与孤立桩柱相互作用的研究进行比较。结果表明:由于桩柱的存在,在波浪-流的传播过程中出现了时间延迟,这就导致了桩柱前后附近自由面水位存在较大的位移差;由于桩柱之间互相干扰,下游桩柱的存在导致上游桩柱附近不能形成非对称的周期性涡脱落,上游桩柱横向拖曳力受到影响最大,顺方的拖曳力也受到明显影响;下游桩柱横向和顺向的拖曳力都受到较大的影响;同时涡的脱落会把其附近的自由面吸入。