
[Sponsors] 
January 15, 2011, 19:59 
Odd units on UDS for Isotropic Diffusion

#1 
New Member
Adam
Join Date: Jan 2011
Posts: 4
Rep Power: 7 
In the ANSYS FLUENT user's guide http://my.fit.edu/itresources/manual...ug/node348.htm , the units for diffusivity are kg/(m*s) as opposed to the conventional units for diffusivity (m^2/s). Can anyone help me to understand this?
Thanks! 

January 17, 2011, 06:31 

#2 
New Member
Kristian Etienne Einarsrud
Join Date: Oct 2010
Location: Trondheim
Posts: 29
Rep Power: 7 
Hi,
The UDS equation is actually solving , being the fluid density and being your scalar. Multiplying your "expected" units for diffusivity with density gives the units shown in the user guide. Density is included "explicitly" in the transient and convective terms, and all other terms (i.e. diffusion and source terms) must have units kg/(m^3 s). Cheers! 

January 17, 2011, 08:35 
Thanks!

#3 
New Member
Adam
Join Date: Jan 2011
Posts: 4
Rep Power: 7 
Thanks for your help!


January 28, 2011, 10:04 

#4 
Member
Neil Duffy
Join Date: Mar 2010
Posts: 34
Rep Power: 8 
Hi,
My question is quite similar. I am solving a number of UDS's, mainly for solid mass and solid temp in a single phase simulation (the multiphase options do not allow me to select the solid in the porous media as the second phase). For the temperature UDS, the diffusivity is the thermal diffusivity and is related to the solid density. My concern is that this is not the case because all UDS diffusivities are hooked into the mixture materials. The manuals are not very clear either. Can anyone shed some light on this? Thanks, Neil 

July 11, 2011, 13:11 
question about UDS diff

#5 
Senior Member
MASOUD
Join Date: Mar 2009
Posts: 102
Rep Power: 9 
Quote:
How about if the scalar something other than 'mass fraction'? Precisely, if the scalar is electronic conductivity (1/ohm.meter), should it be multiplied by fluid density? Thanks 

July 11, 2011, 13:37 

#6 
New Member
Kristian Etienne Einarsrud
Join Date: Oct 2010
Location: Trondheim
Posts: 29
Rep Power: 7 
Hi Masoud,
I assume from your question that you are interested in modelling an electrical potential? The electrical potential is then your "scalar", while the electrical conductivity is the "diffusivity". Anyway, ss long as you solve for UDS without modifying this, it is on the form posted previously, i.e. density weighted. As long as all your values (i.e. fields, conductivity and sourceterms) are consistent, this should be ok. However, for clarity, I would take the time to rewrite the UDS equations (thorugh the DEFINE_UDS_UNSTEADY and DEFINE_UDS_FLUX) so that you can specify the electrical conductivity directly. This is especially crucial in a multiphase setting, where your density varies, but not necessarily the conductivity. Good luck! KE 

July 11, 2011, 14:26 

#7  
Senior Member
MASOUD
Join Date: Mar 2009
Posts: 102
Rep Power: 9 
Quote:
That's right. It's a single phase fuel cell model with two UDSs for the 'electronic potential' and 'ionic potential'. The scalars unit is (Volt) and the diffusivity is (1/(ohm.m)). Also, I am going to model species transport using 5 UDS for the 5 existing species, instead of Fluent species model. I tried to contact Ansys customer support but didn't get any clear answer for these questions: 1. In the material properties panel, Fluent asks for the diffusivity in terms of (Kg/m.s). Is this applicable to all scalars? 2. If the answer is yes, then have a look at the Eq. 9.1.1 of the user guide please. By ignoring the transient and convective terms (that's the case in my simulation), I'll have just diffusive and source terms. The source term unit is (Amp/m3). To be consistent, diffusivity must be (1/(ohm.m)) which is the realistic unit. Then why should I multiply the UDS diffusivity by the fluid density which makes the equation inconsistent? 3. And for the species transport modeling through UDS, again, should I multiply it by density? I appreciate your time answering these questions in advance. Masoud 

July 13, 2011, 03:08 

#8 
New Member
Kristian Etienne Einarsrud
Join Date: Oct 2010
Location: Trondheim
Posts: 29
Rep Power: 7 
Hello again,
Your attached equation (9.11) is the generic UDS equation solved by FLUENT, if no modifications are made. Hence, by dimensional analysis, the diffusivity is the massdiffusivity with dimensions kg/m.s. So, as long as you do not make use of the DEFINE_UDS_UNSTEADY and DEFINE_UDS_FLUXmacros, this will apply to all scalars. This means that you will be solving an equation for , where is the electrical potential. Consequently, for an unmodified UDS, the source term sould have units kg*(scalar unit)/m3.s. However, if you disable the unsteady and convection term, you can specify the electrical conductivity directly as the diffusivity (ignore the units given in the materials panel in this case). The source term for electrical potential should then be in its expected units, i.e. Volt/ohm.m3. To be sure that there isn't any "hidden" weighting within FLUENT, I would recommend that you run two (simple) simulations with different densities, and check that your electrical potential doesn't change. When modelling species transport with an equation of the form of 7.11, assuming that Fick's law is applicable to your system, the diffusive mass flux should be , where is the massdiffusivity, with dimensions kg/m.s. Good luck! 

July 14, 2011, 00:48 

#9 
Senior Member
MASOUD
Join Date: Mar 2009
Posts: 102
Rep Power: 9 
Quote:
Just a quick question; this is NOT about UDS, it's species model! I need to impose a convective mass transfer B.C. to an external boundary. As you may know, Fluent doesn't allow to specify species mass flux on WALL/MASS FLOW INLET/VELOCITY INLET. We have to either use zeroflux or specified mass fraction for each species, which is not what I need. How would you deal with this? Here is the B.C.: DY=h*(Y0Y) I know the other option is to use UDS instead of species model and then specify UDS flux but it's a big headache! 

July 14, 2011, 04:44 

#10 
New Member
Kristian Etienne Einarsrud
Join Date: Oct 2010
Location: Trondheim
Posts: 29
Rep Power: 7 
Glad I could help!
Regarding the implementation for a species massflux at your boundaries, I believe that this is challenging. Could you maybe use Gausstheorem to rewrite your diffusive flux as a volumetric source/sink term present only in the cells adjacent to the wall in question? As far as I can see (from dimensional considerations), the source term should then take the form where is the area of the cell through which the flux occurs. In your UDF for the source you will need to check that your cell is adjacent to the boundary in question and activate the source term only in this region. Note that this just is an idea, I am not certain that it will work as expected, but it might be worth trying out? 

July 14, 2011, 17:05 

#11 
Senior Member
MASOUD
Join Date: Mar 2009
Posts: 102
Rep Power: 9 
Quote:
Meanwhile, isn't (Cell_Area)/(Cell_Volume)=1 in 2D case? Also, is there any other setting I should care about or we're all set just with imposing the sink/sources to the adjacent cells? Maybe setting a velocity profile? Thanks. 

Tags 
diffusion, isotropic, kg/(m*s), uds, units 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
How to define turbulent diffusion coeff. for UDS  John  FLUENT  4  October 16, 2015 12:38 
Moving mesh  Niklas Wikstrom (Wikstrom)  OpenFOAM Running, Solving & CFD  122  June 15, 2014 06:20 
Help! How to calculate the diffusion flux of a UDS  lzw  FLUENT  0  June 18, 2008 03:50 
UDS source term units  Ale  FLUENT  2  July 12, 2002 13:01 
ODD: no UDS flux across interior boundary  Bob  FLUENT  1  February 25, 2002 20:15 