[jpbettenc]

Hi all,

I tried to compile ihFOAM with OF3.0.1 but it gave me an error due to some missing header files:

[jbettencourt@dhcp-892b8198 ihFOAM]$ of301
ModuleCmd_Load.c(213):ERROR:105: Unable to locate a modulefile for 'openmpi-x86_64'
[jbettencourt@dhcp-892b8198 ihFOAM]$ cd solvers/
[jbettencourt@dhcp-892b8198 solvers]$ ls
ihFoamFE30_31  ihFoamOF16ext  ihFoamOF171  ihFoamOF211  ihFoamOF222  ihFoamOF230  ihFoamOF240
[jbettencourt@dhcp-892b8198 solvers]$ cd ihFoam240
bash: cd: ihFoam240: No such file or directory
[jbettencourt@dhcp-892b8198 solvers]$ cd ihFoamOF240
[jbettencourt@dhcp-892b8198 ihFoamOF240]$ ls
allMake           alphaEqn.H          correctPhi.H    createPorosity.H  ihFoam.C  pEqn.H       UEqn.H
alphaCourantNo.H  alphaEqnSubCycle.H  createFields.H  ihDyMFoam         Make      setDeltaT.H
[jbettencourt@dhcp-892b8198 ihFoamOF240]$ ./allMake
Making dependency list for source file ihFoam.C
could not open file turbulenceModel.H for source file ihFoam.C due to No such file or directory
could not open file createPrghCorrTypes.H for source file ihFoam.C due to No such file or directory
g++ -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=64 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -O3  -DNoRepository -ftemplate-depth-100 -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/transportModels/twoPhaseMixture/lnInclude -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/transportModels -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/transportModels/incompressible/lnInclude -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/transportModels/interfaceProperties/lnInclude -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/turbulenceModels/incompressible/turbulenceModel -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/transportModels/immiscibleIncompressibleTwoPhaseMixture/lnInclude -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/finiteVolume/lnInclude -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/fvOptions/lnInclude -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/meshTools/lnInclude -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/sampling/lnInclude -IlnInclude -I. -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/OpenFOAM/lnInclude -I/home/jbettencourt/OpenFOAM/OpenFOAM-3.0.1/src/OSspecific/POSIX/lnInclude   -fPIC -c ihFoam.C -o Make/linux64GccDPInt64Opt/ihFoam.o
ihFoam.C:71:29: fatal error: turbulenceModel.H: No such file or directory
 #include "turbulenceModel.H"
                             ^
compilation terminated.
make: *** [Make/linux64GccDPInt64Opt/ihFoam.o] Error 1
ihFoam compilation failed

I know that OF 3.0.1 is not supported by ihFoam, but several people seem to have been able to compile it for that version. However, I have been having this problem.

Do you know of any solution to this error?

Thanks in advance,
Joao

[nctiwary]

Originally Posted by nctiwary

hi Pablo,
I am using IHFOAM to calculate the force on a verical cylinder under waves.I have imported my msh file from ansys fluent.Everything works fine but when i am trying to run the solver(IHFOAM) after some iteration it shows a Floating point exception (core dumped) error.Can u pls help me in finding the problem in my case set up.I am attaching my case set with this.

https://www.dropbox.com/s/z5gdxy94zn...er.tar.gz?dl=0

Hi Pablo
thank u for ur quick reply.Can u tell me what else i have to do in my mesh??Do i have to refine it more??

[Phicau]

Hi jpbettenc,

I am not sure whether 3.0.1 is supported (3.0.0 is, though), but compiling the version 2.4.0 solver with OpenFOAM 3.0.1 may be the problem. Please be more careful in the future.

Hi Nirbhay,

the resolution seems not adequate, it uses tetrahedra (not a problem per se, but not good either), the mesh is very narrow, the initial setFields yields a bumpy free surface...

[nctiwary]

Originally Posted by Phicau

Hi jpbettenc,

I am not sure whether 3.0.1 is supported (3.0.0 is, though), but compiling the version 2.4.0 solver with OpenFOAM 3.0.1 may be the problem. Please be more careful in the future.

Hi Nirbhay,

the resolution seems not adequate, it uses tetrahedra (not a problem per se, but not good either), the mesh is very narrow, the initial setFields yields a bumpy free surface...

hi pablo
thanks for ur suggestion.Should i refine mesh around the free surface level??What should i use instead of tetrahedra??

[jpbettenc]

Originally Posted by Phicau

Hi jpbettenc,

I am not sure whether 3.0.1 is supported (3.0.0 is, though), but compiling the version 2.4.0 solver with OpenFOAM 3.0.1 may be the problem. Please be more careful in the future.

Of course. I was not sure if it would work, but in the 'solvers' folder the highest solver version is indeed 2.4.0.

[jbettencourt@dhcp-892b8198 solvers]$ ls
ihFoamFE30_31  ihFoamOF171  ihFoamOF222  ihFoamOF240
ihFoamOF16ext  ihFoamOF211  ihFoamOF230

Which version of the solver should I try then?

Thanks,
Joao

[Phicau]

Hi Joao,

that is because you cloned it from the wrong place, the guys at IH Cantabria changed all the links in the wiki and IHFOAM site when I left the place, but they are not updating the code. Take a look here for the latest version, compatible with 3.0.0:

https://github.com/phicau/IHFOAM

Best,

Pablo

[jpbettenc]

Originally Posted by Phicau

Hi Joao,

that is because you cloned it from the wrong place, the guys at IH Cantabria changed all the links in the wiki and IHFOAM site when I left the place, but they are not updating the code. Take a look here for the latest version, compatible with 3.0.0:

https://github.com/phicau/IHFOAM

Best,

Pablo

Hi Pablo,

Thanks for clarifying this. I cloned the lastest version from https://github.com/phicau/IHFOAM and it compiled without problems with OF 3.0.1.

And thanks also for making this code available for everybody.

Best,
Joao

[MrChows]

Hi, I was trying to modify the code to allow custom user inputs for eta and the velocity vector at every patch. I have spent some time looking through the code. As far as I know paddleEta_ does not allow us to specify different eta for different y positions or paddles. Thus I have taken some useful code from multiPistonMovement.C.

Eventually, what I have in my calculatedLevelEtaCustom.H:

... CODE ...

    forAll( timeSeries_, contador1 )
    {
        if ( timeSeries_[contador1] >= currTime )
        {
            indexF = contador1;
            break;
        }
    }

    for(int i=0; i<nPaddles_; i++)
    {
        if ( indexF == 0)
        {
            etaCustomInterp[i] = timeMult * etaCustom_[i][0];
        }
        else
        {
            etaCustomInterp[i] = timeMult * 
	    interpolation(timeSeries_[indexF-1], timeSeries_[indexF], etaCustom_[i][indexF-1], etaCustom_[i][indexF], currTime);
        }
    
        calculatedLevel[i] = waterDepth_ + timeMult * etaCustomInterp[i];
    }

I am getting an error which I have no idea about:

*** Error in `ihFoam': free(): invalid next size (fast): 0x00000000218f01a0 ***
Aborted (core dumped)

Note that this error occurs before IH_Waves_InletVelocityFvPatchVectorField, knowing what went wrong here I should be able to fix the rest.

Can anyone help?

[dneves]

Hi Pablo,

What is exactly the Pseudo-Dynamic Pressure?

Thank you for your attention,

Best regards,

Diogo

[Phicau]

Hi Diogo,

take a look at my thesis, it is explained there.

http://www.tesisenred.net/bitstream/...pdf?sequence=1

Best,

Pablo

[dneves]

Thank you very much Pablo,

This type of questions is not allways well received and you promptly answered mine

My best regards,

Diogo

[Phicau]

Dear all,

I would like to make you all aware that yesterday I published OLAFOAM, the evolution of IHFOAM:

http://www.cfd-online.com/Forums/ope...on-ihfoam.html

This has some implications. All my efforts will be put in developing, maintaining and promoting the new code. This does not mean that I am going to drop support for IHFOAM here in this thread and at this stage (eventually I may), but it will certainly have a lower priority.

I will soon open an OLAFOAM support thread.

Best,

Pablo

[kasra karimi]

Hi Pablo
Thanks a lot for development of IHFoam. I am using ihFoam to generate waves in a simple 2D flume (6m length*0.4m hight* 0.1m width) and The waves parameter are (H=0.065 T=1.53 d=0.3) which makes it follow the StokesII wave theory.
The problem is that the surface elevation-time plot gives a good agreement with theoretical answers, but when it comes to horizontal orbital velocity, the simulation data have unacceptable underestimation in compared to theoretical.

I utilized both quadratic and triangular mesh in two different cases but my problem still remained the same.

heres my case files including some plots of horizontal orbital velocity (Ux) Vs time, for the elevation of -0.1 below surface & also UxMax profile-depth from Matlab code and OpenFOAM.
I would be so thankful if you could help me.
http://s000.tinyupload.com/?file_id=...16525276522753

[mo_na]

Hey guys,
I need to calculate the lift and drag coefficients for a pitching floater that is hinged at one point and is moving due to incoming waves.
What I am wondering about is the parameter magUInf that needs to be specified in the controlDict. As I have "only" a wave and not really a free stream velocity I am not sure what value to put there.
Any ideas?
Cheers,
Mona

[Jung hoo]

Hello pablo! I'm using IHFOAM well for my academic purpose and I'm really appreciative of your effort.

I'm simulating regular waves passing though simple vertical breakwater by using IHFOAM.

I want to get the overtopping flow rate by time above the breakwater and I think that it is possible through paraview using 'plot over line', 'integrating values' and 'plot selection over time' but I'm not sure it is working...

I read your paper 'Three-dimensional interaction of waves and porous coastal structures using OpenFOAM.Part 2 : Application' and I will be very thankful to you if you give me the way you plot 'Fig. 14. Space-time series of instantaneous overtopping discharge rate'.

Thank you so much in advance.. your any comment will very helpful to me..!!

[Phicau]

HI Jung Hoo,

that graph was not performed in paraview. I used OpenFOAM sampling tools to get the flow rate over the crown wall and then postprocessed them with matlab to obtain the final plot.

Best,

Pablo

[Jung hoo]

Thank you for your reply Pablo!

I have used 'sampleDict' recently, and now I can calculate some values like pressure and water surface level by time using 'sampleDict' and 'postSensPres.py', 'postSensVOF.py'.

but I can't understand that overtopping flow rate can be calculated by using 'sampleDict'... it must be due to my poor knowledge.. ....

Do you use 'alpha.water' as filedname in 'sampleDict' ??

If so, calculating overtopping flow discharge above a breakwater through 'sampleDict' seems to be so complicated to me..

If it doesn't hurt your feelings, could you show me the contents of 'sampleDict' that can be used for calculating overtopping flow discharge??

If the answer of my question is so complicated in writing, please just post any references about it!

And I'm very appreciative of your previous kind reply about 'p_rgh'!.. it helped me to understand 'pseudo-dynamic pressure'.

Thank you in advance and I'm always thankful for your effort you made until now in openfoam!!

[Phicau]

Probably you will be better off cutting with a plane in paraview and exporting alpha and velocity. There is plenty of information on how to do that in the forum. Then you can obtain overtopping with your favourite programming tool.

Best,

Pablo

[Jung hoo]

I'll try it using paraview! Thanks a lot!!

[Jung hoo]

Hi, I'm simulating regular waves passing through a breakwater using IHFOAM.

I applied absorbing condition on the outlet, but not on the inlet(only wave generating condition).

In this case, How can I absorb reflected waves induced by breakwater toward inlet?

When I was using waves2foam, I could do that by setting absorbing condition in front of the inlet.

Any comments can help me a lot. Thanks!