CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Viscoplastic Fluid Flow with OpenFOAM

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 2, 2021, 12:44
Default Viscoplastic Fluid Flow with OpenFOAM
  #1
New Member
 
Elliott Sutton
Join Date: Dec 2019
Posts: 5
Rep Power: 6
Ells96 is on a distinguished road
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;
    }
}
Ells96 is offline   Reply With Quote

Reply

Tags
convergence, herschel-bulkley, simplefoam, steady state, yield stress

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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 05:29
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 05:21
How to contribute to the community of OpenFOAM users and to the OpenFOAM technology wyldckat OpenFOAM 17 November 10, 2017 15:54
Multiphase flow - incorrect velocity on inlet Mike_Tom CFX 6 September 29, 2016 01:27
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44


All times are GMT -4. The time now is 16:56.