CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Waves2Foam Related Topics

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

Like Tree49Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   April 11, 2013, 10:05
Default
  #341
New Member
 
Join Date: Apr 2011
Posts: 28
Rep Power: 5
pepe.aero is on a distinguished road
Quote:
Originally Posted by ngj View Post
Hi Andrea,

Welcome to the waves2Foam part of cfd-online

Thanks for reporting the bug, I have made a note of it, though to be honest it might take me some time to figure out a solution, so if you come up with a solution please tell me.

Good luck

Niels

P.S. To begin with you can merely remove the addition of the internal structure in the Allrun script, i.e. the three lines related to the creation of the internal wall. In this way you can at least consider the waves in a 3 dimensional domain.

Hi Niels,

thanks for your very quicky reply!
As you suggested I have removed from Allrun the 3 lines regarding the internal walls but a critical segmentation fault is reported when I tryed to execute setWaveParameters command:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time


Reading g

Reading waveProperties


Constructing: stokesFirstProperties
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigSegv::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 std::basic_string<char, std::char_traits<char>, std::allocator<char> >::basic_string(std::string const&) in "/usr/lib/x86_64-linux-gnu/libstdc++.so.6"
#4 Foam::setWaveProperties::writeBeginning(Foam::Ostr eam&) in "/home/andrea/OpenFOAM/andrea-2.2.0/platforms/linux64GccDPOpt/lib/libwaves2Foam.so"
#5 Foam::stokesFirstProperties::set(Foam::Ostream&) in "/home/andrea/OpenFOAM/andrea-2.2.0/platforms/linux64GccDPOpt/lib/libwaves2Foam.so"
#6
in "/home/andrea/OpenFOAM/andrea-2.2.0/platforms/linux64GccDPOpt/bin/setWaveParameters"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8
in "/home/andrea/OpenFOAM/andrea-2.2.0/platforms/linux64GccDPOpt/bin/setWaveParameters"
Segmentation fault


Do you have any idea about the reason of this error looking to such an output?
Thanks for your support!

Regards,

Andrea
pepe.aero is offline   Reply With Quote

Old   April 11, 2013, 11:03
Default
  #342
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,552
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Andrea

Please see above. The error should be corrected in the most recent and static version of waves2Foam.

/ Niels
ngj is offline   Reply With Quote

Old   April 11, 2013, 12:54
Default
  #343
New Member
 
Join Date: Apr 2011
Posts: 28
Rep Power: 5
pepe.aero is on a distinguished road
Quote:
Originally Posted by ngj View Post
Hi Andrea

Please see above. The error should be corrected in the most recent and static version of waves2Foam.

/ Niels

Thanks Niels,

as was for Silvan setWaveParameters is now working well.
Sorry for not having noticed the thread before. Anyway, thanks for the support.

Best,

Andrea
pepe.aero is offline   Reply With Quote

Old   April 15, 2013, 06:32
Default
  #344
Member
 
Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 5
gregjunqua is on a distinguished road
Hi there i am beginning with with waves2Foam.
I got to do the compilation by myself. But when i do it it seems to don't have the repertory for the different script needed

Quote:
could not open file createFields.H for source file waveDyMFoam.C
could not open file alphaCourantNo.H for source file waveDyMFoam.C
could not open file alphaEqnSubCycle.H for source file waveDyMFoam.C
could not open file UEqn.H for source file waveDyMFoam.C
i guess it's in the options
so here is my options

Quote:
EXE_INC = \
-I.. \
-I$(LIB_SRC)/transportModels \
-I$(LIB_SRC)/transportModels/incompressible/lnInclude \
-I$(LIB_SRC)/transportModels/interfaceProperties/lnInclude \
-I$(LIB_SRC)/turbulenceModels/incompressible/turbulenceModel \
-I$(LIB_SRC)/finiteVolume/lnInclude \
-I$(LIB_SRC)/dynamicMesh/lnInclude \
-I$(LIB_SRC)/dynamicFvMesh/lnInclude \
-I$./home/gregoire/Program_file/openfoam/waves2Foam/src/lnInclude/ \
-I$(LIB_SRC)/meshTools/lnInclude \
-I$(LIB_SRC)/fvOptions/lnInclude \
-I$(LIB_SRC)/sampling/lnInclude \
-I$./home/gregoire/Program_file/openfoam/waves2Foam/src \
-lwaves2Foam

EXE_LIBS = \
-linterfaceProperties \
-ltwoPhaseInterfaceProperties \
-lincompressibleTransportModels \
-lincompressibleTurbulenceModel \
-lincompressibleRASModels \
-lincompressibleLESModels \
-lfiniteVolume \
-ldynamicMesh \
-ldynamicFvMesh \
-ltopoChangerFvMesh \
-lmeshTools \
-lfvOptions \
-lsampling
I know it's the bolted line how make trouble i tried several form (I must tell you that i am not a pro in wmake). If you could help me
-I./../../../../../src/lnInclude

don't look like to work
gregjunqua is offline   Reply With Quote

Old   April 15, 2013, 11:27
Default
  #345
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,552
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi,

You would need to provide us with some more information, e.g.:

- OF-version number?
- what do you want to achieve (waves2Foam is distributed with waveFoam)?
- where have you place waveDyMFoam relative to waveFoam?
- did the supported (downloaded) part of waves2Foam compile correctly?

Kind regards

Niels
ngj is offline   Reply With Quote

Old   April 16, 2013, 04:19
Default
  #346
Member
 
Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 5
gregjunqua is on a distinguished road
Hi, Thank you for your fast answer.
I am working on wind wave interactions, and now i am searching some tools for doing it with openFoam. I had done some experiment with fixed bottom but i would like to have a better fit with true condition.
The original allwmake don't work i use oF 2.2.0
I used the tutorial located at
http://openfoamwiki.net/index.php/Contrib/waves2Foam



waveDyMFoam is at this location :
/home/gregoire/Program_file/openfoam/waves2Foam/applications/solvers/solvers220/waveFoam
the src of waves2foam is here :
/home/gregoire/Program_file/openfoam/waves2Foam/src/lnInclude

thank you

Last edited by wyldckat; December 28, 2013 at 08:57. Reason: fixed broken link
gregjunqua is offline   Reply With Quote

Old   April 16, 2013, 05:21
Default
  #347
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,552
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
OK, first thing first.

Other users above have successfully compiled waves2Foam on OF220 with the standard Allwmake script, so could you elaborate on your problems?

Secondly, the line in your options file for waveDyMFoam, which looks like this

Code:
-I.. \
is the relative path for e.g. createFields.H, which is placed in the waveFoam directory. Make sure that this relative path is correctly defined.

Kind regards

Niels
ngj is offline   Reply With Quote

Old   April 17, 2013, 02:37
Default
  #348
Member
 
Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 5
gregjunqua is on a distinguished road
Hi Niels

It appear that at the first try to install the library
libgsl0-dev
wasn't correctly installed.
Thank you for your help and time
gregjunqua is offline   Reply With Quote

Old   April 17, 2013, 03:59
Default
  #349
Member
 
YS
Join Date: Jan 2010
Posts: 61
Rep Power: 6
Ya_Squall2010 is on a distinguished road
I am using waveFoam to resolve the divergent and transverse waves behind ship hull moving above a steady potential current (U=1 m/s). By following the tutorial case "squarepile", I used a circular domain with cylinder relaxationZone. As you can see from the static picture below, there's strong reflextion of waves travelling from wake to the water incoming direction. I need you kind help on how to reduce the unphysical reflexed waves. Many thanks!


reflextion.jpg
Ya_Squall2010 is offline   Reply With Quote

Old   April 17, 2013, 05:44
Default
  #350
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,552
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi,

It does indeed seem that you have too much reflection. The best explanation I can come up with is the fact that your waves have a large angle to the relaxation zone. Both active absorption and passive relaxation zone techniques have problems with damping oblique waves.

I know it would be sort of cheating, but what if you combined a couple of relaxation shapes (you might need to add a new shape-type), so your domain or at least relaxation layout looks something like the attached figure. This would give you a normal incident wave field and might limit the amount of reflection.

Kind regards,

Niels

P.S. Your figure does nevertheless look good, and it was especially for your type of application that I made the shape run-time selectable; though I have never tested it myself.
Attached Images
File Type: png compDomain.png (5.6 KB, 7 views)
ngj is offline   Reply With Quote

Old   April 22, 2013, 06:32
Default
  #351
Member
 
YS
Join Date: Jan 2010
Posts: 61
Rep Power: 6
Ya_Squall2010 is on a distinguished road
Hi, ngj,

Here's a quick update.

I managed to superimpose two circular relaxation zones as you suggested. Things are much better. I assume if impose an elliptical one at the wake region, the visible weak reflection in my simulation should be gone completely.

sigma.png

vof.jpg
ngj and sail like this.
Ya_Squall2010 is offline   Reply With Quote

Old   April 24, 2013, 03:36
Default Monopile
  #352
New Member
 
Jonas Kastrup
Join Date: Apr 2013
Posts: 3
Rep Power: 3
josk is on a distinguished road
Hi Niels.

I am analyzing the force and overturning moment on a monopile at different wave configurations.

The wave setup is done with stream function theory with the settings:
  • Wave height 13 m.
  • Water depth 21 m.
  • Period 10 s.
  • Eulerian drift velocity 1.5 m/s.

The free surface has been resolved with 2x2x2m cells.

The attached picture of the wave is at time 48 s where there has been 4 waveruns through the domain. It flows from left to right and the white boxes show the location of the relaxation zones.
As you can see on the picture the initial stream function wave in the relaxation zone starts to break and lose energy when it has traveled one to two wavelengths.
Is the physics of this behavior correct?
I understand that we are close to the range of shallow water, but shouldn't it only be possible for the stream function to generate the wave if the wave height to water depth where below the wave breaking point?

I have measured the force and overturning moment on the monopile (see attached image). For info the monopile is fully rigid and not moving. There is these short duration spikes in the force and moment which arise due to the wave breaking on the monopile. Are these physically correct or is the second spike with larger duration a more valid maximum force/moment?

I hope you can help answer the questions.

/Jonas
Attached Images
File Type: jpg wave.jpg (14.2 KB, 48 views)
File Type: jpg force_moment.jpg (49.5 KB, 41 views)
josk is offline   Reply With Quote

Old   April 24, 2013, 04:55
Default
  #353
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,552
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Jonas,

Not considering the current, you are more or less on the stability limit for waves, see e.g. any diagram for the validity range for mono-chromatic waves, which could be Fenton (1990), Nonlinear Wave Theories, chapter 1, in "From the Sea. Ideas and Observations on the Progress in the Study of the Seas". Wiley.

A co-directed current would probably stabilise the wave against breaking, whereas an opposing current would make the wave steeper, and thus lower the breaking criterion. But I am not sure.

Have you successfully run with smaller waves? And also, your domain is quite high compared to the "small" water depth.

Kind regards

Niels
ngj is offline   Reply With Quote

Old   April 24, 2013, 06:07
Default
  #354
New Member
 
Jonas Kastrup
Join Date: Apr 2013
Posts: 3
Rep Power: 3
josk is on a distinguished road
Hi Niels.

My domain is high because I did not change it after I ran with higher waves.
Just for info, the current is co-directed.

I think I have to run some tests on smaller waveheights to find the breaking limit for this waterdepth combined with the current.

Do you have any idea why the force/moment peak is so high?
josk is offline   Reply With Quote

Old   April 24, 2013, 14:01
Default
  #355
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,552
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Jonas,

With respect to the peaks in the force history, then I suppose it is related to the wave breaking. I have seen almost sinusoidal force histories under non-breaking waves.

Do recall that the near surface kinematics under breaking waves do have very large (in excess of the propagation speed) velocities in the surface roller, so this will give rise to large forces of short duration relative to the wave period. In your case the propagation speed is some 15 m/s.

Kind regards

Niels
ngj is offline   Reply With Quote

Old   April 24, 2013, 16:25
Default
  #356
Senior Member
 
kilroy's Avatar
 
Join Date: Mar 2013
Location: USA
Posts: 119
Rep Power: 3
kilroy is on a distinguished road
Niels,

Thank you very much for the waveFoam solver. I am not an expert on Linux and I am having some issues while modifying "interFoam" to "waveFoam". If you could help me with that, I would appreciate it.

I followed exactly the same steps described in the wiki and made the necessary changes to the waveFoam.C and other files (http://openfoamwiki.net/index.php/Co...m#Applications).

I am using OF v2.2.0. My "Make/files" file looks like below:

waveFoam.C

EXE = $(FOAM_USER_APPBIN)/waveFoam


And my "Make/options" file looks like below:

EXE_INC = \
-I$(LIB_SRC)/transportModels \
-I$(LIB_SRC)/transportModels/incompressible/lnInclude \
-I$(LIB_SRC)/transportModels/interfaceProperties/lnInclude \
-I$(LIB_SRC)/turbulenceModels/incompressible/turbulenceModel \
-I$(LIB_SRC)/finiteVolume/lnInclude \
-I$(LIB_SRC)/meshTools/lnInclude \
-I$(LIB_SRC)/fvOptions/lnInclude \
-I$(LIB_SRC)/sampling/lnInclude \
-DOFVERSION=220 \
-I./../../../../src/lnInclude

EXE_LIBS = \
-linterfaceProperties \
-ltwoPhaseInterfaceProperties \
-lincompressibleTransportModels \
-lincompressibleTurbulenceModel \
-lincompressibleRASModels \
-lincompressibleLESModels \
-lfiniteVolume \
-lmeshTools \
-lfvOptions \
-lsampling \
-L$(FOAM_USER_LIBBIN) \
-lwaves2Foam


The path of the solver is below:

/opt/openfoam220/applications/solvers/waves2Foam/applications/solvers/solvers220/waveFoam

After all those modifications that I described above, I go the directory where "waveFoam.C" is located and run "wmake". After that I get the error below:

SOURCE=waveFoam.C ; g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -I/opt/openfoam220/src/transportModels -I/opt/openfoam220/src/transportModels/incompressible/lnInclude -I/opt/openfoam220/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam220/src/turbulenceModels/incompressible/turbulenceModel -I/opt/openfoam220/src/finiteVolume/lnInclude -I/opt/openfoam220/src/meshTools/lnInclude -I/opt/openfoam220/src/fvOptions/lnInclude -I/opt/openfoam220/src/sampling/lnInclude -DOFVERSION=220 -I./../../../../src/lnInclude -IlnInclude -I. -I/opt/openfoam220/src/OpenFOAM/lnInclude -I/opt/openfoam220/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/waveFoam.o
waveFoam.C:48:28: fatal error: relaxationZone.H: No such file or directory
compilation terminated.
make: *** [Make/linux64GccDPOpt/waveFoam.o] Error 1


Can you please help me to solve that problem?

Thanks,

Kilroy
kilroy is offline   Reply With Quote

Old   April 24, 2013, 16:34
Default
  #357
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,552
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Kilroy,

Since you are using waves2Foam with 2.2.0., then waveFoam is distributed in the waves2Foam package. If you by change have an older version, delete it and download the most recent here:

http://www.student.dtu.dk/~ngja/waves2Foam.tar.gz

Once you execute the Allwmake script in the waves2Foam directory, everything should be compiled correctly (as long as you have the needed dependencies).

Kind regards and good luck

Niels
ngj is offline   Reply With Quote

Old   April 24, 2013, 17:31
Default
  #358
Senior Member
 
kilroy's Avatar
 
Join Date: Mar 2013
Location: USA
Posts: 119
Rep Power: 3
kilroy is on a distinguished road
Niels,

Thank you very much for the reply.

I did as you have told me and removed the older versions of "waves2Foam" and downloaded the most recent one. I ran "Allwmake" and everything is completed without any errors.

But when I try to run the tutorials, I see a message like below in the terminal before the run of "blockMesh":

sh: ../../../bin/prepareCase.sh: No such file or directory

Then again i can not see any errors but the tutorials complete the analysis less then a second and I can not see any results. In the tutorials folder there is only the folder "0" after the analysis. It appears like the mesh is generated but the solver is not running for some reason.

Thank you very much for your time,

Kilroy
kilroy is offline   Reply With Quote

Old   April 25, 2013, 02:36
Default
  #359
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,552
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Good morning Kilroy,

Good that you now have a compiled waves2Foam.

With respect to the tutorial, which one are you running? Or is the problem related to all of the tutorials?

If you are standing in the non-working tutorial and do

Code:
cd ../../../bin/
then you should be able to locate prepareCase.sh, which your system claims not to be able to find. If this is working, we probably need a bit of diagnostics, i.e. what is your operating system, what version, etc.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   April 25, 2013, 11:04
Default
  #360
Senior Member
 
kilroy's Avatar
 
Join Date: Mar 2013
Location: USA
Posts: 119
Rep Power: 3
kilroy is on a distinguished road
Good morning Niels,

I searched for my "prepareCase.sh" file and found it. Then I applied the changes in there manually and now all my tutorial cases are working.

But I have a question. I am sorry about my lack of knowledge of Linux. My "prepareCase.sh" file is located in:

/opt/openfoam220/applications/solvers/waves2Foam/bin

To be able to make the "prepareCase.sh" file work automatically how should I modify the line below?

cd ../../../bin/

What does "../" mean? And how many of "../"s should I put in front of the "/bin" to make "prepareCase.sh" file work automatically?

By the way, my operating system is ubuntu 12.04 LTS

Thank you very much,

Kilroy

Last edited by kilroy; April 25, 2013 at 11:24.
kilroy is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Other Topics at OpenFOAM Workshop Milan 2008 hjasak OpenFOAM 2 October 26, 2013 05:33
Sections / Topics in CFD Wiki Roberthealy1 CFD-Wiki 6 August 23, 2007 18:58
CFD Related Educational Programmes Jonas Larsson Main CFD Forum 3 February 9, 2007 11:11
project topics vivekanand CFX 0 October 27, 2004 06:17
Advanced Topics in Aerodynamics Antonio Filippone Main CFD Forum 0 August 28, 1999 13:16


All times are GMT -4. The time now is 03:20.