共查询到18条相似文献,搜索用时 125 毫秒
1.
利用双曲守恒律的Hamilton-Jacobi方程形式,应用Taylor公式与Galerkin有限元给出了求解双曲守恒律的计算方法。采用TVD差分格式的构造思想,对数值通量作修正,在等距网格情形下有限元方法得到的计算格式满足TVD性质,并给出了数值例子。 相似文献
2.
为了减小TVD格式的数值粘性,提高TVD格式的分辨度,提出了一种具有较小数值耗散的Harten型TVD格式.数值模拟结果表明,这种Harten型TVD格式比原来的Harten格式[1]及其改进形式Harten-Yee[2]的TVD格式能更好地模拟出边界层及涡流. 相似文献
3.
本文对存在汽、液平衡态相变的两相流动建立了完全欧拉坐标系统下的数理模型。通过引入考虑真实流体性质的数值算法,并直接从IAPWS水及水蒸汽性质数据库中获取流体工质的性质,使数值计算的精度得到显著提高。采用包括LU-SGS-GE隐式格式和改良型高精度、高分辨率的MUSCL TVD格式的时间推进算法求解平衡态两相流动控制方程组以及低Reynolds数双方程湍流模型,对某汽轮机末级静叶进行了数值模拟,计算结果表明本文采用的模型及方法在某些条件下可以对叶栅主要性能参数进行准确的预测。 相似文献
4.
三维不可压N-S方程的多重网格求解 总被引:2,自引:0,他引:2
应用全近似存储(Full Approximation Storage,FAS)多重网格法和人工压缩性方法求解了三维不可压Navi-er-Stokes方程.在解粗网格差分方程时,对Neumann边界条件采用增量形式进行更新,离散方程用对角化形式的近似隐式因子分解格式求解,其中空间无粘项分别用MUSCL格式和对称TVD格式进行离散.对90°弯曲的方截面管道流动和4:1椭球体层流绕流的数值模拟表明,多重网格的计算时间比单重网格节省一半以上,且无限制函数的MUSCL格式比TVD格式对流动结构有更好的分辨能力. 相似文献
5.
平面叶栅气膜冷却流动的数值模拟 总被引:1,自引:0,他引:1
为了能够准确地对透平叶栅气膜冷却效率进行数值预测,本文采用了FNM形式的结构化网格,对一个平面叶栅中的气膜冷却流场进行了数值模拟。计算中采用了包括LU-SGS-GE隐式格式和改良型高精度、高分辨率的MUSCL TVD格式的时间推进算法求解三维RANS方程以及低Reynolds数q-ω双方程湍流模型。计算结果表明本文采用的模型及方法在低吹风比的条件下可以较准确地对气膜冷却效率进行数值预测。 相似文献
6.
高分辨率有限体积LW格式 总被引:2,自引:0,他引:2
一、引言 采用通量差分分裂技术或通量分裂技术设计出的高分辨率格式,推广应用到双曲型气动方程组时,编程要比常规的中心差分格式复杂很多。使用校正通量技术设计的TVD格式,集中了TVD格式分辨率高,中心差分格式编程容易的优点。在中心差分格式的基础上,根据TVD理论添加新的数值耗散项,通过感知流场的变化,耗散项的大小能够自动地得以调整和控制,使基恰好能抑制住流场中出现的振荡。与常规的人工耗散(AD)不同,新的数值耗散具有感受一反馈机制,不含和问题相关的自由参数,为了和人工耗散相区别,不妨称之为智能耗散(ID)。 相似文献
7.
将守恒型单调迎风格式(MUSCL)和双曲正切界面捕捉(THINC)算法组合到一起进行数值重构,并通过重构建立新的算子。重构的算子主要作用是尽可能降低网格的边界差商。在MUSCL和THINC算法下,网格左右边界处的状态可以导出一个优化算法的差商判定机制。再按照这一边界差商优化值最小(BVD)的原则,对MUSCL和THINC进行二选一。在BVD重构算子的作用下,数值断面变得平滑,并且有效地限制振荡。BVD重构算子的计算性能将通过一维和二维数值算例来进行检验。 相似文献
8.
9.
10.
采用双曲函数展开法得到Modified Benjamin-Bona-Mahony(mBBM)方程的一类扭结-反扭结状的双扭结孤立波解,在不同的极限情况下,此孤立波分别退化为mBBM方程的扭结状和钟状孤立波解.对双扭结型单孤子的结构特征进行分析,构造有限差分格式对其动力学稳定性进行数值研究.有限差分格式为两层隐式格式,在线性化意义下无条件稳定.数值结果表明mBBM方程的双扭结型单孤子在不同类型的扰动下均具有很强的稳定性.对双孤立波的碰撞进行数值模拟,发现既存在弹性碰撞也存在非弹性碰撞. 相似文献
11.
12.
二维交错网格的GAUSS型格式 总被引:2,自引:0,他引:2
利用Gauss型求积公式在交错网格的情况下构造了一类不需解Riemann问题的求解二维双曲守恒律的二阶显式Gauss型差分格式,该格式在CFL条件限制下为MmB格式.并将格式推广到二维方程组,进行了数值试验. 相似文献
13.
A large time step (LTS) TVD scheme originally proposed by Harten is modified and further developed in the present paper and applied to Euler equations in multidimensional problems. By firstly revealing the drawbacks of Harten’s original LTS TVD scheme, and reasoning the occurrence of the spurious oscillations, a modified formulation of its characteristic transformation is proposed and a high resolution, strongly robust LTS TVD scheme is formulated. The modified scheme is proven to be capable of taking larger number of time steps than the original one. Following the modified strategy, the LTS TVD schemes for Yee’s upwind TVD scheme and Yee–Roe–Davis’s symmetric TVD scheme are constructed. The family of the LTS schemes is then extended to multidimensional by time splitting procedure, and the associated boundary condition treatment suitable for the LTS scheme is also imposed. The numerical experiments on Sod’s shock tube problem, inviscid flows over NACA0012 airfoil and ONERA M6 wing are performed to validate the developed schemes. Computational efficiencies for the respective schemes under different CFL numbers are also evaluated and compared. The results reveal that the improvement is sizable as compared to the respective single time step schemes, especially for the CFL number ranging from 1.0 to 4.0. 相似文献
14.
15.
The broken dam problem flow is tested to check accuracy of different procedures for gas-liquid interface resolution based
on solution of the additional equation for the volume fraction of liquid phase. The study is focused on the numerical schemes
used to approximate advection fluxes of this equation. In particular, the MUSCL scheme with QUICK interpolants and compressive
minmod TVD limiters with the slope modification technique for the volume fraction fluxes is applied, as well as the upwind-downwind
donor acceptor procedure designed in the VOF method. As the first stage, the quite simple and explicit procedure adopting
the artificial compressibility method is used to solve the velocity and pressure equations. Computations are initially performed
with a careful grid and time step independence studies. Importance of the wall boundary condition is also discussed. To present
free surface motion, results of numerical investigation are shown in terms of contour plots for the volume fraction at successive
times, as well as surge front and column height positions versus time.
The work was financially supported by the President of the Russian Federation (NSh No. 454.2008.1), Russian Foundation for
Basic Research (Grant No. 06-01-00724), Integration Project of SB RAS (No. 2-16), and National Scientific Council of Taiwan
(NSC, R.O.C., contract NSC-92-2212-E006-102). 相似文献
16.
双曲型守恒律的一种高精度TVD差分格式 总被引:3,自引:0,他引:3
构造了一维双曲型守恒律方程的一个高精度高分辨率的守恒型TVD差分格式.其主要思想是:首先将计算区域划分为互不重叠的小单元,且每个小单元再根据希望的精度阶数分为细小单元;其次,根据流动方向将通量分裂为正、负通量,并通过小单元上的高阶插值逼近得到了细小单元边界上的正、负数值通量,为避免由高阶插值产生的数值振荡,进一步根据流向对其进行TVD校正;再利用高阶Runge KuttaTVD离散方法对时间进行离散,得到了高阶全离散方法.进一步推广到一维方程组情形.最后对一维欧拉方程组计算了几个算例. 相似文献
17.
设计一种基于三单元具有六阶精度的修正Hermite-ENO格式(CHENO),求解一维双曲守恒律问题.CHENO格式利用有限体积法进行空间离散,在空间层上,使用ENO格式中的Newton差商法自适应选择模板.在重构半节点处的函数值及其一阶导数值时,利用Taylor展开给出修正Hermite插值使其提高到六阶精度,并设计了间断识别法与相应的处理方法以抑制间断处的虚假振荡;在时间层上采用三阶TVD Runge-Kutta法进行函数值及一阶导数值的推进.其主要优点是在达到高阶精度的同时具有紧致性.数值实验表明对一维双曲守恒律问题的求解达到了理论分析结果,是有效可行的. 相似文献
18.
Four explicit type time marching methods, including one proposed by the authors, are examined. The TVD conditions of this method are analyzed with the linear conservation law as the model equation. Performance of these methods when applied to the Euler equations are numerically tested. Seven examples are tested, the main concern is the performance of the methods when discontinuities with different strengths are encountered. When the discontinuity is getting stronger, spurious oscillation shows up for three existing methods, while the method proposed by the authors always gives the results with satisfaction. The effect of the limiter is also investigated. To put these methods in the same basis for the comparison the same spatial discretization is used. Roe's solver is used to evaluate the fluxes at the cell interface; spatially second-order accuracy is achieved by the MUSCL reconstruction. 相似文献