Initial Data for the 3D Scalar Wave Equation
This thorn provides different initial data for the 3D scalar wave equation.
The general spherically symmetric solution can be written \begin {equation} \Psi (r,t) = \frac {1}{r}\left (f(r+t)+g(r-t)\right ) \end {equation} where the functions \(f\) and \(g\) can be freely chosen.
Making the additional requirement of time symmetry at \(t=0\), forces \begin {equation} f(r)=g(r) \end {equation} Thus if the solution at t=0 is given by \(\phi (r)\), the general solution will be \begin {equation} \Psi (r,t) = \frac {1}{2r}\left ( (r+t)\phi (r+t)+(r-t)\phi (r-t) \right ) \end {equation}
The gaussian solution is spherically symmetric about the origin of the Cartesian coordinate system, and is time symmetric. The initial profile is \begin {equation} \phi (r) = A \exp (- r^2/\sigma ) \end {equation} with the solution at the origin being \begin {equation} \Psi (r=0,t) = \left (1-2\frac {t^2}{\sigma }\right )\exp (-t^2/\sigma ) \end {equation}
The Gaussian solution is set with the parameters
amplitude = \(A\)
sigma = \(\sigma \)
amplitude | Scope: private | REAL |
Description: The amplitude of the waves
| ||
Range | Default: 1.0 | |
*:* | No restriction
| |
kx | Scope: private | REAL |
Description: The wave number in the x-direction
| ||
Range | Default: 4.0 | |
*:* | No restriction
| |
ky | Scope: private | REAL |
Description: The wave number in the y-direction
| ||
Range | Default: 0.0 | |
*:* | No restriction
| |
kz | Scope: private | REAL |
Description: The wave number in the z-direction
| ||
Range | Default: 0.0 | |
*:* | No restriction
| |
radius | Scope: private | REAL |
Description: The radius of the gaussian wave
| ||
Range | Default: 0.0 | |
0:* | Positive
| |
sigma | Scope: private | REAL |
Description: The sigma for the gaussian wave
| ||
Range | Default: 0.1 | |
0:* | Positive
| |
initial_data | Scope: restricted | KEYWORD |
Description: Type of initial data
| ||
Range | Default: gaussian | |
plane | Plane wave
| |
gaussian | Gaussian wave
| |
box | Box wave
| |
none | No initial data, zero phi
| |
type | Scope: shared from GRID | KEYWORD |
Implements:
idscalarwave
Inherits:
wavetoy
grid
This section lists all the variables which are assigned storage by thorn CactusWave/IDScalarWave. Storage can either last for the duration of the run (Always means that if this thorn is activated storage will be assigned, Conditional means that if this thorn is activated storage will be assigned for the duration of the run if some condition is met), or can be turned on for the duration of a schedule function.
NONE
CCTK_PARAMCHECK
idscalarwave_checkparameters
check parameters
Language: | fortran | |
Type: | function | |
CCTK_INITIAL
wavetoy_initialdata
execute initial data for wavetoy
Type: | group | |
WaveToy_InitialData
idscalarwave_initialdata
initial data for 3d wave equation
Language: | fortran | |
Storage: | wavetoy::scalarevolve[3] | |
Type: | function | |