Related papers: AstroGK: Astrophysical Gyrokinetics Code
We present a new massively parallel code for N-body and cosmological hydrodynamical simulations of modified gravity models. The code employs a multigrid-accelerated Newton-Gauss-Seidel relaxation solver on an adaptive mesh to efficiently…
High order accurate and explicit time-stable solvers are well suited for hyperbolic wave propagation problems. As a result of the complexities of real geometries, internal interfaces and nonlinear boundary and interface conditions,…
We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy…
We incorporate a time-independent gravitational field into the BGK scheme for numerical hydrodynamics. In the BGK scheme the gas evolves via an approximation to the collisional Boltzmann equation, namely the Bhatnagar-Gross-Krook (BGK)…
We present a numerical code for multi-component simulation of the galactic evolution. Our code includes the following parts: $N$-body is used to evolve dark matter, stellar dynamics and dust grains, gas dynamics is based on TVD-MUSCL scheme…
We discuss the development, analysis, implementation, and numerical assessment of a spectral method for the numerical simulation of the three-dimensional Vlasov-Maxwell equations. The method is based on a spectral expansion of the velocity…
In this paper, a unified framework for particulate two-phase flow will be presented with a wide range of solid-particle concentration from dilute to dense limit. The two phase flow is simulated by two coupled flow solvers, i.e., the…
Particle acceleration by turbulence plays a role in many astrophysical environments. The non- linear evolution of the underlying cosmic-ray spectrum is complex and can be described by a Fokker-Planck equation, which in general has to be…
Recent development of real-time equilibrium code Equinox [1] using a fixed-point algorithm [2] allow major plasma magnetic parameters to be identified in real-time, using rigorous analytical method. The code relies on the boundary flux code…
The DG algorithm is a powerful method for solving pdes, especially for evolution equations in conservation form. Since the algorithm involves integration over volume elements, it is not immediately obvious that it will generalize easily to…
We present a multiscale continuous Galerkin (MSCG) method for the fast and accurate stochastic simulation and optimization of time-harmonic wave propagation through photonic crystals. The MSCG method exploits repeated patterns in the…
Upcoming large scale telescope projects such as the Square Kilometre Array (SKA) will see high data rates and large data volumes; requiring tools that can analyse telescope event data quickly and accurately. In modern radio telescopes,…
Direct observations of compact objects, in the form of radiation spectra, gravitational waves from VIRGO/LIGO, and forthcoming direct imaging, are currently one of the primary source of information on the physics of plasmas in extreme…
The formation of structure in the Universe at large scales is dominated by gravity, with baryonic physics becoming significant at $\sim{\rm Mpc}$ scales. To capture the impact of baryonic physics, cosmological simulations must model gas…
A considerable amount of attention has been given to discontinuous Galerkin methods for hyperbolic problems in numerical relativity, showing potential advantages of the methods in dealing with hydrodynamical shocks and other…
We present COSMOCR, a numerical code for the investigation of cosmic ray related studies in computational cosmology. The code follows the diffusive shock acceleration, the mechanical and radiative energy losses and the spatial transport of…
Gyrokinetics is a rich and rewarding playground to study some of the mysteries of modern physics. In this thesis I present work, motivated by the quest for fusion energy, which seeks to uncover some of the inner workings of turbulence in…
Accurate modeling of core instabilities in tokamak plasmas is essential to understand the underlying physical mechanisms and their impact on plasma confinement. The ideal stability of the internal kink mode and the m = 1 collisionless…
In this paper, a unified algorithm will be proposed for the study of gas-solid particle multiphase flow. The gas-kinetic scheme (GKS) is used to simulate the continuum gas phase and the multiscale unified gas-kinetic wave-particle (UGKWP)…
The numerical simulation of turbulence and flows in almost ideal, large-Reynolds-number astrophysical plasmas motivates the implementation of almost conservative MHD computer codes. They should efficiently calculate, use highly parallelized…