求解Euler／Navier—Stokes方程的有限体积龙格库塔方法

提出了一个求解定常Ｅｕｌｅｒ／Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的计算方法,求解的主管方程是积分形式的Ｅｕｌｅｒ方程或雷诺平均的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,计算方法是有限体积多步龙格库塔方法,通过在格式中加入自适应的二阶加四阶人工粘性项以保证数值计算的稳定性,用空间变时间步长、焓阻尼来加速定常问题的收敛速度。     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

三维Euler/N-S方程的多重网格求解方法

1引言 求解飞行器复杂绕流的三维Euler/N-S方程,需要大量的机时,过去求解全机N-S方程,计算一个状态往往需要几天,甚至一周以上的时间,这大大地制约了CFD技术在实际工程问题中的应用.即使在计算机技术迅速发展的今天,提高计算效率,缩短计算时间,对CFD技术应用于工程实际仍然具有重要意义.多重网格技术,是提高计算效率的一种非常有效手段,据国外的经验,应用多重网格技术,可成倍,几倍甚至十倍地减少计算机时.但国内对这方面的研究还很不够,能成功应用于已开发的软件系统的则更为缺乏.本文的研究,就是为解决这一问题提供技术储备.     

求解Burgers方程的特征中心型有限体积法

以Burgers方程为试验模型,提出一种新的有限体积法,空间上的离散采用中心型加权基本无振荡重构,时间上的离散采用数值积分。其中积分节点上的数值通量由特征理论回溯求解,保持了物理量沿特征方向传输的特性,计算量相对Runge-Kutta法明显减少。数值结果表明了方法的有效性和稳定性。     

空间推进方法求解抛物化Navier-Stokes方程及其验证

采用伪时间空间推进方法求解完全气体PNS方程。PNS方程采用有限体积方法离散,在推进方向上采用一阶精度迎风格式,在推进面的展向和法向采用三阶精度的MUSCL插值和AUSM类通量构造格式;在推进面上的伪时间推进采用二维LU-SGS迭代方法。通过算例证明,空间推进方法能得到正确的压力、摩阻和热流分布,且计算时间比时间迭代方法快一个数量级以上。     

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

曲边单元间断有限元方法求解二维Euler方程

考虑到间断有限元方法对边界的敏感性,采用基于八节点曲边四边形单元的间断有限元方法求解了Euler方程的圆柱绕流问题。详细推导了八节点四边形单元的变换关系,给出了Jacobi矩阵行列式的具体表达式。对比直边单元和曲边单元的计算结果,采用曲边单元后,计算结果符合Euler方程的无粘假设,得出了八节点四边形单元间断有限元方法求解Euler方程是合适的结论。     

极坐标系下谱元方法求解不可压缩Navier-Stokes方程

本文利用有限元的思想并结合谱方法的精度提出求解偏微分方程的谱元方法,并将谱元方法应用到极坐标系中;详细推导了在极坐标系下的谱元方法的具体计算公式,求解了极坐标系下的简单椭圆型二阶偏微分方程;并结合时间分裂方法应用于同心旋转圆筒间流体的流动问题,具体求解了原始变量速度和压力的不可压缩NavierStokes方程,均取得了满意的结果.     

基于三阶龙格库塔法的铣削稳定性半解析法预测

针对使用再生颤振理论建立的铣削动力学模型,提出了一种基于三阶龙格库塔法用于预测铣削稳定性的半解析方法.首先,以状态空间方程的形式表示动力学微分方程;其次,利用三阶龙格库塔法推导出传递矩阵;最后利用Floquet理论判断特定切削状态下的稳定性,进而获得铣削稳定性叶瓣图.通过与半离散法的仿真结果进行对比发现,基于三阶龙格库塔法的铣削稳定性求解方法具有更高的预测精度和计算效率.     

Navier-Stokes方程两层稳定有限元算法分析

分析了定常Navier-Stokes方程的两种两层稳定有限元算法。它们将局部Gauss积分稳定化技术和两层算法的思想充分结合,采用低次等阶有限元P1-P1或Q1-Q1对N-S方程进行数值求解。误差分析和数值算例结果表明,当粗、细网格尺度H=O（h1/2）时,它们与在细网格上的单层有限元算法具有相同的收敛速度,而两层算法却节省了大量的计算时间。相比之下,Simple算法具有更高的计算效率。而且进一步发现Oseen算法能够对小粘性系数N-S方程进行有效求解。     

非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

二维可压缩Euler方程组的SUPG有限元数值解

本文在文献[1]的基础上,构造了二维非定常可压缩Euler方程组的SUPG变分方程组。文中对由文献[1]提出的预估——校正算法作了分析,并用该算法对变分方程组进行解算,所做算例分别为:无粘性激波在固壁上的反射、零攻角10％抛物翼型亚、跨声速绕流、带攻角NACA 0012翼型跨声速绕流,取得了较好的数值结果。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

跨声速三维Euler方程的显式／隐式格式有限元解法

本文采用有限体积Ｇａｌｅｒｋｉｎ法和非结构网格数值求解了跨声速三维Ｅｕｌｅｒ方程,通量的求解采用ＶａｎＬｅｅｒ－阶逆风矢能量发裂方式,并通过外插使其上升为二阶精度,分别采用了显式及隐式两种格式求解。     

跨声速欧拉方程的一种隐式有限元解法

本文用有限体积Ｇａｌｅｒｋｉｎ法求解了跨声速二维欧拉方程。求通量运用了Ｏｓｈｅｒ格式，通过外插把原来只有一阶精度的逆风格式改进成二阶精度。并采用基于一性化流量的隐式格式和非结构化网格，使程序具有较高的收敛效率。     

非定常不可压粘性/无粘性耦合方程的一种分步分解方法

给出了数值求解初始变量不可压Navier-Stokes／Euler耦合方程的一种分步块LU分解方法。与传统的时间分裂法不同,该法无需压力中介边条件,从而避免了传统时间分裂法要求的复杂的压力中介边条件逼近。分步块LU分解方法可看做经典的Uzawa算法的改进,后者曾被成功应用于不可压Navier-Stokes／Euler耦合方程的求解。但本文显示分步块LU分解法比经典的Uzawa方法更经济。分析显示该法具有良好的稳定性和高精度,数值结果支持这一理论分析。