|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Anais Ducrot
Join Date: Jul 2024
Posts: 3
Rep Power: 2 ![]() |
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: ![]() [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: ![]() [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" } // ************************************************** *********************** // |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
M
Join Date: Dec 2017
Posts: 692
Rep Power: 12 ![]() |
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. |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Anais Ducrot
Join Date: Jul 2024
Posts: 3
Rep Power: 2 ![]() |
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; |
|
![]() |
![]() |
![]() |
![]() |
#4 |
New Member
Anais Ducrot
Join Date: Jul 2024
Posts: 3
Rep Power: 2 ![]() |
I forgot to attach the picture, here it is
|
|
![]() |
![]() |
![]() |
![]() |
#5 |
Senior Member
M
Join Date: Dec 2017
Posts: 692
Rep Power: 12 ![]() |
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. |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[OpenFOAM.org] compile error in dynamicMesh and thermophysicalModels libraries | NickG | OpenFOAM Installation | 3 | December 30, 2019 00:21 |
[OpenFOAM.org] Compile OF 2.3 on Mac OS X .... the patch | gschaider | OpenFOAM Installation | 225 | August 25, 2015 19:43 |
OpenFOAM without MPI | kokizzu | OpenFOAM Installation | 4 | May 26, 2014 09:17 |
[General] 2 datas on one plot | Akuji | ParaView | 46 | December 1, 2013 14:06 |
[Gmsh] Gmsh and samplesurface | touf | OpenFOAM Meshing & Mesh Conversion | 2 | December 10, 2007 02:27 |