[Behrang]

Hi again,

I have also defined three different geometries in order to model different types of wave breakers (spilling, plunging and surging).
I want to compare the simulation wave run-up against theoretical run-up using the equation of hunt and battjes R2%=H*C*(surf similarity)
Do you have any suggestion for how to get the simulation results?

Best,
Behrang

[ngj]

Hallo Behrang,

With respect to the dissipation, have you remembered to change the water depth in the computational domain?

Secondly, with respect to mesh requirements to capture correct breaking and avoid immature breaking, I will kindly refer you to our paper (Jacobsen et al, 2012).

Kind regards,

Niels

[p.roberto]

Hi

I installed successfully waves2foam in OF221, and now i'm trying to compile waveDyMFoam but when i type wmake it gives me the error below:

SOURCE=waveDyMFoam.C ; g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -I.. -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/transportModels/twoPhaseMixture/lnInclude -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/transportModels -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/transportModels/incompressible/lnInclude -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/transportModels/interfaceProperties/lnInclude -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/turbulenceModels/incompressible/turbulenceModel -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/finiteVolume/lnInclude -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/dynamicMesh/lnInclude -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/dynamicFvMesh/lnInclude -I./../../../../../src/lnInclude -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/meshTools/lnInclude -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/fvOptions/lnInclude -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/sampling/lnInclude -IlnInclude -I. -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/OpenFOAM/lnInclude -I/home/tiago/OpenFOAM/OpenFOAM-2.2.1/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/waveDyMFoam.o
waveDyMFoam.C:44:28: erro fatal: relaxationZone.H: Ficheiro ou directoria inexistente
compilação terminada.
make: ** [Make/linux64GccDPOpt/waveDyMFoam.o] Erro 1

can anyone help me please?

[ngj]

Hi Roberto,

You are using a not updated guide on the wiki, which is also stated explicitly.

This means that your Make/options file is not correct. Consult e.g. the corresponding file for the waveFoam solver.

Kind regards

Niels

[p.roberto]

Hi Niels, thanks for your reply

I've try to compare the corresponding file for the waveFoam solver and i know that the problem is in:
-I./../../../../../src/lnInclude

but i don't know how i have to change that line

My Options file is:

EXE_INC = \
-I.. \
-I$(LIB_SRC)/transportModels/twoPhaseMixture/lnInclude \
-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./../../../../../src/lnInclude \
-I$(LIB_SRC)/meshTools/lnInclude \
-I$(LIB_SRC)/fvOptions/lnInclude \
-I$(LIB_SRC)/sampling/lnInclude

EXE_LIBS = \
-ltwoPhaseMixture \
-linterfaceProperties \
-ltwoPhaseProperties \
-lincompressibleTransportModels \
-lincompressibleTurbulenceModel \
-lincompressibleRASModels \
-lincompressibleLESModels \
-lfiniteVolume \
-ldynamicMesh \
-ldynamicFvMesh \
-ltopoChangerFvMesh \
-lmeshTools \
-lfvOptions \
-lsampling \
-L$(WAVES_LIBBIN) \
-lwaves2Foam

[ngj]

Hi Roberto,

For OF221, you want to add the following to the EXE_INC part:

    -DOFVERSION=221 \
    -DEXTBRANCH=0 \
    -DXVERSION=$(WAVES_XVERSION) \
    -I$(WAVES_SRC)/waves2Foam/lnInclude \
    -I$(WAVES_SRC)/waves2FoamSampling/lnInclude

and the following to the EXE_LIBS part:

    -L$(WAVES_LIBBIN) \
    -lwaves2Foam \
    -lwaves2FoamSampling

That should do it, but remember to add the needed newline characters '\'.

Kind regards,

Niels

[Nick_civ]

Hi Neils,

I am a new linux and OF user and i am trying to install waves2foam on OF 2.2.2.
I followed wiki instruction (svn co http://svn.code.sf.net/p/openfoam-extend/svn/trunk/Breeder_1.6/other/waves2Foam ) and a waves2foam folder has been created outside the OFfolder.
Then I followed the instructions for interfoam modifications and created waveFoam and waveDymfoam for 2.2.2. according to your 2.2.1 solver and wiki instructions.

-DOFVERSION=222 \ -DEXTBRANCH=0 \ -DXVERSION=$(WAVES_XVERSION) \ -I$(WAVES_SRC)/waves2Foam/lnInclude \ -I$(WAVES_SRC)/waves2FoamSampling/lnIncludewhen i apply wmake it shows me error although I tried
export WAVES_DIR=$HOME/waves2Foam
nikos@nikos:~$
nikos@nikos:~$ cd waves2Foam/
nikos@nikos:~/waves2Foam$ ./Allwmake

=====================================
ENVIRONMENTAL VARIABLES
=====================================
EXTBRANCH=0
FOAMEXTENDPROJECT=0
WAVES_APPBIN=/home/nikos/OpenFOAM/nikos-2.2.2/platforms/linux64GccDPOpt/bin
WAVES_DIR=/home/nikos/OpenFOAM/nikos-2.2.2/applications/utilities/waves2Foam
WAVES_GSL_INCLUDE=/usr/include
WAVES_GSL_LIB=/usr/lib64
WAVES_LIBBIN=/home/nikos/OpenFOAM/nikos-2.2.2/platforms/linux64GccDPOpt/lib
WAVES_POST=/home/nikos/OpenFOAM/nikos-2.2.2/applications/utilities/waves2Foam/applications/utilities/postProcessing
WAVES_PRE=/home/nikos/OpenFOAM/nikos-2.2.2/applications/utilities/waves2Foam/applications/utilities/preProcessing
WAVES_SOL=/home/nikos/OpenFOAM/nikos-2.2.2/applications/utilities/waves2Foam/applications/solvers/solvers222
WAVES_SRC=/home/nikos/OpenFOAM/nikos-2.2.2/applications/utilities/waves2Foam/src
WAVES_TUT=/home/nikos/OpenFOAM/nikos-2.2.2/applications/utilities/waves2Foam/tutorials
WAVES_UTIL=/home/nikos/OpenFOAM/nikos-2.2.2/applications/utilities/waves2Foam/applications/utilities
WAVES_XVERSION=0
WM_PROJECT_VERSION_NUMBER=222


FATAL ERROR.
The directory path /home/nikos/OpenFOAM/nikos-2.2.2/applications/utilities/waves2Foam does not exist.
Correct the path in bin/bashrc

This means that you should correct the environmental variable $WAVES_DIR
such that it is pointed (in absolute sense) to the path, where waves2Foam
is located.

Please note that once bin/bashrc is created, bin/bashrc.org is only an
inactive file. The latter is also the only of the two files, which is
updated through the SVN-repository.

EXITING

nikos@nikos:~/waves2Foam$ export WAVES_DIR=$HOME/waves2Foam
nikos@nikos:~/waves2Foam$ ./Allwmake

FATAL ERROR.
The directory path /home/nikos/OpenFOAM/nikos-2.2.2/applications/utilities/waves2Foam does not exist.
Correct the path in bin/bashrc

This means that you should correct the environmental variable $WAVES_DIR
such that it is pointed (in absolute sense) to the path, where waves2Foam
is located.

Please note that once bin/bashrc is created, bin/bashrc.org is only an
inactive file. The latter is also the only of the two files, which is
updated through the SVN-repository.

EXITING

nikos@nikos:~/waves2Foam$


Thank you in advance.

[ngj]

Hi Nick,

Please do read the error message carefully. I cannot say it more clearly myself.
If still in doubt, please address the posts above (e.g. use google), where the same question has been asked and resovled a number of times.

Kind regards,

Niels

[Nick_civ]

Hi Niels,

There was some confusion about directories. I finally compiled it succesfully.
Thank you for your time.

NK

[eternityboy]

hello, i would like to set potentialCurrent on both sides. from the inlet side and the outlet

but when I do so, i got strange freeSurface behavior after the relaxation zone

inlet2Coeffs
{
waveType potentialCurrent;
U (2 0 0);
Tsoft 0;
relaxationZone
{
relaxationScheme Spatial;
relaxationShape Rectangular;
beachType Empty;

relaxType INLET;
startX (-8 2.7 0.1868);
endX (-1.5 -2.7 0.1868);
orientation (1.0 0.0 0.0);
}
};

outletCoeffs
{
waveType potentialCurrent;
U (0 0 0);
Tsoft 0;

relaxationZone
{
relaxationScheme Spatial;
relaxationShape Rectangular;
beachType Empty;

relaxType OUTLET;
startX (8 2.7 0.1868);
endX (10 -2.7 0.1868);
orientation (1.0 0.0 0.0);
}
};

[ngj]

It is very difficult to help, when you are not stating what type of weird behaviour you see. But, I will make a guess that one problem could be that you impose 2.0 m/s from one side and have 0.0 m/s on the other.

Kind regards,

Niels

[Behrang]

Hi Neils,

In waveFlume tutorial, I modified the wave theory in waveProperties file, in order to simulate wave within the wave flume based on Stokes second and fifth order wave theory. By visualizing the results, I can see that waves are generated, but How can I support it? (I mean, is there any experimental data available to compare the results and validate them?)
Please advise.

Best,
Behrang

[ngj]

Good evening Behrang,

The literature is very rich on analytical expressions and experimental (field and laboratory) data, so I would recommend you to find data sets, which are closely related to your topic at hand.

Another source of information could be any of the articles, which have already been applying waves2Foam: http://openfoamwiki.net/index.php/Co...in_These_Works

Kind regards,

Niels

[kiang]

Hi Dear Niels,

Thanks a lot for your waveFoam solver. I am going to make a new waveTheory but I did not find a tutorial about it. Would you guide me how can I make a new wave based on existing wave theories in waves2Foam?

Regards,

Kiang

[ngj]

Hallo Kiang,

Have a look in the folder:

waves2Foam/doc/templateWaveTheory

If you execute the script "makeNewWaveTheory", it should create the necessary folder for you, so all you need to modify is the algebraic expressions. The script is heavily commented, so it should be possible to read through it step-by-step.

Example of execution:

./makeNewWaveTheory myNewWaveTheory

The output to screen is then:

A directory called 'myNewWaveTheory' with [C,H]-files is created.

Steps to finalise:
1. Copy the directory to the correct wave type in waves2Foam/src/waveTheories
2. Insert the needed variables and expressions for eta, U, etc.
3. Add a line in waves2Foam/src/Make/files to include the new theory in the compilation.
4. Done.

Kind regards,

Niels

[luigi21]

Hello everybody..

I'm trying to run a waveflume case in order to compare some results with experimental data, I've created a flume with (40mx0.5mx0.1m) (x,y,z), and the resolution in each direction is (0.02, 0.0025, 1), I've got the error below:

ecopetrol@ecopetrol-OptiPlex-755:~/OpenFOAM/OpenFOAM_simulations/coimbra_cluster_ext$ [0] [1] ##0 0 Foam::error:rintStack(Foam::Ostream&)Foam::error :rintStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[0] #2 at ??:?
[1] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) at ??:?
[1] #4 void Foam::MULES::limit<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double, int, bool) at ??:?
[0] #4 void Foam::MULES::limit<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double, int, bool) at ??:?
[1] #5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) at ??:?
[0] #5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) at ??:?
[1] #6 at ??:?
[0] #6

[1] at ??:?
[1] #7 __libc_start_main[0] at ??:?
[0] #7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #8 in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #8

[0] at ??:?
[ecopetrol-OptiPlex-755:07467] *** Process received signal ***
[ecopetrol-OptiPlex-755:07467] Signal: Floating point exception (8)
[ecopetrol-OptiPlex-755:07467] Signal code: (-6)
[ecopetrol-OptiPlex-755:07467] Failing at address: 0x3e800001d2b
[1] at ??:?
[ecopetrol-OptiPlex-755:07468] *** Process received signal ***
[ecopetrol-OptiPlex-755:07468] Signal: Floating point exception (8)
[ecopetrol-OptiPlex-755:07468] Signal code: (-6)
[ecopetrol-OptiPlex-755:07468] Failing at address: 0x3e800001d2c
[ecopetrol-OptiPlex-755:07467] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x370b0) [0x7f882b4fb0b0]
[ecopetrol-OptiPlex-755:07467] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7f882b4fb037]
[ecopetrol-OptiPlex-755:07467] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x370b0) [0x7f882b4fb0b0]
[ecopetrol-OptiPlex-755:07467] [ 3] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geo metricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEER KT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7vo lMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMes hEEESK_RKT0_RKT1_ddi+0xe3c) [0x7f882d66be1c]
[ecopetrol-OptiPlex-755:07467] [ 4] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES5limitINS_17geome tricOneFieldENS_9zeroFieldES3_EEvRKT_RKNS_14Geomet ricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS7_Id NS_13fvsPatchFieldENS_11surfaceMeshEEERSF_RKT0_RKT 1_ddib+0x244) [0x7f882d66fad4]
[ecopetrol-OptiPlex-755:07467] [ 5] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveER NS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMesh EEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEER S8_dd+0x79) [0x7f882d660ba9]
[ecopetrol-OptiPlex-755:07467] [ 6] waveFoam() [0x42fd9b]
[ecopetrol-OptiPlex-755:07467] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7f882b4e5ea5]
[ecopetrol-OptiPlex-755:07468] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x370b0) [0x7f37f121d0b0]
[ecopetrol-OptiPlex-755:07467] [ 8] waveFoam() [0x430f35]
[ecopetrol-OptiPlex-755:07467] *** End of error message ***
[ecopetrol-OptiPlex-755:07468] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7f37f121d037]
[ecopetrol-OptiPlex-755:07468] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x370b0) [0x7f37f121d0b0]
[ecopetrol-OptiPlex-755:07468] [ 3] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geo metricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEER KT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7vo lMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMes hEEESK_RKT0_RKT1_ddi+0xe6e) [0x7f37f338de4e]
[ecopetrol-OptiPlex-755:07468] [ 4] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES5limitINS_17geome tricOneFieldENS_9zeroFieldES3_EEvRKT_RKNS_14Geomet ricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS7_Id NS_13fvsPatchFieldENS_11surfaceMeshEEERSF_RKT0_RKT 1_ddib+0x244) [0x7f37f3391ad4]
[ecopetrol-OptiPlex-755:07468] [ 5] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveER NS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMesh EEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEER S8_dd+0x79) [0x7f37f3382ba9]
[ecopetrol-OptiPlex-755:07468] [ 6] waveFoam() [0x42fd9b]
[ecopetrol-OptiPlex-755:07468] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7f37f1207ea5]
[ecopetrol-OptiPlex-755:07468] [ 8] waveFoam() [0x430f35]
[ecopetrol-OptiPlex-755:07468] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 7467 on node ecopetrol-OptiPlex-755 exited on signal 8 (Floating point exception).



I attached the whole case looking forward someone could catch some error, I'd appreciate it so much, thank you guys.

coimbra_case.zip

[weezer579]

Helllo Niels,

Thank you very much for this great toolbox. It really saves me tons of time in my research.

I am trying to simulate a 2nd order Stokes wave (using waveFlume from tutorial) with the following properties:

stokesSecond
wave height: 0.3
wave number: 1.70048
omega: 3.1415
wave depth: 2
relaxationZone: default

The result looks very nice in the relaxationZone (Inlet). However, right after passing the this region, the wave crest seems to be pushed further to the right and making the wavelength longer, wave height shorter. I am not sure if this could be due to the my wrong set up with relaxationZone or something else. I have attached the surface elevation to this post. I'd really appreciate if you could take a quick look.

Thank you again.

[ngj]

Good morning Weezer,

Your specifications of water depth, cyclic frequency and wave number does not match the linear dispersion relation. Use the setWaveParameters utility to correct this issue.

Kind regards,

Niels

[weezer579]

Good morning Neils,

Thank you very much for your reply. I did try your suggestion and changed my omega value to 4.2129 to satisfy the dispersion relation. The result looks a lot better and matches the theoretical shape very well. However, the amplitude is still smaller than prediction. I've attached the results to this post. Would you mind taking a look again?

Thank you again for your time and help.

[aisy]

Hi Neils,
Can you tell me how to set the debug in waveProperties for stokesSecond ?
Is it right if I set debug as on ?
Thanks!