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

Boundary Numerical Instability

Fernando Barreto, modified 6 Years ago.

Boundary Numerical Instability

Padawan Posts: 74 Join Date: 6/21/12 Recent Posts
I'm doing a simulation of a meso scale region to understand the baroclinic circulation of this place. My model contains 32 Z-layers and 3 bounds (north, east, south). To the boundary configuration I chose to impose current in the top and east boundary, and elevation in the south boundary. At this configuration I got bizarre instabilities in the boundary. It was generated a velocity jet in the interface between boundary types (elevation and current)
In another configuration I imposed current in north bound., and in the half of the east boundary; the remaning sections was elevation type (south and the rest of east). In this configuration was also generated a jet in the interface between bound types (fig). The jet can be seen in the east boundary.
Someone knows how to solve this. Can this be related with the z-model? if not, any advice?
SA
Steven Ayres, modified 6 Years ago.

RE: Boundary Numerical Instability

Padawan Posts: 33 Join Date: 4/21/11 Recent Posts
Fernando,

Were you able to resolve the boundary problems in your model? I have something similar going on in my 18-level Z-model with salinity. The model is creating velocity jets in the middle of an open stage boundary. I assume this has something to do with the salinity constituent but I'm not sure how to go about fixing the problem.

Thanks,
Steve
Fernando Barreto, modified 6 Years ago.

RE: Boundary Numerical Instability (Answer)

Padawan Posts: 74 Join Date: 6/21/12 Recent Posts
Dear Ayres,
my problem was the input of different boundary conditions (velocity and wl) in adjacent boundaries. This was resolved when I used riemann boundary conditions. What's your model configuration? boundary type, dimension?

regards,
Fernando
SA
Steven Ayres, modified 6 Years ago.

RE: Boundary Numerical Instability

Padawan Posts: 33 Join Date: 4/21/11 Recent Posts
I have water level boundaries for the tide boundaries and an inflowing river total discharge boundary. I think the problem is short wave reflections in the water level boundaries. I raised the alpha coefficient to 1000 and this helped extend the simulation time to a few days, but I still seem to get re-circulation around the open boundary. I again raised the alpha coefficient to 10,000 and am testing now.

A formula is given in the manual to determine the alpha coefficient -

Equation 9.76, Water level boundary: alpha = Td * sqrt(H/g), where Td is the time it takes for a free surface wave to travel from the left boundary to the right boundary of the model area.

Since the model is around 30 miles (48 km) across the shortest open tide boundary, and assuming the tidal wave travels around 30 miles/hour or 13.4 m/s, Td would equal 3600 seconds. Using a typical depth of 200 meters at the open boundary, alpha would be on the order of 16,000. Does that sound right?
Fernando Barreto, modified 6 Years ago.

RE: Boundary Numerical Instability

Padawan Posts: 74 Join Date: 6/21/12 Recent Posts
Hi Steven,

When I use WL boundaries I usually put alpha as 1000. I didn't know how to calculate alpha, it's really to know. In Riemman we didn't use alpha, Riemann is a weaky reflective boundary condition, so I don't know a lot about alpha.
There are many possibilities to cause your error. Could you attach your model ?, I think it's easier to help.

Regards,

Fernando Barreto
SA
Steven Ayres, modified 6 Years ago.

RE: Boundary Numerical Instability

Padawan Posts: 33 Join Date: 4/21/11 Recent Posts
The problem was definitely related to boundary reflections. When I set alpha to 10000, the boundary velocity jets disappear.