[OldUnixGuy49]

I am modeling a large piece of hot equipment that will be mounted above a roof to determine temperatures at roofing and waterproofing membranes. I've used Openfoam a little previously, but I'm not expert. I've researched the error I'm having and resolved several problems on my own including having too large a model for my computer. My mesh seems to be fine now, and solution progresses up to Time = 1 where it crashes as follows:

multiRegion multiRegion:
Spawned additional object(s):
multiRegion:Air1
multiRegion:Air2
multiRegion:Air3
multiRegion:Concrete1
multiRegion:Concrete2
multiRegion:Insulation1
multiRegion:Insulation2
multiRegion:Insulation3
multiRegion:Insulation4
multiRegion:Steel1
multiRegion:Steel2

Time = 1


Solving for fluid region Air1
#0 Foam::error:rintStack(Foam::Ostream&) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#1 Foam::sigFpe::sigHandler(int) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#4 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#5 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#6 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#7 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#8 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#9 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam
#10 ? in /lib/x86_64-linux-gnu/libc.so.6
#11 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#12 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam
Floating point exception (core dumped)

I assume this is trying to tell me that some starting value I have is zero and causing a divide by zero error. Further the problem is in TurbulenceModels for Air1, though it may exist elsewhere. I'm currently using laminar in turbulenceProperties, but will likely need to change that if I need to simulate forced air circulation.

I would appreciate any guidance on which inputs to review/modify to get past this roadblock. Here is my changeDictionary for Air1 if that will help:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Website: https://openfoam.org |
| \ / A nd | Version: v2212 & Ned-Foam 2.0 November 2024 |
| \/ M anipulation | Source: AirSpaceCFD02-01
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object changeDictionaryDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

U
{
internalField uniform (0 0 0);

boundaryField
{
min-Y
{
type calculated;
value uniform (0 0 0);
}
max-Y
{
type calculated;
value uniform (0 0 0);
}
min-X
{
type calculated;
value uniform (0 0 0);
}
max-X
{
type calculated;
value uniform (0 0 0);
}
Air1_to_Steel1
{
type calculated;
value uniform (0 0 0);
}
Air1_to_Concrete1
{
type calculated;
value uniform (0 0 0);
}
Air1_to_Steel2
{
type calculated;
value uniform (0 0 0);
}
Air1_to_Insulation2
{
type calculated;
value uniform (0 0 0);
}
Air1_to_Insulation1
{
type calculated;
value uniform (0 0 0);
}
Air1_to_Insulation3
{
type calculated;
value uniform (0 0 0);
}
}
}

T
{
internalField uniform 322.04;

boundaryField
{
min-Y
{
type calculated;
value uniform 322.04;
}
max-Y
{
type calculated;
value uniform 322.04;
}
min-X
{
type calculated;
value uniform 322.04;
}
max-X
{
type calculated;
value uniform 322.04;
}
Air1_to_Steel1
{
type calculated;
value uniform 322.04;
}
Air1_to_Concrete1
{
type calculated;
value uniform 322.04;
}
Air1_to_Steel2
{
type calculated;
value uniform 322.04;
}
Air1_to_Insulation2
{
type calculated;
value uniform 322.04;
}
Air1_to_Insulation1
{
type calculated;
value uniform 322.04;
}
Air1_to_Insulation3
{
type calculated;
value uniform 322.04;
}
}
}

epsilon
{
internalField uniform 0.01;

boundaryField
{
min-Y
{
type calculated;
value uniform 0.01;
}
max-Y
{
type calculated;
value uniform 0.01;
}
min-X
{
type calculated;
value uniform 0.01;
}
max-X
{
type calculated;
value uniform 0.01;
}
Air1_to_Steel1
{
type calculated;
value uniform 0;
}
Air1_to_Concrete1
{
type calculated;
value uniform 0;
}
Air1_to_Steel2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation1
{
type calculated;
value uniform 0;
}
Air1_to_Insulation3
{
type calculated;
value uniform 0;
}
}
}

k
{
internalField uniform 0.1;

boundaryField
{
min-Y
{
type calculated;
value uniform 0.1;
}
max-Y
{
type calculated;
value uniform 0.1;
}
min-X
{
type calculated;
value uniform 0.1;
}
max-X
{
type calculated;
value uniform 0.1;
}
Air1_to_Steel1
{
type calculated;
value uniform 0;
}
Air1_to_Concrete1
{
type calculated;
value uniform 0;
}
Air1_to_Steel2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation1
{
type calculated;
value uniform 0;
}
Air1_to_Insulation3
{
type calculated;
value uniform 0;
}
}
}

p_rgh
{
internalField uniform 101325;

boundaryField
{
min-Y
{
type calculated;
value uniform 101325;
}
max-Y
{
type calculated;
value uniform 101325;
}
min-X
{
type calculated;
value uniform 101325;
}
max-X
{
type calculated;
value uniform 101325;
}
Air1_to_Steel1
{
type calculated;
value uniform 101325;
}
Air1_to_Concrete1
{
type calculated;
value uniform 101325;
}
Air1_to_Steel2
{
type calculated;
value uniform 101325;
}
Air1_to_Insulation2
{
type calculated;
value uniform 101325;
}
Air1_to_Insulation1
{
type calculated;
value uniform 101325;
}
Air1_to_Insulation3
{
type calculated;
value uniform 101325;
}
}
}

p
{
internalField uniform 101325;

boundaryField
{
min-Y
{
type calculated;
value uniform 0;
}
max-Y
{
type calculated;
value uniform 0;
}
min-X
{
type calculated;
value uniform 0;
}
max-X
{
type calculated;
value uniform 0;
}
Air1_to_Steel1
{
type calculated;
value uniform 0;
}
Air1_to_Concrete1
{
type calculated;
value uniform 0;
}
Air1_to_Steel2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation1
{
type calculated;
value uniform 0;
}
Air1_to_Insulation3
{
type calculated;
value uniform 0;
}
}
}

qr
{
internalField uniform 0;

boundaryField
{
min-Y
{
type calculated;
value uniform 0;
}
max-Y
{
type calculated;
value uniform 0;
}
min-X
{
type calculated;
value uniform 0;
}
max-X
{
type calculated;
value uniform 0;
}
Air1_to_Steel1
{
type calculated;
value uniform 0;
}
Air1_to_Concrete1
{
type calculated;
value uniform 0;
}
Air1_to_Steel2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation1
{
type calculated;
value uniform 0;
}
Air1_to_Insulation3
{
type calculated;
value uniform 0;
}
}
}

G
{
internalField uniform 0;

boundaryField
{
min-Y
{
type calculated;
value uniform 0;
}
max-Y
{
type calculated;
value uniform 0;
}
min-X
{
type calculated;
value uniform 0;
}
max-X
{
type calculated;
value uniform 0;
}
Air1_to_Steel1
{
type calculated;
value uniform 0;
}
Air1_to_Concrete1
{
type calculated;
value uniform 0;
}
Air1_to_Steel2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation1
{
type calculated;
value uniform 0;
}
Air1_to_Insulation3
{
type calculated;
value uniform 0;
}
}
}

IDefault
{
internalField uniform 0;

boundaryField
{
min-Y
{
type calculated;
value uniform 0;
}
max-Y
{
type calculated;
value uniform 0;
}
min-X
{
type calculated;
value uniform 0;
}
max-X
{
type calculated;
value uniform 0;
}
Air1_to_Steel1
{
type calculated;
value uniform 0;
}
Air1_to_Concrete1
{
type calculated;
value uniform 0;
}
Air1_to_Steel2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation2
{
type calculated;
value uniform 0;
}
Air1_to_Insulation1
{
type calculated;
value uniform 0;
}
Air1_to_Insulation3
{
type calculated;
value uniform 0;
}
}
}

alphat
{
internalField uniform 0;

boundaryField
{
min-Y
{
type zeroGradient;
}
max-Y
{
type zeroGradient;
}
min-X
{
type zeroGradient;
}
max-X
{
type zeroGradient;
}
Air1_to_Steel1
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
Air1_to_Concrete1
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
Air1_to_Steel2
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
Air1_to_Insulation2
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
Air1_to_Insulation1
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
Air1_to_Insulation3
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
}
}

nut
{
internalField uniform 0;

boundaryField
{
min-Y
{
type calculated;
value uniform 0;
}
max-Y
{
type calculated;
value uniform 0;
}
min-X
{
type calculated;
value uniform 0;
}
max-X
{
type calculated;
value uniform 0;
}
Air1_to_Steel1
{
type nutUSpaldingWallFunction;
value uniform 0;
}
Air1_to_Concrete1
{
type nutUSpaldingWallFunction;
value uniform 0;
}
Air1_to_Steel2
{
type nutUSpaldingWallFunction;
value uniform 0;
}
Air1_to_Insulation2
{
type nutUSpaldingWallFunction;
value uniform 0;
}
Air1_to_Insulation1
{
type nutUSpaldingWallFunction;
value uniform 0;
}
Air1_to_Insulation3
{
type nutUSpaldingWallFunction;
value uniform 0;
}
}
}

// ************************************************** *********************** //


Thanks in advance.

[OldUnixGuy49]

After a great deal of trial and error I was able to get this model off the ground. Basic issues were not having boundaries setup correctly and needing mappedPatch types for internal surface intersections. The model now runs but has other issues which I will address in a separate post.