CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Conjugate heat transfer giving floating point exception error (https://www.cfd-online.com/Forums/openfoam-solving/242030-conjugate-heat-transfer-giving-floating-point-exception-error.html)

B_R_Khan April 1, 2022 09:05

Conjugate heat transfer giving floating point exception error
 
2 Attachment(s)
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 ??:?
#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 April 5, 2022 14:57

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 April 8, 2022 05:38

Quote:

Originally Posted by Bloerb (Post 825568)
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


All times are GMT -4. The time now is 05:06.