|
[Sponsors] |
July 2, 2021, 13:44 |
Viscoplastic Fluid Flow with OpenFOAM
|
#1 |
New Member
Elliott Sutton
Join Date: Dec 2019
Posts: 5
Rep Power: 7 |
Hi all,
I'm trying to solve some steady pressure-driven pipe flows of Herschel-Bulkley fluids using simpleFOAM as a steady state solver. I'm finding that reaching convergence takes an incredibly long time. I've optimised the schemes and solvers so that, within each outer iteration, U converges in <10 iterations and p takes 1 or 2 iterations. Despite this, the residuals decrease at a very slow rate. Relaxation parameters of 0.7 for both U and p are the largest that I can use without the simulation blowing up. Does anyone have experience with these types of flows using OpenFOAM and has advice on how to improve convergence? I will post the fvSchemes and fvSolution files below. If any more info is needed let me know. Thanks in advance. Code:
ddtSchemes { default steadyState; } gradSchemes { default cellMDLimited Gauss linear 1; } divSchemes { default none; div(phi,U) bounded Gauss linearUpwindV grad(U); div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } Code:
solvers { p { solver GAMG; tolerance 1e-06; relTol 0.01; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration true; agglomerator faceAreaPair; nCellsInCoarsestLevel 10; mergeLevels 1; } U { solver smoothSolver; smoother GaussSeidel; tolerance 1e-06; relTol 0.01; nSweeps 1; } } SIMPLE { nNonOrthogonalCorrectors 1; consistent yes; /*pRefCell 0; pRefValue 0;*/ } relaxationFactors { fields { p 0.7; } equations { U 0.7; } } |
|
Tags |
convergence, herschel-bulkley, simplefoam, steady state, yield stress |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Map of the OpenFOAM Forum - Understanding where to post your questions! | wyldckat | OpenFOAM | 10 | September 2, 2021 06:29 |
mass flow in is not equal to mass flow out | saii | CFX | 12 | March 19, 2018 06:21 |
How to contribute to the community of OpenFOAM users and to the OpenFOAM technology | wyldckat | OpenFOAM | 17 | November 10, 2017 16:54 |
Multiphase flow - incorrect velocity on inlet | Mike_Tom | CFX | 6 | September 29, 2016 02:27 |
An error has occurred in cfx5solve: | volo87 | CFX | 5 | June 14, 2013 18:44 |