wiki:AthenaDocsSMRGeom

Version 4 (modified by jmstone, 4 years ago) (diff)

--

Geometry of Levels

Wiki/Documentation/UserGuide to SMR?/Geometry

There are some important restrictions on the geometry of levels, as described below. At the beginning of execution, the code tests whether the Mesh hierarchy violates any of these restrictions: the code should quit with a warning message if a problem is found.

Restrictions on Domains touching the boundary of the root Domain

Fine (level>0) Domains may touch the boundaries of the root (level=0) Domain (see the examples below). In such cases, the physical boundary conditions applied to the root Domain will be applied to the finer level Domains as well.

However, with periodic BCs, the finer level Domains must touch both sides of the root Domain, or none at all. Thus, the first two examples above are not allowed with periodic BCs, while the third is allowed.

Space is required between Domains at different levels

The edges of Domains at different levels cannot touch, unless they also touch the boundary of the root Domain. In addition, there must be at least nghost (generally 4) cells on the finest level between the boundaries of Domains at different levels. This space is needed so that the ghost zones on the finer level can be set by a prolongation of active cells (not the ghost zones) on the coarser level. The first example below is allowed, and the following two are not allowed.

Domains at the same level cannot touch

In order to provide space for the ghost zones of each Domain to be set by prolongation, Domains at the same level cannot touch, or be closer than nghost (usually 4) cells on that Domain. If Domains touch, then the ghost zones can be set by mapping them to the corresponding active cells on the neighboring Domain, however this logic is complicated with face-centered fields in MHD and has not been implemented (yet). All of the examples below are not allowed.

No image "SMRDoms9.png" attached to AthenaDocsSMRGeom

Attachments