Numerical Solution of Advection Diffusion Equation Using Semi-Discretization Scheme

Numerical diffusion and oscillatory behavior characteristics are averted applying numerical solutions of advection-diffusion equation are themselves immensely sophisticated. In this paper, two numerical methods have been used to solve the advection diffusion equation. We use an explicit finite difference scheme for the advection diffusion equation and semi-discretization on the spatial variable for advection-diffusion equation yields a system of ordinary differential equations solved by Euler’s method. Numerical assessment has been executed with specified initial and boundary conditions, for which the exact solution is known. We compare the solutions of the advection diffusion equation as well as error analysis for both schemes.

Finite element multigrid method for multi-term time fractional advection diffusion equations

In this paper,a class of multi-term time fractional advection diffusion equations(MTFADEs)is considered.By finite difference method in temporal direction and finite element method in spatial direction,two fully discrete schemes of MTFADEs with different definitions on multi-term time fractional derivative are obtained.The stability and convergence of these numerical schemes are discussed.Next,a V-cycle multigrid method is proposed to solve the resulting linear systems.The convergence of the multigrid method is investigated.Finally,some numerical examples are given for verification of our theoretical analysis.

A BFG model for calculation of tidal current and diffusion of pollutants in nearshore areas

This study presents a boundary-fitted grid (BFG) numerical model with an aim to simulate the tidal currents and diffusion of pollutants in complicated nearshore areas. To suit the general model to any curvilinear grids, generalized 2-D shallow sea dynamic equations and the advection diffusion equation are derived in curvilinear coordinates, and the contravariant components of the velocity vector are adopted for easily realizing boundary conditions and making the equations conservational. As the generalized equations are not limited by a speCific coordinate transformation. a self-adaptive grid generation method is then proposed conveniently to generate a boundary-fitted and varying SPacing grid.The calculation in the Yangpu Bay and the Xinying Bay shows that this is an effective model for calculating tidal currents and diffusion of pollutants in the more complicated nearshore areas.

A Numerical Simulation of Air Flow in the Human Respiratory System Based on Lung Model

The lung is an important organ that takes part in the gas exchange process. In the study of gas transport and exchange in the human respiratory system, the complicated process of advection and diffusion (AD) in airways of human lungs is considered. The basis of a lumped parameter model or a transport equation is modeled during the inspiration process, when oxygen enters into the human lung channel. The quantitative measurements of oxygen are detached and the model equation is solved numerically by explicit finite difference schemes. Numerical simulations were made for natural breathing conditions or normal breathing conditions. The respiratory flow results for the resting conditions are found strongly dependent on the AD effect with some contribution of the unsteadiness effect. The contour of the flow rate region is labeled and AD effects are compared with the variation of small intervals of time for a constant velocity when breathing is interrupted for a negligible moment.

A Novel Method for Solving Time-Dependent 2D Advection-Diffusion-Reaction Equations to Model Transfer in Nonlinear Anisotropic Media

This paper presents a new numerical technique for solving initial and bound-ary value problems with unsteady strongly nonlinear advection diffusion reaction(ADR)equations.The method is based on the use of the radial basis functions(RBF)for the approximation space of the solution.The Crank-Nicolson scheme is used for approximation in time.This results in a sequence of stationary nonlinear ADR equations.The equations are solved sequentially at each time step using the proposed semi-analytical technique based on the RBFs.The approximate solution is sought in the form of the analytical expansion over basis functions and contains free parameters.The basis functions are constructed in such a way that the expansion satisfies the boundary conditions of the problem for any choice of the free parameters.The free parameters are determined by substitution of the expansion in the equation and collocation in the solution domain.In the case of a nonlinear equation,we use the well-known procedure of quasilinearization.This transforms the original equation into a sequence of the linear ones on each time layer.The numerical examples confirm the high accuracy and robustness of the proposed numerical scheme.

AN OPERATOR-SPLITTING ALGORITHM FOR ADVECTION-DIFFUSION-REACTION EQUATION

An operator-splitting algorithm for three-dimensional advection-diffusion-reaction equation is presented.The method of characteristics is adopted for the pure advection operator, the explicit difference scheme is used for diffusion,and a prediction-correction scheme is em- ployed for reaction.The condition for stability of the algorithm is analysed.Severall inear and nonlinear examples are illustrated to test the convergence and accuracy of the numerical proce- dure,and satisfactory agreements between computed and analytical solutions are achieved.Due to its simplicity,stability,and validity for both one-and two-dimensional problems,the success- ful algorithm can be used to numerical simulations of viscous fluid flows,the transport of pollu- tants and sedimentations in reservoirs,lakes,rivers,estuaries and other environments,cooling- problems in heat or nuclear power plants,etc.

THREE HIGH-ORDER SPLITTING SCHEMES FOR 3D TRANSPORT EQUATION

Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by these two schemes, only three computational grid points were needed in each direction but the accuracy reaches the spatial fourth-order. The third scheme proposed is based on the classical ADI scheme and the accuracy of the advection term of it can reach the spatial fourth-order. Finally, two typical numerical experiments show that the solutions of these three schemes compare well with that given by the analytical solution when the Peclet number is not bigger than 5.

CONVERGENCE RESULTS FOR NON-OVERLAP SCHWARZ WAVEFORM RELAXATION ALGORITHM WITH CHANGING TRANSMISSION CONDITIONS

In this paper,we establish a new algorithm to the non-overlapping Schwarz domain decomposition methods with changing transmission conditions for solving one dimensional advection reaction diffusion problem.More precisely,we first describe the new algorithm and prove the convergence results under several natural assumptions on the sequences of parameters which determine the transmission conditions.Then we give a simple method to estimate the new value of parameters in each iteration.The interesting advantage of our method is that one may update the better parameters in each iteration to save the computational cost for optimizing the parameters after many steps.Finally some numerical experiments are performed to show the behavior of the convergence rate for the new method.

A NEW APPROACH TO THE NUMERICAL SIMULATION OF MASS TRANSPORT

A refined numerical method, based upon time-line interpolation, for the simulation of advection and diffusion has been tentatively explored. A complete set of temporal reachback numerical scheme in applying the method of characteristics has been derived, and the favorable accuracy of the method demonstrated. The use of interpolations in time, rather than the more widely used interpolations in space, demonstrates that it generates a much smaller numerical error.