绕圆柱荷电气固两相流动的数值解法

提出了一种求解绕圆柱荷电气固两相流动的数值解法，即用时间推进法求定常解，用有限差分方法离散控制方程，用人工粘性法处理固相间断面，计算结果表明，此数值解法能够处理强固相间断面，并且可用于含强间断面的悬浮体两相淳动的数值计算。     

微重力条件下气固两相界面耦合燃烧的数值模拟

在建立气固两相界面耦合燃烧计算模型的基础上，对微重力条件下火焰在固相可燃物表面的蔓延过程进行了三维非稳态数值模拟．考虑气相燃烧过程的同时，对固相可燃材料也用有限差分法研究其内部的传热和热解过程，给出了体现两相耦合作用的相界面处理方法．还研究了强迫流动对微重力火蔓延的影响．为了分析对比，进行了正常重力条件下的计算．数值模拟得了微重力燃烧的一些基本特点，结果令人满意．     

液体－固体粒子两相流流动特性的数值模拟及计算

运用液体－固体粒子两相流理论及计算流体动力学ＳＩＭＰＬＥ方法，对液－固两相流动特性进行了分析和研究，在此基础上建立了物理模型和数学模型，编制了计算程序，并以空心体在压力筒内产生破裂的现象为实例进行了数值模拟计算．结果表明，本文的理论研究方法、数值计算方法及计算程序是正确可靠的，不仅为液－固两相流动的研究提供了手段，还进一步为研究三相流动打下了良好的基础．     

微重力条件下气固两相界面耦合燃的数值模拟

在建立气固两相界面耦合燃烧计算模型的基础上，对微重力条件下火焰在固相可燃物表面的蔓延过程进行了三维非稳态数值模拟，考虑气相燃烧过程的同时，对固相可燃材料也用有限差分法研究其内的传热和热解过程，给出体现两相耦合作用的相界面处理方法，还研究了强迫流动对微重力火蔓延的影响，为了分析对比，进行了正常重力条件下的计算，数值模拟得了重力燃烧的一些基本特点，结果令人满意。     

复式断面明渠中污染物混合输移的数值模拟(Ⅰ)——复式断面二次流的模拟

采用L-Y雷诺应力模型，对水面和固壁边界条件进行有效的处理，结合混合有限分析法数值模拟了复式断面明渠二次流流动．和实测资料的对比验证了数值计算和边界处理技术的合理可靠性，并进一步分析复式断面明渠二次流运动的流动结构和特性．发现二次流主要产生在固壁拐角和近岸区域，其中最突出的位置是滩槽结合部，存在明显倾斜指向主槽自由水面的二次流；主流速度等值线被二次流强烈影响凸向二次流指向．     

激波卷扬附壁煤尘的湍流模型

本文建立了激波卷扬煤尘的气固湍流边界层方程,其中考虑了气固两相间的耦合效应和气固相湍流效应。模型认为颗粒上扬是激波后气固流动的剪切升力和湍流效应的结果。数值计算结果与文献实验结果较一致。     

垂直管道内液-固两相流动的压降和相分布特征

为得到垂直管道内液-固两相的流动特征,采用电阻层析成像方法对流动的相分布和含率等进行测量,并采用无量纲参数对流动的压降进行分析。研究表明,垂直管道内液-固两相流动中固相颗粒均呈非均匀分布,且浓度呈近似轴对称分布;不同颗粒粒径的液-固两相流动中的相间速度滑移程度不同;液-固两相流动中,粒径较小固相颗粒的掺入可有效降低流动的摩擦压降,具有湍流减阻的作用,而由于粒径较大固相颗粒间存在频繁碰撞,会增加流动压降;对于相间速度滑移较小的液-固两相流动,可采用均相流模型对流动的压降等参数进行精确计算,为准确预测液-固两相流动的压降,合理设计海洋资源生产及输送系统提供理论依据。     

耗散粒子动力学在流动数值模拟中的应用

对耗散粒子动力学方法在流动问题中的应用进行了理论及数值研究，给出了耗散流体粒子运动的控制方程组、边界条件、数值方法及随机数据统计方法等．提出用固壁粒子层结合流体粒子反弹运动的方法来处理固壁边界有滑移流动边界条件，有效地消除了目前耗散粒子动力学及分子动力学模型在固壁附近的粒子密度波动问题．利用该方法编程计算了2类流动问题：方腔驱动流动与泊肃叶流动，计算结果与现有的理论及数值结果比较吻合，验证了所采用的数学模型、计算方法在流动数值模拟中的可行性及潜在优势．     

浓密液固两相流动的数值研究与理论分析

在ＩＰＳＡＲ算法基础上,推导了适用于考虑浓度变化影响的浓密液固两相流动数值计算的ＤＩＰＳＡＲ算法,并考虑了颗粒湍流强度梯度传输对颗粒浓度分布的影响,对竖直上升管中浓密液固两相流动分别采用ＩＰＳＡＲ和ＤＩＰＳＡＲ算法进行了数值计算和理论分析,结果表明,ＤＩＰＳＡＲ算法能更有效地预测浓密液固两相流动．     

液体—固体粒子两相流流动特性的数值模拟及计算

运用液体－固体粒子两相流理论及计算流体动力学ＳＩＭＰＬＥ方法，对液－固两相流动特性进行了分析和研究，在此基础上建立了物理模型和数学模型，编制了计算程序，并以空心全在压力筒内产生破裂的现象为实例进行了数值模拟计算，结果表明，本文的理论研究方法，数值计算方法及计算程序是正确可靠的，不仅为液－固两相流动的研究提供了了手段，还进一步研究三相流动打下了良好的基础。     

基于嵌入式离散裂缝和扩展有限元的裂缝性页岩油藏流固耦合高效数值模拟方法

通过耦合嵌入式离散裂缝模型和扩展有限元建立了一种适用于裂缝性页岩油藏的流固耦合高效数值模拟方法。该方法创新性在于:通过对裂缝进行显式降维处理,可以准确刻画水力裂缝对渗流场和应力场的影响,同时可以模拟水力裂缝开度的动态变化过程以及水力裂缝内流体压力对裂缝变形的影响;在划分网格时仅需对基岩系统进行正交网格剖分,不必考虑裂缝的几何形态,渗流场和应力场均选取一套网格,使得网格划分时的复杂度大幅降低。对于渗流方程求解采取模拟有限差分方法,满足局部质量守恒并且可以有效处理全张量形式渗透率;对于岩石力学方程求解采取扩展有限元方法,可以准确处理裂缝两侧的位移间断性。通过算例分析,本文方法的正确性与优势得以验证。     

新型高分辨率格式及其在CFD的应用

 利用小波奇异分析构造了一种新型的高分辨率高效计算格式,借助小波奇异分析所找到的流场中奇异点区域和光滑区域,合理选择高分辨率、高精度WENO格式和高精度中心差分格式进行流场的离散求解;利用该格式对超声速二维前台阶流动问题、跨声速RAE 2822典型翼型的二维外流场、跨声速VKI-LS59涡轮叶栅的二维流场以及NASA Rotor 37跨声速压气机转子和NASA Rotor 67跨声速风扇转子的三维流场进行数值模拟,计算结果表明,采用这种新型格式可以用于内外流复杂流场计算。数值计算还发现,在采用相同网格的情况下,本文格式不仅具有通常WENO格式高分辨率的特点,而且计算效率较WENO格式高出3～5倍。     

SPH方法在两相流动问题中的典型应用

应用光滑粒子法(SPH)对具有间断面的不同介质的流动问题进行模拟分析.在现有的SPH方法的理论基础上进行修正,通过对密度近似方程进行正则化处理,在理想气体状态方程中引入Van Der Waals修正项以及对物质点运动速度进行修正,模拟了密度相差较大的两相流动问题,并且使用图解分析了两相交界面在模拟过程中的清晰程度.通过气(液)泡上升以及两相溃坝问题对SPH方法在两相流动问题模拟中的可行性进行了分析.最后应用修正后的SPH方法有效再现了两相流动的物理过程,并且能清晰分辨不同介质之间的交界面,从而说明SPH方法能有效对两相流动进行模拟.     

激波诱导圆柱运动数值模拟

该文采用基于含有运动边界自适应非结构网格点的显式有限体积方法求解Euler方程，提出了一种运动边界的数值处理方法，对激波诱导的圆柱运动进行了数值模拟。由于物体的运动，物体边界附近的网格将发生变形。对运动边界附近严重变形的网格采用局部网格重新生成的方法进行处理。采用HLLC方法来处理运动边界问题，计算并比较了气流流过静止圆柱和圆柱在静止流场中运动2种情形。结果表明，应用HLLC方法处理运动边界是可行的。     

波导不连续性的加速收敛算法研究

论文提出了利用新的基函数系列求解波导不连续问题的加速收敛算法。结合不连续处的金属楔近旁的电磁场边缘条件,构造了变形的三角基函数用以展开波导不连续处的电磁场,利用里兹算法构造一个线性方程组,求解此方程组可以获得波导不连续的特性参数。计算表明,该算法具有极好的数值收敛特性;与耦合积分技术（CIET）不同,对于一个不连续性,本算法的方程组的阶数与基函数展开所用的模式数相等,从而保持了较低的阶数,利用求解。     

聚合物熔体轴对称共挤出流动的有限元计算

建立数学模型，利用有限元法作数值解，确定界面，速度场及压力场，利用界面处法向速度为零，通过迭代求解界面位置，在界面处采用双结点法以保证速度和应力的连续性并满足压力的不连续性，使用ｅｎｎｅｒ迭代法求解圆形通道中两个共轴层流动的非线性方程。求出共挤出轴对称棒材的流量压降关系，计算结果与解析解比较，有很好的一致性。该数值方法可以运用于无法得到解析解的复杂流道。     

基于光流场模型的医学图像配准

基于最大化互信息的配准方法因其复杂性高而不能满足临床的实时性要求;基于光流场模型的图像配准方法计算简单快速,但已有的采用光流场进行配准的方法都直接采用原始的Horn模型,该模型简单地采用速度场的模作为光滑性约束,而导致配准过程不能保持图像的不连续性,使得图像严重模糊以致无法应用,该模型还存在不能处理大位移的问题.为了能在配准过程中尽可能地保持图像特征,采用了Lucas等人提出的窗口化光流场模型,通过选用小的窗口减轻对图像造成的模糊;对于大位移问题,可以采用金字塔结构来有效处理.实验结果证明,Lucas光流场模型应用于医学图像配准能够得到较准确的配准结果,具有很高的运行效率.     

模拟两相流动问题的修正扩展有限元方法

基于SUPG/PSPG稳定化方法,在处理两相流不连续问题时,引入修正扩展有限元方法使混合单元附加形函数满足单位分解原理,给出一种模拟两相流体流动问题的扩展有限元方法.在模拟两相流体流动过程中采用水平集方法捕捉流体运动界面.利用编制的计算程序模拟液体自由晃动问题,数值模拟结果与理论解一致.进一步采用本文方法和有限元方法模拟了溃坝流动问题,并将数值结果与实验结果进行比较,结果表明本文方法更有效.本文方法能够准确捕捉流体运动过程中两相界面的变化,且具有在计算过程中无需进行网格重构的优点.     

基于光流场模型的医学图像配准

基于最大化互信息的配准方法因其复杂性高而不能满足临床的实时性要求；基于光流场模型的图像配准方法计算简单快速，但已有的采用光流场进行配准的方法都直接采用原始的Horn模型，该模型简单地采用速度场的模作为光滑性约束，而导致配准过程不能保持图像的不连续性，使得图像严重模糊以致无法应用，该模型还存在不能处理大位移的问题．为了能在配准过程中尽可能地保持图像特征，采用了Lucas等人提出的窗口化光流场模型，通过选用小的窗口减轻对图像造成的模糊；对于大位移问题，可以采用金字塔结构来有效处理．实验结果证明，Lucas光流场模型应用于医学图像配准能够得到较准确的配准结果，具有很高的运行效率．     

低耗散TVD格式的构造

目的　为进一步提高双曲型方程数值解的精度,减少数值耗散,同时保证稳定性,处理有强间断的流动问题; 方法　从总变差减小（ＴＶＤ）定义出发,分析保证数值解ＴＶＤ的条件,找出满足研究目的的可能性和途径; 结果　首次提出和证明双曲型方程数值解ＴＶＤ的充分必要条件,并据此简化导出适用于时间显示格式的ＴＶＤ充分条件,并据此构造了一种低耗散ＴＶＤ 格式; 结论　进行的数值实验证明,高分辨格式可以进一步减小数值耗散,提高数值解精度;