[ayu12]

I am just trying to simulate the compressible flow of air within two hotplates maintained at a temperature of 1800 K. It is a 2D Simulation. I want to understand how heat transfer affects the Velocity and temperature of Air flowing within the Hot Plates. There are 2 Solid regions and 1 Fluid region. I am getting a floating-point error in the beginning itself.

To resolve the issue, i tried many cases like,

1.) I tried to run a dummy case, in order to understand if the solver is running.

PIMPLE
{
}

PIMPLE algorithm finds no convergences defined for all regions and It shifts to PISO. This is agreeable. And finally a Floating Point error. I updated the convergence criteria as well.

2.)nI rechecked my boundary conditions. For P, it is Calculated and T Zero Gradient. I kept my maxCo as 0.5.

I do understand this error is because of division by 0 but not able to point out the problem.

Can anyone help me understand what can be the issue? I am very new to OpenFoam and CFD and I am stuck at this point and not able to resolve this issue.

My Log File

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 7-1ff648926f77
Exec   : chtMultiRegionFoam
Date   : May 28 2020
Time   : 17:40:26
Host   : "test-VBox"
PID    : 1699
I/O    : uncollated
Case   : /home/test/OpenFOAM/hotpipe2D
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create fluid mesh for region Air for time = 0

Create solid mesh for region solid1 for time = 0

Create solid mesh for region solid2 for time = 0

*** Reading fluid mesh thermophysical properties for region Air

    Adding to thermoFluid

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

    Adding to rhoFluid

    Adding to UFluid

    Adding to phiFluid

    Adding to gFluid

    Adding to hRefFluid

    Adding to ghFluid

    Adding to ghfFluid

    Adding to turbulenceFluid

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
    Adding to reactionFluid

Combustion model not active: combustionProperties not found
Selecting combustion model none
    Adding to radiationFluid

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding to KFluid

    Adding to dpdtFluid

    Adding to fieldsFluid

    Adding MRF

No MRF models present

    Adding fvOptions

*** Reading solid mesh thermophysical properties for region solid1

    Adding to thermos

Selecting thermodynamics package 
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

*** Reading solid mesh thermophysical properties for region solid2

    Adding to thermos

Selecting thermodynamics package 
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions


PIMPLE: Region Air
PIMPLE: Convergence criteria found
        p: tolerance 1e-05
        U: tolerance 1e-05
        T: tolerance 1e-05


PIMPLE: Region solid1
PIMPLE: Convergence criteria found
        T: tolerance 1e-05


PIMPLE: Region solid2
PIMPLE: Convergence criteria found
        T: tolerance 1e-05


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode


#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::operator/(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#5  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam"
#6  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam"
Floating point exception (core dumped)

My case File is attached for reference.

[lonewanderer]

In p_rgh file give zeroGradient at inlet and fixed value at outlet