[kaleekun]

I have read all the previous threads on source term addition.

My problem is this....

I have to add a source term in the mass conservation equation. This source term is a function of local 'y' co-ordinate and also local 'Uy' (y-component of U)
Lets say, for example....
ddt(rho) + div(phi) == -(rho*Uy)/y


How to add this source term on the right hand side

[gschaider]

For the correct addition of a massSource have a look at the dieselFoam-sources.

[kaleekun]

hi, i have taken a look at it...
but still i have the problem when I multiply local Uy and local y co-ordinate....

I think i face problem while accessing the local Uy and y-coordinate..

[kaleekun]

What I meant is, I still couldn't figure out how to get the y co-ordinates in an array, so that i can manipulate to the right of the equation....

Can any one elaborate on this. Awaiting a reply

[su_junwei]

Hi Kun

try
y=mesh.C().component(1);

Su Junwei

[ash]

So same can be done with Uy component also... right

I am extrememly happy

Uy=U().component(2);

or should i add like this

Uy=mesh.U().component(2)

[ngj]

Hi Aurun

Please note that the components run from 0, thus:

x:= component(0)
y:= component(1)
z:= component(2)

with := meaning defined as.

Best regards,

Niels