[Phicau]

Dear colleagues,


it is my pleasure to present a new numerical model based on the popular and widely used OpenFOAM: OLAFOAM.

OLAFOAM is the evolution of IHFOAM, a well-known three-dimensional Navier-Stokes solver for free surface flows, and is born with the spirit of providing the latest advances for the simulation of wave dynamics to the OpenFOAM and FOAM-extend communities. New implementations and updates will be readily available, bug fixes will be a priority and reasonable support will be given through the appropriate channels.

OLAFOAM provides cutting-edge technologies as state-of-the-art active wave generation and absorption techniques working at the boundaries (i.e. without increasing the computational cost), two phase flow through porous media or the newly introduced moving multi-paddle piston-type wavemakers, which also perform active wave absorption.


The release includes the complete source code, pre- and post-processing tools, 6 tutorials and additional reference materials. You can find all the details in the project web site:

https://sites.google.com/site/olafoamcfd/



Best regards,


Pablo Higuera

PhD in Civil Engineering, MSc

Research Fellow

National University of Singapore

[dneves]

If you keep with the same development dynamics as IHFOAM I will be certainly be a OLAFOAM user

[ashkan]

Hi Pablo,
Many thanks for your great work on IHFoam. I will download and start working wth your new olaFoam as well.

I found your wavemaker option very interesting and wanted to use that along with time series from wave tank. However, I was wondering if it is possible to model flap type wave makers as well or it only work with piston type wave makers?

Many thanks again
Ashkan

[Phicau]

Thanks Ashkan,

currently only piston type wavemakers are supported, but flap type is under development and I don't expect that it will take me long to polish it and release it.

I will keep you updated here: http://www.cfd-online.com/Forums/ope...am-thread.html

Best,

Pablo

[ashkan]

Originally Posted by Phicau

Thanks Ashkan,

currently only piston type wavemakers are supported, but flap type is under development and I don't expect that it will take me long to polish it and release it.

I will keep you updated here: http://www.cfd-online.com/Forums/ope...am-thread.html

Best,

Pablo

Thanks Pablo. Again thanks for your great work.
Ashkan

[smhashemi]

Hi Pablo,
First, I'd like to thank you for your admirable work and support.
Actually I'm going to work with a cluster that OpenFOAM 3.0.1 is installed on it.
Unfortunately for now there is not any IHFOAM or OLAFOAM solver for this version of OpenFOAM.
With all respect, I wanted to know if there is any hope for the solver of version 3.0.1 to be released in the near future?

Best Regards,

Mahdi

[Phicau]

Hi Mahdi,

try to compile the solver version 3.0.0. If I am not wrong, it should work. Otherwise, just let me know.

Best,

Pablo

[smhashemi]

Originally Posted by Phicau

Hi Mahdi,

try to compile the solver version 3.0.0. If I am not wrong, it should work. Otherwise, just let me know.

Best,

Pablo

Thank you Pablo,
I tried your opinion but the solver version 3.0.0 was not compiled on OF3.0.1.
I look forward to have your valuable comments.

Best regards,

Mahdi

[Phicau]

Hello again Mahdi,

there is a problem with your OpenFOAM installation, check/ask for a solution at the installation subforum. I have just tried and olaFoam compiles with no problems for OpenFOAM 3.0.1 (build 119cac7e8750).

Best,

Pablo

[JulienS]

Hi Pablo,

First, I would like to thank you for the work you did on IHFoam (which I was trained to in Santander, good memories!) and now OlaFoam.

I recently installed OlaFoam on my cluster, with OpenFOAM v3.0+. All compiles nicely. But, it seems that there is a problem with IOobject constructor. "NO_READ" and "READ_IF_PRESENT" seems not to work anymore.

-> Error with NO_READ
--> FOAM FATAL ERROR:
NO_READ specified for read-constructor of object porosity of class IOobject

From function Foam::Istream& Foam::regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 46.


-> Error with READ_IF_PRESENT
--> FOAM FATAL IO ERROR:
cannot find file

file: [..]/run/Tutorials/olafoam/baseWaveFlume/0/porosityIndex at line 0.

Have you already encountered this problem?
Best regards

Julien

[Phicau]

Salut Julien,

thanks for your words and your bug report.

Since I don't work with v3.0+ I had not noticed this error before, I have just pushed the fix to the repository, it should work now. It seems that the constructor is less permissive in this version and you need to specify the value and dimensions.

If you encounter any other flaws let me know.

Best regards,

Pablo

[JulienS]

Hi Pablo,

thank you for your quick response, it seems to work fine now.

Best regards

Julien

[theogrm]

Hi all,

I've been trying to simulate bichromatic wave groups propagating over a sloping beach using the irregular waveType. The problem I am having is the water level that is constantly increasing. I have in the waveDict both absGen 1 and absDir 0. Tried different combination with no success. Any hint will be appreciated.

Best Regards,

Theo

[kenneth0811]

Hi all

I am OLAFOAM user. I'd like to thank you for your admirable work and support. Through OLAFOAM i can learn a lot about the CFD analysis.
I used LES model in 2Dcase(RegularWave Stokes II). and I used blockMesh and SnappyHexMesh
but When I execute OLAFOAM, terminal show

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model kEqn
Selecting LES delta type smooth
Selecting LES delta type cubeRootVol
--> FOAM Warning :
From function cubeRootVolDelta::calcDelta()
in file LES/LESdeltas/cubeRooVolDelta/cubeRootVolDelta.C at line 55
Case is 2D, LES is not strictly applicable

Is it correct OLAFOAM exection? Although this FOAM Warning, OLAFOAM activated to the lasttime.

This is turbulenceProperties file , blockMeshDirt file in my case and ParaFoamView(using SnappyHexMesh refine).

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

vertices
(
(0 0 0)
(30 0 0)
(51.5 0 0)
(51.6 0 0)
(81.6 0 0)
(81.6 0 4)
(51.6 0 4)
(51.5 0 4)
(30 0 4)
(0 0 4)
(0 0 4.074)
(30 0 4.074)
(51.5 0 4.074)
(51.6 0 4.074)
(81.6 0 4.074)
(81.6 0 15.24089655)
(51.6 0 15.24089655)
(51.5 0 15.24089655)
(30 0 15.24089655)
(0 0 15.24089655)

(0 0.125 0)
(30 0.125 0)
(51.5 0.125 0)
(51.6 0.125 0)
(81.6 0.125 0)
(81.6 0.125 4)
(51.6 0.125 4)
(51.5 0.125 4)
(30 0.125 4)
(0 0.125 4)
(0 0.125 4.074)
(30 0.125 4.074)
(51.5 0.125 4.074)
(51.6 0.125 4.074)
(81.6 0.125 4.074)
(81.6 0.125 15.24089655)
(51.6 0.125 15.24089655)
(51.5 0.125 15.24089655)
(30 0.125 15.24089655)
(0 0.125 15.24089655)

);

blocks
(
hex (0 1 21 20 9 8 28 29) (240 1 32) simpleGrading (1 1 1)
hex (1 2 22 21 8 7 27 28) (172 1 32) simpleGrading (1 1 1)
hex (2 3 23 22 7 6 26 27) (1 1 32) simpleGrading (1 1 1)
hex (3 4 24 23 6 5 25 26) (240 1 32) simpleGrading (1 1 1)
hex (9 8 28 29 10 11 31 30) (240 1 1) simpleGrading (1 1 1)
hex (8 7 27 28 11 12 32 31) (172 1 1) simpleGrading (1 1 1)
hex (7 6 26 27 12 13 33 32) (1 1 1) simpleGrading (1 1 1)
hex (6 5 25 26 13 14 34 33) (240 1 1) simpleGrading (1 1 1)
hex (10 11 31 30 19 18 38 39) (240 1 88) simpleGrading (1 1 1)
hex (11 12 32 31 18 17 37 38) (172 1 88) simpleGrading (1 1 1)
hex (12 13 33 32 17 16 36 37) (1 1 88) simpleGrading (1 1 1)
hex (13 14 34 33 16 15 35 36) (240 1 88) simpleGrading (1 1 1)
// hex (24 25 37 36 31 30 42 43) (125 1 17) simpleGrading (1 1 1)
// hex (25 26 38 37 30 29 41 42) (80 1 17) simpleGrading (1 1 1)
// hex (26 27 39 38 29 28 40 41) (125 1 17) simpleGrading (1 1 1)
// hex (31 30 42 43 32 33 45 44) (125 1 31) simpleGrading (1 1 1)
// hex (30 29 41 42 33 34 46 45) (80 1 31) simpleGrading (1 1 1)
// hex (29 28 40 41 34 35 47 46) (125 1 31) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
inlet
{
type patch;
faces
(
(0 20 29 9)
(9 29 30 10)
(10 30 39 19)
);
}
outlet
{
type patch;
faces
(
(4 24 25 5)
(5 25 34 14)
(14 34 35 15)
);
}
bottom
{
type wall;
faces
(
(0 1 21 20)
(1 2 22 21)
(2 3 23 22)
(3 4 24 23)
);
}
atmosphere
{
type patch;
faces
(
(19 18 38 39)
(18 17 37 38)
(17 16 36 37)
(16 15 35 36)
);
}
);

mergePatchPairs
(
);

// ************************************************** *********************** //

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType LES;

LES
{
LESModel kEqn;

turbulence on;

printCoeffs on;

delta smooth;

cubeRootVolCoeffs
{
deltaCoeff 1;
}

PrandtlCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

maxDeltaRatio 1.1;
}

Cdelta 0.158;
}

vanDriestCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

maxDeltaRatio 1.1;
}

Aplus 26;
Cdelta 0.158;
}

smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

maxDeltaRatio 1.1;
}
}

// ************************************************** *********************** //

// ************************************************** *********************** //




I do not English well. My post is quite a chaos. Sorry.

[Phicau]

Hi Ju Hyun.

this is not exactly the best place for support, there is a dedicated post for that:

http://www.cfd-online.com/Forums/ope...am-thread.html

You are applying LES in 2D, which given its assumptions is not strictly applicable. For LES to be fully correct it needs to be applied in 3D.

Best,

Pablo

[kenneth0811]

Originally Posted by Phicau

Hi Ju Hyun.

this is not exactly the best place for support, there is a dedicated post for that:

http://www.cfd-online.com/Forums/ope...am-thread.html

You are applying LES in 2D, which given its assumptions is not strictly applicable. For LES to be fully correct it needs to be applied in 3D.

Best,

Pablo

Hi Pablo.

I didn't know 'The OLAFOAM Thread'. I communicate by 'The OLAFOAM Thread' now.
and It is grateful that you gave an answer to my question.

Ju Hyun

[Ahmed Elhanafi]

Hi Pablo,

First of all, thank you for your support to OpenFOAM users. Could you please let me know what are the main differences between OLAFOAM and IHFOAM?

Kind regards,
Ahmed

[Phicau]

Hi Ahmed,

thanks.

Some facts (as of 31/07/2017):

https://github.com/phicau/IHFOAM
First commit: Jul 8, 2014
Last commit: Feb 11, 2016
Number of commits: 21
Last solver included: OpenFOAM 3.0.0
Watch: 4 Star: 2 Fork: 6

https://github.com/IHFOAM/ihFOAM
First commit: Sept 30, 2015
Last commit: Jan 18, 2017
Number of commits: 8
Last solver included: OpenFOAM v1612+
Watch: 1 Star: 2 Fork: 3
Note: My name has been deleted from the code and has been replaced by Javier Lopez Lara and Gabriel Barajas

https://github.com/IHCantabria/ihFOAM
First commit: Sept 30, 2015
Last commit: Mar 21, 2017
Number of commits: 10
Last solver included: OpenFOAM v1612+
Watch: 2 Star: 0 Fork: 0
Note: My name has been deleted from the code and has been replaced by Javier Lopez Lara and Gabriel Barajas
Note: Exactly the same as https://github.com/IHFOAM/ihFOAM, but with a new logo

For the record, they have integrated ihFoam in OpenFOAM+, with a fancier programming style, but they seem to have removed some wave theories and added others.

https://github.com/phicau/OLAFOAM
First commit: Mar 2, 2016
Last commit: Jul 3, 2017
Number of commits: 21
Last solver included: OpenFOAM v1706+
Watch: 10 Star: 11 Fork: 5
Note: The first commit is exactly the same as the last commit of https://github.com/phicau/IHFOAM

Since I started developing olaFoam I have fixed bugs, updated the compatibility for new versions, added wave theories and the moving boundary wave generation module.

Best,

Pablo

[Ahmed Elhanafi]

Hi Pablo,

Thank you for the detailed explanation. It is very much appreciated.

I have another subject to discuss. I have recently attended the ISOPE conference. During this event, a participant has mentioned in his paper that when using the inlet boundary conditions of IHFOAM a damping in wave height (wave height reduction) was observed when a turbulence model such as k-epsilon was implemented compared to a laminar model. I have observed the same damping issue with the k-epsilon turbulence model, but I found no wave damping when utilising SST k-w model. However, a recent paper published in Coastal Engineering Journal "Application of a buoyancy-modified k-ω SST turbulence model to simulate wave run-up around a monopile subjected to regular waves using OpenFOAM", found that there is a damping in wave height when using SST k-w model, and therefore the authors provided a solution of this by including a buoyancy term in the turbulent kinetic energy equation. I have emailed the authors of this paper asking about the turbulence model settings they have used to report such damping events, but they have not responded yet. I'm wondering if you have experienced this issue before with either IHFOAM or OLAFOAM.

Kind regards,
Ahmed

[Phicau]

Hi Ahmed,

it is well known that turbulence modelling does not work perfectly with transient two-phase flows in OpenFOAM. It was developed for other purposes and, therefore, you will sometimes produce a turbulence buildup that will damp the waves. The problem is not related with olaFoam or ihFoam, and it has been widely discussed in the forum, just make a quick search for it.

I suggest that if you have further questions, you can ask them in the support thread, as this one is quite hidden because of the dig-policy of this section of the forum.

https://www.cfd-online.com/Forums/op...am-thread.html

Best,

Pablo