[ngj]

Hi Mohamad,

It appears that the solvers in waves2Foam has not been compiled correctly. Otherwise, the executable should be known.

Try to recompile waves2Foam and identify, whether any errors are reported.

The error using waves2Foam is related to the fact that the boundary condition waveAlpha is distributed along with libwaves2Foam.so, which interFoam does not know.

Kind regards,

Niels

P.S. I have just seen that you are using 2.0.1. No solvers are distributed with waves2Foam for that version. Please follow the guidelines on the wiki for creating the waveFoam solver.

Link: http://openfoamwiki.net/index.php/Contrib/waves2Foam

[CFD-Palma]

Dear Niels,

I am working with a modification of the square pile tutorial. The idea is replacing the pile with a boat stl, so I have added the snappy and the dynamic mesh, so I use waveDyMFoam. I have tried many diferent meshes, the last a millon + cels that checkMesh OK.
No wonder what I do it always ends very soon with FPE. I wander if you could give me a clue as to where to look for the problem, as the courant nº is ok, tried different speeds, etc.
I attach the log file and the terminal output.
Thanks in advance,
Carlos.

[ngj]

Hi Carlos,

Since it crashes so soon, it is probably related to the mesh motion and VOF coupling. Have you tried running it without waves? I would guess that is still crashes.

Kind regards,

Niels

[CFD-Palma]

Hi Niels, thanks for answering so soon.

The case has only current, no other wave type but potential current.
Attached a jpg of a run that lasted longer but ended the same way. Since it runs with no mesh motion, it could be related to some issue with the mesh, it may need more room on top to deform freely, I do not know.
Thanks for your attention,
Carlos.

[ngj]

You say it runs without mesh motion, but the ship you have depicted, isn't that being moved? Or have you meshed it in that skewed position?

- Niels

[CFD-Palma]

Hi Neils, The picture is the case "with mesh motion" that has run for a short while, allowing me to have a look at what happens. (Short in real time, but it took hour running)
The goal of the simulation is to predict the keel angle and position on the generated wave. I am not much interested in the forces (by now) there fore laminar should be adequate.
Regards,
Carlos.

[ngj]

OK. Mesh motion can produce some results from time to time, which are really hard to debug.

A couple of thoughts, which might help:

1. Use a really long ramping time (Tsoft), so the ship motion is not subject to shocks.

2. Use a solution from a potential flow solver with free surface capabilities (do not know of any, though) to initialise the flow field around the ship. Again, avoiding too many shocks.

3. Try with a simple geometry first, e.g. a simplified ship being made of a box and a prism with a triangular cross section. Try to begin with a stationary geometry, i.e. no mesh motion.

4. Look up shipFoam in the forum and couple it with waves2Foam. Unfortunately, the svn/git for the OpenFoam-extend branch is still down, but maybe you can find another source.

5. Some time ago navalPack was released. I have never used it, but the name suggests that it can prove useful.

Good luck,

Niels

[CFD-Palma]

Hi Niels,

Thank you very much for your suggestions. I will try ship-foam, that I have downloaded some time ago but was reluctant to try because it works on OF 1.6 which has some differences with 2.1 I have been using.
In any case I will try some of your ideas to learn and how knows, it may work.

Thank you again for your time and help.
Best regards,
Carlos.

[ngj]

I have seen a thread, where a modification of shipFoam to 2.1 was discussed. Try search for it, i.e. use Google on the cfd-online domain

Good luck,

Niels

[ngj]

Dear all,

We are currently working on an extension to waves2Foam (hopefully coming out in the late spring/summer of 2013). The key extension is a large post-processing utility, so one does not need to transform the output data into matlab, python or whatever, but you can work with the post-processing directly in OpenFoam.

I have a small question to you, namely: What type of post-processing are you doing on your data?
This is for us to be inspired of which tools to put into the extension.

To give you a flavour, currently the following type of post-processing is implemented and tested:

1. Read probes/surface elevation gauges into the post-processing utility (e.g. velocity, surfaceElevation, forces and moments, void fraction). This includes a concatenation of data from restarted simulations.

2. Perform interpolation to equidistant time series if needed (e.g. for spectral analysis).

3. Output the data in easily accessible ASCII format, i.e. pure numbers in a rectangular matrix format.

4. Spectral analysis (two versions with one for regular and one of irregular waves). Can be carried out on both scalar and vector quantities.

5. Reflection analysis based on the surface elevation signal. Works currently in 2D with two different methods for regular and irregular waves.

6. Ensemble averaging of regular wave quantities (scalar and vector)

7. Trapezoidal time integration and cumulative trapezoidal time integration (scalar and vector).

Any suggestions beyond these methods are highly appreciated.

If extensive implementation is needed, please also provide me with a test case, a mathematical formulation of the problem and a set of results from the analysis to use as validation. Furthermore, a transparent implementation in any language would be nice.

Kind regards,

Niels

[CFD-Palma]

Hi Niels,

Very interesting, will keep on waiting!!

I have been trying to apply wave2foam to fast speed boat analysis, and finally manage to make it stable.
Now, on Fr nº over 1, the shock of the initial acceleration is to hight so I have been wandering if there is a possibility of making the speed of the current wave as function of time, at least in two steps.

In your previous answer to my first questions, you pointed this up and suggested a hight tsoft. I could not get results with this or did not do it properly, because I get an elevation of the water surface in the inlet and the simulation crashes. If the default is 0, what could be a number to try?

Thanks you for your work and help!
Carlos.

[ngj]

Hi Carlos,

You can implement any inlet velocity of your liking. Look into the waveTheory folder in the source, and you will quickly be able to make a time varying current.

With respect to the soft start/ramping time, maybe you will need another shape function in the relaxation zone, but I honestly do not know, since I have never tried such types of simulations.

Kind regards,

Niels

[Sagun]

Hello Niels,

I have been trying to install waves2Foam on my personal computer and whenever I try to obtain the source code via SVN I get this error message:

svn: Could not open the requested SVN filesystem

Kindly advise.

Regards,
Sagun

[ngj]

Hi Sagun,

The SVN is currently down (http://openfoamwiki.net/index.php/Main_Page), so meanwhile you can find a static version of the most recent revision here:

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

Kind regards,

Niels

[lucaBonfiglio]

Originally Posted by jordi.muela

Hi,

I've used a constraint plane XY and moments of inertia very highs in the three axis (that's the reason of the small pitch), to assure the stability in this simulation test.

Hi foamers,
I've been using waveDyMFoam with sixDoFRigidBodyDisplacement for simulate a ship in waves with zero speed. I set up the inertia and the mass without constraints, but the angular velocities become very high (especially for roll).
As waves come from the bow, I shouldn't have roll.
I tried to constraint roll and yaw, but the constraint moments grow up very quickly blowing up the simulation. I also tried to increase the inertia but I only delay the blowing up of the simulation.
Does anybody have any suggestion?
Thank you

Luca

[lucaBonfiglio]

Hi Niels,

I've been using waveDyMFoam with sixDoFRigidBodyDisplacement for simulate a ship in waves with zero speed. I set up the inertia and the mass without constraints, but the angular velocities become very high (especially for roll).
As waves come from the bow, I shouldn't have roll.
I tried to constraint roll and yaw, but the constraint moments grow up very quickly blowing up the simulation. I also tried to increase the inertia but I only delay the blowing up of the simulation.
Do you have any suggestion?
Thank you

Luca

[rosswin]

Hi Niels,

I have successfully run the waveFlume tutorial. I have openfoam 160 and ubuntu 12.10.

I am trying to implement a gmsh mesh into the waveFlume tutorial. I have managed to edit all the files required to the best of my knowledge, i.e. the 0 time folder and boundary file.

I'm not really sure what to do once I have used the command: gmshToFoam.
I have checked the mesh using: checkMesh and that says the mesh is ok.

Should I then compile as normal using setWaveField and waveFoam or are there more steps that I need to do.

Here is the output I get when I try run setWaveFields:

Create time

Create mesh for time = 0



cannot open file

file: /home/ross/OpenFOAM/waves2Foam/tutorials/waveFoam/waveFlumeFoil3/system/fvSchemes at line 0.

    From function regIOobject::readStream()
    in file db/regIOobject/regIOobjectRead.C at line 62.

FOAM exiting

Thank you very much for your time and help.
Ross.

[ngj]

Hi Ross,

It seems that you are missing your fvSchemes file. This could be due to the fact that you rely on one of the tutorials, this also means that you have a lot of fvSchemes.* and fvSolution.* files available, though, you need to choose the one called "16". All of these things are in the tutorials taken care of in a small script, which you can find here:

<some path>/waves2Foam/bin/prepareCase.sh

Good luck

Niels

[Ya_Squall2010]

Originally Posted by winden

Hi Niels and Dominic. Just to report some findings, I tried the modified stokesFirst on a very fine mesh and found a discrepancy in amplitude of <1% and a discrepancy in encounter period of <0.01%.

combinedWaves with the same setup gives a 13.5% discrepancy in amplitude and 5% discrepancy in encounter period.

It would be nice if the same could be achieved with combinedWaves though since it would allow for forward speed simulations with all the other wave types as well as Stoke's first without having to derive new expressions for all of them.

I suspect that the problem might be that U and alpha are corrected separately and while combinedWaves adds the velocity contributions, it doesn't add the translation of the free surface. This means that the "wrong" omega hast to be set to ensure that the surface behaves correctly. This in turn means that the waves get the wrong velocity distribution though which may be why they grow when set free from the relaxation.

I may just be rambling here since I'm also confused by this so please correct me if i got something wrong. That would mean that the problem could be solved by adding a translation of eta in potentialCurrent somehow so that this behaviour in the modified stokesFirst

scalar eta = H_ / 2.0 * Foam::cos(omega_ * time - (k_ & x) + (k_ & current_) * time 
+ phi_) * factor(time) + seaLevel_;

is replicated when adding a Stoke's first wave with a potential current.

//Björn

Hello everybody,

I've been playing around with stokesFifth and stokesFirstwCurrent for a while, and wondering what is the correlation between the stokes drift velocity (scalar Q_) in the former and the current (vector current_) in the latter. I made a comparison run by setting Q_ and current_ to values of the same magnitude (-0.5 m/s) and the results are attached below and the wave properties for respective waves are also given for your reference:

stokesFifthCoeffs
{
waveType stokesFifth;
height 0.1;
period 2;
depth 0.4;
stokesDrift -0.5;
direction ( 1 0 0 );
Tsoft 2;
phi 0;
waveNumber (2.54201 0 0);
waveLength 2.47174;
omega 3.14159;
}

stokesFirstwCurrentCoeffs
{
waveType stokesFirstwCurrent;
Tsoft 2;
depth 0.4;
period 2;
direction ( 1 0 0 );
phi 0;
height 0.1;
waveNumber (1.70048 0 0);
waveLength 3.69495;
omega 3.14159;
current (-0.5 0 0);
}

surfaceElevation_wave_current.jpg

velocity_wave-current.jpg

Can anyone advise which one gives better result and what is the difference between Q_ and current_?

BTW, I have created a setProperty class for the stokesFirstwCurrent and attached here for your convenience.

stokesFirstwCurrentProperties.tar.gz

Many thanks!

[JGadelho]

Hello, everyone.
I'm new in OpenFOAM and waves2foam and I've been doing the tutorials of waves2foam successfully in the OF 2.1.1 version.

Today I installed OpenFOAM 2.2.0 and tried to install waves2foam on it running the ./Allwmake command, but gives me a lot of errors:

SOURCE=sampling/sampledSurfaceElevation/sampledSurfaceElevation.C ; g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -DOFVERSION=220 -I/opt/openfoam220/src/finiteVolume/lnInclude -I/opt/openfoam220/src/meshTools/lnInclude -I/opt/openfoam220/src/sampling/lnInclude -I/opt/openfoam220/src/lagrangian/basic/lnInclude -I/usr/local/include -I/usr/include -IlnInclude -I. -I/opt/openfoam220/src/OpenFOAM/lnInclude -I/opt/openfoam220/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/sampledSurfaceElevation.o
In file included from sampling/sampledSurfaceElevation/sampledSurfaceElevation.H:51:0,
from sampling/sampledSurfaceElevation/sampledSurfaceElevation.C:27:
/opt/openfoam220/src/sampling/lnInclude/sampledSet.H:47:22: fatal error: coordSet.H: No such file or directory
compilation terminated.
make: *** [Make/linux64GccDPOpt/sampledSurfaceElevation.o] Error 1
./Allwmake: line 47: cd: applications/solvers/solvers220: No such file or directory
make[1]: Entering directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/misc'
make[2]: Entering directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/misc/matlab'
make[3]: Entering directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/misc/matlab/postprocessing'
make[3]: Nothing to be done for `application'.
make[3]: Leaving directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/misc/matlab/postprocessing'
make[3]: Entering directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/misc/matlab/preprocessing'
make[3]: Nothing to be done for `application'.
make[3]: Leaving directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/misc/matlab/preprocessing'
make[2]: Leaving directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/misc/matlab'
make[1]: Leaving directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/misc'
make[1]: Entering directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/postProcessing'
make[2]: Entering directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/postProcessing/surfaceElevation'
SOURCE=surfaceElevation.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/finiteVolume/lnInclude -I/opt/openfoam220/src/meshTools/lnInclude -I/opt/openfoam220/src/sampling/lnInclude -I/opt/openfoam220/src/lagrangian/basic/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/surfaceElevation.o
In file included from ./../../../../src/lnInclude/sampledSurfaceElevation.H:51:0,
from ./../../../../src/lnInclude/IOsampledSurfaceElevation.H:36,
from surfaceElevation.C:89:
/opt/openfoam220/src/sampling/lnInclude/sampledSet.H:47:22: fatal error: coordSet.H: No such file or directory
compilation terminated.
make[2]: *** [Make/linux64GccDPOpt/surfaceElevation.o] Error 1
make[2]: Target `/home/jorge/OpenFOAM/jorge-2.2.0/platforms/linux64GccDPOpt/bin/surfaceElevation' not remade because of errors.
make[2]: Leaving directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/postProcessing/surfaceElevation'
make[1]: *** [surfaceElevation] Error 2
make[1]: Target `application' not remade because of errors.
make[1]: Leaving directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/postProcessing'
make: *** [postProcessing] Error 2
make[1]: Entering directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/preProcessing'
make[2]: Entering directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/preProcessing/relaxationZoneLayout'
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/finiteVolume/lnInclude -DOFVERSION=220 -I./../../../../src/lnInclude -IlnInclude -I. -I/opt/openfoam220/src/OpenFOAM/lnInclude -I/opt/openfoam220/src/OSspecific/POSIX/lnInclude -fPIC -Xlinker --add-needed -Xlinker --no-as-needed Make/linux64GccDPOpt/relaxationZoneLayout.o -L/opt/openfoam220/platforms/linux64GccDPOpt/lib \
-lfiniteVolume -L/home/jorge/OpenFOAM/jorge-2.2.0/platforms/linux64GccDPOpt/lib -lwaves2Foam -lOpenFOAM -ldl -lm -o /home/jorge/OpenFOAM/jorge-2.2.0/platforms/linux64GccDPOpt/bin/relaxationZoneLayout
/usr/bin/ld: cannot find -lwaves2Foam
collect2: error: ld returned 1 exit status
make[2]: *** [/home/jorge/OpenFOAM/jorge-2.2.0/platforms/linux64GccDPOpt/bin/relaxationZoneLayout] Error 1
make[2]: Leaving directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/preProcessing/relaxationZoneLayout'
make[1]: *** [relaxationZoneLayout] Error 2
make[2]: Entering directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/preProcessing/setWaveField'
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/finiteVolume/lnInclude -DOFVERSION=220 -I./../../../../src/lnInclude -IlnInclude -I. -I/opt/openfoam220/src/OpenFOAM/lnInclude -I/opt/openfoam220/src/OSspecific/POSIX/lnInclude -fPIC -Xlinker --add-needed -Xlinker --no-as-needed Make/linux64GccDPOpt/setWaveField.o -L/opt/openfoam220/platforms/linux64GccDPOpt/lib \
-lfiniteVolume -L/home/jorge/OpenFOAM/jorge-2.2.0/platforms/linux64GccDPOpt/lib -lwaves2Foam -lOpenFOAM -ldl -lm -o /home/jorge/OpenFOAM/jorge-2.2.0/platforms/linux64GccDPOpt/bin/setWaveField
/usr/bin/ld: cannot find -lwaves2Foam
collect2: error: ld returned 1 exit status
make[2]: *** [/home/jorge/OpenFOAM/jorge-2.2.0/platforms/linux64GccDPOpt/bin/setWaveField] Error 1
make[2]: Leaving directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/preProcessing/setWaveField'
make[1]: *** [setWaveField] Error 2
make[2]: Entering directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/preProcessing/setWaveParameters'
g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -DOFVERSION=220 -I/opt/openfoam220/src/finiteVolume/lnInclude -I./../../../../src/lnInclude -I./../../../../src/lnInclude -I/usr/local/include -I/include -IlnInclude -I. -I/opt/openfoam220/src/OpenFOAM/lnInclude -I/opt/openfoam220/src/OSspecific/POSIX/lnInclude -fPIC -Xlinker --add-needed -Xlinker --no-as-needed Make/linux64GccDPOpt/setWaveParameters.o -L/opt/openfoam220/platforms/linux64GccDPOpt/lib \
-lfiniteVolume -lgsl -lgslcblas -L/home/jorge/OpenFOAM/jorge-2.2.0/platforms/linux64GccDPOpt/lib -lwaves2Foam -lOpenFOAM -ldl -lm -o /home/jorge/OpenFOAM/jorge-2.2.0/platforms/linux64GccDPOpt/bin/setWaveParameters
/usr/bin/ld: cannot find -lwaves2Foam
collect2: error: ld returned 1 exit status
make[2]: *** [/home/jorge/OpenFOAM/jorge-2.2.0/platforms/linux64GccDPOpt/bin/setWaveParameters] Error 1
make[2]: Leaving directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/preProcessing/setWaveParameters'
make[1]: *** [setWaveParameters] Error 2
make[1]: Target `application' not remade because of errors.
make[1]: Leaving directory `/home/jorge/Downloads/Software/waves2Foam/applications/utilities/preProcessing'
make: *** [preProcessing] Error 2
make: Target `application' not remade because of errors.

Am I missing anything? Does anyone tried to install waves2foam on OF 2.2.0?

Thank you!