Related papers: Implementation of Sink Particles in the Athena Cod…
We describe the first grid-based simulations of the polar alignment of a circumbinary disk. We simulate the evolution of an inclined disk around an eccentric binary using the grid-based code ATHENA++. The use of a grid-based numerical code…
We compare here several modern versions of SPH with a particular focus on the impact of gradient accuracy. We examine specifically an approximation to the "linearly exact" gradients (aLE) with standard SPH kernel gradients and with linearly…
We study critical phenomena in the gravitational collapse of a radiation fluid. We perform numerical simulations in both spherical symmetry and axisymmetry, and observe critical scaling in both supercritical evolutions, which lead to the…
We study the computation of equilibrium points of electrostatic potentials: locations in space where the electrostatic force arising from a collection of charged particles vanishes. This is a novel scenario of optimization in which…
An overset grid method was used to investigate the interaction between a particle-laden flow and a circular cylinder. The overset grid method was implemented in the Pencil Code , a high-order finite-difference code for compressible flow…
We present the exact solutions for the collapse of a spherically-symmetric, cold (i.e., pressureless) cloud under its own self-gravity, valid for arbitrary initial density profiles and not restricted to the realm of self-similarity. These…
Tectonic faults are commonly modelled as Volterra or Somigliana dislocations in an elastic medium. Various solution methods exist for this problem. However, the methods used in practice are often limiting, motivated by reasons of…
We study the consistency and convergence of smoothed particle hydrodynamics (SPH), as a function of the interpolation parameters, namely the number of particles $N$, the number of neighbors $n$, and the smoothing length $h$, using…
We present a new methodology for simulating self-gravitating general-relativistic fluids. In our approach the fluid is modelled by means of Lagrangian particles in the framework of a general-relativistic (GR) Smooth Particle Hydrodynamics…
We describe similarity solutions that characterize the collapse of collisional gas onto scalefree perturbations in an Einstein-de Sitter universe. We consider the effects of radiative cooling and derive self-similar solutions under the…
The spherical collapse scenario has great importance in cosmology since it captures several crucial aspects of structure formation. The presence of self-similar solutions in the Einstein-de Sitter (EdS) model greatly simplifies its…
Dispersion of low-density rigid particles with complex geometries is ubiquitous in both natural and industrial environments. We show that while explicit methods for coupling the incompressible Navier-Stokes equations and Newton's equations…
Acoustic levitation provides a unique method for manipulating small particles as it completely evades effects from gravity, container walls, or physical handling. These advantages make it a tantalizing platform for studying complex…
A quantum annealer exploits quantum effects to solve a particular type of optimization problem. The advantage of this specialized hardware is that it effectively considers all possible solutions in parallel, thereby potentially…
Hydrodynamic simulations have been used to study accretion disks consisting of counterrotating components with an intervening shear layer(s). Configurations of this type can arise from the accretion of newly supplied counterrotating matter…
Dust configurations are the simplest models for astrophysical objects. Here we examine the gravitational collapse of an infinite cylinder of dust and give an analytic interior solution. Surprisingly, starting with a cylindrically symmetric…
We perform a systematic study of the dynamics of dust particles in protoplanetary disks with embedded planets using global 2-D and 3-D inviscid hydrodynamic simulations. Lagrangian particles have been implemented into magnetohydrodynamic…
Our goal is to present an elementary approach to the analysis and programming of sparse grid finite element methods. This family of schemes can compute accurate solutions to partial differential equations, but using far fewer degrees of…
We study mass outflows driven from accretion discs by radiation pressure due to spectral lines. To investigate non-axisymmetric effects, we use the Athena++ code and develop a new module to account for radiation pressure driving. In 2D, our…
We propose a special-purpose class of compression algorithms for efficient compression of Prolog programs. It is a dictionary-based compression method, specially designed for the compression of Prolog code, and therefore we name it PCA…