blob: a18f46b69e67060088ec906383de444e493e78a5 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
|
% If you are using CVS use this line to give version information
% $Header$
\documentclass{article}
% Use the Cactus ThornGuide style file
% (Automatically used from Cactus distribution, if you have a
% thorn without the Cactus Flesh download this from the Cactus
% homepage at www.cactuscode.org)
\usepackage{../../../../doc/ThornGuide/cactus}
\begin{document}
\title{IDScalarWave}
\author{Gabrielle Allen \\ Horst Beyer}
\date{$ $Date$ $}
\maketitle
% Do not delete next line
% START CACTUS THORNGUIDE
\begin{abstract}
Initial Data for the 3D Scalar Wave Equation
\end{abstract}
\section{Purpose}
This thorn provides different initial data for the 3D scalar wave
equation.
\section{Spherically Symmetric Solutions}
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}
\section{Gaussian}
The gaussian solution is {\it spherically symmetric} about the
origin of the Cartesian coordinate system, and is {\it 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
\begin{itemize}
\item {\tt amplitude} = $A$
\item {\tt sigma} = $\sigma$
\end{itemize}
% Do not delete next line
% END CACTUS THORNGUIDE
\end{document}
|