Numerical approximation of systems modeling large structure formation in cosmology

M.Colombeau

DOI: https://doi.org/10.48550/arXiv.0905.0230

2009-05-04

Abstract:Numerical approximations of two classical fluid dynamics systems modelling large structure formation in cosmology are proposed. These systems model nonrelativistic and relativistic fluids submitted to self-gravitation in an expanding background. They are obtained by an adaptation of an extension of the Godunov method, using delta wave projections, which was first introduced in Le Roux & al. for the system of pressureless fluid dynamics and which adapts to various systems of fluid dynamics.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: How to simulate the formation of large - scale structures in cosmology by numerical methods. Specifically, the paper focuses on the dynamic system modeling of non - relativistic and relativistic fluids driven by self - gravity in the context of cosmic expansion, and proposes an improved Godunov method to handle the numerical approximation of these systems.

### Main problem decomposition:

1. **Proposal of numerical method**:
- A numerical scheme based on the Godunov method is proposed, using the delta - wave projection technique to solve the classical fluid mechanics equations that describe the formation of large - scale structures.
- This method can handle situations where the velocity sign changes frequently, which is a common physical phenomenon in cosmology.

2. **Specific application of the model**:
- This method is applied to single - fluid and multi - fluid systems in Newtonian and relativistic frameworks.
- The interactions between dark matter, baryonic matter, and radiation are simulated, especially their behaviors during the equivalence and decoupling periods.

3. **Verification and testing**:
- Numerical tests are carried out through a series of two - dimensional toy models to verify the effectiveness of the proposed method.
- The test contents include the structure formation of super - horizon perturbations and sub - horizon perturbations at different expansion rates, as well as the Meszaros effect (the phenomenon that structures are frozen by rapid expansion).

### Summary of mathematical formulas:

- **Continuity equation** (Newtonian domain):
\[
\frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{u}) + 3H(t)\rho=0
\]

- **Euler equation** (Newtonian domain):
\[
\frac{\partial (\rho\mathbf{u})}{\partial t}+\nabla\cdot(\rho\mathbf{u}\otimes\mathbf{u})+\nabla p+\rho\nabla\Phi + 4H(t)\rho\mathbf{u}=0
\]

- **Poisson equation** (Newtonian domain):
\[
\Delta\Phi = 4\pi G\rho
\]

- **Equations in the relativistic framework**:
\[
\frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{u})+3H(t)(\rho + p/c^{2})=0
\]
\[
\frac{\partial (\rho\mathbf{u})}{\partial t}+\nabla\cdot(\rho\mathbf{u}\otimes\mathbf{u})+\nabla p+\rho\nabla\Phi + 4H(t)\rho\mathbf{u}=0
\]
\[
\Delta\Phi = 4\pi G(\rho + 3p/c^{2})
\]

### Key conclusions:

- Through the proposed numerical method, the formation process of large - scale structures in cosmology can be effectively simulated, especially in cases where different expansion rates and different components (such as dark matter, baryonic matter, radiation) are mixed.
- The numerical results are consistent with the classical cosmological perturbation theory, indicating that this method has high accuracy and reliability.

In conclusion, this paper solves the complex dynamic problem of large - scale structure formation in the context of cosmic expansion by proposing and verifying a new numerical method.