Search

This paper describes a new morpho-dynamic 3D surf-zone model. The hydrodynamic model used for this purpose is the MIKE 3 Wave FM Model. This model solves the 3D Navier-Stokes equations on a 2D horizontal unstructured mesh with a sigma transformation describing the vertical coordinate. A recent solution to the well-known over-production of turbulence outside the surf-zone in RANS models has been applied and is shown to improve the predicted undertow profiles. The sediment transport is modelled using an advection-diffusion equation for the suspended load, together with a bed-load formulation. In the present paper the predicted surf zone hydrodynamics and sea bed evolutions are compared against experimental laboratory data.

Motivation In the pursuit of higher resolution simulation models that use all seismic, geological, and dynamic reservoir data and to make use of modern parallel computing architectures we consider alternative numerical methods to solve the system of equations governing subsurface porous media flow. It is standard in conventional techniques to use a global linearization in a Newton-type method to solve the strongly nonlinear system of equations arising from the spatial and temporal discretization of the governing system of PDEs. Consequently, the memory requirement to store the sparse Jacobian is significant. Such very large linear systems result in the linear solver component to constitute more than 70% of the computation time in reservoir simulators. Iterative linear solvers depend on effective preconditioners, which can be hard to parallelize to the extent required by many-core simulations. In a first step, we investigate feasibility of using the locally linearizing nonlinear multigrid method Full Approximation Scheme (FAS) in serial to establish algorithmic performance.