Sandia National Laboratories Tramonto
printLeftMenu(); ?>

Mesh Parameters

These parameters set the size and mesh spacing of the computational domain. They also define boundary conditions.

  • Ndim(integer): The number of spatial dimensions to be solved numerically in the problem of interest.
  • Size_x[idim](real vector): The size of the computational domain. There is one entry for each dimension.
  • Esize_x[idim](real vector): The mesh spacing. There is one entry for each dimension (they need not be the same).
  • Type_bc[0,1][idim](integer array): There are three lines of input associated with Type_bc. Each line should contain two integers. The different lines designate the boundary conditions in the different dimensions (some may not be used). The two entries on a given line specify boundary conditions for the two sides of the domain in a given dimension. The options are:
    • -1: IN_WALL: semi-infinite surface : set density=0 when integrating beyond boundary.
    • 0: IN_BULK: constant bulk fluid: set density equal to the bulk density when integrating beyond boundary.
    • 1: PERIODIC: periodic boundary: set density at node N+k equal to the density at node k (Nodes run from 0 to N).
    • 2: REFLECT: reflective boundary: set density at node N+k equal to the density at node N-k.
    • 3: LAST_NODE: continuation boundary: set density at node N+k equal to the density at node N.
    • 4: LAST_NODE_RESTART: continuation boundary: set density at node N+k equal to the density at node N, but keep solution values obtained from a restart file for these nodes. This can be useful for pinning a free interface away from some kind of surface or inclusion.