Related papers: Numerical approaches for multidimensional simulati…
A numerical method is presented for first-principle simulations of charged colloidal dispersions in electrolyte solutions. Utilizing a smoothed profile for colloid-solvent boundaries, efficient mesoscopic simulations are enabled for…
With the advent of large spectroscopic surveys the amount of high quality chemo-dynamical data in the Milky Way (MW) increased tremendously. Accurately and correctly capturing and explaining the detailed features in the high-quality…
The gravitational collapse, bounce, the explosion of an iron core of an 11.2 $M_{\odot}$ star is simulated by two-dimensional neutrino-radiation hydrodynamic code. The explosion is driven by the neutrino heating aided by multi-dimensional…
In this paper, we develop a numerical scheme to handle interfaces across computational domains in multi-block schemes for the approximation of systems of conservation laws. We are interested in transmitting shock discontinuities without…
There are many astrophysical and laboratory scenarios where kinetic effects play an important role. These range from astrophysical shocks and plasma shell collisions, to high intensity laser-plasma interactions, with applications to fast…
Numerical schemes for wave-like systems with small dissipation are often inaccurate and unstable due to truncation errors and numerical roundoff errors. Hence, numerical simulations of wave-like systems lacking proper handling of these…
A major weakness in one-dimensional (1D) stellar structure and evolution modeling is the simplified treatment of convection, which leads to erroneous near-surface stratification and considerable uncertainties in predicted effective…
We here focus on the behavior of supernovae that technically explode in 1D (spherical symmetry). When simulated in 3D, however, the outcomes of representative progenitors of this class are quite different in almost all relevant quantities.…
We present preliminary results from two dimensional numerical studies of pair instability supernova (PSN). We study nuclear burning, hydrodynamic instabilities and explosion of very massive stars. Use a new radiation-hydrodynamics code,…
Numerical studies of core-collapse supernovae have demonstrated the importance of non-radial motions in pre-collapse progenitors on the explosion outcome. We use the CHIMERA neutrino radiation hydrodynamics code running seven…
We describe numerical methods for incorporating gas dynamics into cosmological simulations and present illustrative applications to the cold dark matter (CDM) scenario. Our evolution code, a version of TreeSPH (Hernquist \& Katz 1989)…
Computational astrochemical models are essential for helping us interpret and understand the observations of different astrophysical environments. In the age of high-resolution telescopes such as JWST and ALMA, the substructure of many…
Numerical simulations have become a necessary tool to describe the complex interactions among the different processes involved in galaxy formation and evolution, unfeasible via an analytic approach. The last decade has seen a great effort…
A new energy and enstrophy conserving scheme is evaluated using a suite of test cases over the global spherical domain or bounded domains. The evaluation is organized around a set of pre-defined properties: accuracy of individual opeartors,…
Hydrodynamical simulations are the most accurate way to model structure formation in the universe, but they often involve a large number of astrophysical parameters modeling subgrid physics, in addition to cosmological parameters. This…
Simulations inform all aspects of modern astrophysical research, ranging in scale from 1D and 2D test problems that can run in seconds on an astronomer's laptop all the way to large-scale 3D calculations that run on the largest…
It is well known that numerical errors grow exponentially in $N$-body simulations of gravitational bound stellar systems, but it is not well understood how the accuracy parameters of algorithms affect the physical evolution in simulations.…
We perform two- (2D) and three-dimensional (3D) hydrodynamics simulations of convective oxygen shell-burning that takes place deep inside a massive progenitor star of a core-collapse supernova. Using one dimensional (1D) stellar evolution…
A multi-domain spectral method for computing very high precision 3-D stellar models is presented. The boundary of each domain is chosen in order to coincide with a physical discontinuity (e.g. the star's surface). In addition, a…
We investigate how the accuracy and stability of numerical relativity simulations of 1D colliding plane waves depends on choices of equation formulations, gauge conditions, boundary conditions, and numerical methods, all in the context of a…