基于体积分数形式的多介质流动数值模拟

把非守恒形式的体积分数方程与Euler方程组耦合，用交错中心型格式求解流体力学方程组。对一维情况下两种气体的相互作用和二维气泡的变形及塌陷问题进行了数值模拟，并把模拟结果与γ-model方法和NND（Non—oscillatory and non—free—parameter dissipation difference scheme，无振荡无自由参数耗散）格式的模拟结果进行比较，达到了抑制界面两侧非物理振荡和保持体积分数正性的目的。     

运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

求解织物热湿耦合方程的控制体-时域递归展开算法

为求解织物热湿传递耦合方程,给出了一种基于控制体积法和时域递归展开的求解织物热湿传递耦合方程的算法。首先,在时域对方程变量和参数进行级数展开,然后使用控制体积法对方程空间域进行离散,从而化连续性的非线性微分方程组为一系列的递归形式的线性代数方程组。文中给出了问题的求解步骤和算例。结论表明该方法的预测结果具有不依赖时间步长的特点。     

求解织物热湿耦合方程的控制体一时域递归展开算法

为求解织物热湿传递耦合方程,给出了一种基于控制体积法和时域递归展开的求解织物热湿传递耦合方 程的算法.首先,在时域对方程变量和参数进行级数展开,然后使用控制体积法对方程空间域进行离散,从而化 连续性的非线性微分方程组为一系列的递归形式的线性代数方程组.文中给出了问题的求解步骤和算例.结论 表明该方法的预测结果具有不依赖时间步长的特点.     

力限振动试验技术理论基础研究

为了深入理解双控方程,由运载火箭有限元模型的多自由度子结构方法,导出缩聚在星/箭连接界面处的全箭动力学方程,定义相应参数后,导出星/箭耦合中的双控方程,给出界面加速度、有源机械系统的自由加速度、界面力、紧固力、源阻抗和负载阻抗表达式,从而加深了对星/箭耦合双控方程的理解,进而促进力限试验技术研究.同时,从多自由度模型简...     

基于流体速度势的非线性晃动问题的ALE有限元分析

基于流体速度势，本导出了非线性晃动问题的ALE有限元格式。通过引入ALE运动学描述，使网格的运动独立于流体的运动，通过采用与贮箱固连的非惯性参考系建立控制方程，使网格只须相对贮箱更新，因此运动边界和自由液面的跟踪算法不仅非常简单，而且不会引起网格畸变．本分别采用有限元法和有限差分法进行空间和时间离散，中的算例分别模拟了贮箱受横向激励和俯仰激励时流体的非线性晃动问题，所得结果证实了本方法的有效性与可靠性。中的方法还可进一步推广用于流体—结构耦合问题．     

有效载荷／运载火箭耦合分析的子结构模态综合法

为了制定有效载荷设计和验证所需要的低频振环境条件，需要通过耦合动力分析给出有效载荷与运载火箭连接界面的加速度响应。本文利用自由界面模态综合技术研究有效载荷／运载火箭的耦合分析问题，给出了连接界面分别为刚性连接和弹性连接民政部下计算有效载荷界面加速度响应的子结构模态综合方法的两种形式。     

改进的CFD/CSD耦合系统设计方法及其计算

论述了界面上的数据交换方法和时间推进格式在流体动力学（CFD）和结构动力学（CSD）耦合求解非线性气动弹性问题中的关键作用;改进了CFD／CSD耦合系统．包括基于边界元方法设计了CFD和CSD耦合界面的数据转换方法,该方法可在同一映射矩阵处理结构响应和非定常气动载荷转换,保证了耦合边界上的能量守恒;改进了一种松耦合方法流程,在保持模块化基础上提高了计算精度和效率。最后将本文方法应用于二维位移外插算例和A—GARD445．6机翼的动响应分析中,结果表明本文方法具有较高精度。     

一种充分局部化的子结构并行分析方法

机群系统是一种低成本的松散耦合型分布式并行平台,它要求并行算法设计时遵循"分而治之"的原则,尽量降低节点任务之间的相关性。本文在FETI方法和A-FETI方法的基础上,直接从力学概念出发,提出了一种充分局部化的FETI方法。该方法在进行子域界面处理时引入三重变量:界面节点位移、界面节点力、分区框架上耦合节点位移,由此得到一组近似解耦的界面方程,使得各个子域的计算相对于经典的FETI算法更加独立。对得到的界面方程采用预处理共轭投影梯度法(PCPG)并行求解,所采用的预处理算子为局部化的集中型Dirichlet算子。分别在自建的两套4节点PC机群上进行了两组算例的验算,结果表明,本文方法具有很好的计算精度和收敛速度。并行加速比达到3.76。     

二维非定常化学非平衡流动高精度数值模拟及其与实验的比较

本文成功地将二阶迎风 TVD 数值格式运用到弱电离、化学非平衡斜激波反射流场的数值计算中。气体力学方程与组元质量守恒方程组全耦合求解是目前非平衡流场数值计算中的新方向,文中详细给出了气体力学方程和组元质量守恒方程组全耦合求解过程,最后给出了三种不同情况斜激波反射流场数值结果并进行了详细讨论,本文所得结果与文献〔4〕提供的实验结果和文献〔1〕提供的数值结果相符合。     

多流体碰撞复杂界面的界面追踪法计算

为实现对复杂运动界面的精确描述,本文对多流体碰撞复杂界面问题进行了深入的研究.以欧拉方程的黎曼解为基础,采用MUSCL格式,在现有的FronTier软件包中加入了新的函数和变量,以达到完成复杂界面的计算.FronTier软件包是基于界面追踪方法编写的,可以完成界面不稳定性问题在简单界面情况下的计算.文中的方法成功实现了多流体碰撞复杂界面的计算,计算结果与实际碰撞过程相符.     

非确定结构系统区间分析的直接优化法

工程中的非确定性问题可以用区间分析、随机理论或模糊集理论进行求解。本文采用区间分析法来处理结构静力分析和设计中的不确定性问题。将结构系统中的不确定性参数用区间数来表示，用有限元法建立系统的控制方程。该控制方程是线性区间方程组。本文讨论了一些求解区间线性方程组的方法并提出了一种直接优化法。将方程组中的所有区间数都作为设计变量，区间量的变化区间作为相应的设计变量的边界约束，运用约束优化法求出方程组解的各元素的最大值和最小值。文中给出了两个算例，列出了本文算法与其他算法的结果比较     

二维多介质流动问题的界面处理方法

就二维可压缩多介质流动问题的数值模拟,给出了一种新的界面处理方法。通过在界面处构造Riemann问题,利用Riemann问题的解分别定义界面两边流体的边界条件,由于Riemann问题的解准确地描述了界面处流体的流动状态,因此得到了更加准确的界面边界条件。本文将由Riemann问题的解得到的界面速度外推到整个流场,重新定义速度场,避免了由于速度的大梯度变化而导致的Level—Set等值线相互交错,因而得到了更加精确的界面位置。利用该方法对水下激波与柱型气泡相互作用问题进行数值模拟,结果表明该方法能准确地捕捉各种物理现象。     

星/箭耦合动态响应分析新方法研究

介绍了采用约束边界模态展开表达式进行星/箭耦合响应分析的方法,导出缩聚在星/箭连接界面处的全箭动力学方程.应用这种方法通过解析推导可以将整个星/箭耦合结构多自由度方程减缩为仅含星/箭连接界面处自由度的方程,形成星/箭耦合响应分析新方法,内部自由度的影响化为固定界面模态有效质量的贡献,它们的贡献累加于界面自由度之上,使动...     

振动中的逆问题

本文考虑有限范围内无阻尼弹性系统无限小自由振动问题。述评由一些文献组成。这些文献涉及到从固有频率数据中唯一地重新构造一个振动系统。文献分为两组——一组涉及离散系统,该系统的逆问题与矩阵逆特征值问题紧密相关;另一组则涉及连续系统,该系统由Sturm—Liouville方程组中的某个方程所控制,或者由细梁弯曲振动的Euler—Bernoulli方程控制。     

毛细流动聚焦的实验方法及过程控制

作为一种基于毛细流动制备微纳尺度液滴的技术，毛细流动聚焦（Capillary flow focusing）在工程领域具有重要应用。回顾了基于吹气式和吸气式核心装置的毛细流动聚焦实验方法，展示了完整的测试平台，介绍了施加外部激励控制流动聚焦射流破碎的实验方法。同时，给出了同轴界面的拍摄方法，可获得清晰的内外层流体界面图像。在对流体锥形收缩阶段的研究中，探讨了几何参数与流动控制参数对流体锥形的形态与稳定性的影响。在稳定锥形下，研究了流动控制参数对液体射流直径、扰动波长及复合射流界面耦合的影响，并基于光的折射定律，对复合射流外层界面透镜效应所导致的内界面失真进行了修正。在对激励作用下射流破碎的研究中，考察了射流长度随振幅的变化，建立了尺度率关系，探讨了频率对生成液滴的单分散性及粒径的影响规律，为在实际应用中可控制备单分散性微液滴提供了理论与技术支持。     

圆柱贮箱类耦合系统的平动响应分析

本文以一种新方法建立了水平激励下圆柱贮箱类液固耦合系统的动力学方程。利用两类完全不同的Ｌａｇｒａｎｇｅ函数来描述液体子系统和结构子系统的运动，联立可得耦合系统的动力学方程。通过数值仿真发现，在一定的激励幅值和频率下，该耦合系统会出现零位漂移、分叉等非线性动力学现象。     

两层液体Benard—Marangoni对流的界面张力效应

采用PIV技术对两层流体Benard—Marangoni对流进行了实验研究,得到了不同厚度比下的速度场,分析了界面张力在各种对流模式中的作用。和Rayleig—Benard对流的三种临界对流模式相比,由于界面张力对上下层浮力对流的不同作用效果,使得更复杂的临界对流模式在Benard—Marangoni对流体系中出现;在实验中首次观测到三层涡胞结构等现象。     

气液两相分层流剪切应力的不确定度分析

基于水平气液两相分层流的压力梯度、波状液层和气壁剪切应力实验数据,对两流体剪切应力进行不确定度计算,分析各不确定度分量对液壁和气液界面剪切应力结果的影响.由剪切应力可以直接计算分层流动量平衡方程的本构参数,对本构参数的关联作了合理的解释.气液界面剪切应力受液层高度和压力梯度的准确性影响较大,而液壁剪切应力却相反,和气液界面剪切应力最大程度相关.通过和实验数据对比发现,液壁摩擦因子无法像气壁摩擦因子那样可以用单相管流的指数型关系式来描述,而是应该结合气液界面摩擦因子建立一个基于两相各流体特征参数的有效关联式.     

不同密度梯度的多层流体界面上的Richt myer-Meshkov不稳定性研究

在气-液界面Richtmyer-Meshkov(R-M)不稳定性研究的基础上,利用不同流体间密度的差异,构造了空气-硅油-水、空气-酒精-硅油两种流体粘度方向迥异的气-液-液三相界面,实验中以氮气作为高压驱动气体,在不同激波马赫数下对这两种界面R-M不稳定性后期尖钉与气泡区发展进行了测量与统计,通过对实验数据的分析,得出了相关规律,并用这些规律与已有的单层气-液界面R-M不稳定性研究成果作比较,得出了异同点.同时,还研究了两种三相界面在R-M不稳定发展中的差异.实验结果表明:当流体的粘度梯度方向(从小到大的方向)与激波方向一致时,界面失稳更加明显,湍流混合更为显著.