[wavemaster]

Dear All,
I am probably missing some obvious information. What is the easiest way to generate a given timetrace of surface elevation in IHFOam? I tried to use the wavemaker, but that seems to require paddel information which is not what I want.
Any comments are much appreciated and thanks for making this interesting tool available!
Cheers,
Pal

[rafa13]

Originally Posted by Phicau

Hi Haze,

thanks for your code that extracts free surface as a 3D surface along the whole domain. It has been around other threads, but it can be handy for users, as sometimes things are difficult to look for when you don't know that they exist.

There is a way to run this function object after the case has finished, to work like the sample tool. You just need to run execFlowFunctionObjects.

You are right, alpha = 0.5 is just an arbitrary convention, but that cell is really half air and half water, so there is plenty of water above it yet. In our papers we found out that to estimate run-up it is reasonable to consider alpha = 0.9. Moreover, there are lots of ways to obtain free surface elevation, we cover 4 or 5 in the training courses.

Did you already register in http://ihfoam.ihcantabria.com/source-download/ ? We do provide tools to obtain free surface elevation at any location in the free materials that you can download when you register. Check the breakwater case. There you can find a procedure to sample the pressure and free surface elevation, and some python tools to convert the raw data into a time series. You also have a python tool to plot it.

Best,

Pablo

Hi Pablo,


I am want to ask if you can indicate, in which Paper you about the tuning of the indicator alpha? I want to know a little bit more about this calibration.

I wrote 3 of your paper but the are nothing about the tuning.


thanks a lot,

Greets Rafael Marques

[Phicau]

Hi Pal,

probably the best way is to decompose your signal into components by means of an FFT and set up the irregular waves.

Hi Rafa,

there is no calibration/tuning for that, check "Simulating coastal engineering processes with OpenFOAM": An alpha1 value of 0.1 has been used to consider
a cell as inundated. This 10% threshold translates in an effective
water height on those cells of 1.25 mm, which is a convenient value
regarding the resolution of experimental devices.

Best,

Pablo

[rafa13]

Hi Pablo,

Thanks for the quick response. I will check this paper.


Thanks again you are helping me a lot!

Greets Rafael Marques

[HosseinB]

Hi Pablo,

I have received training in ihFoam in May and have a series of questions and would be grateful if you could help me with it.

I am modelling stokesV waves and they are close to target at the wave-generating boundary, however, when waves travel inside the domain, their trough and crest locations change and make it challenging to correctly predict resultants that depend on the water column elevation. My wave height is 0.512 m and cubic cell size is 0.028 m. I appreciate it if you could help me with this.

Cheers,
Hossein

[jdttl]

Hi, I currently downloaded the IHFoam extension. During compilation It has been a success boundary conditions:

cd genAbs
./allMake

However, I encountered an error for the second compilation, which is for the solvers:

cd genAbs
./allMake

Here is the error on in my terminal:

user@user:~/ihFOAM/solvers/ihFoamFE30_31$ ./allMake
Making dependency list for source file ihFoam.C
could not open file cyclicAMILduInterface.H for source file ihFoam.C due to No such file or directory
could not open file cyclicAMIPolyPatch.H for source file ihFoam.C due to No such file or directory
could not open file twoPhaseMixture.H for source file ihFoam.C due to No such file or directory
could not open file turbulenceModel.H for source file ihFoam.C due to No such file or directory
could not open file readPIMPLEControls.H for source file ihFoam.C due to No such file or directory
g++ -m32 -Dlinux -DWM_ARCH_OPTION=32 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam30/src/transportModels -I/opt/openfoam30/src/transportModels/incompressible/lnInclude -I/opt/openfoam30/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam30/src/turbulenceModels/incompressible/turbulenceModel -I/opt/openfoam30/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam30/src/OpenFOAM/lnInclude -I/opt/openfoam30/src/OSspecific/POSIX/lnInclude   -fPIC -c ihFoam.C -o Make/linuxGccDPInt32Opt/ihFoam.o
In file included from /opt/openfoam30/src/finiteVolume/lnInclude/ddtScheme.C:30:0,
                 from /opt/openfoam30/src/finiteVolume/lnInclude/ddtScheme.H:337,
                 from /opt/openfoam30/src/finiteVolume/lnInclude/fvcDdt.C:28,
                 from /opt/openfoam30/src/finiteVolume/lnInclude/fvcDdt.H:199,
                 from /opt/openfoam30/src/finiteVolume/lnInclude/fvc.H:44,
                 from /opt/openfoam30/src/finiteVolume/lnInclude/fvCFD.H:8,
                 from ihFoam.C:64:
/opt/openfoam30/src/finiteVolume/lnInclude/cyclicAMIFvPatch.H:39:35: fatal error: cyclicAMILduInterface.H: No such file or directory
 #include "cyclicAMILduInterface.H"
                                   ^
compilation terminated.
make: *** [Make/linuxGccDPInt32Opt/ihFoam.o] Error 1
ihFoam compilation failed

Is there something I should have done first before running the second compilation? I just followed the instructions as stated in the installation process. Any reply would be greatly appreciated. Thank you very much.

[Phicau]

Hi Joseph,

you ere trying to compile Foam-extend version with OpenFOAM 3.0. I haven't had the time to prepare an OpenFOAM 3.0 solver yet. I guess you'll need to wait until I find some time or create it yourself.

Best,

Pablo

[jdttl]

Originally Posted by jdttl

Hi, I currently downloaded the IHFoam extension. During compilation It has been a success boundary conditions:

cd genAbs
./allMake

However, I encountered an error for the second compilation, which is for the solvers:

cd genAbs
./allMake

I have tried compilation of ihFOAM in OpenFOAM 2.4.0, both compilation went successful. The problem i had posted above was encountered in OpenFOAM 3.0.0. It seems that the compilation of solvers of ihFOAM for OpenFOAM 3.0.0 still have a minimal problem.

[jdttl]

Originally Posted by Phicau

Hi Joseph,

you ere trying to compile Foam-extend version with OpenFOAM 3.0. I haven't had the time to prepare an OpenFOAM 3.0 solver yet. I guess you'll need to wait until I find some time or create it yourself.

Best,

Pablo

Hi Pablo,

Thank you very much for your time. For the meantime, I will be using OpenFoam 2.4.0 since everything works fine for this version.
By the way, I am new to OpenFoam so as to ihFoam, forgive me if I might ask some noob questions. Like I am trying now to try the tutorials, breakwater, specifically. Am I right to use "blockMesh" to create the Mesh? It gave no errors, however, after as I input to the terminal "ihFoam" for simulation, I encountered this error:

--> FOAM FATAL IO ERROR: 
cannot find file

file: /home/user/OpenFOAM/user-2.4.0/tutorials/OF230_240/breakwater/0/p_rgh at line 0.

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

FOAM exiting

[jdttl]

Hi Pablo,

How are you? I managed to run now the blockMesh and ihFoam commands. I am using OpenFoam 2.4. However, I notice some questionable error that prompted in my terminal after I run the "ihFoam" :

user@user:~/OpenFOAM/user-2.4.0/tutorials/OF230_240/breakwater$ ihFoam 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.4.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.4.0-f0842aea0e77
Exec   : ihFoam
Date   : Nov 22 2015
Time   : 13:41:18
Host   : "user"
PID    : 31674
Case   : /home/user/OpenFOAM/user-2.4.0/tutorials/OF230_240/breakwater
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading field porosityIndex

Porosity NOT activated

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
--> Upgrading k to employ run-time selectable wall functions
    Backup original k to k.old
    Writing updated k
--> Upgrading epsilon to employ run-time selectable wall functions
    Backup original epsilon to epsilon.old
    Writing updated epsilon
--> Creating nut to employ run-time selectable wall functions
    Writing new nut
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    sigmaEps        1.3;
}


Reading g
Calculating field g.h

No finite volume options present

time step continuity errors : sum local = 0, global = 0, cumulative = 0
DICPCG:  Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
deltaT = 0.00119048
Time = 0.00119048

PIMPLE: iteration 1
alpha.water BC on patch inlet
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/opt/openfoam240/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#3  cnoidalFun::calculations(double, double, double, double*, double*) at ??:?
#4  Foam::IH_Waves_InletAlphaFvPatchScalarField::updateCoeffs() at ??:?
#5  Foam::fvPatchField<double>::evaluate(Foam::UPstream::commsTypes) at ??:?
#6  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() at ??:?
#7  ? at ??:?
#8  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#9  ? at ??:?
Floating point exception (core dumped)

I haven't tried editing the corresponding files inside breakWater directory though.

I then tried to run ParaView to see the simulation:

user@user:~/OpenFOAM/user-2.4.0/tutorials/OF230_240/breakwater$ paraFoam
created temporary 'breakwater.OpenFOAM'

ParaView opens but when I clicked Apply, and tries to simulate visually, nothing happens.

Is there something wrong with my process?

Regards,

Joseph D.

[Phicau]

Hi Joseph,

have you taken a look at the wiki? It's a nice place to start:

https://openfoamwiki.net/index.php/Contrib/IHFOAM

Each case has an runCase with all the steps needed to run the case, take a look at them.

Best,

Pablo

[Phicau]

Dear all,

I have updated IHFOAM to include the solvers and tutorials for the latest release: OpenFOAM 3.0.

https://github.com/phicau/IHFOAM

Best,

Pablo

[jdttl]

Hi Pablo,

Thank you very much for your work and time developing ihFoam. I had tried compiling the new solver that you've added for OpenFoam lastest release 3.0. However, it gave an error :

user@user:~/IHFOAM-master/solvers/ihFoamOF300$ ./allMake
Making dependency list for source file ihFoam.C
could not open file turbulentTransportModel.H for source file ihFoam.C due to No such file or directory
could not open file CorrectPhi.H for source file ihFoam.C due to No such file or directory
could not open file createTimeControls.H for source file ihFoam.C due to No such file or directory
could not open file createRDeltaT.H for source file ihFoam.C due to No such file or directory
could not open file readhRef.H for source file ihFoam.C due to No such file or directory
could not open file gh.H for source file ihFoam.C due to No such file or directory
could not open file createMRF.H for source file ihFoam.C due to No such file or directory
SOURCE=ihFoam.C ;  g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam240/src/transportModels/twoPhaseMixture/lnInclude -I/opt/openfoam240/src/transportModels -I/opt/openfoam240/src/transportModels/incompressible/lnInclude -I/opt/openfoam240/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam240/src/TurbulenceModels/turbulenceModels/lnInclude -I/opt/openfoam240/src/TurbulenceModels/incompressible/lnInclude -I/opt/openfoam240/src/transportModels/immiscibleIncompressibleTwoPhaseMixture/lnInclude -I/opt/openfoam240/src/finiteVolume/lnInclude -I/opt/openfoam240/src/fvOptions/lnInclude -I/opt/openfoam240/src/meshTools/lnInclude -I/opt/openfoam240/src/sampling/lnInclude -IlnInclude -I. -I/opt/openfoam240/src/OpenFOAM/lnInclude -I/opt/openfoam240/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linuxGccDPOpt/ihFoam.o
ihFoam.C:69:37: fatal error: turbulentTransportModel.H: No such file or directory
 #include "turbulentTransportModel.H"
                                     ^
compilation terminated.
make: *** [Make/linuxGccDPOpt/ihFoam.o] Error 1
ihFoam compilation failed

[Phicau]

Hi Joseph,

please be more careful and double check always, you are trying to compile the version 3.0 solver with version 2.4.

Best,

Pablo

[jdttl]

Hi Pablo,

Thank you very much for your reply. As soon as you posted the update for OF 3.0, i uninstalled my OF 2.4 and replaced it with 3.0. Then I tried the updated solver.

I tried to simulate the breakwater tutorial but got the same error that I encountered using OF 2.4.


Regards,

Joseph D.

[jdttl]

Hi Pablo,

Good Day! Thank you very much Pablo for your work. I reinstalled OpenFoam 3.0.0 along with Paraview 4.4. I completed the compilation for boundary conditions and solvers. Afterwards, I tried the tutorial inside OF300 directory, by running:

user@user:~/Downloads/IHFOAM/tutorials/OF300/baseWaveFlume$ ./runCase

The simulation was successful. And opened the baseWaveFlume.foam through Paraview, and played the simulation. I have another question, what is the unit of p_rgh? is it in Pascal? And I have an idea of adding some wind factor, is it possible to add wind velocity parameter?

Regards,

Joseph D.

[nctiwary]

hi Dmitrjs;
i want ur case setup for your above problem....i think i have the same kind of problem.the link u have attached is not working.

[haze_1986]

Hi, does anyone knows how do we decide on the wavePhase? When is it 3pi/2? I am getting problems with the boundary not being entirely correct.

[Gizela Zucula]

Hello, everyone. I'm new to IHFoam. I wanna know what are the unities for p_rgh? After running IhFoam, when it comes to read the results we have p and p_rgh, what's the difference between them? Is there anyway of calculating hydrostatic forces on a surface through using some tools in paraview or IHfoam solver? Thanks for your attention. Im waiting for any answer. Regards Gizela Zucula

[nctiwary]

Greatings Pablo,

I am a new user of the IHFOAM 2.0.

I am very impressed and pleased with the way the model behaves, I'm exploring it intensively.
i am also facing the same problem what u were facing for evaluating the wave-induced forces on a horizontal cylinder.(in "The IHFOAM Thread" post no.30).Can i know how u solve this problem.....can u send me the correct case file setup???
further do u have any case setup for vertical cylinder???....it will be very helpful for my experimental work validation.
do u have any idea about how to couple my "fully non-linear potential theory" code which is generating a list of velocity ordinates , with IHFoam???
thanx in advance
nirbhay
Capture.JPG