Related papers: Fully-Coupled Simulation of Cosmic Reionization. I…
We propose a locally one dimensional (LOD) finite difference method for multidimensional Riesz fractional diffusion equation with variable coefficients on a finite domain. The numerical method is second-order convergent in both space and…
Non-thermal (relativistic) particles are injected into the cosmos by structure formation shock waves, active galactic nuclei and stellar explosions. We present a suite of unigrid cosmological simulations (up to $2048^3$) using a two-fluid…
Reduced graphene oxide (rGO) exhibits strong anisotropic light absorption and high compatibility with photonic integrated chips, making it a promising material for implementing high performance onchip polarization selective devices. The…
Modeling and simulation of High Power Microwave (HPM) breakdown, a multiscale phenomenon, is computationally expensive and requires solving Maxwell's equations (EM solver) coupled with a plasma continuity equation (plasma solver). In this…
We present an interface between the (magneto-) hydrodynamics code PLUTO and the plasma simulation and spectral synthesis code CLOUDY. By combining these codes, we constructed a new photoionization hydrodynamics solver: The PLUTO-CLOUDY…
We perform three-dimensional simulations of homogeneous and inhomogeneous cosmologies via the coupling of a numerical relativity code for spacetime evolution and smoothed particle hydrodynamics (SPH) code. Evolution of a flat dust and…
A hybrid Maxwell solver for fully relativistic and electromagnetic (EM) particle-in-cell (PIC) codes is described. In this solver, the EM fields are solved in $k$ space by performing an FFT in one direction, while using finite difference…
The emergence of three-dimensional magneto-hydrodynamic (MHD) simulations of stellar atmospheres has sparked a need for efficient radiative transfer codes to calculate detailed synthetic spectra. We present RH 1.5D, a massively parallel…
Cosmological simulations of reionization often treat radiative transfer by solving for the monopole and dipoles of the intensity field and by making ansatz for the quadrupole moments to close the system of equations. We investigate the…
We present the Technicolor Dawn simulations, a suite of cosmological radiation-hydrodynamic simulations of the first 1.2 billion years. By modeling a spatially-inhomogeneous UVB on-the-fly with 24 frequencies and resolving dark matter halos…
The linear perturbation theory of inhomogeneous reionization (LPTR) has been developed as an analytical tool for predicting the global ionized fraction and large-scale power spectrum of ionized density fluctuations during reionization. In…
Disc winds play a crucial role in many accreting astrophysical systems across all scales. In accreting white dwarfs (AWDs) and active galactic nuclei (AGN), radiation pressure on spectral lines is a promising wind-driving mechanism.…
We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the intergalactic medium, to detail the growth of HII…
The morphology of the total gas and the neutral hydrogen (H I) distributions during the cosmological epoch of reionization can be quantified with Minkowski Functionals (MFs) of isodensity surfaces. We compute the MFs from the output of a…
We combine observational data on a dozen independent cosmic properties at high-$z$ with the information on reionization drawn from the spectra of distant luminous sources and the cosmic microwave background (CMB) to constrain the…
We have developed a method of solving for multi-species chemical reaction flows in non--equilibrium and self--consistently with the hydrodynamic equations in an expanding FLRW universe. The method is based on a backward differencing scheme…
The explicit quasi-monotonic conservative TVD scheme and numerical method for the solution of the gravitational MHD equations are developed. The 2D numerical code for the simulation of multidimensional selfgravitating MHD flows on the…
The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics.…
We present a new three-dimensional general-relativistic hydrodynamic evolution scheme coupled to dynamical spacetime evolutions which is capable of efficiently simulating stellar collapse, isolated neutron stars, black hole formation, and…
Understanding the detailed process of cosmic reionization is one of the remaining problems in astrophysics and cosmology. Here we construct a model of cosmic reionization that includes contribution from high-$z$ galaxies and active galactic…