Related papers: Lagrangian versus Eulerian Methods for Toroidally-…
In this paper, we construct a novel Eulerian-Lagrangian finite volume (ELFV) method for nonlinear scalar hyperbolic equations in one space dimension. It is well known that the exact solutions to such problems may contain shocks though the…
In this work, we present a novel family of high order accurate numerical schemes for the solution of hyperbolic partial differential equations (PDEs) which combines several geometrical and physical structure preserving properties. First, we…
We perform direct numerical simulations (DNS) of passive heavy inertial particles (dust) in homogeneous and isotropic two-dimensional turbulent flows (gas) for a range of Stokes number, ${\rm St} < 1$, using both Lagrangian and Eulerian…
Global numerical simulations of protoplanetary disk formation and evolution were conducted in thin-disk limit, where the model included magnetically layered disk structure, a self-consistent treatment for the infall from cloud core as well…
Turbulent accretion disks threaded by a large-scale vertical field near equipartition can drive tenuous and fast self-confined jets. Self-similar solutions of these jet-emitting disks (JEDs) have been known for a long time and provide the…
This paper presents a series of global three dimensional accretion disk simulations carried out in the cylindrical limit in which the vertical component of the gravitational field is neglected. The simulations use a cylindrical…
We present an adaptive multiresolution method for the numerical simulation of ideal magnetohydrodynamics in two space dimensions. The discretization uses a finite volume scheme based on a Cartesian mesh and an explicit compact Rung-Kutta…
This paper proposes a semi-implicit arbitrary Lagrangian-Eulerian (ALE) method for the solution of the unified Godunov-Peshkov-Romenski (GPR) model of continuum mechanics. To handle the curl free involutions arising in the solid limit of…
Turbulence in protoplanetary disks, when present, plays a critical role in transporting dust particles embedded in the gaseous disk component. When using a field description of dust dynamics, a diffusion approach is traditionally used to…
Recent work demonstrated that breaking the up-down symmetry of tokamak flux surfaces removes a constraint that limits intrinsic momentum transport, and hence toroidal rotation, to be small. We show, through MHD analysis, that ellipticity is…
In this work, the benefits of the phase fitting technique are embedded in high order discrete Lagrangian integrators. The proposed methodology creates integrators with zero phase lag in a test Lagrangian in a similar way used in phase…
We present results of axisymmetric magnetohydrodynamic (MHD) simulations investigating the launching of jets and outflows from a magnetically diffusive accretion disk. The time evolution of the disk structure is self-consistently taken into…
We investigate accretion onto a central star, with the size, rotation rate, and magnetic dipole of a young stellar object, to study the flow pattern (velocity and density) of the fluid within and outside of the disc. We perform resistive…
Recent spectropolarimetric observations suggest that young low-mass stars such as classical T Tauri stars (CTTSs) possess relatively strong (~kG) magnetic field. This supports a scenario in which the final accretion onto the stellar surface…
The error of smoothed particle hydrodynamics (SPH) using kernel for particle-based approximation mainly comes from smoothing and integration errors. The choice of kernels has a significant impact on the numerical accuracy, stability and…
Using 2D MHD numerical simulations performed with two different finite difference Eulerian codes, we analyze the effect that a toroidal magnetic field has on low mass planet migration in nonturbulent protoplanetary discs. The presence of…
In this paper we present a conservative cell-centered Lagrangian finite volume scheme for the solution of the hyper-elasticity equations on unstructured multidimensional grids. The starting point of the new method is the Eucclhyd scheme,…
The Accretion-Ejection Instability (AEI) described by Tagger & Pellat (1999) is explored numerically using a global 2d model of the inner region of a magnetised accretion disk. The disk is initially currentless but threaded by a vertical…
Recently, Lobb and Nijhoff initiated the study of variational (Lagrangian) structure of discrete integrable systems from the perspective of multi-dimensional consistency. In the present work, we follow this line of research and develop a…
It has recently been established that the evolution of protoplanetary disks is primarily driven by magnetized disk winds, requiring large-scale magnetic flux threading the disks. The size of such disks is expected to shrink in time, as…