Related papers: A Godunov Method for Multidimensional Radiation Ma…
We present a discontinuous Galerkin-finite-difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics. The hybrid method is conceptually quite simple.…
Magnetohydrodynamics (MHD) describes the interaction between electrically conducting fluids and electromagnetic fields. We propose and analyze a symplectic, second-order algorithm for the evolutionary MHD system in Els\"asser variables. We…
This article concerns a scalar multidimensional conservation law where the flux is of Panov type and may contain spatial discontinuities. We define a notion of entropy solution and prove that entropy solutions are unique. We propose a…
We construct a new Godunov type relativistic hydrodynamics code in Milne coordinates, using a Riemann solver based on the two-shock approximation which is stable under the existence of large shock waves. We check the correctness of the…
We consider the phenomenon of magnetic reconnection, namely a change in the topology of magnetic lines, for sufficiently regular solutions of the three-dimensional periodic magnetohydrodynamic (MHD) equations. We provide examples where…
We present a virtual element method (VEM) for the numerical approximation of the electromagnetics subsystem of the resistive magnetohydrodynamics (MHD) model in two spatial dimensions. The major advantages of the virtual element method…
We present a new numerical scheme for solving the advection equation and its application to the Vlasov simulation. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent…
We describe a photon-conserving radiative transfer algorithm, using a spatially-adaptive ray tracing scheme, and its parallel implementation into the adaptive mesh refinement (AMR) cosmological hydrodynamics code, Enzo. By coupling the…
High-dimensional fractional reaction-diffusion equations have numerous applications in the fields of biology, chemistry, and physics, and exhibit a range of rich phenomena. While classical algorithms have an exponential complexity in the…
This paper studies the shallow Ritz method for solving the one-dimensional diffusion problem. It is shown that the shallow Ritz method improves the order of approximation dramatically for non-smooth problems. To realize this optimal or…
"Generalized Hydrodynamics" (GHD) stands for a model that describes one-dimensional \textit{integrable} systems in quantum physics, such as ultra-cold atoms or spin chains. Mathematically, GHD corresponds to nonlinear equations of kinetic…
A high-performance shooting algorithm is developed to compute time-periodic solutions of the free-surface Euler equations with spectral accuracy in double and quadruple precision. The method is used to study resonance and its effect on…
Collisionless plasmas, mostly present in astrophysical and space environments, often require a kinetic treatment as given by the Vlasov equation. Unfortunately, the six-dimensional Vlasov equation can only be solved on very small parts of…
Diffusive Radiation is a new type of radiation predicted to occur in randomly inhomogeneous media due to the multiple scattering of pseudophotons. This theoretical effect is now observed experimentally. The radiation is generated by the…
We develop a new, efficient, and accurate method to simulate frequency-domain borehole electromagnetic (EM) measurements acquired in the presence of three-dimensional (3D) variations of the anisotropic subsurface conductivity. The method is…
We present a robust and accurate numerical method for the anisotropic diffusion equation in curvilinear coordinates. This study extends the recent work [Muir et al., Computer Physics Communications, 2025] for solving the anisotropic…
In this paper we study the three-dimensional two-phase magnetohydrodynamic interface problem in a bounded domain. The two incompressible fluids are both Newtonian and the surface tension is considered. We shall use the Galerkin method to…
We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer…
Radiative transfer plays a major role in high-energy astrophysics. In multiple scenarios and in a broad range of energy scales, the coupling between matter and radiation is essential to understand the interplay between theory, observations…
We present an efficient and accurate grid method for solving the time-dependent Schr\"{o}dinger equation of atomic systems interacting with intense laser pulses. As usual, the angular part of the wave function is expanded in terms of…