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

Compressible subsonic 2d airfoil investigation

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 22, 2018, 11:59
Default Compressible subsonic 2d airfoil investigation
  #1
New Member
 
Max
Join Date: Feb 2018
Posts: 1
Rep Power: 0
happyFoamer is on a distinguished road
Hi Foamer-Community,

I'm trying to simulate an Airfoil 2D with OpenFoam and run into divergence problems. I started with 100 m/s (AoA=2) and increased the velocity step by step to now 180 m/s. For 200 m/s the simulation crashes after reaching about 800 iterations. The residuals of k and epsilon increase and they are bounded.
I use rhoPimpleFoam and the mesh should be fine. There are no warnings when running checkMesh. As turbulence model kEpsilon is activated. TimeStep is set to 1e-4. I'm running openFoam in a virtualbox.

The domain consists of one disk with one cell thickness. The airfoil is cut out in the center of the disk. So there are 3 boundaries: Airfoil(wall); Farfield(inlet and outlet to produce 360 polars); Symmetry_Plane(front and back side of the disk)

It would be great to get some advices. I tried a lot but nothing succeeded.

Here are some specifications of my case:

Code:
fvSolution:

solvers
{
    "(p)"
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          0.001;
    }

    "(p)Final"
    {
        $p;
        relTol          0;
    }

    "(rho|U|e)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-08;
        relTol          0.001;
    }

    "(rho|U|e)Final"
    {
        $U;
        relTol          0;
    }

    "(k|nTilda|epsilon)"
    {
        solver          GAMG;
        smoother        DILUGaussSeidel;
        tolerance       1e-08;
        relTol          0.001;
    }

    "(k|nTilda|epsilon)Final"
    {
        $U;
        relTol          0;
    }

}

PIMPLE
{
    momentumPredictor   yes;
    transonic           no;
    nOuterCorrectors    50;
    nCorrectors         1;
    nNonOrthogonalCorrectors 0;
    consistent          yes;
    SIMPLErho           yes;

    residualControl
    {
        "(U|k|epsilon)"
        {
            relTol          0;
            tolerance       0.0001;
        }
    }

    turbOnFinalIterOnly no;
}

relaxationFactors
{
    fields
    {
        "p.*"           1;
        "rho.*"         1;
    }
    equations
    {
        "U.*"           0.9;
        "(e|h)"         0.7;
        "(k|epsilon|omega).*" 0.8;
    }
}

thermophysicalProperties:
thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       sutherland;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture
{
    specie
    {
        molWeight   28.97;
    }
    thermodynamics
    {
        Cp          1007;
        Hf          0;
    }
    transport
    {
        As          1.458e-06;
        Ts          110.4;

    }
}

boundaries:

alphat:
dimensions      [1 -1 -1 0 0 0 0];

internalField   uniform 1e-3;

boundaryField
{
    Airfoil
    {
        type            compressible::alphatWallFunction;
        value           $internalField;
    }
    Symmetry
    {
        type            empty;

    }
    Farfield
    {
        type            calculated;
        value           $internalField;
    }
}

epsilon:
dimensions      [0 2 -3 0 0 0 0];

internalField   uniform 200;

boundaryField
{
    Airfoil
    {
        type            epsilonWallFunction;
        value           uniform 200;
    }
    Symmetry
    {
    type        empty;

    }
    Farfield
    {
        type            inletOutlet;
        inletValue      uniform 200;
        value           uniform 200;
    }
}

k:
dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 1;

boundaryField
{
    Airfoil
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
    Symmetry
    {
    type        empty;

    }
    Farfield
    {
        type            inletOutlet;
        inletValue      uniform 1;
        value           uniform 1;
    }
}

p:
dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 1.0e5;

boundaryField
{
    Airfoil
    {
        type            zeroGradient;
    }
    Symmetry
    {
        type            empty;
    }
    Farfield
    {
        type            fixedValue;
        value           $internalField;
    }
}

T: 
dimensions      [0 0 0 1 0 0 0];

internalField   uniform 288.15;

boundaryField
{

    Airfoil
    {
        type            zeroGradient;
    }
    Symmetry
    {
    type        empty;

    }
    Farfield
    {
        type            inletOutlet;
        value           $internalField;
        inletValue      $internalField;
    }
}

U:
dimensions [ 0 1 -1 0 0 0 0 ];

internalField uniform (99.94 3.49 0);

boundaryField
{
  Airfoil
  {
    type noSlip;
  }
  Symmetry_Plane
  {
    type empty;
  }
  Farfield
  {
    type freestream;
    freestreamValue uniform (99.939083 3.489950 0);
  }
}
Error:

Code:
[1] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[1] #1  Foam::sigFpe::sigHandler(int) at ??:?
[1] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
[1] #4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
[1] #5  ? at ??:?
[1] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #7  ? at ??:?
[openfoam-VirtualBox:16350] *** Process received signal ***
[openfoam-VirtualBox:16350] Signal: Floating point exception (8)
[openfoam-VirtualBox:16350] Signal code:  (-6)
[openfoam-VirtualBox:16350] Failing at address: 0x3e800003fde
[openfoam-VirtualBox:16350] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fb0575734b0]
[openfoam-VirtualBox:16350] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fb057573428]
[openfoam-VirtualBox:16350] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fb0575734b0]
[openfoam-VirtualBox:16350] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEE9calculateEv+0x184)[0x7fb05cd58de4]
[openfoam-VirtualBox:16350] [ 4] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEE7correctEv+0x2a)[0x7fb05cd6865a]
[openfoam-VirtualBox:16350] [ 5] rhoPimpleFoam[0x42ab96]
[openfoam-VirtualBox:16350] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fb05755e830]
[openfoam-VirtualBox:16350] [ 7] rhoPimpleFoam[0x4328f9]
[openfoam-VirtualBox:16350] *** End of error message ***
Thanks a lot in advance! :-)
happyFoamer is offline   Reply With Quote

Reply

Tags
airfoil 2d, compressible, rhopimplefoam, transient 2d

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
Ffd_control_point_2d feiyi SU2 4 September 30, 2019 13:42
High drag for airfoil compared to XFOIL and wind tunnel data Ry10 SU2 15 October 30, 2016 18:27
2D FFD Optimization RLangtry SU2 2 August 5, 2014 10:48
SU2 AOA optimization 454514566@qq.com SU2 8 July 7, 2014 08:01
Problem with airfoil shape optimization robyTKD SU2 Shape Design 6 June 13, 2013 20:13


All times are GMT -4. The time now is 11:48.