|
[Sponsors] |
July 13, 2005, 22: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 cfx-pre 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, 22: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 CFX-Pre 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, 14: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, cfx-pre 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, 14:56 |
Re: Unit for a source term
|
#4 |
Guest
Posts: n/a
|
Hi,
That's what I thought, CFX-Pre 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, 15: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 |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
momentum source term | zwdi | FLUENT | 14 | June 27, 2017 16:40 |
[swak4Foam] swak4foam building problem | GGerber | OpenFOAM Community Contributions | 54 | April 24, 2015 17:02 |
DxFoam reader update | hjasak | OpenFOAM Post-Processing | 69 | April 24, 2008 02:24 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 18:51 |
UDFs for Scalar Eqn - Fluid/Solid HT | Greg Perkins | FLUENT | 0 | October 11, 2000 04:43 |