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

Film cooling problem

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 26, 2020, 10:08
Default Film cooling problem
  #1
Member
 
Arun subramanian
Join Date: Jun 2016
Location: Florence,Italy
Posts: 48
Rep Power: 10
arun1994 is on a distinguished road
Hello,


I have a film cooling problem of a flat plate with RANS (k-ep). The geometry is simple as in figure attached. I am using rhoSimpleFoam to solve the problem but I have a floating point exception always I have tried to play with the boundary conditions quite a bit but I still cannot find where the problem is.



Code:
[25] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[25] #1  Foam::sigFpe::sigHandler(int) at ??:?
[25] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[25] #3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
[25] #4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
[25] #5  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
[25] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[25] #7  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
[ARUBUNTU:10065] *** Process received signal ***
[ARUBUNTU:10065] Signal: Floating point exception (8)
[ARUBUNTU:10065] Signal code:  (-6)
[ARUBUNTU:10065] Failing at address: 0x3e800002751
[ARUBUNTU:10065] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3efd0)[0x7fc201030fd0]
[ARUBUNTU:10065] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7fc201030f47]
[ARUBUNTU:10065] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3efd0)[0x7fc201030fd0]
[ARUBUNTU:10065] [ 3] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEE9calculateEv+0x1e2)[0x7fc206c0eb02]
[ARUBUNTU:10065] [ 4] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEE7correctEv+0x2a)[0x7fc206c4458a]
[ARUBUNTU:10065] [ 5] rhoSimpleFoam(+0x2889f)[0x562fd890489f]
[ARUBUNTU:10065] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7fc201013b97]
[ARUBUNTU:10065] [ 7] rhoSimpleFoam(+0x2d1aa)[0x562fd89091aa]
[ARUBUNTU:10065] *** End of error message ***
My 0 files look like this :
----------------------------------------------------------------------------------

alphat



Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    Body.outlet
    {
        type            calculated;
        value           uniform 0;
    }
    Body.wall_base
    {
        type            compressible::alphatWallFunction;
        Prt             0.85;
        value           uniform 0;
    }
    Body.wall_side2
    {
        type            compressible::alphatWallFunction;
        Prt             0.85;
        value           uniform 0;
    }
    Body.topwall
    {
        type            compressible::alphatWallFunction;
        Prt             0.85;
        value           uniform 0;
    }
    Body.wall_side
    {
        type            compressible::alphatWallFunction;
        Prt             0.85;
        value           uniform 0;
    }
    Body.inlet_coolant
    {
        type            calculated;
        value           uniform 0;
    }
    Body.inlet
    {
        type            calculated;
        value           uniform 0;
    }
    Body.wall_coolant
    {
        type            compressible::alphatWallFunction;
        Prt             0.85;
        value           uniform 0;
    }
}
----------------------------------------------------------------------------------------------------
velocity U



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

internalField   uniform (0 0 0);

boundaryField
{
    Body.outlet
    {
        type            inletOutlet;
        value           uniform (0 0 0);
        inletValue      uniform (0 0 0);
    }
    Body.wall_base
    {
        type            noSlip;
    }
    Body.wall_side2
    {
        type            noSlip;
    }
    Body.topwall
    {
        type            noSlip;
    }
    Body.wall_side
    {
        type            noSlip;
    }
    Body.inlet_coolant
    {
    type surfaceNormalFixedValue;
    refValue uniform -10;
        rhoInlet        1.177;    // Guess for rho
    }
    Body.inlet
    {
    type surfaceNormalFixedValue;
    refValue uniform -10;
        rhoInlet        1;    // Guess for rho
    }
    Body.wall_coolant
    {
        type            noSlip;
    }
}
----------------------------------------------------------------------------------------------





Temperature


Code:
dimensions      [0 0 0 1 0 0 0];

internalField   uniform 300;

boundaryField
{
    Body.outlet
    {

        type            inletOutlet;
        value           $internalField;
        inletValue      $internalField;
    }
    Body.wall_base
    {
        type            zeroGradient;
    }
    Body.wall_side2
    {
        type            zeroGradient;
    }
    Body.topwall
    {
        type            zeroGradient;
    }
    Body.wall_side
    {
        type            zeroGradient;
    }
    Body.inlet_coolant
    {
        type            fixedValue;
        value           uniform 313;
    }
    Body.inlet
    {
        type            fixedValue;
        value           uniform 293;
    }
    Body.wall_coolant
    {
        type            zeroGradient;
    }
}
----------------------------------------------------------------------------------------------
Pressure



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

internalField   uniform 101325;

boundaryField
{
    Body.outlet
    {
        type            fixedValue;
        value           uniform 101325;

        // type            mixed;
        // refValue        uniform 110000;
        // refGradient     uniform 0;
        // valueFraction   uniform 1;
        // type            transonicOutletPressure;
        // U               U;
        // phi             phi;
        // gamma           1.4;
        // psi             psi;
        // pInf            uniform 110000;
    }
    Body.wall_base
    {
        type            zeroGradient;
    }
    Body.wall_side2
    {
        type            zeroGradient;
    }
    Body.topwall
    {
        type            zeroGradient;
    }
    Body.wall_side
    {
        type            zeroGradient;
    }
    Body.inlet_coolant
    {
        type            zeroGradient;
        // type            mixed;
        refValue        uniform 101325;
        refGradient     uniform 0;
        valueFraction   uniform 0.3;
    }
    Body.inlet
    {
        type            zeroGradient;
        // type            mixed;
        refValue        uniform 101325;
        refGradient     uniform 0;
        valueFraction   uniform 0.3;
    }
    Body.wall_coolant
    {
        type            zeroGradient;
    }
}





----------------------------------------------------------------------------------------------
nut
Code:
dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    Body.outlet
    {
        type            calculated;
        value           uniform 0;
    }
    Body.wall_base
    {
        type            nutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0;
    }
    Body.wall_side2
    {
        type            nutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0;
    }
    Body.topwall
    {
        type            nutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0;
    }
    Body.wall_side
    {
        type            nutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0;
    }
    Body.inlet_coolant
    {
        type            calculated;
        value           uniform 0;
    }
    Body.inlet
    {
        type            calculated;
        value           uniform 0;
    }
    Body.wall_coolant
    {
        type            nutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0;
    }
}

----------------------------------------------------------------------------------------------


k


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

internalField   uniform 1;

boundaryField
{
    Body.outlet
    {
        type            inletOutlet;
        inletValue      uniform 1;
        value           uniform 1;
    }
    Body.wall_base
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
    Body.wall_side2
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
    Body.topwall
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
    Body.wall_side
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
    Body.inlet_coolant
    {
        type            turbulentIntensityKineticEnergyInlet;
        intensity       0.05;
        value           uniform 1;
    }
    Body.inlet
    {
        type            turbulentIntensityKineticEnergyInlet;
        intensity       0.05;
        value           uniform 1;
    }
    Body.wall_coolant
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
}



----------------------------------------------------------------------------------------------
epsilon





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

internalField   uniform 200;

boundaryField
{
    Body.outlet
    {
        type            inletOutlet;
        inletValue      uniform 200;
        value           uniform 200;
    }
    Body.wall_base
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 200;
    }
    Body.wall_side2
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 200;
    }
    Body.topwall
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 200;
    }
    Body.wall_side
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 200;
    }
    Body.inlet_coolant
    {
        type            turbulentMixingLengthDissipationRateInlet;
        mixingLength    0.005;
        value           uniform 200;
    }
    Body.inlet
    {
        type            turbulentMixingLengthDissipationRateInlet;
        mixingLength    0.005;
        value           uniform 200;
    }
    Body.wall_coolant
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 200;
    }
}
----------------------------------------------------------------------------------------------
I am posting the BCs because I think this is where the problem should be.



If I am missing some details let me know please

Thanks for your time!
Arun
Attached Images
File Type: png CFDonline.png (34.2 KB, 12 views)
arun1994 is offline   Reply With Quote

Reply

Tags
film cooling, openfaom-7, rhosimplefoam

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
[interCondensatingEvaporatingFoam] Film boiling - mixing problem bouhoura OpenFOAM Running, Solving & CFD 11 January 7, 2023 16:13
SU2-7.0.1 on ubuntu 18.04 hyunko SU2 Installation 7 March 16, 2020 05:37
BuoyantBoussinesqSimpleFoam_Facing problem Mondal131211 OpenFOAM Running, Solving & CFD 1 April 10, 2019 20:41
Problem on Fluent Tutorial: Horizontal Film Boilig Feng FLUENT 2 April 13, 2013 06:34
Supersonic film cooling Solver selection Balasubramanian.S FLUENT 1 May 5, 2006 11:20


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