|
[Sponsors] |
October 21, 2012, 10:24 |
|
#41 |
Senior Member
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16 |
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 |
|
October 21, 2012, 13:51 |
|
#42 |
New Member
yakouna
Join Date: Oct 2012
Posts: 2
Rep Power: 0 |
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 |
|
October 26, 2012, 06:40 |
|
#43 | |
Senior Member
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 240
Rep Power: 17 |
Quote:
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? |
||
October 26, 2012, 09:48 |
|
#44 | |
Senior Member
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16 |
Quote:
|
||
November 20, 2012, 00:29 |
|
#45 |
New Member
Tsun-Hua Yang
Join Date: Nov 2012
Posts: 4
Rep Power: 14 |
Hey,
I wonder did you solve your problem yet? I am kind of running into a same issue as yours. Thanks. Josh |
|
November 21, 2012, 06:27 |
|
#46 | |
Senior Member
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 240
Rep Power: 17 |
Quote:
|
||
March 29, 2013, 13:41 |
|
#47 |
New Member
Join Date: Feb 2010
Posts: 14
Rep Power: 16 |
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. |
|
March 29, 2013, 20:38 |
|
#48 |
Senior Member
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16 |
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; } Matt |
|
August 9, 2013, 08:59 |
Instability in Open channel flow with interFoam
|
#49 |
New Member
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 14 |
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 |
|
August 9, 2013, 10:10 |
|
#50 |
New Member
Join Date: Feb 2010
Posts: 14
Rep Power: 16 |
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. |
|
August 9, 2013, 11:15 |
fvSolushion
|
#51 |
New Member
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 14 |
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 |
|
December 6, 2017, 13:28 |
|
#52 |
Member
Honza Höll
Join Date: Mar 2016
Location: Brno, CZ
Posts: 37
Rep Power: 10 |
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. |
|
December 7, 2017, 08:34 |
|
#53 |
Senior Member
Saideep
Join Date: Apr 2015
Location: INDIA
Posts: 203
Rep Power: 12 |
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!! |
|
December 7, 2017, 09:13 |
|
#54 |
Member
Honza Höll
Join Date: Mar 2016
Location: Brno, CZ
Posts: 37
Rep Power: 10 |
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. Last edited by indy07cz; December 7, 2017 at 11:01. Reason: Idea |
|
December 19, 2017, 04:03 |
|
#55 |
Member
Honza Höll
Join Date: Mar 2016
Location: Brno, CZ
Posts: 37
Rep Power: 10 |
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; } |
|
November 2, 2018, 21:20 |
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 |
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 |
|
Tags |
interfoam, spillways |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |