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

[OLAFLOW] The OLAFOAM Thread

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

Like Tree16Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 30, 2020, 17:04
Default Reflections from the structure
  #241
New Member
 
paaresh
Join Date: Mar 2009
Posts: 19
Rep Power: 15
sri31049 is on a distinguished road
Hi Pablo


I am simulating the wave interaction with rectangular structure using olaFlow. In my simulation generated waves are spoiling due to wave reflections from the structure. Is there any ways to suppress those reflections. I am using cnoidal theory for 1 m wave height with 10 seconds period at 6m water depth.



Many thanks




regards
Sri
Attached Images
File Type: png figure_1.png (65.9 KB, 28 views)
sri31049 is offline   Reply With Quote

Old   May 4, 2020, 07:43
Default
  #242
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 18
Phicau is on a distinguished road
Hi Sri,
I do not have enough information to even understand what or even if something is wrong. I suggest that you take a look at this thread before posting: How to give enough info to get help

What I do notice is that the waves are at 2.7 m water depth instead of at 6 m, so they are quite nonlinear.



Best,
Pablo
__________________
Check out my new project: olaFlow --> The olaFlow Support Thread
Phicau is offline   Reply With Quote

Old   September 7, 2020, 10:39
Default Tracking surface
  #243
New Member
 
Sergio Croquer
Join Date: Jan 2011
Posts: 15
Rep Power: 14
sercro is on a distinguished road
Hello all,

I was wondering how do you track the interface height at specific locations?

I'm doing a test in OF2006 using interfaceHeight (which worked well in OF5) but it does not work when running in parallel, i.e. the output file does not have any values (it seems to work fine in serial)

I used to track the interface profiles along vertical lines and then interpolate the surface position in post-processing, but this approach is less accurate and generates a lot of files so it is not practical for big cases.

Thanks
sercro is offline   Reply With Quote

Old   September 7, 2020, 19:47
Default
  #244
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 18
Phicau is on a distinguished road
Hi Sergio,
I use the vertical lines method: https://github.com/phicau/olaFlow/tr...mpling/python2
I agree with you that the number of files it outputs is enormous, but not in that it is less accurate. It is the most flexible method available, allowing you to calculate the interface location in different ways (threshold[s], compressing volume...) and providing you with full control and possibility to obtain subgrid resolution.
What I usually do is include another small piece of code to the python script to delete those files that have already been processed.



Best,
Pablo
__________________
Check out my new project: olaFlow --> The olaFlow Support Thread
Phicau is offline   Reply With Quote

Old   September 7, 2020, 20:33
Default
  #245
New Member
 
Sergio Croquer
Join Date: Jan 2011
Posts: 15
Rep Power: 14
sercro is on a distinguished road
Quote:
Originally Posted by Phicau View Post
Hi Sergio,
I use the vertical lines method: https://github.com/phicau/olaFlow/tr...mpling/python2
I agree with you that the number of files it outputs is enormous, but not in that it is less accurate. It is the most flexible method available, allowing you to calculate the interface location in different ways (threshold[s], compressing volume...) and providing you with full control and possibility to obtain subgrid resolution.
What I usually do is include another small piece of code to the python script to delete those files that have already been processed.



Best,
Pablo
Hello Pablo,

Thanks for the quick reply. I said it's less accurate because the resulting profile was not as smooth as using interfaceHeight on the same case. Perhaps the difference came from my matlab script for the interpolation.

Definitely the script you shared will come handy. I'll try to tune so it deletes files regularly.

Thanks and have a nice day.

Sergio-
sercro is offline   Reply With Quote

Old   September 9, 2020, 12:20
Default
  #246
New Member
 
Sergio Croquer
Join Date: Jan 2011
Posts: 15
Rep Power: 14
sercro is on a distinguished road
Hello Pablo,

I'm going through your postSensVOF.py script and maybe I'm confused but I don't get why you integrate alpha in the z direction (line 81 onwards).

If I want to know the interface position, shouldn't I look for the position at which alpha changes from 0 to 1 (using some threshold)?
sercro is offline   Reply With Quote

Old   September 10, 2020, 18:46
Default
  #247
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 18
Phicau is on a distinguished road
Hi Sergio,
integrating the water column is one of the options to calculate the free surface elevation. It does not differ much from looking for the highest alpha = 0.5 location in most cases, but is especially useful to prevent jumps (for example when you get splashes) or when you get overturning waves (it is debatable where the free surface should be).
That is why I was suggesting that saving all the VOF data throughout the water column is the most complete approach, because then you can process it as you wish.


Best,
Pablo
__________________
Check out my new project: olaFlow --> The olaFlow Support Thread
Phicau is offline   Reply With Quote

Old   November 13, 2020, 10:43
Default wave-current interaction
  #248
Member
 
Callum Guy
Join Date: Dec 2019
Location: Scotland
Posts: 44
Rep Power: 5
CallumG is on a distinguished road
Hi Pablo,

I just wanted some more information regarding the superposition of waves and currents. As olaFlow is now, do orbital trajectories from surface waves propagate down into the water column?

And is there any documentation currently describing this wave-current superposition?

Great work on the toolbox!

All the best,
Callum
CallumG is offline   Reply With Quote

Old   November 15, 2020, 15:10
Default
  #249
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 18
Phicau is on a distinguished road
Hi Callum,
thanks!

I am not really sure on what you mean by the orbital trajectories propagating down into the water column. olaflow generates waves by setting the wave kinematics at the boundary, in that sense deep water waves have the expected decay of orbital velocities near the bottom, while shallow water waves have a constant profile throughout the water depth. Accordingly, intermediate water depth waves can also be generated.
At this stage the public release of wave-current superposition is a linear superposition but if you have any particular questions feel free to drop me an email.


Best,
Pablo
__________________
Check out my new project: olaFlow --> The olaFlow Support Thread
Phicau is offline   Reply With Quote

Old   January 11, 2021, 08:45
Default Problem when running setFields in parallel
  #250
Member
 
Haoran Zhou
Join Date: Nov 2019
Posts: 45
Rep Power: 5
Stan Zhou is on a distinguished road
Hi all,

As the porous media area which I wanted to set in my case is pretty large, I tried to run setFields in parallel but a problem occurred.

The commands I executed are as follows:
1. blockMesh
2. cp -r 0.org 0
3. decomposePar
4. mpirun -np 4 snappyHexMesh -parallel -overwrite
5. mpirun -np 4 setFields -parallel

During step 5, I got the FOAM FATAL IO ERROR:
Cannot find patchField entry for caisson
The error directs to the alpha.water files in every processor. However, I've already added the boundaryField 'caisson' in the alpha.water file of 0 folder before executing decomposePar.

Do I miss any steps before running setFields in parallel? Is there any methods to deal with this problem? Thanks in advance!

Best regards,
Stan
Stan Zhou is offline   Reply With Quote

Old   January 25, 2021, 21:22
Default
  #251
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 18
Phicau is on a distinguished road
Hi Stan,
I believe this issue is OpenFOAM related rather than due to olaFlow, so probably you will have more chances of it being solved by posting your own thread, if the solution below does not work.


I have not used snappy and setFields in parallel extensively, but I would not do step 2. Instead I would just create an empty zero folder (mkdir 0) and continue with 3 and 4. Then perform step 2 and try running "decomposePar -fields", which should only decompose the fields. Finally, end with step 5.


Best,
Pablo
__________________
Check out my new project: olaFlow --> The olaFlow Support Thread
Phicau is offline   Reply With Quote

Old   June 16, 2021, 21:54
Default Error-Tutorial- baseWaveFlumeNewAbs
  #252
New Member
 
Join Date: Sep 2020
Posts: 19
Rep Power: 4
paulathikalam is on a distinguished road
Hii all,

Could you help me out to fix this error.

I am trying to run the tutorial case of baseWaveFlumeNewAbs. I am running it in OpenFoam6. I have not modified anything in the tutorial case, just started the case by executing the runCase script.

The script runs fine till setfields, but when olaflow is executed it pops up an error, i think - related to active wave absoprtion
-----
Active wave absorption BC on patch outlet
"Initial water depths for absorption" 1( 0.4 )

Selecting active wave absorption model extendedRangeAWA
AWA model: patch outlet
Theory: extendedRangeAWA
Number of paddles: 1
Reference water depth: 1(0.4)
Wave period: 3
Updating extendedRangeAWA absorption model for patch outlet
----
full olaflow logfile is also attached.

Is this something to do with the OF version that i am using to run ? or is there something that i should edit ?

Thanks

Paul
Attached Files
File Type: txt olaFlow_log.txt (3.0 KB, 6 views)
paulathikalam is offline   Reply With Quote

Old   September 7, 2021, 04:41
Default Tide simulation
  #253
New Member
 
George
Join Date: Jun 2020
Location: Taiwan
Posts: 3
Rep Power: 4
CC_Joe is on a distinguished road
Hello, Phicau,

I want to simulate the interaction of tides on pore structures.

Only mean water level Change with time,without making wave.

Which boundary condition should I use ?

Maybe I just need to adjust setFieldDict or waveDict file to Reach this simulation ?

Any insight or knowledge would be greatly appreciated.
CC_Joe is offline   Reply With Quote

Old   September 7, 2021, 20:42
Default
  #254
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 18
Phicau is on a distinguished road
Hi George,
What is a tide but an extremely long wave? Set your wave height and wave period based on your local conditions and use Stokes I theory.

Best,
Pablo
__________________
Check out my new project: olaFlow --> The olaFlow Support Thread
Phicau is offline   Reply With Quote

Old   October 7, 2021, 06:31
Default kEpsilon Multiphase turbulence model
  #255
New Member
 
Rosangela
Join Date: Jul 2021
Posts: 4
Rep Power: 3
Rosy is on a distinguished road
Hello everyone!

I used Olaflow to simulate the current and the wave motion that impacts two vertical rigid cylinders partially immersed in the water, obtaining good results with a laminar simulation mode. Now I want to implement a similar model, but using kEpsilon multiphase turbulence mode.
I downloaded and compiled the Olaflow supplementary turbulenceMultiphaseLibrary and, taking the tutorial as a reference, I changed the boundary conditions of my case.
The model runs, but fails to create even the first time step.

Could someone help me figure out what am I wrong? I am attaching controlDict and other files.

Thank you.
Attached Files
File Type: txt controlDict.txt (1.4 KB, 2 views)
File Type: txt epsilon.txt (1.9 KB, 1 views)
File Type: txt k.txt (1.9 KB, 1 views)
File Type: txt nut.txt (1.6 KB, 1 views)
File Type: txt U_forw.txt (2.0 KB, 2 views)
Rosy is offline   Reply With Quote

Old   October 11, 2021, 21:04
Default
  #256
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 18
Phicau is on a distinguished road
Hi Rosy,
Try to save a snapshot before any instabilities arise and observe it in ParaView to see where they are generated. Perhaps is due to the mesh. Double check the numerical schemes and you can also try other seeding (initial) values.

Best,
Pablo
__________________
Check out my new project: olaFlow --> The olaFlow Support Thread
Phicau is offline   Reply With Quote

Old   October 17, 2021, 08:13
Post Question of focusing wave generation
  #257
New Member
 
yujia Wei
Join Date: Mar 2020
Posts: 1
Rep Power: 0
yujia111 is on a distinguished road
Hi Pablo,

Could I ask how to set wave phase focusing at a given time and a given location by using the static boundary wave generation bc?

The paper I tried to validate is "A blind comparative study of focused wave interactions with a fixed FPSO-like structure". It says a member function has to be added in IHFOAM to generate phase-focues wave. Could you provide any hints on how to modify the solver?


Many thanks,
Yujia
yujia111 is offline   Reply With Quote

Old   October 17, 2021, 21:19
Default
  #258
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 18
Phicau is on a distinguished road
Hi Yujia,
I have no idea regarding your question, this thread is on olaFlow.

In general you can calculate it using linear wave theory for a first order approximation, or use the methodology outlined in https://ascelibrary.org/doi/10.1061/...3-5460.0000420

Best,
Pablo
__________________
Check out my new project: olaFlow --> The olaFlow Support Thread
Phicau is offline   Reply With Quote

Old   November 12, 2021, 13:58
Default
  #259
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 5
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
Quote:
Originally Posted by Phicau View Post
Hi mleary29,

thanks for the report. You are right, at this point wave and current are just generated by linear superposition.

Please let me refer you to previous posts: The OLAFOAM Thread and The OLAFOAM Thread

Best,

Pablo
Hi Pablo,

I am so confused about the setting for currentWaveFlume (here current only considered without wave) since I tried many conditions regarding uCurrent (, , ) and setFieldsDict. i.g. 1st condition: keep outlet boundary condition of U zero and only give a value uCurrent in waveDict; 2rd condition: reversed 1st condition; 3rd conditon: set a uCurrent for both but zero in setFields. All situations I saw the movement of solid. Does "forward" means wave following wave and "back" means wave traveling against current? Why we need to set uCurrent in outlet insteady of inlet? Could you give me some tips? Thanks.
qi.yang@polimi.it is offline   Reply With Quote

Old   November 14, 2021, 12:32
Default
  #260
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 5
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
Hi, I am simulating wave-current flume and then I have to extract the dynamic wave pressure at the bottom. So, as you can see from the attached plot in which the pressure induced by wave only is reasonable whereas that of wave-current induced is so weird. What is your opinion? I do not have any ideas. Hope somebody could give me some tips, thanks.
Attached Images
File Type: jpg plot.jpg (57.3 KB, 9 views)
qi.yang@polimi.it is offline   Reply With Quote

Reply

Tags
generation, ihfoam, olafoam, waves

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
Divergence detected in AMG solver: k when udf loaded google9002 Fluent UDF and Scheme Programming 3 November 8, 2019 00:34
udf problem jane Fluent UDF and Scheme Programming 37 February 20, 2018 05:17
UDF velocity profile willroca Fluent UDF and Scheme Programming 2 January 10, 2016 04:13
Error messages atg enGrid 7 August 30, 2013 12:16
Phase locked average in run time panara OpenFOAM 2 February 20, 2008 15:37


All times are GMT -4. The time now is 05:30.