[Users] Problem with CarpetRegrid2/AMR
Hal Finkel
hfinkel at anl.gov
Tue Aug 30 14:37:20 CDT 2011
I am trying to setup a simple simulation using CarpetRegrid2's AMR
capability. As Erik has already pointed out to me (off list), the
adaptive-refinement triggering code is still under development.
Nevertheless, I would really like to get it working.
In my current test case, I am using Kranc to solve a Klein-Gordon
equation (with a non-linear potential). With the provided initial
conditions, the system has an expanding bubble. This works well when
refinement is turned off. To be clear, the evolution equation is:
evolCalc =
{
Name -> "SFBubble_Evolve",
Schedule -> {"in MoL_CalcRHS"},
Where -> Interior,
Shorthands -> {},
Equations -> {
dot[chi] -> chiM,
dot[chiM] -> PD[chi,li,lj] Euc[ui,uj] - Vchi[chi]
}
};
where Vchi[chi_] := V'[chi] (V is the potential). The AMR level_mask is
set by:
lvlmskSet =
{
Name -> "SFBubble_SetLevelMask",
Schedule -> {"AT preregrid"},
Where -> Interior,
Shorthands -> { adchix, adchiy, adchiz },
Equations ->
{
adchix -> fabs[dx PD[chi,1]],
adchiy -> fabs[dy PD[chi,2]],
adchiz -> fabs[dz PD[chi,3]],
"level_mask" -> (
Max[
adchix, adchiy, adchiz
]/dchimax
) /. Max[a_, b_, c_] -> fmax[a, fmax[b, c]]
}
};
(in words: the value of chi should not change by more than dchimax over
the length of one grid cell)
I am using the Periodic thorn to provide periodic boundary conditions.
I've attached six pictures. They are at t = 0, the first level 1
intermediate time step, and the first coarse time step. The first three
pictures are the value of chi on a slice through the middle of the box.
The second three pictures are the value of level_mask. As can be seen in
the second chi image, the interpolation into the refined region seems to
be incorrect. The "boundary" seems to have appeared in the middle of the
block. I have also attached the block descriptions. Is there something I
could be doing incorrectly which would cause this? Otherwise, does
anyone have any ideas on where to look for the problem?
I am using the trunk version as of this morning.
Thank you in advance,
Hal
--
Hal Finkel
Postdoctoral Appointee
Leadership Computing Facility
Argonne National Laboratory
-------------- next part --------------
A non-text attachment was scrubbed...
Name: visit0011.png
Type: image/png
Size: 92315 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20110830/820e1008/attachment-0006.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: visit0010.png
Type: image/png
Size: 107243 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20110830/820e1008/attachment-0007.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: visit0009.png
Type: image/png
Size: 25568 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20110830/820e1008/attachment-0008.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: visit0008.png
Type: image/png
Size: 88644 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20110830/820e1008/attachment-0009.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: visit0007.png
Type: image/png
Size: 86746 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20110830/820e1008/attachment-0010.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: visit0006.png
Type: image/png
Size: 75319 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20110830/820e1008/attachment-0011.png
-------------- next part --------------
# grid coordinates
# format: map reflevel region mglevel bounding-box
iteration 0
maps 1
0 mglevels 1
0 0 reflevels 1
0 0 0 regions 1
0 0 0 0 ([-24,-24,-24]:[216,216,216]:[8,8,8]/[-3,-3,-3]:[27,27,27]/[31,31,31]/29791)
iteration 1
maps 1
0 mglevels 1
0 0 reflevels 2
0 0 0 regions 1
0 0 0 0 ([-24,-24,-24]:[216,216,216]:[8,8,8]/[-3,-3,-3]:[27,27,27]/[31,31,31]/29791)
0 0 1 regions 8
0 0 1 0 ([64,164,64]:[128,208,128]:[4,4,4]/[16,41,16]:[32,52,32]/[17,12,17]/3468)
0 0 1 1 ([64,64,164]:[128,128,208]:[4,4,4]/[16,16,41]:[32,32,52]/[17,17,12]/3468)
0 0 1 2 ([64,32,32]:[160,60,60]:[4,4,4]/[16,8,8]:[40,15,15]/[25,8,8]/1600)
0 0 1 3 ([32,132,64]:[160,160,128]:[4,4,4]/[8,33,16]:[40,40,32]/[33,8,17]/4488)
0 0 1 4 ([32,32,64]:[160,60,128]:[4,4,4]/[8,8,16]:[40,15,32]/[33,8,17]/4488)
0 0 1 5 ([32,64,64]:[208,128,128]:[4,4,4]/[8,16,16]:[52,32,32]/[45,17,17]/13005)
0 0 1 6 ([32,64,32]:[160,160,60]:[4,4,4]/[8,16,8]:[40,40,15]/[33,25,8]/6600)
0 0 1 7 ([32,32,132]:[160,160,160]:[4,4,4]/[8,8,33]:[40,40,40]/[33,33,8]/8712)
-------------- next part --------------
# grid structure
# format: map reflevel component mglevel processor bounding-box is-outer-boundary
iteration 0
maps 1
0 mglevels 1
0 0 reflevels 1
0 0 0 components 1
0 0 0 0 0 ([0,0,0]:[15360,15360,15360]:[512,512,512]/[0,0,0]:[30,30,30]/[31,31,31]/29791) [[1,1,1],[1,1,1]]
iteration 1
maps 1
0 mglevels 1
0 0 reflevels 2
0 0 0 components 1
0 0 0 0 0 ([0,0,0]:[15360,15360,15360]:[512,512,512]/[0,0,0]:[30,30,30]/[31,31,31]/29791) [[1,1,1],[1,1,1]]
0 0 1 components 8
0 0 1 0 0 ([5632,12032,5632]:[9728,14848,9728]:[256,256,256]/[22,47,22]:[38,58,38]/[17,12,17]/3468) [[0,0,0],[0,1,0]]
0 0 1 1 0 ([5632,5632,12032]:[9728,9728,14848]:[256,256,256]/[22,22,47]:[38,38,58]/[17,17,12]/3468) [[0,0,0],[0,0,1]]
0 0 1 2 0 ([5632,3584,3584]:[11776,5376,5376]:[256,256,256]/[22,14,14]:[46,21,21]/[25,8,8]/1600) [[0,0,0],[0,0,0]]
0 0 1 3 0 ([3584,9984,5632]:[11776,11776,9728]:[256,256,256]/[14,39,22]:[46,46,38]/[33,8,17]/4488) [[0,0,0],[0,0,0]]
0 0 1 4 0 ([3584,3584,5632]:[11776,5376,9728]:[256,256,256]/[14,14,22]:[46,21,38]/[33,8,17]/4488) [[0,0,0],[0,0,0]]
0 0 1 5 0 ([3584,5632,5632]:[14848,9728,9728]:[256,256,256]/[14,22,22]:[58,38,38]/[45,17,17]/13005) [[0,0,0],[1,0,0]]
0 0 1 6 0 ([3584,5632,3584]:[11776,11776,5376]:[256,256,256]/[14,22,14]:[46,46,21]/[33,25,8]/6600) [[0,0,0],[0,0,0]]
0 0 1 7 0 ([3584,3584,9984]:[11776,11776,11776]:[256,256,256]/[14,14,39]:[46,46,46]/[33,33,8]/8712) [[0,0,0],[0,0,0]]
More information about the Users
mailing list