结构网格下求解非规则区域流场的一种差分法

对于粘性不可压缩流场的计算，提出一种标识边界的差分方法，对不规则区域的流场边界做标记，应用SIMPLE算法在直角坐标系的结构网格下离散求解N-s方程；对圆柱绕流，不规则腔体流等问题进行数值模拟，取得了较好的效果。     

不可压缩流动的无网格Galerkin方法模拟

Galerkin型无网格方法在求解不可压缩流动问题时,会遇到对流占优、速度-压力失耦等问题,本文基于CBS有限元方法的基本思想,提出了无网格CBS方法来解决上述问题.通过对平面Poiseuille流动的计算表明:无网格CBS方法在采用压力速度等线性基近似的情况下,当时间步长大于某个临界值时可很好地解决速度-压力失耦问题,且具有相当高的计算精度.     

基于嵌入式迭代的高精度浸入边界法

针对传统浸入边界法施加边界条件精度不高的情况,将附加体积力的求解内嵌到压强泊松方程的迭代当中,提出了基于嵌入式迭代的浸入边界法.将所提出的算法应用到圆柱绕流数值模拟中,并与他人的实验值和数值解进行了对比,验证了所提算法具有较高的计算精度.分析了插值函数是导致浸入边界法高估圆柱阻力系数的原因,并对阻力系数进行了修正,修正值与采用贴体网格得到的结果非常吻合.比较了本文算法和传统算法的圆柱迎流侧驻点附近的速度矢量和圆柱周围的流线,证实了本文算法在施加边界条件方面的优越性.比较了两种算法的计算效率,迭代计算附加体积力所导致的额外耗时不显著.     

基函数法在黏性可压缩流动中的研究与应用

研究基函数法在二维和三维轴对称情形下可压缩黏性流动中的应用, 并利用此法数值地计算了二维圆柱的超音速黏性绕流和三维球头的轴对称超音速黏性绕流两个算例。取一阶三角函数为基函数, 构造出导数的中心格式和迎风格式。对于可压缩N-S 方程中的对流项, 采用通量分裂法及中心格式与迎风格式相结合的技术;对于可压缩N-S方程中的黏性项, 则采用中心格式进行处理。由此, 构造出了数值求解黏性可压缩流动的一阶三角函数类型的基函数格式。两个算例的数值计算结果表明, 基函数法不仅在处理无黏可压缩流动时, 是一种高精度、高分辨率的新型计算方法, 而且在处理黏性可压缩流动问题时也一样行之有效。     

低速湍流中气动噪声预测方法

 为了研究气动噪声的产生机理和传播过程,在考虑介质黏性的影响情况下,采用分解法并结合湍流模型,对低速湍流流动中气动噪声问题在时域上进行数值计算。基于非结构化同位网格和有限体积法,把可压缩N-S方程分解成不可压N-S方程、含黏性项的声扰动方程;为了考虑质点振速和声压耦合,采用SimpleC算法来同步求解不可压N-S方程和声扰动方程。进出口远场边界采用以渐近解为基础的无反射边界条件,并采用与内部区域相对应的有限体积法、时间隐式格式对其进行求解。利用所编制的程序进行层流状态下圆柱绕流气动噪声仿真验证,并与文献结果进行对比,检验本方法的正确性;并结合湍流模型将数值解法推广到湍流状态下气动噪声数值模拟中。结果表明该方法能够很好地反映流场和声场的形态,无反射边界能很好地抑制声波在边界处的反射,适合低速流气动噪声问题模拟,为实际工程中的降噪工作提供预测信息。     

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

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

基于四叉树切削网格的N-S方程求解方法

以四叉树非结构化网格为基础,提出了背景区域采用正方形四叉树网格、边界区域采用切削网格的一种可以表达复杂几何形状的网格生成方法,该网格具有生成过程简单,正交性好等优点.在这种网格的基础上,采用非结构网格有限体积法进行离散得到了多种形状切削网格并存时Navier-Stokes（N-S）方程的求解算法,并以顶盖驱动斜方腔流和方腔内热圆柱自然对流为例,应用上述算法实现了网格生成和流动数值模拟,与基准解进行了比较,一致性较好.计算结果表明这种切削网格方法及其N-S方程求解方法具有可靠性和应用前景.     

应用滑移界面方法求解动/静相干叶排内非定常流动

基于滑移界面技术,发展了求解二维动／静相干叶排内非定常粘性流动的计算方法和软件,该方法克服了混合平面方法所引入的方法误差．压力与速度耦合采用非迭代ＰＩＳＯ算法;计算区域的离散采用非交错斜交网格．计算了二维不可压缩Ｃｏｕｅｔｅ流动及无限长相对运动动／静叶排内的非定常流动,验证了用该方法和软件计算动／静相干叶排内非定常流动的有效性     

基于预处理的不可压缩N-S方程拟压缩数值算法研究

通过研究求解不可压缩Navier—Stokes方程的拟压缩方法与加速刚性双曲型方程时间推进的预处理技术,推导了一般曲线坐标系下带预处理的拟压缩Navier—Stokes方程特征系统,并结合有限差分法,建立了适用于不可压缩黏性流动计算的拟压缩快速算法．通过对不可压缩无黏圆柱绕流、平板层流流动、低Reynolds数定常圆柱绕流问题的数值模拟研究,得到了与相关理论与实验测试相吻合的结果,验证了所建立数值方法的快速可靠性．较系统地研究了预处理引入的参数和拟压缩因子的选择对收敛特性的影响．结果表明,Roe格式相对于二阶中心差分格式得到的结果更令人满意;对拟压缩Navier—Stokes方程进行预处理能有效提高数值计算的收敛速度;自适应的拟压缩因子取值能在很大程度上改善数值解的收敛特性,且不需要根据具体流动问题进行人工调节．最后将本文发展的数值方法用于低Reynolds数非定常圆柱绕流的数值模拟,所得结果亦和实验观测结果及其他文献的计算吻合很好．     

基于SMAC方法的90°弯管内部流场数值模拟

基于SMAC(simplified marker and cell)方法和交错网格技术,发展了一套在贴体坐标系下求解三维不可压缩流动的数值算法,并应用该方法求解以逆变速度和压力为未知变量的Navier-Stokes方程,对弯曲度为2.3的三维90°方形截面弯管内部不可压缩流场进行了数值模拟.结果表明:文中提出的SMAC方法对强曲率弯管内流有较强的预测能力,能够较准确地模拟其主流和两次流的运动规律,从而证明了本算法的可靠性和有效性.     

激波诱导变形液滴外流场数值模拟

基于Delaunay非结构网格生成方法和局部阵面推进重新生成网格的方法,应用边界边剖分与合并技术,对二维含自由表面流动中运动边界网格的生成方法进行了研究,改善了局部网格重新生成过程中的网格质量;对激波诱导变形水液滴在高速气流作用下的受力情况做了数值模拟.结果表明:该文在处理二维自由表面的动边界问题中采用的方法是可行的.     

基于动态嵌套网格方法的摆动翼型粘性绕流数值模拟

该文根据摆动翼型的流场特性运用动态嵌套网格技术将其划分为不同的子区,在不同的子区中使用不同的流动控制方程,以提高计算效率.对翼型生成一个随其摆动的相对较小较密的贴体网格,在它上面采用Navier Stokes方程;对整个流场生成了一个大的相对较稀的静止背景网格,在它上面采用Euler方程.两个网格间的流场信息交换通过动态嵌套网格技术来实现.通过与单块翼型网格的计算过程和计算结果的对比分析,表明该方法在不损失计算精度的前提下计算非定常问题时具有相对较高的计算效率.     

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

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

二维移动网格方法中边界网格点移动的数值研究

本文对二维移动网格方法中边界网格点移动处理方法进行了探讨,一种做法是把边界上的网格点的移动速度等同于和它相邻的内部点的移动速度,另一种做法是将边界上的网格移动看作一维情形,结合数值算例分析和比较了两种处理方法的效果。     

回转状态下导管螺旋桨水动力特性实用数值模拟方法

以数值模拟手段研究作为水下机器人主要姿态与轨迹控制机构的导管螺旋桨在回转状态下的水动力特性。首先根据所要计算的导管螺旋桨的几何要素构造导管和螺旋桨的三维几何模型,在此基础上采用动网格技术以计算流体力学方法借助计算流体力学软件Fluent,采用有限体积法在导管和螺旋桨所在的流域内求解其N-S方程以此对在作回转运动的水下机器人流场作用下导管螺旋桨以一定转速和一定来流方向下所产生的推力特性进行数值模拟分析,观察在水下机器人主体流场干扰下螺旋桨推进器的水动力现象。为验证所采用的预测导管螺旋桨水动力特性数值方法的有效性,在Ka 4-70/19A标准导管螺旋桨已有实验数据的基础上,将数值模拟结果与试验结果进行比较,对比结果表明本文所采用的数值方法是可靠的。     

薄平板气动导数的数值计算

用基于一般曲线坐标系和交错网格的差分法求解原始变量二维不可压粘性流体的Ｎ－Ｓ方程,计算了宽高比为２２．５、两端圆滑的薄平板的气动导数．将薄平板当作刚性物体,首先分别计算物体做竖向强迫振动和扭转强迫振动的气动力,雷诺数取为１０５,物体运动和气流的耦合作用用动网格法考虑,Ｎ－Ｓ方程用二阶投影法和多重网格法求解．然后用最小二乘法确定Ｓｉｍｉｕ和Ｓｃａｎｌａｎ给出的气动力表达式中的８个气动导数．绝大部分气动导数的计算结果和Ｔｈｅｏｄｏｒｓｅｎ的理论值及实验值都十分吻合     

偏格子中的格子Boltzmann方程

给出偏格子的格子Boltzmann模型. 应用Chapman-Enskog展开和多重尺度技术, 通过 选择平衡态分布函数的高阶矩和偏张量的形式, 得到偏网格上的Euler方程.     

基于动网格的液压阀阀芯启闭中的液动力分析

为了深入研究液压阀阀芯启闭过程中阀芯上的压强分布和液动力,采用动网格技术,数值模拟了液压球阀阀芯开启与闭合过程中的非定常流场.将阀芯作为运动实体,阀芯的运动由其所受的流体力、弹簧力及自身重力确定,随着阀芯的运动,相应流场计算区域的边界也随着发生变化,获得了球阀启闭过程中不同瞬时阀腔内的压力场.在此基础上,分析了液压阀在...     

自适应的移动单位分解法求解边界层问题

提出一种新的无网格方法,即自适应移动单位分解法来处理边界层问题,该方法的主要思想是将移动网格技巧与单位分解法结合起来.利用等弧长分配自适应地移动节点分布,在此基础上使用单位分解格式求解一维对流扩散模型,并给出方法的理论基础与算法流程.文中数值实验结果验证了该自适应移动方法的有效性和优越性.     

二阶双曲型方程变网格有限元方法的后验误差分析

讨论双曲型方程变网格有限元方法，给出了它的后验误差估计，实际计算中的局部网格调整，并对后验估计进行合理性分析．