[anaduc]

Hi!

I am having the following error (see below) when running my case with simpleFoam. I am trying to simulate the wind flow around a 3D object. The wind is uniform and unidirectional. So it should be a pretty simple case.

From what I understood by checking online, this floating point exception means that there is an arithmetic problem, probably a division by 0. Although I don't understand how to correct this mistake.

I tried checking my BCs and ICs as mentioned online, but they seem correct to me. I think that either I didn't input my IC and BC correctly, or the mistake might come from elsewhere.

I attached under the error message the U, p, omega, nut and k files.

Thank you for your help!


[8] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[8] #1 Foam::sigFpe::sigHandler(int) at ??:?
[8] #2 ? in "/usr/lib/x86_64-linux-gnu/libc.so.6"
[8] #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
[8] #4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[8] #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[8] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[8] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[8] #8 Foam::fvMatrix<double>::solve() at ??:?
[8] #9 Foam::solvers::incompressibleFluid::correctPressur e() at ??:?
[8] #10 Foam::solvers::incompressibleFluid:ressureCorrec tor() at ??:?
[8] #11 ? in "/opt/openfoam12/platforms/linux64GccDPInt32Opt/bin/foamRun"
[8] #12 ? in "/usr/lib/x86_64-linux-gnu/libc.so.6"
[8] #13 __libc_start_main in "/usr/lib/x86_64-linux-gnu/libc.so.6"
[8] #14 ? in "/opt/openfoam12/platforms/linux64GccDPInt32Opt/bin/foamRun"
[HKPW09ZCHX:15801] *** Process received signal ***
[HKPW09ZCHX:15801] Signal: Floating point exception (8)
[HKPW09ZCHX:15801] Signal code: (-6)
[HKPW09ZCHX:15801] Failing at address: 0x3e800003db9
[HKPW09ZCHX:15801] [ 0] /usr/lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7f9360924520]
[HKPW09ZCHX:15801] [ 1] /usr/lib/x86_64-linux-gnu/libc.so.6(pthread_kill+0x12c)[0x7f93609789fc]
[HKPW09ZCHX:15801] [ 2] /usr/lib/x86_64-linux-gnu/libc.so.6(raise+0x16)[0x7f9360924476]
[HKPW09ZCHX:15801] [ 3] /usr/lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7f9360924520]
[HKPW09ZCHX:15801] [ 4] /opt/openfoam12/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5Fi eldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEE RKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0x 230)[0x7f9361433f50]
[HKPW09ZCHX:15801] [ 5] /opt/openfoam12/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7 PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS 8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x825)[0x7f9361437875]
[HKPW09ZCHX:15801] [ 6] /opt/openfoam12/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5Fi eldIdEERKS2_h+0x4c8)[0x7f9361439bb8]
[HKPW09ZCHX:15801] [ 7] /opt/openfoam12/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegr egatedERKNS_10dictionaryE+0x26b)[0x7f93627cdbeb]
[HKPW09ZCHX:15801] [ 8] /opt/openfoam12/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveERKNS _10dictionaryE+0x23b)[0x7f936229830b]
[HKPW09ZCHX:15801] [ 9] /opt/openfoam12/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveEv+0x 8d)[0x7f936229860d]
[HKPW09ZCHX:15801] [10] /opt/openfoam12/platforms/linux64GccDPInt32Opt/lib/libincompressibleFluid.so(_ZN4Foam7solvers19incomp ressibleFluid15correctPressureEv+0x46b)[0x7f935c105f7b]
[HKPW09ZCHX:15801] [11] /opt/openfoam12/platforms/linux64GccDPInt32Opt/lib/libincompressibleFluid.so(_ZN4Foam7solvers19incomp ressibleFluid17pressureCorrectorEv+0x2c)[0x7f935c12bcec]
[HKPW09ZCHX:15801] [12] foamRun(+0x3ff6)[0x55750e07fff6]
[HKPW09ZCHX:15801] [13] /usr/lib/x86_64-linux-gnu/libc.so.6(+0x29d90)[0x7f936090bd90]
[HKPW09ZCHX:15801] [14] /usr/lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0x80)[0x7f936090be40]
[HKPW09ZCHX:15801] [15] foamRun(+0x4795)[0x55750e080795]
[HKPW09ZCHX:15801] *** End of error message ***
--------------------------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpiexec noticed that process rank 8 with PID 0 on node HKPW09ZCHX exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------


U:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include "include/initialConditions"

dimensions [0 1 -1 0 0 0 0];

internalField uniform $flowVelocity;

boundaryField
{
//- Set patchGroups for constraint patches
// #includeEtc "caseDicts/setConstraintTypes"

#include "include/fixedInlet"

inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type zeroGradient;
}

lowerWall
{
type noSlip;
}

walls
{
type noSlip;
}



#include "include/frontBackUpperPatches"
}


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


p:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include "include/initialConditions"

dimensions [0 2 -2 0 0 0 0];

internalField uniform $pressure;

boundaryField
{
//- Set patchGroups for constraint patches
// #includeEtc "caseDicts/setConstraintTypes"

inlet
{
type zeroGradient;
}

outlet
{
type fixedValue;
value $internalField;
}

lowerWall
{
type zeroGradient;
}

walls
{
type zeroGradient;
}

#include "include/frontBackUpperPatches"
}

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



omega:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include "include/initialConditions"

dimensions [0 0 -1 0 0 0 0];

internalField uniform $turbulentOmega;

boundaryField
{
//- Set patchGroups for constraint patches
// #includeEtc "caseDicts/setConstraintTypes"

#include "include/fixedInlet"

outlet
{
type inletOutlet;
inletValue $internalField;
value $internalField;
}

lowerWall
{
type omegaWallFunction;
value $internalField;
}

walls
{
type omegaWallFunction;
value $internalField;
}

#include "include/frontBackUpperPatches"
}


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

nut :
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
//- Set patchGroups for constraint patches
// #includeEtc "caseDicts/setConstraintTypes"

frontAndBack
{
type calculated;
value uniform 0;
}

inlet
{
type calculated;
value uniform 0;
}

outlet
{
type calculated;
value uniform 0;
}

lowerWall
{
type nutkWallFunction;
value uniform 0;
}

upperWall
{
type calculated;
value uniform 0;
}

walls
{
type nutkWallFunction;
value uniform 0;
}
}


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


k:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include "include/initialConditions"

dimensions [0 2 -2 0 0 0 0];

internalField uniform $turbulentKE;

boundaryField
{
//- Set patchGroups for constraint patches
// #includeEtc "caseDicts/setConstraintTypes"

//- Define inlet conditions
#include "include/fixedInlet"

outlet
{
type inletOutlet;
inletValue $internalField;
value $internalField;
}

lowerWall
{
type kqRWallFunction;
value $internalField;
}

walls
{
type kqRWallFunction;
value $internalField;
}

#include "include/frontBackUpperPatches"
}


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

[AtoHM]

Is that the first iteration or were some iterations completed before the error occured? Best to post a complete log. Error states a Floating Point Exception, so probably it diverged. Have you checked your mesh thoroughly?

Have you defined BCs according to one of the tutorial cases? You can try pressureInletOutletVelocity for U @ Outlet.

[anaduc]

Hi AtoHM, thank you for your help.

You can find below a log file that I outputted when running simpleFoam. Some iterations were completed. (I took out some parts of the logfile and replaced them with "..." as my message was too long. Let me know if I should send the logfile again, in a seperate message)

My BC are defined according to a university professor, but I will try using the tutorial you mentioned and see if I have different results.

I used checkMesh to check my mesh (output is : mesh ok), from my understanding it is not enough, but I also used a model from a professor which leads me to think the problem shouldn't arise from there. Attached you can find a picture of the meh preview (ParaView). I also attached the blockMesh file and the snappyHexMesh file.


Log file :

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 12
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 12-c59e5ce98388
Exec : foamRun -solver incompressibleFluid -parallel
Date : Jul 24 2024
Time : 10:28:22
Host : "HKPW09ZCHX"
PID : 15786
I/O : uncollated
Case : /home/anaisducrot/canopy-short1box
nProcs : 10
Slaves :
9
(
"HKPW09ZCHX.15787"
"HKPW09ZCHX.15789"
"HKPW09ZCHX.15792"
"HKPW09ZCHX.15794"
"HKPW09ZCHX.15796"
"HKPW09ZCHX.15797"
"HKPW09ZCHX.15800"
"HKPW09ZCHX.15801"
"HKPW09ZCHX.15804"
)

Pstream initialised with:
floatTransfer : false
nProcsSimpleSum : 0
commsType : nonBlocking
polling iterations : 0
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 mesh for time = 0

Selecting solver incompressibleFluid
Selecting viscosity model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
RAS
{
RASModel kOmegaSST;
turbulence on;
printCoeffs on;
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.555556;
gamma2 0.44;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}

No MRF models present

Selecting patchDistMethod meshWave

PIMPLE: No convergence criteria found


PIMPLE: Corrector convergence criteria found
p: tolerance 5e-05, relTol 0
k: tolerance 5e-05, relTol 0
omega: tolerance 5e-05, relTol 0
U: tolerance 5e-06, relTol 0
Calculations will do 10 corrections if the convergence criteria are not met


PIMPLE: Operating solver in steady-state mode with 10 outer correctors



Starting time loop

forcesBase forces1:
Not including porosity effects
forces forces1:
Not including porosity effects
wallShearStress wallShearStress:
processing all wall patches

forces forces1 write:
sum of forces:
pressure : (0 0 0)
viscous : (-1.20607e-07 106.164 1.55307e-05)
porous : (0 0 0)
sum of moments:
pressure : (0 0 0)
viscous : (-530.882 -5.44985e-06 -0.0171853)
porous : (0 0 0)

Time = 1s

PIMPLE: Iteration 1
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0792024, No Iterations 10
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0999848, No Iterations 5
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0858495, No Iterations 12
GAMG: Solving for p, Initial residual = 1, Final residual = 0.00976087, No Iterations 6
GAMG: Solving for p, Initial residual = 0.191862, Final residual = 0.00127844, No Iterations 3
time step continuity errors : sum local = 6.30086e-05, global = -8.59509e-06
GAMG: Solving for p, Initial residual = 0.158256, Final residual = 0.0014351, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0299484, Final residual = 0.000207886, No Iterations 4
time step continuity errors : sum local = 1.13169e-05, global = -5.49845e-07
PIMPLE: Iteration 2
smoothSolver: Solving for Ux, Initial residual = 0.427448, Final residual = 0.0307296, No Iterations 7
smoothSolver: Solving for Uy, Initial residual = 0.112374, Final residual = 0.00840946, No Iterations 8
smoothSolver: Solving for Uz, Initial residual = 0.418704, Final residual = 0.0326524, No Iterations 8
GAMG: Solving for p, Initial residual = 0.432929, Final residual = 0.00199944, No Iterations 5
GAMG: Solving for p, Initial residual = 0.0124887, Final residual = 8.0103e-05, No Iterations 3
time step continuity errors : sum local = 4.24614e-05, global = 2.89485e-06
GAMG: Solving for p, Initial residual = 0.0194392, Final residual = 0.00010519, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00205136, Final residual = 9.41258e-06, No Iterations 4
time step continuity errors : sum local = 5.22223e-06, global = -3.85455e-07
...
time step continuity errors : sum local = 3.21654e-06, global = 6.92728e-07
PIMPLE: Iteration 5
smoothSolver: Solving for Ux, Initial residual = 0.143607, Final residual = 0.0118188, No Iterations 6
smoothSolver: Solving for Uy, Initial residual = 0.022934, Final residual = 0.00176465, No Iterations 7
smoothSolver: Solving for Uz, Initial residual = 0.171491, Final residual = 0.0157175, No Iterations 6
GAMG: Solving for p, Initial residual = 0.0142665, Final residual = 8.30776e-05, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00209389, Final residual = 1.47656e-05, No Iterations 4
time step continuity errors : sum local = 1.05518e-05, global = 1.71176e-06
GAMG: Solving for p, Initial residual = 0.00691267, Final residual = 2.98678e-05, No Iterations 4
GAMG: Solving for p, Initial residual = 0.000967677, Final residual = 5.48831e-06, No Iterations 5
time step continuity errors : sum local = 3.90245e-06, global = 1.12953e-06
PIMPLE: Iteration 6
smoothSolver: Solving for Ux, Initial residual = 0.171827, Final residual = 0.0131641, No Iterations 5
...
time step continuity errors : sum local = 0.000228021, global = 4.32945e-05
GAMG: Solving for p, Initial residual = 0.00975068, Final residual = 4.50914e-05, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00653201, Final residual = 6.28742e-05, No Iterations 4
time step continuity errors : sum local = 0.000217799, global = 2.97954e-05, cumulative = 2.97954e-05
smoothSolver: Solving for omega, Initial residual = 0.111799, Final residual = 0.0106602, No Iterations 2
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0832829, No Iterations 2
PIMPLE: Not converged within 10 iterations
ExecutionTime = 3.40957 s ClockTime = 4 s

forces forces1 write:
sum of forces:
pressure : (-10216.1 9057.88 -171950)
viscous : (-1.24085 110.425 -4.58193)
porous : (0 0 0)
sum of moments:
pressure : (-286215 -64215.2 16995.8)
viscous : (-557.082 -6.08914 0.352403)
porous : (0 0 0)

Time = 2s

PIMPLE: Iteration 1
smoothSolver: Solving for Ux, Initial residual = 0.704672, Final residual = 0.0445682, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.352361, Final residual = 0.0275593, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.786524, Final residual = 0.0498515, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0673504, Final residual = 0.000615409, No Iterations 12
GAMG: Solving for p, Initial residual = 0.00879416, Final residual = 7.91831e-05, No Iterations 8
time step continuity errors : sum local = 0.00145827, global = 4.32979e-05
GAMG: Solving for p, Initial residual = 0.0119511, Final residual = 9.71019e-05, No
...
Iterations 4
GAMG: Solving for p, Initial residual = 0.0166182, Final residual = 0.000161263, No Iterations 8
time step continuity errors : sum local = 0.00259233, global = 0.000436285
PIMPLE: Iteration 7
smoothSolver: Solving for Ux, Initial residual = 0.664621, Final residual = 0.0471121, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.51013, Final residual = 0.0303934, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.73645, Final residual = 0.0541652, No Iterations 4
GAMG: Solving for p, Initial residual = 0.114971, Final residual = 0.00090848, No Iterations 6
GAMG: Solving for p, Initial residual = 0.00116456, Final residual = 1.12998e-05, No Iterations 9
time step continuity errors : sum local = 0.00750699, global = 0.00134886
GAMG: Solving for p, Initial residual = 0.00418543, Final residual = 3.58571e-05, No Iterations 5
GAMG: Solving for p, Initial residual = 0.00284843, Final residual = 2.68966e-05, No Iterations 9
time step continuity errors : sum local = 0.00409727, global = 0.000518198
PIMPLE: Iteration 8
smoothSolver: Solving for Ux, Initial residual = 0.684966, Final residual = 0.0463303, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.539349, Final residual = 0.0392157, No Iterations 5
smoothSolver: Solving for Uz, Initial residual = 0.786623, Final residual = 0.0565337, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0794518, Final residual = 0.000696336, No Iterations 7
GAMG: Solving for p, Initial residual = 0.000689193, Final residual = 6.35372e-06, No Iterations 12
time step continuity errors : sum local = 0.0085613, global = 0.00205354
GAMG: Solving for p, Initial residual = 0.00320712, Final residual = 2.39184e-05, No Iterations 5
GAMG: Solving for p, Initial residual = 0.00158733, Final residual = 1.5083e-05, No Iterations 11
time step continuity errors : sum local = 0.00526274, global = 0.000998422
PIMPLE: Iteration 9
smoothSolver: Solving for Ux, Initial residual = 0.671852, Final residual = 0.0462837, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.558497, Final residual = 0.0386459, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.741289, Final residual = 0.0528999, No Iterations 4
GAMG: Solving for p, Initial residual = 0.043713, Final residual = 0.000337943, No Iterations 7
GAMG: Solving for p, Initial residual = 0.000565967, Final residual = 5.64039e-06, No Iterations 11
time step continuity errors : sum local = 0.0152731, global = 0.00343629
GAMG: Solving for p, Initial residual = 0.00241708, Final residual = 1.59423e-05, No Iterations 5
GAMG: Solving for p, Initial residual = 0.00107662, Final residual = 9.48259e-06, No Iterations 9
time step continuity errors : sum local = 0.00748813, global = 0.00167025
PIMPLE: Iteration 10
smoothSolver: Solving for Ux, Initial residual = 0.683772, Final residual = 0.0500167, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.537289, Final residual = 0.0390661, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.766006, Final residual = 0.0524462, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0304014, Final residual = 0.000296393, No Iterations 6
GAMG: Solving for p, Initial residual = 0.000367517, Final residual = 3.66684e-06, No Iterations 13
time step continuity errors : sum local = 0.0209913, global = 0.00649364
GAMG: Solving for p, Initial residual = 0.00175619, Final residual = 1.14259e-05, No Iterations 5
GAMG: Solving for p, Initial residual = 0.000873487, Final residual = 8.72338e-06, No Iterations 10
time step continuity errors : sum local = 0.0115466, global = 0.00306665, cumulative = 0.00309645
smoothSolver: Solving for omega, Initial residual = 0.0129123, Final residual = 0.00126264, No Iterations 3
smoothSolver: Solving for k, Initial residual = 0.235247, Final residual = 0.0201433, No Iterations 3
PIMPLE: Not converged within 10 iterations
ExecutionTime = 7.796 s ClockTime = 8 s

forces forces1 write:
sum of forces:
pressure : (-3.18227e+07 -552255 -4.87339e+08)
viscous : (33.197 -705.829 -212.958)
porous : (0 0 0)
sum of moments:
pressure : (-8.52121e+08 -1.45813e+08 5.08253e+07)
viscous : (3037.44 240.222 -3.01841)
porous : (0 0 0)

Time = 3s

PIMPLE: Iteration 1
smoothSolver: Solving for Ux, Initial residual = 0.670706, Final residual = 0.0426299, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.560932, Final residual = 0.039382, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.759684, Final residual = 0.0489113, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0267526, Final residual = 0.000203579, No Iterations 7
GAMG: Solving for p, Initial residual = 0.000206444, Final residual = 1.97147e-06, No Iterations 12
time step continuity errors : sum local = 0.0258458, global = 0.00773113
GAMG: Solving for p, Initial residual = 0.00123992, Final residual = 1.10343e-05, No Iterations 4
GAMG: Solving for p, Initial residual = 0.000442311, Final residual = 4.06279e-06, No Iterations 11
time step continuity errors : sum local = 0.0152685, global = 0.00564694
PIMPLE: Iteration 2
smoothSolver: Solving for Ux, Initial residual = 0.677924, Final residual = 0.0490353, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.542447, Final residual = 0.0401997, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.762623, Final residual = 0.0526566, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0192802, Final residual = 0.000125666, No Iterations 5
GAMG: Solving for p, Initial residual = 0.000189406, Final residual = 1.78461e-06, No Iterations 12
time step continuity errors : sum local = 0.038365, global = 0.00600849
GAMG: Solving for p, Initial residual = 0.00115666, Final residual = 1.02544e-05, No Iterations 4
GAMG: Solving for p, Initial residual = 0.000413854, Final residual = 4.12383e-06, No Iterations 9
time step continuity errors : sum local = 0.025111, global = 0.00376544
PIMPLE: Iteration 3
...
GAMG: Solving for p, Initial residual = 1, Final residual = 0.340115, No Iterations 1000
GAMG: Solving for p, Initial residual = 1, Final residual = 3.9186e-05, No Iterations 1
time step continuity errors : sum local = 1.39615e+25, global = 2.88728e+16
GAMG: Solving for p, Initial residual = 0.972745, Final residual = 0.30648, No Iterations 1000
GAMG: Solving for p, Initial residual = 1, Final residual = 1.25823e-05, No Iterations 1
time step continuity errors : sum local = 3.94821e+35, global = 1.77325e+28
PIMPLE: Iteration 8
smoothSolver: Solving for Ux, Initial residual = 0.934004, Final residual = 3.42843e+10, No Iterations 1000
smoothSolver: Solving for Uy, Initial residual = 0.911895, Final residual = 2.42188e+11, No Iterations 1000
smoothSolver: Solving for Uz, Initial residual = 0.949511, Final residual = 1.29109e+11, No Iterations 1000




blockMeshDict:

scale 1;

vertices
(
( 30 -20 -0)
( 30 50 -0)
(-30 50 -0)
(-30 -20 -0)
( 30 -20 20)
( 30 50 20)
(-30 50 20)
(-30 -20 20)
);

blocks
(
hex (0 1 2 3 4 5 6 7) (60 70 20) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
frontAndBack
{
type patch;
faces
(
(1 5 4 0)
(3 7 6 2)
);
}
inlet
{
type patch;
faces
(
(0 4 7 3)
);
}
outlet
{
type patch;
faces
(
(2 6 5 1)
);
}
lowerWall
{
type wall;
faces
(
(0 3 2 1)
);
}
upperWall
{
type patch;
faces
(
(4 5 6 7)
);
}
);




snappyHexMeshDict:


// Which of the steps to run
castellatedMesh true;
snap true;
addLayers true;


geometry
{
canopy.stl
{
type triSurfaceMesh;
file "canopy.stl";
name walls
}

box_1
{
type searchableBox;
min (-20 -12 0);
max (20 38 14);
}

}



// Settings for the castellatedMesh generation.
castellatedMeshControls
{

// Refinement parameters
// ~~~~~~~~~~~~~~~~~~~~~

maxLocalCells 2000000;

maxGlobalCells 10000000;


minRefinementCells 10;

maxLoadUnbalance 0.10;

nCellsBetweenLevels 3;



// Explicit feature edge refinement
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


features
(
{
file "canopy.eMesh";
level 4;
}
);



// Surface based refinement
// ~~~~~~~~~~~~~~~~~~~~~~~~

refinementSurfaces
{
walls
{

level (4 4);

patchInfo
{
type wall;
}
}
}


resolveFeatureAngle 30;


// Region-wise refinement
// ~~~~~~~~~~~~~~~~~~~~~~

refinementRegions
{
walls
{
mode distance;
levels ((0.3 4));
}

box_1
{
mode inside;
levels ((1e15 1));
}


}


// Mesh selection
// ~~~~~~~~~~~~~~


locationInMesh (24.84 33.44 12.77);

allowFreeStandingZoneFaces true;
}



// Settings for the snapping.
snapControls
{

nSmoothPatch 3;

tolerance 2.0;

nSolveIter 30;

nRelaxIter 5;

// Feature snapping

nFeatureSnapIter 20;

implicitFeatureSnap false;

explicitFeatureSnap true;

multiRegionFeatureSnap true;
}



// Settings for the layer addition.
addLayersControls
{

relativeSizes false;

// Per final patch (so not geometry!) the layer information
layers
{
"walls*"
{
nSurfaceLayers 4;
}
}

// Expansion factor for layer mesh
expansionRatio 1.7;


firstLayerThickness 9.86e-04;

// 1x at 1e-2 M3
// 1x at 5e-3 M4
// 1x at 2e-3 M5
// 1x at 1e-3 M6
// 1x at 5e-4 M7
// 1x at 3.8e-4 M8

// first layer at 1e-2 x 1 with 1.8 then 4x at 4e-4

minThickness 1e-05;


nGrow 0;

// Advanced settings

featureAngle 130;

nRelaxIter 3;

nSmoothSurfaceNormals 1;


nSmoothNormals 3;


nSmoothThickness 10;


maxFaceThicknessRatio 0.5;


maxThicknessToMedialRatio 0.3;


minMedialAxisAngle 90;

nBufferCellsNoExtrude 0;

nLayerIter 50;
}



// Generic mesh quality settings. At any undoable phase these determine
// where to undo.
meshQualityControls
{
#include "meshQualityDict"


// Advanced


nSmoothScale 4;

errorReduction 0.75;
}


// Advanced

// Write flags
writeFlags
(
scalarLevels
layerSets
layerFields // write volScalarField for layer coverage
);


mergeTolerance 1e-7;

[anaduc]

I forgot to attach the picture, here it is

[AtoHM]

Checking the mesh means you have to look if the quality is sufficient. One cannot see that from the dictionaries used to generate the mesh. Use the checkMesh utility for this and check also visually.


I see you have prescribed U in the whole domain and then reference this at the inlet. Its probably more stable to start with (0 0 0) in the domain and set the inlet velocity at the inlet patch.


Good thing to do if a solver runs and crashes shortly after is to run again and set the writeInterval such, that you get an output written directly in the first/each iteration. This way you can check visually in ParaView what was calculated and might see obvious mistakes. It looks like p does not converge already in the first iteration which is a hint that the BCs are not correct (given the mesh is fine). As I said, its always good to start with a tutorial and work from there.