BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Dinshaw S. Balsara

Dinshaw S. Balsara contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.comp-phastro-ph.HEastro-phastro-ph.GAastro-ph.IMmath.NAastro-ph.SRgr-qcNumerical Analysisphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - Emerging
16works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 published item(s)

preprint2016arXiv

A High-Order Relativistic Two-Fluid Electrodynamic Scheme with Consistent Reconstruction of Electromagnetic Fields and a Multidimensional Riemann Solver for Electromagnetism

In various astrophysics settings it is common to have a two-fluid relativistic plasma that interacts with the electromagnetic field. While it is common to ignore the displacement current in the ideal, classical magnetohydrodynamic limit, when the flows become relativistic this approximation is less than absolutely well-justified. In such a situation, it is more natural to consider a positively charged fluid made up of positrons or protons interacting with a negatively charged fluid made up of electrons. The two fluids interact collectively with the full set of Maxwell's equations. As a result, a solution strategy for that coupled system of equations is sought and found here. Our strategy extends to higher orders, providing increasing accuracy. Three important innovations are reported here. In our first innovation, the magnetic field within each zone is reconstructed in a divergence-free fashion while the electric field within each zone is reconstructed in a form that is consistent with Gauss' law. In our second innovation, a multidimensionally upwinded strategy is presented which ensures that the magnetic field can be updated via a discrete interpretation of Faraday's law and the electric field can be updated via a discrete interpretation of the generalized Ampere's law. Our third innovation consists of an efficient design for several popular IMEX schemes so that they provide strong coupling between the finite-volume-based fluid solver and the electromagnetic fields at high order. Several accuracy analyses are presented showing that our method meets its design accuracy in the MHD limit as well as in the limit of electromagnetic wave propagation. Several stringent test problems are also presented. We also present a relativistic version of the GEM problem, which shows that our algorithm can successfully adapt to challenging problems in high energy astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Subluminal Relativistic Magnetohydrodynamics Scheme with ADER-WENO Predictor and Multidimensional Riemann Solver-Based Corrector

The relativistic magnetohydrodynamics (RMHD) set of equations has recently seen increased use in astrophysical computations. Even so, RMHD codes remain fragile. The reconstruction can sometimes yield superluminal velocities in certain parts of the mesh. In this paper we present a reconstruction strategy that overcomes this problem by making a single conservative to primitive transformation per cell followed by higher order WENO reconstruction on a carefully chosen set of primitives that guarantee subluminal reconstruction of the flow variables. For temporal evolution via a predictor step we also present second, third and fourth order accurate ADER methods that keep the velocity subluminal during the predictor step. The RMHD system also requires the magnetic field to be evolved in a divergence-free fashion. In the treatment of classical numerical MHD the analogous issue has seen much recent progress with the advent of multidimensional Riemann solvers. By developing multidimensional Riemann solvers for RMHD, we show that similar advances extend to RMHD. As a result, the face-centered magnetic fields can be evolved much more accurately using the edge-centered electric fields in the corrector step. Those edge-centered electric fields come from a multidimensional Riemann solver for RMHD which we present in this paper. In this paper we also develop several new test problems for RMHD. We show that RMHD vortices can be designed that propagate on the computational mesh as self-preserving structures. These RMHD vortex test problems provide a means to do truly multidimensional accuracy testing for RMHD codes. Several other stringent test problems are presented. We show the importance of resolution in certain test problems. Our tests include a demonstration that RMHD vortices are stable when they interact with shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Linear Stability of Sheared and Magnetized Jets Without Current Sheets - Non-Relativistic Case

In a prior paper (Kim et al. 2015) we considered the linear stability of magnetized jets that carry no net electric current and do not have current sheets. In this paper, in addition to physically well-motivated magnetic field structures, we also include the effects of jet shear. The jets we study have finite thermal pressure in addition to having realistic magnetic field structures and velocity shear. We find that shear has a strongly stabilizing effect on various modes of jet instability. Increasing shear stabilizes the fundamental pinch modes at long wavelengths and short wavelengths. Increasing shear also stabilizes the first reflection pinch modes at short wavelengths. Increasing shear has only a very modest stabilizing effect on the fundamental kink modes at long wavelengths; however, increasing shear does have a strong stabilizing effect on the fundamental kink modes at short wavelengths. The first reflection kink modes are strongly stabilized by increasing shear at shorter wavelengths. Overall, we find that the combined effect of magnetic field and shear stabilizes jets more than shear alone. In addition to the results from a formal linear stability analysis, we present a novel way of visualizing and understanding jet stability. This gives us a deeper understanding of the enhanced stability of sheared, magnetized jets. We also emphasize the value of our numerical approach in understanding the linear stability of jets with realistic structure.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement

We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method. Due to the one-step nature of the underlying scheme, the resulting algorithm is particularly well suited for an AMR strategy on space-time adaptive meshes, i.e.with time-accurate local time stepping. The AMR property has been implemented 'cell-by-cell', with a standard tree-type algorithm, while the scheme has been parallelized via the Message Passing Interface (MPI) paradigm. The new scheme has been tested over a wide range of examples for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations of compressible gas dynamics and the equations of magnetohydrodynamics (MHD). High order in space and time have been confirmed via a numerical convergence study and a detailed analysis of the computational speed-up with respect to highly refined uniform meshes is also presented. We also show test problems where the presented high order AMR scheme behaves clearly better than traditional second order AMR methods. The proposed scheme that combines for the first time high order ADER methods with space--time adaptive grids in two and three space dimensions is likely to become a useful tool in several fields of computational physics, applied mathematics and mechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Linear Stability of Magnetized Jets Without Current Sheets: Non-Relativistic Case

In this paper we consider stability of magnetized jets that carry no net electric current and do not have current sheets (Gourgouliatos et al. 2012). The non-relativistic MHD equations are linearized around the background velocity and the magnetic field structure of the jet. The resulting linear equations are solved numerically inside the jet. We find that introduction of current-sheet-free magnetic field significantly improves jet stability relative to unmagnetized jets or magnetized jets with current sheets at their surface. This particularly applies to the fundamental pinch and kink modes - they become completely suppressed in a wide range of long wavelengths that are known to become most pernicious to jet stability when the evolution enters the non-linear regime. The reflection modes, both for the pinch and kink instability, also become progressively more stable with increased magnetization.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lagrangian ADER-WENO Finite Volume Schemes on Unstructured Triangular Meshes Based On Genuinely Multidimensional HLL Riemann Solvers

In this paper we use the genuinely multidimensional HLL Riemann solvers recently developed by Balsara et al. to construct a new class of computationally efficient high order Lagrangian ADER-WENO one-step ALE finite volume schemes on unstructured triangular meshes. A nonlinear WENO reconstruction operator allows the algorithm to achieve high order of accuracy in space, while high order of accuracy in time is obtained by the use of an ADER time-stepping technique based on a local space-time Galerkin predictor. The multidimensional HLL and HLLC Riemann solvers operate at each vertex of the grid, considering the entire Voronoi neighborhood of each node and allows for larger time steps than conventional one-dimensional Riemann solvers. The results produced by the multidimensional Riemann solver are then used twice in our one-step ALE algorithm: first, as a node solver that assigns a unique velocity vector to each vertex, in order to preserve the continuity of the computational mesh; second, as a building block for genuinely multidimensional numerical flux evaluation that allows the scheme to run with larger time steps compared to conventional finite volume schemes that use classical one-dimensional Riemann solvers in normal direction. A rezoning step may be necessary in order to overcome element overlapping or crossing-over. We apply the method presented in this article to two systems of hyperbolic conservation laws, namely the Euler equations of compressible gas dynamics and the equations of ideal classical magneto-hydrodynamics (MHD). Convergence studies up to fourth order of accuracy in space and time have been carried out. Several numerical test problems have been solved to validate the new approach.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Observational Diagnostics for Two-Fluid Turbulence in Molecular Clouds As Suggested by Simulations

We present high resolution simulations of two-fluid (ion-neutral) MHD turbulence with resolutions as large as 512^3. The simulations are supersonic and mildly sub-Alfvenic, in keeping with the conditions present in molecular clouds. Such turbulence is thought to influence star formation processes in molecular clouds because typical cores form on length scales that are comparable to the dissipation scales of this turbulence in the ions. The simulations are motivated by the fact that recent studies of isophotologue lines in molecular clouds have found significant differences in the linewidth-size relationship for neutral and ion species. The goals of this paper are to explain those observations using simulations and analytic theory, present a new set of density-based diagnostics by drawing on similar diagnostics that have been obtained by studying single-fluid turbulence, and show that our two-fluid simulations play a vital role in reconciling alternative models of star formation. The velocity-dependent diagnostics display a complementarity with the density-dependent diagnostics. We find that the linewidth-size relationships show a prominent difference between ions and neutrals when the line of sight is orthogonal to the mean field. This is because the MHD waves in the ions differ from hydrodynamic waves in the neutrals. We also find that the density probability distribution functions (PDFs) show prominent differences between the ions and neutrals when the line of sight is parallel to the mean field. This is because the velocity fluctuations in the ions tend to produce column density fluctuations along field lines. When the magnetic field makes an angle to the line of sight, both observable differences should be visible. These diagnostics should be easy for observers to test. This analysis assumes optically thin lines and a mean magnetic field that is uniform in direction within a cloud.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Two-dimensional HLLC Riemann Solver for Conservation Laws : Application to Euler and MHD Flows

In this paper we present a genuinely two-dimensional HLLC Riemann solver. On logically rectangular meshes, it accepts four input states that come together at an edge and outputs the multi-dimensionally upwinded fluxes in both directions. This work builds on, and improves, our prior work on two-dimensional HLL Riemann solvers. The HLL Riemann solver presented here achieves its stabilization by introducing a constant state in the region of strong interaction, where four one-dimensional Riemann problems interact vigorously with one another. A robust version of the HLL Riemann solver is presented here along with a strategy for introducing sub-structure in the strongly-interacting state. Introducing sub-structure turns the two-dimensional HLL Riemann solver into a two-dimensional HLLC Riemann solver. The sub-structure that we introduce represents a contact discontinuity which can be oriented in any direction relative to the mesh. The Riemann solver presented here is general and can work with any system of conservation laws. We also present a second order accurate Godunov scheme that works in three dimensions and is entirely based on the present multidimensional HLLC Riemann solver technology. The methods presented are cost-competitive with traditional higher order Godunov schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Direct Evidence for Two-Fluid Effects in Molecular Clouds

We present a combination of theoretical and simulation-based examinations of the role of two-fluid ambipolar drift on molecular line widths. The dissipation provided by ion-neutral interactions can produce a significant difference between the widths of neutral molecules and the widths of ionic species, comparable to the sound speed. We demonstrate that Alfven waves and certain families of magnetosonic waves become strongly damped on scales comparable to the ambipolar diffusion scale. Using the RIEMANN code, we simulate two-fluid turbulence with ionization fractions ranging from 10^{-2} to 10^{-6}. We show that the wave damping causes the power spectrum of the ion velocity to drop below that of the neutral velocity when measured on a relative basis. Following a set of motivational observations by Li & Houde (2008), we produce synthetic line width-size relations that shows a difference between the ion and neutral line widths, illustrating that two-fluid effects can have an observationally detectable role in modifying the MHD turbulence in the clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Efficient Implementation of ADER Schemes for Euler and Magnetohydrodynamical Flows on Structured Meshes -- Comparison with Runge-Kutta Methods

ADER (Arbitrary DERivative in space and time) methods for the time-evolution of hyperbolic conservation laws have recently generated a fair bit of interest. The ADER time update can be carried out in a single step, which is desirable in many applications. However, prior papers have focused on the theory while downplaying implementation details. The purpose of the present paper is to make ADER schemes accessible by providing two useful formulations of the method as well as their implementation details on three-dimensional structured meshes. We therefore provide a detailed formulation of ADER schemes for conservation laws with non-stiff source terms in nodal as well as modal space along with useful implementation-related detail. We also provide details for the efficient use of ADER schemes in obtaining the numerical flux for conservation laws as well as electric fields for divergence-free magnetohydrodynamics. An efficient WENO-based strategy for obtaining zone-averaged magnetic fields from face-centered magnetic fields in MHD is also presented. The schemes catalogued here have been implemented in the first author's RIEMANN code. The speed of ADER schemes is shown to be almost twice as fast as that of strong stability preserving Runge-Kutta time stepping schemes for all the orders of accuracy that we tested.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evaluating the Magnetorotational Instability's Dependence on Numerical Algorithms and Resolution

We have studied saturated, MRI-driven turbulence using three-dimensional, isothermal simulations with resolutions that extend from 64 to 192 zones in each direction. The simulations were performed with several higher order Godunov algorithms. A variety of reconstruction strategies as well as a variety of Riemann solvers are tried. We show that the details of the isothermal MRI-driven turbulence depend principally on the Riemann solver and secondarily on the reconstruction strategy. Furthermore, we find that the effective viscosity parameter parameter tends to show progressively smaller decrements with increasing resolution when the best reconstruction strategy (WENO) and the best Riemann solver (linearized)are used. We attribute this result to the more sophisticated dissipation mechanisms that are used in higher-order Godunov schemes. Spectral analysis and transfer functions have been used to quantify the dissipative processes in these higher-order Godunov schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Multidimensional HLLE Riemann solver; Application to Euler and Magnetohydrodynamic Flows

In this work we present a general strategy for constructing multidimensional Riemann solvers with a single intermediate state, with particular attention paid to detailing the two-dimensional Riemann solver. This is accomplished by introducing a constant resolved state between the states being considered, which introduces sufficient dissipation for systems of conservation laws. Closed form expressions for the resolved fluxes are also provided to facilitate numerical implementation. The Riemann solver is proved to be positively conservative for the density variable; the positivity of the pressure variable has been demonstrated for Euler flows when the divergence in the fluid velocities is suitably restricted so as to prevent the formation of cavitation in the flow. We also focus on the construction of multidimensionally upwinded electric fields for divergence-free magnetohydrodynamical flows. A robust and efficient second order accurate numerical scheme for two and three dimensional Euler and magnetohydrodynamic flows is presented. The scheme is built on the current multidimensional Riemann solver. The number of zones updated per second by this scheme on a modern processor is shown to be cost competitive with schemes that are based on a one-dimensional Riemann solver. However, the present scheme permits larger timesteps.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Divergence-Free Reconstruction of Magnetic Fields and WENO Schemes for Magnetohydrodynamics

Balsara (2001, J. Comput. Phys., 174, 614) showed the importance of divergence-free reconstruction in adaptive mesh refinement problems for magnetohydrodynamics (MHD) and the importance of the same for designing robust second order schemes for MHD was shown in Balsara (2004, ApJS, 151, 149). Second order accurate divergence-free schemes for MHD have shown themselves to be very useful in several areas of science and engineering. However, certain computational MHD problems would be much benefited if the schemes had third and higher orders of accuracy. In this paper we show that the reconstruction of divergence-free vector fields can be carried out with better than second order accuracy. As a result, we design divergence-free weighted essentially non-oscillatory (WENO) schemes for MHD that have order of accuracy better than second. A multistage Runge-Kutta time integration is used to ensure that the temporal accuracy matches the spatial accuracy. Accuracy analysis is carried out and it is shown that the schemes meet their design accuracy for smooth problems. Stringent tests are also presented showing that the schemes perform well on those tests.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Dust Settling in Magnetorotationally-Driven Turbulent Discs I: Numerical Methods and Evidence for a Vigorous Streaming Instability

(Abridged) In this paper we have used the RIEMANN code for computational astrophysics to study the interaction of a realistic distribution of dust grains with gas in a vertically stratified protostellar accretion disc. The disc was modeled to have the density and temperature of a minimum mass solar nebula, and was driven to a fully-developed turbulence via the magnetorotational instability (MRI). We find that the inclusion of standard dust to gas ratios does not have any significant effect on the MRI even when the dust sediments to the midplane of the accretion disc. The density distribution of the dust reaches a Gaussian profile, and the scale heights for the dust that we derive are shown to be proportional to the reciprocal of the square root of the dust radius. The largest dust shows a strong tendency to settle to the midplane of the accretion disc, and tends to organize itself into elongated clumps of high density. The dynamics of these clumps is shown to be consistent with a streaming instability. The streaming instability is seen to be very vigorous and persistent once it forms. Each stream of high density dust displays a reduced RMS velocity dispersion, and the densest clumpings of large dust are shown to form where the streams intersect. We have also shown that the mean free path and collision time for the dust that participates in the streaming instability is reduced by almost two orders of magnitude relative to the average mean free paths and collision times. We show that some of the large dust in our 10 au simulations should have a propensity for grain coalescence.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Simulating the boundary layer between a white dwarf and its accretion disc

We describe the results of numerical simulations of the dynamics of the boundary layer (BL) between the accretion disk and the surface of a non-magnetic white dwarf (WD) for different viscosities which corresponds to different stages for dwarf novae burst cycles. The simulations cover the inner part of the accretion disk, the BL, and the upper atmosphere of the star. The high viscosity case, which corresponds to a dwarf nova in outburst, shows a optically thick BL which after one Keplerian rotation period (t_K=19s) extends more than 30 degrees to either side of the disk plane. The BL is optically thick and thus occludes part of the star. The low viscosity case, which corresponds to a dwarf nova in quiescense, also shows a BL, but it is optically thin.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

The accretion and spreading of matter on white dwarfs

For a slowly rotating non-magnetized white dwarf the accretion disk extends all the way to the star. Here the matter impacts and spreads towards the poles as new matter continuously piles up behind it. We have solved the 3d compressible Navier-Stokes equations on an axisymmetric grid to determine the structure of this boundary layer for different viscosities corresponding to different accretion rates. The high viscosity cases show a spreading BL which sets off a gravity wave in the surface matter. The accretion flow moves supersonically over the cusp making it susceptible to the rapid development of gravity wave and/or Kelvin-Helmholtz instabilities. This BL is optically thick and extends more than 30 degrees to either side of the disk plane after 3/4 of a Keplerian rotation period (t=19s). The low viscosity cases also show a spreading BL, but here the accretion flow does not set off gravity waves and it is optically thin.

0 citations0 reviewsNo code linkedOpen access