[thomas030]

Hello,

I try to model a block with two fluid regions and one solid heat source as test case. I ran the multiRegionHeaterRadiation tutorial and it works. Now I try to understand how the solver works and for that I rebuild the described testcase. I work with version 2.1.1.

In the first Time step I get an error message for my solid region. The setup for the solid is the same as is the tutorial case.

Time = 1


Solving for fluid region volume_1
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.02505913, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.03564643, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.0262261, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.05032147, No Iterations 2

Decomposing C matrix...

LU Back substitute C matrix..
Min/max T:299.8731 300.0793
GAMG: Solving for p_rgh, Initial residual = 0.9244711, Final residual = 0.004468303, No Iterations 5
time step continuity errors : sum local = 0.6615856, global = -4.182891e-17, cumulative = -4.182891e-17
Min/max rho:1.158317 1.159113

Solving for fluid region volume_2
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.0254788, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.03565059, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.02496372, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.05084451, No Iterations 2

Decomposing C matrix...

LU Back substitute C matrix..
Min/max T:299.8266 300.0985
GAMG: Solving for p_rgh, Initial residual = 0.8836019, Final residual = 0.003833346, No Iterations 5
time step continuity errors : sum local = 0.605208, global = 4.839946e-18, cumulative = -3.698896e-17
Min/max rho:1.158242 1.159293

Solving for solid region heater


--> FOAM FATAL ERROR:
Attempt to cast type zeroGradient to type compressible::turbulentTemperatureRadCoupledMixed

From function refCast<To>(From&)
in file /usr/local/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude/typeInfo.H at line 114.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/lib/libOpenFOAM.so"
#2 Foam::Ostream& Foam:perator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/bin/chtMultiRegionSimpleFoam"
#3 Foam::compressible::turbulentTemperatureRadCoupled MixedFvPatchScalarField::updateCoeffs() in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/lib/libcompressibleTurbulenceModel.so"
#4 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/bin/chtMultiRegionSimpleFoam"
#5 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricFi eld<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/lib/libfiniteVolume.so"
#6 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/lib/libfiniteVolume.so"
#7 Foam::fv::laplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/lib/libfiniteVolume.so"
#8 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&, Foam::word const&) in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/bin/chtMultiRegionSimpleFoam"
#9 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 "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/bin/chtMultiRegionSimpleFoam"
#10 main in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/bin/chtMultiRegionSimpleFoam"
#11 __libc_start_main in "/lib64/libc.so.6"
#12 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/usr/local/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64Gcc44DPOpt/bin/chtMultiRegionSimpleFoam"

Has anybody an idea whats the problem here?

Thank you for your help.

Regards

Thomas