berkozaydin |
June 22, 2021 07:48 |
buoyantSimpleFoam Error
Hi,
I try to solve a steady, turbulent, incompressible and heat transfer problem. In a rectangular shape, there are 13 batteries. And they have a heat flux value. In all heat transfer models, I have to change the compressible file to incompressible. Therefore, I always make mistakes.
In this study, I used buoyantSimpleFoam. There are some problems while writing the necessary input. And I got this error.
What should I do? Maybe, you can suggest another solver (probably you will say chtMultiRegion) or we can talk about changing compressible file to incompressible (I don't think it's such an easy matter that I can replace density only.)
And also, I want to add thermophysical properties file. Maybe, I have to change something in this file.
Thanks!
Code:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 8
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}
mixture
{
specie
{
molWeight 28.96;
}
thermodynamics
{
Cp 1004.4;
Hf 0;
}
transport
{
mu 1.831e-05;
Pr 0.705;
}
}
// ************************************************************************* //
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
SIMPLE: Convergence criteria found
p_rgh: tolerance 0.001
U: tolerance 0.0001
h: tolerance 0.0001
"(k|epsilon|omega)": tolerance 0.005
Reading thermophysical properties
Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}
Reading field U
Reading/calculating face flux field phi
Creating turbulence model
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
model kEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}
Creating thermophysical transport model
Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model eddyDiffusivity
Reading g
Reading hRef
Calculating field g.h
Reading pRef
Reading field p_rgh
No MRF models present
Radiation model not active: radiationProperties not found
Selecting radiationModel none
No finite volume options present
Starting time loop
Time = 1
DILUPBiCGStab: Solving for h, Initial residual = 1, Final residual = 0.00536049, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.999989, Final residual = 0.0037206, No Iterations 7
time step continuity errors : sum local = 0.269408, global = -0.00834916, cumulative = -0.00834916
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.0321705, Final residual = 7.80071e-06, No Iterations 1
DILUPBiCGStab: Solving for k, Initial residual = 1, Final residual = 0.000295899, No Iterations 1
ExecutionTime = 0.46 s ClockTime = 0 s
Time = 2
DILUPBiCGStab: Solving for h, Initial residual = 0.026843, Final residual = 0.000173081, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.892057, Final residual = 0.00336125, No Iterations 3
time step continuity errors : sum local = 766.721, global = -4.69339, cumulative = -4.70174
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.740648, Final residual = 0.000894831, No Iterations 1
bounding epsilon, min: -335.513 max: 60702.7 average: 189.468
DILUPBiCGStab: Solving for k, Initial residual = 0.983617, Final residual = 0.000904046, No Iterations 1
bounding k, min: -6.35689 max: 628.844 average: 10.4528
ExecutionTime = 0.59 s ClockTime = 0 s
Time = 3
DILUPBiCGStab: Solving for h, Initial residual = 0.0810137, Final residual = 0.000642639, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.931238, Final residual = 0.00480081, No Iterations 5
time step continuity errors : sum local = 929.728, global = -13.246, cumulative = -17.9478
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.854529, Final residual = 0.000398156, No Iterations 1
bounding epsilon, min: -1834.41 max: 381686 average: 1660.7
DILUPBiCGStab: Solving for k, Initial residual = 0.63389, Final residual = 0.000277833, No Iterations 1
bounding k, min: -113.523 max: 7561.63 average: 67.972
ExecutionTime = 0.71 s ClockTime = 1 s
Time = 4
DILUPBiCGStab: Solving for h, Initial residual = 0.47617, Final residual = 0.00421619, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
#8 Foam::fvMatrix<double>::solve() in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
#9 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
Floating point exception (core dumped)
|