PLUTO
4.4-patch2
|
Set boundary conditions. More...
#include "pluto.h"
Functions | |
static void | FlipSign (int, int, int *) |
void | Boundary (const Data *d, int idim, Grid *grid) |
void | OutflowBoundary (double ***q, RBox *box, int side) |
void | PeriodicBoundary (double ***q, RBox *box, int side) |
void | ReflectiveBoundary (double ***q, int s, int stag, RBox *box, int side) |
void | PolarAxisBoundary (const Data *d, RBox *box, int side) |
The Boundary() function sets both internal and physical boundary conditions on one or more sides of the computational domain. It is used to fill ghost zones of both cell-centered and face-centered data arrays.
The type of boundary conditions at the leftmost or rightmost side of a given grid is specified by the integers grid.lbound[dir]
or grid.rbound[dir]
, respectively. When this value is different from zero, the local processor borders the physical boundary and the admissible values for lbound
or rbound
are OUTFLOW, REFLECTIVE, AXISYMMETRIC, EQTSYMMETRIC, PERIODIC, SHEARING or USERDEF. Conversely, when this value is zero (internal boundary), two neighboring processors that share the same side need to fill ghost zones by exchanging data values. This step is done here only for parallel computations on static grids.
Predefined physical boundary conditions are handled by the following functions:
-->Before: proc #0 proc #1 proc #2 proc #3 |.......| |.......| |.......| |.......| A B C D E F G H --> PeriodicBoundary() --> After: proc #0 proc #1 proc #2 proc #3 |.......| |.......| |.......| |.......| H B B D D F F H
Set boundary conditions on one or more sides of the computational domain.
[in,out] | d | pointer to PLUTO Data structure containing the solution array (including centered and staggered fields) |
[in] | idim | specifies on which side(s) of the computational domain boundary conditions must be set. Possible values are
|
[in] | grid | pointer to grid structure. |
|
static |
Reverse the sign of vector components with respect to axis side. Depending on type, one needs to symmetrize or anti-symmetrize:
where (n) is the normal components, (p) and (t) are the transverse (or poloidal and toroidal for cylindrical and spherical coordinates) components.
If the radiation module is used, the transformation law for the radiation flux is the same as for the velocities.
[in] | side | boundary side |
[in] | type | boundary condition type |
[out] | vsign | an array of values (+1 or -1) giving the sign |
void OutflowBoundary | ( | double *** | q, |
RBox * | box, | ||
int | side | ||
) |
Impose zero-gradient boundary conditions on 'q' on the boundary side specified by 'side'. The input array 'q' must not represent the normal component of a staggered magnetic fied.
[in,out] | q | a 3D array requiring ghost zone filling |
[in] | box | pointer to a RBox structure defining the extent of the boundary region |
[in] | side | the side of the computational domain. |
void PeriodicBoundary | ( | double *** | q, |
RBox * | box, | ||
int | side | ||
) |
Implements periodic boundary conditions in serial mode or when one processor only handle the periodic direction.
[in,out] | q | a 3D array requiring ghost zone filling |
[in] | box | pointer to a RBox structure defining the extent of the boundary region |
[in] | side | the side of the computational domain. |
Boundary conditions on singular axis.
void ReflectiveBoundary | ( | double *** | q, |
int | s, | ||
int | stag, | ||
RBox * | box, | ||
int | side | ||
) |
Make symmetric (s = 1) or anti-symmetric (s = -1) profiles with respect to the boundary plane specified by box->side. The sign is set by the FlipSign() function.
[in,out] | q | a 3D array requiring ghost zone filling |
[in] | s | an integer taking only the values +1 (symmetric profile) or -1 (antisymmetric profile) |
[in] | stag | an integer taking the values 0 (centered with respect to the boundary) or 1 (staggered with respect to the boundary). |
[in] | box | pointer to a RBox structure defining the extent of the boundary region |
[in] | side | the side of the computational domain. |