[B_R_Khan]

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.

      "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

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib/x86_64-linux-gnu/libpthread.so.0
#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

Please help me solve this error!


Thanks!

[Bloerb]

    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]

Originally Posted by Bloerb

    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