BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

J. B. Bell

J. B. Bell appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEmath.NAphysics.flu-dynastro-ph.IMastro-ph.SRcond-mat.mes-hallmath-phmath.MPastro-ph.COphysics.chem-phphysics.comp-ph

Catalog footprint

What is connected

22works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

22 published item(s)

preprint2019arXiv

Improved Coupling of Hydrodynamics and Nuclear Reactions via Spectral Deferred Corrections

Simulations in stellar astrophysics involve the coupling of hydrodynamics and nuclear reactions under a wide variety of conditions, from simmering convective flows to explosive nucleosynthesis. Numerical techniques such as operator splitting (most notably Strang splitting) are usually employed to couple the physical processes, but this can affect the accuracy of the simulation, particularly when the burning is vigorous. Furthermore, Strang splitting does not have a straightforward extension to higher-order integration in time. We present a new temporal integration strategy based on spectral deferred corrections and describe the second- and fourth-order implementations in the open-source, finite-volume, compressible hydrodynamics code Castro. One notable advantage to these schemes is that they combine standard low-order discretizations for individual physical processes in a way that achieves an arbitrarily high order of accuracy. We demonstrate the improved accuracy of the new methods on several test problems of increasing complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Comparisons of Two- and Three-Dimensional Convection in Type I X-ray Bursts

We perform the first detailed three-dimensional simulation of low Mach number convection preceding thermonuclear ignition in a mixed H/He X-ray burst. Our simulations include a moderate-sized, approximate network that captures hydrogen and helium burning up through rp-process breakout. We look in detail at the difference between two- and three-dimensional convective fields, including the details of the turbulent convection.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fluctuating hydrodynamics of multispecies reactive mixtures

We formulate and study computationally the fluctuating compressible Navier-Stokes equations for reactive multi-species fluid mixtures. We contrast two different expressions for the covariance of the stochastic chemical production rate in the Langevin formulation of stochastic chemistry, and compare both of them to predictions of the chemical Master Equation for homogeneous well-mixed systems close to and far from thermodynamic equilibrium. We develop a numerical scheme for inhomogeneous reactive flows, based on our previous methods for non-reactive mixtures [K. Balakrishnan, A. L. Garcia, A. Donev and J. B. Bell, Phys. Rev. E 89:013017, 2014]. We study the suppression of non-equilibrium long-ranged correlations of concentration fluctuations by chemical reactions, as well as the enhancement of pattern formation by spontaneous fluctuations. Good agreement with available theory demonstrates that the formulation is robust and a useful tool in the study of fluctuations in reactive multi-species fluids. At the same time, several problems with Langevin formulations of stochastic chemistry are identified, suggesting that future work should examine combining Langevin and Master Equation descriptions of hydrodynamic and chemical fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low Mach Number Fluctuating Hydrodynamics of Binary Liquid Mixtures

Continuing on our previous work [ArXiv:1212.2644], we develop semi-implicit numerical methods for solving low Mach number fluctuating hydrodynamic equations appropriate for modeling diffusive mixing in isothermal mixtures of fluids with different densities and transport coefficients. We treat viscous dissipation implicitly using a recently-developed variable-coefficient Stokes solver [ArXiv:1308.4605]. This allows us to increase the time step size significantly compared to the earlier explicit temporal integrator. For viscous-dominated flows, such as flows at small scales, we develop a scheme for integrating the overdamped limit of the low Mach equations, in which inertia vanishes and the fluid motion can be described by a steady Stokes equation. We also describe how to incorporate advanced higher-order Godunov advection schemes in the numerical method, allowing for the treatment of fluids with high Schmidt number including the vanishing mass diffusion coefficient limit. We incorporate thermal fluctuations in the description in both the inertial and overdamped regimes. We apply our algorithms to model the development of giant concentration fluctuations during the diffusive mixing of water and glycerol, and compare numerical results with experimental measurements. We find good agreement between the two, and observe propagative (non-diffusive) modes at small wavenumbers (large spatial scales), not reported in published experimental measurements of concentration fluctuations in fluid mixtures. Our work forms the foundation for developing low Mach number fluctuating hydrodynamics methods for miscible multi-species mixtures of chemically reacting fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low Mach Number Fluctuating Hydrodynamics of Multispecies Liquid Mixtures

We develop a low Mach number formulation of the hydrodynamic equations describing transport of mass and momentum in a multispecies mixture of incompressible miscible liquids at specified temperature and pressure that generalizes our prior work on ideal mixtures of ideal gases and binary liquid mixtures. In this formulation we combine and extend a number of existing descriptions of multispecies transport available in the literature. The formulation applies to non-ideal mixtures of arbitrary number of species, without the need to single out a 'solvent' species, and includes contributions to the diffusive mass flux due to gradients of composition, temperature and pressure. Momentum transport and advective mass transport are handled using a low Mach number approach that eliminates fast sound waves (pressure fluctuations) from the full compressible system of equations and leads to a quasi-incompressible formulation. Thermal fluctuations are included in our fluctuating hydrodynamics description following the principles of nonequilibrium thermodynamics. We extend the semi-implicit staggered-grid finite-volume numerical method developed in our prior work on binary liquid mixtures, and use it to study the development of giant nonequilibrium concentration fluctuations in a ternary mixture subjected to a steady concentration gradient. We also numerically study the development of diffusion-driven gravitational instabilities in a ternary mixture, and compare our numerical results to recent experimental measurements in a Hele-Shaw cell. We find that giant nonequilibrium fluctuations can trigger the instability but are eventually dominated by the deterministic growth of the unstable mode, in both quasi two-dimensional (Hele-Shaw), and fully three-dimensional geometries used in typical shadowgraph experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmological Fluid Mechanics with Adaptively Refined Large Eddy Simulations

We investigate turbulence generated by cosmological structure formation by means of large eddy simulations using adaptive mesh refinement. In contrast to the widely used implicit large eddy simulations, which resolve a limited range of length scales and treat the effect of turbulent velocity fluctuations below the grid scale solely by numerical dissipation, we apply a subgrid-scale model for the numerically unresolved fraction of the turbulence energy. For simulations with adaptive mesh refinement, we utilize a new methodology that allows us to adjust the scale-dependent energy variables in such a way that the sum of resolved and unresolved energies is globally conserved. We test our approach in simulations of randomly forced turbulence, a gravitationally bound cloud in a wind, and the Santa Barbara cluster. To treat inhomogeneous turbulence, we introduce an adaptive Kalman filtering technique that separates turbulent velocity fluctuations on resolved length scales from the non-turbulent bulk flow. From the magnitude of the fluctuating component and the subgrid-scale turbulence energy, a total turbulent velocity dispersion of several 100 km/s is obtained for the Santa Barbara cluster, while the low-density gas outside the accretion shocks is nearly devoid of turbulence. The energy flux through the turbulent cascade and the dissipation rate predicted by the subgrid-scale model correspond to dynamical time scales around 5 Gyr, independent of numerical resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficient Variable-Coefficient Finite-Volume Stokes Solvers

We investigate several robust preconditioners for solving the saddle-point linear systems that arise from spatial discretization of unsteady and steady variable-coefficient Stokes equations on a uniform staggered grid. Building on the success of using the classical projection method as a preconditioner for the coupled velocity-pressure system [B. E. Griffith, J. Comp. Phys., 228 (2009), pp. 75657595], as well as established techniques for steady and unsteady Stokes flow in the finite-element literature, we construct preconditioners that employ independent generalized Helmholtz and Poisson solvers for the velocity and pressure subproblems. We demonstrate that only a single cycle of a standard geometric multigrid algorithm serves as an effective inexact solver for each of these subproblems. Contrary to traditional wisdom, we find that the Stokes problem can be solved nearly as efficiently as the independent pressure and velocity subproblems, making the overall cost of solving the Stokes system comparable to the cost of classical projection or fractional step methods for incompressible flow, even for steady flow and in the presence of large density and viscosity contrasts. Two of the five preconditioners considered here are found to be robust to GMRES restarts and to increasing problem size, making them suitable for large-scale problems. Our work opens many possibilities for constructing novel unsplit temporal integrators for finite-volume spatial discretizations of the equations of low Mach and incompressible flow dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-Order Algorithms for Compressible Reacting Flow with Complex Chemistry

In this paper we describe a numerical algorithm for integrating the multicomponent, reacting, compressible Navier-Stokes equations, targeted for direct numerical simulation of combustion phenomena. The algorithm addresses two shortcomings of previous methods. First, it incorporates an eighth-order narrow stencil approximation of diffusive terms that reduces the communication compared to existing methods and removes the need to use a filtering algorithm to remove Nyquist frequency oscillations that are not damped with traditional approaches. The methodology also incorporates a multirate temporal integration strategy that provides an efficient mechanism for treating chemical mechanisms that are stiff relative to fluid dynamical time scales. The overall methodology is eighth order in space with option for fourth order to eighth order in time. The implementation uses a hybrid programming model designed for effective utilization of many-core architectures. We present numerical results demonstrating the convergence properties of the algorithm with realistic chemical kinetics and illustrating its performance characteristics. We also present a validation example showing that the algorithm matches detailed results obtained with an established low Mach number solver.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Low Mach Number Fluctuating Hydrodynamics of Diffusively Mixing Fluids

We formulate low Mach number fluctuating hydrodynamic equations appropriate for modeling diffusive mixing in isothermal mixtures of fluids with different density and transport coefficients. These equations eliminate the fluctuations in pressure associated with the propagation of sound waves by replacing the equation of state with a local thermodynamic constraint. We demonstrate that the low Mach number model preserves the spatio-temporal spectrum of the slower diffusive fluctuations. We develop a strictly conservative finite-volume spatial discretization of the low Mach number fluctuating equations in both two and three dimensions and construct several explicit Runge-Kutta temporal integrators that strictly maintain the equation of state constraint. The resulting spatio-temporal discretization is second-order accurate deterministically and maintains fluctuation-dissipation balance in the linearized stochastic equations. We apply our algorithms to model the development of giant concentration fluctuations in the presence of concentration gradients, and investigate the validity of common simplifications such as neglecting the spatial non-homogeneity of density and transport properties. We perform simulations of diffusive mixing of two fluids of different densities in two dimensions and compare the results of low Mach number continuum simulations to hard-disk molecular dynamics simulations. Excellent agreement is observed between the particle and continuum simulations of giant fluctuations during time-dependent diffusive mixing.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multidimensional Modeling of Type I X-ray Bursts. II. Two-Dimensional Convection in a Mixed H/He Accretor

Type I X-ray Bursts (XRBs) are thermonuclear explosions of accreted material on the surfaces of a neutron stars in low mass X-ray binaries. Prior to the ignition of a subsonic burning front, runaway burning at the base of the accreted layer drives convection that mixes fuel and heavy-element ashes. In this second paper in a series, we explore the behavior of this low Mach number convection in mixed hydrogen/helium layers on the surface of a neutron star using two-dimensional simulations with the Maestro code. Maestro takes advantage of the highly subsonic flow field by filtering dynamically unimportant sound waves while retaining local compressibility effects, such as those due to stratification and energy release from nuclear reactions. In these preliminary calculations, we find that the rp-process approximate network creates a convective region that is split into two layers. While this splitting appears artificial due to the approximations of the network regarding nuclear flow out of the breakout reaction 18Ne(a,p)21Na, these calculations hint at further simplifications and improvements of the burning treatment for use in subsequent calculations in three dimensions for a future paper.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Carbon Deflagration in Type Ia Supernova: I. Centrally Ignited Models

A leading model for Type Ia supernovae (SNe Ia) begins with a white dwarf near the Chandrasekhar mass that ignites a degenerate thermonuclear runaway close to its center and explodes. In a series of papers, we shall explore the consequences of ignition at several locations within such dwarfs. Here we assume central ignition, which has been explored before, however, the problem is worth revisiting, if only to validate those previous studies and to further elucidate the relevant physics for future work. A perturbed sphere of hot iron ash with a radius of ~100 km is initialized at the middle of the star. The subsequent explosion is followed in several simulations using a thickened flame model in which the flame speed is either fixed --- within the range expected from turbulent combustion --- or based on the local turbulent intensity. Global results, including the explosion energy and bulk nucleosynthesis (e.g. 56Ni of 0.48--0.56 $\Msun$) turn out to be insensitive to this speed. In all completed runs, the energy released by the nuclear burning is adequate to unbind the star, but not enough to give the energy and brightness of typical SNe Ia. As found previously, the chemical stratification observed in typical events is not reproduced. These models produce a large amount of unburned carbon and oxygen in central low velocity regions, which is inconsistent with spectroscopic observations, and the intermediate mass elements and iron group elements are strongly mixed during the explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Deflagration Stage of Chandrasekhar Mass Models For Type Ia Supernovae: I. Early Evolution

We present high-resolution, full-star simulations of the post-ignition phase of Type Ia supernovae using the compressible hydrodynamics code Castro. Initial conditions, including the turbulent velocity field and ignition site, are imported directly from a simulation of the last few hours of presupernova convection using a low Mach number code, Maestro. Adaptive mesh refinement allows the initial burning front to be modeled with an effective resolution of 36,864^3 zones (~136 m/zone). The initial rise and expansion of the deflagration front are tracked until burning reaches the star's edge and the role of the background turbulence on the flame is investigated. The effect of artificially moving the ignition location closer to the star's center is explored. The degree to which turbulence affects the burning front decreases with increasing ignition radius since the buoyancy force is stronger at larger radii. Even central ignition --- in the presence of a background convective flow field --- is rapidly carried off-center as the flame is carried by the flow field. We compare our results to analytic models for burning thermals, and find that they reproduce the general trends of the bubble's size and mass, but underpredict the amount of buoyant acceleration due to simplifying assumptions of the bubble's properties. Overall, we find that the amount of mass that burns prior to flame break out is small, consistent with a "gravitationally confined detonation" occurring at a later epoch, but additional burning will occur following breakout that may modify this conclusion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Low Mach Number Modeling of Convection in Helium Shells on Sub-Chandrasekhar White Dwarfs. I. Methodology

We assess the robustness of a low Mach number hydrodynamics algorithm for modeling helium shell convection on the surface of a white dwarf in the context of the sub-Chandrasekhar model for Type Ia supernovae. We use the low Mach number stellar hydrodynamics code, MAESTRO, to perform three-dimensional, spatially-adaptive simulations of convection leading up to the point of the ignition of a burning front. We show that the low Mach number hydrodynamics model provides a robust description of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Staggered Schemes for Fluctuating Hydrodynamics

We develop numerical schemes for solving the isothermal compressible and incompressible equations of fluctuating hydrodynamics on a grid with staggered momenta. We develop a second-order accurate spatial discretization of the diffusive, advective and stochastic fluxes that satisfies a discrete fluctuation-dissipation balance, and construct temporal discretizations that are at least second-order accurate in time deterministically and in a weak sense. Specifically, the methods reproduce the correct equilibrium covariances of the fluctuating fields to third (compressible) and second (incompressible) order in the time step, as we verify numerically. We apply our techniques to model recent experimental measurements of giant fluctuations in diffusively mixing fluids in a micro-gravity environment [A. Vailati et. al., Nature Communications 2:290, 2011]. Numerical results for the static spectrum of non-equilibrium concentration fluctuations are in excellent agreement between the compressible and incompressible simulations, and in good agreement with experimental results for all measured wavenumbers.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Diffusive Transport Enhanced by Thermal Velocity Fluctuations

We study the contribution of advection by thermal velocity fluctuations to the effective diffusion coefficient in a mixture of two indistinguishable fluids. The enhancement of the diffusive transport depends on the system size L and grows as \ln(L/L_0) in quasi two-dimensional systems, while in three dimensions it scales as L_0^{-1}-L^{-1}, where L_0 is a reference length. The predictions of a simple fluctuating hydrodynamics theory are compared to results from particle simulations and a finite-volume solver and excellent agreement is observed. Our results conclusively demonstrate that the nonlinear advective terms need to be retained in the equations of fluctuating hydrodynamics when modeling transport in small-scale finite systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Enhancement of Diffusive Transport by Nonequilibrium Thermal Fluctuations

We study the contribution of advection by thermal velocity fluctuations to the effective diffusion coefficient in a mixture of two identical fluids. The steady-state diffusive flux in a finite system subject to a concentration gradient is enhanced because of long-range correlations between concentration fluctuations and fluctuations of the velocity parallel to the concentration gradient. The enhancement of the diffusive transport depends on the system size L and grows as \ln(L/L_{0}) in quasi-two dimensional systems, while in three dimensions it grows as L_{0}^{-1}-L^{-1}, where L_{0} is a reference length. The predictions of a simple fluctuating hydrodynamics theory, closely related to second-order mode-mode coupling analysis, are compared to results from particle simulations and a finite-volume solver and excellent agreement is observed. We elucidate the direct connection to the long-time tail of the velocity autocorrelation function in finite systems, as well as finite-size corrections employed in molecular dynamics calculations. Our results conclusively demonstrate that the nonlinear advective terms need to be retained in the equations of fluctuating hydrodynamics when modeling transport in small-scale finite systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

From Convection to Explosion: End-to-End Simulation of Type Ia Supernovae

We present our end-to-end capability for computing the convective phase through the explosion phase of Type Ia supernovae. We compute the convective phase up to the time of ignition using our low Mach number code, MAESTRO, and the subsequent explosion phase using our compressible code, CASTRO. Both codes share the same BoxLib software framework and use finite-volume, block-structured adaptive mesh refinement (AMR) to enable high-resolution, three-dimensional full-star simulations that scale to 100,000+ cores. We present preliminary results from the first-ever simulations of convection preceding ignition using MAESTRO with AMR. We also demonstrate our ability to initialize a compressible simulation of the explosion phase in CASTRO using data obtained directly from MAESTRO just before ignition. Some care must be taken during this initialization procedure when interpreting the size and distribution of hot spots.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High-Resolution Simulations of Convection Preceding Ignition in Type Ia Supernovae Using Adaptive Mesh Refinement

We extend our previous three-dimensional, full-star simulations of the final hours of convection preceding ignition in Type Ia supernovae to higher resolution using the adaptive mesh refinement capability of our low Mach number code, MAESTRO. We report the statistics of the ignition of the first flame at an effective 4.34 km resolution, and general flow field properties at an effective 2.17 km resolution. We find that off-center ignition is likely, with radius of 50 km most favored and a likely range of 40 to 75 km. This is consistent with our previous coarser (8.68 km resolution) simulations, implying that we have achieved sufficient resolution in our determination of likely ignition radii. The dynamics of the last few hot spots preceding ignition suggest that a multiple ignition scenario is not likely. With improved resolution, we can more clearly see the general flow pattern in the convective region, characterized by a strong outward plume with a lower speed recirculation. We show that the convective core is turbulent with a Kolmogorov spectrum and has a lower turbulent intensity and larger integral length scale than previously thought (on the order of 16 km s$^{-1}$ and 200 km, respectively), and we discuss the potential consequences for the first flames.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

CASTRO: A New Compressible Astrophysical Solver. I. Hydrodynamics and Self-Gravity

We present a new code, CASTRO, that solves the multicomponent compressible hydrodynamic equations for astrophysical flows including self-gravity, nuclear reactions and radiation. CASTRO uses an Eulerian grid and incorporates adaptive mesh refinement (AMR). Our approach to AMR uses a nested hierarchy of logically-rectangular grids with simultaneous refinement in both space and time. The radiation component of CASTRO will be described in detail in the next paper, Part II, of this series.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

MAESTRO: An Adaptive Low Mach Number Hydrodynamics Algorithm for Stellar Flows

Many astrophysical phenomena are highly subsonic, requiring specialized numerical methods suitable for long-time integration. In a series of earlier papers we described the development of MAESTRO, a low Mach number stellar hydrodynamics code that can be used to simulate long-time, low-speed flows that would be prohibitively expensive to model using traditional compressible codes. MAESTRO is based on an equation set derived using low Mach number asymptotics; this equation set does not explicitly track acoustic waves and thus allows a significant increase in the time step. MAESTRO is suitable for two- and three-dimensional local atmospheric flows as well as three-dimensional full-star flows. Here, we continue the development of MAESTRO by incorporating adaptive mesh refinement (AMR). The primary difference between MAESTRO and other structured grid AMR approaches for incompressible and low Mach number flows is the presence of the time-dependent base state, whose evolution is coupled to the evolution of the full solution. We also describe how to incorporate the expansion of the base state for full-star flows, which involves a novel mapping technique between the one-dimensional base state and the Cartesian grid, as well as a number of overall improvements to the algorithm. We examine the efficiency and accuracy of our adaptive code, and demonstrate that it is suitable for further study of our initial scientific application, the convective phase of Type Ia supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multidimensional Modeling of Type I X-ray Bursts. I. Two-Dimensional Convection Prior to the Outburst of a Pure Helium Accretor

We present multidimensional simulations of the early convective phase preceding ignition in a Type I X-ray burst using the low Mach number hydrodynamics code, MAESTRO. A low Mach number approach is necessary in order to perform long-time integration required to study such phenomena. Using MAESTRO, we are able to capture the expansion of the atmosphere due to large-scale heating while capturing local compressibility effects such as those due to reactions and thermal diffusion. We also discuss the preparation of one-dimensional initial models and the subsequent mapping into our multidimensional framework. Our method of initial model generation differs from that used in previous multidimensional studies, which evolved a system through multiple bursts in one dimension before mapping onto a multidimensional grid. In our multidimensional simulations, we find that the resolution necessary to properly resolve the burning layer is an order of magnitude greater than that used in the earlier studies mentioned above. We characterize the convective patterns that form and discuss their resulting influence on the state of the convective region, which is important in modeling the outburst itself.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Low Mach Number Modeling of Type Ia Supernovae. IV. White Dwarf Convection

We present the first three-dimensional, full-star simulations of convection in a white dwarf preceding a Type Ia supernova, specifically the last few hours before ignition. For these long-time calculations we use our low Mach number hydrodynamics code, MAESTRO, which we have further developed to treat spherical stars centered in a three-dimensional Cartesian geometry. The main change required is a procedure to map the one-dimensional radial base state to and from the Cartesian grid. Our models recover the dipole structure of the flow seen in previous calculations, but our long-time integration shows that the orientation of the dipole changes with time. Furthermore, we show the development of gravity waves in the outer, stable portion of the star. Finally, we evolve several calculations to the point of ignition and discuss the range of ignition radii.

0 citations0 reviewsNo code linkedOpen access