一维定常型对流占优扩散方程的一类迎风有限体积格式

本文针对一维定常型对流占优扩散方程提出了一类迎风有限体积格式 .该格式对对流项具有二阶精度 ,对扩散项保持一阶精度 ,符合对流占优扩散问题强对流、弱扩散的特点 .     

Numerical solution of 2D and 3D backward facing step flows

The work deals with numerical modelling of flow through 2-dimensional (2D) and 3-dimensional (3D) backward facing step. In laminar case, we apply several higher order upwind and central discretizations and compare numerical results with measurements. The turbulent regime is considered in 2D as well as in 3D and influence of secondary flow is observed. Different modifications of low-Re two equation turbulence models and an explicit algebraic Reynolds stress model (EARSM) are considered. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Combined higher order finite volume and finite element scheme for double porosity and non-linear adsorption of transport problem in porous media

In this work, a dual porosity model of reactive solute transport in porous media is presented. This model consists of a nonlinear-degenerate advection-diffusion equation including equilibrium adsorption to the reaction combined with a first-order equation for the non-equilibrium adsorption interaction processes. The numerical scheme for solving this model involves a combined high order finite volume and finite element scheme for approximation of the advection-diffusion part and relaxation-regularized algorithm for nonlinearity-degeneracy. The combined finite volume-finite element scheme is based on a new formulation developed by Eymard et al. (2010) [10]. This formulation treats the advection and diffusion separately. The advection is approximated by a second-order local maximum principle preserving cell-vertex finite volume scheme that has been recently proposed whereas the diffusion is approximated by a finite element method. The result is a conservative, accurate and very flexible algorithm which allows the use of different mesh types such as unstructured meshes and is able to solve difficult problems. Robustness and accuracy of the method have been evaluated, particularly error analysis and the rate of convergence, by comparing the analytical and numerical solutions for first and second order upwind approaches. We also illustrate the performance of the discretization scheme through a variety of practical numerical examples. The discrete maximum principle has been proved.     

Interaction of the Filter and Spatial Discretisation Operators for Large-Eddy Simulation Using the Approximate Deconvolution Model

The Approximate Deconvolution Model (ADM) for Large–Eddy Simulation is an approach for the computation of turbulent flows. The main idea of ADM is to approximate the unfiltered data in the filtered Navier–Stokes equations by deconvolution of filtered values. This is achieved via repeated filtering. Thereby problems may arise for cases where the order of the spatial discretisation scheme is too low or lower than the order of the filter operation. The interaction of the discretisation and the filtering was investigated for one–dimensional forced Burgers turbulence. For this study, both the order of the filter operator and the order of the spatial discretisation of the derivatives were varied independently. It was found that the use of ADM with low–order discretisation schemes can not be recommended. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulent flow computations on 3D unstructured grids

An edge-based finite element method is presented for the simulation of compressible turbulent flows on unstructured tetrahedral grids. A two equation k–ω turbulence model is employed and the standard Galerkin approach is used for spatial discretisation. Stabilisation of the resulting procedure is achieved by the addition of an appropriate diffusion. An explicit multistage time-stepping scheme is used to advance the solution in time to steady state. The performance of the algorithm is demonstrated for the simulation of a high Reynolds number transonic separated flow over a wing.     

无结构网格有限体积方法及在热对流中的应用研究

首次将无结构三角网格的有限体积方法和压强连接半隐式算法相结合 ,用于求解非平行壁管道中的热对流问题 .并由此分析了化学汽相淀积薄膜生长的均匀性问题 .计算结果对于分析一类管道中热和动量输运现象均有普遍指导意义     

Multidimensional upwind convection schemes for flow in porous media on structured and unstructured quadrilateral grids

Standard reservoir simulation schemes employ first order upwind schemes for approximation of the convective fluxes when multiple phases or components are present. These convective flux schemes rely upon upwind information that is determined according to grid geometry. As a consequence directional diffusion is introduced into the solution that is grid dependent. The effect can be particularly important for cases where the flow is across grid coordinate lines and is known as cross-wind diffusion.Truly higher dimensional upwind schemes that minimize cross-wind diffusion are presented for convective flow approximation on quadrilateral unstructured grids. The schemes are locally conservative and yield improved results that are essentially free of spurious oscillations. The higher dimensional schemes are coupled with full tensor Darcy flux approximations.The benefits of the resulting schemes are demonstrated for classical test problems in reservoir simulation including cases with full tensor permeability fields. The test cases involve a range of structured and unstructured grids with variations in orientation and permeability that lead to flow fields that are poorly resolved by standard simulation methods. The higher dimensional formulations are shown to effectively reduce the numerical cross-wind diffusion effect, leading to improved resolution of concentration and saturation fronts.     

An adaptive refinement method for convection‐diffusion problems

An adaptive finite difference method for singularly perturbed convection‐diffusion problems is presented. The method is introduced using a first‐order upwind scheme and a suitable error estimator based on the first derivatives. To obtain the grid structure needed for the cross stencil a special refinement strategy is considered. To avoid the slave points we change the stencil at the interface points from a cross to a skew one. After the convergence of the refinement algorithm we use a combination of a first order upwind and a second order central schemes to achieve higher order of convergence. Several numerical examples show the efficiency of our treatment. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Finite element approximation of spatially extended predator–prey interactions with the Holling type II functional response

We study the numerical approximation of the solutions of a class of nonlinear reaction–diffusion systems modelling predator–prey interactions, where the local growth of prey is logistic and the predator displays the Holling type II functional response. The fully discrete scheme results from a finite element discretisation in space (with lumped mass) and a semi-implicit discretisation in time. We establish a priori estimates and error bounds for the semi discrete and fully discrete finite element approximations. Numerical results illustrating the theoretical results and spatiotemporal phenomena are presented in one and two space dimensions. The class of problems studied in this paper are real experimental systems where the parameters are associated with real kinetics, expressed in nondimensional form. The theoretical techniques were adapted from a previous study of an idealised reaction–diffusion system (Garvie and Blowey in Eur J Appl Math 16(5):621–646, 2005).     

Numerical simulation of unsteady flows with transient regimes

Self-oscillatory flows in aerodynamics and astrophysics are studied. The two-dimensional compressible gas equations are solved using the implicit Runge-Kutta scheme of the third order with respect to the inviscid terms and of the second order with respect to the viscous terms. An algebraic Cebeci-Smith turbulence model is used. Weakly unsteady and strongly unsteady flow regimes are observed. The former occur in a supersonic flow past a cylinder with a front projection and in the heliosphere. Such flows became stable when the turbulent diffusion is taken into account. The latter flows occur when a supersonic jet meets an obstacle and when such a jet penetrates a cavity. In these flows, the amplitude of oscillations slightly decreases when the turbulent diffusion is taken into account.     

A new combined finite element‐upwind finite volume method for convection‐dominated diffusion problems

In this article, we develop a combined finite element‐weighted upwind finite volume method for convection‐dominated diffusion problems in two dimensions, which discretizes the diffusion term with the standard finite element scheme, and the convection and source terms with the weighted upwind finite volume scheme. The developed method leads to a totally new scheme for convection‐dominated problems, which overcomes numerical oscillation, avoids numerical dispersion, and has high‐order accuracy. Stability analyses of the scheme are given for the problems with constant coefficients. Numerical experiments are presented to illustrate the stability and optimal convergence of our proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 799–818, 2016     

Simulation of planar ionization wave front propagation on an unstructured adaptive grid

The paper deals with the numerical solution of a basic 2D model of the propagation of an ionization wave. The system of equations describing this propagation consists of a coupled set of reaction–diffusion-convection equations and a Poissons equation. The transport equations are solved by a finite volume method on an unstructured triangular adaptive grid. The upwind scheme and the diamond scheme are used for the discretization of the convection and diffusion fluxes, respectively. The Poisson equation is also discretized by the diamond scheme. Numerical results are presented. We deal in more detail with numerical tests of the grid adaptation technique and its influence on the numerical results. An original behavior is observed. The grid refinement is not sufficient to obtain accurate results for this particular phenomenon. Using a second order scheme for convection is necessary.     

Upwind Based Parameter Uniform Convergence Analysis for Two Parametric Parabolic Convection Diffusion Problems by Moving Mesh Methods

A parabolic two parametric convection-diffusion reaction problem is considered for the moving mesh error analysis. The continuous problem is discretized by the first order upwind scheme on a non uniform mesh. A curvature based error monitor function is proposed to generate the layer adapted mesh. It is proved that the numerical solution converges to the exact solution on the mesh obtained by the equidistribution of the proposed monitor function. The convergence is first order accurate. The present analysis generalizes the results obtained in earlier publications [8,9]. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of turbulent magnetic reconnection in the small-scale solar wind

Some observational examples for the possible occurrence of the turbulent magnetic reconnection in the solar wind are found by analysing Helios spacecraft's high resolution data. The phenomena of turbulent magnetic reconnections in small scale solar wind are simulated by introducing a third order accuracy upwind compact difference scheme to the compressible two_dimensional MHD flow. Numerical results verify that the turbulent magnetic reconnection process could occur in small scale interplanetary solar wind, which is a basic feature characterizing the magnetic reconnection in high_magnetic Reynolds number (RM=2 000-10 000) solar wind. The configurations of the magnetic reconnection could evolve from a single X_line to a multiple X-line reconnection, exhibiting a complex picture of the formation, merging and evolution of magnetic islands, and finally the magnetic reconnection would evolve into a low_energy state. Its life_span of evolution is about one hour order of magnitude. Various magnetic and flow signatures are recorded in the numerical test for different evolution stages and along different crossing paths, which could in principle explain and confirm the observational samples from the Helios spacecraft. These results are helpful for revealing the basic physical processes in the solar wind turbulence.     

Second order blended multidimensional upwind residual distribution scheme for steady and unsteady computations

Multidimensional upwind residual distribution (RD) schemes have become an appealing alternative to more widespread finite volume and finite element methods (FEM) for solving compressible fluid flows. The RD approach allows to construct nonlinear second order and non-oscillatory methods at the same time. They are routinely used for steady state calculations of the complex flow problems, e.g., 3D turbulent transonic industrial-type simulations [H. Deconinck, K. Sermeus, R. Abgrall, Status of multidimensional upwind residual distribution schemes and applications in aeronautics, AAIA Paper 2000-2328, AIAA, 2000; K. Sermeus, H. Deconinck, Drag prediction validation of a multi-dimensional upwind solver, CFD-based aircraft drag prediction and reduction, VKI Lecture Series 2003-02, Von Karman Institute for Fluid Dynamics, Chausée do Waterloo 72, B-1640 Rhode Saint Genèse, Belgium, 2003].     

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

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

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

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

High order Hybrid central—WENO finite difference scheme for conservation laws

In this article we present a high resolution hybrid central finite difference—WENO scheme for the solution of conservation laws, in particular, those related to shock–turbulence interaction problems. A sixth order central finite difference scheme is conjugated with a fifth order weighted essentially non-oscillatory WENO scheme in a grid-based adaptive way. High order multi-resolution analysis is used to detect the high gradients regions of the numerical solution in order to capture the shocks with the WENO scheme while the smooth regions are computed with the more efficient and accurate central finite difference scheme. The application of high order filtering to mitigate the dispersion error of central finite difference schemes is also discussed. Numerical experiments with the 1D compressible Euler equations are shown.     

Upwind compact finite difference scheme for time-accurate solution of the incompressible Navier-Stokes equations

This article presents a time-accurate numerical method using high-order accurate compact finite difference scheme for the incompressible Navier-Stokes equations. The method relies on the artificial compressibility formulation, which endows the governing equations a hyperbolic-parabolic nature. The convective terms are discretized with a third-order upwind compact scheme based on flux-difference splitting, and the viscous terms are approximated with a fourth-order central compact scheme. Dual-time stepping is implemented for time-accurate calculation in conjunction with Beam-Warming approximate factorization scheme. The present compact scheme is compared with an established non-compact scheme via analysis in a model equation and numerical tests in four benchmark flow problems. Comparisons demonstrate that the present third-order upwind compact scheme is more accurate than the non-compact scheme while having the same computational cost as the latter.