Related papers: A Fast Numerical Scheme for Causal Relativistic Hy…
Brownian Dynamics simulations are an important tool for modeling the dynamics of soft matter. However, accurate and rapid computations of the hydrodynamic interactions between suspended, microscopic components in a soft material is a…
A recent paper by Lucas-Serrano et al. indicates that a high-resolution central (HRC) scheme is robust enough to yield accurate hydrodynamical simulations of special relativistic flows in the presence of ultrarelativistic speeds and strong…
We address a numerical instability that arises in the directionally split computation of hydrodynamic flows when shock fronts are parallel to a grid plane. Transverse oscillations in pressure, density and temperature are produced that are…
We present the publicly available moving-mesh hydrodynamics code Sprout. Sprout solves the equations of ideal hydrodynamics on an expanding Cartesian mesh. The expanding mesh can follow fluid outflows for several orders of magnitude with…
A numerical scheme of relaxation type is proposed to approximate hyperbolic conservation laws in canal networks. Physical conditions at the junction are given and a novel strategy based on [Briani, Natalini, Ribot, 2025] is introduced to…
We introduce a coupled Cahn-Hilliard Navier-Stokes model that governs the two-phase dynamics of a system that consists of a fluid and a solid phase and prove its thermodynamic consistency. Moreover, we present an associated fully-discrete…
Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and…
We present a numerical method to solve the equations of general relativistic hydrodynamics in a given external gravitational field. The method is based on a generalization of Roe's approximate Riemann solver for the non relativistic Euler…
We present a convergence analysis of an unconditionally energy-stable first-order semi-discrete numerical scheme designed for a hydrodynamic Q-tensor model, the so-called Beris-Edwards system, based on the Invariant Energy Quadratization…
An approach for constructing a low-dissipation numerical method is described. The method is based on a combination of the operator-splitting method, Godunov method, and piecewise-parabolic method on the local stencil. Numerical method was…
In this paper, we propose and analyze a positivity-preserving, energy stable numerical scheme for certain type reaction-diffusion systems involving the Law of Mass Action with the detailed balance condition. The numerical scheme is…
In this paper, we propose and analyze a second order accurate (in both time and space) numerical scheme for the Poisson-Nernst-Planck-Navier-Stokes system, which describes the ion electro-diffusion in fluids. In particular, the…
In this paper, we present a computationally efficient method for including fluid-solid interactions into direct numerical simulations of the Navier-Stokes equations. This method is found to be as powerful as our earlier formulation [J.…
Fully non-linear equations of motion for dissipative general relativistic multi-fluids can be obtained from an action principle involving the explicit use of lower dimensional matter spaces. More traditional strategies for incorporating…
This paper detailedly discusses the locally one-dimensional numerical methods for efficiently solving the three-dimensional fractional partial differential equations, including fractional advection diffusion equation and Riesz fractional…
Dissipative Particle Dynamics (DPD) is a popular simulation model for investigating hydrodynamic behavior of systems with non-negligible equilibrium thermal fluctuations. DPD employs soft core repulsive interactions between the system…
We present an alternative method to estimate the numerical viscosity in simulations of astrophysical objects, which is based in the damping of fluid oscillations. We apply the method to general relativistic hydrodynamic simulations using a…
Many problems in stellar astrophysics feature flows at low Mach numbers. Conventional compressible hydrodynamics schemes frequently used in the field have been developed for the transonic regime and exhibit excessive numerical dissipation…
A Riemann problem with prescribed initial conditions will produce one of three possible wave patterns corresponding to the propagation of the different discontinuities that will be produced once the system is allowed to relax. In general,…
In this paper, we develop the numerical inverse scattering transform (NIST) for solving the derivative nonlinear Schrodinger (DNLS) equation. The key technique involves formulating a Riemann-Hilbert problem (RHP) that is associated with the…