Simulation does exit with error

[Sidemountyucatan]

Hi all,

OS: Windows 10
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.18.4 (GitTag)
Build type: Release
Branch: releases/FreeCAD-0-18
Hash: 980bf9060e28555fecd9e3462f68ca74007b70f8
Python version: 3.6.6
Qt version: 5.6.2
Coin version: 4.0.0a
OCC version: 7.3.0
Locale: German/Switzerland (de_CH)

I tried to solve a simple case (flow channel / air stream) and there are problems when solving the case. I have attached the file and the error message below (blue text).

Does somebody whatthe Problem is ? I am no OpenFoam Insider regarding programming therefore the messages says nothing to me.

Start to write case to folder C:\Users\marcu\AppData\Local\Temp
Matching boundary conditions ...
Populating createPatchDict to update BC names
Successfully wrote case to folder C:\Users\marcu\AppData\Local\Temp
Solver run command: C:/BLUECF~1/OPENFO~1.X\..\msys64\usr\bin\bash --login -O expand_aliases -c cd "/c/Users/marcu/AppData/Local/Temp/case" && ./Allrun
C:/BLUECF~1/OPENFO~1.X\..\msys64\usr\bin\bash --login -O expand_aliases -c cd "/c/Users/marcu/AppData/Local/Temp/case" && ./Allrun
Raw command: ['C:/Program Files/FreeCAD 0.18/bin\\python.exe', '-u', 'C:\\Users\\marcu\\AppData\\Roaming\\FreeCAD\\Mod\\CfdOF\\WindowsRunWrapper.py', 'C:/BLUECF~1/OPENFO~1.X\\..\\msys64\\usr\\bin\\bash', '--login', '-O', 'expand_aliases', '-c', 'cd "/c/Users/marcu/AppData/Local/Temp/case" && ./Allrun']
Setting environment for OpenFOAM 5.x mingw-w64 Double Precision (of5-64), using MSMPI71 - please wait...
Environment is now ready. Notes:
- You can change between installed versions by running: ofmenu
- You can change to other predefined versions by running: ofmenuNew
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
/* Windows 32 and 64 bit porting by blueCAPE: http://www.bluecape.com.pt *\
| Based on Windows porting (2.0.x v4) by Symscape: http://www.symscape.com |
\*---------------------------------------------------------------------------*/
Build : 5.x-963176928289
Exec : C:/BLUECF~1/OpenFOAM-5.x/platforms/mingw_w64GccDPInt32Opt/bin/createPatch.exe -overwrite
Date : Feb 09 2020
Time : 16:20:14
Host : "LAPTOP-RPFT5VS1"
PID : 6404
I/O : uncollated
Case : C:/Users/marcu/AppData/Local/Temp/case
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 polyMesh for time = 0

Reading createPatchDict

Adding new patch wall as patch 14 from
{
type wall;
}

Adding new patch outlet as patch 15 from
{
type patch;
}

Adding new patch inlet as patch 16 from
{
type patch;
}


Moving faces from patch face0 to patch 14
Moving faces from patch face1 to patch 14
Moving faces from patch face11 to patch 14
Moving faces from patch face12 to patch 14
Moving faces from patch face13 to patch 14
Moving faces from patch face2 to patch 14
Moving faces from patch face3 to patch 14
Moving faces from patch face4 to patch 14
Moving faces from patch face5 to patch 14
Moving faces from patch face6 to patch 14
Moving faces from patch face7 to patch 14
Moving faces from patch face9 to patch 14
Moving faces from patch face8 to patch 15
Moving faces from patch face10 to patch 16

Doing topology modification to order faces.

Not synchronising points.

Removing patches with no faces in them.

Removing zero-sized patch face0 type wall at position 0
Removing zero-sized patch face1 type wall at position 1
Removing zero-sized patch face10 type wall at position 2
Removing zero-sized patch face11 type wall at position 3
Removing zero-sized patch face12 type wall at position 4
Removing zero-sized patch face13 type wall at position 5
Removing zero-sized patch face2 type wall at position 6
Removing zero-sized patch face3 type wall at position 7
Removing zero-sized patch face4 type wall at position 8
Removing zero-sized patch face5 type wall at position 9
Removing zero-sized patch face6 type wall at position 10
Removing zero-sized patch face7 type wall at position 11
Removing zero-sized patch face8 type wall at position 12
Removing zero-sized patch face9 type wall at position 13
Removing patches.
Writing repatched mesh to 0

End

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
/* Windows 32 and 64 bit porting by blueCAPE: http://www.bluecape.com.pt *\
| Based on Windows porting (2.0.x v4) by Symscape: http://www.symscape.com |
\*---------------------------------------------------------------------------*/
Build : 5.x-963176928289
Exec : C:/BLUECF~1/OpenFOAM-5.x/platforms/mingw_w64GccDPInt32Opt/bin/potentialFoam.exe -initialiseUBCs -writep
Date : Feb 09 2020
Time : 16:20:18
Host : "LAPTOP-RPFT5VS1"
PID : 23272
I/O : uncollated
Case : C:/Users/marcu/AppData/Local/Temp/case
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 mesh for time = 0


potentialFlow: Operating solver in PISO mode

Reading velocity field U

Constructing pressure field p

Constructing velocity potential field Phi

No MRF models present


Calculating potential flow
GAMG: Solving for Phi, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for Phi, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for Phi, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for Phi, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for Phi, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for Phi, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for Phi, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for Phi, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for Phi, Initial residual = 0, Final residual = 0, No Iterations 0
Continuity error = 0
Interpolated velocity error = 0

Calculating approximate pressure field
Generating stack trace...


Backtrace:
ZN10StackTraceC1Ev [0x705c1465+0x25]
module: C:\BLUECF~1\ThirdParty-5.x\platforms\mingw_w64GccDPInt32\lib\libstack_trace.dll
ZN4Foam5error10printStackERNS_7OstreamE [0x11e1c88+0x218]
module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
ZN4Foam6sigFpe13sigFpeHandlerEi [0x11e2af3+0x33]
module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
(No symbol) [0x403ccd]
module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\potentialFoam.exe
_C_specific_handler [0x7ff8aaaf8048+0x98]
module: C:\Windows\System32\msvcrt.dll
0_chkstk [0x7ff8ab7811ff+0x11f]
module: C:\Windows\SYSTEM32\ntdll.dll
RtlRaiseException [0x7ff8ab74a289+0x399]
module: C:\Windows\SYSTEM32\ntdll.dll
KiUserExceptionDispatcher [0x7ff8ab77fe6e+0x2e]
module: C:\Windows\SYSTEM32\ntdll.dll
(No symbol) [0x439323]
module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\potentialFoam.exe
(No symbol) [0x44146e]
module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\potentialFoam.exe
(No symbol) [0x461313]
module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\potentialFoam.exe
(No symbol) [0x4013f7]
module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\potentialFoam.exe
(No symbol) [0x40152b]
module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\potentialFoam.exe
BaseThreadInitThunk [0x7ff8a9fc7bd4+0x14]
module: C:\Windows\System32\KERNEL32.DLL
RtlUserThreadStart [0x7ff8ab74ced1+0x21]
module: C:\Windows\SYSTEM32\ntdll.dll

[Jee-Bee]

please edit your question and use code brackets
use the button what looks like </>

[thschrader]

When using "Potential flow" for initializing, I get the same error.
With "Specify values" it works. Pressure difference is 10 kPa from inlet to outlet.
My testcase:

pressure_difference.FCStd

(19.12 KiB) Downloaded 59 times

initialize_field.JPG (32.83 KiB) Viewed 2626 times

OS: Windows 10 (10.0)
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.19.19478 (Git)
Build type: Release
Branch: master
Hash: e6a29986eb7caf569b9353e997f310c1c2efbec8
Python version: 3.6.8
Qt version: 5.12.1
Coin version: 4.0.0a
OCC version: 7.3.0
Locale: German/Germany (de_DE)

[Sidemountyucatan]

Hi thschrader,

you are right ! With this approach it works. But what is the reason ? Simply because it is no potential flow ?
(not friction free, not rotation free) ? Or is it a code issue ?

Do you have any clue ?

Thx

[thschrader]

Hi thschrader,
you are right ! With this approach it works. But what is the reason ? Simply because it is no potential flow ?
(not friction free, not rotation free) ? Or is it a code issue ?
Do you have any clue ?
Thx

I am not a cfd-freak, but when I understand some comments on potentialFoam correct,
the solver only solves the U-field (as a starting point for another solver). Using only p as bondary
with zero velocity is maybe not correct for potentialFoam.
https://www.openfoam.com/documentation/ ... lFoam.html

See comments from Tobias Holzmann:
https://ww3.cad.de/foren/ubb/Forum527/HTML/000775.shtml

[Sidemountyucatan]

Thx a lot for your feedback.

in my solution all residuals are going towards 10^-4. Except the one for the pressure. Do you have any idea
why ?

[thschrader]

Thx a lot for your feedback.
in my solution all residuals are going towards 10^-4. Except the one for the pressure. Do you have any idea
why ?

I can only guess: coarse mesh, no mesh refinement zones, no boundary layers for turbulence model.
picture: my run with simplefoam, laminar approach. (v_max 170 km/h???)
BTW: what do you want to simulate?

simpleFoam_result.JPG (64.85 KiB) Viewed 2576 times

[Sidemountyucatan]

Constructive discussion ! I appriciate that a lot ! Basically I want see the stream in / around the two structures (red circles in picture).
(...has to do with stream in clean room service, particle separation etc.)
Regarding the velocity: I have to check but it should be 0.2 m/s =approx. 0.8 km/h. Maybe I did some mistake during input. (simulation is
currently running so I can just check later on)

Regarding residual of p: yes, you are right there is no mesh refinement. I will change that in the next run. But what do you mean with boundary consitions for turbulence models ? There is just limited choice in the menu for boundary conditions.

Which one you are reflecting for ?

[thschrader]

Freecad cfd-tutorials (Tutorial 3 is on turbulence):
https://www.opensim.co.za/training.html
cfd-tools:
https://www.cfd-online.com/Tools/

You can add boundary-layers/surface-mesh refinement after pressing the mesh-refinement button.
Boundary-layers work for cfmesh only.
But: at v=0,2 m/s and a step height of 100 mm I get Re=1330. This is very low and should be done
with a laminar calculation in the first step.

boundary_layers.JPG (84.96 KiB) Viewed 2536 times