Crank-Nicolson Quasi-Compact Scheme for the Nonlinear Two-Sided Spatial Fractional Advection-Diffusion Equations

The higher-order numerical scheme of nonlinear advection-diffusion equations is studied in this article, where the space fractional derivatives are evaluated by using weighted and shifted Grünwald difference operators and combining the compact technique, in the time direction is discretized by the Crank-Nicolson method. Through the energy method, the stability and convergence of the numerical scheme in the sense of L2-norm are proved, and the convergence order is . Some examples are given to show that our numerical scheme is effective.

High Order IMEX Stochastic Galerkin Schemes for Linear Transport Equation with Random Inputs and Diffusive Scalings

In this paper,we consider the high order method for solving the linear transport equations under diffusive scaling and with random inputs.To tackle the randomness in the problem,the stochastic Galerkin method of the generalized polynomial chaos approach has been employed.Besides,the high order implicit-explicit scheme under the micro-macro decomposition framework and the discontinuous Galerkin method have been employed.We provide several numerical experiments to validate the accuracy and the stochastic asymptotic-preserving property.

Investigation of Acoustomagnetoelectric Effect in Bandgap Graphene by the Boltzmann Transport Equation

We study the acoustomagnetoelectric (AME) effect in two-dimensional graphene with an energy bandgap using the semiclassical Boltzmann transport equation within the hypersound regime, (where represents the acoustic wavenumber and is the mean free path of the electron). The Boltzmann transport equation and other relevant equations were solved analytically to obtain an expression for the AME current density, consisting of longitudinal and Hall components. Our numerical results indicate that both components of the AME current densities display oscillatory behaviour. Furthermore, geometric resonances and Weiss oscillations were each defined using the relationship between the current density and Surface Acoustic Wave (SAW) frequency and the inverse of the applied magnetic field, respectively. Our results show that the AME current density of bandgap graphene, which can be controlled to suit a particular electronic device application, is smaller than that of (gapless) graphene and is therefore, more suited for nanophotonic device applications.

Analysis of anomalous transport with temporal fractional transport equations in a bounded domain

Anomalous transport in magnetically confined plasmas is investigated using temporal fractional transport equations.The use of temporal fractional transport equations means that the order of the partial derivative with respect to time is a fraction. In this case, the Caputo fractional derivative relative to time is utilized, because it preserves the form of the initial conditions. A numerical calculation reveals that the fractional order of the temporal derivative α(α ∈(0, 1), sub-diffusive regime) controls the diffusion rate. The temporal fractional derivative is related to the fact that the evolution of a physical quantity is affected by its past history, depending on what are termed memory effects. The magnitude of α is a measure of such memory effects. When α decreases, so does the rate of particle diffusion due to memory effects. As a result,if a system initially has a density profile without a source, then the smaller the α is, the more slowly the density profile approaches zero. When a source is added, due to the balance of the diffusion and fueling processes, the system reaches a steady state and the density profile does not evolve. As α decreases, the time required for the system to reach a steady state increases. In magnetically confined plasmas, the temporal fractional transport model can be applied to off-axis heating processes. Moreover, it is found that the memory effects reduce the rate of energy conduction and hollow temperature profiles can be sustained for a longer time in sub-diffusion processes than in ordinary diffusion processes.

TRANSPORTATION COST-INFORMATION INEQUALITY FOR A STOCHASTIC HEAT EQUATION DRIVEN BY FRACTIONAL-COLORED NOISE

In this paper,we prove Talagrand’s T2 transportation cost-information inequality for the law of stochastic heat equation driven by Gaussian noise,which is fractional for a time variable with the Hurst index H∈(1/2,1),and is correlated for the spatial variable.The Girsanov theorem for fractional-colored Gaussian noise plays an important role in the proof.

Comparison between Non-Gaussian Puff Model and a Model Based on a Time-Dependent Solution of Advection-Diffusion Equation

A comparison between a non-Gaussian puff model and an advanced time-dependent model to simulate the pollutant dispersion in the Planetary Boundary Layer is presented. The puff model is based on a general technique for solving the K-equation, using the truncated Gram-Charlier expansion (type A) of the concentration field and finite set equations for the corresponding moments. The other model (named ADMM: Analytical Dispersion Multilayers Model) is an semi- analytical solution to the time-dependent two-dimensional advection-diffusion equation based on a discretization of the PBL in N sub-layers;in each sub-layers the advection-diffusion equation is solved by the Laplace transform technique, considering an average value for eddy diffusivity and the wind speed. A preliminary performance evaluation is shown in the case of continuous emission from an elevated source in a variable boundary layer. Both models were able to correctly reproduce the concentration field measured and so to be used as operative air pollution models.

Air Pollution Steady-State Advection-Diffusion Equation: The General Three-Dimensional Solution

Atmospheric air pollution turbulent fluxes can be assumed to be proportional to the mean concentration gradient. This assumption, along with the equation of continuity, leads to the advection-diffusion equation. Many models simulating air pollution dispersion are based upon the solution (numerical or analytical) of the advection-diffusion equation as- suming turbulence parameterization for realistic physical scenarios. We present the general steady three-dimensional solution of the advection-diffusion equation considering a vertically inhomogeneous atmospheric boundary layer for arbitrary vertical profiles of wind and eddy-diffusion coefficients. Numerical results and comparison with experimental data are shown.

Solution of Nonlinear Advection-Diffusion Equations via Linear Fractional Map Type Nonlinear QCA

Linear fractional map type (LFMT) nonlinear QCA (NLQCA), one of the simplest reversible NLQCA is studied analytically as well as numerically. Linear advection equation or Time Dependent Schrödinger Equation (TDSE) is obtained from the continuum limit of linear QCA. Similarly it is found that some nonlinear advection-diffusion equations including inviscid Burgers equation and porous-medium equation are obtained from LFMT NLQCA.

Two Modified QUICK Schemes for Advection-Diffusion Equation of Pollutants on Unstructured Grids

In this paper, two modified QUICK schemes, namely Q-QUICK and UQ-QUICK, for improving the preci-sion of convective flux approximation are verified in advection-diffusion equation of pollutants on unstruc-tured grids. The constructed auxiliary nodes for Q-QUICK/UQ-QUICK are composed of two neighboring nodes plus the next upwind node, the later node is generated from intersection of the line of current neighboring nodes and their corresponding interfaces. 2D unsteady advection-diffusion equation of pollut-ants is conducted for their verifications on unstructured grids. The numerical results show that Q-QUICK and UQ-QUICK have similar computational accuracy to the central difference scheme and similar numerical stability to upwind difference scheme after applying the deferred correction method. In addition, their corre-sponding CPU times are approximately equivalent to those of traditional difference schemes and their abili-ties for adapting high grid deformation are robust.

Discontinuous Legendre Wavelet Galerkin Method for One-Dimensional Advection-Diffusion Equation

This paper presents discontinuous Legendre wavelet Galerkin (DLWG) approach for solving one-dimensional advection-diffusion equation (ADE). Variational formulation of this type equation and corresponding numerical fluxes are devised by utilizing the advantages of both the Legendre wavelet bases and discontinuous Galerkin (DG) method. The distinctive features of the proposed method are its simple applicability for a variety of boundary conditions and able to effectively approximate the solution of PDEs with less storage space and execution. The results of a numerical experiment are provided to verify the efficiency of the designed new technique.

A Numerical Algorithm for the Caputo Tempered Fractional Advection-Diffusion Equation

By transforming the Caputo tempered fractional advection-diffusion equation into the Riemann–Liouville tempered fractional advection-diffusion equation,and then using the fractional-compact Grünwald–Letnikov tempered difference operator to approximate the Riemann–Liouville tempered fractional partial derivative,the fractional central difference operator to discritize the space Riesz fractional partial derivative,and the classical central difference formula to discretize the advection term,a numerical algorithm is constructed for solving the Caputo tempered fractional advection-diffusion equation.The stability and the convergence analysis of the numerical method are given.Numerical experiments show that the numerical method is effective.

A Third-Order Scheme for Numerical Fluxes to Guarantee Non-Negative Coefficients for Advection-Diffusion Equations

According to Godunov theorem for numerical calculations of advection equations, there exist no high-er-order schemes with constant positive difference coefficients in a family of polynomial schemes with an accuracy exceeding the first-order. In case of advection-diffusion equations, so far there have been not found stable schemes with positive difference coefficients in a family of numerical schemes exceeding the second-order accuracy. We propose a third-order computational scheme for numerical fluxes to guarantee the non-negative difference coefficients of resulting finite difference equations for advection-diffusion equations. The present scheme is optimized so as to minimize truncation errors for the numerical fluxes while fulfilling the positivity condition of the difference coefficients which are variable depending on the local Courant number and diffusion number. The feature of the present optimized scheme consists in keeping the third-order accuracy anywhere without any numerical flux limiter by using the same stencil number as convemtional third-order shemes such as KAWAMURA and UTOPIA schemes. We extend the present method into multi-dimensional equations. Numerical experiments for linear and nonlinear advection-diffusion equations were performed and the present scheme’s applicability to nonlinear Burger’s equation was confirmed.

Alternating Group Explicit Iterative Methods for One-Dimensional Advection-Diffusion Equation

The finite difference method such as alternating group iterative methods is useful in numerical method for evolutionary equations and this is the standard approach taken in this paper. Alternating group explicit (AGE) iterative methods for one-dimensional convection diffusion equations problems are given. The stability and convergence are analyzed by the linear method. Numerical results of the model problem are taken. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show that the behavior of the method with emphasis on treatment of boundary conditions is valuable.

SPIRAL SOLUTION TO THE TWO-DIMENSIONAL TRANSPORT EQUATIONS

The existence of spiral solution for the two-dimensional transport equations is considered in the present paper. Based on the notion of generalized solutions in the sense of Lebesgue-stieltjes integral, the global weak solution of transport equations which includes δ-shocks and vacuum is constructed for some special initial data.

WEAK SOLUTIONS OF MONGE-AMPRE TYPE EQUATIONS IN OPTIMAL TRANSPORTATION

This paper concerns the weak solutions of some Monge-Amp^re type equa- tions in the optimal transportation theory. The relationship between the Aleksandrov solutions and the viscosity solutions of the Monge-Ampere type equations is discussed. A uniform estimate for solution of the Dirichlet problem with homogeneous boundary value is obtained.

NUMERICAL SIMULATION OF SINGLE-GROUP, STEADY STATE AND ISOTROPIC NEUTRON TRANSPORT EQUATION IN DIFFUSIVE REGIMES

We present an algorithm for numerical solution of transport equation in diffusive regimes, in which the transport equation is nearly singular and its solution becomes a solution of a diffusion equation. This algorithm, which is based on the Least-squares FEM in combination with a scaling transformation, presents a good approximation of a diffusion operator in diffusive regimes and guarantees an accurate discrete solution. The numerical experiments in 2D and 3D case are given, and the numerical results show that this algorithm is correct and efficient.

Spherical harmonics method for neutron transport equation based on unstructured-meshes

Based on a new second-order neutron transport equation, self-adjoint angular flux (SAAF) equation, the spherical harmonics (PN) method for neutron transport equation on unstructured-meshes is derived. The spherical harmonics function is used to expand the angular flux. A set of differential equations about the spatial variable, which are coupled with each other, can be obtained. They are solved iteratively by using the finite element method on un- structured-meshes. A two-dimension transport calculation program is coded according to the model. The numerical results of some benchmark problems demonstrate that this method can give high precision results and avoid the ray effect very well.

A MODIFIED CRANK－NICOLSON SCHEME FOR THE INITIAL－BOUNDARY VALUE PROBLEM OF SUPERTHERMAL ELECTRON TRANSPORT EQUATION

The superthermal electron transport equation is a degenerate andnon-local evolutionary equation. In this paper, a modified Crank-Nicolsonscheme is constructed for the numerical solution. It is proved that the schemeis uniquely solvable and unconditionally

Nanofluids Transport Model Based on Fokker-Planck Equation and the Convection Heat Transfer Calculation

In current research about nanofluid convection heat transfer, random motion of nanoparticles in the liquid distribution problem mostly was not considered. In order to study on the distribution of nanoparticles in liquid, nanofluid transport model in pipe is established by using the continuity equation, momentum equation and Fokker-Planck equation. The velocity distribution and the nanoparticles distribution in liquid are obtained by numerical calculation, and the effect of particle size and particle volume fraction on convection heat transfer coefficient of nanofluids is analyzed. The result shows that in high volume fraction （ 0 _-- 0.8% ）, the velocity distribution of nanofluids characterizes as a ＂cork-shaped＂ structure, which is significantly different from viscous fluid with a parabolic distribution. The convection heat transfer coefficient increases while the particle size of nanoparticle in nanofluids decreases. And the convection heat transfer coefficient of nanofluids is in good agreement with the experimental result both in low （0 ~〈 0.1% ） and high （ q = 0.6% ） volume fractions. In presented model, Brown motion, the effect of interactions between nanoparticles and fluid coupling, is also considered, but any phenomenological parameter is not introduced. Nanoparticles in liquid transport distribution can be quantitatively calculated by this model.