求解N-S方程的多边形网格与三角形网格对比

本文介绍了基于边的数据结果离散对流扩散方程的方法,编制了适用于任意多边形单元网格的有限容积法计算程序,利用方腔自然对流、外掠后台阶和移动顶盖驱等算例考察了三角形网格与多边形网格在数值计算方面的性能差别.结果表明,采用三角形网格收敛较慢且不容易获得网格无关的解,而采用与三角形网格对偶的多边形网格则收敛较快且易获得网格无关的解.     

非结构网格下涡轮级三维非定常N-S方程的数值解

本文在非结构网格下,采用时间上二阶精度、空间上不低于二阶的高分辨率格式求解涡轮级三维流场,得到了与实验数据大体接近的N-S方程数值解。计算表明:本文所构造的双时间步迭代格式具有稳定、高效的特征。实践显示:本文用FORTRAN与C语言混合编制的N-S方程源程序,可用于涡轮级三维非定常流场计算。     

利用拉格朗日第二类方程求正交曲线坐标系下的加速度

在正交曲线坐标系下计算加速度，因为涉及变矢量的求导问题，一般计算比较复杂．即使用张量方法计算也不简单．本文从第二类拉格朗日方程出发，推出了一个计算正交曲线坐标系下求加速度的公式．用此法计算便捷，出错的可能性减小．     

非结构混合网格上的NS方程求解方法

提出了一套较为通用的,完全自动化的非结构混合网格生成方法.在物面粘性作用区,采用一种改进的推进层方法生成三棱柱形和金字塔形网格;在其他流动区域采用阵面推进方法生成四面体网格.采用一种改进精度的格心有限体积法对三维NS方程进行了求解,在加速收敛措施方面,提出了一种新的当地时间步长取定方法来减小质量较差的网格单元对流场计算稳定性和收敛速度的不利影响.以M6机翼和DLR/F4翼身组合体外形的粘性流场作为数值算例,验证了上述网格生成和流场求解方法的正确性和实用性.     

二维不可压流函数N-S方程的多重网格方法

通过对二维不可压缩N-S方程的涡量-流函数方程组消去涡量而得到仅以流函数为求解变量的控制方程,从而 使不可压N-S方程的求解个数减到最少。求解方法采用本文提出的二阶精度的九节点紧致差分格式,因此无须对靠近边 界的网格点作特殊处理。为了加快迭代收敛速度,采用多重网格加速技术。数值实验结果验证了方法的精确性和可靠性。     

Euler及N-S方程的二维四边形/叉树混合网格自适应算法

描述了一种基于直角叉树网格的Euler和N-S方程自适应算法.由于考虑了粘性的作用,提出并使用了四边形叉树混合网格的方法,在几何表面附近生成贴体的四边形网格,外流场使用直角叉树网格.采用中心有限体积法,对Euler及N-S方程进行数值求解,对N-S方程的计算中加入了B-L代数湍流模型.在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离等特性的捕捉和分辨能力.采用上述方法,数值分析了单段和多段翼型的绕流问题.     

非交错网格下三维准稳态激光重熔熔池数值模拟

对三维维稳态激光重熔问题进行数值模拟,模型考虑了熔池内金属液体的流动和换热。在非交错网格下,采用修正的动量插值技术和ＳＩＭＰＬＥ算法求解动量方程。数值结果表明表面张力温度系数δγ／δＴ对熔池内金属液体的流动方式有重要影响。当δγ／δＴ为－３．５×１０＾－６时,熔池中环。当δγ／δＴ为＋３．５×１０＾－４时,与负的表面温度张力系数情形相反,熔池中只存在沿逆时针方向流动的主对流环。通过灰铸铁激光相变硬     

三维非匹配网格上的扩散方程数值求解

将子网格剖分的支撑算子方法,拓展应用于三维非匹配网格上的扩散方程求解.算例表明该方法在正交非匹配网格上能够精确获得线性解;在一般非匹配网格上可以达到二阶精度;在求解曲面网格和节点不共面网格时,精度比平面近似的方法要高,也可以达到2阶精度,同时也适合求解含有物质界面的混合介质网格.     

正交曲线坐标系中加速度的拉氏方程求法

采用拉格朗日方程给出一种求解正交曲线坐标系中加速度的简捷方法。     

复杂流道内N-S方程的非交错网格数值求解

本文试图用非交错网格求解二维曲线坐际系下的N-S方程,文中对压力修正方程进行了改造,导出了曲线速度分量的修正表达式,作为计算方法的检验,本文针对90°弯曲管道和波瓣喷管混合器三维流动进行了两个算例分析,结果表明该方法是稳定收敛的。     

求正交曲线坐标系中速度和加速度的解析方法

引入“广义速度”和“广义加速度”的概念，把正交曲线坐标系或的速度加速度简化为广义坐标及其微商的偏导。     

色散方程的一类具任意稳定性的显格式

本文对色散方程u_r=au_xxx构造了一类中间层包含四个结点,带两个参数m和θ的三层显式差分格式。当m和θ满足一定的关系时,其稳定性条件为|γ|≤(m+1)/(4(m-1))(|m|>1),从而当取m充分接近1时,可得到任意大的稳定性条件,并且保持截断误差阶不变。数值例子验证了理论分析的结果。     

用非正交曲线坐标系数值计算航空发动机加力燃烧室无化学反应的湍流流场

本文采用非正交曲线坐标系下非交错网格的SIMPLE方法对航空发动机加力燃烧室无化学反应的湍流流场进行了数值计算，湍流模型采用k-ε双方程模型．差分网格采用分区方法生成，计算时对整个流场进行分区迭代直至得到收敛结果．加力燃烧室湍流流场数值计算结果合理．     

弯曲时片压气机时栅内二次流的数值研究

本文应用Ｂｅａｍ－Ｗａｒｍｉｎｇ近似隐式因子分解格式以及ＭＭＬ代数湍流模型,采用拟压缩性方法求解雷诺平均拟压缩Ｎ－Ｓ方程组,对弯曲叶片压气机叶栅内的三维粘性流场进行了数值研究．结果表明,正弯曲叶栅内通道涡较直叶栅强,其诱导产生的壁面涡较弱且位于吸力面与壁端的角区内,与壁角涡接近．在正弯曲叶栅出口处,通道涡开始处于破裂前期,叶栅总损失增加。反弯曲叶栅通道涡较弱,其诱导的壁面涡很强,位于通道涡左上方,壁面涡和通道涡的有利作用,使通道涡更稳定,叶栅总损失比直叶栅和正弯曲叶栅要小．     

二维交错网格的GAUSS型格式

利用Gauss型求积公式在交错网格的情况下构造了一类不需解Riemann问题的求解二维双曲守恒律的二阶显式Gauss型差分格式,该格式在CFL条件限制下为MmB格式.并将格式推广到二维方程组,进行了数值试验.     

不规则形状热结构计算中的任意坐标变换和代数网格生成技术

航天器可重复使用防热结构通常需要采用数值模拟,进行外流场和结构热响应的耦合计算。本文针对前缘类热结构的差分计算中,关于不规则（不能用初等函数表示）区域的热响应计算,应用流场计算的思路,探讨了任意坐标变换和代数网格的生成。给出了有关控制方程、边界条件的处理方法,为进一步耦合计算提供了有关基础。     

复代数方程组几种迭代法的比较

本文从二维非线性Schrödinger方程出发,推导出五对角的复代数方程组,并应用高斯-赛德尔迭代法、SOR迭代法、复双共轭梯度法以及预处理复双共轭梯度法等对求解的计算量进行了比较。同时,又将复代数方程组化成七对角的实代数方程组,用高斯-赛德尔迭代法、SOR迭代法以及PCG法(预处理共轭梯度法)等进行了比较。结果表明,PCG法在上述几种方法中是最有效的。本文还对SOR松弛因子的选择进行了讨论。     

自旋为5/2的Bargmann-Wigner方程的严格解

在坐标表象中严格求解自旋为5/2的Bargmann-Wigner方程,导出了自旋为5/2的场(m≠0)的相对论性方程和动量表象波函数 关键词： 自旋为5/2的场(m≠0) Bargmann-Wigner方程 坐标表象 严格解     

多孔介质分形结构热方程的多重双尺度概率模型与数值模拟

多孔介质分形结构既含有大量微米量级甚至纳米量级的微孔隙,也含有大量厘米甚至米量级的宏观孔洞。如果又假设复相材料和孔隙的就位性及形状大小服从某一概率分布。此时无论是理论分析或是数值计算均非常困难。本文对此问题,初步建立多重双尺度概率模型,井给出数值实验结果。