[wolfindark]

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:

    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:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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)

[jherb]

I guess you made the same mistake as me:
http://www.openfoam.org/mantisbt/view.php?id=806

You used

        alphaEff        kappat;

instead of

        alphaEff        kappaEff;

[wolfindark]

Thanks jherb,
you were right.

[S.M.H]

hi

im modeling the same problem
did this solver work for you?
how did you get the result?

thanks