固体火箭发动机燃烧室层流流场的数值模拟

该文发展了一种在任意曲线坐标系上求解层流Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的数值方法，该算法以ＳＩＭＰＬＥ为基础，采用了非交错网格，因而对原始算法中的压力修正方程进行了改进。用准定常方法数值模拟了固体火箭发动机燃烧室内的二维轴对称流场，计算结果能够反映流场内的旋涡与各参数的分布。计算表明，压力与速度等参数的分布明显受旋涡存在的影响，比传统的一维流场复杂得多。     

用一个改善的标志网格方法对二维方柱涡脱落流的计算（英文）

对处于均匀流中高Ｒｅ数的二维方柱的气动力行为进行了数值研究．使用了二维不可压粘性ＮａｖｉｅｒＳｔｏｋｅｓ 方程和三维不可压ＮａｖｉｅｒＳｔｏｋｅｓ 方程两个模型,没有采用湍流模型,在非均匀网格系统下进行了计算,对对流项采用三阶逆风格式,求解压力方程组使用多重网格,大大提高了传统的ＭＡＣ方法的效率．对均值和脉动气动力以及Ｓｔｒｏｕｈａｌ数、功率谱等项与实验进行了比较,取得了比较满意的结果     

Navier-Stokes方程的最佳非线性谱Galerkin算法

提出了求解二维非定常ＮａｖｉｅｒＳｔｏｋｅｓ方程的最佳非线性谱Ｇａｌｅｒｋｉｎ算法,分析了近似解的一致收敛速度．和标准的谱Ｇａｌｅｒｋｉｎ算法与非线性谱Ｇａｌｅｒｋｉｎ算法相比,该算法具有节省计算量的优点．     

二维方腔流动的有限元模拟

采用了Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法离散Ｎ－Ｓ方程，该法保留了高阶的空间离散精度，隐含了流线迎风的耗散作用，采用了压力校正法求解流场中各原始变量，推导了压力修正Ｐｏｉｓｓｏｎ方程的有限元离散格式，最后给出了二维不可压缩方腔流动计算结果。     

求解不可压缩流动的一种新算法

在Jameson有限体积法的基础上发展了一种模拟不可压流动的数值方法。该算法用四步Runge-Kutta时间推进法直接求解非线性的动量方程,充分利用了良好的非线性收敛特性,并从连续方程出发推导了求解压力的方程,解决了压力与速度的耦合问题。最后用该方法在同位网格上,对二维平板边界层流、二维顶盖驱动方腔层流及二维扩压器管道湍流进行了数值模拟。计算结果与理论值、标准解及试验数据进行了比较分析,从而证实了该方法对模拟不可压流动的适用性,并显示出了编程简单、易于推广、计算效率高等优点。     

不可压缩流动的TG有限元模拟

采用三阶Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ（ＴＧ）有限元法离散Ｎ－Ｓ方程，并对离散格式中由高阶时间导数转化而来的附加项进行了缩并运算修正，使其起到了人工粘性的作用。由于该法隐含了流线迎风的耗散作用，可较好地应用于较高Ｒｅｙｎｏｌｄｓ数流动的计算。同时采用压力校正法求解不可压缩流场的各原始变量，并讨论了压力校正Ｐｏｉｓｓｏｎ方程的边界条件；最后给出了二维方腔拖带流的计算结果，证明了算法的正确性。     

N—S方程的非线性Galerkin算法的稳定性

提出了求解二维非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（Ｎ－Ｓ）方程的３种全离散非线性Ｇａｌｅｒｋｉｎ算法，其中空间离散由谱元非线性Ｇａｌｅｒｋｉｎ算法进行，时间离散分别由Ｅｕ－ｌｅｒ法，修正的Ｇｒａｎｋ－Ｎｉｃｏｌｓｏｎ法和两步法进行．此外还分析了这些算法的有界性和稳定性，并给了时间步长的约束条件     

SIMPLE算法中压力修正方程边界条件确定新探

不可压缩粘性流体Ｎａｖｅｒ－Ｓｔｏｋｅｓ方程ＳＩＭＰＬＥ算法求解中压力修正方程起着很关键的作用。但其边界处理的特殊性常给正确应用带来的一定困难。提出从连续性方程本质出发确定压力修正方程边界条件的新方法。     

人工压缩法数值模拟非定常不可压缩热对流流场

利用人工压缩方法数值模拟了由于强热浮力不稳定性引起的非定常不可压缩热对流流场．在非定常动量、能量方程中引入了一伪时间项,连续性方程中加入与伪时间层相关的人工压缩项,使连续性方程与动量方程自动耦合,得到一沿物理时间准确推进的沿伪时间层采用ＧａｕｓｓＳｉｄｌｅ内迭代的推进迭代算法．在每一时间层进行内迭代直到连续性方程满足为止．对流项用四阶中心差分格式,时间项用三阶向后差分．利用此算法数值研究了二维ＰｏｉｓｅｕｉｌｅＢｅｎａｒｄ流场和底部加热的非定常水平射流问题,成功地模拟了流场中涡旋的生成及发展过程．结果表明,用该算法模拟非定常热对流问题具有较高的效率和准确性     

一维和二维溃坝波的混合有限分析解

针对水深平均二维浅水方程高分辨率模型,采用非交错网格下的二维对流扩散方程的混合有限分析法,运用SIMPLE算法结合Rhie and Chow的动量插值技术计算二维溃坝水流．一维理论解和其他格式的二维数值结果符合良好,表明非交错网格下的混合有限分析法能正确有效地求解溃坝问题．进而对部分溃决情况下的溃坝波传播和绕流问题进行了数值模拟,给出了精细的数值结果,揭示了复杂的运动特性,进一步表明混合有限分析法具有较强的溃坝波的捕捉能力,是求解溃坝流动的有效方法之一．     

圆柱绕流的非交错多层笛卡尔网格混合有限分析法

采用非交错多层笛卡尔网格离散计算区域，应用混合有限分析法(HFAM)离散控制方程，利用“壁法向插值法”处理不规则边界，以改进的动量插值法(MIM)解决压力锯齿波问题，立了一套完整的计算二维复杂边界流动的数值方法．通过低雷诺数圆柱绕流资料进行计算检验，计算结果与资料基本吻合，显示了本模型在不规则区域计算中的实用价值．     

四阶紧致格式有限体积法湍流大涡模拟

为准确预测不可压复杂湍流,提出了一种应用于大涡模拟的高精度有限体积法。该方法在非交错网格上数值求解大涡模拟方程,空间离散采用有限体四阶紧致格式,时间推进用四阶Runge-Kutta法,压力速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动和振荡平板上方的流动,证实了该方法具有近四阶精度;并在此基础上,采用动力Smagorinsky亚格子应力模式,成功地完成了充分发展槽道湍流的大涡模拟计算,所得结果与直接数值模拟结果吻合良好。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

基于混合有限分析插值的非交错网格方法

在求原始变量N—S方程的数值解时，如果采用常规方法建立网格，将速度和压力放在同一套网格节点上，计算中可能出现不合理的压力锯齿波现象．为解决这一问题，应用Rhie和Chow首先提出的动量插值思想，建立了一种基于混合有限分析插值的非交错网格方法，并通过计算实例对该方法进行了检验．检验结果表明，该方法具有良好的数值效应．     

三维抛物型流动问题的非错列网格计算方法

用控制容积有限差分法计算三维抛物型流动，通常要把三个速度分量分别定义在三套相互错列的网格上，以消除不皮形场对解的真空性所产生的影响。网格的错列使得计算过程和程度变得相当复杂，造成了应用上的困难。本文利用虚拟网格和插值技术，在网格不错列的的前提下，建立了三维抛物型流动的新的计算方法。对典型流动问题的计算结果与实验数据、近似分析解以及用错列网格方法计算的结果一致。这种方法对简化述化床悬浮空间内的气－…     

柱体绕流的自适应数值模拟方法

发展了基于四叉树数据结构的网格生成和流动的Navier-Stokes方程数值求解器。采用压力梯度或者密度梯度的绝对值作为网格自适应的控制参量,同时采用基于最小二乘法的无网格方法处理对于一般Cartesian网格难于处理的物面边界条件。采取了绕方柱流动和绕圆柱流动的经典二维定常和非定常层流算例对所发展的方法进行了验证。计算的结果验证了所发展的方法在处理绕流流动时的合理性和有效性。从而也为将来数值模拟具有较复杂几何外形的流动提供了一种网格布局合理、高效,边界处理简单易行的新思路。     

多段翼型襟翼滑动非定常粘性流数值模拟

运用动态嵌套网格技术和双时间推进算法,对多段翼型襟翼滑动的非定常粘性流进行数值模拟和研究.计算定常情况下的压力分布,并与实验数据进行对比,验证数值方法的实用性.同时对非定常情况进行数值研究,分析襟翼在不同的滑动速度下,幅度和频率对多段翼型的阻力、升力和力矩产生的影响.得到的结果对襟翼运动形式的设计具有一定的借鉴意义.     

求解不可压缩流动的分步有限元格式

提出求解不可压缩 Navier- Stokes方程的分步有限元格式 ,该格式没有高阶微分项产生、程序编制简单 ,适用于非线性的多维复杂流动。应用该方法实际模拟了二维圆柱绕流的旋涡形成与脱落过程 ,得出了不同 Re情况下圆柱绕流的流速分布。计算得到的不同 Re下的旋涡脱落频率(Strouhal数 )与前人已有的经典解答符合良好     

HLLC方法在二维非结构动网格上的应用

以自适应非结构网格的显式有限体积法为基础，采用格心方法以及基于近似黎曼解的Godunov一阶精度方法求解Euler方程，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，空间和时间都是一阶精度，对含有运动边界的多体干扰流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步进行多体分离数值模拟打下了基础。     

一种低雷诺数有回流区的流场的数值计算

本文利用改进的SIMPLER方法求解二维N-S原始方程,并以突扩管道在突扩附近的流场为例进行了数值模拟,计算了雷诺数分别为180和500的两种流动,与用其它方法计算的结果对比,表明本文所采取的方法是可行的。     

二维非定常有涡流动的数值模拟方法

为有效地模拟二维有分离现象的粘性流动,发展了求解涡量流函数方程的数值方法。对涡量输运方程的时间项运用四阶Runge-Kutta法离散,从而将方程分解为4个计算步分别求解。在每一计算步对方程进行Steger-Warming近似因式分解,从而使涡量的计算在空间的两个方向上分别进行。对流项采用Chakravaythy-Oscher总变差减少(TVD)格式离散。流函数的Poisson方程运用Tschebyscheff SLOR方法交替方向迭代求解。运用该方法对突然起动圆柱和高雷诺数时弯曲薄翼的非定常有涡流动进行了数值模拟,并将计算结果与其它方法的计算结果和试验结果进行了对比,结果表明本方法具有精度高、数值稳定性好和计算效率高的优点。