LES with CfdOF ?

A subforum specific to the development of the OpenFoam-based workbenches ( Cfd https://github.com/qingfengxia/Cfd and CfdOF https://github.com/jaheyns/CfdOF )

Moderator: oliveroxtoby

shIxx
Posts: 55
Joined: Wed Mar 28, 2018 10:13 am
Location: Bavaria (Germany)

LES with CfdOF ?

Postby shIxx » Fri Feb 05, 2021 8:43 am

Hallo community,

I was watching this video (https://www.youtube.com/watch?v=l8jDBPN4rXo&t=982s) and was thinking, why not write this case with CfdOF.
So I modeled a similar body and assigned the boundary conditions. Then I wrote a single phase, incompressible, turbulent case. I looked up the "controlDict" and saw that pimpleFoam is set up als solver. To not mess up the written case I chose a tutorial case that is also incompressible and run with pimpleFoam (C:\blueCFD-Core-2017\OpenFOAM-5.x\tutorials\incompressible\pimpleFoam\LES\channel395). It's the only LES case in pimpleFoam tutorial folder btw. So I copy and paste the "fvSchemes" file from the tutorial in my system folder. I have also replaced the (k | omega | ...) parts with the "(U | k | nuTilda)" ones in the fvSolution file. The "decomposeParDict" method scotch should be fine. From the constant folder I copy and paste "fvOptions", "postChannelDict" and "turbulenceProperties"

The simulation is running but it exit with error after 10 iterations. I think thats because deltaT is 0.001 and writeInterval is 0.01

What I'm doing wrong?

Error:

Code: Select all

Backtrace:
	ZN10StackTraceC1Ev [0x705c1465+0x25]
		 module: C:\BLUECF~1\ThirdParty-5.x\platforms\mingw_w64GccDPInt32\lib\libstack_trace.dll
	ZN4Foam5error10printStackERNS_7OstreamE [0xfa1c88+0x218]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
	ZN4Foam6sigFpe13sigFpeHandlerEi [0xfa2af3+0x33]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
	(No symbol) [0x4041cd]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe
	_C_specific_handler [0x7ff9bbea8048+0x98]
		 module: C:\WINDOWS\System32\msvcrt.dll
	0_chkstk [0x7ff9bdba184f+0x11f]
		 module: C:\WINDOWS\SYSTEM32\ntdll.dll
	RtlRaiseException [0x7ff9bdb6a889+0x399]
		 module: C:\WINDOWS\SYSTEM32\ntdll.dll
	KiUserExceptionDispatcher [0x7ff9bdba04be+0x2e]
		 module: C:\WINDOWS\SYSTEM32\ntdll.dll
	ZN4Foam4pow3ERNS_5FieldIdEERKNS_5UListIdEE [0xf4822d+0x2d]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
	ZN4Foam4pow3INS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_ [0x7083b2a6+0x36]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libincompressibleTurbulenceModels.dll
	ZN4Foam4pow3INS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_ [0x7083afc8+0x1b8]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libincompressibleTurbulenceModels.dll
	ZNK4Foam9LESModels4WALEINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE1kERKNS_14GeometricFieldINS_6TensorIdEENS_12fvPatchFieldENS_7volMeshEEE [0x708e263a+0x2da]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libincompressibleTurbulenceModels.dll
	ZN4Foam9LESModels4WALEINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE10correctNutEv [0x70863cc0+0x40]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libincompressibleTurbulenceModels.dll
	(No symbol) [0x44f79e]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe
	(No symbol) [0x4013f7]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe
	(No symbol) [0x40152b]
		 module: C:\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe
	BaseThreadInitThunk [0x7ff9bda17c24+0x14]
		 module: C:\WINDOWS\System32\KERNEL32.DLL
	RtlUserThreadStart [0x7ff9bdb6d4d1+0x21]
		 module: C:\WINDOWS\SYSTEM32\ntdll.dll
Case file:
https://drive.google.com/file/d/1rjHUCD ... sp=sharing


OS: Windows 10 (10.0)
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.19.23756 (Git)
Build type: Release
Branch: master
Hash: 9c6e9184930a52b165a0b7274e3a45d1006bfe67
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)

BlueCFD 2017
shIxx
Posts: 55
Joined: Wed Mar 28, 2018 10:13 am
Location: Bavaria (Germany)

Re: LES with CfdOF ?

Postby shIxx » Sun Apr 11, 2021 3:32 pm

Hallo community,

I managed to modify the case so that a large eddy simulation is possible with the help of CfdOF. It doesn't converge and I don't think the result is accurate because I have too little knowledge with openFOAM and all the parameters, but it is possible. You can run it also with the PIMPLE algorithm just don't do the last 3 edits but PISO is faster on my machine. Now it is time to study openFOAM more deeply.

https://youtu.be/r1Jt3vmwdo4

OS: Windows 10 (10.0)
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.19.23756 (Git)
Build type: Release
Branch: master
Hash: 9c6e9184930a52b165a0b7274e3a45d1006bfe67
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)
thschrader
Posts: 2146
Joined: Sat May 20, 2017 12:06 pm
Location: Germany

Re: LES with CfdOF ?

Postby thschrader » Sun Apr 11, 2021 3:56 pm

Thanks for producing the video, excellent job!
I will try your approach in the next days.

To become an expert in cfd the whole live is not enough...
https://www.youtube.com/user/kehoste/videos
https://www.youtube.com/channel/UCcqQi9 ... oUu8eYaEkg
thschrader
Posts: 2146
Joined: Sat May 20, 2017 12:06 pm
Location: Germany

Re: LES with CfdOF ?

Postby thschrader » Fri Apr 16, 2021 11:45 am

thschrader wrote: Sun Apr 11, 2021 3:56 pm I will try your approach in the next days.
Done. For me it works.
Lets see how my cell-phone antennas works with LES.
Thanks for your tip.
les_test_cfdof.JPG
les_test_cfdof.JPG (61.07 KiB) Viewed 622 times
shIxx
Posts: 55
Joined: Wed Mar 28, 2018 10:13 am
Location: Bavaria (Germany)

Re: LES with CfdOF ?

Postby shIxx » Wed Apr 21, 2021 10:16 am

Thank you for the information. Now it's getting exciting.
The next question is which parameters have to be optimized in order to achieve convergence.