Migration status

  • Home Page: Activity map.

intro story DELWAQ

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.

** PLEASE TAG YOUR POST! **

 

 

Sub groups
D-Flow Flexible Mesh
DELWAQ

Cohesive sediments & muddy systems

 


Message Boards

Tidal water level variation inside an estuary

SS
Saeed Shaeri, modified 5 Years ago.

Tidal water level variation inside an estuary

Youngling Posts: 6 Join Date: 10/7/12 Recent Posts
Dear forum members,

My question is about tidal water level variation inside an estuary.

To make my question simple, at a location some kilometres upstream of a tidal inlet, water level is needed. The flow domain is forced by a boundary condition which is a variable water level, some kilometres offshore of the inlet.

My understanding is that the changes of the boundary, make the water level variation all across the modelling domain. Therefore, it should take some time steps which the very far cells from the variable boundary are influenced by such changes. I believe that is the case for a wave modelling; however, in a flow model I experienced different.

Hence, my question/confusion is that if I should observe water level variation at all cells exactly from the initial time step. This would also show no time-lag between any two cells of the domain.

Looking forward to your kind comments,
Regards,
Saeed
PS
Phil Shepperd, modified 5 Years ago.

RE: Tidal water level variation inside an estuary

Padawan Posts: 33 Join Date: 5/31/11 Recent Posts
If you think of your model not as a computational one, but as a physical one, made of fibreglass, filled with water ,sitting in a laboratory. You can visualise that it will take some time for a water level change at the boundary to propagate up through the estuary to your point of interest.
Your FLOW model works in essentially the same way - a water change at the boundary will move through the model grids cells, one timestep at a time, until it gets to the point of interest. So if you were to look at some trim* output, you would see the tidal wave propagating just like it would in the physical model, or the real world.
PS
Phil Shepperd, modified 5 Years ago.

RE: Tidal water level variation inside an estuary

Padawan Posts: 33 Join Date: 5/31/11 Recent Posts
Possibly.

As a quick check, the speed of a wave (tidal or otherwise) in shallow water is sqrt (gd) where g is the acceleration due to gravity (9.81 m/s) and d is the local water depth (m).
Giordano Lipari, modified 5 Years ago.

RE: Tidal water level variation inside an estuary

Youngling Posts: 12 Join Date: 3/23/11 Recent Posts
Hi Saeed.

The intuition that there's a propagating long wave into the estuary is correct. For this to be an effective start of a numerical simulation, you need to consider what your initial water level is and the phase of the tide at simulated time zero.

You can keep in mind the common distinction in numerical modelling between a 'cold start' and a 'hot start'.

In the cold-start case, you begin from fictitious conditions such as a quiescent basin that, then, is acted upon by the relevant forcing (tide, discharge wind and all that applies). A cold state hints at no motion here. The motion in the water body, then, develops until reaching a dynamical equilibrium between the applied forcing and the resisting forces (bottom and internal friction). You recognize such equilibrium because it's either a steady state or an unsteady process that has lost memory of the fictitious initial conditions (e.g. a proper tide in the estuary).

The hot start is when you start, actually restart, a simulation from a snapshot situation taken from a previous (correct) simulation. So, at simulated time zero, velocity and water levels are variable all over the domain, but the balance of forces already holds in each cell of the domain.

By contrast, at time zero in the cold start, you are unlikely to have balance between the initial condition in the inner domain and the boundary condition. That's I guess what was puzzling you. Well, in that case you run the simulation until you observe regular tidal oscillations everywhere in the basin. This is the so-called spin-up time. As long as your domain's description is realistic, the simulation will take the right course (see caveats below though). To judge what a regular tidal oscillation is, it's better to refer to measurements but common sense will certainly guide you ruling out oddities. There's some more physics specific to estuarine circulation to take into account here but, as a first pass, it'd be alright.

The final touch is a good match of initial condition and boundary conditions. At the boundary you have the tide which is principally synchronised with the calendar. Start then a simulation at a calendar time when the predicted water level at the boundary is not largely different from the water level set in your initial condition (or the other way round). It seems preferable to me that you also start at flood so that the water comes into the estuary and fills it (but I may be overlooking other aspects now, so happy to be integrated).
Either way, there's one more feature in the numerical parameters that has been devised to keep issues at bay. You can ask the solver not to apply the boundary condition at full value from the outset, but to increase the forcing gradually within a limited period (that you choose) so as to avoid a jump start (which, numerically speaking, can be obnoxious in cases).

I don't know how large your estuary and tidal range are, which determines the order of magnitude of it all. As a guess out of my hips, I would say that with a ramp-up time of 1-2 day and discarding the first week of simulated time you should be sorted or getting there. With the cold start, you always need to take into account a buffer time of sorts though.

Happy to be corrected/integrated. The no-slip condition applies to water, not to water experts. Hope this helps.

Giordano
www.watermotion.eu