2,415 research outputs found
RAM: A Relativistic Adaptive Mesh Refinement Hydrodynamics Code
We have developed a new computer code, RAM, to solve the conservative
equations of special relativistic hydrodynamics (SRHD) using adaptive mesh
refinement (AMR) on parallel computers. We have implemented a
characteristic-wise, finite difference, weighted essentially non-oscillatory
(WENO) scheme using the full characteristic decomposition of the SRHD equations
to achieve fifth-order accuracy in space. For time integration we use the
method of lines with a third-order total variation diminishing (TVD)
Runge-Kutta scheme. We have also implemented fourth and fifth order Runge-Kutta
time integration schemes for comparison. The implementation of AMR and
parallelization is based on the FLASH code. RAM is modular and includes the
capability to easily swap hydrodynamics solvers, reconstruction methods and
physics modules. In addition to WENO we have implemented a finite volume module
with the piecewise parabolic method (PPM) for reconstruction and the modified
Marquina approximate Riemann solver to work with TVD Runge-Kutta time
integration. We examine the difficulty of accurately simulating shear flows in
numerical relativistic hydrodynamics codes. We show that under-resolved
simulations of simple test problems with transverse velocity components produce
incorrect results and demonstrate the ability of RAM to correctly solve these
problems. RAM has been tested in one, two and three dimensions and in
Cartesian, cylindrical and spherical coordinates. We have demonstrated
fifth-order accuracy for WENO in one and two dimensions and performed detailed
comparison with other schemes for which we show significantly lower convergence
rates. Extensive testing is presented demonstrating the ability of RAM to
address challenging open questions in relativistic astrophysics.Comment: ApJS in press, 21 pages including 18 figures (6 color figures
On Validating an Astrophysical Simulation Code
We present a case study of validating an astrophysical simulation code. Our
study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code
for studying the compressible, reactive flows found in many astrophysical
environments. We describe the astrophysics problems of interest and the
challenges associated with simulating these problems. We describe methodology
and discuss solutions to difficulties encountered in verification and
validation. We describe verification tests regularly administered to the code,
present the results of new verification tests, and outline a method for testing
general equations of state. We present the results of two validation tests in
which we compared simulations to experimental data. The first is of a
laser-driven shock propagating through a multi-layer target, a configuration
subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second
test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported
against the force of gravity by a light fluid. Our simulations of the
multi-layer target experiments showed good agreement with the experimental
results, but our simulations of the Rayleigh-Taylor instability did not agree
well with the experimental results. We discuss our findings and present results
of additional simulations undertaken to further investigate the Rayleigh-Taylor
instability.Comment: 76 pages, 26 figures (3 color), Accepted for publication in the ApJ
Simulating radiative shocks in nozzle shock tubes
We use the recently developed Center for Radiative Shock Hydrodynamics
(CRASH) code to numerically simulate laser-driven radiative shock experiments.
These shocks are launched by an ablated beryllium disk and are driven down
xenon-filled plastic tubes. The simulations are initialized by the
two-dimensional version of the Lagrangian Hyades code which is used to evaluate
the laser energy deposition during the first 1.1ns. The later times are
calculated with the CRASH code. This code solves for the multi-material
hydrodynamics with separate electron and ion temperatures on an Eulerian
block-adaptive-mesh and includes a multi-group flux-limited radiation diffusion
and electron thermal heat conduction. The goal of the present paper is to
demonstrate the capability to simulate radiative shocks of essentially
three-dimensional experimental configurations, such as circular and elliptical
nozzles. We show that the compound shock structure of the primary and wall
shock is captured and verify that the shock properties are consistent with
order-of-magnitude estimates. The produced synthetic radiographs can be used
for comparison with future nozzle experiments at high-energy-density laser
facilities.Comment: submitted to High Energy Density Physic
A study of planar Richtmyer-Meshkov instability in fluids with Mie-GrĂĽneisen equations of state
We present a numerical comparison study of planar Richtmyer-Meshkov instability with the intention of exposing the role of the equation of state. Results for Richtmyer-Meshkov instability in fluids with Mie-Grüneisen equations of state derived from a linear shock-particle speed Hugoniot relationship (Jeanloz, J. Geophys. Res. 94, 5873, 1989; McQueen et al., High Velocity Impact Phenomena (1970), pp. 294–417; Menikoff and Plohr, Rev. Mod. Phys. 61(1), 75 1989) are compared to those from perfect gases under nondimensionally matched initial conditions at room temperature and pressure. The study was performed using Caltech’s Adaptive Mesh Refinement, Object-oriented C++ (AMROC) (Deiterding, Adaptive Mesh Refinement: Theory and Applications (2005), Vol. 41, pp. 361–372; Deiterding, “Parallel adaptive simulation of multi-dimensional detonation structures,” Ph.D. thesis (Brandenburgische Technische Universität Cottbus, September 2003)) framework with a low-dissipation, hybrid, center-difference, limiter patch solver (Ward and Pullin, J. Comput. Phys. 229, 2999 (2010)). Results for single and triple mode planar Richtmyer-Meshkov instability when a reflected shock wave occurs are first examined for mid-ocean ridge basalt (MORB) and molybdenum modeled by Mie-Grüneisen equations of state. The single mode case is examined for incident shock Mach numbers of 1.5 and 2.5. The planar triple mode case is studied using a single incident Mach number of 2.5 with initial corrugation wavenumbers related by k_1 = k_2+k_3. Comparison is then drawn to Richtmyer-Meshkov instability in perfect gases with matched nondimensional pressure jump across the incident shock, post-shock Atwood ratio, post-shock amplitude-to-wavelength ratio, and time nondimensionalized by Richtmyer’s linear growth time constant prediction. Differences in start-up time and growth rate oscillations are observed across equations of state. Growth rate oscillation frequency is seen to correlate directly to the oscillation frequency for the transmitted and reflected shocks. For the single mode cases, further comparison is given for vorticity distribution and corrugation centerline shortly after shock interaction. Additionally, we examine single mode Richtmyer-Meshkov instability when a reflected expansion wave is present for incident Mach numbers of 1.5 and 2.5. Comparison to perfect gas solutions in such cases yields a higher degree of similarity in start-up time and growth rate oscillations. The formation of incipient weak waves in the heavy fluid driven by waves emanating from the perturbed transmitted shock is observed when an expansion wave is reflected
A model of Mira's cometary head/tail entering the Local Bubble
We model the cometary structure around Mira as the interaction of an AGB wind
from Mira A, and a streaming environment. Our simulations introduce the
following new element: we assume that after 200 kyr of evolution in a dense
environment Mira entered the Local Bubble (low density coronal gas). As Mira
enters the bubble, the head of the comet expands quite rapidly, while the tail
remains well collimated for a 100 kyr timescale. The result is a
broad-head/narrow-tail structure that resembles the observed morphology of
Mira's comet. The simulations were carried out with our new adaptive grid code
WALICXE, which is described in detail.Comment: 12 pages, 8 figures (4 in color). Accepted for publication in The
Astrophysical Journa
Computational Models of Material Interfaces for the Study of Extracorporeal Shock Wave Therapy
Extracorporeal Shock Wave Therapy (ESWT) is a noninvasive treatment for a
variety of musculoskeletal ailments. A shock wave is generated in water and
then focused using an acoustic lens or reflector so the energy of the wave is
concentrated in a small treatment region where mechanical stimulation enhances
healing. In this work we have computationally investigated shock wave
propagation in ESWT by solving a Lagrangian form of the isentropic Euler
equations in the fluid and linear elasticity in the bone using high-resolution
finite volume methods. We solve a full three-dimensional system of equations
and use adaptive mesh refinement to concentrate grid cells near the propagating
shock. We can model complex bone geometries, the reflection and mode conversion
at interfaces, and the the propagation of the resulting shear stresses
generated within the bone. We discuss the validity of our simplified model and
present results validating this approach
Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics
In this paper we present a full-fledged scheme for the second order accurate,
divergence-free evolution of vector fields on an adaptive mesh refinement (AMR)
hierarchy. We focus here on adaptive mesh MHD. The scheme is based on making a
significant advance in the divergence-free reconstruction of vector fields. In
that sense, it complements the earlier work of Balsara and Spicer (1999) where
we discussed the divergence-free time-update of vector fields which satisfy
Stoke's law type evolution equations. Our advance in divergence-free
reconstruction of vector fields is such that it reduces to the total variation
diminishing (TVD) property for one-dimensional evolution and yet goes beyond it
in multiple dimensions. Divergence-free restriction is also discussed. An
electric field correction strategy is presented for use on AMR meshes. The
electric field correction strategy helps preserve the divergence-free evolution
of the magnetic field even when the time steps are sub-cycled on refined
meshes. The above-mentioned innovations have been implemented in Balsara's
RIEMANN framework for parallel, self-adaptive computational astrophysics which
supports both non-relativistic and relativistic MHD. Several rigorous, three
dimensional AMR-MHD test problems with strong discontinuities have been run
with the RIEMANN framework showing that the strategy works very well.Comment: J.C.P., figures of reduced qualit
- …