Related papers: A simple cure for numerical shock instability in H…
The off-lattice Boltzmann (OLB) method consists of numerical schemes which are used to solve the discrete Boltzmann equation. Unlike the commonly used lattice Boltzmann method, the spatial and time steps are uncoupled in the OLB method. In…
In this paper, a class of high order numerical schemes is proposed for solving Hamilton-Jacobi (H-J) equations. This work is regarded as an extension of our previous work for nonlinear degenerate parabolic equations, see Christlieb et al.…
Accurately reconstructing and forecasting high-resolution (HR) states from computationally cheap low-resolution (LR) observations is central to estimation-and-control of spatio-temporal PDE systems. We develop a unified superresolution…
An unsolved issue in widely used methods such as Support Vector Data Description (SVDD) and Small Sphere and Large Margin SVM (SSLM) for anomaly detection is their nonconvexity, which hampers the analysis of optimal solutions in a manner…
I study a class of global, causal geodesically complete solutions to the spherically symmetric Einstein scalar field (SSESF) system . Extending results of Luk-Oh (Quantitative Decay Rates for Dispersive Solutions to the Einstein-Scalar…
High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its…
A high-frequency recovered fully discrete low-regularity integrator is constructed to approximate rough and possibly discontinuous solutions of the semilinear wave equation. The proposed method, with high-frequency recovery techniques, can…
We present a new numerical method of special relativistic resistive magnetohydrodynamics with scalar resistivity that can treat a range of phenomena, from nonrelativistic to relativistic (shock, contact discontinuity, and Alfv\'en wave).…
The most essential concept in concurrent multiscale methods involving atomistic-continuum coupling is how to define the relation between atomistic and continuum regions. A well-known coupling method that has been frequently employed in…
Compositional simulation is challenging, because of highly nonlinear couplings between multi-component flow in porous media with thermodynamic phase behavior. The coupled nonlinear system is commonly solved by the fully-implicit scheme.…
We study a deformation of the nonlinear Schr\"odinger equation recently derived in the context of deformation of hierarchies of integrable systems. This systematic method also led to known integrable equations such as the Camassa-Holm…
In this paper, we develop a new method for magnetohydrodynamics (MHD) using smoothed particle hydrodynamics (SPH). To describe MHD shocks accurately, the Godunov method is applied to SPH instead of artificial dissipation terms. In the…
This paper studies the nonlinear evolution of magnetic field turbulence in proximity of steady ideal MHD configurations characterized by a small electric current, a small plasma flow, and approximate flux surfaces, a physical setting that…
Molecular Dynamic (MD) approach is applied to study the converging cylindrical shock waves in a dense Lennard-Jones (LJ) fluid. MD method is based on tracking of the atom motions and hence it has an fundamental advantages over hydrodynamic…
We introduce new adaptive schemes for the one- and two-dimensional hyperbolic systems of conservation laws. Our schemes are based on an adaption strategy recently introduced in [{\sc S. Chu, A. Kurganov, and I. Menshov}, Appl. Numer. Math.,…
We implement advanced Riemann solvers HLLC and HLLD \cite{Mignone:2005ft,MUB:2009} together with an advanced constrained transport scheme \cite{Gardiner:2007nc} in a numerical-relativity neutrino-radiation magnetohydrodynamics code. We…
Statistical solutions of incompressible Euler describe turbulent dynamics as time-parameterized laws on $L^2$ whose multi-point correlations satisfy an infinite hierarchy of weak identities. Modern generative samplers for PDE forecasting…
A main disadvantage of many high-order methods for hyperbolic conservation laws lies in the famous Gibbs-Wilbraham phenomenon, once discontinuities appear in the solution. Due to the Gibbs-Wilbraham phenomenon, the numerical approximation…
Symplectic integrators offer vastly superior performance over traditional numerical techniques for conservative dynamical systems, but their application to \emph{dissipative} systems is inherently difficult due to dissipative systems' lack…
This work establishes nonlinear orbital asymptotic stability of scalar radiative shock profiles, namely, traveling wave solutions to the simplified model system of radiating gas \cite{Hm}, consisting of a scalar conservation law coupled…