共查询到20条相似文献,搜索用时 15 毫秒
1.
This paper makes some numerical comparisons of time–space iterative method and spatial iterative methods for solving the stationary Navier–Stokes equations. The time–space iterative method consists in solving the nonstationary Stokes equations based on the time–space discretization by the Euler implicit/explicit scheme under a weak uniqueness condition (A2). The spatial iterative methods consist in solving the stationary Stokes scheme, Newton scheme, Oseen scheme based on the spatial discretization under some strong uniqueness assumptions. We compare the stability and convergence conditions of the time–space iterative method and the spatial iterative methods. Moreover, the numerical tests show that the time–space iterative method is the more simple than the spatial iterative methods for solving the stationary Navier–Stokes problem. Furthermore, the time–space iterative method can solve the stationary Navier–Stokes equations with some small viscosity and the spatial iterative methods can only solve the stationary Navier–Stokes equations with some large viscosities. 相似文献
2.
双曲型守恒律的一种高精度TVD差分格式 总被引:3,自引:0,他引:3
构造了一维双曲型守恒律方程的一个高精度高分辨率的守恒型TVD差分格式.其主要思想是:首先将计算区域划分为互不重叠的小单元,且每个小单元再根据希望的精度阶数分为细小单元;其次,根据流动方向将通量分裂为正、负通量,并通过小单元上的高阶插值逼近得到了细小单元边界上的正、负数值通量,为避免由高阶插值产生的数值振荡,进一步根据流向对其进行TVD校正;再利用高阶Runge KuttaTVD离散方法对时间进行离散,得到了高阶全离散方法.进一步推广到一维方程组情形.最后对一维欧拉方程组计算了几个算例. 相似文献
3.
4.
Pierre-Henri Maire 《Journal of computational physics》2009,228(18):6882-6915
The goal of this paper is to present high-order cell-centered schemes for solving the equations of Lagrangian gas dynamics written in cylindrical geometry. A node-based discretization of the numerical fluxes is obtained through the computation of the time rate of change of the cell volume. It allows to derive finite volume numerical schemes that are compatible with the geometric conservation law (GCL). Two discretizations of the momentum equations are proposed depending on the form of the discrete gradient operator. The first one corresponds to the control volume scheme while the second one corresponds to the so-called area-weighted scheme. Both formulations share the same discretization for the total energy equation. In both schemes, fluxes are computed using the same nodal solver which can be viewed as a two-dimensional extension of an approximate Riemann solver. The control volume scheme is conservative for momentum, total energy and satisfies a local entropy inequality in its first-order semi-discrete form. However, it does not preserve spherical symmetry. On the other hand, the area-weighted scheme is conservative for total energy and preserves spherical symmetry for one-dimensional spherical flow on equi-angular polar grid. The two-dimensional high-order extensions of these two schemes are constructed employing the generalized Riemann problem (GRP) in the acoustic approximation. Many numerical tests are presented in order to assess these new schemes. The results obtained for various representative configurations of one and two-dimensional compressible fluid flows show the robustness and the accuracy of our new schemes. 相似文献
5.
本文给出了一个模拟叶栅内准三维定常和非定常粘性流动的数值方法。对于定常流动,采用TVD Lax-Wendroff格式和代数湍流模型求解雷诺平均Navier-Stokes方程,使用当地时间步长和多网格技术使计算加速收敛到定常状态;对于非定常流动,使用双时间步长和全隐式离散,采用与求解定常流动相似的多网格方法求解隐式离散方程。文中给出了VKI透平叶栅内的定常流结果和1.5级透平叶栅内的非定常数值结果。 相似文献
6.
Samet Y. Kadioglu Rupert Klein Michael L. Minion 《Journal of computational physics》2008,227(3):2012-2043
A fourth-order numerical method for the zero-Mach-number limit of the equations for compressible flow is presented. The method is formed by discretizing a new auxiliary variable formulation of the conservation equations, which is a variable density analog to the impulse or gauge formulation of the incompressible Euler equations. An auxiliary variable projection method is applied to this formulation, and accuracy is achieved by combining a fourth-order finite-volume spatial discretization with a fourth-order temporal scheme based on spectral deferred corrections. Numerical results are included which demonstrate fourth-order spatial and temporal accuracy for non-trivial flows in simple geometries. 相似文献
7.
Floraine Cordier Pierre Degond Anela Kumbaro 《Journal of computational physics》2012,231(17):5685-5704
We present an Asymptotic-Preserving ‘all-speed’ scheme for the simulation of compressible flows valid at all Mach-numbers ranging from very small to order unity. The scheme is based on a semi-implicit discretization which treats the acoustic part implicitly and the convective and diffusive parts explicitly. This discretization, which is the key to the Asymptotic-Preserving property, provides a consistent approximation of both the hyperbolic compressible regime and the elliptic incompressible regime. The divergence-free condition on the velocity in the incompressible regime is respected, and an the pressure is computed via an elliptic equation resulting from a suitable combination of the momentum and energy equations. The implicit treatment of the acoustic part allows the time-step to be independent of the Mach number. The scheme is conservative and applies to steady or unsteady flows and to general equations of state. One and two-dimensional numerical results provide a validation of the Asymptotic-Preserving ‘all-speed’ properties. 相似文献
8.
Pierre-Henri Maire 《Journal of computational physics》2009,228(7):2391-2425
We present a high-order cell-centered Lagrangian scheme for solving the two-dimensional gas dynamics equations on unstructured meshes. A node-based discretization of the numerical fluxes for the physical conservation laws allows to derive a scheme that is compatible with the geometric conservation law (GCL). Fluxes are computed using a nodal solver which can be viewed as a two-dimensional extension of an approximate Riemann solver. The first-order scheme is conservative for momentum and total energy, and satisfies a local entropy inequality in its semi-discrete form. The two-dimensional high-order extension is constructed employing the generalized Riemann problem (GRP) in the acoustic approximation. Many numerical tests are presented in order to assess this new scheme. The results obtained for various representative configurations of one and two-dimensional compressible fluid flows show the robustness and the accuracy of our new scheme. 相似文献
9.
F. Bassi C. De Bartolo R. Hartmann A. Nigro 《Journal of computational physics》2009,228(11):3996-4011
In this work we extend the high-order discontinuous Galerkin (DG) finite element method to inviscid low Mach number flows. The method here presented is designed to improve the accuracy and efficiency of the solution at low Mach numbers using both explicit and implicit schemes for the temporal discretization of the compressible Euler equations. The algorithm is based on a classical preconditioning technique that in general entails modifying both the instationary term of the governing equations and the dissipative term of the numerical flux function (full preconditioning approach). In the paper we show that full preconditioning is beneficial for explicit time integration while the implicit scheme turns out to be efficient and accurate using just the modified numerical flux function. Thus the implicit scheme could also be used for time accurate computations. The performance of the method is demonstrated by solving an inviscid flow past a NACA0012 airfoil at different low Mach numbers using various degrees of polynomial approximations. Computations with and without preconditioning are performed on different grid topologies to analyze the influence of the spatial discretization on the accuracy of the DG solutions at low Mach numbers. 相似文献
10.
Huajun Zhu Songhe Song & Yaming Chen 《advances in applied mathematics and mechanics.》2011,3(6):663-688
In this paper, we develop a multi-symplectic wavelet collocation method for
three-dimensional (3-D) Maxwell's equations. For the multi-symplectic formulation
of the equations, wavelet collocation method based on autocorrelation functions
is applied for spatial discretization and appropriate symplectic scheme is employed
for time integration. Theoretical analysis shows that the proposed method is
multi-symplectic, unconditionally stable and energy-preserving under periodic
boundary conditions. The numerical dispersion relation is investigated. Combined
with splitting scheme, an explicit splitting symplectic wavelet collocation method
is also constructed. Numerical experiments illustrate that the proposed methods are
efficient, have high spatial accuracy and can preserve energy conservation laws exactly. 相似文献
11.
The compact implicit integration factor (cIIF) method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear Schro¨dinger (NLS) equation and the complex Ginzburg–Landau (GL) equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient. 相似文献
12.
The proposed work concerns the numerical approximations of the shallow-water equations with varying topography. The main objective is to introduce an easy and systematic technique to enforce the well-balance property and to make the scheme able to deal with dry areas. To access such an issue, the derived numerical method is obtained by involving the free surface instead of the water height and this produces the scheme well-balanced independently from the numerical flux function associated with the homogeneous problem. As a consequence, we obtain an easy well-balanced scheme which preserves non-negative water height. When compared with the well-known hydrostatic reconstruction, the presented topography discretization does not involve any max function known to introduce some numerical errors as soon as the topography admits very strong variations or discontinuities. A second-order MUSCL accurate reconstruction is adopted. The proposed hydrostatic upwind scheme is next extended for considering 2D simulations performed over unstructured meshes. Several 1D and 2D numerical experiments are performed to exhibit the relevance of the scheme. 相似文献
13.
研究二维散乱点集上数值求解非线性扩散方程的有限方向差分方法。利用五个邻点信息构造具有最小模板的离散格式,并且离散系数具有显式表达式。另外,利用五点公式获得了间断问题物质界面的离散格式,该格式对界面流的计算具有近似二阶精度。不同计算区域及不同类型的离散点集上的计算结果验证了方法的有效性。 相似文献
14.
15.
16.
In this paper, a semi-implicit finite element method is presented for the coupled Cahn–Hilliard and Navier–Stokes equations with the generalized Navier boundary condition for the moving contact line problems. In our method, the system is solved in a decoupled way. For the Cahn–Hilliard equations, a convex splitting scheme is used along with a P1-P1 finite element discretization. The scheme is unconditionally stable. A linearized semi-implicit P2-P0 mixed finite element method is employed to solve the Navier–Stokes equations. With our method, the generalized Navier boundary condition is extended to handle the moving contact line problems with complex boundary in a very natural way. The efficiency and capacity of the present method are well demonstrated with several numerical examples. 相似文献
17.
Grid convergence studies for subsonic and transonic flows over airfoils are presented in order to compare the accuracy of several spatial discretizations for the compressible Navier–Stokes equations. The discretizations include the following schemes for the inviscid fluxes: (1) second-order-accurate centered differences with third-order matrix numerical dissipation, (2) the second-order convective upstream split pressure scheme (CUSP), (3) third-order upwind-biased differencing with Roe's flux-difference splitting, and (4) fourth-order centered differences with third-order matrix numerical dissipation. The first three are combined with second-order differencing for the grid metrics and viscous terms. The fourth discretization uses fourth-order differencing for the grid metrics and viscous terms, as well as higher-order approximations near boundaries and for the numerical integration used to calculate forces and moments. The results indicate that the discretization using higher-order approximations for all terms is substantially more accurate than the others, producing less than two percent numerical error in lift and drag components on grids with less than 13,000 nodes for subsonic cases and less than 18,000 nodes for transonic cases. Since the cost per grid node of all of the discretizations studied is comparable, the higher-order discretization produces solutions of a given accuracy much more efficiently than the others. 相似文献
18.
In recent years multigrid methods have been proven to be very efficient for solving large systems of linear equations resulting from the discretization of positive definite differential equations by either the finite difference method or theh-version of the finite element method. In this paper an iterative method of the multiple level type is proposed for solving systems of algebraic equations which arise from thep-version of the finite element analysis applied to indefinite problems. A two-levelV-cycle algorithm has been implemented and studied with a Gauss–Seidel iterative scheme used as a smoother. The convergence of the method has been investigated, and numerical results for a number of numerical examples are presented. 相似文献
19.
John W. Barrett Mark Galassi Warner A. Miller Rafael D. Sorkin Philip A. Tuckey Ruth M. Williams 《International Journal of Theoretical Physics》1997,36(4):815-839
The role of Regge calculus as a tool for numerical relativity is discussed, and a parallelizable implicit evolution scheme
described. Because of the structure of the Regge equations, it is possible to advance the vertices of a triangulated spacelike
hypersurface in isolation, solving at each vertex a purely local system of implicit equations for the new edge lengths involved.
(In particular, equations of global “elliptic type” do not arise.) Consequently, there exists a parallel evolution scheme
which divides the vertices into families of nonadjacent elements and advances all the vertices of a family simultaneously.
The relation between the structure of the equations of motion and the Bianchi identities is also considered. The method is
illustrated by a preliminary application to a 600-cell Friedmann cosmology. The parallelizable evolution algorithm described
in this paper should enable Regge calculus to be a viable discretization technique in numerical relativity. 相似文献
20.
Margarete O. Domingues Sônia M. Gomes Olivier Roussel Kai Schneider 《Journal of computational physics》2008,227(8):3758-3780
We present a fully adaptive numerical scheme for evolutionary PDEs in Cartesian geometry based on a second-order finite volume discretization. A multiresolution strategy allows local grid refinement while controlling the approximation error in space. For time discretization we use an explicit Runge–Kutta scheme of second-order with a scale-dependent time step. On the finest scale the size of the time step is imposed by the stability condition of the explicit scheme. On larger scales, the time step can be increased without violating the stability requirement of the explicit scheme. The implementation uses a dynamic tree data structure. Numerical validations for test problems in one space dimension demonstrate the efficiency and accuracy of the local time-stepping scheme with respect to both multiresolution scheme with global time stepping and finite volume scheme on a regular grid. Fully adaptive three-dimensional computations for reaction–diffusion equations illustrate the additional speed-up of the local time stepping for a thermo-diffusive flame instability. 相似文献