CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Adaptive mesh refinement multiphase flow crashes

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By Venky_94

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 7, 2021, 12:47
Default Adaptive mesh refinement multiphase flow crashes
  #1
New Member
 
Join Date: Mar 2020
Posts: 11
Rep Power: 6
GV871 is on a distinguished road
Hello everyone,


I am working with multiphaseEulerFoam on a simplified case of droplet deposition on a surface. On a static mesh everything works fine.

If I try to activate the AMR, I get a couple of time steps where the refinement works fine and then the solver crashes.


Code:
Selected 0 cells for refinement out of 18000.
Selected 0 split points out of a possible 0.
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.982667 -2.29925e-22 1
water fraction, min, max = 0.0173333 0 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.982667 -2.82047e-22 1
water fraction, min, max = 0.0173333 0 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
air min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 0.000940122, Final residual = 3.51079e-08, No Iterations 1000
ExecutionTime = 30.83 s


Courant Number mean: 0.00129843 max: 0.00139923
deltaT = 5.39279e-06
Time = 3.52865e-05

Selected 44 cells for refinement out of 18000.
Refined from 18000 to 18308 cells.
Selected 0 split points out of a possible 44.
GAMG:  Solving for pcorr, Initial residual = 0.999998, Final residual = 6.6057e-09, No Iterations 22
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.98267 -1.22775e-08 1
water fraction, min, max = 0.0173299 -4.67376e-21 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.98267 -8.65858e-09 1
water fraction, min, max = 0.0173299 -1.05202e-24 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
air min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 0.00598428, Final residual = 3.25813e-08, No Iterations 1000
ExecutionTime = 32.31 s


Courant Number mean: 0.00176672 max: 0.00190832
deltaT = 6.47135e-06
Time = 4.17578e-05

Selected 524 cells for refinement out of 18308.
Refined from 18308 to 21976 cells.
Selected 0 split points out of a possible 524.
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 6.67404e-09, No Iterations 19
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.981933 -8.84395e-18 1.00001
water fraction, min, max = 0.0180674 -7.57791e-06 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.981932 -9.83074e-08 1.00001
water fraction, min, max = 0.0180676 -7.636e-06 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
air min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 0.010646, Final residual = 1.35044e-08, No Iterations 1000
ExecutionTime = 34.2 s


Courant Number mean: 0.00251206 max: 0.00946315
deltaT = 7.28027e-06
Time = 4.90381e-05

Selected 599 cells for refinement out of 21976.
Refined from 21976 to 26169 cells.
Selected 234 split points out of a possible 1061.
Unrefined from 26169 to 24531 cells.
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 5.30496e-09, No Iterations 20
[0] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::phaseSystem::solve(Foam::PtrList<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::PtrList<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) at ??:?
[0] #4  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
[0] #5  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #6  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
[GV-Linux:2495036] *** Process received signal ***
[GV-Linux:2495036] Signal: Floating point exception (8)
[GV-Linux:2495036] Signal code:  (-6)
[GV-Linux:2495036] Failing at address: 0x3e80026123c
[GV-Linux:2495036] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x46210)[0x7f9f5307e210]
[GV-Linux:2495036] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7f9f5307e18b]
[GV-Linux:2495036] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x46210)[0x7f9f5307e210]
[GV-Linux:2495036] [ 3] /opt/openfoam9/platforms/linux64GccDPInt32Opt/lib/libphaseSystem.so(_ZN4Foam11phaseSystem5solveERKNS_7PtrListINS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEEEERKNS1_INS2_IdNS_13fvsPatchFieldENS_11surfaceMeshEEEEE+0x134d)[0x7f9f5654f50d]
[GV-Linux:2495036] [ 4] multiphaseEulerFoam(+0x3eaa1)[0x557c4e8f9aa1]
[GV-Linux:2495036] [ 5] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f9f5305f0b3]
[GV-Linux:2495036] [ 6] multiphaseEulerFoam(+0x4b72e)[0x557c4e90672e]
[GV-Linux:2495036] *** End of error message ***
Whatever suggestion you might have, I will greatly appreciate. The case files are attached, if you want to have fun .
Thanks
Attached Files
File Type: zip drop_amr.zip (67.1 KB, 13 views)
GV871 is offline   Reply With Quote

Old   March 29, 2022, 06:56
Default
  #2
Senior Member
 
Nguyen Duy Trong
Join Date: Apr 2014
Posts: 124
Rep Power: 12
ndtrong is on a distinguished road
Quote:
Originally Posted by GV871 View Post
Hello everyone,


I am working with multiphaseEulerFoam on a simplified case of droplet deposition on a surface. On a static mesh everything works fine.

If I try to activate the AMR, I get a couple of time steps where the refinement works fine and then the solver crashes.


Code:
Selected 0 cells for refinement out of 18000.
Selected 0 split points out of a possible 0.
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.982667 -2.29925e-22 1
water fraction, min, max = 0.0173333 0 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.982667 -2.82047e-22 1
water fraction, min, max = 0.0173333 0 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
air min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 0.000940122, Final residual = 3.51079e-08, No Iterations 1000
ExecutionTime = 30.83 s


Courant Number mean: 0.00129843 max: 0.00139923
deltaT = 5.39279e-06
Time = 3.52865e-05

Selected 44 cells for refinement out of 18000.
Refined from 18000 to 18308 cells.
Selected 0 split points out of a possible 44.
GAMG:  Solving for pcorr, Initial residual = 0.999998, Final residual = 6.6057e-09, No Iterations 22
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.98267 -1.22775e-08 1
water fraction, min, max = 0.0173299 -4.67376e-21 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.98267 -8.65858e-09 1
water fraction, min, max = 0.0173299 -1.05202e-24 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
air min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 0.00598428, Final residual = 3.25813e-08, No Iterations 1000
ExecutionTime = 32.31 s


Courant Number mean: 0.00176672 max: 0.00190832
deltaT = 6.47135e-06
Time = 4.17578e-05

Selected 524 cells for refinement out of 18308.
Refined from 18308 to 21976 cells.
Selected 0 split points out of a possible 524.
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 6.67404e-09, No Iterations 19
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.981933 -8.84395e-18 1.00001
water fraction, min, max = 0.0180674 -7.57791e-06 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.981932 -9.83074e-08 1.00001
water fraction, min, max = 0.0180676 -7.636e-06 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
air min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 0.010646, Final residual = 1.35044e-08, No Iterations 1000
ExecutionTime = 34.2 s


Courant Number mean: 0.00251206 max: 0.00946315
deltaT = 7.28027e-06
Time = 4.90381e-05

Selected 599 cells for refinement out of 21976.
Refined from 21976 to 26169 cells.
Selected 234 split points out of a possible 1061.
Unrefined from 26169 to 24531 cells.
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 5.30496e-09, No Iterations 20
[0] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::phaseSystem::solve(Foam::PtrList<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::PtrList<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) at ??:?
[0] #4  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
[0] #5  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #6  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
[GV-Linux:2495036] *** Process received signal ***
[GV-Linux:2495036] Signal: Floating point exception (8)
[GV-Linux:2495036] Signal code:  (-6)
[GV-Linux:2495036] Failing at address: 0x3e80026123c
[GV-Linux:2495036] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x46210)[0x7f9f5307e210]
[GV-Linux:2495036] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7f9f5307e18b]
[GV-Linux:2495036] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x46210)[0x7f9f5307e210]
[GV-Linux:2495036] [ 3] /opt/openfoam9/platforms/linux64GccDPInt32Opt/lib/libphaseSystem.so(_ZN4Foam11phaseSystem5solveERKNS_7PtrListINS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEEEERKNS1_INS2_IdNS_13fvsPatchFieldENS_11surfaceMeshEEEEE+0x134d)[0x7f9f5654f50d]
[GV-Linux:2495036] [ 4] multiphaseEulerFoam(+0x3eaa1)[0x557c4e8f9aa1]
[GV-Linux:2495036] [ 5] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f9f5305f0b3]
[GV-Linux:2495036] [ 6] multiphaseEulerFoam(+0x4b72e)[0x557c4e90672e]
[GV-Linux:2495036] *** End of error message ***
Whatever suggestion you might have, I will greatly appreciate. The case files are attached, if you want to have fun .
Thanks
Hi,

Did you solve this issue, I am also having a similar problem with AMR technique.
Could you please share with me your solution?

Thanks
ndtrong is offline   Reply With Quote

Old   March 29, 2022, 18:04
Default
  #3
Member
 
Venkat Ganesh
Join Date: May 2020
Location: Cincinnati, Ohio
Posts: 49
Rep Power: 6
Venky_94 is on a distinguished road
Quote:
Originally Posted by GV871 View Post
Hello everyone,


I am working with multiphaseEulerFoam on a simplified case of droplet deposition on a surface. On a static mesh everything works fine.

If I try to activate the AMR, I get a couple of time steps where the refinement works fine and then the solver crashes.


Code:
Selected 0 cells for refinement out of 18000.
Selected 0 split points out of a possible 0.
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.982667 -2.29925e-22 1
water fraction, min, max = 0.0173333 0 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.982667 -2.82047e-22 1
water fraction, min, max = 0.0173333 0 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
air min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 0.000940122, Final residual = 3.51079e-08, No Iterations 1000
ExecutionTime = 30.83 s


Courant Number mean: 0.00129843 max: 0.00139923
deltaT = 5.39279e-06
Time = 3.52865e-05

Selected 44 cells for refinement out of 18000.
Refined from 18000 to 18308 cells.
Selected 0 split points out of a possible 44.
GAMG:  Solving for pcorr, Initial residual = 0.999998, Final residual = 6.6057e-09, No Iterations 22
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.98267 -1.22775e-08 1
water fraction, min, max = 0.0173299 -4.67376e-21 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.98267 -8.65858e-09 1
water fraction, min, max = 0.0173299 -1.05202e-24 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
air min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 0.00598428, Final residual = 3.25813e-08, No Iterations 1000
ExecutionTime = 32.31 s


Courant Number mean: 0.00176672 max: 0.00190832
deltaT = 6.47135e-06
Time = 4.17578e-05

Selected 524 cells for refinement out of 18308.
Refined from 18308 to 21976 cells.
Selected 0 split points out of a possible 524.
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 6.67404e-09, No Iterations 19
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.981933 -8.84395e-18 1.00001
water fraction, min, max = 0.0180674 -7.57791e-06 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.981932 -9.83074e-08 1.00001
water fraction, min, max = 0.0180676 -7.636e-06 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
air min/max T 300 - 300
water min/max T 300 - 300
GAMG:  Solving for p_rgh, Initial residual = 0.010646, Final residual = 1.35044e-08, No Iterations 1000
ExecutionTime = 34.2 s


Courant Number mean: 0.00251206 max: 0.00946315
deltaT = 7.28027e-06
Time = 4.90381e-05

Selected 599 cells for refinement out of 21976.
Refined from 21976 to 26169 cells.
Selected 234 split points out of a possible 1061.
Unrefined from 26169 to 24531 cells.
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 5.30496e-09, No Iterations 20
[0] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::phaseSystem::solve(Foam::PtrList<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::PtrList<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) at ??:?
[0] #4  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
[0] #5  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #6  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
[GV-Linux:2495036] *** Process received signal ***
[GV-Linux:2495036] Signal: Floating point exception (8)
[GV-Linux:2495036] Signal code:  (-6)
[GV-Linux:2495036] Failing at address: 0x3e80026123c
[GV-Linux:2495036] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x46210)[0x7f9f5307e210]
[GV-Linux:2495036] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7f9f5307e18b]
[GV-Linux:2495036] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x46210)[0x7f9f5307e210]
[GV-Linux:2495036] [ 3] /opt/openfoam9/platforms/linux64GccDPInt32Opt/lib/libphaseSystem.so(_ZN4Foam11phaseSystem5solveERKNS_7PtrListINS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEEEERKNS1_INS2_IdNS_13fvsPatchFieldENS_11surfaceMeshEEEEE+0x134d)[0x7f9f5654f50d]
[GV-Linux:2495036] [ 4] multiphaseEulerFoam(+0x3eaa1)[0x557c4e8f9aa1]
[GV-Linux:2495036] [ 5] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f9f5305f0b3]
[GV-Linux:2495036] [ 6] multiphaseEulerFoam(+0x4b72e)[0x557c4e90672e]
[GV-Linux:2495036] *** End of error message ***
Whatever suggestion you might have, I will greatly appreciate. The case files are attached, if you want to have fun .
Thanks
I faced a similar issue with my bubble growth simulation. Playing around with the unrefinelevel helped get past the initial hurdle. I'm however still facing the issue at a much later time step. I still haven't been able to completely resolve the issue. Please do leave a comment if you find a fix.
Sid@17 likes this.
Venky_94 is offline   Reply With Quote

Old   December 29, 2023, 13:21
Default
  #4
New Member
 
Siddhesh Pujari
Join Date: Feb 2022
Posts: 2
Rep Power: 0
Sid@17 is on a distinguished road
Quote:
Originally Posted by Venky_94 View Post
I faced a similar issue with my bubble growth simulation. Playing around with the unrefinelevel helped get past the initial hurdle. I'm however still facing the issue at a much later time step. I still haven't been able to completely resolve the issue. Please do leave a comment if you find a fix.
Hello Venky_94,

I also tried playing with an unrefined level and my simulation crashed after few time-steps. I was wondering, did you get passed this step?
Sid@17 is offline   Reply With Quote

Old   January 5, 2024, 15:55
Default
  #5
Senior Member
 
Saeed Jamshidi
Join Date: Aug 2019
Posts: 214
Rep Power: 8
saeed jamshidi is on a distinguished road
Please try this and see what happens; otherwise, share your fvsolutionDict and controlDict files.

Code:
// How often to refine
refineInterval  2;
saeed jamshidi is offline   Reply With Quote

Reply

« Previous Thread | Next Thread »

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts
BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] non uniform mesh near the stl object vava10 OpenFOAM Meshing & Mesh Conversion 0 January 31, 2021 15:41
problem with parallel calculations OF8 mpirun otaolafr OpenFOAM Programming & Development 0 January 5, 2021 08:33
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
Flow separation point and adaptive mesh refinement yhaomin2007 OpenFOAM 0 July 25, 2013 11:59
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55


All times are GMT -4. The time now is 00:15.

Contact Us - CFD Online - Privacy Statement - Top