CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > REEF3D

Coupling between FNPF and CFD

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 8, 2023, 07:43
Smile Coupling between FNPF and CFD
  #1
New Member
 
Mar
Join Date: Nov 2022
Posts: 5
Rep Power: 4
Mar_Vayreda is on a distinguished road
Good afternoon!

I'm trying to run a simple case of waves in a 2D tank, running first the FNPF model and at the end of the domain running CFD (coupling). The first simulation went well, I got the REEF3D_FNPF_State and added it to the CFD folder. Then I did the DiveMESH for the CFD and got the REEF3D_HDC_INPUT. I attached a photo of paraview where you can see the simulation done by FNPF and at the end of the domain, the mesh created by the CFD.

The problem appears when I run the mpirun command for the CFD. It only runs the first iteration and this message appears:

EMERGENCY STOP -- velocities exceeding critical value N 61

Also, I saw in this forum an example where you have to add the wave generation for the second simulation (CFD). So I added dirichlet wave generation but I don't really understand it, because it should be taking as input the output of the first simulation.

If someone can help me understand why you have to add the wave generation in the second simulation and why the emergency stop appears it would be of great help. I also attached the ctrl and control files.

Thank you!!

Screenshot from 2023-02-08 15-04-05.jpg

control_cfd.txt

ctrl_cfd.txt

ctrl_fnpf.txt

control_fnpf.txt
Mar_Vayreda is offline   Reply With Quote

Old   February 9, 2023, 18:05
Default
  #2
Super Moderator
 
Hans Bihs
Join Date: Jun 2009
Location: Trondheim, Norway
Posts: 403
Rep Power: 19
valgrinda is on a distinguished road
B 99 4 doesn't work for HDC, because no wave theory input is know prior to the simulation (at least the way it is implemented now). So please use B 99 3 (which is the one to use anyway).

Another comment: you are creating standing waves with FNPF, because your waves are too steep for AWA. So best to use B 99 1 or B 99 2 for FNPF with the HDC overlap somewhat away from the beach.

And one question: for your case, why do you use HDC? Why not generate the waves directly in CFD? I think that would be easier.
__________________
Hans Bihs
Team REEF3D
www.reef3d.com
valgrinda is offline   Reply With Quote

Old   February 15, 2023, 09:50
Smile Coupling between FNPF and CFD
  #3
New Member
 
Mar
Join Date: Nov 2022
Posts: 5
Rep Power: 4
Mar_Vayreda is on a distinguished road
Hi Hans!

Sorry for the late reply and thank you so much for your suggestions. I tried changing with B 99 3 and also adding the relaxation zone for the CFD and the numerical beach for the FNPF and it started running. The problem is that when I visualize it with paraview, I can see that the wave is badly generated (I attached a photo and ctrl and control files). Do you have any other suggestion for why this is happening?

Answer to your question: The purpose of this test, is being able to run more complex cases in the future in a faster way, so in the CFD area is where the development would be added.

ctrl_fnfp.txt

control_fnfp.txt

control_cfd.txt

ctrl_cfd.txt

MicrosoftTeams-image.jpg
Mar_Vayreda is offline   Reply With Quote

Old   February 18, 2023, 07:50
Default
  #4
Senior Member
 
Arun Kamath
Join Date: Nov 2014
Location: Trondheim, Norway
Posts: 265
Rep Power: 14
kamath is on a distinguished road
Unless you have surface piercing structures z_max=1.0 in control.txt for FNPF.
Also ten cells in the z-direction are sufficient.

Z-grids in FNPF need to be stretched to get a good dispersion. Check this paper by Pakozdi et al

A good start (and mostly works) is:
B 103 5
B 113 2.0
B 116 1.0

Check the tutorial files for FNPF in the folder you download from github.

In control.txt for the coupled CFD simulation:
Since you use a truncated domain in B 10, you don't need to specify this with H 21-22-23.

(And this is the cause of your "bad waves". I suspect if you check the console output, you have velocities of the order 0.000x)
You also have H 31-32 where H 32 points to a non existent location. Your FNPF simulation ended at 30s whereas you ask for decomp to 60s.
Divemesh seems to handle this without throwing an error, but you should be doing for clarity's sake atleast.

In ctrl.txt for CFD, while B 98 2 works, imo it defeats the purpose of the coupling when you hold off a part of the domain for wave generation and in addition introduce a modulated wave rather than the actual output at the coupling point.
Here I would recommend B 98 3.
Set I 241 to 20 s assuming you want to start the sim at 20 s.



control-fnpf:
C 11 6
C 12 3
C 13 3
C 14 7
C 15 21
C 16 3

B 2 1000 1 10
B 10 0.0 230.0 0.0 0.2 0.0 1.0

B 103 5
B 113 2.0
B 116 1.0


M 20 2
M 10 8


control-cfd: H 10 4, H 31 20.0 H 32 30.0

ctrl-cfd: I 241 20.0

This will work.
__________________
Arun
X years with REEF3D

Last edited by kamath; February 18, 2023 at 07:57. Reason: fixed some errors
kamath is offline   Reply With Quote

Old   February 21, 2023, 02:21
Smile Coupling between FNPF and CFD
  #5
New Member
 
Mar
Join Date: Nov 2022
Posts: 5
Rep Power: 4
Mar_Vayreda is on a distinguished road
Thankyou so much for your answer!

Your corrections did work. I only had to add a numerical beach at the end of CFD to avoid reflection.
Mar_Vayreda is offline   Reply With Quote

Reply

Tags
cfd, coupling, fnpf

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
REEF3D 23.01 & DIVEMesh 23.01 valgrinda REEF3D News & Announcements 0 February 1, 2023 09:33
REEF3D 22.11 & DIVEMesh 22.11 valgrinda REEF3D News & Announcements 1 November 11, 2022 08:30
Hydrodynamic coupling between FNPF and CFD modules in REEF3D FelixD REEF3D 17 September 20, 2022 18:51
REEF3D 22.03 & DIVEMesh 22.03 valgrinda REEF3D News & Announcements 0 March 26, 2022 08:50
Module Coupling Felix_Sp REEF3D 5 January 14, 2022 04:02


All times are GMT -4. The time now is 15:42.