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

Poor convergence for Energy e (rhoSimpleFoam, Steady-State)

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 26, 2024, 10:05
Default Poor convergence for Energy e (rhoSimpleFoam, Steady-State)
  #1
New Member
 
S03r3n
Join Date: Feb 2024
Posts: 15
Rep Power: 2
S03r3n is on a distinguished road
Hello together,

I have the issue that my simulation (rhoSimpleFoam, steady-state) shows good convergence for velocities U and pressure p but poor convergence for energy e.
What can be the reason? Is it solver related?

I have attached the convergence plot.

Below you can find the thermoType I'm using and the chosen solver settings.

thermophysicalProperties
Code:
thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture
{
    specie
    {
        nMoles      1;
        molWeight   600000;
    }
    thermodynamics
    {
        Cp          1850;
        Hf          0;
    }
    transport
    {
        mu          40000;
        Pr          370000000;
    }
    equationOfState
    {
        rho         1154;
    } 
}
controlDict
Code:
application     rhoSimpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         300;

deltaT          1;

writeControl    timeStep;

writeInterval   100;

purgeWrite      1;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;
fvSolution
Code:
solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-10;
        relTol          0.1;
        smoother        GaussSeidel;
    }

    "(U|e)"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        tolerance       1e-10;
        relTol          0.1;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 2;
}

relaxationFactors
{
    fields
    {
        p               0.3;
        rho             0.3;   
    }
    equations
    {
        p               0.3;
        U               0.7;
        e               0.7;
    }
}
fvSchemes
Code:
ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      bounded Gauss linearUpwind grad(U);
    div(((rho*nuEff)*dev2(T(grad(U)))))      Gauss linear; 
    
    div(phi,e)          bounded Gauss linearUpwind grad(e)

    div(phid,p)         Gauss upwind;
    div(phi,Ekp)        bounded Gauss upwind;
    div((phi|interpolate(rho)),p)  Gauss upwind;  
}

laplacianSchemes
{
    default         Gauss linear limited 0.33;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.33;
}
Attached Images
File Type: jpg rhoSimpleFoam_poorEnergyConvergence.jpg (76.1 KB, 15 views)
S03r3n is offline   Reply With Quote

Old   June 26, 2024, 11:54
Default
  #2
Senior Member
 
Join Date: Dec 2021
Posts: 251
Rep Power: 5
Alczem is on a distinguished road
Hey!


Your thermophysicalProperties values seem weird to me, what kind of fluid has such molecular weight, viscosity and prandtl? Are you sure about these properties?


Apart from that, you could maybe add a limiter to the gradient, but everything looks ok to me
Alczem is offline   Reply With Quote

Old   June 26, 2024, 12:49
Default
  #3
New Member
 
S03r3n
Join Date: Feb 2024
Posts: 15
Rep Power: 2
S03r3n is on a distinguished road
Quote:
Originally Posted by Alczem View Post
Hey!


Your thermophysicalProperties values seem weird to me, what kind of fluid has such molecular weight, viscosity and prandtl? Are you sure about these properties?


Apart from that, you could maybe add a limiter to the gradient, but everything looks ok to me

Thanks for your message.
I try to simulate uncured rubber.
Do you have some suggestions?
S03r3n is offline   Reply With Quote

Old   June 28, 2024, 18:57
Default
  #4
Senior Member
 
Join Date: Apr 2020
Location: UK
Posts: 742
Rep Power: 14
Tobermory will become famous soon enough
Just reiterating Alczem's earlier post ... lead has a molecular weight of 207 kg/kmol, and you are trying to apply 6x10^5 kg/kmol! Your density is approx the density of lead and your viscosity is off the scale. None of this seems very physical (although I appreciate that uncured rubber might be very viscous).

So again - as per Alczem's post, can I also suggest you check your inputs.
Tobermory is offline   Reply With Quote

Old   July 1, 2024, 10:51
Default
  #5
New Member
 
S03r3n
Join Date: Feb 2024
Posts: 15
Rep Power: 2
S03r3n is on a distinguished road
I found out that the convergence behavior of energy e highly depends on Prandtl number, thermal conductivity respectively.

I have changed the transport model to polynomial in order to specify the thermal conductivity (kappa).

However, I get bad e-convergence for this case:
Code:
    transport
    {
        muCoeffs<8>     ( 15000.0 0 0 0 0 0 0 0 );
        kappaCoeffs<8>  ( 0.15 0 0 0 0 0 0 0 );
    }
And good e-convergence for this case:
Code:
    transport
    {
        muCoeffs<8>     ( 15000.0 0 0 0 0 0 0 0 );
        kappaCoeffs<8>  ( 150 0 0 0 0 0 0 0 );
    }
I thought that the unit of kappa should be W/(m*K), so the bad case should have the correct value of 0.15 W/(m*K)?
Any suggestions?

Best regards
S03r3n is offline   Reply With Quote

Reply

Tags
convergence, energy, openfoam

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
Setting the height of the stream in the free channel kevinmccartin CFX 12 October 13, 2022 22:43
Domain Reference Pressure and mass flow inlet boundary AdidaKK CFX 75 August 20, 2018 06:37
Convergence in steady state simulations vs transient ones cardioCFD CFX 5 January 21, 2018 11:59
Constant velocity of the material Sas CFX 15 July 13, 2010 09:56
About the difference between steady and unsteady problems Lisa Main CFD Forum 11 July 5, 2000 15:37


All times are GMT -4. The time now is 01:20.