Source author record

D. S. Balsara

D. S. Balsara appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.IM astro-ph.EP astro-ph.GA physics.comp-ph

Catalog footprint

What is connected

3works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

Technologies for supporting high-order geodesic mesh frameworks for computational astrophysics and space sciences

Many important problems in astrophysics, space physics, and geophysics involve flows of (possibly ionized) gases in the vicinity of a spherical object, such as a star or planet. The geometry of such a system naturally favors numerical schemes based on a spherical mesh. Despite its orthogonality property, the polar (latitude-longitude) mesh is ill suited for computation because of the singularity on the polar axis, leading to a highly non-uniform distribution of zone sizes. The consequences are (a) loss of accuracy due to large variations in zone aspect ratios, and (b) poor computational efficiency from a severe limitations on the time stepping. Geodesic meshes, based on a central projection using a Platonic solid as a template, solve the anisotropy problem, but increase the complexity of the resulting computer code. We describe a new finite volume implementation of Euler and MHD systems of equations on a triangular geodesic mesh (TGM) that is accurate up to fourth order in space and time and conserves the divergence of magnetic field to machine precision. The paper discusses in detail the generation of a TGM, the domain decomposition techniques, three-dimensional conservative reconstruction, and time stepping.

preprint2013arXiv

MHD modeling on geodesic grids

This report describes a new magnetohydrodynamic numerical model based on a hexagonal spherical geodesic grid. The model is designed to simulate astrophysical flows of partially ionized plasmas around a central compact object, such as a star or a planet with a magnetic field. The geodesic grid, produced by a recursive subdivision of a base platonic solid (an icosahedron), is free from control volume singularities inherent in spherical polar grids. Multiple populations of plasma and neutral particles, coupled via charge-exchange interactions, can be simulated simultaneously with this model. Our numerical scheme uses piecewise linear reconstruction on a surface of a sphere in a local two-dimensional "Cartesian" frame. The code employs HLL-type approximate Riemann solvers and includes facilities to control the divergence of magnetic field and maintain pressure positivity. Several test solutions are discussed, including a problem of an interaction between the solar wind and the local interstellar medium, and a simulation of Earth's magnetosphere.

preprint2009arXiv

Dust Settling in Magnetorotationally-Driven Turbulent Discs II: The Pervasiveness of the Streaming Instability and its Consequences

We present a series of simulations of turbulent stratified protostellar discs with the goal of characterizing the settling of dust throughout a minimum-mass solar nebula. We compare the evolution of both compact spherical grains, as well as highly fractal grains. Our simulations use a shearing-box formulation to study the evolution of dust grains locally within the disc, and collectively our simulations span the entire extent of a typical accretion disc. The dust is stirred by gas that undergoes MRI-driven turbulence. This establishes a steady state scale height for the dust that is different for dust of different sizes. This sedimentation of dust is an important first step in planet formation and we predict that ALMA should be able to observationally verify its existence. When significant sedimentation occurs, the dust will participate in a streaming instability that significantly enhances the dust density. We show that the streaming instability is pervasive in the outer disc. We characterize the scale heights of dust whose size ranges from a few microns to a few centimeters. We find that for spherical grains, a power-law relationship develops for the scale height with grain size, with a slope that is slightly steeper than -1/2. The sedimentation is strongest in the outer disc and increases for large grains. The results presented here show that direct measurements of grain settling can be made by ALMA and we present favorable conditions for observability. The streaming instability should also be directly observable and we provide conditions for directly observing it. We calculate collision rates and growth rates for the dust grains in our simulations of various sizes colliding with other grains, and find that these rates are significantly enhanced through the density enhancement arising from the streaming instability.