An operator splitting algorithm for the three-dimensional advection–diffusion equation

Operator splitting algorithms are frequently used for solving the advection–diffusion equation, especially to deal with advection dominated transport problems. In this paper an operator splitting algorithm for the three-dimensional advection–diffusion equation is presented. The algorithm represents a second-order-accurate adaptation of the Holly and Preissmann scheme for three-dimensional problems. The governing equation is split into an advection equation and a diffusion equation, and they are solved by a backward method of characteristics and a finite element method, respectively. The Hermite interpolation function is used for interpolation of concentration in the advection step. The spatial gradients of concentration in the Hermite interpolation are obtained by solving equations for concentration gradients in the advection step. To make the composite algorithm efficient, only three equations for first-order concentration derivatives are solved in the diffusion step of computation. The higher-order spatial concentration gradients, necessary to advance the solution in a computational cycle, are obtained by numerical differentiations based on the available information. The simulation characteristics and accuracy of the proposed algorithm are demonstrated by several advection dominated transport problems. © 1998 John Wiley & Sons, Ltd.     

三维对流问题的拟协调六面体单元解法

寻找一种高精度的空间单元插值模式是数值求解三维对流问题的关键。在前人研究的基础上,探讨了一种任意空间六面体的拟协凋单元,保证节点上的物理量函数及其一阶导数连续。算例表明,该方法具有良好的计算稳定性和低数值阻尼的优点,且计算工作量大大小于协调单元法,有利于推广应用于对流扩散方程的数值求解。     

Investigation of use of reach-back characteristics method for 2D dispersion equation

The use of the Holly-Preissmann two-point scheme has been very popular for the calculation of the dispersion equation. The key to this scheme is to use the characteristics method incorporating the Hermite cubic interpolation technique to approximate the trajectory foot of the characteristics. This method can avoid the excessive numerical damping and oscillation associated with most finite difference schemes for advection computation. On the basis of the fundamental idea of the Holly-Preissmann two-point scheme, a new technique is introduced herein for the computation of the two-dimensional dispersion equation. This new scheme allows the characteristics projecting back several time steps to fall on the spatial or temporal axis, while the characteristics foot is still solved by the Holly-Preissmann two-point method. The diffusion portion of the dispersion equation is solved by the commonly used Crank-Nicholson method. The calculation for these two processes consisting of advection and diffusion is carried out separately but consecutively in one time step, a method known as the split operator algorithm. A hypothetical model was constructed to demonstrate the applicability of this new technique for the calculation of the pure advection and dispersion equation in two dimensions.     

A WAVE EQUATION MODEL TO SOLVE THE MULTIDIMENSIONAL TRANSPORT EQUATION

The wave equation model, originally developed to solve the advection–diffusion equation, is extended to the multidimensional transport equation in which the advection velocities vary in space and time. The size of the advection term with respect to the diffusion term is arbitrary. An operator-splitting method is adopted to solve the transport equation. The advection and diffusion equations are solved separate ly at each time step. During the advection phase the advection equation is solved using the wave equation model. Consistency of the first-order advection equation and the second-order wave equation is established. A finite element method with mass lumping is employed to calculate the three-dimensional advection of both a Gaussian cylinder and sphere in both translational and rotational flow fields. The numerical solutions are accurate in comparison with the exact solutions. The numerical results indicate that (i) the wave equation model introduces minimal numerical oscillation, (ii) mass lumping reduces the computational costs and does not significantly degrade the numerical solutions and (iii) the solution accuracy is relatively independent of the Courant number provided that a stability constraint is satisfied. © 1997 by John Wiley & Sons, Ltd.     

On the use of the reach-back characteristics method for calculation of dispersion

The Holly-Preissmann two-point finite difference scheme (HP method) has been popularly used for solving the advection equation. The key idea of this scheme is to solve the dependent variable (i.e. the concentration for the pollutant transport problem) by the method of characteristics with the use of cubic interpolation on the spatial axis. The interpolating polynomials of higher order are constructed by use of the dependent variable and its derivatives at two adjacent grid points. In this paper a new interpolating technique is introduced for incorporation with the Holly-Preissmann two-point method. The new method is denoted herein as the Holly-Preissmann reach-back method (HPRB) and allows the characteristics to project back several time steps beyond the present time level. Through stability analyses it has been observed that the increase of the reach-back time step numbers for the characteristics indeed reduces the numerical damping and dispersive phenomena. A schematic model has been constructed to demonstrate the merits of this new technique for the calculation of the pure advection and dispersion equations. Numerical experiments and comparisons with analytical solutions which support and demonstrate this new technique are presented.     

An accurate moving grid Eulerian Lagrangian localized adjoint method for solving the one‐dimensional variable‐coefficient ADE

An accurate finite‐volume Eulerian Lagrangian localized adjoint method (ELLAM) is presented for solving the one‐dimensional variable coefficients advection dispersion equation that governs transport of solute in porous medium. The method uses a moving grid to define the solution and test functions. Consequently, the need for spatial interpolation, or equivalently numerical integration, which is a major issue in conventional ELLAM formulations, is avoided. After reviewing the one‐dimensional method of ELLAM, we present our strategy and detailed calculations for both saturated and unsaturated porous medium. Numerical results for a constant‐coefficient problem and a variable‐coefficient problem are very close to analytical and fine‐grid solutions, respectively. The strength of the developed method is shown for a large range of CFL and grid Peclet numbers. Copyright 2004 John Wiley & Sons, Ltd.     

A simple high‐resolution advection scheme

A simple, robust, mass‐conserving numerical scheme for solving the linear advection equation is described. The scheme can estimate peak solution values accurately even in regions where spatial gradients are high. Such situations present a severe challenge to classical numerical algorithms. Attention is restricted to the case of pure advection in one and two dimensions since this is where past numerical problems have arisen. The authors' scheme is of the Godunov type and is second‐order in space and time. The required cell interface fluxes are obtained by MUSCL interpolation and the exact solution of a degenerate Riemann problem. Second‐order accuracy in time is achieved via a Runge–Kutta predictor–corrector sequence. The scheme is explicit and expressed in finite volume form for ease of implementation on a boundary‐conforming grid. Benchmark test problems in one and two dimensions are used to illustrate the high‐spatial accuracy of the method and its applicability to non‐uniform grids. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical modelling of shear-dependent mass transfer in large arteries

A numerical scheme for the simulation of blood flow and transport processes in large arteries is presented. Blood flow is described by the unsteady 3D incompressible Navier–Stokes equations for Newtonian fluids; solute transport is modelled by the advection–diffusion equation. The resistance of the arterial wall to transmural transport is described by a shear-dependent wall permeability model. The finite element formulation of the Navier–Stokes equations is based on an operator-splitting method and implicit time discretization. The streamline upwind/Petrov–Galerkin (SUPG) method is applied for stabilization of the advective terms in the transport equation and in the flow equations. A numerical simulation is carried out for pulsatile mass transport in a 3D arterial bend to demonstrate the influence of arterial flow patterns on wall permeability characteristics and transmural mass transfer. The main result is a substantial wall flux reduction at the inner side of the curved region. © 1997 John Wiley & Sons, Ltd.     

The perturbation finite element method for solving problems with nonlinear materials

The perturbation method is one of the effective methods for so-lving problems in nonlinear continuum mechanics.It has been de-veloped on the basis of the linear analytical solutions for the o-riginal problems.If a simple analytical solution cannot be ob-tained.we would encounter difficulties in applying this method tosolving certain complicated nonlinear problems.The finite ele-ment method appears to be in its turn a very useful means for sol-ving nonlinear problems,but generally it takes too much time incomputation.In the present paper a mixed approach,namely,theperturbation finite element method,is introduced,which incorpo-rates the advantages of the two above-mentioned methods and enablesus to solve more complicated nonlinear problems with great savingin computing time.Problems in the elastoplastic region have been discussed anda numerical solution for a plate with a central hole under tensionis given in this paper.     

Dispersion error reduction for acoustic problems using the smoothed finite element method (SFEM)

The smoothed finite element method (SFEM), which was recently introduced for solving the mechanics and acoustic problems, uses the gradient smoothing technique to operate over the cell‐based smoothing domains. On the basis of the previous work, this paper reports a detailed analysis on the numerical dispersion error in solving two‐dimensional acoustic problems governed by the Helmholtz equation using the SFEM, in comparison with the standard finite element method. Owing to the proper softening effects provided naturally by the cell‐based gradient smoothing operations, the SFEM model behaves much softer than the standard finite element method model. Therefore, the SFEM can significantly reduce the dispersion error in the numerical solution. Results of both theoretical and numerical experiments will support these important findings. It is shown clearly that the SFEM suits ideally well for solving acoustic problems, because of the crucial effectiveness in reducing the dispersion error. Copyright © 2015 John Wiley & Sons, Ltd.     

A higher-order eulerian scheme for coupled advection-diffusion transport

A new accurate high-order numerical method is presented for the coupled transport of a passive scalar (concentration) by advection and diffusion. Following the method of characteristics, the pure advection problem is first investigated. Interpolation of the concentration and its first derivative at the foot of the characteristic is carried out with a fifth-degree polynomial. The latter is constructed by using as information the concentration and its first and second derivatives at computational points on current time level t in Eulerian co-ordinates. The first derivative involved in the polynomial is transported by advection along the characteristic towards time level t + Δt in the same way as is the concentration itself. Second derivatives are obtained at the new time level t + Δt by solving a system of linear equations defined only by the concentrations and their derivatives at grid nodes, with the assumption that the third-order derivatives are continuous. The approximation of the method is of sixth order. The results are extended to coupled transport by advection and diffusion. Diffusion of the concentration takes place in parallel with advection along the characteristic. The applicability and precision of the method are demonstrated for the case of a Gaussian initial distribution of concentrations as well as for the case of a steep advancing concentration front. The results of the simulations are compared with analytical solutions and some existing methods.     

Particle transport method for convection problems with reaction and diffusion

The paper is devoted to the further development of the particle transport method for the convection problems with diffusion and reaction. Here, the particle transport method for a convection–reaction problem is combined with an Eulerian finite‐element method for diffusion in the framework of the operator‐splitting approach. The technique possesses a special spatial adaptivity to resolve solution singularities possible due to convection and reaction terms. A monotone projection technique is used to transfer the solution of the convection–reaction subproblem from a moving set of particles onto a fixed grid to initialize the diffusion subproblem. The proposed approach exhibits good mass conservation and works with structured and unstructured meshes. The performance of the presented algorithm is tested on one‐ and two‐dimensional benchmark problems. The numerical results confirm that the method demonstrates good accuracy for the convection‐dominated as well as for convection–diffusion problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Conservative monotone streamline upwind formulation using simplex elements

A streamline upwind formulation is presented for the treatment of the advection terms in the general transport equation. The formulation is monotone and conservative and is based on the discontinuous nature of the advection mechanism. The results of there benchmark test cases for the full range of flow Peclet numbers are presented. The new formulation is shown to accurately model the advection phenomenon with significantly smaller numerical diffusion than the existing methods. The results are also free of all spatial oscillations. Considerable savings in computer storage and execution time have been achieved by employing the three-noded triangular element for which exact integrations exist. The formulation is straightforward and can be readily incorporated into any finite element code using the conventional Galerkin approach.     

饱和-非饱和土壤中污染物运移过程的数值模拟

本文提出了一个模拟饱和 非饱和土壤中溶和污染物运移过程的数值模型．模拟的控制污染物运移的物理 化学现象包括：对流，机械逸散，分子弥散，吸附，蜕变，不动水效应．发展了一个修正的特征线Ｇａｌｅｒｋｉｎ方法以离散污染物运移过程的控制方程并导出了一个用于有限元方程求解的显式算法．数值例题结果表明所提出模型和算法的功能     

A non‐iterative implicit algorithm for the solution of advection–diffusion equation on a sphere

A numerical algorithm for the solution of advection–diffusion equation on the surface of a sphere is suggested. The velocity field on a sphere is assumed to be known and non‐divergent. The discretization of advection–diffusion equation in space is carried out with the help of the finite volume method, and the Gauss theorem is applied to each grid cell. For the discretization in time, the symmetrized double‐cycle componentwise splitting method and the Crank–Nicolson scheme are used. The numerical scheme is of second order approximation in space and time, correctly describes the balance of mass of substance in the forced and dissipative discrete system and is unconditionally stable. In the absence of external forcing and dissipation, the total mass and L₂‐norm of solution of discrete system is conserved in time. The one‐dimensional periodic problems arising at splitting in the longitudinal direction are solved with Sherman–Morrison's formula and Thomas's algorithm. The one‐dimensional problems arising at splitting in the latitudinal direction are solved by the bordering method that requires a prior determination of the solution at the poles. The resulting linear systems have tridiagonal matrices and are solved by Thomas's algorithm. The suggested method is direct (without iterations) and rapid in realization. It can also be applied to linear and nonlinear diffusion problems, some elliptic problems and adjoint advection–diffusion problems on a sphere. Copyright © 2015 John Wiley & Sons, Ltd.     

SOLUTION OF THE ADVECTION–DIFFUSION EQUATION USING A COMBINATION OF DISCONTINUOUS AND MIXED FINITE ELEMENTS

When transport is advection-dominated, classical numerical methods introduce excessive artificial diffusion and spurious oscillations. Special methods are required to overcome these phenomena. To solve the advection‒diffusion equation, a numerical method is developed using a discontinuous finite element method for the discretization of the advective terms. At the discontinuities of the approximate solution, numerical advective fluxes are calculated using one-dimensional approximate Riemann solvers. The method is stabilized with a multidimensional slope limiter which introduces small amounts of numerical diffusion when sharp concentration fronts occur. In addition, the diffusive term is discretized using a mixed hybrid finite element method. With this approach, numerical oscillations are completely avoided for a full range of cell Peclet numbers. The combination of discontinuous and mixed finite elements can be easily applied to 2D and 3D models using various types of elements in regular and irregular meshes. Numerical tests show good agreement with 1D and 2D analytical solutions. This approach is compared at the same time with two different numerical methods, a standard mixed finite method and a finite volume approach with high-resolution upwind terms. Regular and irregular meshes are used for the numerical tests to study the mesh effects on the numerical results. Our data show that in all cases this approach performs well. © 1997 by John Wiley & Sons, Ltd.     

Numerical modeling of cohesive sediments dynamics in estuaries: Part I—Description of the model and simulations in the Po River Estuary

The present work contributes to the numerical modeling of complex turbulent multiphasic fluid flows occurring in estuarine channels. This research finds its motivation in the increasing need for efficient management of estuaries by taking into account the complex turbulent stratified flows encountered in estuaries and costal zones. A time‐dependent, 3D finite element model of suspended sediment transport taking into account the effects of cohesiveness between sediments is presented. The model estuary is the forced time‐dependent winds, time elevation at open boundaries and river discharge. To cope with the stiffness problems a decoupling method is employed to solve the shallow‐water equations of mass conservation, momentum and suspended sediment transport with the conventional hydrostatic pressure. The decoupling method partitions a time step into three subcycles according to the physical phenomena. In the first sub‐cycle the pure hydrodynamics including the k–ε turbulence model is solved, followed by the advection–diffusion equations for pollutants (salinity, temperature, suspended sediment concentration, (SSC)), and finally the bed evolution is solved. The model uses a mass‐preserving method based on the so‐called Raviart–Thomas finite element on the unstructured mesh in the horizontal plane, while the multi‐layers system is adopted in vertical with the conventional conforming finite element method, with the advantage that the lowermost and uppermost layers of variable height allow a faithful representation of the time‐varying bed and free surface, respectively. The model has been applied to investigate the SSC and seabed evolution in Po River Estuary (PRE) in Italy. The computed results mimic the field data well. Copyright © 2007 John Wiley & Sons, Ltd.     

Augmented upwind numerical schemes for the groundwater transport advection‐dispersion equation with local operators

When solute transport is advection‐dominated, the advection‐dispersion equation approximates to a hyperbolic‐type partial differential equation, and finite difference and finite element numerical approximation methods become prone to artificial oscillations. The upwind scheme serves to correct these responses to produce a more realistic solution. The upwind scheme is reviewed and then applied to the advection‐dispersion equation with local operators for the first‐order upwinding numerical approximation scheme. The traditional explicit and implicit schemes, as well as the Crank‐Nicolson scheme, are developed and analyzed for numerical stability to form a comparison base. Two new numerical approximation schemes are then proposed, namely, upwind–Crank‐Nicolson scheme, where only for the advection term is applied, and weighted upwind‐downwind scheme. These newly developed schemes are analyzed for numerical stability and compared to the traditional schemes. It was found that an upwind–Crank‐Nicolson scheme is appropriate if the Crank‐Nicolson scheme is only applied to the advection term of the advection‐dispersion equation. Furthermore, the proposed explicit weighted upwind‐downwind finite difference numerical scheme is an improvement on the traditional explicit first‐order upwind scheme, whereas the implicit weighted first‐order upwind‐downwind finite difference numerical scheme is stable under all assumptions when the appropriate weighting factor (θ) is assigned.     

Solution of the Navier–Stokes equations in velocity–vorticity form using a Eulerian–Lagrangian boundary element method

This paper describes the Eulerian–Lagrangian boundary element model for the solution of incompressible viscous flow problems using velocity–vorticity variables. A Eulerian–Lagrangian boundary element method (ELBEM) is proposed by the combination of the Eulerian–Lagrangian method and the boundary element method (BEM). ELBEM overcomes the limitation of the traditional BEM, which is incapable of dealing with the arbitrary velocity field in advection‐dominated flow problems. The present ELBEM model involves the solution of the vorticity transport equation for vorticity whose solenoidal vorticity components are obtained iteratively by solving velocity Poisson equations involving the velocity and vorticity components. The velocity Poisson equations are solved using a boundary integral scheme and the vorticity transport equation is solved using the ELBEM. Here the results of two‐dimensional Navier–Stokes problems with low–medium Reynolds numbers in a typical cavity flow are presented and compared with a series solution and other numerical models. The ELBEM model has been found to be feasible and satisfactory. Copyright © 2000 John Wiley & Sons, Ltd.     

A least-squares procedure for the solution of transport problems

In this paper a least-squares formulation associated with a conjugate gradient algorithm is proposed for the solution of transport problems. In this procedure the advection–diffusion equation is first discretized in time using an implicit scheme. At each time step the resulting partial differential equation is replaced by an optimal control problem. This minimization problem involves the minimization of a functional defined via a state equation. This functional is chosen in order to force the numerical solution of the advection–diffusion equation to be equal to the hyperbolic advective part of this equation. The effectiveness of the method is shown through a one-dimensional example involving advective and diffusive transport. No oscillation and high accuracy have been obtained for the entire range of Peclet numbers with a Courant number well in excess of unity.