[thanhphatvt]

Dear Adrian,
I appreciate your help! I found it.
Have a nice weekend!

[thanhphatvt]

Dear Bukreev,
Thanks for your reply, I’m trying to use it for my simulation now.
And I tried with the Shearsmagorinsky model but I have this error, could you help me to see it?
[SuperGeometryStatistics3D] materialNumber=0; count=280032; minPhysR=(-0.2,0.2,0.2); maxPhysR=(80,11.2,11.2)
[SuperGeometryStatistics3D] materialNumber=1; count=891332; minPhysR=(0.4,0.4,0.4); maxPhysR=(79.6,11,11)
[SuperGeometryStatistics3D] materialNumber=2; count=86880; minPhysR=(0.2,0.2,0.2); maxPhysR=(79.8,11.2,11.2)
[SuperGeometryStatistics3D] materialNumber=3; count=80; minPhysR=(0.2,4.8,4.8); maxPhysR=(0.2,6.6,6.6)
[SuperGeometryStatistics3D] materialNumber=4; count=2348; minPhysR=(79.8,0.4,0.4); maxPhysR=(79.8,11,11)
[SuperGeometryStatistics3D] countTotal[1e6]=1.26067
[prepareGeometry] Prepare Geometry … OK
[prepareLattice] Prepare Lattice …
[prepareLattice] Prepare Lattice … OK
[VortexMethod] inletArea=6.4
terminate called after throwing an instance of ‘thrust::THRUST_200400_600_NS::system::system_error’
terminate called after throwing an instance of ‘thrust::THRUST_200400_600_NS::system::system_error’
what(): CUDA free failed: cudaErrorIllegalAddress: an illegal memory access was encountered what(): CUDA free failed: cudaErrorIllegalAddress: an illegal memory access was encountered
[acmt-gpu:77057] *** Process received signal ***
[acmt-gpu:77057] Signal: Aborted (6)
[acmt-gpu:77057] Signal code: (-6)
[acmt-gpu:77057] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7f5172219520]
[acmt-gpu:77057] [ 1] /lib/x86_64-linux-gnu/libc.so.6(pthread_kill+0x12c)[0x7f517226d9fc]
[acmt-gpu:77057] [ 2] /lib/x86_64-linux-gnu/libc.so.6(raise+0x16)[0x7f5172219476]
[acmt-gpu:77057] [ 3] /lib/x86_64-linux-gnu/libc.so.6(abort+0xd3)[0x7f51721ff7f3]
[acmt-gpu:77057] [ 4] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xa2b9e)[0x7f5173c76b9e]
[acmt-gpu:77057] [ 5] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xae20c)[0x7f5173c8220c]
[acmt-gpu:77057] [ 6] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xae277)[0x7f5173c82277]
[acmt-gpu:77057] [ 7] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xae4d8)[0x7f5173c824d8]
[acmt-gpu:77057] [ 8] /home/phat/openlb/olb_nvhpc/olb-1.7r0/examples/turbulence/nozzle3d/libolbcuda.so(+0xf0d19)[0x7f5176de9d19]
[acmt-gpu:77057] [ 9] /home/phat/openlb/olb_nvhpc/olb-1.7r0/examples/turbulence/nozzle3d/libolbcuda.so(+0x18c99c)[0x7f5176e8599c]
[acmt-gpu:77057] [10] /home/phat/openlb/olb_nvhpc/olb-1.7r0/examples/turbulence/nozzle3d/libolbcuda.so(+0x113d38)[0x7f5176e0cd38]
[acmt-gpu:77057] [11] /home/phat/openlb/olb_nvhpc/olb-1.7r0/examples/turbulence/nozzle3d/libolbcuda.so(_ZN3olb23ConcreteBlockCollisionOIfNS_11descriptors5D3Q19IJEEELNS_8PlatformE2ENS_8dynamics5TupleIfS3_NS_7momenta5TupleINS7_12FixedDensityENS7_12ZeroMomentumENS7_10ZeroStressENS7_16DefineSeparatelyEEENS_10equilibria11SecondOrderENS_9collision6RevertENS5_18DefaultCombinationEEEE15applyIndividualERNS_20ConcreteBlockLatticeIfS3_LS4_2EEERNS_17ConcreteBlockMaskIfLS4_2EEE+0x937)[0x7f5176e9e917]
[acmt-gpu:77057] [12] ./nozzle3d[0x5192cc]
[acmt-gpu:77057] [13] ./nozzle3d[0x4ac401]
[acmt-gpu:77057] [14] ./nozzle3d[0x455c3b]
[acmt-gpu:77057] [15] /lib/x86_64-linux-gnu/libc.so.6(+0x29d90)[0x7f5172200d90]
[acmt-gpu:77057] [16] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0x80)[0x7f5172200e40]
[acmt-gpu:77057] [17] ./nozzle3d[0x44fc25]
[acmt-gpu:77057] *** End of error message ***

[acmt-gpu:77056] *** Process received signal ***
[acmt-gpu:77056] Signal: Aborted (6)
[acmt-gpu:77056] Signal code: (-6)
[acmt-gpu:77056] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7f88a6419520]
[acmt-gpu:77056] [ 1] /lib/x86_64-linux-gnu/libc.so.6(pthread_kill+0x12c)[0x7f88a646d9fc]
[acmt-gpu:77056] [ 2] /lib/x86_64-linux-gnu/libc.so.6(raise+0x16)[0x7f88a6419476]
[acmt-gpu:77056] [ 3] /lib/x86_64-linux-gnu/libc.so.6(abort+0xd3)[0x7f88a63ff7f3]
[acmt-gpu:77056] [ 4] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xa2b9e)[0x7f88a7e76b9e]
[acmt-gpu:77056] [ 5] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xae20c)[0x7f88a7e8220c]
[acmt-gpu:77056] [ 6] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xae277)[0x7f88a7e82277]
[acmt-gpu:77056] [ 7] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xae4d8)[0x7f88a7e824d8]
[acmt-gpu:77056] [ 8] /home/phat/openlb/olb_nvhpc/olb-1.7r0/examples/turbulence/nozzle3d/libolbcuda.so(+0xf0d19)[0x7f88ab110d19]
[acmt-gpu:77056] [ 9] /home/phat/openlb/olb_nvhpc/olb-1.7r0/examples/turbulence/nozzle3d/libolbcuda.so(+0x18c99c)[0x7f88ab1ac99c]
[acmt-gpu:77056] [10] /home/phat/openlb/olb_nvhpc/olb-1.7r0/examples/turbulence/nozzle3d/libolbcuda.so(+0x113d38)[0x7f88ab133d38]
[acmt-gpu:77056] [11] /home/phat/openlb/olb_nvhpc/olb-1.7r0/examples/turbulence/nozzle3d/libolbcuda.so(_ZN3olb23ConcreteBlockCollisionOIfNS_11descriptors5D3Q19IJEEELNS_8PlatformE2ENS_8dynamics5TupleIfS3_NS_7momenta5TupleINS7_12FixedDensityENS7_12ZeroMomentumENS7_10ZeroStressENS7_16DefineSeparatelyEEENS_10equilibria11SecondOrderENS_9collision6RevertENS5_18DefaultCombinationEEEE15applyIndividualERNS_20ConcreteBlockLatticeIfS3_LS4_2EEERNS_17ConcreteBlockMaskIfLS4_2EEE+0x937)[0x7f88ab1c5917]
[acmt-gpu:77056] [12] ./nozzle3d[0x5192cc]
[acmt-gpu:77056] [13] ./nozzle3d[0x4ac401]
[acmt-gpu:77056] [14] ./nozzle3d[0x455c3b]
[acmt-gpu:77056] [15] /lib/x86_64-linux-gnu/libc.so.6(+0x29d90)[0x7f88a6400d90]
[acmt-gpu:77056] [16] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0x80)[0x7f88a6400e40]
[acmt-gpu:77056] [17] ./nozzle3d[0x44fc25]
[acmt-gpu:77056] *** 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.
————————————————————————–
————————————————————————–
mpirun noticed that process rank 0 with PID 77056 on node acmt-gpu exited on signal 6 (Aborted).
————————————————————————–
[acmt-gpu:77052] 1 more process has sent help message help-mpi-btl-base.txt / btl:no-nics
[acmt-gpu:77052] Set MCA parameter “orte_base_help_aggregate” to 0 to see all help / error messages

Thank you!

[thanhphatvt]

I used the stlreader to load the file but I can not define the zone, like it’s a circle or rectangle, I can use the other indicator for defining it, but it’s the mouth shape so I have some difficulty doing it.
Or I can just import the mouth shape like a surface and create the cubic for the domain in the code? I think I misunderstood your first comment, right?

[thanhphatvt]

Yes, I tried the STL file based on your comment. I need to use the “Indicator” for defining the mouth zone, right? So I can not define the zone of mouth in my code

[thanhphatvt]

I have already designed the 3D model using Solidworks. I created the mouth like the pictures in Solidworks. I used the intersection for making the mouth at the inlet surface. The red zone is the mouth, and the blue zone is a wall. I have some difficulty defining the shape of the inlet at the mouth.
https://drive.google.com/file/d/1VrTQ_HDwt8gSIVI3M7B7Sz-et_lnxjjK/view?usp=sharing

[thanhphatvt]

https://drive.google.com/file/d/1uz5KvoUdt2aNpHLotPHCYLdf-1PicYvg/view?usp=drive_link
Here is my model.

[thanhphatvt]

Dear arpiuk,
Thank you so much! Now my results have the 5-7% of discrepancy. I think it’s ok now. I’m trying to simulate the cough flow with the aerosol particle with open LB.
Thanks

[thanhphatvt]

Dear Adrian and Mathias,
With the resolution is over 400, the results showed the good agreement with the references.
Thank for your help!

[thanhphatvt]

Dear Adrian,
I checked the cavity3dBenchmarks: the MLUPs is 252.7
Our simulation: the MLUPs is 216.3
Thanks

[thanhphatvt]

Dear Adrian,
I used K80, the MLUPs is 250. I disabled the GNUplot. I don’t know how to check the GPU bandwidth.
Thank you!

[thanhphatvt]

Dear Mathias and Adrian,
I already reduced the VTK export. I’m running the lid-driven cavity with Smagorinski model at Re = 18000 with resolution 300 in 1000 second. I think my simulation is too slow. It’s about 1 days for 20% of sim.
Do you have any suggestions for me to improve the computation speed?
Thank you so much!

[thanhphatvt]

Dear Adrian,
Yes, I’m using ./cavity3d for run this case.
Thanks

[thanhphatvt]

Dear Adrian,
I used gpu_openmpi.mk config. I will reduce the VTK output.
I’m running the case. I will update it after finishing.
Thanks for your help!
Have a nice day!

[thanhphatvt]

Dear Mathias,
Got it! I’ll do the mesh independence test. It runs quite slow on my desktop. I use K80 card.
Thanks

[thanhphatvt]

Dear Mathias,
I tried with resolution 200x200x200 but the result is quite the same. Should I increase the resolution?
Thank you so much!

[Author]

Posts