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

Change the heat exchange model "Ocean" for ice

YD
yafei Duan, modified 2 Years ago.

Change the heat exchange model "Ocean" for ice

Youngling Posts: 10 Join Date: 9/28/14 Recent Posts
Hi,
The ice temperature is -1.8 ?;
I want to change heat exchange model "Ocean"
Since this model get the only water temperature above zero degree. However, for salty water, Sea, the water temperaure could be lower than 0.
I want to change it in the heatu.f90
However
when I change the code
from
"
else
if (qtotk > 0.0_fp) then
sour(nm, k0, ltem) = sour(nm, k0, ltem) + qtotk/(thick(k0)*h0old)
elseif (r0(nm, k0, ltem) >0) then
sink(nm, k0, ltem) = sink(nm, k0, ltem) - qtotk/(thick(k0)*h0new*r0(nm, k0, ltem))
elseif (r0(nm, k0, ltem) > 0.0_fp .and. r0(nm, k0, ltem) < 0.01_fp) then
!
! No addition to sink when the water temperature is lower than 0.01 degree.
!
else
msgcount = msgcount + 1
endif
endif
"
to
"
else
if (qtotk > 0.0_fp) then
sour(nm, k0, ltem) = sour(nm, k0, ltem) + qtotk/(thick(k0)*h0old)
elseif (r0(nm, k0, ltem) >-1.8) then
sink(nm, k0, ltem) = sink(nm, k0, ltem) - qtotk/(thick(k0)*h0new*r0(nm, k0, ltem))
!elseif (r0(nm, k0, ltem) > 0.0_fp .and. r0(nm, k0, ltem) < 0.01_fp) then
!
! No addition to sink when the water temperature is lower than 0.01 degree.
!
else
msgcount = msgcount + 1
endif
endif
"

There is no different result. and the water still could not decrease any more after dwon to 0.

Why?

Many thanks.
VM
Vahid Mohammadian, modified 1 Year ago.

RE: Change the heat exchange model "Ocean" for ice

Youngling Posts: 5 Join Date: 1/29/16 Recent Posts
Hi yafei

I faced a similar issue,

The temperature is about -2 in the domain that I am gonna simulate (because of highly saline sea water). But the model does not accept negative values for water temperature.

Did you find any solution?

Thanks
Vahid