|
[Sponsors] |
buoyantBoussinesqPimpleFoam & turbulentHeatFluxTemperature |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 11, 2013, 00:42 |
buoyantBoussinesqPimpleFoam & turbulentHeatFluxTemperature
|
#1 |
Member
Neilson Whit
Join Date: Aug 2011
Posts: 74
Rep Power: 14 |
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) |
|
April 5, 2013, 05:29 |
|
#2 |
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 21 |
I guess you made the same mistake as me:
http://www.openfoam.org/mantisbt/view.php?id=806 You used Code:
alphaEff kappat; Code:
alphaEff kappaEff; |
|
April 7, 2013, 21:22 |
|
#3 |
Member
Neilson Whit
Join Date: Aug 2011
Posts: 74
Rep Power: 14 |
Thanks jherb,
you were right. |
|
September 25, 2015, 16:26 |
|
#4 |
New Member
seyyed
Join Date: Jun 2014
Posts: 7
Rep Power: 11 |
hi
im modeling the same problem did this solver work for you? how did you get the result? thanks |
|
|
|