11 research outputs found
Solving Vlasov Equations Using NRxx Method
In this paper, we propose a moment method to numerically solve the Vlasov
equations using the framework of the NRxx method developed in [6, 8, 7] for the
Boltzmann equation. Due to the same convection term of the Boltzmann equation
and the Vlasov equation, it is very convenient to use the moment expansion in
the NRxx method to approximate the distribution function in the Vlasov
equations. The moment closure recently presented in [5] is applied to achieve
the globally hyperbolicity so that the local well-posedness of the moment
system is attained. This makes our simulations using high order moment
expansion accessible in the case of the distribution far away from the
equilibrium which appears very often in the solution of the Vlasov equations.
With the moment expansion of the distribution function, the acceleration in the
velocity space results in an ordinary differential system of the macroscopic
velocity, thus is easy to be handled. The numerical method we developed can
keep both the mass and the momentum conserved. We carry out the simulations of
both the Vlasov-Poisson equations and the Vlasov-Poisson-BGK equations to study
the linear Landau damping. The numerical convergence is exhibited in terms of
the moment number and the spatial grid size, respectively. The variation of
discretized energy as well as the dependence of the recurrence time on moment
order is investigated. The linear Landau damping is well captured for different
wave numbers and collision frequencies. We find that the Landau damping rate
linearly and monotonically converges in the spatial grid size. The results are
in perfect agreement with the theoretic data in the collisionless case
A Nonlinear Multigrid Steady-State Solver for Microflow
We develop a nonlinear multigrid method to solve the steady state of
microflow, which is modeled by the high order moment system derived recently
for the steady-state Boltzmann equation with ES-BGK collision term. The solver
adopts a symmetric Gauss-Seidel iterative scheme nested by a local Newton
iteration on grid cell level as its smoother. Numerical examples show that the
solver is insensitive to the parameters in the implementation thus is quite
robust. It is demonstrated that expected efficiency improvement is achieved by
the proposed method in comparison with the direct time-stepping scheme
Physics-based adaptivity of a spectral method for the Vlasov-Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space
We propose a spectral method for the 1D-1V Vlasov-Poisson system where the
discretization in velocity space is based on asymmetrically-weighted Hermite
functions, dynamically adapted via a scaling and shifting of the
velocity variable. Specifically, at each time instant an adaptivity criterion
selects new values of and based on the numerical solution of the
discrete Vlasov-Poisson system obtained at that time step. Once the new values
of the Hermite parameters and are fixed, the Hermite expansion is
updated and the discrete system is further evolved for the next time step. The
procedure is applied iteratively over the desired temporal interval. The key
aspects of the adaptive algorithm are: the map between approximation spaces
associated with different values of the Hermite parameters that preserves total
mass, momentum and energy; and the adaptivity criterion to update and
based on physics considerations relating the Hermite parameters to the
average velocity and temperature of each plasma species. For the discretization
of the spatial coordinate, we rely on Fourier functions and use the implicit
midpoint rule for time stepping. The resulting numerical method possesses
intrinsically the property of fluid-kinetic coupling, where the low-order terms
of the expansion are akin to the fluid moments of a macroscopic description of
the plasma, while kinetic physics is retained by adding more spectral terms.
Moreover, the scheme features conservation of total mass, momentum and energy
associated in the discrete, for periodic boundary conditions. A set of
numerical experiments confirms that the adaptive method outperforms the
non-adaptive one in terms of accuracy and stability of the numerical solution
Quantum Hydrodynamic Model by Moment Closure of Wigner Equation
In this paper, we derive the quantum hydrodynamics models based on the moment
closure of the Wigner equation. The moment expansion adopted is of the Grad
type firstly proposed in \cite{Grad}. The Grad's moment method was originally
developed for the Boltzmann equation. In \cite{Fan_new}, a regularization
method for the Grad's moment system of the Boltzmann equation was proposed to
achieve the globally hyperbolicity so that the local well-posedness of the
moment system is attained. With the moment expansion of the Wigner function,
the drift term in the Wigner equation has exactly the same moment
representation as in the Boltzmann equation, thus the regularization in
\cite{Fan_new} applies. The moment expansion of the nonlocal Wigner potential
term in the Wigner equation is turned to be a linear source term, which can
only induce very mild growth of the solution. As the result, the local
well-posedness of the regularized moment system for the Wigner equation remains
as for the Boltzmann equation
Recent Advances in Industrial and Applied Mathematics
This open access book contains review papers authored by thirteen plenary invited speakers to the 9th International Congress on Industrial and Applied Mathematics (Valencia, July 15-19, 2019). Written by top-level scientists recognized worldwide, the scientific contributions cover a wide range of cutting-edge topics of industrial and applied mathematics: mathematical modeling, industrial and environmental mathematics, mathematical biology and medicine, reduced-order modeling and cryptography. The book also includes an introductory chapter summarizing the main features of the congress. This is the first volume of a thematic series dedicated to research results presented at ICIAM 2019-Valencia Congress