Related papers: Numerical approaches for multidimensional simulati…
Turbulence cascade has been modeled using various methods; the one we have used applies to a more exact representation of turbulence where people use the multifractal representation. The nature of the energy dissipation is usually governed…
We present a new numerical scheme for solving the advection equation and its application to Vlasov simulations. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent…
We model the structure and evolution of black hole accretion disks, and their neighboring regions, using numerical simulations. The numerics is governed by the equations of general relativistic magneto-hydrodynamics (GRMHD). In particular,…
A method for symbolically computing conservation laws of nonlinear partial differential equations (PDEs) in multiple space dimensions is presented in the language of variational calculus and linear algebra. The steps of the method are…
Real-world applications of computational fluid dynamics often involve the evaluation of quantities of interest for several distinct geometries that define the computational domain or are embedded inside it. For example, design optimization…
The early evolution of dense star clusters is possibly dominated by close interactions between stars, and physical collisions between stars may occur quite frequently. Simulating a stellar collision event can be an intensive numerical task,…
Radiation hydrodynamics simulations based on the one-fluid two-temperature model may violate the law of energy conservation because the governing equations are expressed in a nonconservative formulation. Here, we maintain the important…
We propose efficient numerical algorithms for approximating statistical solutions of scalar conservation laws. The proposed algorithms combine finite volume spatio-temporal approximations with Monte Carlo and multi-level Monte Carlo…
Aligning theoretical atomistic structural models of materials with available experimental data presents a significant challenge for disordered systems. The configurational space to navigate is vast, and faithful realizations require large…
We present new two-dimensional (2D) axisymmetric neutrino radiation/hydrodynamic models of core-collapse supernova (CCSN) cores. We use the CASTRO code, which incorporates truly multi-dimensional, multi-group, flux-limited diffusion (MGFLD)…
An Eulerian TVD code and a Lagrangian SPH code are used to simulate the off-axis collision of equal-mass main sequence stars in order to address the question of whether stellar mergers can produce a remnant star where the interior has been…
We have developed a new simulation code aimed at studying the stellar dynamics of a galactic central star cluster surrounding a massive black hole. In order to include all the relevant physical ingredients (2-body relaxation, stellar mass…
In the present work we compare reliability of several most widely used reduced detailed chemical kinetic schemes for hydrogen-air and hydrogen-oxygen combustible mixtures. The validation of the schemes includes detailed analysis of 0D and…
Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially…
The set-up of the initial conditions in cosmological N-body simulations is usually implemented by rescaling the desired low-redshift linear power spectrum to the required starting redshift consistently with the Newtonian evolution of the…
The failure of heterogeneous materials with microstructures is a complex process of damage nucleation, growth and localisation. This process spans multiple length scales and is challenging to simulate numerically due to its high…
Using the new state-of-the-art core-collapse supernova (CCSN) code F{\sc{ornax}}, we have simulated the three-dimensional dynamical evolution of the cores of 9-, 10-, 11-, 12-, and 13-M$_{\odot}$ stars from the onset of collapse. Stars from…
We present new fully self-consistent models of the formation and evolution of isolated dwarf galaxies. We have used the publicly available N-body/SPH code HYDRA, to which we have added a set of star formation criteria, and prescriptions for…
We present a technique for optimizing hundreds of thousands of variational parameters in variational quantum Monte Carlo. By introducing iterative Krylov subspace solvers and by multiplying by the Hamiltonian and overlap matrices as they…
Neutrino transport in spherically symmetric models of stellar core collapse and bounce has achieved a technically complete level, rewarded by the agreement among independent groups that a multi-dimensional treatment of the…