CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[IHFOAM] The IHFOAM Thread

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree57Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 17, 2016, 06:18
Default ihFOAM fails to compile with OpenFOAM 3.0.1
  #201
New Member
 
Join Date: Jan 2016
Posts: 3
Rep Power: 10
jpbettenc is on a distinguished road
Hi all,

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

Code:
[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
jpbettenc is offline   Reply With Quote

Old   February 17, 2016, 06:59
Default
  #202
New Member
 
nirbhay tiwary
Join Date: Aug 2015
Posts: 24
Rep Power: 10
nctiwary is on a distinguished road
Quote:
Originally Posted by nctiwary View Post
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??
nctiwary is offline   Reply With Quote

Old   February 17, 2016, 07:08
Default
  #203
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
Phicau is on a distinguished road
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...
Phicau is offline   Reply With Quote

Old   February 17, 2016, 07:28
Default
  #204
New Member
 
nirbhay tiwary
Join Date: Aug 2015
Posts: 24
Rep Power: 10
nctiwary is on a distinguished road
Quote:
Originally Posted by Phicau View Post
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??
nctiwary is offline   Reply With Quote

Old   February 17, 2016, 08:23
Default
  #205
New Member
 
Join Date: Jan 2016
Posts: 3
Rep Power: 10
jpbettenc is on a distinguished road
Quote:
Originally Posted by Phicau View Post
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.

Code:
[jbettencourt@dhcp-892b8198 solvers]$ ls
ihFoamFE30_31  ihFoamOF171  ihFoamOF222  ihFoamOF240
ihFoamOF16ext  ihFoamOF211  ihFoamOF230
Which version of the solver should I try then?

Thanks,
Joao
jpbettenc is offline   Reply With Quote

Old   February 17, 2016, 21:22
Default
  #206
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
Phicau is on a distinguished road
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
Phicau is offline   Reply With Quote

Old   February 18, 2016, 03:07
Default
  #207
New Member
 
Join Date: Jan 2016
Posts: 3
Rep Power: 10
jpbettenc is on a distinguished road
Quote:
Originally Posted by Phicau View Post
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
jpbettenc is offline   Reply With Quote

Old   February 22, 2016, 03:05
Default IHFoam custom eta at different paddles at wave generation
  #208
New Member
 
Chow
Join Date: Feb 2016
Posts: 1
Rep Power: 0
MrChows is on a distinguished road
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:
... 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:
Code:
*** 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?
MrChows is offline   Reply With Quote

Old   March 1, 2016, 09:48
Default
  #209
New Member
 
Diogo R.C.B. Neves
Join Date: Aug 2014
Posts: 16
Rep Power: 11
dneves is on a distinguished road
Hi Pablo,

What is exactly the Pseudo-Dynamic Pressure?

Thank you for your attention,

Best regards,

Diogo
dneves is offline   Reply With Quote

Old   March 1, 2016, 20:46
Default
  #210
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
Phicau is on a distinguished road
Hi Diogo,

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

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

Best,

Pablo
Phicau is offline   Reply With Quote

Old   March 2, 2016, 04:12
Default
  #211
New Member
 
Diogo R.C.B. Neves
Join Date: Aug 2014
Posts: 16
Rep Power: 11
dneves is on a distinguished road
Thank you very much Pablo,

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

My best regards,

Diogo
dneves is offline   Reply With Quote

Old   March 3, 2016, 21:01
Default
  #212
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
Phicau is on a distinguished road
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
Phicau is offline   Reply With Quote

Old   April 11, 2016, 15:31
Default problem with StokesII-orbital velocity
  #213
New Member
 
kasra karimi
Join Date: Jan 2016
Location: Tehran
Posts: 10
Rep Power: 10
kasra karimi is on a distinguished road
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
kasra karimi is offline   Reply With Quote

Old   May 25, 2016, 10:24
Default
  #214
Member
 
mo_na's Avatar
 
Mona
Join Date: Mar 2016
Location: Berlin
Posts: 49
Rep Power: 10
mo_na is on a distinguished road
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
mo_na is offline   Reply With Quote

Old   May 30, 2016, 00:59
Default Calculating overtopping flow rate
  #215
Member
 
Lee Jung Hoo
Join Date: Dec 2015
Posts: 37
Rep Power: 10
Jung hoo is on a distinguished road
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..!!
Jung hoo is offline   Reply With Quote

Old   June 1, 2016, 04:05
Default
  #216
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
Phicau is on a distinguished road
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
Phicau is offline   Reply With Quote

Old   June 1, 2016, 08:29
Default
  #217
Member
 
Lee Jung Hoo
Join Date: Dec 2015
Posts: 37
Rep Power: 10
Jung hoo is on a distinguished road
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!!
Jung hoo is offline   Reply With Quote

Old   June 1, 2016, 22:17
Default
  #218
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
Phicau is on a distinguished road
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
Phicau is offline   Reply With Quote

Old   June 2, 2016, 00:38
Default
  #219
Member
 
Lee Jung Hoo
Join Date: Dec 2015
Posts: 37
Rep Power: 10
Jung hoo is on a distinguished road
I'll try it using paraview! Thanks a lot!!
Jung hoo is offline   Reply With Quote

Old   June 5, 2016, 03:29
Default Absorption at inlet
  #220
Member
 
Lee Jung Hoo
Join Date: Dec 2015
Posts: 37
Rep Power: 10
Jung hoo is on a distinguished road
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!
Jung hoo is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Divergence detected in AMG solver: k when udf loaded google9002 Fluent UDF and Scheme Programming 3 November 7, 2019 23:34
udf problem jane Fluent UDF and Scheme Programming 37 February 20, 2018 04:17
UDF velocity profile willroca Fluent UDF and Scheme Programming 2 January 10, 2016 03:13
Error messages atg enGrid 7 August 30, 2013 11:16
Phase locked average in run time panara OpenFOAM 2 February 20, 2008 14:37


All times are GMT -4. The time now is 12:57.