CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Close-packing spheres natural convection diverging

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   December 17, 2021, 16:01
Post Close-packing spheres natural convection diverging
Join Date: Feb 2020
Location: Türkiye
Posts: 61
Rep Power: 6
saidc. is on a distinguished road

I'm trying to solve a laminar flow natural convection case that close-packed geometry. The cylinder's wall adiabatic. The mesh generated with this steps:
  1. Generate sphere coordinates with desired porosity by DEM code (cemfDEM).
  2. Create spheres and neck bridges for contact points approach with SALOME and fuse all together.
  3. Generate 2D mesh with NETGEN 1D-2D algorithm with SALOME and export as a STL file (Images are attached).
  4. Create background cylinder with blockMesh and create 3D grid with snappyHexMesh with desired features such as max-nonorthogonality < 55.
Necks image:
2D mesh image:

checkMesh results attached. For quick answer:
  • 43 million cells (snapping start with 1 millon background cells)
  • Max-nonorthogonality < 55 and average is 17
  • Some (10-20) skewed faces.
buoyantPimpleFoam Solver settings are:
  • 1e-6 p_rgh and 1e-8 U tolerance
  • 1 nNonOrthogonalCorrectors
  • 3 nCorrectors
  • 1000 nOuterCorrectors
  • relaxationFactors p_rgh --> 0.3 and others -->0.7
Initial conditions are:
  • U (0 0 0)
  • p_rgh (1 atm)
  • T (fluid: 293K and spheres: 393K)
The problem is, pressure diverges and time step contuinity error is increases. After 50th inner iteration solver blows up. I tried change schemes (increased most robust to accurate ones), reducing tolerances, decreasing time step (1e-12 to 1e-21) but non of them helped. If I try to decreasing tolerances for control TSCE, pressure iterations reaches max and diverges.

Boundary conditions seem right to me but I know that if TSCE increases this is mostly related to BCs.

I'm open to suggestions. Files are attached.

Kind regards,
Attached Files
File Type: zip (25.4 KB, 3 views)

Last edited by saidc.; December 17, 2021 at 18:00.
saidc. is offline   Reply With Quote

Old   December 21, 2021, 06:12
Join Date: Feb 2020
Location: Türkiye
Posts: 61
Rep Power: 6
saidc. is on a distinguished road

Instead of using the boussinesq approximation and other equationOfStates in thermophysicalProperties, I changed the solver to buoyantBoussinesqPimpleFoam and it runs smooth for now. Really interesting.

But I don't think the results will be good because in this case the temperature difference is 100 K. Let's see what happens. I'll keep update.

Kind regards,

Last edited by saidc.; December 22, 2021 at 09:18.
saidc. is offline   Reply With Quote


buoyantpimplefoam, natural convection, packed bed, sphere

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Very simple natural convection problem Naseem FLUENT 19 December 17, 2020 16:00
Natural Convection Problem - Helium marzoa STAR-CCM+ 0 April 18, 2011 14:12
CFX Natural Convection, please help. dbecker CFX 5 October 13, 2010 19:07
Natural convection - Inlet boundary condition max91 CFX 1 July 29, 2008 20:28
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29

All times are GMT -4. The time now is 17:59.