Related papers: Self-Consistent Solution of Cosmological Radiation…
In this paper, we develop a computational multiscale to solve the parabolic wave approximation with heterogeneous and variable media. Parabolic wave approximation is a technique to approximate the full wave equation. One benefit of the…
The axisymmetric form of the hydrodynamic equations within the smoothed particle hydrodynamics (SPH) formalism is presented and checked using idealized scenarios taken from astrophysics (free fall collapse, implosion and further pulsation…
We introduce a framework for solving a class of parabolic partial differential equations on triangle mesh surfaces, including the Hamilton-Jacobi equation and the Fokker-Planck equation. PDE in this class often have nonlinear or stiff terms…
We present a 3D finite element solver for the nonlinear Poisson-Nernst-Planck (PNP) equations for electrodiffusion, coupled to the Stokes system of fluid dynamics. The model serves as a building block for the simulation of macromolecule…
This paper is devoted to studying the global and finite convergence of the semi-smooth Newton method for solving a piecewise linear system that arises in cone-constrained quadratic programming problems and absolute value equations. We first…
The smoothed particle hydrodynamics (SPH) method is a useful numerical tool for the study of a variety of astrophysical and planetlogical problems. However, it turned out that the standard SPH algorithm has problems in dealing with…
Perfectly matched layers are a very efficient and accurate way to absorb waves in media. We present a stable convolutional unsplit perfectly matched formulation designed for the linearized stratified Euler equations. However, the technique…
We extend the results of the FBSDE theory in order to construct a probabilistic representation of a viscosity solution to the Cauchy problem for a system of quasilinear parabolic equations. We derive a BSDE associated with a class of…
Our study is dedicated to the probabilistic representation and numerical approximation of solutions to coupled systems of variational inequalities. The dynamics of each component of the solution is driven by a different linear parabolic…
In this work numerical methods for solving Einstein's equations are developed and applied to the study of inhomogeneous cosmological models. A two-dimensional computer code is described which implements two advanced numerical methods:…
A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is…
In this paper we introduce Enzo, a 3D MPI-parallel Eulerian block-structured adaptive mesh refinement cosmology code. Enzo is designed to simulate cosmological structure formation, but can also be used to simulate a wide range of…
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…
The gas of the interacted electrons is usually described within Kohn-Sham approximation by the set of Poisson and Schr\"{o}dinger equations with an effective potential for the single-particle wave functions. The solution of these equations…
In this paper a simple, robust, and general purpose approach to implement the Incompressible Smoothed Particle Hydrodynamics (ISPH) method is proposed. This approach is well suited for implementation on CPUs and GPUs. The method is…
Solving optimization problems with transient PDE-constraints is computationally costly due to the number of nonlinear iterations and the cost of solving large-scale KKT matrices. These matrices scale with the size of the spatial…
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…
During the past decades the study of strongly interacting fluids experienced a tremendous progress. In the relativistic heavy ion accelerators, specially the RHIC and LHC colliders, it became possible to study not only fluids made of…
We introduce a simple, rigorous, and unified framework for solving nonlinear partial differential equations (PDEs), and for solving inverse problems (IPs) involving the identification of parameters in PDEs, using the framework of Gaussian…
Shape Dynamics is a 3D conformally invariant theory of gravity which possesses a large set of solutions in common with General Relativity. When looked closely, these solutions are found to behave in surprising ways, so in order to probe the…