
[Sponsors] 
July 13, 2005, 21:01 
Unit for a source term

#1 
Guest
Posts: n/a

Hi, all
I have just got a problem about the unit of a source term, hope someone can help me figure it out. To solve a bacteria decay problem, I set two volumetric additional variables, bacteria A(m^3) and chemical B(kg m^3). The source term of A (decay rate)is 0.6*A*(k*B^0.7)^1.4, which is based on experimental study. K is a constant with a unit of s^1. My question is how to set the unit for this source term? The cfxpre always remind me that I should set it as dimentionless, but the sover showed error and say it should a unit of [m^3*s^1]. Any comments are welcome. Best regards! JP 

July 13, 2005, 21:39 
Re: Unit for a source term

#2 
Guest
Posts: n/a

Hi,
The transport equation for an Aditional Variable is: d(phi)/dt + ... = ... + S(phi)So, the source term S(phi) has to have the units of phi devided by s. This is the case when you chose the option Source. However if you chose the option Total Source, CFX interprets S(phi) as being the Total Source devided by volume, and therefore the Total Source term has to have the units of phi multiplied by m^3 and devided by s. But as it seems you're using the option Source, the source term of your variable A must have units [m^3 s^1]. So I don't understand how CFXPre reminds you to set the source term as dimensionless!!!If you try to introduce a value instead of the expression, what are the units you see? 

July 14, 2005, 13:08 
Re: Unit for a source term

#3 
Guest
Posts: n/a

Hi, Rui,
You are right. I understand that in my case the source term of variable A must have units [m^3 s^1]. However, since the source term of A ( SA=0.6*A*(k*B^0.7)^1.4 )is based on experimental study, cfxpre will take the unit as (m^3*(kg*m^3)^0.7)^1.4. The error message is as follows: CEL error: Error in setting: "SA" via the expression: 0.6*A*(k*B^0.7)^1.4 [m^3*s^1]; Error detected at the end of the expression.; Exponent must be dimensionless; So I guess I should change my question as: how to make the unit of source term independent of the variables in its expression. Thanks a lot for your comments! Best regards! JP 

July 14, 2005, 13:56 
Re: Unit for a source term

#4 
Guest
Posts: n/a

Hi,
That's what I thought, CFXPre isn't telling you that the source term must be dimensionless, it is telling you that the exponent 1.4 has to be dimensionless.When you type 1.4 [m^3 s^1], this means 1.4 m^3*s^1. If A has units [m^3], B has units [kg m^3] and if you have defined K with units [s^1], you'll have to do: SA = 0.6*A*1 [s^1]*(K*1 [s]*(B*1 [kg^1 m^3])^0.7)^1.4K*1 [s] is dimensionless, and B*1 [kg^1 m^3] is also dimensionless, so SA will have units [m^3 s^1]. Note that when you type, for example, K*1 [s], it means that K * 1s, and not that K*1 has units [s] Regards, Rui 

July 14, 2005, 14:29 
Re: Unit for a source term

#5 
Guest
Posts: n/a

Hi, Rui
I just used a similar approach solved my problem, but your method is much more straightfoward. Your explaination also helps me a lot in understanding the definition of the units in CFX. Thanks a lot! Best regards! JP 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
momentum source term  zwdi  FLUENT  14  June 27, 2017 15:40 
[swak4Foam] swak4foam building problem  GGerber  OpenFOAM Community Contributions  54  April 24, 2015 16:02 
DxFoam reader update  hjasak  OpenFOAM PostProcessing  69  April 24, 2008 01:24 
DecomposePar links against liblamso0 with OpenMPI  jens_klostermann  OpenFOAM Bugs  11  June 28, 2007 17:51 
UDFs for Scalar Eqn  Fluid/Solid HT  Greg Perkins  FLUENT  0  October 11, 2000 03:43 