1. Simulation input, run, performance

RE: 3D D-Flow FM model coupled with D-Waves ?

2021-08-23T15:40:39Z

Hi Dmitry,

It is correct that it is not possible to couple D-Flow FM 3D model with D-Waves. Unfortunately, it is unclear when it will be made available. In 2D, it is possible to couple a D-Flow FM model with D-Waves. You can find all the information in the User Manual (http://content.oss.deltares.nl/delft3d/manuals/D-Flow_FM_User_Manual.pdf) in chapter 17 "Coupling with D-WAVES".

Best regards,
Tjitske

3D D-Flow FM model coupled with D-Waves ?

2021-08-20T07:25:17Z

Hi all,

D-Flow FM User manual states that there is no coupling available for 3D FM model with D-Waves yet.
Is it available as a beta functionality? (alpha will be OK too).

Regards,
Dmytro

sediment/morphology + tracers + sourcesink

2021-06-17T07:56:41Z

Hi all,
I'm doing some numerical exercises with sediments and passive tracers. First I ran simulations with sediments only and had no problems. When I added the tracers in a source discharge point (passive tracer, neither salinity nor temperature), error messages appeared and the model did not run (diag message: “the model needs “n” columns and you are providing “n+1” columns). When I manually removed one column directly from the “source-sink.tim” file, the model has shown the same message and re-generated that “such extra column”. I noticed that by removing the “bedload” fraction from the model setup, the problem has gone. So I think it's a bug: bedload and tracers. Has anyone experienced a similar issue? I have not run simulation with bedload and temperature/salinity altogether. So, I can not say if the issue would persist in a simulation (bedload + temperature/salinity) or if it is related to the source-sink feature.

Best regards,
Hemerson

(DelftFM 2021.04 Version: 1.2.110.68456M)

Nonconservative salinity, temperature in DFM

2021-05-11T20:29:32Z

I am developing a DFM domain for a small coastal lagoon, and encountering issues with salinity and temperature. Notably, I get salinity values which are greater than the initial condition or any boundary conditions. It seems to be a substantial mass, not limited to very shallow cells at the wetting/drying front. My initial condition is 0ppt, with a 34 ppt ocean boundary, and I see salinity values of 40 or more. In a 12 hour run, I summed the mass of salt exceeding 34ppt (i.e. sum( FlowElem_bac*waterdepth*max(sa1-34,0)) over the full domain). In some cases this showed as much as 1500 kg of "extra" salt (for reference, the volume of the domain is about 4e5 m3). Temperature also has output values that are sometimes greater than the initial condition or boundary conditions (with no additional heat fluxes).

Things I have tried to eliminate or constrain the source of this:
1. make initial condition and boundary conditions for salt and temperature constant and equal. In this case no new maxima are generated.
2. Compare running in 2D (Kmx=0) and 3D, one layer (Kmx=1): Different results but qualitatively the same behavior. Ultimately we want to run this in 3D, but currently just making sure it runs okay in 2D.
3. No change with increasing Epshu, changing limtypsa, CFLmax, changing TransportTimestepping
4. The domain is in UTM, and high resolution (dx ~ 3-10m). Tried translating the domain to be centered around 0 in case floating point truncation was a problem, but no change.
5. Mostly testing with 2021.03, but also tested with 1.6.1 with no change in behavior. All runs on Linux, using MPI (though it doesn't appear that the cells with excess salt are in the vicinity of inter-processor ghost cells).

The one setting that seems to have a big effect is BedLevType. The worst behavior seems to be with BedLevType=4. Type 3 is slightly better, and type 6 is the best (excess salt mass down to a max of 60kg).

Are there other settings in the MDU we should consider (aside from using Salimax and Tempmax)?

Thanks!
Rusty

RE: D-Flow tracer time series inside grid

2020-10-09T20:55:34Z

Hi Dmitry,
Thank you very much for your help.
You understood me very well. That's what I want! Apologies for my rough English.
I had already tried using a restart file (to get the hydrodynamic conditions at the "spill" instant) and at the same time space-varying initial conditions for tracer concentration. However, it seems that restart file "prevails" over space-varying initial conditions and I don't get tracer cloud at the specific part of the domain.
I will try editing tracer concentrations in restart file (using a netcdf4 Python module) prior launching the simulation.
By the way, do you know how many tracers are supported? I used both GUI: the 1.2.2.36603 version as well as the Beta Suite 1.3.2.37884.
Any other suggestion is welcome.
Thanks for your time and for the contribution to the forum.
Regards,
Sebastian

RE: D-Flow tracer time series inside grid

2020-10-09T07:27:24Z

Hi Sebastian,

If I understand you correctly, you want to put instantaneous sources of tracer to the model area?

This is best simulated via space-varying initial conditions or edited restart files.

In the (older) Delft3D-FLOW model a few matlab functions were available for this, which made possible to create space-varying initial conditions and restart files, namely delft3d_io_ini and trirst.

You can create your own initial conditions, which contain i.e. a tracer cloud at the specific part of the domain, using Delta Shell spatial editor (see appendices D and H of D-Flow FM User Manual). Then you just have to start your simulation from those initial conditions and investigate the behavior of the tracer.

Hope that helps,

Regards,

Dmitry

D-Flow tracer time series inside grid

2020-10-07T17:05:04Z

It is possible in D-FLOW to set up a tracer time series discharge within the domain?
I want to simulate and instantaneous release of a tracer (like a plume, a patch/cloud of pulluted water that cover a number cells). This is the case when a sensor located in an estuary registers an abrupt increase/decrease of water quality parameter. The idea is to forecast the transport of the plume as it moves downstream, with a 2D simulation.One way is using many "point source" and assinging a low salinity value (at this point is a fresh water estuary and dont want to affect density). However I also want to "realease" several peaks of different tracer at certain interval time in order to know how the transport time (from this point to another place) is changing with de meteorological condition.
Hope to be clear. Thanks.
Sebastian.

RE: Fourier Analysis

2020-02-20T19:25:15Z

Hi Nill,I believe the easy way to make .fou work is
1. to look up the mdu file, and
2. check the Tunit, DtUser, and Tstop settings.
3. you can use the number for the start and stop time in .fou file.
4. if the model works and produces fou output, then you can try to change the start time (in case you need spin up time in the model)

The default time step is 5 minutes, so if the Tunit is M; it should be the multiplication of 5 or if it is S; it should be the multiplication of 300. For template, I usually copy the number to excel and copy it back to .txt file so it looks better and ensure the tabs/spaces.Please let me know if this tips is work or not.Cheers! Good luck!

RE: Fourier Analysis

2020-02-05T08:06:26Z

Thank you Munawir. That works for me after changing the start time to 1200 instead of 2880 (because 2880 is not a multiple of user defined time step).Is there a reason why that my file do not work? Formatting aside, the difference is only about having the statistic start time to be slightly after the model start time (original statistic start time is 0).Or does the number of space between input matter? That does not seem to be included in the manual though.

RE: Fourier Analysis

2020-01-28T14:18:42Z

Hi Nill,
Please try this file for your model.
It works for me, it is for curr velocity magnitude.
Cheers,

Fourier Analysis

2020-01-21T07:05:49Z

I am testing out the Fourier analysis feature in a very simplified model. Yet it does not work out for me. Here are what I have tried: -
(1) Model run without temperature and Fourier analysis - model run normally
(2) Model run without temperature and with Fourier analysis for water level (wl) - model run normally, Fourier map results include statistics of water level
(3) Model run without temperature and with Fourier analysis for current velocity (ux, uy) - model run cannot start
(4) Model run without temperature and with Fourier analysis for current velocity (uv, this is the old keyword used in Delft3D Classic)- model run normally, but without map output of Fourier analysis for current
(5) Model run with temperature and without Fourier analysis (ct) - model run cannot start
(6) Model run with temperature with Fourier analysis period start in the middle of the model - model run stop progressing in the middle of the modelling period
From (2)-(5) I can deduce the keywords wl, ux, uy and ct are probably right as stated in the DFM manual. And from (3), (5) & (6) I can deduce the only reason with model not working is about the Fourier analysis for current and temperature.
The relevant setup are attached if anyone is interested to have a look. There is nothing confidential there and the test run take only 20s on my PC. I would love to hear from other users on the issue. My DFM version is 2018.04

WES3.m Input

2019-05-02T19:58:59Z

Hi,

I am trying to set up a tropical cyclone model in matlab and I want to use WES3.m to generate. But I can't find the wes_read_spw_input.m file. Where can I find it.

Thanks a lot,

Sven

RE: Use Integrated model in Delft3D FM suite 2019.01(1.5.1.41875)

2019-04-21T15:25:07Z

Another very simple question, Can we use irrigular grid in D-Wave?

thank's

Use Integrated model in Delft3D FM suite 2019.01(1.5.1.41875)

2019-04-17T13:08:21Z

Hello everyone, I am a new user of the model, I already realized modeling with model Delft 3D (coupled mode), but this time I try to make an integrated model, my question is the following :
is it possible to use an irrigular mesh in an integrated model, particularly in wave model?

any tutorial to explore integrated model in Delft3D FM suite ?

thank's

Example input files for integrating D-Water Quality Process in D-Flow FM

2019-04-05T22:40:55Z

Hello!

I'm trying to set up an integrated water quality and dfm simulation as opposed to a coupled DFM-DELWAQ run.

I've reviewed chapter 19 of the manual (Version 1.5.0) and still can't seem to figure out how to set up an integrated run.

Does anyone have any example files of an integrated water quality and dfm simulation?

Greatly appreciated!

Lily

NetCDF error code -40: Index exceeds dimension bound.

2019-02-09T00:15:02Z

Hello everyone!
I am a beginner and am trying to run my first model using Delft3D FM HMWQ (1.5.1.41875). It is a Wave-Flow-Mor model. After performing the configuration, I tried to run the model and found the following:

ERROR DelftTools.Shell.Core.Workflow.ActivityRunner.RunActivity - Finishing of Integrated Model has failed.
2019-02-03 16: 08: 05,169 [5] ERROR DelftTools.Shell.Core.Workflow.ModelBase.ChangeState - NetCDF error code -40: Index exceeds dimension bound. Consider cloning the file into NetCDF 4 format to enable data extdending

Has anyone ever had an error like this? Anyone have any clues on how to solve it?

Thanks.

Vinicius Pessanha

Exporting Quickin file (.dep-file) for a following simulation

2019-01-29T17:25:27Z

Hello Everyone,

I am working on a project for which I have to simulate three years of hydrodynamical and morphological development in an open river flow. Since some geometrical data/barriers (implemented via a *.bar-file) change every couple of weeks i am working with multiple simulations using the restart-file of the previous simulation. However, I checked the morphology development and it seems to me that the restart file does not contain any bathymetry information of the previous simulation.

With Quickplot I can export the 'bed level in water level points' at the end of the simulation and use the exported file as the bathymetry file for the next simulation. But since the bathymetry for the computation is defined negative in the upper direction I need to change the orientation of the bathymetry data. I can either do that using Quickin (operations -> combine depth and uniform value -> multiply depth and uniform value -> -1) OR I can directly export a negative Quickin file. However, when I directly export a negative Quickin file (-Quickin) the default value of -999 is also changed to +999. Thus, using the exported file as the bathymetry file for the next computation doesn't work because when I try to run the simulation +999 is not regarded as a default value but as bathymetry data, which is false and obvioulsy does not work with the simulation.

This doesn't happen when I manually change the orientation of the file. But it would be resonably faster only having to export the file and using it right away.

Is there any way around the manually changing method? Any way I can directly use the exported *.dep-file for the next simulation?

Thank you for your help.

Regards, Rebekka

Delft3dFM + D-RTC error

2018-11-14T03:04:41Z

Hi Everyone,

I am trying to set up an integrated delft3dFM+ RTC model. I would like to use dead band trigger for pump control. However, the current GUI does not seem to support this feature. So I modified the rtcToolsConfig.xml and add the following:

<trigger>
           <deadBand id="test">
               <conditionOn>
           <x1Series ref="EXPLICIT">input_ObservationPoint01_water_level</x1Series>
                   <relationalOperator>Greater</relationalOperator>
                   <x2Value>2</x2Value>
               </conditionOn>
               <conditionOff>
           <x1Series ref="EXPLICIT">input_ObservationPoint01_water_level</x1Series>
                   <relationalOperator>Less</relationalOperator>
                   <x2Value>1.8</x2Value>
               </conditionOff>
               <true>
               <trigger>
               <ruleReference>run1</ruleReference>
               </trigger>
               </true>

               <false>
               <trigger>
               <ruleReference>run2</ruleReference>
               </trigger>
               </false>
               <output>
               <status>output1</status>
               </output>
           </deadBand>
</trigger>

xmlspy did not complain but when I run the model it gave me the following error:

<line description="rtcToolsRuntime::rtcToolsRuntime(...) - error - trigger* schematisation::getTrigger(TriggerComplexType triggerType) - error - trigger is not implemented" level="1"/>

Can anyone please let me know what may have caused above error? Many Thanks.

Yafei

RE: Z-Model Set up

2018-11-06T11:15:40Z

Hi everyone,

I am also trying to implement several Z-layers in my model (in Delft3D Fm 2018.03). I have added the suggested part in the *.mdu file:

Kmx               = 5
Layertype         = 2
StretchType       = 1
stretchCoef       = 30 20 20 20 10

However, when i try to run the model in the Deltashell it doesn't seem like it's using the given layers (e.g. It doesn't show the velocity in a certain layer when looking at the results). I am thinking it could not be working because when starting the model in the Deltashell it says it is first exporting the *.mdu file - which would not contain any layers since I can't add than in the Deltashell. However, opening the *.mdu file after running the model still shows the layers defined like above. So it doesn't overwrite the .mdu file after all. How is this possible? Any ideas what might be wrong?

Thanks! Regards, Rebekka

RE: Access Dflow FM model results during model run

2018-11-06T10:16:10Z

Hi everyone,

i am still confused about how to actually look at the map results. I am working with the beta version Delft3D FM 2018.03 in Windows 7. As soon as i start the model run the Gui is blocked by the running process window. I can locate the *.his and *.map files in the output files. However - how can i look at them in the Gui?

Thank you for your help.

Regards, Rebekka