相关论文: A physical-constraints-preserving genuinely multid…
In this paper, we introduce a hyperbolic model for entropy dissipative system of viscous conservation laws via a flux relaxation approach. We develop numerical schemes for the resulting hyperbolic relaxation system by employing the…
This is the first of two papers examining the critical collapse of spherically symmetric perfect fluids with the equation of state P = (Gamma -1)rho. Here we present the equations of motion and describe a computer code capable of simulating…
We consider solutions of two-dimensional $m \times m$ systems hyperbolic conservation laws that are constant in time and along rays starting at the origin. The solutions are assumed to be small $L^\infty$ perturbations of a constant state…
We present and discuss a novel approach to deal with conservation properties for the simulation of nonlinear complex porous media flows in the presence of: 1) multiscale heterogeneity structures appearing in the elliptic-pressure-velocity…
We report on the development of a test-particle based kinetic Monte Carlo code for large systems and its application to simulate matter in the continuum regime. Our code combines advantages of the Direct Simulation Monte Carlo and the…
A high-order numerical method is developed for solving the Cahn-Hilliard-Navier-Stokes equations with the Flory-Huggins potential. The scheme is based on the $Q_k$ finite element with mass lumping on rectangular grids, the second-order…
The hydrodynamic description of probabilistic ballistic annihilation, for which no conservation laws hold, is an intricate problem with hard sphere-like dynamics for which no exact solution exists. We consequently focus on simplified…
The complete integrability of a generalized Riemann type hydrodynamic system is studied by means of symplectic and differential-algebraic tools. A compatible pair of polynomial Poissonian structures, Lax type representation and related…
We present a finite volume method that is applicable to hyperbolic PDEs including spatially varying and semilinear nonconservative systems. The spatial discretization, like that of the well-known Clawpack software, is based on solving…
An efficient algorithm for solving Poisson's equation in two and three spatial dimensions is discussed. The algorithm, which is described in detail, is based on the integral form of Poisson's equation and utilizes spherical coordinates and…
Methods to solve the relativistic hydrodynamic equations are a key computational kernel in a large number of astrophysics simulations and are crucial to understanding the electromagnetic signals that originate from the merger of…
Vortices, turbulence, and unsteady non-laminar flows are likely both prominent and dynamically important features of astrophysical disks. Such strongly nonlinear phenomena are often difficult, however, to simulate accurately, and are…
This paper proposes and analyzes a robust and efficient second-order positivity-preserving constrained transport (PPCT) scheme for ideal magnetohydrodynamics (MHD) on non-staggered Cartesian meshes. The PPCT scheme ensures two critical…
This work provides a comparative assessment of several low-dissipation numerical schemes for hyperbolic conservation laws, highlighting their performance relative to the classical Harten-Lax-van Leer (HLL) schemes. The schemes under…
We present for astrophysical use a multi-dimensional numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference method on an Eulerian grid, called the Total Variation Diminishing…
Hamiltonian systems of ordinary and partial differential equations are fundamental mathematical models spanning virtually all physical scales. A critical property for the robustness and stability of computational methods in such systems is…
We consider the design of structure-preserving discretization methods for the solution of systems of boundary controlled Partial Differential Equations (PDEs) thanks to the port-Hamiltonian formalism. We first provide a novel general…
Central-upwind (CU) schemes are Riemann-problem-solver-free finite-volume methods widely applied to a variety of hyperbolic systems of PDEs. Exact solutions of these systems typically satisfy certain bounds, and it is highly desirable or…
We report on the development of a computational framework for the parallel, mesh-adaptive solution of systems of hyperbolic conservation laws like the time-dependent Euler equations in compressible gas dynamics or Magneto-Hydrodynamics…
We numerically solve fully (3+1)-dimensional relativistic hydrodynamical equation coupled with the baryon number conservation law without spatial symmetry. We discuss the effect of transverse expansion based on the deviation our numerical…