A Global Spectral Element Model for Poisson Equations and Advective Flow over a Sphere

A global spherical Fourier-Legendre spectral element method is proposed to solve Poisson equations and advective flow over a sphere. In the meridional direction, Legendre polynomials are used and the region is divided into several elements. In order to avoid coordinate singularities at the north and south poles in the meridional direction, Legendre-Gauss-Radau points are chosen at the elements involving the two poles. Fourier polynomials are applied in the zonal direction for its periodicity, with only one element. Then, the partial differential equations are solved on the longitude-latitude meshes without coordinate transformation between spherical and Cartesian coordinates. For verification of the proposed method, a few Poisson equations and advective flows are tested. Firstly, the method is found to be valid for test cases with smooth solution. The results of the Poisson equations demonstrate that the present method exhibits high accuracy and exponential convergence. High- precision solutions are also obtained with near negligible numerical diffusion during the time evolution for advective flow with smooth shape. Secondly, the results of advective flow with non-smooth shape and deformational flow are also shown to be reasonable and effective. As a result, the present method is proved to be capable of solving flow through different types of elements, and thereby a desirable method with reliability and high accuracy for solving partial differential equations over a sphere.

RETROSPECTIVE TIME INTEGRAL SCHEME AND ITS APPLICATIONS TO THE ADVECTION EQUATION

To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme.

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.

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.

Solution of the Time Dependent Schrodinger Equation and the Advection Equation via Quantum Walk with Variable Parameters

We propose a solution method of Time Dependent Schr?dinger Equation (TDSE) and the advection equation by quantum walk/quantum cellular automaton with spatially or temporally variable parameters. Using numerical method, we establish the quantitative relation between the quantum walk with the space dependent parameters and the “Time Dependent Schr?dinger Equation with a space dependent imaginary diffusion coefficient” or “the advection equation with space dependent velocity fields”. Using the 4-point-averaging manipulation in the solution of advection equation by quantum walk, we find that only one component can be extracted out of two components of left-moving and right-moving solutions. In general it is not so easy to solve an advection equation without numerical diffusion, but this method provides perfectly diffusion free solution by virtue of its unitarity. Moreover our findings provide a clue to find more general space dependent formalisms such as solution method of TDSE with space dependent resolution by quantum walk.

A Comparative Study of Two Spatial Discretization Schemes for Advection Equation

In this paper, we describe a comparison of two spatial discretization schemes for the advection equation, namely the first finite difference method and the method of lines. The stability of the methods has been studied by Von Neumann method and with the matrix analysis. The methods are applied to a number of test problems to compare the accuracy and computational efficiency. We show that both discretization techniques approximate correctly solution of advection equation and compare their accuracy and performance.

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 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.

CONVERGENCE OF AN IMMERSED INTERFACE UPWIND SCHEME FOR LINEAR ADVECTION EQUATIONS WITH PIECEWISE CONSTANT COEFFICIENTS I:L^1-ERROR ESTIMATES

We study the L^l-error estimates for the upwind scheme to the linear advection equations with a piecewise constant coefficients modeling linear waves crossing interfaces. Here the interface condition is immersed into the upwind scheme. We prove that, for initial data with a bounded variation, the numerical solution of the immersed interface upwind scheme converges in L^l-norm to the differential equation with the corresponding interface condition. We derive the one-halfth order L^l-error bounds with explicit coefficients following a technique used in [25]. We also use some inequalities on binomial coefficients proved in a consecutive paper [32].

CONVERGENCE OF AN IMMERSED INTERFACE UPWIND SCHEME FOR LINEAR ADVECTION EQUATIONS WITH PIECEWISE CONSTANT COEFFICIENTS Ⅱ: SOME RELATED BINOMIAL COEFFICIENT INEQUALITIES

In this paper we give proof of three binomial coefficient inequalities. These inequalities are key ingredients in [Wen and Jin, J. Comput. Math. 26, （2008）, 1-22] to establish the L^1-error estimates for the upwind difference scheme to the linear advection equations with a piecewise constant wave speed and a general interface condition, which were further used to establish the L^1-error estimates for a Hamiltonian-preserving scheme developed in [Jin and Wen, Commun. Math. Sci. 3, （2005）, 285-315] to the Liouville equation with piecewise constant potentials [Wen and Jin, SIAM J. Numer. Anal. 46, （2008）, 2688-2714].

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.

Comparison of the performance of traditional advection-dispersion equation and mobile-immobile model for simulating solute transport in heterogeneous soils

The traditional advection-dispersion equation(ADE) and the mobile-immobile model(MIM) are widely used to describe solute transport in heterogeneous porous media. However, the fitness of the two models is casedependent. In this paper, the transport of conservative,adsorbing and degradable solutes through a 1 m heterogeneous soil column under steady flow condition was simulated by ADE and MIM, and sensitivity analysis was conducted. Results show that MIM tends to prolong the breakthrough process and decrease peak concentration for all three solutes, and tailing and skewness are more pronounced with increasing dispersivity. Breakthrough curves of the adsorbing solute simulated by MIM are less sensitive to the retardation factor compared with the results simulated by ADE. The breakthrough curves of degradable solute obtained by MIM and ADE nearly overlap with a high degradation rate coefficient, indicating that MIM and ADE perform similarly for simulating degradable solute transport when biochemical degradation prevails over the mass exchange between mobile and immobile zones. The results suggest that the physical significance of dispersivity should be carefully considered when MIM is applied to simulate the degradable solute transport and/or ADE is applied to simulate the adsorbing solute transport in highly dispersive soils.

A unified theory for gas dynamics and aeroacoustics in viscous compressible flows.PartⅡ.Sources on solid boundary

This work attempts to extend the fundamental theory for classic gas dynamics to viscous compressible flow,of which aeroacoustics will naturally be a special branch.As a continuation of Part I.Unbounded fluid(Mao et al.,2022),this paper studies the source of longitudinal field at solid boundary,caused by the on-wall kinematic and viscous dynamic coupling of longitudinal and transverse processes.We find that at this situation the easiest choice for the two independent thermodynamic variables is the dimensionless pressure P and temperature T.The two-level structure of boundary dynamics of longitudinal field is obtained by applying the continuity equation and its normal derivative to the surface.We show that the boundary dilatation flux represents faithfully the boundary production of vortex sound and entropy sound,and the mutual generation mechanism of the longitudinal and transverse fields on the boundary does not occur symmetrically"at the samc level,but appears along a zigzag route.At the first level,it is the pressure gradient that generates vorticity unidirectionally;while at the second level,it is the vorticity that generates dilatation unidirectionally.

A unified theory for gas dynamics and aeroacoustics in viscous compressible flows.Part I.Unbounded fluid

This paper presents a deep reflection on the advective wave equations for velocity vector and dilatation discovered in the past decade.We show that these equations can form the theoretical basis of modern gas dynamics,because they dominate not only various complex viscous and heat-conducting gas flows but also their associated longitudinal waves,including aero-generated sound.Current aeroacoustics theory has been developing in a manner quite independently of gas dynamics;it is based on the advective wave equations for thermodynamic variables,say the exact Phillips equation of relative disturbance pressure as a representative one.However,these equations do not cover the fluid flow that generates and propagates sound waves.In using them,one has to assume simplified base-flow models,which we argue is the main theoretical obstacle to identifying sound source and achieving effective noise control.Instead,we show that the Phillips equation and alike is nothing but the first integral of the dilatation equation that also governs the longitudinal part of the flow field.Therefore,we conclude that modern aeroacoustics should merge back into the general unsteady gas dynamics as a special branch of it,with dilatation of multiple sources being a new additional and sharper sound variable.

Time-Harmonic Acoustic Scattering in a Complex Flow:A Full Coupling Between Acoustics and Hydrodynamics

For the numerical simulation of time harmonic acoustic scattering in a complex geometry,in presence of an arbitrary mean flow,the main difficulty is the coexistence and the coupling of two very different phenomena:acoustic propagation and convection of vortices.We consider a linearized formulation coupling an augmented Galbrun equation(for the perturbation of displacement)with a time harmonic convection equation(for the vortices).We first establish the well-posedness of this time harmonic convection equation in the appropriatemathematical framework.Then the complete problem,with Perfectly Matched Layers at the artificial boundaries,is proved to be coercive+compact,and a hybrid numerical method for the solution is proposed,coupling finite elements for the Galbrun equation and a Discontinuous Galerkin scheme for the convection equation.Finally a 2D numerical result shows the efficiency of the method.