Related papers: A Fast Potential and Self-Gravity Solver for Non-A…
Multigrid solvers are among the most efficient methods for solving the Poisson equation, which is ubiquitous in computational physics. For example, in the context of incompressible flows, it is typically the costliest operation. The present…
The calculation of the underground density field from measured gravity data has been done by a variety of methods of varied types. The use of the vector gravity components is here addressed in order to develop one accurate gravity inversion…
Accretion disks around stars, or other central massive bodies, can support long-lived, slowly precessing $m=1$ disturbances in which the fluid motion is nearly Keplerian with non-zero eccentricity. We study such `slow modes' in disks that…
A Continuous Galerkin method-based approach is presented to compute the seismic normal modes of rotating planets. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core using a polynomial…
We present a Mathematica package for performing algebraic and numerical computations in cosmological models based on supersymmetric theories. The programs allow for (I) evaluation and study of the properties of a scalar potential in a large…
Aims. Linear magnetohydrostatic (MHS) models of solar magnetic fields balance plasma pressure gradients, gravity and Lorentz forces where the current density is composed of a linear force-free component and a cross-field component that…
The short--wave asymptotics (WKB) of spiral density waves in self-gravitating stellar discs is well suited for the study of the dynamics of tightly--wound wavepackets. But the textbook WKB theory is not well adapted to the study of the…
We give an exact solution of the quadratic gravity in D dimensions. The solution is a plane fronted wave metric with a cosmological constant. This metric solves not only the full quadratic gravity field equations but also the linearized…
We present Superbox, a particle-mesh code with high resolution sub-grids and an NGP (nearest grid point) force-calculation scheme based on the second derivatives of the potential. Superbox implements a fast low-storage FFT-algorithm, giving…
In this paper, we consider the three dimensional Vlasov equation with an inhomogeneous, varying direction, strong magnetic field. Whenever the magnetic field has constant intensity, the oscillations generated by the stiff term are periodic.…
As hydrostatic equilibrium of gaseous discs is partly governed by the gravity field, we have estimated the component caused by a vertically homogeneous disc, with a special attention for the outer regions where self-gravity classically…
Gravity is scale free. Thus gravity may form similar structures in self-gravitating systems on different scales. Indeed, observations of the interstellar medium, spiral disks and cosmic structures, reveal similar characteristics. The…
In \emph{Wang et al., A Shifted Laplace Rational Filter for Large-Scale Eigenvalue Problems}, the SLRF method was proposed to compute all eigenvalues of a symmetric definite generalized eigenvalue problem lying in an interval on the real…
For simulating incompressible flows by projection methods. it is generally accepted that the pressure-correction stage is the most time-consuming part of the flow solver. The objective of the present work is to develop a fast hybrid…
The simulation of fracture using continuum ductile damage models attains a pathological discretization dependence caused by strain localization, after loss of ellipticity of the problem, in regions whose size is connected to the spatial…
In the framework of metric $f(R)$ gravity, we find the dispersion relation for the propagation of tightly wound spiral density waves in the surface of rotating, self-gravitating disks. Also, new Toomre-like stability criteria for…
We describe a newly developed code for the extrapolation of nonlinear force-free coronal magnetic fields in spherical coordinates. The program uses measured vector magnetograms on the solar photosphere as input and solves the force-free…
This paper introduces a parallel directional fast multipole method (FMM) for solving N-body problems with highly oscillatory kernels, with a focus on the Helmholtz kernel in three dimensions. This class of oscillatory kernels requires a…
Self-gravitating discs are believed to play an important role in astrophysics in particular regarding the star and planet formation process. In this context, discs subject to an idealized cooling process, characterized by a cooling…
Massive circumstellar disks are prone to gravitational instabilities, which trigger the formation of spiral arms that can fragment into bound clumps under the right conditions. Two dimensional simulations of self-gravitating disks are…