[daineythomas]

Dear all,
I have been trying to model the scouring near seawall toe for 15 hours but unfortunately, some error occurs. After 500 seconds accretion takes place faster than erosion.
Could anyone please help me to resolve the issue?
I do not have any idea about incorporating seabed slope. I would be grateful if anyone could help.
Input files and output is attached.

Thanks
Dainey

Control.txt

C 11 6
C 12 3
C 13 3
C 14 7
C 15 21
C 16 3

B 1 0.05

B 10 0.0 30 0.0 0.05 -0.3 0.6

S 10 29 30 0.0 0.05 -0.3 0.6

M 10 4

M 20 2

ctrl.txt

B 10 1
B 50 0.0001

B 90 1
B 92 4
B 93 0.28 2.68
b 96 6 12
B 98 3
B 99 0

D 10 4
D 20 2
D 30 1

F 30 3
F 40 3
F 42 0.7
F 60 0.3

I 12 1

N 40 3
N 46 225000
N 47 0.2

M 10 4

P 10 1
P 27 1
P 30 20.0

P 52 0.5
P 52 1.0
P 53 1
P 55 20.0

P 77 1
P 78 1
P 79 1

P 122 1

P 123 -0.3
P 123 -0.2
P 123 -0.1
P 123 0

P 126 1

S 10 1
S 11 1
S 12 0
S 13 0.18
S 14 0.2

S 16 1
S 19 360000.0
S 20 0.000024
S 22 2700
S 30 0.047
S 32 4
S 37 2
S 41 3
S 42 3
S 47 1
S 48 0.005

S 57 0.0
S 71 25
S 72 30
S 80 2
S 90 1

T 10 2
T 36 2

W 22 -9.81



seabed profile after 500 s.jpg

14 hours.jpg

MAXI EROSION.jpg

[pjoynt]

Hi Dainey,

I tried running your simulation using the ctrl and control files you provided. Potentially the waves aren't behaving as you are expecting? After 2.5s the wave profile looks like this in paraview (attached).

I would first ensure you are getting your waves correct before jumping into sediment transport.

Cheers.

[kamath]

Pjoynt is right (Thank you!)

Dainey, the tank is too low.

You have:
zmax=0.6,
d=0.30, H=0.28

So about 0.58/0.6 (not including 0.3 m bed) is filled with water, with 0.16 m for air. This in you case is 3 cells.
This is insufficient.

I suggest starting with zmax=1.2 m atleast.

[daineythomas]

Dear Dr Kamath and Pjoynt,
Thank you so much or helping

[raul_xusernamex_]

Hi how are you?.
You know with your original input files "ctl" and "control" the iteration is stopped with the message "Stop Emergy", since for the established conditions, the flow is very high for the numerical domain. Now, changing the "B 99 0" parameter to "B 99 1" achieves a greater number of iterations but returns to displaying the "Stop Emergy" critical speed message.
I would like to be able to help but first I need to solve this and know in what conditions they have to be given to avoid this message in the execution.

[kamath]

Hi Hernest,
The "Emergency stop" error appears due to unphysically large velocities in the domain. So, you would need to check where exactly things are going wrong in the tank.

B 99 0 stands for not having a beach. B 99 1 introduces a numerical beach.
One way to look at it, if your simulation works for the former and not for the latter, then may be there are extremely large reflections that then cause other problems?
You need to look at the output (vtu files) velocity fields to locate the source of the high velocities to determine what the problem is.

[raul_xusernamex_]

Originally Posted by kamath

Hi Hernest,
The "Emergency stop" error appears due to unphysically large velocities in the domain. So, you would need to check where exactly things are going wrong in the tank.

B 99 0 stands for not having a beach. B 99 1 introduces a numerical beach.
One way to look at it, if your simulation works for the former and not for the latter, then may be there are extremely large reflections that then cause other problems?
You need to look at the output (vtu files) velocity fields to locate the source of the high velocities to determine what the problem is.

Thaks so much for the orientation.
Thank you very much!