Scaling relations between numerical simulations and physical systems they represent

Abstract

The dynamical equations describing the evolution of a physical system generally have a freedom in the choice of units, where different choices correspond to different physical systems that are described by the same equations. Since there are three basic physical units, of mass, length and time, there are up to three free parameters in such a rescaling of the units, $N_f\le 3$. In Newtonian hydrodynamics, e.g., there are indeed usually three free parameters, $N_f=3$. If, however, the dynamical equations contain a universal dimensional constant, such as the speed of light in vacuum $c$ or the gravitational constant $G$, then the requirement that its value remains the same imposes a constraint on the rescaling, which reduces its number of free parameters by one, to $N_f=2$. This is the case, for example, in magneto-hydrodynamics (MHD) or special relativistic hydrodynamics, where $c$ appears in the dynamical equations and forces the length and time units to scale by the same factor, or in Newtonian gravity where the gravitational constant $G$ appears in the equations. More generally, when there are $N_{udc}$ independent (in terms of their units) universal dimensional constants, then the number of free parameters is $N_f=max(0,3-N_{udc})$. When both gravity and relativity are included, there is only one free parameter ($N_f=1$, as both $G$ and $c$ appear in the equations so that $N_{udc}=2$), and the units of mass, length and time must all scale by the same factor. The explicit rescalings for different types of systems are discussed and summarized here. Such rescalings of the units also hold for discrete particles, e.g. in N-body or particle in cell simulations. They are very useful when numerically investigating a large parameter space or when attempting to fit particular experimental results, by significantly reducing the required number of simulations.

6 pages, 2 tables, accepted to MNRAS (expanded discussion of the general context in the introduction)

References in corpus (6)

• Particle acceleration in relativistic collisionless shocks: Fermi process at last?
• The Dynamics and Afterglow Radiation of Gamma-Ray Bursts. I. Constant Density Medium
• Gamma-ray burst afterglow broadband fitting based directly on hydrodynamics simulations
• Deceleration of arbitrarily magnetized GRB ejecta: the complete evolution
• Simulations of GRB Jets in a Stratified External Medium: Dynamics, Afterglow Lightcurves, Jet Breaks and Radio Calorimetry
• Gamma-Ray Burst Dynamics and Afterglow Radiation from Adaptive Mesh Refinement, Special Relativistic Hydrodynamic Simulations

Cited by in corpus (15)

• A cocoon shock breakout as the origin of the $\gamma$-ray emission in GW170817
• From $\gamma$ to Radio - The Electromagnetic Counterpart of GW 170817
• Off-Axis Emission of Short GRB Jets from Double Neutron Star Mergers and GRB 170817A
• Rayleigh-Taylor Instability in a Relativistic Fireball on a Moving Computational Grid
• Photospheric emission throughout GRB 100507 detected by Fermi
• Off-axis afterglow light curves and images from 2D hydrodynamic simulations of double-sided GRB jets in a stratified external medium
• PIC methods in astrophysics: Simulations of relativistic jets and kinetic physics in astrophysical systems
• On the opening angle of magnetised jets from neutron-star mergers: the case of GRB170817A
• What determines the structure of short gamma-ray burst jets?
• Radio limits on off-axis GRB afterglows and VLBI observations of SN 2003gk
• Numerical Simulations of an Initially Top-Hat Jet and the Afterglow of GW170817/GRB170817A
• Detectability of neutron star merger afterglows
• On the Deceleration and Spreading of Relativistic Jets I: Jet Dynamics
• Gamma-ray burst afterglow blast waves
• Temporal Evolution of Prompt GRB Polarization