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

buoyantBoussinesqPimpleFoam & turbulentHeatFluxTemperature

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 11, 2013, 01:42
Question buoyantBoussinesqPimpleFoam & turbulentHeatFluxTemperature
  #1
Member
 
Neilson Whit
Join Date: Aug 2011
Location: USA
Posts: 31
Rep Power: 5
wolfindark is on a distinguished road
Dear Foamers

I try to simulate natural convection and temporal variation of temperature distribution inside a water pool.
with:
OpenFOAM version: 2.1.1.
Solver : buoyantBoussinesqPimpleFoam
RAS : kEps turbulent

I use turbulentHeatFluxTemperature for a constant heat input from a surface with following block in 0/T file:

Code:
    PIPEOUT
    {
           type            turbulentHeatFluxTemperature;
           heatSource      flux;        // power [W]; flux [W/m2]
            q               uniform 10;  // heat power or flux
            alphaEff        kappat;    // alphaEff field name;
                                         // alphaEff in [kg/m/s]
            Cp              Cp;          // Cp field name; Cp in [J/kg/K]
            value           uniform 300; // initial temperature value

    }

Code runs but it gives following error in the second time step. I would appreciate if you could give me an idea where the error is originated?

Error:

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : buoyantBoussinesqPimpleFoam
Date   : Mar 11 2013
Time   : 14:37:12
Host   : "ubun"
PID    : 8648
Case   : /home/neo/OpenFOAM/neo-2.1.1/run/exp_tr
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Reading thermophysical properties

Reading field T

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Creating turbulence model

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    sigmaEps        1.3;
}

Reading field kappat

Calculating field g.h

Courant Number mean: 0 max: 0

PIMPLE: Operating solver in PISO mode


Starting time loop

Time = 1

Courant Number mean: 0 max: 0
DILUPBiCG:  Solving for T, Initial residual = 1, Final residual = 9.88958e-07, No Iterations 208
DICPCG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.0093319, No Iterations 160
time step continuity errors : sum local = 1.97326e-09, global = 1.3291e-21, cumulative = 1.3291e-21
DICPCG:  Solving for p_rgh, Initial residual = 0.00686217, Final residual = 9.33739e-09, No Iterations 212
time step continuity errors : sum local = 4.3591e-14, global = 3.96046e-14, cumulative = 3.96046e-14
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0714385, Final residual = 8.4943e-07, No Iterations 50
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 9.70453e-07, No Iterations 100
ExecutionTime = 5.43 s  ClockTime = 5 s

Time = 2

Courant Number mean: 5.30969e-08 max: 3.99768e-07
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  Foam::operator/(Foam::UList<double> const&, Foam::tmp<Foam::Field<double> > const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5  Foam::incompressible::turbulentHeatFluxTemperatureFvPatchScalarField::updateCoeffs() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#6   at gaussLaplacianSchemes.C:0
#7  Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8  Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#9  Foam::fv::laplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#10  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantBoussinesqPimpleFoam"
#11  Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantBoussinesqPimpleFoam"
#12  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantBoussinesqPimpleFoam"
#13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantBoussinesqPimpleFoam"
 Floating point exception (core dumped)
wolfindark is offline   Reply With Quote

Old   April 5, 2013, 05:29
Default
  #2
Senior Member
 
Joachim Herb
Join Date: Sep 2010
Posts: 206
Rep Power: 9
jherb is on a distinguished road
I guess you made the same mistake as me:
http://www.openfoam.org/mantisbt/view.php?id=806

You used
Code:
        alphaEff        kappat;
instead of
Code:
        alphaEff        kappaEff;
jherb is offline   Reply With Quote

Old   April 7, 2013, 21:22
Default
  #3
Member
 
Neilson Whit
Join Date: Aug 2011
Location: USA
Posts: 31
Rep Power: 5
wolfindark is on a distinguished road
Thanks jherb,
you were right.
wolfindark is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 21:38.