2.2.1 Waves and Coastal simulations

Domain Dries Up Due to Overflow

2023-06-27T10:00:48Z

Hi,

I am simulating a non-hydrostatic model with a simple, low-angle-slope bathymetry using 100 z layers and 0.01 mins timestep. I also impose unidirectional 4 cm/s current as boundary condition and uniform-and-constant zero water level boundary on the other end. However, the model dries out (ie. the water elevation becomes equal to the water depth) after 7 mins (a half of smoothing time, 14 mins). The velocity profile shows extremely high value, up to 1000 m/s.
Has somebody encountered this situation before? Any lead is much appreciated.
Please find attached the model set up files and the velocity profiles.

Thanks,
Rushdi

Delft3D stack overflow error

2023-06-26T18:56:18Z

Hi,

I'm new to Delft3D and am running the Wave module and getting a stack overflow error:

.
.
.

+SWAN is processing output request 1
1 file(s) copied.

C:\D3D\Southern California>rem exit

>>...End of SWAN run
*** ERROR: file 'norm_end' expected to signal a correct SWAN calculation

forrtl: severe (170): Program Exception - stack overflow
Image PC Routine Line Source
swan_4072ABCDE_de 00000001401D7D87 Unknown Unknown Unknown
swan_4072ABCDE_de 00000001400D00AE SWOEXF 4403 swanout1.for
swan_4072ABCDE_de 00000001400C068C SWOUTP 383 swanout1.for
swan_4072ABCDE_de 000000014002CE1A SWMAIN 598 swanmain.for
swan_4072ABCDE_de 000000014002B764 MAIN__ 110 swanmain.for
swan_4072ABCDE_de 00000001401D4FBC Unknown Unknown Unknown
swan_4072ABCDE_de 00000001401B2AE5 Unknown Unknown Unknown
KERNEL32.DLL 00007FFD0AC47614 Unknown Unknown Unknown
ntdll.dll 00007FFD0B7226F1 Unknown Unknown Unknown

I'm running on:

Windows 10 Enterprise
Delft3D 4.03.01
Delft3D-WAVE Version 3.05.04.7545M, Sep 13 2017, 20:57:19

It appears the model runs through the wave computations okay, but crashes when it's trying to do the output:

use of rescaling in 100.00 % of wet points with maximum in spectral space = 247
accuracy OK in 99.86 % of wet grid points ( 98.00 % required)

++ trace subr: HSOBND
++ trace subr: SWOUTP
+SWAN is processing output request 1
++ trace subr: SWORDC
++ trace subr: SWODDC
++ trace subr: SWOEXC
++ trace subr: EVALF
++ trace subr: EVALF
++ trace subr: SWOBSTO
++ trace subr: SWOEXD
++ trace subr: SWIPOL
++ trace subr: SWIPOL
++ trace subr: SWOEXA
++ trace subr: SWOINA
++ trace subr: ANGDEG
++ trace subr: SWOINA
++ trace subr: ANGDEG
1 file(s) copied.
1 file(s) copied.

My domain is pretty large with:

Gridresolution : MXC 465 MYC 1126
: MCGRD 407901

I watched the system memory usage as it runs and it's < 15%. CPU usage is 100% during the iterations, but it drops down to ~5% when it crashes processing the output request.

I derefined the grid to:

Gridresolution : MXC 233 MYC 563
: MCGRD 102115

and it runs fine.

Any help would be greatly appreciated.

RE: Bed stress coefficient ERROR

2023-06-14T12:35:17Z

I am making a spatially varying roughness file and everytime I generate this in QUCIKIN it produces negative values even though I have specified the values I wanted. Have you had this issue?

Errfile Terminating error: Input file missing

2023-06-08T13:45:29Z

The model runs for around half of time steps then shows

file 'norm_end' expected to signal a correct SWAN calculation

The process cannot access the file because it is being used by another process.

Problem in Sediment Transport with Discharge Input Only in Flow

2023-05-17T12:26:21Z

Dear Delft3D Community,

I am currently working on Delft3D-FLOW, utilizing a 3D model to observe sediment transport under different discharge conditions. While the simulation itself concludes without any errors or warnings, the results are not in line with expectations. Despite the velocity distribution appearing reasonable across the intake and discharge locations, some suspended sediment transport results in too far areas which are illogical, and cumulative erosion/deposition distribution is not calculated, seeming like an instability problem??

In FLOW, my model has 10 layers vertically, where the layer thickness decreases from top to bottom (bed layer). I have specified each intake section in 15 grid cells (5 horizontal*3 vertical) and the discharge section in 35 grid cells (5 horizontal*7 vertical). These discharge grids are specified for the 2nd, 3rd, and 4th layers from the top. For each grid of intakes, the flow is determined as -5.55 m3/s with no additional sediment fraction and for each grid of the discharge, the flow is determined as 9.52 m3/s with no additional sediment fraction. Only sand (as a non-cohesive sediment) has been considered through the simulation process, where the simulation duration is 12 hours and the time step is 0.1 min. The only input is the negative (pumping) discharge and the positive discharge where each of the discharge types is defined as “Normal” and the interpolation method is chosen as “Linear”.

The simulation results are expectable in terms of water level and velocity but there is suspended sediment transport in too far areas from the input (discharge) locations where expected transport does not appear near the input area. The cumulative erosion/sedimentation results are also all zero.

Do you have any experience or ideas related to such a problem?

Consequently, I suspect that the issue may stem from the specified discharges; however, I am unable to comprehend why the sediment transport fails to correlate with the observed velocity behavior. Could you provide insights into the possible reasons behind this problem?

With best regards,

Hydraulic structures in Delft3d

2023-05-15T05:24:39Z

Good Day Delft community!

Our group is currently conducting a design project through delft3d. The objective of the project is to reduce storm surge induced coastal flooding in a specific area using hydraulic structures (Seawall, Breakwater and Dyke).

We have already done the initial water level simulation in the area and we would like to know if it is possible to insert hydraulic structures using delft3d. We plan on simulating and recording the water levels during a typhoon's occurrence with and without a hydraulic structure in order to quantify its effectivity. Is this method possible using delft? If yes, is there a specific procedure that we can follow in order to do it?

Thanks in advance!

RE: Simple Wave Run Error

2023-03-29T01:50:02Z

Hi! The model work just fine when I remove the wave diffraction and setup.

Thank you very much!

RE: Simple Wave Run Error

2023-03-24T12:23:29Z

Hi Sofia,
I think this warning is quite general.
I have a quick look at the MDW files. I would suggest running the model without wave setup and wave diffraction. Would it be run normally?
Cheers, Nawir

Simple Wave Run Error

2023-03-24T00:34:10Z

Hi!

I am having trouble running my .mdw file. I am suspecting the bathymetry hinders it to run normally. Can I ask for help what are the possible reasons to this? The run also generates a message like this:

*** ERROR: file 'norm_end' expected to signal a correct SWAN calculation
ERROR: *** ERROR: file 'norm_end' expected to signal a correct SWAN calculation
Wave exited abnormally

forrtl: severe (172): Program Exception - exception code = 0x7363 (29539)
Image PC Routine Line Source
KERNELBASE.dll 00007FF9B4E78B9C Unknown Unknown Unknown
VCRUNTIME140.dll 00007FF9B0E166C0 Unknown Unknown Unknown
wave.exe 00007FF725813701 Unknown Unknown Unknown
wave.exe 00007FF725600EDA Unknown Unknown Unknown
wave.exe 00007FF7256ACB5D Unknown Unknown Unknown
wave.exe 00007FF7255B9552 Unknown Unknown Unknown
wave.exe 00007FF725551615 Unknown Unknown Unknown
wave.exe 00007FF7257F8552 Unknown Unknown Unknown
wave.exe 00007FF7258203D4 Unknown Unknown Unknown
KERNEL32.DLL 00007FF9B5AC13D2 Unknown Unknown Unknown
ntdll.dll 00007FF9B7B45454 Unknown Unknown Unknown

Thank you!

Calculating vorticity from FM output

2023-03-22T21:09:30Z

Hello,

I was wondering if anyone has been able to calculate vorticity from the map outputs of a delft3D FM model. I can figure it out for a regular grid, however am struggling to do so for a triangular, irregular grid.

Thanks

Unable to deactivate refraction in WAVE module

2023-02-03T15:27:16Z

Dear Delft3D community,

in my installation of Delft3D it seems to be impossible to deactivate the refraction process in the WAVE module / SWAN.
Comparing WAVE runs (tried both coupled and standalone) showed that I get exactly the same results regardless of the value (true / false) of the refraction keyword in the mdw file. The SWAN INPUT file does contain the correct keyword OFF REFRAC and no error / warning message is given in SWAN or D3D output.
I observed the same behaviour when trying to disable frequency shift. For wind growth it worked, however.
The same problem occurs on both my local Windows installation (Delft3D-WAVE Version 3.07.00.62434) and the Linux server (Delft3D-WAVE Version 3.07.00.64862) I'm using.
My colleagues (who use the same program version on their computers) encounter the same behaviour when they try to disable refraction in their models.

It seems to me that the problem is not directly related to the WAVE version or the model setup. Nevertheless, I have attached a simple model setup that yields the same simulation results regardless of the refraction keyword.

Any suggestion would be much appreciated.

Cheers,
Arne

Combing Sea Surface Height (SSH) with Tidal forcing

2023-01-18T00:02:29Z

I am currently trying to add SSH to an astronomic tidal signal in Delft3D4.
I have a single offshore boundary and within that, I have included the two water level forcing signals. One is driven by Astronomic constituents. The other is a time series gathered from CMEMS.

When I run the simulation, however, the results do not seem to be the combination of these two signals. The tidal signal does not seem represented at all (See below).

My assumption was that the two signals would be added together and the resulting signal would be forced at the boundary.

Any ideas on why this might not be working?

I've included the PLI and BC files below.

FLOW ERROR：NEFIS message 6030: On writing variable pointer table, table 0

2022-12-28T14:29:08Z

Hi,

In the FLOW simulation, when I try a short time case, the model is working fine. But when I try small step and long time simulation (large amount of data), the model is always shutting down normally without finishing the run, see figure 1. I checked the log and it shows "ERROR: NEFIS message 6030: On writing variable pointer table, table 0", as in Figure 2.

What is the cause of this? Is it related to the settings of my computer system or the hardware of the machine itself? Or is it due to boundary conditions or mesh quality? My model has a total of 7000 meshes.

With best regards
Zhen Wei

RE: Bed stress coefficient ERROR

2022-10-31T01:47:14Z

UPDATE: I was able to solve the problem by checking the rgh file for zeros or negative values generated after exporting before copy and pasting the values inside the file

The xyz file seems not to be plotted in accordance with our grid area

2022-10-13T09:25:33Z

Good Day! We would like to ask a question with regard to dealing with the software Delft3D Flow. I hope you we're able to help us and make our study more explicit in its scope. So we have this reference that we got from youtube, and it's very useful for us though it does not have audio. And currently, we're okay with simulating not until we encounter this trouble in using the QUICKIN part in the Delft3D. Well, the xyz file seems not to be plotted in accordance with our grid area. As in the reference video, the xyz file plot very accurately in its grid area. May we know what might error we did in having this issue? Kindly go to 12:07, as we're in that part. Reference Video Link: https://www.youtube.com/watch?v=ZLOmFUC9IbA Thank you so much for the invaluable help and looking forward to communicating with you for better understanding.

Weird crash after few days of simulation

2022-10-07T09:20:04Z

Hi all

I'm trying to implement a simulation of the Strait of Gibraltar (South of Spain), with a domain spanning from the westernmost limit of the Iberian peninsula (San Vicente Cape) to approximately half of the southern Spain Coast length (Cartagena). It is few hundreds of km, hence a quite big domain.
The grid is curvilinear, with spherical coordinates, and resolution is of order of 1-2 km in the large areas and few hundreds of meters in the very Strait.
I interpolated bathymetry from a variety of sources and then I smoothed it by iteratively computing the slope factor, so to guarantee it is never higher than abs(0.2).
The vertical coordinates are sigma, with 30 uneven layers.
My intention is to improve horizontal and vertical resolution, and run it as a decomposition domain, but first I want to verify that it works adequately with only one domain.
The model is forced with astronomic water level conditions at both open boundaries (the western and the eastern one). Yes, I know that the same boundary condition is not recommended, but the domain should be big enough to not be an issue. For both boundaries I defined two conditions, the southern and the northern extremes, with tidal constants of 8 constituents, interpolated by the TPXO9-atlas model. M2 amplitude is of order 1 m and 2 cm at western and eastern boundary, respectively, while the other constituents (only semidiurnal and diurnal) are somehow smaller.
The initial condition is interpolated from the MedSea Reanalysis Model from Copernicus.
In a first instance I only configured sea level and current, nor salinity neither temperature, just to make the most simple approximation.

The run goes very well during the first 6 days, and then suddenly crashes. The more struggling fact is that the problem originates only at the northernmost cell of the western boundary and after 6 days of smooth and expected behavior. At that cell, despite the water level matches very well with the one imposed by the astronomic forcing, the current velocity behaves quite weird. It is not harmonic and, after 6 days it suddenly (with no previous changes of the sea level) starts increasing up to completely unreal values (of order of tens of m/s). Obviously, after few iteration, the CFL condition starts failing, and the model ends crashing.

I tried different configurations of viscosity, roughness and 3D turbulence scheme, but nothing works.
Any idea of what is wrong?

I attach the .bca, and .mdf files and few screenshots of sea level and current of the western boundary. The blue line is the northern cell and the red one is the southern one.

Thank you in advance.
Best regards
Simone

Bed stress coefficient ERROR

2022-09-23T10:22:40Z

Hi everyone. I have encountered an error about bed stress. The simulation stopped at part five. It says "ERROR Zero or negative value found for Bed stress coefficient" in the diagnosis file. How can I fix this? This happened after I created a space-varying roughness file. I have attached the diagnosis file, MDF and roughness file for your reference

Sincerely,
James Paul

RE: missing libnetcdf.so.7

2022-09-09T13:48:26Z

I sorted it out myself, just install the netcdf-c-4.3.3 version and add the lib into LD_LIBRARY_PATH environment.

missing libnetcdf.so.7

2022-09-09T11:35:44Z

hi,

I am using the 65936 version and all examples run normally except the one in dflowfm/e100..._FM_WAVE/, below is the error message. looks like this library is not available in the netcdf-c-4.6.1 I installed, which version of netcdf I should use?

RE: Inputting .xyz file for samples not working

2022-08-30T05:22:55Z

Hi Charles!

Can you share your *.xyz file to check it and try to find the source problem?

Regards

Gabriel