Computational wave propagation
The aim is to model sound propagation in waveguides consisting of water, soft sediment, hard sediment or bedrock layers.
Our solution approach is to couple the solution of the (Helmholtz) differential equations (DEs) for the inhomogeneous fluid (water and soft sediment) layers to the solution of boundary integral formulations of the Helmholtz or Navier equations (BIEs) for the homogeneous fluid or solid layers. The gist of this hybrid method is to radically decrease the number of unknowns, since discretizations of BIEs involve only the boundary in contrast to DEs.
The figure shows how sound is damped (typically in the summer) in a slice of inhomogeneous water from the Baltic Sea. The computation has been performed using our Helmholtz solver, which is based on tensor notation and implemented in an object-oriented style.