radarfdtd.man

a full 3D simulation of electromagnetic waves with efficient absorbing boundary a full 3D simulation

.rn '' }`
''' $RCSfile$$Revision$$Date$
'''
''' $Log$
'''
.de Sh
.br
.if t .Sp
.ne 5
.PP
\fB\\$1\fR
.PP
..
.de Sp
.if t .sp .5v
.if n .sp
..
.de Ip
.br
.ie \\n(.$>=3 .ne \\$3
.el .ne 3
.IP "\\$1" \\$2
..
.de Vb
.ft CW
.nf
.ne \\$1
..
.de Ve
.ft R

.fi
..
'''
'''
'''     Set up \*(-- to give an unbreakable dash;
'''     string Tr holds user defined translation string.
'''     Bell System Logo is used as a dummy character.
'''
.tr \(*W-|\(bv\*(Tr
.ie n \{\
.ds -- \(*W-
.ds PI pi
.if (\n(.H=4u)&(1m=24u) .ds -- \(*W\h'-12u'\(*W\h'-12u'-\" diablo 10 pitch
.if (\n(.H=4u)&(1m=20u) .ds -- \(*W\h'-12u'\(*W\h'-8u'-\" diablo 12 pitch
.ds L" ""
.ds R" ""
'''   \*(M", \*(S", \*(N" and \*(T" are the equivalent of
'''   \*(L" and \*(R", except that they are used on ".xx" lines,
'''   such as .IP and .SH, which do another additional levels of
'''   double-quote interpretation
.ds M" """
.ds S" """
.ds N" """""
.ds T" """""
.ds L' '
.ds R' '
.ds M' '
.ds S' '
.ds N' '
.ds T' '
'br\}
.el\{\
.ds -- \(em\|
.tr \*(Tr
.ds L" ``
.ds R" ''
.ds M" ``
.ds S" ''
.ds N" ``
.ds T" ''
.ds L' `
.ds R' '
.ds M' `
.ds S' '
.ds N' `
.ds T' '
.ds PI \(*p
'br\}
.\"	If the F register is turned on, we'll generate
.\"	index entries out stderr for the following things:
.\"		TH	Title 
.\"		SH	Header
.\"		Sh	Subsection 
.\"		Ip	Item
.\"		X<>	Xref  (embedded
.\"	Of course, you have to process the output yourself
.\"	in some meaninful fashion.
.if \nF \{
.de IX
.tm Index:\\$1\t\\n%\t"\\$2"
..
.nr % 0
.rr F
.\}
.TH MANUAL 1 "perl 5.005, patch 03" "4/May/2000" "User Contributed Perl Documentation"
.UC
.if n .hy 0
.if n .na
.ds C+ C\v'-.1v'\h'-1p'\s-2+\h'-1p'+\s0\v'.1v'\h'-1p'
.de CQ          \" put $1 in typewriter font
.ft CW
'if n "\c
'if t \\&\\$1\c
'if n \\&\\$1\c
'if n \&"
\\&\\$2 \\$3 \\$4 \\$5 \\$6 \\$7
'.ft R
..
.\" @(#)ms.acc 1.5 88/02/08 SMI; from UCB 4.2
.	\" AM - accent mark definitions
.bd B 3
.	\" fudge factors for nroff and troff
.if n \{\
.	ds #H 0
.	ds #V .8m
.	ds #F .3m
.	ds #[ \f1
.	ds #] \fP
.\}
.if t \{\
.	ds #H ((1u-(\\\\n(.fu%2u))*.13m)
.	ds #V .6m
.	ds #F 0
.	ds #[ \&
.	ds #] \&
.\}
.	\" simple accents for nroff and troff
.if n \{\
.	ds ' \&
.	ds ` \&
.	ds ^ \&
.	ds , \&
.	ds ~ ~
.	ds ? ?
.	ds ! !
.	ds /
.	ds q
.\}
.if t \{\
.	ds ' \\k:\h'-(\\n(.wu*8/10-\*(#H)'\'\h"|\\n:u"
.	ds ` \\k:\h'-(\\n(.wu*8/10-\*(#H)'\`\h'|\\n:u'
.	ds ^ \\k:\h'-(\\n(.wu*10/11-\*(#H)'^\h'|\\n:u'
.	ds , \\k:\h'-(\\n(.wu*8/10)',\h'|\\n:u'
.	ds ~ \\k:\h'-(\\n(.wu-\*(#H-.1m)'~\h'|\\n:u'
.	ds ? \s-2c\h'-\w'c'u*7/10'\u\h'\*(#H'\zi\d\s+2\h'\w'c'u*8/10'
.	ds ! \s-2\(or\s+2\h'-\w'\(or'u'\v'-.8m'.\v'.8m'
.	ds / \\k:\h'-(\\n(.wu*8/10-\*(#H)'\z\(sl\h'|\\n:u'
.	ds q o\h'-\w'o'u*8/10'\s-4\v'.4m'\z\(*i\v'-.4m'\s+4\h'\w'o'u*8/10'
.\}
.	\" troff and (daisy-wheel) nroff accents
.ds : \\k:\h'-(\\n(.wu*8/10-\*(#H+.1m+\*(#F)'\v'-\*(#V'\z.\h'.2m+\*(#F'.\h'|\\n:u'\v'\*(#V'
.ds 8 \h'\*(#H'\(*b\h'-\*(#H'
.ds v \\k:\h'-(\\n(.wu*9/10-\*(#H)'\v'-\*(#V'\*(#[\s-4v\s0\v'\*(#V'\h'|\\n:u'\*(#]
.ds _ \\k:\h'-(\\n(.wu*9/10-\*(#H+(\*(#F*2/3))'\v'-.4m'\z\(hy\v'.4m'\h'|\\n:u'
.ds . \\k:\h'-(\\n(.wu*8/10)'\v'\*(#V*4/10'\z.\v'-\*(#V*4/10'\h'|\\n:u'
.ds 3 \*(#[\v'.2m'\s-2\&3\s0\v'-.2m'\*(#]
.ds o \\k:\h'-(\\n(.wu+\w'\(de'u-\*(#H)/2u'\v'-.3n'\*(#[\z\(de\v'.3n'\h'|\\n:u'\*(#]
.ds d- \h'\*(#H'\(pd\h'-\w'~'u'\v'-.25m'\f2\(hy\fP\v'.25m'\h'-\*(#H'
.ds D- D\\k:\h'-\w'D'u'\v'-.11m'\z\(hy\v'.11m'\h'|\\n:u'
.ds th \*(#[\v'.3m'\s+1I\s-1\v'-.3m'\h'-(\w'I'u*2/3)'\s-1o\s+1\*(#]
.ds Th \*(#[\s+2I\s-2\h'-\w'I'u*3/5'\v'-.3m'o\v'.3m'\*(#]
.ds ae a\h'-(\w'a'u*4/10)'e
.ds Ae A\h'-(\w'A'u*4/10)'E
.ds oe o\h'-(\w'o'u*4/10)'e
.ds Oe O\h'-(\w'O'u*4/10)'E
.	\" corrections for vroff
.if v .ds ~ \\k:\h'-(\\n(.wu*9/10-\*(#H)'\s-2\u~\d\s+2\h'|\\n:u'
.if v .ds ^ \\k:\h'-(\\n(.wu*10/11-\*(#H)'\v'-.4m'^\v'.4m'\h'|\\n:u'
.	\" for low resolution devices (crt and lpr)
.if \n(.H>23 .if \n(.V>19 \
\{\
.	ds : e
.	ds 8 ss
.	ds v \h'-1'\o'\(aa\(ga'
.	ds _ \h'-1'^
.	ds . \h'-1'.
.	ds 3 3
.	ds o a
.	ds d- d\h'-1'\(ga
.	ds D- D\h'-1'\(hy
.	ds th \o'bp'
.	ds Th \o'LP'
.	ds ae ae
.	ds Ae AE
.	ds oe oe
.	ds Oe OE
.\}
.rm #[ #] #H #V #F C
.SH "NAME"
radarFDTD \- a full 3D simulation of electromagnetic waves with efficient absorbing boundary conditions
.SH "SYNOPSIS"
radarFDTD directory [configuration file]
.SH "DESCRIPTION"
When started the program needs a directory as a command line parameter. If there is no second parameter defined it assumes the file is named sim.cfg. After starting it will output quite a few informations about the simulation it is going to start. The most important part is the section about memory, do not try to run the program in the swap area, since this will surely slow it down and may risk your hard disk which may start to smoke after a few weeks of hard work :\-(.
.SH "CONFIGURATION FILE"
In this section there will be a short description of all sections within the file sim.cfg. Feel free to experiment with those settings, but do not blame me, if something goes entirely wrong... ;\-)
.PP
Since I am not a native English speaker, I hope you are not offended if I left some mistakes in here.
.PP
Last thing, before we start: all lines which are empty or have a leading pound sign \*(L'#\*(R' are ignored (BTW, can someone explain to me why this double cross is called pound sign?).
.Sh "\s-1VERSION\s0 \s-1NUMBER\s0"
.PP
.Vb 2
\&  # version number is...
\&  3
.Ve
This just states the version number this configuration file was written for. This is checked by the program and if it does not match the internal layout, it will not use it.
.Sh "\s-1MEMORY\s0"
.PP
.Vb 2
\&  # maximum number of bytes the program may allocate
\&  100000000
.Ve
This is the number of bytes the program is allowed to allocate. Please change this according to your own need (or the memory limits of your computer; remember: this limit should be less than the physical amount of memory installed)
.Sh "\s-1MAXIMUM\s0 \s-1PML\s0 \s-1CONDUCTIVITY\s0"
.PP
.Vb 2
\&  # insert maximum pml conductivity here [S/m]
\&  0.025
.Ve
There is a fixed equation which rules the rise of conductivity. If you are interested in a more detailed explanation, have a look at the paper of Fang and Wu from 1996 (citation can be found in the \s-1README\s0 file). As a good first guess try something between 10 and 50 mS/m, which should be reasonable good.
.Sh "\s-1MAXIMUM\s0 \s-1STRETCHING\s0"
.PP
.Vb 2
\&  # insert maximum pml-stretching here [1]
\&  4.0
.Ve
Stretching should prevent evanescent waves to keep on moving within the perfectly matched layers, do not use too high values here, 2 to 4 seems to be quite reasonable.
.Sh "\s-1STEEPNESS\s0"
.PP
.Vb 2
\&  # insert steepness here
\&  2.0
.Ve
The higher the value for steepness is, the less reflection will occur at the boundary between the simulation space and the absorbing layer which is generally a good thing. But unfortunately a high value here means that there will be a high reflection before a wave manages to reach the layer, which is a \s-1PEC\s0 (perfectly conducting layer)
.Sh "\s-1SIMULATION\s0 \s-1LENGTH\s0"
.PP
.Vb 2
\&  # Length of Simulation [s]
\&  30.0e-9
.Ve
This just states the overall time duration of the simulation in seconds. Do not try to run a simulation too long, i.e. for a few seconds if your simulation space is only a few meters wide... it will take ages
.Sh "\s-1SIZE\s0 \s-1OF\s0 \s-1SIMULATION\s0 \s-1SPACE\s0"
.PP
.Vb 2
\&  # size of simulation space (number of boxes in x-, y- and z-direction, respectively
\&  41 41 41
.Ve
These three numbers specify how man cells are used in each direction (let us call them x,y and z for the time being). I usually use odd numbers, so I can place the transmitter in the very middle of the simulation box, but handle this as you like to
.Sh "\s-1CELL\s0 \s-1SIZES\s0"
.PP
.Vb 2
\&  # real size of a box [m] (x,y,z)
\&  0.025 0.025 0.025
.Ve
Here you specify how big (or small) a single Yee-cell should be, in this example this is about 2.5 cm (which is roughly an inch ;\-)). Be careful with this selection, since you should have at least 10 boxes per main wavelength you want to simulate, otherwise you might get strange wavy things out of the simulation.
.Sh "\s-1TIME\s0 \s-1STEP\s0 \s-1SIZE\s0"
.PP
.Vb 4
\&  # this number is a bit special: if positive it is the size of a time steps in seconds,
\&  # if negative, a matching time step is computed via the Courant-criterion
\&  # and finally divided by this number, i.e.