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.

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.

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.

Asymptotic Behavior of Solutions for a Periodic Transport Equation with Delay

This paper analyzes a class of first order partial differential equations with delay (a model for the blood production system). The asymptotic behavior of solutions are studied.

NONLINEAR SEMIGROUP APPROACH TOTRANSPORT EQUATIONS WITH DELAYED NEUTRONS

This paper deal with a nonlinear transport equation with delayed neutron andgeneral boundary conditions. We establish, via the nonlinear semigroups approach, the exis-tence and uniqueness of the mild solution, weak solution, strong solution and local solutionon LP-spaces （1 ≤ p 〈 ＋∞）. Local and non local evolution problems are discussed.

ON THE TRANSPORT EQUATION WITH INTEGRAL BOUNDARY CONDITIONS AND ITS CONSERVAIVE LAW

The well-posedness for the time-dependent neutron transport equationwith integral boundary conditions is established in L^1 space. Some spectral propertiesof the transport operator are discussed, the dominant eigenvalue is proved existing,and furthermore, the conservative law of migrating particle numbers is established.

Global Lipschitz Stability in an Inverse Problem for the Non-stationary Transport Equation

An inverse problem of determining the source term of the non-stationary transport equation was considered.The extra lateral boundary condition was chosen as the observation data.Based on a Carleman estimate on the non-stationary transport equation,a global Lipschitz stability for the inverse problem was proved.

Analysis of High Concentration Spikes Appearing in Mass Plume in Nearly Homogeneous Turbulent Flow Based on the PDF Transport Equation

To evaluate the pollutant dispersion in background turbulent flows, most researches focus on statistical variation of concentrations or its fluctuations. However, those time-averaged quantities may be insufficient for risk assessment, because there emerge many high-intensity pollutant areas in the instantaneous concentration field. In this study, we tried to estimate the frequency of appearance of the high concentration areas in a turbulent flow based on the Probability Density Function (PDF) of concentration. The high concentration area was recognized by two conditions based on the concentration and the concentration gradient values. We considered that the estimation equation for the frequency of appearance of the recognized areas consisted of two terms based on each condition. In order to represent the two terms with physical quantities of velocity and concentration fields, simultaneous PIV (Particle Image Velocimetry) and PLIF (Planar Laser-Induced Fluorescence) measurement and PLIF time-serial measurement were performed in a quasi-homogeneous turbulent flow. According to the experimental results, one of the terms, related to the condition of the concentration, was found to be represented by the concentration PDF, while the other term, by the streamwise mean velocity and the integral length scale of the turbulent flow. Based on the results, we developed an estimation equation including the concentration PDF and the flow features of mean velocity and integral scale of turbulence. In the area where the concentration PDF was a Gaussian one, the difference between the frequencies of appearance estimated by the equation and calculated from the experimental data was within 25%, which showed good accuracy of our proposed estimation equation. Therefore, our proposed equation is feasible for estimating the frequency of appearance of high concentration areas in a limited area in turbulent mass diffusion.

Ultra-efficient and parameter-free computation of submicron thermal transport with phonon Boltzmann transport equation

Understanding thermal transport at the submicron scale is crucial for engineering applications,especially in the thermal management of electronics and tailoring the thermal conductivity of thermoelectric materials.At the submicron scale,the macroscopic heat diffusion equation is no longer valid and the phonon Boltzmann transport equation(BTE)becomes the governing equation for thermal transport.However,previous thermal simulations based on the phonon BTE have two main limitations:relying on empirical parameters and prohibitive computational costs.Therefore,the phonon BTE is commonly used for qualitatively studying the non-Fourier thermal transport phenomena of toy problems.In this work,we demonstrate an ultra-efficient and parameter-free computational method of the phonon BTE to achieve quantitatively accurate thermal simulation for realistic materials and devices.By properly integrating the phonon properties from first-principles calculations,our method does not rely on empirical material properties input.It can be generally applicable for different materials and the predicted results can match well with experimental results.Moreover,by developing a suitable ensemble of advanced numerical algorithms,our method exhibits superior numerical efficiency.The full-scale(from ballistic to diffusive)thermal simulation of a 3-dimensional fin field-effect transistor with 13 million degrees of freedom,which is prohibitive for existing phonon BTE solvers even on supercomputers,can now be completed within two hours on a single personal computer.Our method makes it possible to achieve the predictive design of realistic nanostructures for the desired thermal conductivity.It also enables accurately resolving the temperature profiles at the transistor level,which helps in better understanding the self-heating effect of electronics.

Convergence of the Modified Discrete Ordinates Method for the Anisotropic Scattering Transport Equation

Criticality problem of nuclear tractors generMly refers to an eigenvalue problem for the transport equations. In this paper, we deal with the eigenvalue of the anisotropic scattering transport equation in slab geometry. We propose a new discrete method which was called modified discrete ordinates method. It is constructed by redeveloping and improving discrete ordinates method in the space of L1（X）. Different from traditional methods, norm convergence of operator approximation is proved theoretically. Furthermore, convergence of eigenvalue approximation and the corresponding error estimation are obtained by analytical tools.

AN UNCONDITIONALLY STABLE HYBRID FE-FD SCHEME FOR SOLVING A 3-D HEAT TRANSPORT EQUATION IN A CYLINDRICAL THIN FILM WITH SUB-MICROSCALE THICKNESS

Heat transport at the microscale is of vital importance in microtechnology applications. The heat transport equation is different from the traditional heat transport equation since a second order derivative of temperature with respect to time and a third-order mixed derivative of temperature with respect to space and time are introduced. In this study, we develop a hybrid finite element-finite difference (FE-FD) scheme with two levels in time for the three dimensional heat transport equation in a cylindrical thin film with sub-microscale thickness. It is shown that the scheme is unconditionally stable. The scheme is then employed to obtain the temperature rise in a sub-microscale cylindrical gold film. The method can be applied to obtain the temperature rise in any thin films with sub-microscale thickness, where the geometry in the planar direction is arbitrary.

Transmission line model of carbon nanotubes:through the Boltzmann transport equation

A transmission line（TL）model of a carbon nanotube（CNT）is analyzed through the Boltzmann transport equation（BTE）.With the help of a numerical solution of the BTE,we study the kinetic inductance（L_K）, quantum capacitance（C_Q）and resistivity（R_S）of a CNT under a high frequency electric field.Values of L_K and C_Q obtained from BTE accord with the theoretical values,and the TL model is verified by transport theory for the first time.Moreover,our results show that the AC resistivity of CNTs deviates from DC,increasing along with shorter electric field wave length.This shows that changes in R_S in the high frequency condition must be considered in the TL model.

CONVERGENCE ANALYSIS FOR SPECTRAL APPROXIMATION TO A SCALAR TRANSPORT EQUATION WITH A RANDOM WAVE SPEED

This paper is concerned with the initial-boundary value problems of scalar transport equations with uncertain transport velocities. It was demonstrated in our earlier works that regularity of the exact solutions in the random spaces （or the parametric spaces） can be determined by the given data. In turn, these regularity results are crucial to convergence analysis for high order numerical methods. In this work, we will prove the spectral conver- gence of the stochastic Galerkin and collocation methods under some regularity results or assumptions. As our primary goal is to investigate the errors introduced by discretizations in the random space, the errors for solving the corresponding deterministic problems will be neglected.

A UNIFORM FIRST-ORDER METHOD FOR THE DISCRETE ORDINATE TRANSPORT EQUATION WITH INTERFACES IN X,Y-GEOMETRY

A uniformly first-order convergent numerical method for the discrete-ordinate transport equation in the rectangle geometry is proposed in this paper. Firstly we approximate the scattering coefficients and source terms by piecewise constants determined by their cell averages. Then for each cell, following the work of De Barros and Larsen [1, 19], the solution at the cell edge is approximated by its average along the edge. As a result, the solution of the system of equations for the cell edge averages in each cell can be obtained analytically. Finally, we piece together the numerical solution with the neighboring cells using the interface conditions. When there is no interface or boundary layer, this method is asymptotic-preserving, which implies that coarse meshes （meshes that do not resolve the mean free path） can be used to obtain good numerical approximations. Moreover, the uniform first-order convergence with respect to the mean free path is shown numerically and the rigorous proof is provided.

Self-similar Solutions for a Transport Equation with Non-local Flux

The authors construct self-similar solutions for an N-dimensional transport equation,where the velocity is given by the Riezs transform.These solutions imply nonuniqueness of weak solution.In addition,self-similar solution for a one-dimensional conservative equation involving the Hilbert transform is obtained.

A SPECTRAL STREAMLINE DIFFUSION FINITE ELEMENT COUPLING METHOD OF UNSTEADY TRANSPORT EQUATION IN THE FIELD OF NEUTRON LOGGING

In this paper, a new numerical method, the coupling method of spherical harmonic function spectral and streamline diffusion finite element for unsteady Boltzmann equation in the neutron logging field, is discussed. The convergence and error estimations of this scheme are proved. Its applications in the field of neutron logging show its effectiveness.

An Inverse Problem of Determining Coefficients in a One-Dimensional Radiative Transport Equation

We consider an inverse problem of determining unknown coefficients for a one-dimensional analogue of radiative transport equation.We show that some combination of the unknown coefficients can be uniquely determined by giving pulse-like inputs at the boundary and observing the corresponding outputs.Our result can be applied for determination of absorption and scattering properties of an optically turbid medium if the radiative transport equation is appropriate for describing the propagation of light in the medium.

Two Uniform Tailored Finite Point Schemes for the Two Dimensional Discrete Ordinates Transport Equations with Boundary and Interface Layers

This paper presents two uniformly convergent numerical schemes for the two dimensional steady state discrete ordinates transport equation in the diffusive regime,which is valid up to the boundary and interface layers.A five-point nodecentered and a four-point cell-centered tailored finite point schemes(TFPS)are introduced.The schemes first approximate the scattering coefficients and sources by piecewise constant functions and then use special solutions to the constant coefficient equation as local basis functions to formulate a discrete linear system.Numerically,both methods can not only capture the diffusion limit,but also exhibit uniform convergence in the diffusive regime,even with boundary layers.Numerical results show that the five-point scheme has first-order accuracy and the four-point scheme has second-order accuracy,uniformly with respect to the mean free path.Therefore a relatively coarse grid can be used to capture the two dimensional boundary and interface layers.