CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   buoyantBoussinesqPimpleFoam & turbulentHeatFluxTemperature (http://www.cfd-online.com/Forums/openfoam-solving/114418-buoyantboussinesqpimplefoam-turbulentheatfluxtemperature.html)

wolfindark March 11, 2013 01:42

buoyantBoussinesqPimpleFoam & turbulentHeatFluxTemperature
 
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)


jherb April 5, 2013 05:29

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;

wolfindark April 7, 2013 21:22

Thanks jherb,
you were right.


All times are GMT -4. The time now is 16:09.