Related papers: A general-purpose timestep criterion for simulatio…
Tests of general relativity (GR) with gravitational waves (GWs) introduce additional deviation parameters in the waveform model. The enlarged parameter space makes inference computationally costly, which has so far limited systematic,…
We present a novel multiscale numerical approach that combines parallel-in-time computation with hybrid domain adaptation for linear collisional kinetic equations in the diffusive regime. The method addresses the computational challenges of…
For the first time we introduce an error estimator for the numerical approximation of the equations describing the dynamics of sea ice. The idea of the estimator is to identify different error contributions coming from spatial and temporal…
The time step criterion plays a crucial role in direct N-body codes. If not chosen carefully, it will cause a secular drift in the energy error. Shared, adaptive time step criteria commonly adopt the minimum pairwise time step, which…
A longstanding problem in galactic simulations is to resolve the dynamical friction (DF) force acting on massive black hole particles when their masses are comparable to or less than the background simulation particles. Many sub-grid models…
Dynamic simulators are computational models governed by differential equations that evolve over time. They are essential for scientific and engineering applications but remain challenging to emulate because of the unpredictable behavior of…
The well-developed separate universe technique enables accurate calibration of the response of any observable to an isotropic long-wavelength density fluctuation. The large-scale environment also hosts tidal modes that perturb all…
A shift-symmetric Galileon model in presence of spacetime torsion has been constructed for the first time. This has been realized by localizing (or, gauging) the Galileon symmetry in flat spacetime in an appropriate manner. We have applied…
Context. High-resolution numerical methods have been developed for nonlinear, discontinuous problems as they appear in simulations of astrophysical objects. One of the strategies applied is the concept of artificial viscosity. Aims.…
The Galileon model is a tensor-scalar theory of gravity which explains the late acceleration of the Universe expansion with no instabilities and recovers General Relativity in the strong field limit. Most constraints obtained so far on…
We present a new three-parameter family of self-consistent equilibrium models for quasi-relaxed stellar systems that are subject to the combined action of external tides and rigid internal rotation. These models provide an idealised…
We introduce a highly-parallelizable architecture for estimating parameters of compact binary coalescence using gravitational-wave data and waveform models. Using a spherical harmonic mode decomposition, the waveform is expressed as a sum…
We derive equations of motion for bivariational wave functions with orthogonal adaptive basis sets and specialize the formalism to the coupled cluster ansatz. The equations are related to the biorthogonal case in a transparent way, and…
Advanced ground-based gravitational-wave detectors are capable of measuring tidal influences in binary neutron-star systems. In this work, we report on the statistical uncertainties in measuring tidal deformability with a full Bayesian…
Fluid dynamical simulations are often performed using cheap macroscopic models like the Euler equations. For rarefied gases under near-equilibrium conditions, however, macroscopic models are not sufficiently accurate and a simulation using…
We suggest a simple adaptive step-size procedure, which does not require any line-search, for a general class of nonlinear optimization methods and prove convergence of a general method under mild assumptions. In particular, the goal…
We present a methodology for simulating dilute suspensions of particles settling under gravity, with the main purpose of overcoming limitations of triply periodic configurations, mainly the strong vertical correlation that hinders the study…
This paper describes a flexible framework for generalized low-rank tensor estimation problems that includes many important instances arising from applications in computational imaging, genomics, and network analysis. The proposed estimator…
Gravitational wave measurements of binary neutron star coalescences offer information about the properties of the extreme matter that comprises the stars. Despite our expectation that all neutron stars in the Universe obey the same equation…
Dynamic tensor data are becoming prevalent in numerous applications. Existing tensor clustering methods either fail to account for the dynamic nature of the data, or are inapplicable to a general-order tensor. Also there is often a gap…