共查询到14条相似文献,搜索用时 125 毫秒
1.
2.
3.
4.
给出了求解一维双曲型守恒律的一种半离散三阶中心迎风格式,并利用逐维进行计算的方法将格式推广到二维守恒律。构造格式时利用了波传播的单侧局部速度,三阶重构方法的引入保证了格式的精度。时间方向的离散采用三阶TVD Runge—Kutta方法。本文格式保持了中心差分格式简单的优点,即不需用Riemann解算器,避免了进行特征分解过程。用该格式对一维和二维守恒律进行了大量的数值试验,结果表明本文格式是高精度、高分辨率的。 相似文献
5.
将Jin's的界面方法应用到求解双曲守恒型方程的半离散中心迎风方法中,给出了一种新的求解浅水波方程的半离散中心迎风差分方法。对于源项,不是采用传统的单元均值而是采用单元界面处的值来近似,使所得格式对稳定态的求解是均衡的。且已证明所给的二阶精度的求解格式保持水深的非负性,这一特性使其能够较好的处理干河床问题。使用该方法产生的数值粘性(与O(Δ2r-1)同阶)要比交错的中心格式小(与O(Δx2r/Δt)同阶),而且由于数值粘性与时间步长无关,从而时间步长可根据稳定性需要尽可能的小,因此适用于稳定态的求解。 相似文献
6.
7.
8.
利用松弛近似,将非线性的凝聚炸药爆轰控制方程转化为线性的松弛方程组,并采用五阶WENO格式和五阶线性多步显隐格式对线性松弛方程组进行空间方向和时间方向的离散,由此建立具有高精度和高分辨率性质的计算凝聚炸药爆轰的松弛方法。建立的松弛方法可以避免求解Riemann问题及计算非线性通量的Jacobi矩阵,同时无需分裂处理反应源项。通过对凝聚炸药的平面一维定常爆轰波结构及球面一维聚心、散心爆轰起爆和传播过程的数值模拟,验证了所建立的松弛方法能够很好地计算凝聚炸药爆轰问题。 相似文献
9.
基于HLL-HLLC的高阶WENO格式及其应用研究 总被引:1,自引:0,他引:1
HLL-HLLC格式能够克服HLLC在强激波附近的激波不稳定现象,并且保持了HLLC的低耗散特性,是一种适合更大马赫数范围的近似黎曼求解器。本文从RANS方程出发,将HLL-HLLC近似黎曼求解器结合五阶WENO重构,实现了对无粘通量的高阶离散;同时,采用完全守恒形式的四阶中心差分格式处理粘性项,建立了RANS 方程的高阶数值求解格式。通过对四个经典算例,钝头体、 ONERA M6机翼、DLR F6-WB翼身组合体和DLR F6-WBNP复杂外形的数值模拟,考察了两种WENO改进格式在复杂流场中的表现,研究了高阶格式的收敛特性;给出了在复杂流动中 WENO自由参数的推荐值,以增强求解的收敛性。算例结果表明,本文构造的高阶格式鲁棒性好,能够显著改善激波位置和激波强度,捕获更丰富的流场细节,满足复杂工程应用需求。 相似文献
10.
HLL-HLLC格式能够克服HLLC在强激波附近的激波不稳定现象,并且保持了HLLC的低耗散特性,是一种适合更大马赫数范围的近似黎曼求解器。本文从RANS方程出发,将HLL-HLLC近似黎曼求解器结合五阶WENO重构,实现了对无粘通量的高阶离散;同时,采用完全守恒形式的四阶中心差分格式处理粘性项,建立了RANS方程的高阶数值求解格式。通过对四个经典算例,钝头体、ONERA M6机翼、DLR F6-WB翼身组合体和DLR F6-WBNP复杂外形的数值模拟,考察了两种WENO改进格式在复杂流场中的表现,研究了高阶格式的收敛特性;给出了在复杂流动中WENO自由参数的推荐值,以增强求解的收敛性。算例结果表明,本文构造的高阶格式鲁棒性好,能够显著改善激波位置和激波强度,捕获更丰富的流场细节,满足复杂工程应用需求。 相似文献
11.
We introduce a new fourth order, semi-discrete, central-upwind scheme for solving systems of hyperbolic conservation laws. The scheme is a combination of a fourth order non-oscillatory reconstruction, a semi-discrete central-upwind numerical flux and the third order TVD Runge-Kutta method. Numerical results suggest that the new scheme achieves a uniformly high order accuracy for smooth solutions and produces non-oscillatory profiles for discontinuities. This is especially so for long time evolution problems. The scheme combines the simplicity of the central schemes and accuracy of the upwind schemes. The advantages of the new scheme will be fully realized when solving various examples. 相似文献
12.
Mohammed M. Babatin 《International Journal of Computational Fluid Dynamics》2013,27(10):723-735
In this article, we develop an adaptive scheme for solving systems of hyperbolic conservation laws. In this scheme nonlinear shock and linear contact waves will be treated differently. The proposed scheme uses the Kurganov central-upwind scheme. Fourth-order non-oscillatory reconstruction is employed near shock only while the unlimited fifth-order reconstruction is used for smooth regions and linear contact waves. To distinguish between the smooth parts and discontinuities, we use an efficient adaptive multi-resolution technique. The advantages of the scheme are high resolution and computationally efficient since limiters are used only for shocks. Numerical experiments with one- and two-dimensional problems are presented which show the robustness of the proposed scheme. 相似文献
13.
The governing equations of shallow water magnetohydrodynamics describe the dynamics of a thin layer of nearly incompressible and electrically conducting fluids for which the evolution is nearly two-dimensional with magnetic equilibrium in the third direction. A high-resolution central-upwind scheme is applied to solve the model equations considering non-flat bottom topography. The suggested method is an upwind biased non-oscillatory finite volume scheme which doées not require a Riemann solver at each time step. To satisfy the divergence-free constraint, the projection method is used. Several case studies are carried out. For validation, a gas kinetic flux vector splitting scheme is also applied to the same model. 相似文献
14.
A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow. 相似文献