Home_intro

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

Wave boundary condition - point source

SM
Samantha Maticka, modified 1 Month ago.

Wave boundary condition - point source

Youngling Posts: 6 Join Date: 10/2/20 Recent Posts
Hello,

I am trying to run Delft3D with just waves. I am using a curvilinear grid and would like to force it with time-varying waves that approach from the southwest. Can anyone advise on the best wave to do this? What I have tried is resulting in a point-source effect of waves rather than a spatially uniform forcing. I included pictures below.

I have field observations from a single location that I would like to force my model with. It is about 20km from my area of interest. Currently, I am using a curvilinear circular domain, and I am applying a non-stationary time-varying boundary forcing by using a TPAR file, specifying a 'segment' type of boundary in the [Boundary] block of the .mdw file, and declaring specific times to perform the calculations as [TimePoint] blocks. Something like this:

[TimePoint]
   Time                 =  0.000000000000000000e+000             
   WaterLevel           =  1.000000000000000000e+000         
   XVeloc               =  0.000000000000000000e+000         
   YVeloc               =  0.000000000000000000e+000         
   WindSpeed            =  0.000000000000000000e+000         
   WindDir              =  0.000000000000000000e+000  
[TimePoint]
   Time                 =  60.0      
   WaterLevel           =  1.000000000000000000e+000         
   XVeloc               =  0.000000000000000000e+000         
   YVeloc               =  0.000000000000000000e+000         
   WindSpeed            =  0.000000000000000000e+000         
   WindDir              =  0.000000000000000000e+000  
...
...[Boundary]
   Name                 = PointP6                  
   Definition           = xy-coordinates
   StartCoordX          = 575571.0
   EndCoordX            = 575571.0
   StartCoordY          = 8177208.0
   EndCoordY            = 8177208.0 
   SpectrumSpec         = from file        
   Spectrum             = spectral.bnd


I tried using a 15-point segment for the boundary (~1/15th of the perimeter), but the accuracy reported in the PRINT file from SWAN suggested that the model did not converge. So I reduced it from 15 points along the encompassing boundary to 2 points. It looked to converge, but both results had somewhat of a 'beam-like' behavior.

Any suggestions would be greatly appreciated!!

Thank you emoticon
Sam
This is what our results look like. For a depth like this (outside the red circle is not part of the domain):
The significant wave height looks like this (1st is 15-point boundary, 2nd is 2-pt):

and 

but this is what we would like them to look like: