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

interFoam simulation blowing up

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

Like Tree12Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 21, 2012, 10:24
Default
  #41
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16
mgdenno is on a distinguished road
Hi ardjouna,

You haven't really provided enough information for me to help you. Perhaps you should have a look at this thread.

From what you said I would guess that maybe either your OF environment variables are not set properly or possibly you are using a recent version of OF but following a tutorial from an older version; recent versions use interFoam and not rasInterFoam. I think rasInterFoam may be from older versions or maybe from the extend version. I am not sure but have seen it referenced in some older posts.

If you type the following at the command:

Code:
echo $FOAM_INST_DIR
what do you see?
mgdenno is offline   Reply With Quote

Old   October 21, 2012, 13:51
Default
  #42
New Member
 
yakouna
Join Date: Oct 2012
Posts: 2
Rep Power: 0
ardjouna is on a distinguished road
HI Mathiew,

Thank you for your rapid reply excuse me for not being clear
error message is:this application is not available i am working with OpenFOAM-1.5.00b-wininst it works correctly with cavity tutorial because icofoam.exe exist in Bin folder but interfoam and also rasinterfoam doesn't exist is there any way to get it for windows ?

best regards
ardjouna is offline   Reply With Quote

Old   October 26, 2012, 06:40
Default
  #43
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 240
Rep Power: 17
vonboett is on a distinguished road
Quote:
Originally Posted by MOHAMMAD67 View Post
Dear Oliver
I did all things except extending the domain. Unfortunately it didn't work and after 1.3 seconds it blows out. I don't know how deal with. I reduced the delta T to 0.00001 and turn the adjustable time step off. It is running now. I will inform you from the result.
Maybe I should make the mesh finer. Whats your opinion. Does it help me to get result.
Kind Regard
Hi Mohammad,

the issue with phase reflection at the outlet appears to several users, and type buoyantPressure; value uniform 0; usually helps. There is a funny effect: could you move your mesh to the negative coordinate quadrant (all vertices have x coordinate < 0) and check what happens at the outlet if you use type outletInlet; outletValue uniform 0; for pressure at the open boundary?
vonboett is offline   Reply With Quote

Old   October 26, 2012, 09:48
Default
  #44
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16
mgdenno is on a distinguished road
Quote:
Originally Posted by ardjouna View Post
HI Mathiew,

Thank you for your rapid reply excuse me for not being clear
error message is:this application is not available i am working with OpenFOAM-1.5.00b-wininst it works correctly with cavity tutorial because icofoam.exe exist in Bin folder but interfoam and also rasinterfoam doesn't exist is there any way to get it for windows ?

best regards
I believe there are a few more recent versions available for Windows. I would suggest Googling it. The other option is to run a virtual machine with Linux on it. This approach ios covered on the OpenFOAM website I believe.
mgdenno is offline   Reply With Quote

Old   November 20, 2012, 00:29
Default
  #45
New Member
 
Tsun-Hua Yang
Join Date: Nov 2012
Posts: 4
Rep Power: 14
Josh Yang is on a distinguished road
Hey,

I wonder did you solve your problem yet?
I am kind of running into a same issue as yours.

Thanks.

Josh
Josh Yang is offline   Reply With Quote

Old   November 21, 2012, 06:27
Default
  #46
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 240
Rep Power: 17
vonboett is on a distinguished road
Quote:
Originally Posted by vonboett View Post
Hi Mohammad,

the issue with phase reflection at the outlet appears to several users, and type buoyantPressure; value uniform 0; usually helps. There is a funny effect: could you move your mesh to the negative coordinate quadrant (all vertices have x coordinate < 0) and check what happens at the outlet if you use type outletInlet; outletValue uniform 0; for pressure at the open boundary?
Actually, zeroGradient for pressure at the outlet works best for me.
vonboett is offline   Reply With Quote

Old   March 29, 2013, 13:41
Default
  #47
cjz
New Member
 
Join Date: Feb 2010
Posts: 14
Rep Power: 16
cjz is on a distinguished road
Hello mgdenno,
I'm wondering if you ever got your simulation to run to completion? I've been working on a similar open channel flow for the past couple of weeks and find my simulation fails part way through. There is a sudden velocity magnitude that is orders of magnitude higher than the previous. I've set cAlpha = 0 but the problem remains. If you've been successful I'd be interested in knowing what combination of fvSolution, fvSchemes, and boundary conditions worked for you.

Thanks for any help you can provide.
cjz is offline   Reply With Quote

Old   March 29, 2013, 20:38
Default
  #48
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16
mgdenno is on a distinguished road
I think what worked in my case was to initially use:

Code:
gradSchemes
{
    default         Gauss linear;
    grad(U)        cellLimited Gauss linear 1; 
}

divSchemes
{
    div(rho*phi,U) Gauss linearUpwindV grad(U);
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;
    div(phi,k)      Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(phi,R)      Gauss upwind;
    div(R)          Gauss linear;
    div(phi,nuTilda) Gauss upwind;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}
Hope that helps.

Matt
lourencosm and jylee4 like this.
mgdenno is offline   Reply With Quote

Old   August 9, 2013, 08:59
Default Instability in Open channel flow with interFoam
  #49
New Member
 
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 14
sorockin is on a distinguished road
Dear Matthew,

I try to simulate flow in open channel with interFoam.
http://www.cfd-online.com/Forums/ope...interfoam.html
My case is 2D and I used BCs from your working case. But I can't overcome instability that arise at inlet. May be you can help me with advice? Unlike your case I turned off turbulence.

Thank you.
Maxim
sorockin is offline   Reply With Quote

Old   August 9, 2013, 10:10
Default
  #50
cjz
New Member
 
Join Date: Feb 2010
Posts: 14
Rep Power: 16
cjz is on a distinguished road
Have you tried adjusting cAlpha?

from:

http://www.openfoam.org/docs/user/damBreak.php

The cAlpha keyword is a factor that controls the compression of the interface where: 0 corresponds to no compression; 1 corresponds to conservative compression; and, anything larger than 1, relates to enhanced compression of the interface. We generally recommend a value of 1.0 which is employed in this example.
cjz is offline   Reply With Quote

Old   August 9, 2013, 11:15
Default fvSolushion
  #51
New Member
 
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 14
sorockin is on a distinguished road
My cAlpha=1 was and now. But now I took fvSolushion setting from Matthew and looks like issue is solved.
http://www.cfd-online.com/Forums/ope...tml#post444738

Thank you Matthew!

Maxim
sorockin is offline   Reply With Quote

Old   December 6, 2017, 13:28
Default
  #52
Member
 
Honza Höll
Join Date: Mar 2016
Location: Brno, CZ
Posts: 37
Rep Power: 10
indy07cz is on a distinguished road
Hello foamers,
I have a little question about interfaceCompression divScheme for phirb,alpha. When I use this scheme values of Max(alpha.water) go high to numbers like 1.9-2.2 during simulation and due to I must drastically decrease maxAlphaCo to achieve Max(alpha.water) equal to 1 (then simulation takes long time to compute). If I use linear scheme Max(alpha.water) is very close to 1. Can anybody confirm that and tell me why is this happening? Is it somehow connected with mesh quality? Thanks.
indy07cz is offline   Reply With Quote

Old   December 7, 2017, 08:34
Default
  #53
Senior Member
 
Saideep
Join Date: Apr 2015
Location: INDIA
Posts: 203
Rep Power: 12
Saideep is on a distinguished road
Not sure about your case study but in short to run interfoam cases you generally need to use very small time steps. Courant number of 0.1 is already too large.

There is "Brackbill number" that determines the time step size for volume of fluid computations. I believe/ never had a problem with "interfaceCompression" for the "artificial compression advection term". Probably your also running into classic "spurious current" problem or "volumeRatio" is too high when you used "snappyHexMesh".

So, need more description from your end!!
BlnPhoenix likes this.
Saideep is offline   Reply With Quote

Old   December 7, 2017, 09:13
Default
  #54
Member
 
Honza Höll
Join Date: Mar 2016
Location: Brno, CZ
Posts: 37
Rep Power: 10
indy07cz is on a distinguished road
Well I simulate flow over stepped spillway with flow over 200 cms. Mesh has about 2mil cells, and yes there could be problem in mesh volume ratio because smallest cell is 4.3E-5 cubic meters (around spillway crest) and largest is 2.2E-2. CheckMesh is OK. I use Co and alphaCo equal to 0.6 because of computing time (60 seconds takes 2 days). I now when I decrease Co to something about 0.2 then alpha.water is lets say fine, but 5 seconds of simulation takes 1-2 days and I don't have too much time.

I attached some pics to illustrate problem. You can see when I clip results in paraview by threshold for alpha 0-1 lots of cells are missing mainly near wall surfaces.

So I think that biggest problem is high Co number for this simulation.

So I made some research, and I suppose that problem is fine mesh on surfaces in combination with high Co number. I need lower Co or make mash more coarse. Interesting is when i use linear scheme my p_rgh residuals are higher than residuals with interfaceCompression scheme but Max(alpha.water) is higher with interfaceCompression.
Attached Images
File Type: png clip.png (121.7 KB, 59 views)
File Type: png threshold.png (63.6 KB, 40 views)

Last edited by indy07cz; December 7, 2017 at 11:01. Reason: Idea
indy07cz is offline   Reply With Quote

Old   December 19, 2017, 04:03
Default
  #55
Member
 
Honza Höll
Join Date: Mar 2016
Location: Brno, CZ
Posts: 37
Rep Power: 10
indy07cz is on a distinguished road
Hi, after some simulations I found solution for my problem around boundening alpha.water solutions. Even if I reduced Co my Max(alpha.water) was still above 1. Then I added MULES correction to alpha.water in fvSolution and now everything is fine and solution is bounded. But without it I think high alpha values is signal for large timestep and simulation needs Co reduction!

Now the code looks like this:
Code:
"alpha.water.*"
    {
        nAlphaCorr      1;
        nAlphaSubCycles 2; 
        cAlpha          1;

        MULESCorr       yes;  
        nLimiterIter    3;

        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-8;
        relTol 0;
    }
blahblahblah likes this.
indy07cz is offline   Reply With Quote

Old   November 2, 2018, 21:20
Default some tutorial or project of an open water channel with transition?
  #56
New Member
 
Jesus Laor Flores
Join Date: Nov 2018
Posts: 1
Rep Power: 0
jesucitolf is on a distinguished road
some tutorial or project of an open water channel with transition
im new in openfoam, i want to learn plss. some normal water channel or with transition pls...thx
jesucitolf is offline   Reply With Quote

Reply

Tags
interfoam, spillways

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
Solar Radiation in OpenFOAM plainstyle OpenFOAM Running, Solving & CFD 15 July 8, 2014 05:43
GUI crash and simulation engine still running RPJones FLOW-3D 2 November 9, 2010 09:18
velocity profile export from a simulation onto another sudhirlv STAR-CCM+ 1 September 12, 2010 19:57
Continuous vs interrupted simulation sega OpenFOAM Running, Solving & CFD 4 November 3, 2008 15:29
strange simulation error Ralf Schmidt FLUENT 2 May 4, 2007 14:02


All times are GMT -4. The time now is 07:18.