Related papers: Higher-order initial conditions for mixed baryon-C…
We develop new perturbation techniques for conducting convergence analysis of various first-order algorithms for a class of nonsmooth optimization problems. We consider the iteration scheme of an algorithm to construct a perturbed…
This is a tutorial and survey paper on Karush-Kuhn-Tucker (KKT) conditions, first-order and second-order numerical optimization, and distributed optimization. After a brief review of history of optimization, we start with some preliminaries…
We investigate the Lagrangian perturbation theory of a homogeneous and isotropic universe in the non-relativistic limit, and derive the solutions up to the fourth order. These solutions are needed for example for the next-to-leading order…
In this paper, we develop a class of high-order conservative methods for simulating non-equilibrium radiation diffusion problems. Numerically, this system poses significant challenges due to strong nonlinearity within the stiff source terms…
Prior to recombination photons, electrons, and atomic nuclei rapidly scattered and behaved, almost, like a single tightly-coupled photon-baryon plasma. We investigate here the accuracy of the tight-coupling approximation commonly used to…
Precise initial conditions (ICs) are crucial for accurate computation in cosmological perturbation theory. We derive the consistent ICs for Horndeski theory in the Effective Field Theory (EFT) approach, assuming constant EFT functions at…
First-order variational equations are widely used in N-body simulations to study how nearby trajectories diverge from one another. These allow for efficient and reliable determinations of chaos indicators such as the Maximal Lyapunov…
The Lagrangian theory of structure formation in cosmological fluids, restricted to the matter model ``dust'', provides successful models of large-scale structure in the Universe in the laminar regime, i.e., where the fluid flow is…
We report simulations of the inspiral and merger of binary neutron stars performed with \texttt{WhiskyTHC}, the first of a new generation of numerical relativity codes employing higher than second-order methods for both the spacetime and…
Most numerical methods developed for solving nonlinear programming problems are designed to find points that satisfy certain optimality conditions. While the Karush-Kuhn-Tucker conditions are well-known, they become invalid when constraint…
We critically examine how the evolution of the matter density field in cosmological simulations is affected by details of setting up initial conditions. We show that it is non-trivial to realise an initial distribution of matter in…
We build a transient multidimensional multiphysical model based on continuum theories, involving the coupled mechanical, thermal and electrochemical phenomena occurring simultaneously in the discharge or charge of lithium-ion batteries. The…
Confronting measurements of the Lyman-$\alpha$ forest with cosmological hydrodynamical simulations has produced stringent constraints on models of particle dark matter and the thermal and ionization state of the intergalactic medium. We…
We review the field of collisionless numerical simulations for the large-scale structure of the Universe. We start by providing the main set of equations solved by these simulations and their connection with General Relativity. We then…
A fluid-structure interaction model in a port-Hamiltonian representation is derived for a classical guitar. We combine the laws of continuum mechanics for solids and fluids within a unified port-Hamiltonian (pH) modeling approach by…
We present novel coupling schemes for partitioned multi-physics simulation that combine four important aspects for strongly coupled problems: implicit coupling per time step, fast and robust acceleration of the corresponding iterative…
We present a new fast and efficient approach to model structure formation with Augmented Lagrangian Perturbation Theory (ALPT). Our method is based on splitting the displacement field into a long and a short-range component. The long-range…
We present and test a new numerical method to determine second-order Lagrangian displacement fields in the context of modified gravity (MG) theories. We start from the extension of Lagrangian Perturbation Theory to a class of MG models that…
We present a simple and intuitive approximation for solving perturbation theory (PT) of small cosmic fluctuations. We consider only the spherically symmetric or monopole contribution to the PT integrals, which yields the exact result for…
We present an improved prediction of the nonlinear perturbation theory (PT) via the Lagrangian picture, which was originally proposed by Matsubara (2008). Based on the relations between the power spectrum in standard PT and that in…