基于旋转坐标系的高阶间断有限元方法非定常湍流数值模拟

基于高阶间断有限元方法(Discontinuous Galerkin method,DGM),对旋转非惯性系下耦合了修正的一方程S-A模型的RANS方程进行了离散求解。为了在稀疏网格上获得更贴近真实的物面形状,使用了多层高阶弯曲网格方法对物面进行拟合。非定常时间推进采用了隐式双时间步方法,每个时间步产生的线性系统采用预处理的方法,即广义最小残差方法(Generalized minimal residual method,GMRES)来求解。计算了旋转圆柱绕流以及经典翼型振荡算例的升力和力矩迟滞曲线,与实验结果以及前人的计算结果对比验证了本文方法的正确性和有效性。     

一种模拟不可压流的节点间断有限元格子波尔兹曼方法

基于节点间断有限元方法,本文提出了一种求解格子波尔兹曼方程(Lattice Boltzmann equation,LBE)的新方法,即节点间断有限元-格子波尔兹曼方法(Nodal discontinuous Galerkin-Lattice Boltzmann method,NDGLBM)。在该方法中,LBE的碰撞过程和迁移过程被拆分成了两步:碰撞过程用LBM多松弛时间(Multiple relaxation time,MRT)模型进行求解,迁移过程则写成对流方程并采用节点间断有限元方法进行求解。其中,空间离散采用了非结构网格,时间离散采用了四阶、五步龙格-库塔格式。为了验证NDG-LBM的可行性,文中模拟了顶盖驱动方腔流、静止单个圆柱绕流、旋转-静止双圆柱绕流以及高雷诺数NACA0012翼型绕流。计算所得的数值结果与其他文献中的结果吻合度很好。     

非定常欠膨胀射流的间断有限元模拟

采用间断有限元方法对非定常欠膨胀射流进行了数值模拟,将二维守恒方程的间断有限元方法发展到轴对称Euler方程,对欠膨胀射流问题进行了数值计算。计算结果表明,采用间断有限元方法能够有效地捕捉到包括激波、滑移线、涡环和多级马赫盘等复杂流场结构,与实验照片反映的流动特征吻合较好。此外,两个马赫盘的位置与其他数值结果和实验测量结果吻合很好,表明该方法具有很好的鲁棒性,对非定常欠膨胀燃气射流的数值模拟是有效的。     

环形激波聚焦流场特性的间断有限元方法(英文)

针对环形激波聚焦时产生的高温、高压特性,采用间断有限元方法对环形激波在同轴圆柱形激波管内的聚焦流场进行了数值模拟。为了验证采用方法的有效性,首先对经典激波管问题进行了求解,计算结果与解析解吻合很好。其次,对入射激波马赫数为2.0到5.0的环形激波聚焦流场进行了计算,获得了环形激波在圆柱形激波管内聚焦过程的气动特性。模拟结果表明,采用间断有限元方法能够有效地捕捉激波聚焦过程形成的绕射、反射和聚焦等主要流动特征,而且在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力。最后比较了网格尺度对聚焦点峰值压力的影响。     

高阶龙格库塔间断有限元方法求解二维谐振腔问题

提出了使用高阶龙格库塔间断有限元方法在时域求解经典的电磁场谐振腔问题,间断有限元方法在空间离散时采用非结构化网格且在时域显式求解方程,这是有限元方法和有限体积方法的最佳结合。该方法通过采用局部高阶多项式插值基函数获得高阶精度。文中使用该方法研究了横磁波在二维谐振腔中的传播情况,以及高阶的拉格朗日基函数。数值实验采用了高阶的二变量拉格朗日多项式基函数,数值计算结果与理论解析解相吻合。文中还讨论了不同阶数多项式插值基函数对计算精度的影响。结果表明,随插值基函数阶数增加,计算精度迅速提高。最后讨论了不同插值基函数阶数对L^2误差的影响,结果显示L^2误差随插值基函数阶数增加呈指数下降。     

Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

基于间断伽辽金法的悬停状态旋翼涡流场低耗散数值模拟

建立间断伽辽金方法并将其应用于直升机旋翼桨尖涡演化和气动特性的精确模拟.首先,选取旋转参考坐标系下的雷诺平均N-S方程作为控制方程,采用嵌套网格方法对旋翼流场进行模拟,在桨叶网格和背景网格上分别采用有限体积法和间断伽辽金法.此外,空间离散采用高分辨率、低耗散的Harten-Lax-Van Leer contact(HL...     

一类矩阵方程的迭代解法

利用Kronecker积、矩阵拉直运算以及线性方程组的迭代解法,给出了解矩阵方程A1X1B1 A2X2B2 … AnXnBn=D的一种新的迭代方法,并给出了它的几种特殊形式.     

一类非线性双曲型方程的Galerkin方法

主要讨论了平面有界凸多角形区域上的一类非线性双曲型方程utt- .( a( x,u) u) =f ( x,u)u( x,0 ) =u0 ( x)ut( x,0 ) =φ( x)u( x,t) =0　　( x,t)∈Ω× [0 ,T]x∈Ωx∈Ω( x,t)∈ Ω× [0 ,T]的 Galerkin有限元方法 ,首先给出了所讨论问题的 Galerkin有限元方法的离散格式 ,其次对所讨论问题的解与其离散问题的解之间的误差进行了分析研究 ,最后利用椭圆投影算子的性质 ,得到了 L2 模和能量模方面的一些误差估计。     

High-Order Discontinuous Galerkin Solution of Compressible Flows with a Hybrid Lattice Boltzmann Flux

A discontinuous Galerkin(DG)-based lattice Boltzmann method is employed to solve the Euler and Navier-Stokes equations.Instead of adopting the widely used local Lax-Friedrichs flux and Roe Flux etc.,a hybrid lattice Boltzmann flux solver(LBFS)is employed to evaluate the inviscid flux across the cell interfaces.The main advantage of the hybrid LBFS is its flexibility for capturing both strong shocks and thin boundary layers through introducing a function which varies from zero to one to control the artificial viscosity.Numerical results indicate that the hybrid lattice Boltzmann flux solver behaves very well combining with the high-order DG method when simulating both inviscid and viscous flows.     

混合网格在求解三维N—S方程中的应用

用三棱柱／四面体组成的混合网格求解三维Navier-Stokes方程，在物面附近采用三陵柱网格，而其他地方则采用四面体网格，与单一类型网格相比，混合网格易于调节粘性层内网点数目，并能减少存储量，提高解的效率，Navier-Stokes方程的求解采用有限体积法，对人工耗散项进行了改进，使得求解的结果更接近真实的物理解，本文对上述方法开展了数值实验研究。结果表明，用三棱柱／四面体组成的混合网格求解N－S方程是非常有效的。     

关于二阶常系数线性椭圆型偏微分方程组解的唯一性的若干结论

本文证明了一个矩阵方面的有用结论,即文中定理２,说明了当条件（Ⅰ）、（Ⅱ）成立时,对于二个自变量、二个未知函数的二阶常系数线性方程组（１）可化为强椭圆型方程组,这一结论也可推广到某些三个未知函数的情形。利用强椭圆型方程组解必定唯一的结论,证明了某些二阶常系数线性椭圆型方程组在有界闭区域内Ｄｉｒｉｃｈｌｅｔ问题解的唯一性。     

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

超声速和高超声速粘性绕流QPNS方程的GLS有限元解

从完全非定常Ｎ－Ｓ方程（ＦＮＳ）中略去流向粘性导数导出拟抛物Ｎ－Ｓ方程（ＱＰＮＳ），ＱＰＮＳ技术结合经典抛物推进和拟时松弛，拟时松弛能有效地抑制解的发散。本文导出的熵变量形式ＱＰＮＳ方程具有对称性和自动满足热力第二定律，这将提高解的稳定性，伽辽金／最小二乘法（ＧＬＳ）被用来构造ＱＰＮＳ的弱解式。对于钝前缘物体，ＱＰＮＳ在前缘区并不适用，本文代之以用ＦＮＳ求解该区，由离散解得到的非对称线性方程组，对     

几种载体表面缝隙对雷达目标特性的影响

对在几种载体(平板和杏仁体)上的不同缝隙所引起的目标散射特性影响进行了研究。首先分析了不同载体目标的雷达散射截面特性;以此为对比基准,对含不同类型缝隙的载体目标开展散射特性计算和系列试验测试,从中获得由缝隙所引起的目标特性变化情况,进而分析缝隙对雷达目标特性的影响。结果显示,直槽会引起比直缝更强的散射增幅,锯齿缝相对直缝则具有明显的减缩作用;杏仁体载体能更有效地模拟弯曲表面上的缝隙散射。结果可为飞行器雷达隐身设计提供依据。