Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe timestep constraint in solving the gravitational wave equations sourced by hydromagnetic turbulence
59
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Hydromagnetic turbulence produced during phase transitions in the early universe can be a powerful source of stochastic gravitational waves (GWs). GWs can be modelled by the linearised spatial part of the Einstein equations sourced by the Reynolds and Maxwell stresses. We have implemented two different GW solvers into the {sc Pencil Code} -- a code which uses a third order timestep and sixth order finite differences. Using direct numerical integration of the GW equations, we study the appearance of a numerical degradation of the GW amplitude at the highest wavenumbers, which depends on the length of the timestep -- even when the Courant--Friedrichs--Lewy condition is ten times below the stability limit. This degradation leads to a numerical error, which is found to scale with the third power of the timestep. A similar degradation is not seen in the magnetic and velocity fields. To mitigate numerical degradation effects, we alternatively use the exact solution of the GW equations under the assumption that the source is constant between subsequent timesteps. This allows us to use a much longer timestep, which cuts the computational cost by a factor of about ten.
rate research
Read More
We study the evolution of kinetic and magnetic energy spectra in magnetohydrodynamic flows in the presence of strong cross helicity. For forced turbulence, we find weak inverse transfer of kinetic energy toward the smallest wavenumber. This is plausibly explained by the finiteness of scale separation between the injection wavenumber and the smallest wavenumber of the domain, which here is a factor of 15. In the decaying case, there is a slight increase at the smallest wavenumber, which is probably explained by the dominance of kinetic energy over magnetic energy at the smallest wavenumbers. Within a range of wavenumbers covering almost an order of magnitude the decay is purely exponential, which is argued to be a consequence of a suppression of nonlinearity due to the presence of strong cross helicity.
We consider a generic dispersive massive gravity theory and numerically study its resulting modified energy and strain spectra of tensor gravitational waves (GWs) sourced by (i) fully developed turbulence during the electroweak phase transition (EWPT) and (ii) forced hydromagnetic turbulence during the QCD phase transition (QCDPT). The GW spectra are then computed in both spatial and temporal Fourier domains. We find, from the spatial spectra, that the slope modifications are weakly dependent on the eddy size at QCDPT, and, from the temporal spectra, that the modifications are pronounced in the $1$--$10{rm nHz}$ range -- the sensitivity range of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) -- for a graviton mass $m_{rm g}$ in the range $2times10^{-23}{rm eV}lesssim m_{rm g}c^2lesssim7times10^{-22}{rm eV}$.
A general Hamiltonian wave system with quartic resonances is considered, in the standard kinetic limit of a continuum of weakly interacting dispersive waves with random phases. The evolution equation for the multimode characteristic function $Z$ is obtained within an interaction representation and a perturbation expansion in the small nonlinearity parameter. A frequency renormalization is performed to remove linear terms that do not appear in the 3-wave case. Feynman-Wyld diagrams are used to average over phases, leading to a first order differential evolution equation for $Z$. A hierarchy of equations, analogous to the Boltzmann hierarchy for low density gases is derived, which preserves in time the property of random phases and amplitudes. This amounts to a general formalism for both the $N$-mode and the 1-mode PDF equations for 4-wave turbulent systems, suitable for numerical simulations and for investigating intermittency.
The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $lesssim 40%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.
We report results of sumulation of wave turbulence. Both inverse and direct cascades are observed. The definition of mesoscopic turbulence is given. This is a regime when the number of modes in a system involved in turbulence is high enough to qualitatively simulate most of the processes but significantly smaller then the threshold which gives us quantitative agreement with the statistical description, such as kinetic equation. Such a regime takes place in numerical simulation, in essentially finite systems, etc.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
