Related papers: An Adaptive Code For Radial Stellar Pulsations
We introduce numerical algorithms for initializing multidimensional simulations of stellar explosions with 1D stellar evolution models. The initial mapping from 1D profiles onto multidimensional grids can generate severe numerical…
We present a new numerical code, PLUTO, for the solution of hypersonic flows in 1, 2 and 3 spatial dimensions and different systems of coordinates. The code provides a multi-physics, multi-algorithm modular environment particularly oriented…
Context: We study the impact of two-dimensional spherical shells on compressible convection. Realistic profiles for density and temperature from a one-dimensional stellar evolution code are used to produce a model of a large stellar…
Stars may be understood as self-gravitating masses of a compressible fluid whose radiative cooling is compensated by nuclear reactions or gravitational contraction. The understanding of their time evolution requires the use of detailed…
An adaptive moving mesh finite element method is proposed for the numerical solution of the regularized long wave (RLW) equation. A moving mesh strategy based on the so-called moving mesh PDE is used to adaptively move the mesh to improve…
The document describes a numerical algorithm for plasms and fluids by the Lagrange method, in which the spatial meshes follow the plasma and fluid behavior. Through the mesh wall the plasma and the fluid do not escape. The 3D Lagrange code…
We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR…
We present a new magnetohydrodynamic (MHD) simulation code with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements…
We describe a new scheme for evolving the equations of force-free electrodynamics, the vanishing-inertia limit of magnetohydrodynamics. This pseudospectral code uses global orthogonal basis function expansions to take accurate spatial…
We present a review of the TORUS radiation transfer and hydrodynamics code. TORUS uses a 1-D, 2-D or3-D adaptive mesh refinement scheme to store and manipulate the state variables, and solves the equation of radiative transfer using Monte…
The study of the interaction between solid objects and magnetohydrodynamic (MHD) fluids is of great importance in physics as consequence of the significant phenomena generated, such as planets interacting with stellar wind produced by their…
Many astrophysical explosions, such as type Ia supernovae, classical novae, and X-ray bursts, are dominated by thermonuclear runaway. To model these processes accurately, one must evolve nuclear reactions concurrently with hydrodynamics. We…
Compressible magnetohydrodynamic (MHD) turbulence is ubiquitous in astrophysical phenomena ranging from the intergalactic to the stellar scales. In studying them, numerical simulations are nearly inescapable, due to the large degree of…
We construct a new relativistic viscous hydrodynamics code optimized in the Milne coordinates. We split the conservation equations into an ideal part and a viscous part, using the Strang spitting method. In the code a Riemann solver based…
A one-dimensional turbulent convection model in the form of a time-dependent diffusion equation for the turbulent energy is incorporated into our numerical pulsation code. The effect of turbulent convection on the structural rearrangement…
We present an analysis of the response of a radiative region to waves generated by a convective region of the star; this wave treatment of the classical problem of ``overshooting'' gives extra mixing relative to the treatment traditionally…
We present the first magnetohydrodynamic model of the stellar chromospheric heating and acceleration of the outer atmospheres of cool evolved stars, using alpha Tau as a case study. We used a 1.5D MHD code with a generalized Ohm's law that…
This paper develops entropy stable (ES) adaptive moving mesh schemes for the 2D and 3D special relativistic hydrodynamic (RHD) equations. They are built on the ES finite volume approximation of the RHD equations in curvilinear coordinates,…
We present an object-oriented geophysical and astrophysical spectral-element adaptive refinement (GASpAR) code for application to turbulent flows. Like most spectral-element codes, GASpAR combines finite-element efficiency with…
We introduce a new computational method for embedding Lagrangian sink particles into an Eulerian calculation. Simulations of gravitational collapse or accretion generally produce regions whose density greatly exceeds the mean density in the…