Due to the large spatial scale of sand waves and the different time scales involved for the tidal flow (hours) and the bed pattern changes (years to decades), modelling the sand wave behaviour from its initial stage up to its final equilibrium shape is a challenge, let alone to include the spatial variation within a sand wave field.
In this paper a numerical solver is presented that describes sand waves from their initial state up to their final equilibrium. It is a 2DV idealized model, based on non-linear stability analysis. Both the model equations and the numerical setup are described and the model is validated against linear stability analysis.
To investigate variations in a field of sand waves and the interaction between individual waves, simulations on large domains are presented and show that this solver is able to describe the development of a realistic sand wave field, with variations, from an initially flat sand bed with small random disturbances.
The solver shows to be a promising tool, it is computational fast due to the efficient numerical algorithms and is easy to extend with other physical processes. Results show good agreement with analytical results in the linear regime and with field data.