
[Sponsors] 
Surface Integrals of Species Mass Flow to Wall Incorrect 

LinkBack  Thread Tools  Search this Thread  Display Modes 
June 17, 2015, 13:15 
Surface Integrals of Species Mass Flow to Wall Incorrect

#1 
New Member
Brian
Join Date: Oct 2012
Posts: 6
Rep Power: 13 
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 multicomponent 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 postprocessing. 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 CFDPost, 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. 

June 25, 2015, 19:08 
Any way to calculate species wall fluxes?

#2 
New Member
Brian
Join Date: Oct 2012
Posts: 6
Rep Power: 13 
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 CFDPost when diffusion is important? Do I have to write a UDF? 

January 9, 2016, 00:00 

#3 
Member
sean read
Join Date: Sep 2014
Location: Melbourne, Australia
Posts: 41
Rep Power: 11 
Did you have any luck with this? I am also trying to compute species flux at the wall in fluent


January 11, 2016, 23:58 

#4 
New Member
Brian
Join Date: Oct 2012
Posts: 6
Rep Power: 13 
Sorry, no. I contacted ANSYS and they confirmed that it is a bug that they will fix at some point.


January 12, 2016, 16:38 

#5 
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,672
Rep Power: 65 
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.


Tags 
mass conservation 
Thread Tools  Search this Thread 
Display Modes  


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 