定常不可压Navier—Stokes方程迎风格式

本文利用齐次定解条件对定常不可压Navier—Stokes方程的非线性项进行处理,给出了相应的一种迎风Galerkin有限元算法;针对这种迎风Galerkin有限元算法,在迎风参数满足一定条件下,利用其三项式具有的一些很好性质,更简单地证明了该问题解的存在唯一性。     

高维抛物方程迎风区域分裂差分方法

结合迎风方法和区域分裂思想,采用一阶迎风、二阶修正迎风法逼近高维抛物方程的对流项.内边界处和子区域分别对应区域分裂显隐格式;并运用极值原理和嵌入定理给出了收敛性分析,最后给出数值试验,说明其实际意义.     

数值求解Euler方程的UCGVC差分格式的注记

１．引言 近年来高精度差分格式的研究引起国内外的普遍重视,目的是更准确地模拟复杂流场的流动．众所周知,传统的二阶ＴＶＤ类格式虽然能较好地捕捉激波,但却存在局部极值点降阶的问题,而且由于一些格式的数值粘性过大,当用该格式计算粘性流特别是高雷诺数问题时,格式本身的数值粘性可能掩盖了流场的物理粘性,从而降低了格式对边界层的分辨率,因而无法正确计算热流值。文献［３］指出,采用高精度格式可适当放松对网格雷诺数的要求,因此发展三阶或三阶以上的格式是需要的。近年来,人们已经发展了一些无伪振荡的高阶格式,如ＥＮ…     

THE UPWIND FINITE ELEMENT SCHEME AND MAXIMUM PRINCIPLE FOR NONLINEAR CONVECTION-DIFFUSION PROBLEM

In this paper, a kind of partial upwind finite element scheme is studied for twodimensional nonlinear convection-diffusion problem. Nonlinear convection term approximated by partial upwind finite element method considered over a mesh dual to the triangular grid, whereas the nonlinear diffusion term approximated by Galerkin method. A linearized partial upwind finite element scheme and a higher order accuracy scheme are constructed respectively. It is shown that the numerical solutions of these schemes preserve discrete maximum principle. The convergence and error estimate are also given for both schemes under some assumptions. The numerical results show that these partial upwind finite element scheme are feasible and accurate.     

Smooth volume integral conservation law and method for problems in Lagrangian coordinates

An integral conservation law is derived for smooth volume in Lagrangian coordinates (a comoving frame). A method for approximation of the integral smooth volume conservation law is discussed. An extension technique is suggested for development of smooth volume schemes. For hyperbolic systems, smooth volume upwind and Godunov schemes with monotonic reconstruction are derived. The schemes are applied to equations of gas dynamics and tested for three gas-dynamics shock tube problems. The solutions are monotonic and precise. This article was submitted by the author in English.     

对流扩散方程的四阶紧凑迎风差分格式

§1.引言 流动和传热传质的基本方程均是对流扩散型的.对流扩散方程的高阶紧凑差分格式,作为提高计算可靠性和节省计算量的一条有效途径,已引起相当的重视.作为该领域的一大进展,新近由Dennis推出的对流扩散方程四阶紧凑格式,在二维情形下呈九点式且勿须引入中间变量,只涉及对流扩散量本身,能在较粗网格下获取较为准确的数值结果.从本质上说,该格式系指数型四阶紧凑格式的多项式型翻版.它与指数型紧凑格     

一类非线性离散系统的指数二分性及其在数值分析和计算中的应用

本文中我们考虑一类二阶非线性常微分方程的边值问题的迎风差分格式.我们运用奇异摄动方法构造了该迎风差分方程解的渐近近似,并利用指数二分性理论证明了有一个低阶方程其解是该迎风方程式的在边界外的一个良好近似.我们还构造了校正项,使校正项与低阶方程的解之和是一个渐近近似.最后一些数值例子用于显示本文方法的应用.     

Minimal dissipation hybrid bicompact schemes for hyperbolic equations

New monotonicity-preserving hybrid schemes are proposed for multidimensional hyperbolic equations. They are convex combinations of high-order accurate central bicompact schemes and upwind schemes of first-order accuracy in time and space. The weighting coefficients in these combinations depend on the local difference between the solutions produced by the high- and low-order accurate schemes at the current space-time point. The bicompact schemes are third-order accurate in time, while having the fourth order of accuracy and the first difference order in space. At every time level, they can be solved by marching in each spatial variable without using spatial splitting. The upwind schemes have minimal dissipation among all monotone schemes constructed on a minimum space-time stencil. The hybrid schemes constructed has been successfully tested as applied to a number of two-dimensional gas dynamics benchmark problems.     

Error analysis of upwind‐discretizations for the steady‐state incompressible Navier–Stokes equations

Within the framework of finite element methods, the paper investigates a general approximation technique for the nonlinear convective term of the Navier–Stokes equations. The approach is based on an upwind method of finite volume type. It is proved that the discrete convective term satisfies a well‐known collection of sufficient conditions for convergence of the finite element solution. This revised version was published online in June 2006 with corrections to the Cover Date.     

An efficient explicit/implicit domain decomposition method for convection‐diffusion equations

A nonoverlapping domain decomposition method for some time‐dependent convection‐diffusion equations is presented. It combines predictor‐corrector technique, modified upwind differences with explicit/implicit coupling to provide intrinsic parallelism, and unconditional stability while improving the accuracy. Both rigorous mathematical analysis and numerical experiments are carried out to illustrate the stability, accuracy, and parallelism. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Modified Upwind Taylor Finite Element Schemes for 1-D Conservation Laws I. A Basic Idea

In this paper, first, modified upwind finite element schemes are presented for two-point value problem, and then a class of modified upwind Taylor finite element schemes are derived for one dimensional linear hyperbolic equation. The main point of the paper is how to consider the upwind property to construct base functions to make the schemes obtained be MmB (or TVD). Numerical experiments are given to show that the method is efficient to solve the discontinuous solutions.     

激波捕捉差分方法研究

在迎风型格式和矢通量分裂技术的基础之上,对捕捉激波方法进行一种新的尝试．该方法首先对原始格式在特征方向上进行投影,然后用限制器对这些特征分量的变化幅值进行限制以抑止非物理波动,最后再把它转换成守恒形式,得到了基本上无振荡的激波捕捉格式．用该方法对两种迎风显示格式(二阶和三阶)和3种迎风紧致格式(三阶、五阶和七阶)进行处理,并在一维和二维的情况下进行了应用测试．通过与高阶WENO、MP、Compact-WENO等格式的比较,表明该方法在光滑捕捉激波的前提下仍有较高精度和分辨率．     

Uniformly convergent finite difference methods for singularly perturbed problems with turning points

In this work one-dimensional singular perturbation problems with turning points are considered. To resolve these problems numerically we consider a family of finite difference schemes, which includes classical methods in literature, such as the upwind method, the Samarskii method and exponential fitting type methods. Once the uniform convergence of the upwind method on irregular meshes has been established, the same property is easily shown on all the elements of the family.Work supported by a grant of the Diputación General de Aragón.     

An upwind center difference parallel method and numerical analysis for the displacement problem with moving boundary

A nonlinear coupled mathematical system of two‐phase seepage flow displacement is discussed in this paper including an elliptic equation for the pressure and a convection‐dominated diffusion equation for the saturation. In fact, the boundary of an underground region where the fluid flows through is nonstationary. So a moving boundary should be considered. The saturation equation is convection‐dominated, therefore the method of upwind finite difference is introduced for the accurate computation. The upwind approximation could eliminate numerical oscillation and strong stability is shown. Since the computational work of saturation is larger than the pressure, the authors apply a parallel method, decomposing the whole domain into several nonoverlapping subdomains, to simplify the computation. A domain decomposition method coupled with upwind differences is presented for the saturation. The pressure equation is discretized by a five‐point center finite difference method. By using a transformation and defining new inner products and norms, error estimates in l² norm is discussed. Finally, two experimental tests are given to illustrate the efficiency and accuracy of the parallel algorithm.     

Une méthode numérique décentrée pour la résolution d'écoulements diphasiques

This Note focuses on the numerical approximation of two-fluid flow models described by six balance equations. We introduce an original splitting technique especially derived to use the approximate Riemann solvers of the usual gas dynamics and to allow for a straightforward extension to various and detailed exchange source terms. When based on suitable kinetic upwind schemes, the whole scheme preserves the positivity of all the thermodynamic variables under a fairly unrestrictive “CFL like” condition. Several stiff numerical teats, are presented including phase separation, in order to highlight the efficiency of the proposed method.     

Construction of locally conservative fluxes for the SUPG method

We consider the construction of locally conservative fluxes by means of a simple postprocessing technique obtained from the finite element solutions of advection diffusion equations. It is known that a naive calculation of fluxes from these solutions yields nonconservative fluxes. We consider two finite element methods: the usual continuous Galerkin finite element method for solving nondominating advection diffusion equations and the streamline upwind/Petrov‐Galerkin method for solving advection dominated problems. We then describe the postprocessing technique for constructing conservative fluxes from the numerical solutions of the general variational formulation. The postprocessing technique requires solving an auxiliary Neumann boundary value problem on each element independently and it produces a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. We provide a convergence analysis for the postprocessing technique. Performance of the technique and the convergence behavior are demonstrated through numerical examples including a set of test problems for advection diffusion equations, advection dominated equations, and drift‐diffusion equations. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1971–1994, 2015     

Numerical Solution of Hamilton-Jacobi-Bellman Equations by an Upwind Finite Volume Method

In this paper we present a finite volume method for solving Hamilton-Jacobi-Bellman(HJB) equations governing a class of optimal feedback control problems. This method is based on a finite volume discretization in state space coupled with an upwind finite difference technique, and on an implicit backward Euler finite differencing in time, which is absolutely stable. It is shown that the system matrix of the resulting discrete equation is an M-matrix. To show the effectiveness of this approach, numerical experiments on test problems with up to three states and two control variables were performed. The numerical results show that the method yields accurate approximate solutions to both the control and the state variables.     

三维两相渗流驱动问题迎风区域分裂显隐差分法

对三维两相渗流驱动问题提出了两种迎风区域分裂显隐差分格式．压力方程采用了七点差分格式,为了能达到实际并行计算的要求,对饱和度方程采用了迎风区域分裂差分法,内边界处和各子区域分别对应显隐格式．得到了离散l2模收敛性分析,最后给出数值试验,支撑了理论分析结果．     

同位非结构和结构网格摄动有限体积法（PFV）求解二维Navier-Stokes方程组

根据NS方程组的一阶迎风和二阶中心有限体积（UFV和CFV）格式,导出NS方程组迎风和中心摄动有限体积(UPFV和CPFV)格式．该格式通过把控制体界面质量通量摄动展开成网格间距的幂级数,并由守恒方程本身求得幂级数系数而获得．迎风和中心摄动有限体积格式使用了与一阶迎风和二阶中心格式相同的基点数和相同的表达形式,宜于计算机编程．顶盖驱动方腔流和驻点流标量输运的数值实验证明,迎风PFV格式比一阶UFV、二阶CFV格式有更高的精度,更高的分辨率．尤其是良好的鲁棒特性．对顶盖驱动方腔流,在Re数从102到104范围内,亚松弛系数可在0.3～0.8任取,收敛性能良好．     

Front tracking for two-phase flow in reservoir simulation by adaptive mesh

In this article, an algorithm for the numerical approximation of two-phase flow in porous media by adaptive mesh is presented. A convergent and conservative finite volume scheme for an elliptic equation is proposed, together with the finite difference schemes, upwind and MUSCL, for a hyperbolic equation on grids with local refinement. Hence, an IMPES method is applied in an adaptive composite grid to track the front of a moving solution. An object-oriented programmation technique is used. The computational results for different examples illustrate the efficiency of the proposed algorithm. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 673–697, 1997