Related papers: Toward an efficient second-order method for comput…
A collection of algorithms is described for numerically computing with smooth functions defined on the unit sphere. Functions are approximated to essentially machine precision by using a structure-preserving iterative variant of Gaussian…
This paper constructs a theoretical framework for calculating the distribution of masses for secondary bodies forming via gravitational instability in the outer regions of circumstellar disks. We show that several alternate ways to specify…
In geophysical fluid dynamics, the screened Poisson equation appears in the shallow-water, quasi geostrophic equations. Recently, many attempts have been made to solve those equations on the sphere using different numerical methods. These…
We present the first application of the Poisson-Wiseman-Anderson method of matched expansions, to compute the self-force acting on a point particle moving in a curved spacetime. The method uses two expansions for the Green function, valid…
We extend the concept of galaxy environment from the local galaxy number density to the gravitational potential and its functions like the shear tensor. For this purpose we examine whether or not one can make an accurate estimation of the…
Thomson problem is a classical problem in physics to study how $n$ number of charged particles distribute themselves on the surface of a sphere of $k$ dimensions. When $k=2$, i.e. a 2-sphere (a circle), the particles appear at equally…
Recently it has been suggested that the fragmentation boundary in Smoothed Particle Hydrodynamic (SPH) and FARGO simulations of self-gravitating accretion discs with beta-cooling do not converge as resolution is increased. Furthermore, this…
We present a new method for generating the gravitational field of thin disks within the Weyl class of static and axially symmetric spacetimes. Such a gravitational field is described by two metric functions: one satisfies the Laplace…
We study the structure and dynamics of the gap created by a protoplanet in an accretion disc. The hydrodynamic equations for a flat, two-dimensional, non-selfgravitating protostellar accretion disc with an embedded, Jupiter sized…
Using a version of the ZEUS code, we carry out two-dimensional simulations of self-gravitating shearing sheets, with application to QSO accretion disks at a few thousand Schwarzschild radii, corresponding to a few hundredths of a parsec for…
An unsolved issue in the standard core accretion model for gaseous planet formation is how kilometre-sized planetesimals form from, initially, micron-sized dust grains. Solid growth beyond metre sizes can be difficult both because the…
We report on the results of the first 3D SPH simulation of massive, gravitationally unstable protoplanetary disks with radiative transfer. We adopt a flux-limited diffusion scheme justified by the high opacity of most of the disk. The…
General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their…
We discuss the ability of the smoothed particle hydrodynamics (SPH) method combined with a grid-based solver for the Poisson equation to model mass accretion onto protostars in gravitationally unstable protostellar discs. We scrutinize…
Self-gravitating accretion discs in a gravitoturbulent state, including radiation and gas pressures, are studied using a set of new analytical solutions. While the Toomre parameter of the disc remains close to its critical value for the…
The Gravitational Instability (GI) is a leading theory for explaining early planet formation in massive discs. In the early 2010s, 3D SPH simulations of GI failed to converge, initially attributed to resolution-dependent viscosity but later…
Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional Smoothed Particle Hydrodynamics (SPH) code by Meru and Bate, have been interpreted as implying that three-dimensional global discs fragment much…
In this paper, we extend our previous analysis (Lodato & Rice 2004) of the transport properties induced by gravitational instabilities in cooling, gaseous accretion discs to the case where the disc mass is comparable to the central object.…
A Green's function based solver for the modified Bessel equation has been developed with the primary motivation of solving the Poisson equation in cylindrical geometries. The method is implemented using a Discrete Hankel Transform and a…
The solution of the Poisson equation is a ubiquitous problem in computational astrophysics. Most notably, the treatment of self-gravitating flows involves the Poisson equation for the gravitational field. In hydrodynamics codes using…