Related papers: High-resolution simulations and visualization of p…
Radiative transfer is an important component of hydrodynamic simulations as it determines the thermal properties of a physical system. It is especially important in cases where heating and cooling regulate significant processes, such as in…
We present a method to simulate fluid flow on evolving surfaces, e.g., an oil film on a water surface. Given an animated surface (e.g., extracted from a particle-based fluid simulation) in three-dimensional space, we add a second simulation…
We investigate exact nonlinear waves on surfaces locally approximating the rotating sphere for two-dimensional inviscid incompressible flow. Our first system corresponds to a beta-plane approximation at the equator and the second to a gamma…
In the core accretion model, planetesimals grow by mutual collisions and engulfing millimeter-to-centimeter particles, i.e., pebbles. Pebble accretion can significantly increase the accretion efficiency and help explain the presence of…
We develop a block-structured solver for high-fidelity simulation of flows in complex geometries, based on overlapping (Chimera) meshes. The key components of the algorithm are a baseline dissipation-free central discretization and…
In this paper a new primal-dual mixed finite element method is introduced, aimed to model multiscale problems with several geometric subregions in the domain of interest. In each of these regions porous media fluid flow takes place, but…
While recent observational progress is converging on the detection of compact regions of thermal emission due to embedded protoplanets, further theoretical predictions are needed to understand the response of a protoplanetary disk to the…
In this paper we propose and analyze a preconditioner for a system arising from a finite element approximation of second order elliptic problems describing processes in highly het- erogeneous media. Our approach uses the technique of…
The dynamics of planetesimals plays an important role in planet formation, because their velocity distribution sets the growth rate to larger bodies. When planetesimals form in protoplanetary discs, their orbits are nearly circular and…
We explore the evolution of the eccentricity of an accretion disc perturbed by an embedded planet whose mass is sufficient to open a large gap in the disc. Various methods for representing the orbit-averaged motion of an eccentric disc are…
There has been disagreement currently about whether cooling in protoplanetary disks can be sufficiently fast to induce the formation of gas giant protoplanets via gravitational instabilities. Simulations by our own group and others indicate…
We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in…
A new discretization approach is presented for the simulation of flow in complex poro-fractured media described by means of the Discrete Fracture and Matrix Model. The method is based on the numerical optimization of a properly defined…
Meshless methods are commonly used to determine numerical solutions to partial differential equations (PDEs) for problems involving free surfaces and/or complex geometries, approximating spatial derivatives at collocation points via local…
The formation and early evolution of protoplanetary disks are governed by a wide variety of physical processes during a gravitational collapse. Observations have begun probing disks in their earliest stages, and have favored the…
A numerical scheme is developed for systems of conservation laws on manifolds which arise in high speed aerodynamics and magneto-aerodynamics. The systems are presented in an arbitrary coordinate system on the manifold and involve source…
Explicit simulations of fluid mixtures of highly size-dispersed particles are constrained by numerical challenges associated with identifying pair-interaction neighbors. Recent algorithmic developments have ameliorated these difficulties to…
We investigate outflows driven by a giant protoplanet using three-dimensional MHD nested grid simulations. We consider a local region around the protoplanet in the protoplanetary disk, and calculate three models: (a) unmagnetized disk…
Recent observations have detected excess H$\alpha$ emission from young stellar systems with an age of several Myr such as PDS 70. One-dimensional radiation-hydrodynamic models of shock-heated flows that we developed previously demonstrate…
Spatial distribution and growth of dust in a clumpy protoplanetary disk subject to vigorous gravitational instability and fragmentation is studied numerically with sub-au resolution using the FEOSAD code. Hydrodynamics equations describing…