[Users] Meeting minutes

Roland Haas roland.haas at physics.gatech.edu
Mon Nov 21 12:06:32 CST 2011


Hello all,

> Constraint transport:
> * unigrid, cell centred implementation seems to work for the balsara
> shock tube tests as long as one is careful picking the parameters for
> the shock to make sure it is perfectly diagonal on the grid
> * Roland will re-run the test and send around the parameter files and
> sample output
> * still unsure about how many ghost zones are required
As promised attached please find parameter files and a plot of the
interesting output for the Balsara shock tube tests.

These are /similar/ but not identical to Bruno Mundom's original
balsara1d.par file (in GRHydro_initData/par). Please note that at least
the diagonal shock tube only runs on a single MPI processor (it produces
wrong results on more than one MPI process since it makes assumptions
about cctk_lsh and the location of the origin). Differences are: use TVD
instead of PPM, use ICN instead of RK2, use "flat" instead of "none",
use 4 ghost zones, put point at the origin, make domain cubical. None of
these should in principle be required (and they might not be, it was
just easier to test with them turned on since they make MoL much simpler
and detach the outer boundaries).

As of now the status is:
* any x,y,z shock has perfect divergence transport (divB=0 exactly)
* the diagonal shock (after checking it is diagonal on the grid),
develops small divergences on the order of divB=4e-12 in step 3. These
seem to be caused by roundoff. The expressions in the FORTRAN code
should give precisely zero. This I have tested with the attached Maple
worksheet which was created by copying and pasting the FORTRAN code then
making it digestible to Maple (in grhydro_cons.mw,
constraint_transport.mw is the beautified version used to check the
expressions first).

My current guess for the non-zero divB would be that (a) I have an
insidious bug, (b) divB is zero because of many numerical cancellations
(there are on the order of 16 terms contributing to it) so what I see is
simply loss of numerical accuracy.

In the attached image the upper half shows the x component of the B
field B^i and the lower half the divergence.

Yours,
Roland

-- 
My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://keys.gnupg.net.
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: balsara1d.par
Url: http://lists.einsteintoolkit.org/pipermail/users/attachments/20111121/45ae91d9/attachment-0004.pl 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: balsara1d0000.png
Type: image/png
Size: 69100 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20111121/45ae91d9/attachment-0001.png 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: balsara1d_x.par
Url: http://lists.einsteintoolkit.org/pipermail/users/attachments/20111121/45ae91d9/attachment-0005.pl 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: balsara1d_y.par
Url: http://lists.einsteintoolkit.org/pipermail/users/attachments/20111121/45ae91d9/attachment-0006.pl 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: balsara1d_z.par
Url: http://lists.einsteintoolkit.org/pipermail/users/attachments/20111121/45ae91d9/attachment-0007.pl 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: grhydro_cons.mw
Type: text/xml
Size: 22771 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20111121/45ae91d9/attachment-0002.xml 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: constraint_transport.mw
Type: text/xml
Size: 29225 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20111121/45ae91d9/attachment-0003.xml 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 189 bytes
Desc: OpenPGP digital signature
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20111121/45ae91d9/attachment-0001.bin 


More information about the Users mailing list