Conjugate heat transfer giving floating point exception error

 Register Blogs Members List Search Today's Posts Mark Forums Read

April 1, 2022, 09:05
Conjugate heat transfer giving floating point exception error
#1
Member

Bushra Rasheed
Join Date: Dec 2020
Posts: 97
Rep Power: 5
Hi!

I am trying to solve a simple problem of conjugate heat transfer involving two regions; a solid rod heated by two plates at the end (I have not defined plates as a separate solid region, just fixed the temperature of fluid boundary) and there is forced convection around the rod for which heat transfer coefficient is given. The temperature profile of at rod interface is to be calculated. I have defined externalWallHeatFluxTemperature boundary condition at the rod-air interface.

Code:
```      "fluid_to_.*"
{
type            externalWallHeatFluxTemperature;
mode            coefficient;
Ta              uniform 298.0;
h               uniform 25;
kappaMethod     fluidThermo;
value           \$internalField;
}```
I have removed U and p_rgh boundaries from the solid region and gave an inlet velocity of 0.5 to fluid region. The simulation gives me the following error:

Solving for fluid region fluid
Code:
```#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in ~/OpenFOAM/OpenFOAM-v2106/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#6  Foam::CompressibleTurbulenceModel<Foam::fluidThermo>::nu() const at ??:?
#7  Foam::laminarModels::Stokes<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > >::nuEff() const at ??:?
#8  Foam::linearViscousStress<Foam::laminarModel<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::divDevRhoReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const at ??:?
#9  ? in ~/OpenFOAM/OpenFOAM-v2106/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam
#10  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#11  ? in ~/OpenFOAM/OpenFOAM-v2106/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam
Floating point exception (core dumped)```
I have tried changing fvSchemes and fvSolutions but I think the problem lies in velocity boundary condition. I am attaching my boundary conditions. I have tried using both chtMultiRegionFoam and chtMultiRegionSimpleFoam

Thanks!
Attached Files
 fluid.zip (3.3 KB, 7 views) solid.zip (2.0 KB, 3 views)

 April 5, 2022, 14:57 #2 Senior Member   Join Date: Sep 2013 Posts: 353 Rep Power: 20 Code: ``` inlet { type fixedValue; value uniform ( 0.5 0 0 ); } outlet { type fixedValue; value uniform ( 0 0 0 ); }``` Your are sending something into the domain but nothing can leave? Your outlet is a wall? B_R_Khan likes this.

April 8, 2022, 05:38
#3
Member

Bushra Rasheed
Join Date: Dec 2020
Posts: 97
Rep Power: 5
Quote:
 Originally Posted by Bloerb Code: ``` inlet { type fixedValue; value uniform ( 0.5 0 0 ); } outlet { type fixedValue; value uniform ( 0 0 0 ); }``` Your are sending something into the domain but nothing can leave? Your outlet is a wall?
Thanks for pointing that out... I realized my boundary conditions were not correct; I also had a faulty p_rgh boundary condition