共查询到17条相似文献,搜索用时 171 毫秒
1.
首先通过理论推导给出了三阶WENO格式(WENO-JS3格式)满足收敛精度的充分条件.采用Taylor(泰勒)级数展开的方法,分析发现传统的三阶WENO-Z格式(WENO-Z3格式)在光滑流场极值点处精度降低.为了提高WENO-Z3格式在极值点处的计算精度,根据收敛精度的充分条件构造一种改进的三阶WENO-Z格式(WENO-NZ3格式),并综合权衡计算精度和计算稳定性确定所构造格式的参数.通过两个典型的精度测试,验证了WENO-NZ3格式在光滑流场极值点区域逼近三阶精度.选用Sod激波管、激波与熵波相互作用、Rayleigh-Taylor不稳定性、二维Riemann(黎曼)问题经典算例,进一步证实了本文提出的WENO-NZ3格式相较其他格式(WENO-JS3、WENOZ3、WENO-N3),不仅提高了计算精度,而且提高了对复杂流场结构的分辨率. 相似文献
2.
两点边值问题基于应力佳点的一类二次有限体积元方法 总被引:1,自引:0,他引:1
本文针对常微分方程两点边值问题提出了一种基于应力佳点的二次有限体积元格式,证明了该格式按离散能量模具有三阶收敛精度.具体算例表明该格式计算效果良好. 相似文献
3.
为了降低经典的三阶加权本质无振荡(WENO)格式的数值耗散,提出了一种新的三阶WENO格式的修正模板近似方法.改进了经典WENO-JS格式中各候选模板上数值通量的一阶多项式逼近,通过加入二次项使模板逼近达到三阶精度.计算了相应的候选通量,并且通过引入可调函数φ(x),使得新的格式具有ENO性质.最后给出了一系列数值算例,证明了该方法的有效性. 相似文献
4.
提出了求解三维抛物型方程的一个高精度显式差分格式.首先,推导了一个特殊节点处一阶偏导数(■u)/(■/t)的一个差分近似表达式,利用待定系数法构造了一个显式差分格式,通过选取适当的参数使格式的截断误差在空间层上达到了四阶精度和在时间层上达到了三阶精度.然后,利用Fourier分析法证明了当r1/6时,差分格式是稳定的.最后,通过数值试验比较了差分格式的解与精确解的区别,结果说明了差分格式的有效性. 相似文献
5.
6.
针对计算流体力学对高精度高分辨率的需求,基于降低经典的三阶加权本质无振荡(WENO)格式的数值耗散特性,该文提出了一种新的参考光滑性指示子.其构造方法与经典的WENO-Z格式不同,它是通过候选子模板上重构多项式的导数的线性组合与整个全局模板上重构多项式的导数的L^(2)范数逼近获得的.采用该计算方法可以得到比WENO-Z格式更高阶的参考光滑性指示子,另外改变自由参数φ的取值,可以获得不同的参考光滑性指示子.该文通过一系列数值算例证明了该参考光滑性指示子的有效性. 相似文献
7.
低耗散、高分辨率激波捕捉格式对含激波流场的数值模拟具有重要意义.在传统三阶WENO格式(WENO-JS3)和三阶WENO Z格式(WENO-Z3)基础上,基于映射函数,给出WENO-M3、WENO MZ3格式.选用Sod激波管、激波与熵波相互作用、双爆轰波碰撞及双Mach(马赫)反射等经典算例,考察上述格式的计算性能.数值结果表明,WENO-MZ3格式相较其他格式具有耗散低、对流场结构分辨率高的特性.为了进一步扩展WENO-MZ3格式的应用范围,采用该格式数值研究封闭方形舱室内柱形高压、高密度气体爆炸波传播过程,波系演化规律以及壁面典型测点压力载荷.数值计算结果表明WENO-MZ3格式能够较好地模拟包含高压比、高密度比的爆炸波且给出数值耗散较小的壁面压力载荷. 相似文献
8.
9.
本文考虑一维单个守恒律方程,对其设计了一个基于熵耗散的非线性守恒型差分格式.本格式的数值流函数是Lax-Freidrichs格式和Lax-Wendroff格式数值流函数的凸组合,凸组合中的系数是由考虑耗散熵来决定的.这样在解的光滑区域内,格式几乎、甚至完全是Lax-Wendroff格式,而在解的间断处,格式几乎、甚至完全是Lax—Freidrichs格式.从而消除了间断附近的非物理振荡,实现了计算的非线性稳定性.理论分析表明本格式在解的非极值点处是二阶精度的,而在解的极值点处至少有一阶精度.数值试验表明格式是有效的. 相似文献
10.
通过一个典型的Bratu问题,研究了小波Galerkin法(WGM)在非线性分岔问题求解方面的应用.首先,利用基于Coiflet的小波Galerkin法,对一维和二维Bratu方程进行离散;然后针对单参数问题,推导了追踪解曲线的伪弧长格式和直接计算极值型分岔点的扩展方程;针对双参数问题,推导了追踪稳定边界的伪弧长格式和求解尖点型分岔点的扩展方程.数值结果表明,基于小波Galerkin法的非线性分岔计算不仅具有更高的计算精度,而且能够有效地捕捉双参数分岔问题的折迭线和尖点突变曲面.该算例展示了基于小波Galerkin法的数值分岔计算的具体过程及其求解多参数分岔问题复杂行为的应用潜力. 相似文献
11.
随着计算流体力学的快速发展,设计精确、高效并且健壮的数值格式变得尤为重要.通过对3种流行的通量分裂方法(AUSM、Zha-Bilgen和Toro-Vázquez)的对流通量和压力通量进行特征分析,构造了一种简单、低耗散并且健壮的通量分裂格式(命名为R-ZB格式).采用Zha-Bilgen分裂方法将Euler方程的通量分裂成对流通量和压力通量,其中对流通量采用迎风方法来计算,压力通量采用低耗散的HLL格式来计算,从而克服了原始的HLL格式不能精确分辨接触间断的缺点.数值实验表明,该文给出的R-ZB格式不仅保留了原始Zha-Bilgen格式简单高效、能够精确分辨接触间断等优点,而且具有更好的健壮性,在计算二维问题时不会出现数值激波不稳定现象. 相似文献
12.
13.
激波捕捉差分方法研究 总被引:1,自引:1,他引:0
在迎风型格式和矢通量分裂技术的基础之上,对捕捉激波方法进行一种新的尝试.该方法首先对原始格式在特征方向上进行投影,然后用限制器对这些特征分量的变化幅值进行限制以抑止非物理波动,最后再把它转换成守恒形式,得到了基本上无振荡的激波捕捉格式.用该方法对两种迎风显示格式(二阶和三阶)和3种迎风紧致格式(三阶、五阶和七阶)进行处理,并在一维和二维的情况下进行了应用测试.通过与高阶WENO、MP、Compact-WENO等格式的比较,表明该方法在光滑捕捉激波的前提下仍有较高精度和分辨率. 相似文献
14.
15.
In this paper Homotopy Analysis Method (HAM) is implemented for obtaining approximate solutions of (2+1)-dimensional Navier-Stokes equations with perturbation terms. The initial approximations are obtained using linear systems of the Navier-Stokes equations;by the iterations formula of HAM,the first approxima-tion solutions and the second approximation solutions are successively obtained and Homotopy Perturbation Method(HPM)is also used to solve these equations;finally, approximate solutions by HAM of (2+1)-dimensional Navier-Stokes equations with-out perturbation terms and with perturbation terms are compared. Because of the freedom of choice the auxiliary parameter of HAM,the results demonstrate that the rapid convergence and the high accuracy of the HAM in solving Navier-Stokes equa-tions;due to the effects of perturbation terms,the 3rd-order approximation solutions by HAM and HPM have great fluctuation. 相似文献
16.
Weizhong Dai 《Numerical Methods for Partial Differential Equations》2010,26(5):1079-1098
We consider a mathematical model for thermal analysis in a 3D N‐carrier system with Neumann boundary conditions, which extends the concept of the well‐known parabolic two‐step model for micro heat transfer. To solve numerically the complex system, we first reduce 3D equations in the model to a succession of 1D equations by using the local one‐dimensional (LOD) method. The obtained 1D equations are then solved using a fourth‐order compact finite difference scheme for the interior points and a second‐order combined compact finite difference scheme for the points next to the boundary, so that the Neumann boundary condition can be applied directly without discretizing. By using matrix analysis, the compact LOD scheme is shown to be unconditionally stable. The accuracy of the solution is tested using two numerical examples. Results show that the solutions obtained by the compact LOD finite difference scheme are more accurate than those obtained by a Crank‐Nicholson LOD scheme, and the convergence rate with respect to spatial variables is about 2.6. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010 相似文献