CFD Online Discussion Forums - Error in running case in chtMultiRegionFoam

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)

- - Error in running case in chtMultiRegionFoam (https://www.cfd-online.com/Forums/openfoam-programming-development/153018-error-running-case-chtmultiregionfoam.html)

Error in running case in chtMultiRegionFoam

Hi,
Greeting all!!!
I want to use the chtMultiRegionFoam for a geometry of air surrounding with solid.
When i was trying to run this case I got the following error .....

Quote:

Create time

Create fluid mesh for region fluid_air for time = 0

Create solid mesh for region solid_insu for time = 0

*** Reading fluid mesh thermophysical properties for region fluid_air

Adding to thermoFluid

Selecting thermodynamics package hRhoThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Adding to rhoFluid

Adding to kappaFluid

Adding to UFluid

Adding to phiFluid

Adding to gFluid

Adding to turbulence

Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
alphaK1 0.85034;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.85616;
Prt 1;
gamma1 0.5532;
gamma2 0.4403;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
c1 10;
}

Adding to ghFluid

Adding to ghfFluid

Selecting radiationModel viewFactor
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting scatterModel constantScatter
Total number of clusters : 115
Insert elements in the matrix...
Smoothing the matrix...
Adding to KFluid

Adding to dpdtFluid

*** Reading solid mesh thermophysical properties for region solid_insu

Adding to thermos

Constructed constSolidThermo with
rho : rho [1 -3 0 0 0 0 0] 8000
Cp : Cp [0 2 -2 -1 0 0 0] 450
K : K [1 1 -3 -1 0 0 0] 80
Hf : Hf [0 2 -2 0 0 0 0] 0
emissivity : emissivity [0 0 0 0 0 0 0] 0.7
kappa : kappa [0 -1 0 0 0 0 0] 0
sigmaS : sigmaS [0 -1 0 0 0 0 0] 0

Region: fluid_air Courant Number mean: 1.61574e-05 max: 8e-05
Region: solid_insu Diffusion Number mean: 4.7133365e-07 max: 1.7777782e-06
Region: fluid_air Courant Number mean: 1.61574e-05 max: 8e-05
Region: solid_insu Diffusion Number mean: 4.7133365e-07 max: 1.7777782e-06
Time = 0.0001

Solving for fluid region fluid_air
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2
at sigaction.c:0
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::operator/(Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::diagonalSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fvMatrix<double>::solve() in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/chtMultiRegionFoam"
#8
in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/chtMultiRegionFoam"
#9 __libc_start_main in "/lib64/libc.so.6"
#10
in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/chtMultiRegionFoam"

Can anyone kindly mention exactly where I have done the mistake??

Thanks in advance
regards,
baran

Problem in high temperature boundary condition in chtMultiregionFOAM

Hello ,
greeting all!!!
I have some simple geometry as mentioned above.Now i want to implement this in openfoam-2.4.0.
Some details are as follows...
geometry- cavity, surrounding this cavity solid, some heating element in cavity
turbulence model- k epsilon
radiation model- viewfactor
boundary condition- constant T (1. 400 K ,2. 600 K)
In the first case of 400 K it works well,but in 600K or higher temperature it diverges after one iteration.....
Can some one please clarify where is the exact problem in my case?

For better understaanding of the problem...I have attached my case file...so please clarify how I resolve my problem??
I have really facing trouble in this...
please help...

https://www.dropbox.com/s/egyyc1vt7tqq1a1/box.rar?dl=0

Regards,
baran

Quote:

Originally Posted by baran_foam (Post 546397)

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2
at sigaction.c:0
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::operator/(Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::diagonalSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

It appears as if you divided something by zero... this could also depend on your mesh (at least I saw this kind of error with mine). Did you run a checkMesh to see if you have issues with some cells?

PS:use [/CODE] tags instead of [/QUOTE]

Greeting all!
@cocone...In the mesh there is no problem...I have some initial condition issue mostly and it create high Courant no ..but I have sort it out...

Regards,
baran