This page here has been created for the latest stable release of Salvus. You have chosen to view the documentation for another Salvus version. Please be aware there might be some small differences you have to account for.

Boundary Conditions

Without additional boundary conditions that constrain the wavefield at the boundaries of the computational domain, the wave equation does not obey a unique solution.

All boundary conditions are defined on so-called side sets of the mesh. A side set is just a way to flag edges (in 2D) or faces (in 3D) of a mesh and to give these a name. This is important, whenever different conditions are applied on different parts of the boundary of the domain.

The example mesh below has 4 side sets, each denoting one side of the mesh. The actual names are arbitrary and mesh dependent.

Common defaults for labelling the boundaries are

  • for Cartesian meshes: x0, x1, y0, y1, and additionally z0, z1 in 3D,
  • for spherical meshes: r0 for the inner radius (if present), r1 for the outer radius,
  • for spherical chunks: r0, r1 p0, p1, t0, t1.

In the following, Γ\Gamma denotes a side set where a certain boundary condition is imposed.

Without explicitly specifying boundary conditions, the spectral-element method will fall back to the natural boundary conditions.