[B6_lacen]

Good morning,

When I try to run the PimpleFoam solver my code crashes, I don't know why, as you can see in photos the dt goes from a value of 1e-8 to 1e-104.

Could you help me please I am new in CFD.

Thanks

Capture du 2022-02-28 10-53-47.jpg
Capture du 2022-02-28 10-54-13.png

[SHUBHAM9595]

Hi Boucetta, I hope that you have already solved ur issue.

If not, please provide more info to the users about your test case. Can you attach the output of checkMesh command.

Also, details about fvSchemes, fvSolution and controldict would be helpful along with ur initial & boundary conditions.

[B6_lacen]

Hello SHUBHAM, Thanks for your replay.

controlDict.png

CheckMesh8.png

fvschemes.jpg

fvsolution_.png

fvsolution2.png

Thas it my fiels, it take a lote of time, and after it's damped. my Re= 50000

Thanks a lot i apreciate.

[B6_lacen]

Hello again
I forgot initial condition, I wonder if in the fvscheme folder i can use the 2nd order schema?

I don't want to use a wall law because my y+=0,3

Thanks

k.png

nut.png

omega.jpg

U.jpg

[SHUBHAM9595]

Now, we have a better idea about your test case .


Few Comments


1) Try with Euler first or maybe Crank - Nicolson 0.5 & then shift to backwards

2) You can definitely opt for one non - orthogonal correctors. Try going for higher number of outerCorrectors like 50 with the residual control options...

PIMPLE
{
    momentumPredictor true;
    transonic       false;
    consistent		true;// no need to use urf with consistent 
    nOuterCorrectors 50;//50 
    nCorrectors     3;
    nNonOrthogonalCorrectors 1;
    
    residualControl
    {

 "(p|p_rgh)"
        {              // If this INITIAL tolerance is reached , leave outer corrector
			    tolerance  1e-03;
                relTol     0;
        }
        
              U
        {              // If this INITIAL tolerance is reached , leave outer corrector
			    tolerance  1e-03;
                relTol     0;
        }    
        

     }
}

3) You can check from ur first image that the FPE is coming from the GAMG solver which is used for the pressure, but we still don't see your initial conditions for pressure.

4) I've no idea about your geometry, but based on the naming conventions....inletOutlet for BOTH ymin & ymax might not be correct...


P.S. No expertise in turbulence so cant comment on the k, omega, & nut. I presume that aspect ratio of 700 is not the culprit for divergence while simulating turbulence