CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT > Fluent Multiphase

Divergence in hopper flow using DDPM and DEM Collision

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 19, 2016, 08:57
Default Divergence in hopper flow using DDPM and DEM Collision
  #1
New Member
 
Evan Johnson
Join Date: Jun 2016
Posts: 6
Rep Power: 9
evanJ is on a distinguished road
Hello CFD community,

I am modeling a hopper type multiphase/granular flow of sand particles and air, driven by gravity, with the DDPM and “DEM Collision” models in Fluent 16.2. As I have read in other research, a transient approach is necessary for a hopper simulation, where first particles are injected into the domain and allowed to come to rest (with hopper outlet closed), and then the outlet is opened and particles/air flow out.

A number of time steps after injecting particles, I always get divergence due to “pressure correction” or “pressure coupling”.

Divergence happens after a number of time steps, but it is my impression that it is when the particles start to collide and interact. For example, if I inject a burst of particles from the right boundary (as shown in the attached image showing Volume Fraction of particles), a logical solution is found at each time step until the particles hit the tapered surface, at which point it immediately diverges. (In this setup, the top, tapered-top, left, and bottom surfaces are walls, and the right surface is a “pressure-inlet” with an injection along the surface. Gravity is pulling particles to the left. )

I have tried and checked many parameters to get convergence:
- Set Under-Relaxation factors to 0.01
- Time step 5e-5
- Particle diameter is definitely less than element size
- Gravity on/off (doesn’t seem to make a difference)
- DPM Collision settings – reduced dramatically down to k = 100 and eta = .2
- Solving Volume Fraction equations only (neglecting fluid) – works, but when turning fluid equations back on it immediately diverges.
- Solvers – PC Simple or Coupled, tried both
- Granular material properties and models – tried various and used recommendations from research
- Mesh Info – Maximum Orthogonal Skew is .013, Minimum Orthogonal Quality is .99, and Maximum Aspect Ratio is 1.9. I believe these all show an adequate mesh.

Any advice on why I have this consistent divergence would be very helpful and appreciated. I have not yet found a good example or tutorial on this type of dense flow using Fluent, but if/when I solve this I will post the details to help the CFD community!
Attached Images
File Type: jpg Mesh.JPG (33.8 KB, 37 views)
File Type: jpg Volume Fraction Sand.JPG (26.8 KB, 51 views)
evanJ is offline   Reply With Quote

Old   September 22, 2016, 11:43
Default
  #2
New Member
 
shahjehan's Avatar
 
Syed Shah Jehan Gillani
Join Date: Oct 2014
Posts: 15
Rep Power: 11
shahjehan is on a distinguished road
Quote:
Originally Posted by evanJ View Post
Hello CFD community,

I am modeling a hopper type multiphase/granular flow of sand particles and air, driven by gravity, with the DDPM and “DEM Collision” models in Fluent 16.2. As I have read in other research, a transient approach is necessary for a hopper simulation, where first particles are injected into the domain and allowed to come to rest (with hopper outlet closed), and then the outlet is opened and particles/air flow out.

A number of time steps after injecting particles, I always get divergence due to “pressure correction” or “pressure coupling”.

Divergence happens after a number of time steps, but it is my impression that it is when the particles start to collide and interact. For example, if I inject a burst of particles from the right boundary (as shown in the attached image showing Volume Fraction of particles), a logical solution is found at each time step until the particles hit the tapered surface, at which point it immediately diverges. (In this setup, the top, tapered-top, left, and bottom surfaces are walls, and the right surface is a “pressure-inlet” with an injection along the surface. Gravity is pulling particles to the left. )

I have tried and checked many parameters to get convergence:
- Set Under-Relaxation factors to 0.01
- Time step 5e-5
- Particle diameter is definitely less than element size
- Gravity on/off (doesn’t seem to make a difference)
- DPM Collision settings – reduced dramatically down to k = 100 and eta = .2
- Solving Volume Fraction equations only (neglecting fluid) – works, but when turning fluid equations back on it immediately diverges.
- Solvers – PC Simple or Coupled, tried both
- Granular material properties and models – tried various and used recommendations from research
- Mesh Info – Maximum Orthogonal Skew is .013, Minimum Orthogonal Quality is .99, and Maximum Aspect Ratio is 1.9. I believe these all show an adequate mesh.

Any advice on why I have this consistent divergence would be very helpful and appreciated. I have not yet found a good example or tutorial on this type of dense flow using Fluent, but if/when I solve this I will post the details to help the CFD community!
can you share your case file. let me try it myself first.
shahjehan is offline   Reply With Quote

Old   September 23, 2016, 15:38
Default
  #3
Member
 
Bhargav Bharathan
Join Date: Jun 2015
Location: Montreal, Canada
Posts: 71
Rep Power: 10
bhargavbharathan is on a distinguished road
Hi,

Just a thought. Before you inject your particles into the hopper, is it filled with air? Did you assign air as the primary phase here? Did you also patch it in the hopper zone?

-BB
bhargavbharathan is offline   Reply With Quote

Old   September 24, 2016, 11:51
Default
  #4
New Member
 
Evan Johnson
Join Date: Jun 2016
Posts: 6
Rep Power: 9
evanJ is on a distinguished road
Thanks for your thoughts and ideas!

I am attaching a cas and data file showing injected particles which have not yet hit the tapered section.

Yes, air is the primary phase, and I initialize the domain to have a volume fraction (of sand) of 0, so it is filled with air only.
evanJ is offline   Reply With Quote

Old   September 24, 2016, 12:07
Default
  #5
New Member
 
Evan Johnson
Join Date: Jun 2016
Posts: 6
Rep Power: 9
evanJ is on a distinguished road
The files are now attached to this post.
Attached Files
File Type: zip DAT file.zip (189.1 KB, 25 views)
File Type: zip CAS file.zip (65.1 KB, 26 views)
evanJ is offline   Reply With Quote

Old   December 27, 2016, 13:20
Default Hello did you get a solution for your case?
  #6
New Member
 
Join Date: Jan 2015
Posts: 11
Rep Power: 11
Dawood Al-Mosuli is on a distinguished road
Quote:
Originally Posted by evanJ View Post
Thanks for your thoughts and ideas!

I am attaching a cas and data file showing injected particles which have not yet hit the tapered section.

Yes, air is the primary phase, and I initialize the domain to have a volume fraction (of sand) of 0, so it is filled with air only.

Hello did you make sure that the parcel size is less than minimum cell size in your mesh?
Dawood Al-Mosuli is offline   Reply With Quote

Old   December 27, 2016, 14:36
Default
  #7
New Member
 
Evan Johnson
Join Date: Jun 2016
Posts: 6
Rep Power: 9
evanJ is on a distinguished road
Hello and thank you for the reply-

Yes, I did indeed check this and tried various parcel sizes, even a single particle per parcel. I did not ever find a real solution, though I will say that the Ansys tutorial "Modeling bubbling fluidized bed using DDPM+DEM" was very helpful.

While not 100% positive, I believe some of my problems were arising from the fact that the width of my geometry is ~10 mm, and the particles are ~.3 mm in diameter. I believe the mesh size should be several times the particle size - let's say 5. This means my geometry can only have 10/(5*.3), or ~seven cells across, which is not very fine mesh at all! The tutorial works fine, but when keeping all else the same and reducing the geometry size I get divergence every time (even after trying all sorts of things).

During this study, I also built a Matlab script to output an injection or ".inj" file, which I will attach in (.txt format), that lets you inject a group of particles given some x, y, and z boundaries. I hope others will find this useful!
Attached Files
File Type: txt particle_filling_script_w_plot.txt (2.3 KB, 40 views)
evanJ is offline   Reply With Quote

Reply

Tags
ddpm, dem, dem collision, granular flow, hopper

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
DDPM (dem) model specifying collisions with boundaries grasingerm Fluent Multiphase 32 January 9, 2024 07:20
DDPM or DEM Simulation for big diameter particle roopesh99 Fluent Multiphase 6 July 15, 2016 05:06
DDPM or DEM Simulation for big diameter particle roopesh99 Mesh Generation & Pre-Processing 0 June 1, 2016 02:28
DDPM or DEM Simulation for big diameter particle roopesh99 FLUENT 0 June 1, 2016 02:18
cyclone analysis error arjun3020 FLUENT 4 May 15, 2014 02:05


All times are GMT -4. The time now is 23:02.