Related papers: Self-Consistent Solution of Cosmological Radiation…
We present a scalable approach to solve a class of elliptic partial differential equation (PDE)-constrained optimization problems with bound constraints. This approach utilizes a robust full-space interior-point (IP)-Gauss-Newton…
In toroidally confined plasmas, the Grad-Shafranov equation, in general a non-linear PDE, describes the hydromagnetic equilibrium of the system. This equation becomes linear when the kinetic pressure is proportional to the poloidal magnetic…
Self similarity allows for analytic or semi-analytic solutions to many hydrodynamics problems. Most of these solutions are one dimensional. Using linear perturbation theory, expanded around such a one-dimensional solution, we find…
Quantum integrable systems, such as the interacting Bose gas in one dimension and the XXZ quantum spin chain, have an extensive number of local conserved quantities that endow them with exotic thermalization and transport properties. We…
We study confined solutions of certain evolutionary partial differential equations (pde) in 1+1 space-time. The pde we study are Lie-Poisson Hamiltonian systems for quadratic Hamiltonians defined on the dual of the Lie algebra of vector…
This paper is devoted to the study of the large time behaviour of viscosity solutions of parabolic equations with Neumann boundary conditions. This work is the sequel of [13] in which a probabilistic method was developped to show that the…
Recent advances in high-resolution imaging techniques and particle-based simulation methods have enabled the precise microscopic characterization of collective dynamics in various biological and engineered active matter systems. In…
We describe a computational framework for simulating suspensions of rigid particles in Newtonian Stokes flow. One central building block is a collision-resolution algorithm that overcomes the numerical constraints arising from particle…
A new method of solution is proposed for solution of the wave equation in one space dimension with continuously-varying coefficients. By considering all paths along which information arrives at a given point, the solution is expressed as an…
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 develop a new numerical scheme for solving the radiative transfer equation in a spherically symmetric system. This scheme does not rely on any kind of diffusion approximation and it is accurate for optically thin, thick, and intermediate…
In this work we introduce a novel semi-implicit structure-preserving finite-volume/finite-difference scheme for the viscous and resistive equations of magnetohydrodynamics (MHD) based on an appropriate 3-split of the governing PDE system,…
A method is developed within an adaptive framework to solve quasilinear diffusion problems with internal and possibly boundary layers starting from a coarse mesh. The solution process is assumed to start on a mesh where the problem is badly…
The dynamics of self-gravitating fluid bodies is described by the Euler-Einstein system of partial differential equations. The break-down of well-posedness on the fluid-vacuum interface remains a challenging open problem, which is…
Exact self-consistent particle-like solutions with spherical and/or cylindrical symmetry to the equations governing the interacting system of scalar, electromagnetic and gravitational fields have been obtained. As a particular case it is…
Ion transport, often described by the Poisson--Nernst--Planck (PNP) equations, is ubiquitous in electrochemical devices and many biological processes of significance. In this work, we develop conservative, positivity-preserving, energy…
We describe an extension of the Enzo code to enable fully-coupled radiation hydrodynamical simulation of inhomogeneous reionization in large $\sim (100 Mpc)^3$ cosmological volumes with thousands to millions of point sources. We solve all…
The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely…
The density functional approach in the Kohn-Sham approximation is widely used to study properties of many-electron systems. Due to the nonlinearity of the Kohn-Sham equations, the general self-consistence searching method involves…
In this work, we investigate the existence and uniqueness properties of a composite structure (multilayered) fluid interaction PDE system which arises in multi-physics problems, and particularly in biofluidic applications related to the…