intro story DELWAQ


DELWAQ is the engine of the D-Water Quality and D-Ecology programmes of the Delft3D suite. It is based on a rich library from which relevant substances and processes can be selected to quickly put water and sediment quality models together.

The processes library covers many aspects of water quality and ecology, from basic tracers, dissolved oxygen, nutrients, organic matter, inorganic suspended matter, heavy metals, bacteria and organic micro-pollutants, to complex algae and macrophyte dynamics. High performance solvers enable the simulation of long periods, often required to capture the full cycles of the processes being modelled.

The finite volume approach underlying DELWAQ allows it to be coupled to both the structured grid hydrodynamics of the current Delft3D-FLOW engine and the upcoming D-Flow Flexible Mesh engine (1D-2D-3D) of the Delft3D Flexible Mesh Suite (or even other models such as TELEMAC).

'DELWAQ in open source' is our invitation to all leading experts to collaborate in further development and research in the field of water quality, ecology and morphology using Delft3D. Feel free to post your DELWAQ related questions or comments in this dedicated forum space. If you are new to DELWAQ, the tutorial (in the user manual) is a good place to start. A list of DELWAQ related publications is available here.




Sub groups
D-Flow Flexible Mesh

Cohesive sediments & muddy systems


Message Boards

How does DELWAQ deal with sediment concentrations in cells that wet & dry?

Ben Williams, modified 7 Years ago.

How does DELWAQ deal with sediment concentrations in cells that wet & dry?

Jedi Knight Posts: 114 Join Date: 3/23/11 Recent Posts

I have what may seem like a silly question (below). I have had a skim through the WAQ manual and process technical reference manual and I haven't found anything that jumps out at me in terms of answering it. Maybe someone on this forum might be able to provide some assistance?

DELWAQ deals with fine suspended sediment transport as a mass balance process. For a simulation dealing with the tidal zone (such as an estuary), there will be times when some cells are either dry or will be in a very shallow depth. Therefore, where the water depth becomes shallow within a grid cell, the concentration (for a given mass) must necessarily increase.

How does DELWAQ deal with calculating IM1, IM2, IM3 (all specified in the results files as gDM/m3) in these circumstances to avoid 'divide by zero' issues or circumstances where the water depth is so shallow that it would lead to artificially high IM concentrations in the water column?

Best regards,

Arjen Markus, modified 7 Years ago.

RE: How does DELWAQ deal with sediment concentrations in cells that wet & d (Answer)

Jedi Knight Posts: 223 Join Date: 1/26/11 Recent Posts
It is not that silly a question, but there are several parts to it:
- Treatment of drying and flooding in general
- Treatment of processes, notably sedimentation, in very shallow segments

To answer the first one:
If a segment dries up, the flow rates at all sides are supposed to become zero. This is an indication for DELWAQ to skip the evaluation of all transport processes. Only when the flow rates become non-zero again. the transport over such sides is calculated again. All the while the volume of the segment remains constant.

This behaviour is part of the "contract" with any hydrodynamic model that DELWAQ can be coupled to. Sometimes that is a trifle difficult to realise - for instance with some of the drying/flooding optiions in Delft3D-FLOW. There are more technical details involved and we introduced a more sophisticated algorithm to detect the drying and flooding some time ago, but this is short story.

To answer the second one:
Segments that are about to become dry are very shallow and that means that there are some numerical problems to be solved. Sediment can sink to the bottom, but within one time step
there should not be more mass settling than is actually present. This could happen because the depth may be so small that within a single time step the sediment could be deposited, that is:
depth < settling velocity * time step.

The routine that takes care of settling uses a limiting condition to avoid this problem - at most the mass that is present can settle. Also if the depth is lower than a given threshold, no settling is assumed to occur at all.

Effectively, if a segment is "dry" (depth below a given threshold), it is assumed no processes are taking place. This is not always the best solution - like with Ecoli bacteria, but it avoids unrealistic concentrations.