共查询到17条相似文献,搜索用时 187 毫秒
1.
2.
3.
p型多重网格间断Galekin有限元方法研究 总被引:1,自引:0,他引:1
在二维非结构网格上,使用p型多重网格间断Galerkin方法求解定常可压缩欧拉方程。p型多重网格方法主要特征是通过对不同阶次多项式的近似解进行递归迭代求解。文中高阶近似(p0)上使用显式格式,在最低阶近似(p=0)上选用隐式格式,而非显式格式,从而在保证精度和占用较小内存的情况下加速收敛到定常解。运用该方法对NACA0012跨音速无粘流动进行数值模拟,数值结果表明:p型多重网格方法同单重显式Runge-Kutta方法相比,收敛速度能够提高6倍左右,并且精度保持不变。 相似文献
4.
在使用间断Galerkin有限元方法的计算过程中,需要构造相应的积分表达式作为数值求解的出发点,继而会引入体积分和面积分。对于这些积分项的值,一般需要通过数值积分的方法获得。当需要使用高阶间断Galerkin有限元方法时,数值积分计算精度的要求会相应地增加,其所需的计算量将变得很大,而数值积分的计算量又在很大程度上决定了间断Galerkin有限元方法的计算效率。针对这一问题,通过建立Lagrange插值多项式基函数和Jacobi正交多项式基函数的一定关系,构造了一种无数值积分的间断Galerkin有限元方法显式半离散格式,并对不同条件下的线性和非线性一维、二维守恒律进行了直接数值模拟,得到了理想的数值结果。该方法不再需要通过数值积分来计算每个单元的积分项,而且有效地达到了间断Galerkin有限元方法的高阶精度,其对于构造更为高效的高阶间断Galerkin有限元计算方法具有非常显著的意义。 相似文献
5.
6.
7.
8.
采用间断有限元方法对雷诺平均Navier-Stokes(RANS)方程进行了数值求解,对Spalart-Allmaras单方程湍流模型进行了部分修正,使得求解器更加鲁棒。构造了分段高次多项式来逼近真实物面,同时物面附近采用多层弯曲网格来避免网格交叉,此外提出了一种快速计算积分点的曲面物面距的方法。采用混合网格对NACA0012翼型以及RAE翼型进行了数值计算,并与实验数据以及前人数据进行了对比。计算结果表明,采用物面弯曲网格结合修正的湍流模型方法在相对稀疏的网格上就能得到比较好的数值解。 相似文献
9.
拓展了二维间断 Galerkin(DG)有限元方法研究,将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元,采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元,更好地表述了真实物面特征;物面边界相邻体网格单元相应构造为高阶体网格单元,其余体网格单元采用八节点六面体单元,以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解,验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外,由于采用了隐式计算方法,仅需较少的时间步就能迭代收敛。 相似文献
10.
无单元Galerkin法需要在背景网格上积分,计算量大,且在求解对流占优问题时会出现非物理的数值伪振荡现象。为此,基于局部Taylor展开思想,采用节点处的局部Taylor展开计算积分,建立了局部Taylor展开积分无单元Galerkin法。该方法同时解决了标准的无单元Galerkin法计算量大和对流占优时会出现数值伪振荡的问题。一维定常对流扩散方程和二维Burgers方程的求解说明了该方法的有效性,且计算效率远高于无单元Galerkin法。 相似文献
11.
在二维非结构网格上,对高阶间断Galerkin方法求解定常二维欧拉方程进行数值研究。为了很好地抑制解在流场间断附近处产生的伪振荡,采用了间断探测器和限制器相结合的方法,并分别对Krivodonova间断探测器和基于激波物理特征的激波探测器进行了比较和研究。数值结果表明:如果流场中只存在激波时,两者探测效果相似,且后者更适合应用于求解激波问题。如果流场更加复杂,即存在多种间断时,前者依然可以比较准确地用来探测间断单元,能够很好抑制流场间断附近处产生的伪振荡,而后者失效。 相似文献
12.
13.
由于航空发动机内高温热端部件是非规则形状的,其中的热辐射过程存在明显的间断遮蔽效应,如何实现该过程的精准模拟是一个挑战。采用间断有限元对辐射传递方程中角度和空间计算域进行离散,并基于分层限制的思想在角度和空间单元上引入巴斯-叶斯帕森限制器,对辐射强度数值解在角度和空间上的非物理振荡进行抑制。通过与文献中有限体积法和蒙特卡洛辐射模型的结果对比,证明了方法的有效性。此外,通过与辐射强度在角度上的解析解对比,验证了该限制器可以抑制辐射间断效应造成的数值振荡,消除辐射强度的非物理解。实现了任意空间位置上辐射强度在角度方向上的三维刻画。 相似文献
14.
《中国航空学报》2016,(6):1496-1505
High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is inves-tigated in combination with discontinuous Galerkin (DG) method for solving the 3D compressible Euler and Navier-Stokes equations. In this method, a fine linear mesh is first generated by standard commercial mesh generation tools. By taking advantage of an agglomeration method, a quadratic high-order mesh is quickly obtained, which is coarse but provides a high-quality geometry represen-tation, thus very suitable for high-order computations. High-order discretizations are performed on the obtained grids with DG method and the discretized system is treated fully implicitly to obtain steady state solutions. Numerical experiments on several flow problems indicate that the agglomer-ated high-order mesh works well with DG method in dealing with flow problems of curved geome-tries. It is also found that with a fully implicit discretized system and a p-sequencing method, the DG method can achieve convergence state within several time steps which shows significant effi-ciency improvements compared to its explicit counterparts. 相似文献
15.
16.
考虑到间断有限元方法对边界的敏感性,采用基于八节点曲边四边形单元的间断有限元方法求解了Euler方程的圆柱绕流问题。详细推导了八节点四边形单元的变换关系,给出了Jacobi矩阵行列式的具体表达式。对比直边单元和曲边单元的计算结果,采用曲边单元后,计算结果符合Euler方程的无粘假设,得出了八节点四边形单元间断有限元方法求解Euler方程是合适的结论。 相似文献