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 Tree76Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   April 2, 2012, 10:09
Default
  #61
New Member
 
Jan Löhrmann
Join Date: Sep 2010
Posts: 21
Rep Power: 6
JanL is on a distinguished road
Hi,

my calculation for StokesFirst within the limits of the wave-theory just finished (wave-height=0.02. The rest of the parameters in waveProperties.org are left as in the standard-case). The results are unfortunately similar to the previous ones (see figure 5). There is still a considerable dissipation of the wave. I used 12 cells per wave-height.

Any clues? Am I missing something?

Regards

Jan
Attached Files
File Type: pdf figure5.pdf (81.0 KB, 47 views)
JanL is offline   Reply With Quote

Old   April 2, 2012, 10:30
Default
  #62
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,607
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi both of you

@Jan: Could you please upload the deep-water test case. It is impossible to say anything without it.

@Björn: A non-zero speed in the air phase is available in the "official" release; see previous discussion in this thread. Furthermore, could you attach your waveProperties file?

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   April 2, 2012, 10:44
Default
  #63
New Member
 
Jan Löhrmann
Join Date: Sep 2010
Posts: 21
Rep Power: 6
JanL is on a distinguished road
Hi Niels,

please find attached the deepwater-case. I used a slight grading in y-direction as the wave-height is already quite small (0.02m) and I wanted to avoid further refinement to achieve 12 cells per wave-height.

Thanks for your help!

Regards

jan
Attached Files
File Type: gz deepWater.tar.gz (8.8 KB, 43 views)
JanL is offline   Reply With Quote

Old   April 2, 2012, 10:45
Default
  #64
Member
 
Björn Windén
Join Date: Feb 2012
Location: University of Southampton UK
Posts: 37
Rep Power: 5
winden is on a distinguished road
Hi Niels.

Thank you for your reply. I think the air velocity option wasn't added when I downloaded it (mid January.) My waveProperties-file is attached.

//Björn
Attached Files
File Type: zip waveProperties.zip (841 Bytes, 53 views)
winden is offline   Reply With Quote

Old   April 2, 2012, 11:10
Default
  #65
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,607
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi

@Jan: You have not changed the wave properties even though you have changed the water depth. This is probably the reason for the bad results in deep water - run setWaveParameters.

@Björn: setWaveParameters does nothing, when it is "combinedWaves", so the wave number for your waves is computed as if there were no current; i.e. it is approximately a factor of two wrong! At present you have to compute the wave number by hand instead - I have no plans of implementing a automatic wave number computation for wave-current, as two solutions exists, thus the user should be aware of which solution is the one to be used.
What you experience in the simulation is the elongation of the waves once they become free of the explicit forcing in the beginning of the relaxation zone.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   April 2, 2012, 11:25
Default
  #66
Member
 
Björn Windén
Join Date: Feb 2012
Location: University of Southampton UK
Posts: 37
Rep Power: 5
winden is on a distinguished road
Ah. I didn't think about double checking the wave number after running setWaveParameters. Thank you very much Niels! This will hopefully solve my problems.

//Björn
winden is offline   Reply With Quote

Old   April 2, 2012, 11:39
Default
  #67
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,607
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Actually, if you use streamFunction as a wave theory, then wave-current interaction is inherently an option, merely choose e.g. stokes drift different from 0.

Currently, however, you can only manually set those wave parameters (using the matlab script in application/utilities/misc), since I have not been able to make a stable solution from within OF.

/ Niels
ngj is offline   Reply With Quote

Old   April 4, 2012, 05:41
Default
  #68
New Member
 
Jan Löhrmann
Join Date: Sep 2010
Posts: 21
Rep Power: 6
JanL is on a distinguished road
Hi Niels,

thanks a lot for your help! I simply ran the Allrun-Script where setWaveParameters isn't included initially. With this command included (and the missing keywords in waveProperties) everything is running just fine!
Awesome!

Regards

Jan
JanL is offline   Reply With Quote

Old   April 5, 2012, 08:51
Default
  #69
New Member
 
Jan Löhrmann
Join Date: Sep 2010
Posts: 21
Rep Power: 6
JanL is on a distinguished road
Sorry to bother again. I just wanted to apply waveFoam to a bit more complex case of mine. setWaveParameters went through fine but when I'm running setWaveField I get this messages:
Code:
Reading waveProperties

Reading g

Reading field alpha

Reading field U

Reading field p

Setting the wave field ...

OBS: noi is larger than 2. Possible error could occur: 3
OBS: noi is larger than 2. Possible error could occur: 3
OBS: noi is larger than 2. Possible error could occur: 4
OBS: noi is larger than 2. Possible error could occur: 4
OBS: noi is larger than 2. Possible error could occur: 3
OBS: noi is larger than 2. Possible error could occur: 4
Is that an waveFoam-related error? What is "noi" standing for?

Best regards

Jan
JanL is offline   Reply With Quote

Old   April 5, 2012, 10:04
Default
  #70
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,607
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Jan

Yes, it is a waves2Foam problem, however, the bug has been fixed as of revision 1944, so please make an
Code:
svn update
and recompile.

Happy Easter,

Niels
ngj is offline   Reply With Quote

Old   April 9, 2012, 21:23
Default
  #71
New Member
 
Nima
Join Date: Feb 2012
Location: Perth, Western Australia
Posts: 13
Rep Power: 5
nima3906m is on a distinguished road
Dear Niels

I am using ubuntu 11.10 coupled with OF1.7 at the moment and read about the previous problems compiling waveFoam, although I have downlowded the latest release of the code I still get an error which means the solver (not the utilities) can not be compiled , the error message is sth like below :

root@ubuntu:/home/nima/waves2Foam# ./Allwmake
'/root/OpenFOAM/root-1.7.0/lib/linux64GccDPOpt/libwaves2Foam.so' is up to date.
make[1]: Entering directory `/home/nima/waves2Foam/applications/solvers/solvers17/waveFoam'
g++ -m64 -Dlinux64 -DWM_DP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-40 -I/opt/openfoam170/src/transportModels -I/opt/openfoam170/src/transportModels/incompressible/lnInclude -I/opt/openfoam170/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam170/src/turbulenceModels/incompressible/turbulenceModel -I/opt/openfoam170/src/finiteVolume/lnInclude -DOFVERSION=17 -I./../../../../src/lnInclude -IlnInclude -I. -I/opt/openfoam170/src/OpenFOAM/lnInclude -I/opt/openfoam170/src/OSspecific/POSIX/lnInclude -fPIC Make/linux64GccDPOpt/waveFoam.o -L/opt/openfoam170/lib/linux64GccDPOpt \
-linterfaceProperties -lincompressibleTransportModels -lincompressibleTurbulenceModel -lincompressibleRASModels -lincompressibleLESModels -lfiniteVolume -llduSolvers -L/root/OpenFOAM/root-1.7.0/lib/linux64GccDPOpt -lwaves2Foam -lOpenFOAM -liberty -ldl -lm -o /root/OpenFOAM/root-1.7.0/applications/bin/linux64GccDPOpt/waveFoam
/usr/bin/ld: cannot find -llduSolvers
collect2: ld returned 1 exit status
make[1]: *** [/root/OpenFOAM/root-1.7.0/applications/bin/linux64GccDPOpt/waveFoam] Error 1
make[1]: Leaving directory `/home/nima/waves2Foam/applications/solvers/solvers17/waveFoam'
make: *** [waveFoam] Error 2
make: Target `application' not remade because of errors.
make[1]: Entering directory `/home/nima/waves2Foam/applications/utilities/misc'
make[2]: Entering directory `/home/nima/waves2Foam/applications/utilities/misc/matlab'
make[2]: Nothing to be done for `application'.
make[2]: Leaving directory `/home/nima/waves2Foam/applications/utilities/misc/matlab'
make[1]: Leaving directory `/home/nima/waves2Foam/applications/utilities/misc'
make[1]: Entering directory `/home/nima/waves2Foam/applications/utilities/preProcessing'
make[2]: Entering directory `/home/nima/waves2Foam/applications/utilities/preProcessing/relaxationZoneLayout'
make[2]: `/root/OpenFOAM/root-1.7.0/applications/bin/linux64GccDPOpt/relaxationZoneLayout' is up to date.
make[2]: Leaving directory `/home/nima/waves2Foam/applications/utilities/preProcessing/relaxationZoneLayout'
make[2]: Entering directory `/home/nima/waves2Foam/applications/utilities/preProcessing/setWaveField'
make[2]: `/root/OpenFOAM/root-1.7.0/applications/bin/linux64GccDPOpt/setWaveField' is up to date.
make[2]: Leaving directory `/home/nima/waves2Foam/applications/utilities/preProcessing/setWaveField'
make[2]: Entering directory `/home/nima/waves2Foam/applications/utilities/preProcessing/setWaveParameters'
make[2]: `/root/OpenFOAM/root-1.7.0/applications/bin/linux64GccDPOpt/setWaveParameters' is up to date.
make[2]: Leaving directory `/home/nima/waves2Foam/applications/utilities/preProcessing/setWaveParameters'
make[1]: Leaving directory `/home/nima/waves2Foam/applications/utilities/preProcessing'

Please help me finding out what should be the reason

many thanks

Nima
nima3906m is offline   Reply With Quote

Old   April 10, 2012, 06:19
Default
  #72
New Member
 
Jan Löhrmann
Join Date: Sep 2010
Posts: 21
Rep Power: 6
JanL is on a distinguished road
Hi Niels,

I updated and recompiled waveFoam. The error I mentioned last time doesn't occur anymore. What confuses me though, is that when running setWaveFields, it doesn't seem to finish the actual setting of the wave field. This command is already running for hours, which doesn't seem to be reasonable to me (although I'm running it on a mesh with 4.1 M cells).
Have you a clue what could be the reason for that?

Best regards

Jan
JanL is offline   Reply With Quote

Old   April 11, 2012, 04:15
Default
  #73
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,607
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi

@Nima: I apologise, but there was a small bug in the Allwmake script - puzzling that it has not been noticed before. The svn is update, merely type
Code:
svn update
in the waves2Foam directory.

@Jan: Check whether or not the setWaveField is stalling in one of the cells or whether it is just extremely slow on your computer - I have set >1M cells in minutes, where most of the time was spent on constructing the mesh.
Stalling: E.g. an Info statement the loop in
Code:
void setWaveField::correct()
in src/setWaveField/setWaveField.C

Best regards,

Niels
ngj is offline   Reply With Quote

Old   April 11, 2012, 15:36
Default
  #74
Member
 
Join Date: May 2009
Posts: 51
Rep Power: 8
gfilip is on a distinguished road
I had a question about the relaxation zones introducing high-frequency content with a modified waveFlume tutorial case where the domain was extended, the depth increased to give kh=4. When I refined the tutorial grid resolution using refineMesh, I saw strange behavior containing high frequencies near the beginning of the outlet relaxation zone. Here is Neils' reply to my question:

Quote:
My only suggestion is that since the relaxation scheme only considers space, then the distance over which, the relaxation occurs, decreases with the increasing number of cells, as the time step goes down, i.e. a higher number of relaxations per cell and the relaxation goes like k^N * value, where k < 1 and value is e.g. the velocity. Thus in the extreme case, the relaxation zone becomes a reflecting wall. Therefore, a Courant based relaxation approach could possibly help you.
To follow up my question to Neils:
Does this mean that coarsening the outlet relaxation zone help with this issue?


Greg
gfilip is offline   Reply With Quote

Old   April 17, 2012, 11:11
Default
  #75
Member
 
Björn Windén
Join Date: Feb 2012
Location: University of Southampton UK
Posts: 37
Rep Power: 5
winden is on a distinguished road
Hi Niels.

Regarding your previous reply to my problems with using waves and current:
Waves2Foam Related Topics (post #65)

I'm still having some trouble with this. Would it be possible to enlighten me as to what would be the correct way of using waveProperties when modelling an object travelling with forward speed in waves.

I have tried setting T and k in the wave properties to those of the wave without current and omega to ωe =ωwave+kU.

This seems to give the correct encounter period and the wave keeps it's wavelength when it leaves the relaxation zone. However; there is still something that is not right in that the wave increases slightly in amplitude when leaving the relaxation zone (see attached image, red line is end of zone and propagation direction is left->right.)

Kind regards
Björn
Attached Images
File Type: jpg elevation.jpg (21.9 KB, 28 views)

Last edited by wyldckat; December 28, 2013 at 08:50. Reason: updated the link to the other post, given the transfer of posts from the news thread
winden is offline   Reply With Quote

Old   April 17, 2012, 16:22
Default
  #76
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,607
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Bjørn

I have never worked with such a problem, so I can hardly give a best practice help, but this question might solve your problems:

What wave theory are you using and is it valid?

You could e.g. test with a streamfunction wave and see, whether the wave height is correctly computed, because that theory must be valid for regular waves/regular waves superimposed on a current.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   April 18, 2012, 01:36
Default Compilation Error
  #77
New Member
 
Nima
Join Date: Feb 2012
Location: Perth, Western Australia
Posts: 13
Rep Power: 5
nima3906m is on a distinguished road
Hi Niels

thanks for updating the Svn but unfortunately I still have the problem compiling the code. I even started from the scratch compiled the OpenFOAM again myself to solve it but it still gives me the same error , I follow the wiki procedure as below , maybe you can find a mistake in what I do :

1. I download the SVN code (waves2foam)
2- create a folder named solver21 in solvers dictionary
3.paste the interFoam source code in it ( a folder named waveFoam)
4.change interFoam according to wiki
5.run ./Allwmake in the Waves2foam

It is really making me sick as I have been struggling with it nearly 2 weeks now and can't solve the problem, please tell me how do you interpret the following error message!:

it is a part of the message which contains the error

-I. -I/home/nima/OpenFOAM/OpenFOAM-2.1.x/src/OpenFOAM/lnInclude -I/home/nima/OpenFOAM/OpenFOAM-2.1.x/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/waveFoam.o
/bin/sh: cannot open The: No such file
make[1]: *** [Make/linux64GccDPOpt/waveFoam.o] Error 2
make[1]: Target `/home/nima/OpenFOAM/nima-2.1.x/platforms/linux64GccDPOpt/bin/waveFoam' not remade because of errors.
nima3906m is offline   Reply With Quote

Old   April 18, 2012, 05:35
Default Waves2Foam Related Topics
  #78
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,607
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
The problem is in your Make/options file. You need to manually add the two first two digits replacing the long string in the brackets.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   April 18, 2012, 05:37
Default
  #79
Senior Member
 
stephane sanchi
Join Date: Mar 2009
Posts: 300
Rep Power: 9
openfoam_user is on a distinguished road
Hi,

I try to redo the 3Dwaves tutorial with OF-2.1.x.

I'm having problem with the createBaffles command.

Regards,
Stephane.
openfoam_user is offline   Reply With Quote

Old   April 19, 2012, 05:54
Default
  #80
Senior Member
 
stephane sanchi
Join Date: Mar 2009
Posts: 300
Rep Power: 9
openfoam_user is on a distinguished road
How to create baffles for the 3Dwaves tutorial with OF-2.1.x.

# Create internal wall
topoSet
setsToZones -noFlipMap
createBaffles -overwrite f0 '(internalWall internalWall)'

As attachment is the topoSetDict file.

Hope it can help.

Regards,
Stephane.
Attached Files
File Type: txt topoSetDict.txt (1.2 KB, 21 views)
openfoam_user 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 04:33
Sections / Topics in CFD Wiki Roberthealy1 CFD-Wiki 6 August 23, 2007 17:58
CFD Related Educational Programmes Jonas Larsson Main CFD Forum 3 February 9, 2007 11:11
project topics vivekanand CFX 0 October 27, 2004 05:17
Advanced Topics in Aerodynamics Antonio Filippone Main CFD Forum 0 August 28, 1999 12:16


All times are GMT -4. The time now is 18:24.