微重力下密闭腔体内火焰传播过程的数值模拟

建立了微重力条件下火焰沿固体可燃物表面蔓延的数学模型．在考虑固相传热和热解过程的基础上，研究气相的燃烧动力学过程，并且给出了体现微重力燃烧两相耦合作用的相界面处理方法．首次用数值模拟的方法研究了密闭腔体内火焰传播的三维非稳态过程，对于微重力大小分别为１０－４ｇ和１０－２ｇ时的计算结果进行了分析对比．数值模拟的结果基本合理．     

二元工质Marangoni对流的实验研究与数值模拟

用实验和数值模拟方法研究了正值表面张力温度系数的二元醇类特殊水溶液Marangoni对流流动. 首先通过实验测量确定正戊醇水溶液表面张力在特定的浓度分布和温度区间内具有明显随温度升高而增加的变化区域, 然后利用特定浓度配比的正戊醇水溶液,采用PIV方法实验观测了矩形液池中二元工质液层在水平温差驱动下的Marangoni对流, 发现了不同于常规的反向热毛细对流流动, 测量的表面速度分布与相同工况的数值模拟结果进行了比较, 发现二者变化趋势一致. 实验观测和理论结果的比较进一步验证了表面张力温度系数为正值时二元工质液层的热毛细流动输运特性.      

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

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

多层流体的Marangoni对流

介绍微重力环境下矩形液池中多层不相混液体的Marangoni对流及其不稳定性.Marangoni对流流动的形成是由于在该系统施加一个与液体交界面相垂直的外加温度梯度.Marangoni对流的线性不稳定性分析是基于无限延长矩形液池内的两层流体系统.应用数值模拟研究了Marangoni对流的流体动力学和热传输特性,将多层流体的Marangoni对流同两层流体中外加温度梯度与流体交界面平行时引起的热毛细对流的主要特征进行了比较.     

圆筒形贮腔中微重力液体非线性晃动的数值模拟

本文讨论低重力液体在圆筒形贮腔中的非线性晃动问题。将ALE（任意的拉格朗日－欧拉）运动学描述关系引入到Navier-Stokes方程中,在时间域上采用分步离散方法中的速度修正格式,利用Galerkinlk加权余量方法推导了系统的有限元数值离散方程;推导了考虑表面张力效应时有限元边界条件的弱积分形式。推导了自由液面上法向矢量的计算公式。模拟了圆筒形贮腔中低重力液体的非线性晃动,并得到了自由液面、波高变化、压力响应等非线性动力特性。揭示了微重力液体非线性晃动的重要特征并将所得结论与现有的实验结果进行了比较。从而证实了本文方法的有效性与正确性。     

板料三维成形过程的有限元数值模拟

基于有限变形理论建立了三维金属板料成形过程的刚塑性有限元数学模拟,该数学模型采用物质坐标系中的Ｕｐｄａｔｅ－Ｌａｇｒａｎｇｅ描述,等向强化假设,考虑了板料的厚向异性,对于金属板料与模具有摩擦采用近似的库仑摩擦定律以改善计算的收敛性。为简经计算采用薄膜单元,并根据此模型编制程序。     

板材成形有限元数值模拟的研究进展

本文以不同时期引起人们广泛注意的标准考题１）为主要线索，对板材成形有限元数值模拟的发展作了简要的综述．并讨论了研究中仍存在的问题和研究的发展方向．     

半浮区液桥热毛细振荡流

采用非定常、三维直接数值模拟方法研究大Pr数半浮区液桥热毛细对流从定常流向振荡流的过渡过程．文中详细描述了热毛细振荡流的起振和振荡特征,给出了液桥横截面上振荡流的流场和温度分布．在地面引力场条件下计算的结果与地面实验的结果进行比较,得出液桥水平截面上的流场和温度分布图样以一定的速度旋转,自由表面固定点处流体的环向流速正、负交替变化的一致结论．     

液封液桥内热毛细对流的数值模拟

建立了液封液桥(不相溶混的双层同轴液柱)内热毛细对流的物理数学模型,采用涡量-流 函数法、利用有限差分格式对微重力条件下液封液桥内热毛细对流进行了数值模拟,得到了 双层液柱主流区的温度场和流场,证实了液封技术能削弱液桥主流区的热毛细对流,从而提 高浮区晶体生长质量,找出了液封厚度以及内、外层流体物性参数比对液桥内热毛细对流的 影响规律.     

浮区Marangoni对流的浓度边界条件

对于浮区热和浓度毛细对流问题,溶化界面处可以有给定浓度或给定浓度通量两类边界条件,它们给出的液桥内浓度分布有很大差异。本文讨论了这两类边界条件所对应的物理过程,以及其物理含意。     

半浮区热毛细对流的不稳定性与转捩

近二十年来,微重力流体开展了半浮区液桥热毛细对流的不稳定性与转捩的研究．文中给出了热毛细振荡对流发生的临界参数,分析了液桥几何位形（尺度比,体积比）、物理参数及传热参数对临界Ｍａｘａｎｇｏｎｉ的影响．报导了有关的地面模拟实验,微重力实验以及本问题的线性稳定性分析、能量分析和数值模拟结果,并介绍了定常轴对称热毛细对流通过非定常振荡热毛细对流到湍流的转捩过程和三种热毛细振荡对流的产生机理．      

半浮区液桥振荡对流的微重力实验

利用日本微重力中心８００ｍ落井装置，完成了半浮区液桥振荡对流的微重力实验，对振荡对流的典型物理量诸如内部温度、流场、自由面边缘变化及表面波进行了综合测量。实验结果给出了振荡对流由地球重力环境向微重力环境的过渡，以及不同几何参数半浮区液桥的振荡特征，并首次获得了微重力环境下热毛细对流的表面波位形及边缘振荡特征．     

半浮区液桥热毛细对流速度场特征

液桥流场的速度分布是液桥表面张力梯度驱动对流的一个重要研究内容。通过地面模拟微重力环境的液桥实验以及数值计算,研究了外加温差、液桥腰径等对流场速度分布的影响,实验结果与数值计算结果相当符合。结果还表明,当液桥腰径减小时,流场的结构发生了变化,涡心由一个转化成两个。     

固体颗粒在热对流下沉降的直接数值模拟研究

在任意拉格朗日-欧拉算法模拟等温惰性颗粒两相流的基础上,增加对能量方程的联立求解,在热对流条件下对固体颗粒在不同雷诺数下的沉降规律进行了直接数值模拟.结果表明:热对流引起了流场流动的变化和不对称,颗粒在热流体中沉降,热对流产生的力加速了冷颗粒的运动,尾部形成了涡脱落;颗粒在冷流体中沉降,热对流产生的力阻碍了冷颗粒的运动,尾部形成了羽流. 随着雷诺数的增大,颗粒经历了稳定沉降、周期性摆动到不稳定、无规则运动3个阶段.      

热毛细对流的液桥几何位形及浮力效应

本文讨论重力对不同高度、直径比液桥的热毛细对流的影响。当液桥高度、直径比增大时,液桥中的等流函数线呈双涡结构,这种流动图样并不必然与热毛细振荡流相联系。在地面热毛细对流实验中模拟空间微重力情况,液桥高度需小于1.5mm。在微重力环境中,液桥内的流场和温度分布介于地面相同参数液桥的上部加热和下部加热两种结果之间。因此,可以用地面实验结果估计空间液桥的对流和热输运情况。     

黄河龙口水利枢纽库坝区渗漏三维有限元分析

应用三维等参数有限单元法,对黄河龙口水库库坝区周围一定范围的初始渗流场、万家寨蓄水后的渗流场及龙口工程蓄水后的渗流场进行了数值模拟,并计算了相应的渗流量。