`specfem::assembly::boundaries_impl::acoustic_free_surface`¶

template<specfem::element::dimension_tag DimensionTag> struct acoustic_free_surface¶

Primary template declaration for acoustic free surface boundary conditions.

Data required to compute acoustic free surface boundary conditions.

General template declaration for acoustic free surface boundary conditions.

Dimension-Specific Implementations¶

template<> struct acoustic_free_surface<specfem::element::dimension_tag::dim2>¶

Data container used to store acoustic free surface boundary information.

This class provides data structures and data access functions required to implement acoustic free surface boundary conditions in 2D spectral element simulations.

Physical Background: An acoustic free surface boundary condition enforces zero traction (pressure) at the boundary. Mathematically, this corresponds to setting the acceleration to zero at quadrature points located on the free surface. This is physically appropriate for:

Air-water interfaces
Fluid-vacuum boundaries
Ocean-atmosphere interfaces in seismic modeling

Implementation Details: The class stores boundary tag information for each quadrature point on elements that contain acoustic free surface boundaries. The boundary tags are used during the assembly process to identify which quadrature points require free surface treatment.

Type Definitions

using BoundaryTagView = Kokkos::View<specfem::element::boundary_tag_container***, Kokkos::LayoutLeft, Kokkos::DefaultExecutionSpace>¶

Kokkos view type for storing boundary tag containers on device.

This view stores boundary tag information for all quadrature points within elements that contain acoustic free surface boundaries. The view has dimensions [nspec, ngllz, ngllx] where nspec is the number of spectral elements with free surface boundaries.

Public Constants

static constexpr auto dimension_tag = specfem::element::dimension_tag::dim2 ¶: Dimension tag indicating this is a 2D implementation.

Data Members

BoundaryTagView quadrature_point_boundary_tag¶

Device-accessible boundary tag information for quadrature points.

This view contains boundary tag containers for every quadrature point within spectral elements that have acoustic free surface boundaries. The boundary tags are used during device kernel execution to identify which points require free surface treatment.

Dimensions: [nspec_acoustic_free_surface, ngllz, ngllx]

nspec_acoustic_free_surface: Number of elements with free surface boundaries
ngllz, ngllx: Number of GLL points in z and x directions respectively

BoundaryTagView::host_mirror_type h_quadrature_point_boundary_tag¶

Host-accessible mirror of boundary tag information.

This is the host mirror of quadrature_point_boundary_tag, used for CPU access during initialization, debugging, and host-side computations. Data is synchronized between host and device as needed.

Constructors

acoustic_free_surface() = default¶

Default constructor.

Creates an empty acoustic free surface boundary condition container. Data members are left uninitialized and must be set up via assignment or the parametrized constructor.

acoustic_free_surface(const int nspec, const int ngllz, const int ngllx, const specfem::mesh::acoustic_free_surface<dimension_tag> &acoustic_free_surface, const specfem::assembly::mesh<dimension_tag> &mesh, const Kokkos::View<int*, Kokkos::HostSpace> &boundary_index_mapping, std::vector<specfem::element::boundary_tag_container> &boundary_tag)¶

Construct acoustic free surface boundary condition data from mesh information.

This constructor processes mesher-supplied boundary information and converts it into per-quadrature-point data required for implementing acoustic free surface boundary conditions during SEM simulations.

Parameters:

nspec – Number of spectral elements with acoustic free surface boundaries
ngllz – Number of GLL points in the z (vertical) direction
ngllx – Number of GLL points in the x (horizontal) direction
acoustic_free_surface – Mesh-level acoustic free surface boundary information containing element lists and edge definitions
mesh – Assembly mesh containing coordinate and connectivity information
boundary_index_mapping – Mapping from global spectral element indices to boundary-specific indices for elements with free surface boundaries
boundary_tag – Vector of boundary tag containers that will be populated with per-element boundary information

Device Data Access Methods

template<typename IndexType, typename std::enable_if_t<specfem::data_access::is_index_type<IndexType>::value && specfem::data_access::is_point<IndexType>::value && IndexType::using_simd == false, int> = 0> inline void load_on_device(const IndexType &index, specfem::point::boundary<boundary_tag, dimension_tag, false> &boundary) const¶

Load acoustic free surface boundary data for a non-SIMD quadrature point on device.

Note

This is an implementation detail and is typically called by a higher-level load_on_device function

Template Parameters:

IndexType – Must be a valid index type (specfem::point::index) with SIMD disabled

Parameters:

index – Quadrature point location specifier containing
boundary – Output boundary object

template<typename IndexType, typename std::enable_if_t<specfem::data_access::is_index_type<IndexType>::value && specfem::data_access::is_point<IndexType>::value && IndexType::using_simd == false, int> = 0> inline void load_on_device(const IndexType &index, specfem::point::boundary<specfem::element::boundary_tag::composite_stacey_dirichlet, dimension_tag, false> &boundary) const¶

Load acoustic free surface data for composite boundary conditions on device.

specfem::assembly::boundaries_impl::acoustic_free_surface¶

Dimension-Specific Implementations¶

`specfem::assembly::boundaries_impl::acoustic_free_surface`¶