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

Surface Integrals of Species Mass Flow to Wall Incorrect

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

Like Tree1Likes
  • 1 Post By koad

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 17, 2015, 13:15
Default Surface Integrals of Species Mass Flow to Wall Incorrect
  #1
New Member
 
Brian
Join Date: Oct 2012
Posts: 6
Rep Power: 11
koad is on a distinguished road
I have a a Fluent simulation of CVD on the inner wall of a pipe. It seems to work, but when I try to calculate the total mass flow to the pipe wall, I get strange answers.

The simulation is fairly simple. 2D Laminar flow with multi-component diffusion in a long, straight pipe. A carrier gas and a reactive species (call it AB) come in the inlet. I'm using a reaction boundary condition at the wall. The only chemical reaction is on the wall of the pipe:

AB(g) -> A(s) + B(g)

Species A deposits on the wall, and B diffuses into the gas. Only a few percent of AB is consumed, so at the outlet there is a mix of the carrier gas, AB, and B.

The problem comes in post-processing. If I use Results->Reports->Surface Integrals->Mass Flow Rate to calculate the total mass flow at each boundary, I get something like

Inlet: 1
Outlet: -0.98
Wall: -0.02

That makes sense. But if I use Results->Reports->Surface Integrals->Flow Rate->Species->Mass Fraction of Species B to calculate the mass flow of B to the boundaries I get something like

Inlet: 0
Outlet: -0.01
Wall: -0.00001

So Fluent is saying that species B is flowing out of the system but not into it, which would violate mass conservation.

Summing the mass flows for the individual species gives the same answer as calculating the total mass flow, so that is consistent. At the wall Fluent achieves this by calculating a net flow of carrier gas to the wall, even though it is not involved in the reaction.

I've tried doing the same calculations in CFD-Post, using the massFlow and massFlowInt functions in a table, and I get the same answers.

Does anyone have any ideas on what I might be doing wrong? Or is this a defect in Fluent? I have tried running the simulation for more iterations to reduce the residuals, and decreasing the mesh size, but neither changes the answer.

On a related note, I tried changing the surface washcoat factor on the wall from 1 to 2. This doubles the mass flow of species B to each boundary, as it should, but mass is still not conserved. Also, the total deposition rate of species A, as computed with Results->Reports->Surface Integrals->Integral->Species->Surface Deposition Rate of A, is independent of the washcoat factor. So apparently Fluent doesn't take the washcoat factor into account when doing surface integrals.
koad is offline   Reply With Quote

Old   June 25, 2015, 19:08
Default Any way to calculate species wall fluxes?
  #2
New Member
 
Brian
Join Date: Oct 2012
Posts: 6
Rep Power: 11
koad is on a distinguished road
OK, so I think I figured out the problem for myself. The problem is that

Results->Reports->Surface Integrals->Flow Rate->Species->Mass Fraction of Species B

does not integrate the mass flow of species B. It calculates the integral of the total mass flow time the mass fraction of species B. Which is NOT the same thing if there is diffusion. With diffusion, the mass flux of species B is

rho_B*v_B = Y_B(rho*v) - rho*D_B*grad(Y_B)

where rho and v are the total mass density and velocity, rho_B and v_B are the mass density and velocity of species B, and Y_B and D_B are the mass fraction of species B. So if D_B = 0 you can get the mass flux of species B by multiplying the total mass flow by Y_B, but otherwise no.

So this begs the question: Is there any way to calculate the mass flux of an individual species in Fluent or CFD-Post when diffusion is important? Do I have to write a UDF?
koad is offline   Reply With Quote

Old   January 9, 2016, 00:00
Default
  #3
Member
 
sean read
Join Date: Sep 2014
Location: Melbourne, Australia
Posts: 41
Rep Power: 9
seanread is on a distinguished road
Did you have any luck with this? I am also trying to compute species flux at the wall in fluent
seanread is offline   Reply With Quote

Old   January 11, 2016, 23:58
Default
  #4
New Member
 
Brian
Join Date: Oct 2012
Posts: 6
Rep Power: 11
koad is on a distinguished road
Sorry, no. I contacted ANSYS and they confirmed that it is a bug that they will fix at some point.
ghost82 likes this.
koad is offline   Reply With Quote

Old   January 12, 2016, 16:38
Default
  #5
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 4,927
Rep Power: 58
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
I don't know what the status of the bug is but a workaround is to use a custom field function for the integrand and integrate the custom field function.
LuckyTran is offline   Reply With Quote

Reply

Tags
mass conservation

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
Target Mass Flow Rate Nitin FLUENT 9 June 17, 2017 10:30
Waterwheel shaped turbine inside a pipe simulation problem mshahed91 CFX 3 January 10, 2015 11:19
[swak4Foam] funkyDoCalc mass flow over stl surface? lourencosm OpenFOAM Community Contributions 0 November 28, 2014 07:50
Difficulty In Setting Boundary Conditions Moinul Haque CFX 4 November 25, 2014 17:30
CFX4.3 -build analysis form Chie Min CFX 5 July 12, 2001 23:19


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