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

steadyState - Euler - different endresults

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 30, 2021, 06:18
Default steadyState - Euler - different endresults
  #1
New Member
 
Join Date: Jan 2020
Posts: 26
Rep Power: 6
PSander is on a distinguished road
I have problems to understand if my results do make sense or not.
I have the geometry of a square, x=6, y=3, z=100 meters. The upper surface works as a source of heat, where all the time 303°K are streaming in. All the other sides are zero gradient, and the internal temperature is 293°K. The thermophysical parameters are similar to sand. (kappa =2, Cv=1000, rho=2200)
I run with chtMultiRegionFoam and OF8.

Inside the square I have 2 tubes filles with water, but for this simulation the temperatures are as well 293°K and the velocity is 0.

My expectacion is, that after a certain time everything should be heated up the surface temperature, since there is no other source of heat than the surface.

If I select steadyState exactly this happens, in the end everything is heated up to 303°K.

If I repeat the same calculation with Euler, it takes around 5 days, while the square is constantly heating up from the upper surface, when the heating process seems to be stopping. Actually its not stopping completly, but the diference for the same cell between day 5 and day 10 is like 0.000001°K, so almost nothing.
But the square is far away from having the upper surface temperature, in the lower regions it is still at 293°K, the starting value.

So I expected, that steadyState and Euler produce in the end the same result, is this not true?

Because if it is true, than I dont know why my calculation is with Euler not coming to that point. There must be some problem in my settings?

Do you have any idea if
1. Euler and steadyState come to the same results or not?
2. If yes, what could be the problem in my settings?


I expect problems in the system-files, so I show you my fvSolution, fvSchemes and controlDict.


fvSolution for region soil:

Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    e
    {
        solver           PCG;
        preconditioner   DIC;
        tolerance        1e-06;
        relTol           0.1;
    }

    eFinal
    {
        $e;
        tolerance        1e-06;
        relTol           0;
    }

}

PIMPLE
{
    nNonOrthogonalCorrectors 0;
    outerCorrectorResidualControl
    {
        e
        {
            tolerance 1e-5;
            relTol 0;
        }
    }
}

relaxationFactors
{
    equations
    {
        e           0.5;
        eFinal        1;
    }
}

// ************************************************************************* //

fvSchemes for region soil:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler; // Euler steadyState
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
}

laplacianSchemes
{
    default             none;
    laplacian(alpha,e)  Gauss linear uncorrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         uncorrected;
}


// ************************************************************************* //


controlDict

Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  8
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     chtMultiRegionFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         864000;

deltaT          1;

writeControl    adjustableRunTime;

writeInterval   3600;

purgeWrite      0;

writeFormat     ascii;

writePrecision  22;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  true; 

maxCo           1;

functions
{

  #includeFunc "residuals_fluid"

  avT_fluid_outlet
  {
      type            surfaceFieldValue;
      libs            ("libfieldFunctionObjects.so");
      log             yes;
      writeControl    runTime;
      writeInterval   3600;	
      writeFields     false;
      surfaceFormat   none;
      regionType      patch;
      name            outlet;
      region 	      fluid;	
      operation       areaAverage;
      fields
      (
          T
      );
  }

  yPlus_fluid
  {
    type            yPlus;
    libs            ("libfieldFunctionObjects.so");
    executeControl  writeTime;
    writeControl    writeTime;
    region  fluid;
  } 

  wallHeatFlux_fluid
  {
    type            wallHeatFlux;
    libs            ("libfieldFunctionObjects.so");
    executeControl  writeTime;
    writeControl    writeTime;
    region  fluid;
  }

  wallShearStress_fluid
  {
    type            wallShearStress;
    libs            ("libfieldFunctionObjects.so");
    executeControl  writeTime;
    writeControl    writeTime;
    region  fluid;
  }

  turbulenceIntensity_fluid
  {
    type            turbulenceIntensity;
    libs            ("libfieldFunctionObjects.so");
    executeControl  writeTime;
    writeControl    writeTime;
    region  fluid;
  } 

}


libs (
"libOpenFOAM.so"
"libsimpleSwakFunctionObjects.so"
"libswakFunctionObjects.so"
"libgroovyBC.so"
);

// ************************************************************************* //

general fvSolution
Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

PIMPLE
{
//    momentumPredictor yes;
    nOuterCorrectors 50;
    nCorrectors     2;
    nNonOrthogonalCorrectors 0;
    consistent      true;
    pRefCell        0;
    pRefValue       0;
}


// ************************************************************************* //
Attached Files
File Type: zip test_01.zip (51.1 KB, 3 views)
PSander is offline   Reply With Quote

Reply

Tags
euler, steadystate

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
When can I use Euler equations? CellZone Main CFD Forum 8 April 22, 2020 09:51
hard to convergent with SU2 4.0 for Euler adjoint Xianguu SU2 1 July 14, 2015 18:09
euler euler multiphase modelling mdsvjpk FLUENT 0 June 14, 2011 12:50
Steadystate Euler solver jelmer OpenFOAM Running, Solving & CFD 1 June 19, 2006 09:24
euler model CN FLUENT 0 July 4, 2005 17:35


All times are GMT -4. The time now is 23:27.