|
[Sponsors] |
October 7, 2016, 02:41 |
rhoPimpleFoam hardship
|
#1 |
New Member
|
https://www.dropbox.com/home?preview...%3A24%3A15.png
https://www.dropbox.com/home?preview...%3A29%3A07.png Hi. My first post, I hope give all necessary information to get a pleasant, good and quickly answer. Thanks in advance. So, trying simulate a jet injection inside a combustion chamber (without combustion of course). Nitrogen gas. The mesh is axisymmetric, he nozzle of injector and the combustion chamber. Split in 3 blocks. blocks ( hex ( 0 2 2 0 1 3 9 8) (20 1 8) simpleGrading (1.00 1.00 0.3) // block number0 hex ( 2 5 5 2 3 6 11 9) (100 1 8) simpleGrading (10.00 1.00 0.3) // block number1 hex ( 3 6 11 9 4 7 12 10) (100 1 25) simpleGrading (10.00 1.00 35.00) // block number2 ); It run correctly, but the animation is not as expected, better said it is very strange. p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 2.14e6; boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value uniform 2.14e6; } wall { type zeroGradient; } top { type zeroGradient; } frontface { type wedge; } backface { type wedge; } axis { type empty; } } --------------------------------------------------------------------------------------------- object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 298; boundaryField { inlet { type fixedValue; value uniform 298; } outlet { type zeroGradient; } wall { type zeroGradient; } top { type zeroGradient; } frontface { type wedge; } backface { type wedge; } axis { type empty; } } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * / object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type fixedValue; value uniform (400 0 0); } outlet { type pressureInletOutletVelocity; value uniform (0 0 0); } wall { type fixedValue; value uniform (0 0 0); } top { type zeroGradient; } frontface { type wedge; } backface { type wedge; } axis { type empty; } } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // object sig; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 -1 0 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type fixedValue; value uniform 0; } outlet { type zeroGradient; } wall { type fixedValue; value uniform 0; } top { type zeroGradient; } frontface { type wedge; } backface { type wedge; } axis { type empty; } } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // object alphat; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -1 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type calculated; value uniform 0; } outlet { type calculated; value uniform 0; } top { type calculated; value uniform 0; } wall { type compressible::alphatWallFunction; Prt 0.85; value uniform 0; } frontface { type wedge; } backface { type wedge; } axis { type empty; } } // ************************************************** *********************** // object epsilon; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -3 0 0 0 0]; internalField uniform 14903; boundaryField { inlet { type fixedValue; value uniform 14903; } outlet { type zeroGradient; } top { type zeroGradient; } wall { type compressible::epsilonWallFunction; value uniform 14903; Cmu 0.09; kappa 0.41; E 9.8; value uniform 14903; } frontface { type wedge; } backface { type wedge; } axis { type empty; } } // ************************************************** *********************** // object k; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 37.5; boundaryField { inlet { type fixedValue; value uniform 37.5; } outlet { type zeroGradient; } top { type zeroGradient; } wall { type compressible::kqRWallFunction; value uniform 37.5; } frontface { type wedge; } backface { type wedge; } axis { type empty; } } // ************************************************** *********************** // object mut; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -1 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type calculated; value uniform 0; } outlet { type calculated; value uniform 0; } top { type calculated; value uniform 0; } wall { type mutkWallFunction; Cmu 0.09; kappa 0.41; E 9.8; value uniform 0; } frontface { type wedge; } backface { type wedge; } axis { type empty; } } // ************************************************** *********************** // object nut; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type calculated; value uniform 0; } outlet { type calculated; value uniform 0; } top { type calculated; value uniform 0; } wall { type nutkWallFunction; Cmu 0.09; kappa 0.41; E 9.8; value uniform 0; } frontface { type wedge; } backface { type wedge; } axis { type empty; } } // ************************************************** *********************** // object omega; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 -1 0 0 0 0]; internalField uniform 1118000; boundaryField { inlet { type fixedValue; value uniform 1118000; } outlet { type zeroGradient; } top { type zeroGradient; } wall { type omegaWallFunction; Cmu 0.09; kappa 0.41; E 9.8; beta1 0.075; value uniform 1118000; } frontface { type wedge; } backface { type wedge; } axis { type empty; } } // ************************************************** *********************** // blockMesh is ok, checkMesh also is ok. I think the problem is in the boundary conditions, but I dunno where. Thanking again any suggestion. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
rhoPimpleFoam Boundary Condition Problem | dancfd | OpenFOAM Pre-Processing | 18 | September 16, 2021 07:43 |
rhoPimpleFoam solver | adrieno | OpenFOAM Running, Solving & CFD | 11 | April 6, 2016 11:01 |
Pressure stair-step behaviour using rhopimplefoam | joegi.geo | OpenFOAM Running, Solving & CFD | 3 | December 12, 2014 12:10 |
rhoPimpleFoam floating point error | dancfd | OpenFOAM Running, Solving & CFD | 6 | January 5, 2014 20:57 |
Questions about buoyantPimpleFoam and rhoPimpleFoam | Mojtaba.a | OpenFOAM Running, Solving & CFD | 6 | August 1, 2012 04:50 |