# set Min and Max values for pressure in pimpleFoam

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

 September 2, 2013, 17:57 set Min and Max values for pressure in pimpleFoam #1 Senior Member   Join Date: Jun 2011 Posts: 151 Rep Power: 7 Hi all in my problem the velocities are calculated correctly but pressure is not. I use pimpleFoam and I want to limit the pressure I add bellow statment in fvSolution/PIMPLE pMin pMin [ 0 2 -2 0 0 0 0 ] -150; pMax pMax [ 0 2 -2 0 0 0 0 ] 150; but it not work any help will be appreciated

 September 3, 2013, 00:13 #2 Senior Member   Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,195 Blog Entries: 1 Rep Power: 16 Hello did you change your solver at all!? then those lines would not reflect anything to solver besides as you use incopressible solver i dont gauss its a good idea to limit the pressure which will be effect on continuity too and it may impair divergence issue __________________ Telegram channel (https://telegram.me/openfoam4Iranian) My Weblog (http://openfoam.blogfa.com/) Training Course on OpenFOAM at (http://www.isme.ir/)

September 3, 2013, 01:05
#3
Senior Member

Join Date: Jun 2011
Posts: 151
Rep Power: 7
Hi Nima

I have added the setting of my test case as follow:

I use my icoFsiFoam solver which used pimpleFoam for solve fluid.
the results is good until time=0.1205 with maxCor=.6 ( deltaT ~1e-5)
at the next state (Time = 0.1205+deltaT) the velocity is calculated correctly, but pressure not.
I check the original pimpleFoam solver to solve fluid only, but the pressure is diverged again. I have attached its contours.

fvSchemes

HTML Code:
```ddtSchemes
{
default Euler;
}

{
default           Gauss linear ;//cellMDLimited leastSquares 1;//extendedLeastSquares 1;
}

divSchemes
{
default         none;
div(phi,U)      Gauss limitedLinearV 1;//Gauss linearUpwindV grad(U);//
div(phi,k)      Gauss upwind;//limitedLinear 1;
div(phi,omega)  Gauss upwind;//limitedLinear 1;
}

laplacianSchemes
{
default         Gauss linear corrected;//Gauss linear limited 0.5;//
}

interpolationSchemes
{
default         linear;
}

{
default         corrected;
}

fluxRequired
{
default         no;
pcorr           ;
p;
}

// ************************************************************************* //```

fvSolutions

HTML Code:
```solvers
{
p
{
solver           GAMG;
tolerance        1e-4;
relTol           0;
smoother         GaussSeidel;
nPreSweeps       1;
nPostSweeps      2;
cacheAgglomeration true;
directSolveCoarsest true;
agglomerator     faceAreaPair;
nCellsInCoarsestLevel 40;
mergeLevels      1;
minIter          1;
}
pFinal
{
solver           GAMG;
tolerance        1e-8;
relTol           0;
smoother         GaussSeidel;
nPreSweeps       1;
nPostSweeps      2;
cacheAgglomeration true;
agglomerator     faceAreaPair;
nCellsInCoarsestLevel 40;
mergeLevels      1;
minIter          1;
}

"(U|k|omega)"
{
solver          PBiCG;
preconditioner  DILU;

tolerance        1e-5;
relTol           0;
nSweeps          1;
minIter          1;
}
"(UFinal|kFinal|omegaFinal)"
{

solver           smoothSolver;
smoother         GaussSeidel;
tolerance        1e-6;
relTol           0.0;
nSweeps          1;
minIter          1;
}
}

PIMPLE
{
nOuterCorrectors    3;
nCorrectors         2;
nNonOrthogonalCorrectors 3;

pRefCell            0;
pRefValue           0;
//    momentumPredictor yes;
ddtPhiCorr no;
correctPhi no;

}

relaxationFactors
{
fields
{
p               0.3;
}
equations
{
U               0.7;
k               0.7;
omega           0.7;
}
}
cache
{
}

// ************************************************************************* //```
Attached Images
 P_1205OK.jpg (21.0 KB, 15 views) P_120501_notOK.jpg (17.0 KB, 12 views) U_120501_OK.jpg (17.8 KB, 17 views)

 September 3, 2013, 02:24 #4 Senior Member   Bernhard Join Date: Sep 2009 Location: Delft Posts: 790 Rep Power: 14 Bounding pressure is certainly not the way to go. Check your velocity in a few more time steps, and you will probably see something weird there as well. In other words, your solution is diverging, please check the residuals of the equations from the log file. It might be wise to study the origin of the pressure issue, which might well be related to the FSI part. nimasam and FrankFlow like this.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post florian_krause OpenFOAM 110 April 21, 2016 11:54 liguifan OpenFOAM Running, Solving & CFD 5 June 3, 2014 02:53 Elise OpenFOAM Native Meshers: snappyHexMesh and Others 1 April 22, 2013 02:32 paean OpenFOAM Running, Solving & CFD 0 November 14, 2008 22:14 Sans CFX 2 May 31, 2008 04:22

All times are GMT -4. The time now is 17:32.