pimpleDyMFoam stability problems
4 Attachment(s)
I'm having some difficulties getting a pimpleDyMFoam calculation to run. It will run for a long simulation time (as in thousands of time steps) but the calculation always seems to suddenly blow up. However, it needs to run longer and with second order accurate schemes. A quick summary of what I'm trying to do (I can't show a picture of the actual geometry):
- running 1.6-ext, pimpleDyMFoam - 2D incompressible flow (V_inf = 10 m/s, L_ref ~= 0.5 m, air) - uniform, horizontal flow with a cylindrical rotating section (omega = 3000 RPM, R_ref ~= 0.053 m => V_tip ~= 2.65 m/s) - turbulent flow, using standard k-e - using a GGI between the non-rotating/rotating sections. - initial condition is an MRFSimpleFoam solution. - running adaptive time-stepping with CFL = 2.0 The mesh is a mix of quad and tri cells (quad in BL regions, tri's elsewhere). Here's the last bit of the checkMesh output: Code:
Checking geometry... I've attached my controlDict, fvSchemes, fvSolution and dynamicMeshDict for completeness but here's a short summary for each: fvSchemes: ddt: Euler grad: cellMDLimited Gauss linear 1.0 div: upwind As you can see, pretty much first order for everything (I'm just trying to get it to run!) fvSolution: p: GAMG (tolerance = 1e-8, relTol = 0) everything else: PBiCG (tolerance = 1e-12, relTol = 0) nOuterCorrectors: 2 nCorrectors: 2 nNonOrthogonalCorrectors: 0 Here's the output from the second last time before it blows up: Code:
Initializing the GGI interpolator between master/shadow patches: outerGGI/innerGGI I've checked the mechanics of everything and all seems to be working (e.g. the ggi mass flux error is small). I looked at several instances of the solution prior to it blowing up with Paraview: the rotating zone is moving correctly and there are no spikes in velocity, pressure or the turbulence quantities. This all leads me to believe it's has to do with the pressure equation. I've tried different # of nCorrectors, outerCorrectorsn NonOrthogonalCorrectors and numerical schemes (I'm using the most stable ones already!) but it doesn't seem help. As it appears to be the pressure equation blowing up, what should I tweak to get better control/convergence of the pressure? Sorry for the long post but I'm running out of things to try so I was hoping some fresh eyes and someone with more knowledgeable about the inner workings of pimple* type solvers might be able to send me down a better rabbit hole. Let me know if you need more information. |
Hello,
I have a similar issue. I run turbDyMFoam in 1.5-dev, and I also use GAMG for my p calculation. This also blows up my number of iterations, and my speed increases very suddenly to unphysically high values. I am running the same case with PBiCG for p as well as for the other variables, and until now it seems stable (but is much slower). Maybe there is a difficulty in combining GGI with grid coarsening? |
Thanks for the tip flowris. I tried that but unfortunately that didn't make a difference for me.
I took a look at how the pimple*Foam works. It appears the outer corrector updates all the equations whereas the inner corrector updates only the pressure equation. Here's the code from pimpleFoam: Code:
for (int oCorr=0; oCorr<nOuterCorr; oCorr++) Futhermore, I've finally been able to switch over to second-order accurate div and ddt schemes. X'ing my fingers ... |
Another issue I've discovered with pimple*Foam is that you still can not use large CFL numbers. At least for my applications I've found anything over a CFL of 2 and I start having stability issues. I suppose it's better than a PISO algorithm but you still have keep the CFL within reason.
|
All times are GMT -4. The time now is 10:31. |