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

Mass imlalance of water in air-droplets simulation using CFX

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

Like Tree1Likes
  • 1 Post By ghorrocks

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 6, 2023, 22:14
Default Mass imlalance of water in air-droplets simulation using CFX
  #1
New Member
 
Mengying
Join Date: Mar 2023
Posts: 5
Rep Power: 3
MengyingShu is on a distinguished road
Hi friends,

I'm using CFX to simulate an air-droplet flow in a radial turbine. The droplet is produced by hydrogen fuel cell stack. Phase change is not considered here. The diameter of water droplets is 0.3mm. The mass ratio of air to water is about 10 to 1. Since the volume fraction of droplet is very low, I am not sure to use dispersed fluid or particle transport fluid. Then I gave the continuous-dispersed model a go.

The simulation ran well, but in the result I found the mass flow of water is evidently imbalanced. About 15g/s water flows into the turbine, but almost none flows out. When I monitor the mass flow of water at exit, it is always zero.

Then I used a bend pipe to test the model, to exclude the influence of turbine geometry and rotating etc.

Some of my setting are as follows:
Air is continuous, and water is dispersed fluid. Free surface model is none. Particle model is used as interphase transfer model. Drag force is set as Schiller Naumann model. Boundary condition at wall is default (Volume Fraction). Total pressure and temperature is given at inlet, and static pressure at outlet. A CCL file of a bend pipe is attached.

Followings are what I've tried:
1. set Free Surface Model to standard:
The outlet mass flow slightly increased. It is about 0.5g/s when converged, but still imbalanced.
2. based on 1, change wall boundary condition to area fraction, and set the value of water as 1.
Nothing changed.
3. change the bend pipe to a straight one:
Mass is balanced.
4. using particle transport fluid, and set wall absorption as 0, both perpendicular and parallel coef. are 1.
Mass is balanced in bend pipe, but only 1g/s flows out when I change the model back to a turbine volute. The droplets just go round and round in the volute, can’t escape.
5. still particle transport fluid, and use a wall film boundary.
Transient simulation can’t converge after about 100 steps. The droplets hit the wall then like “disappear”.

Is the main problem here the wall boundary condition, or some other settings in the continuous-disperse flow? And, if I use LPT, is there any setting that can simulate the movement of droplets along the wall?

Could someone give me some advice? Thanks a lot.
Attached Files
File Type: txt bend pipe air droplets.txt (23.5 KB, 2 views)
MengyingShu is offline   Reply With Quote

Old   March 7, 2023, 08:12
Default
  #2
Senior Member
 
Gert-Jan
Join Date: Oct 2012
Location: Europe
Posts: 1,910
Rep Power: 28
Gert-Jan will become famous soon enough
When running Euler-Euler, then you could read this post:
H-energy imbalance and P-mass imbalance does not converge.
However, where I write energy, you should read massfraction.
Gert-Jan is online now   Reply With Quote

Old   March 7, 2023, 19:20
Default
  #3
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,850
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
I have done quite a lot of work recently on aerosol deposition on walls. My experience is:

When using the Eularian model when particles hit to walls you get several effects occurring:
1) The volume fraction of the aerosol builds up in the element adjacent to the wall. Once it builds up a bit (or maybe even gets to 1.0) then it overflows into adjacent cells and can start spurious flows along the wall and release particles back into the flow.
2) You want aerosol which hit the walls to be eliminated, rather than just building up in the wall adjacent cell. I tried using a source term to destroy any aerosol in the wall adjacent cell, but it lead to non-physical jumps in the volume fraction field.

In short, I have not found a way yet of modelling aerosol deposition on surfaces yet with a Eularian approach.

When you use a Lagrangian approach you have much better control over the wall behaviour, but you will require a very fine wall mesh to accurately predict which droplets hit the wall and which ones miss. Note you need to set the wall restitution coefficients to 0 or the particle will bounce off the wall, and you need a super-fine wall mesh to capture this accurately. Often the mesh required for this is far finer than the mesh required for a turbulent boundary layer.

Looking at your output file, I see two weird things:
1) You are injecting water at 85K. Do you really intend to inject water at almost -200C? Isn't it ice at those conditions?
2) You have left the restitution coefficients for the wall at 1.0. You will need this to be 0.0 if you want particles to hit the wall and stop.

And finally, You will probably need a super fine wall mesh to get mesh independence on the particles depositing on the wall.
karachun likes this.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
ghorrocks is offline   Reply With Quote

Old   March 9, 2023, 01:54
Default
  #4
New Member
 
Mengying
Join Date: Mar 2023
Posts: 5
Rep Power: 3
MengyingShu is on a distinguished road
Thank you so much. I’ll try it in Euler-Euler simulation.
MengyingShu is offline   Reply With Quote

Old   March 9, 2023, 03:13
Default
  #5
New Member
 
Mengying
Join Date: Mar 2023
Posts: 5
Rep Power: 3
MengyingShu is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
Looking at your output file, I see two weird things:
1) You are injecting water at 85K. Do you really intend to inject water at almost -200C? Isn't it ice at those conditions?
2) You have left the restitution coefficients for the wall at 1.0. You will need this to be 0.0 if you want particles to hit the wall and stop.

Thank you for the reply.

Sorry. It's a mistake. The temperature should be 85C, not Kelvin.
I've tried to set restitution coefficients as 0. Then particles stop at the inlet part of the volute, where they hit the wall. No droplets at all in the downstream domains, but we want to see how the particles 'flow' in the whole turbine.
As for the mesh, I'll improve the near wall grid and simulate the volute only to see if it changes anything.

Thanks again.
MengyingShu is offline   Reply With Quote

Old   March 9, 2023, 05:07
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,850
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
If you read my previous post you will see I recommend you do not use the Eularian approach, at least not until you fix the problem I discussed.

The Lagrangian approach models droplet in the air. Are you trying to model flow of the accumulated liquid on the surface as well? This is an entirely different model with different physics. CFX has a wall film model which is meant to do this but I have never used it - but note it was designed to model wall accumulation of fuel in the walls of IC engines, so use of this model in different applications is at the user's risk.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
ghorrocks is offline   Reply With Quote

Old   March 9, 2023, 22:19
Default
  #7
New Member
 
Mengying
Join Date: Mar 2023
Posts: 5
Rep Power: 3
MengyingShu is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
Are you trying to model flow of the accumulated liquid on the surface as well?
No, the accumulated fluid can be neglected.

In our experiment in a pipe, we see small droplets (diameter less than 0.3mm) hit the wall and don't stick on. We also see the water droplets flow out of the turbine with the air, instead of along the wall. I think it is reasonable to ignore the accumulation on the surface.

I want to model the movement of droplets in the turbine rotor that located downstream the volute. How they approach the blades: like velocity, direction, where they hit the blades, and so on. The problem is in the simulation all droplets hit the wall which they first encountered in the volute and stop (restitution coefficient as 0.0).

This result is much far away from what we see in the experiment. That's why I'm confused. I'm hoping to see at least some droplets leave the turbine with air.

In addition, is there any principle to set the 'number of positions'? I'm having 15g/s water flow with diameter 0.3mm. Are 100 particles too less? Would that be the problem?
MengyingShu is offline   Reply With Quote

Old   March 9, 2023, 22:32
Default
  #8
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,850
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
You need to do a sensitivity study on the key simulation settings. I can tell you that 100 droplets is going to be nowhere near enough. But I cannot tell you how many you need, you need to establish that with a sensitivity analysis.

A sensitivity analysis is:
1) Do a baseline simulation (like the 100 particle simulation you have just done)
2) Double or halve the important parameter (so do a 200 particle simulation)
3a) If the results of the two simulations for an output of importance are the same within a tolerance you are happy with then the baseline simulation settings is OK
3b) If the results of the two simulations are different by an unacceptable amount then return to point 1 with the 200 particle simulation as the baseline, and repeat with a 400 particle simulation.
4) Keep doubling the number of particles until you obtain convergence

I would recommend a sensitivity analysis on number of particles and mesh resolution near the wall. So that is two sensitivity studies you will need to run. I am sure you will find you need thousands to millions of particles, and a VERY fine wall mesh.

You will also find that the simulations required to obtain sensitivity convergence are large and take a long time to run. That is why CFD is run on supercomputers.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
ghorrocks is offline   Reply With Quote

Old   March 9, 2023, 23:49
Default
  #9
New Member
 
Mengying
Join Date: Mar 2023
Posts: 5
Rep Power: 3
MengyingShu is on a distinguished road
Thank you Glenn.

Glad to have some idea on this problem. I’ll sun some cases and update my results.
MengyingShu is offline   Reply With Quote

Reply

Tags
cfx, droplet, lpt, multiphase flow, turbine

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
Mass imbalance problem in multiphase water and steam CFX case Antech CFX 1 October 26, 2020 05:03
Mass Transfer Air and Water aminhgn FLUENT 0 May 22, 2019 14:13
air into box filled with water, outflow for water at bottom dieterdanger Fluent Multiphase 0 July 16, 2015 17:11
Radiation interface hinca CFX 15 January 26, 2014 18:11
CFX CHT increased mass flow rate for water matters little for outlet air temperature dingsheng1206 CFX 7 December 4, 2013 21:04


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