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

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

中间梯度法视电阻率ρ_s/ρ_l的边界元法数值解

本文用边界元法的常数元解法与线性元解法,计算了复杂二维地电断面上中间梯度法视电阻率畸变值ρ_s/ρ_1,并相应地进行了物理模拟。该法是以场所满足的控制方程为基础,根据加权剩余法建立边界积分方程。然后,在电性连续体的边界上划分单元,进行数值求解的一种新的计算方法。数值计算与物理模拟结果对比相一致表明:边界元法是研究电法勘探复杂场问题的一种经济有效的数值计算方法。边界元法具有可使所求解的电法位场问题的空间维数降低,输入数据少、计算速度快、工作效率高等优点。     

地下洞室断面形状的优化设计

基于初参数数值解法计算隧洞衬砌的变位和内力可以得到很好的解答,而复形法是优化设计中的常用方法,将这两种方法合二为一,给定内力控制值,对隧洞的断面形状进行优化设计的初步探讨,在工程具有一定的实际意义。     

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

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

焊接熔池浮力－表面张力驱动流数值解的源项与边值处理

从固／液相变问题的单相统一模型控制方程组出发，用焓法‘模糊，处理固／液运动界面和相变潜热，用Ｔａｙｌｏｒ级数展开法处理非线性源项，并将附加源项法扩展应用于处理动量方程的表面张力边界条件，使整个控制方程组的迭代求解局限于内节点上，从而加快了方程组耦合求解的收敛速度，二维轴对称定点熔化焊浮力－表面张力联合驱动流的数值计算结果表明，这种源项与边值的处理方法较非线性源项的显式处理和热边界条件的附加源项处理计算速度提高１５％～２０％，且计算结果相当一致。     

抛物型问题的边界元重叠型区域分解法

边界元法是一种求解偏微分方程数值的计算方法，用边界元法来求解抛物型方程，如采用与时间有关的基本解，较其它方法可以采用较长的时间步长，从而节省计算时间，且计算结果精度高。区域分解法是把计算区域分解成若干子区域来分别求解，由于它将原问题分解，由大化小，由复杂化简单，并且可以并行计算，优越性是显而易见的。将这两种方法结合起来（边界元重叠型区域分解法）求解抛物型方程，利用区域分解法将求解区域划分为两个小的子区域，然后在子区域上用边界元法并行求解方程。数值算例表明边界元重叠型区域分解法行之有效的，数值试验显示这种方法的收敛速度依赖于子区域重叠面积。     

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

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

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

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

直齿外插齿刀变位系数计算机解法

改变传统插齿刀变位系数的解法，利用建立起来的直齿外插齿刀变位系数计算机解法的数学公式，在计算机上用数值法实施了插齿刀变位系数计算机解泊。根据计算运行的结果，进行插齿刀变位系数图表法校验，获得了满意的结果，也大大提高了其计算精度和计算速度。     

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

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

地下水非稳定流计算的有限元方法(Ⅰ)

本文研究地下水资源评价及水位预报的有限元解法,文章分两部分。第一部分研究非稳定流计算的有限元方法,提出了处理井点的“奇点磨光法”,数值例子说明了这种方法的有效性;第二部分将研究地下水计算的反问题,着重讨论处理不适定问题的惩罚函数法。     

水库调洪数值—解析解法

在水库调洪数值解法的基础上，提出复杂情况下调洪数值－解妥法。两种解法优于且可以取代水库调洪的水文传统方法。从理论与实用上满足了水库实时洪水调度计算所需的高速度和高精度的双重要求。     

坝体-库水动力相互作用的数值求解格式

针对目前数值格式求解效率的不足,基于广义Newmark-β法,建立了求解坝体和库水流固耦合动力相互作用问题的2种新的数值求解格式:交错迭代法和直接积分整体求解法.导出了有限元求解公式,讨论了2种数值离散格式的特点,并结合实际工程计算比较了2种求解方法的计算效率.结果表明:2种方法均能有效地求解坝体-库水动力相互作用问题;在相同计算条件下,迭代法的计算规模要小,但花费的总机时长;整体求解法的计算效率更具优势.     

低温贮罐支座的热解析

利用数值解法中较新型的边界元法来计算低温贮罐支座传热问题，它的主要基本原理是通过加权残余法或构造Ｇｒｅｅｎ函数，有助于两点函数表示的基本解及利用分部积分法，将区域上的偏微分方程转换成区域的边界上的积分方程，通过对区域边界离散化对该积分方程进行数值求解。本文对计算误差进行估计，并以计算新型ＬＰＧ船型上的贮罐支座温度场分布为例，对计算结果进行了分析和讨论。     

统计线性化方法及其计算机数值解

本文研究了用统计线性化方法计算一类特殊的多自由度非线性系统在随机激励下的响应问题,并讨论了这种方法的数值解法及在计算机上如何实现这种数值解法。本文以一个用泡沫塑料作减振元件的导弹包装箱为例,进行了减振性能计算,并用Monte—Carlo仿真法的计算结果作对比,说明随机激励的非线性系统用统计线性化求解的精度是相当高的,数值解法是可靠的。 本文的研究结果有助于解决现代工程中用非线性元件作减振器时存在的设计计算困难。     

超导网络临界温度的计算方法研究

详细介绍了用解析解法和数值解法（ 即矩阵方程的多项式解法） 分别求出了基元回路为正方形的“田”字形超导网络在外磁场中的临界温度，结果表明两种计算方法是完全等价的。并对多节点规则网络所采用的矩阵方程本征值解法的意义进行了讨论。     

非平稳随机过程非线性结合的首次超越可靠度

针对非平稳随机过程非线性组合的首次超越可靠度问题，提出了一种数值解法。该法将所考察时段离散化，以数值计算手段克服复杂的数学处理上的困难，所得公式可以方便地退化为平稳随机过程甚至随机变量非线性组合的问题，因而具有普遍的实用价值。     

非平稳随机过程非线性组合的首次超越可靠度

针对非平稳随机过程非线性组合的首次超越可靠度问题，提出了一种数值解法．该法将所考察时段离散化，以数值计算手段克服复杂的数学处理上的困难．所得公式可以方便地退化为平稳随机过程甚至随机变量非线性组合的问题，因而具有普遍的实用价值．     

局部均布荷载作用下四边支承矩形板的内力计算

以矩形板的Navier解为基础，采用带补充项的傅里叶级数作为挠度函数，研究了局部均布荷载作用下四边支承矩形薄板的弯曲问题. 推导了确定待定系数的线性代数方程组，给出了简支边和固支边不同组合条件下的统一计算公式. 讨论了带补充项法级数解的收敛速度，并与叠加法级数解及有限元数值解分别进行了精度和计算量的对比. 结果表明，带补充项法的级数解达到收敛的级数项数约为40项. 带补充项法的级数解与叠加法级数解具有同样的求解精度. 有限元解随网格的细分，计算结果逐渐接近级数法解. 级数解法的计算量与有限元解法相比是微不足道的. 研究成果适于进行构筑物顶板受局部均布荷载作用的结构计算.     

风口布置对空调车室内流动与传热影响的数值研究

采用k-ε紊流模型及贴体坐标，应用整体求解法计算了空调车室内气固耦合传热问题.用Monte Carlo法分析计算了太阳透射辐射在车室内各固体表面引起的附加热流变化，并以此作为能量方程的附加源项，对不同风口布置的空调客车室内流动与传热进行了数值模拟研究．结果表明，送回风口的布置对空调车室内空调效果有较大影响．