Momentum source term problem

kdarc · February 10, 2009, 12:00

Hi. I have a momentum source I need to use (in 2d) that operates with a positive value in the x-direction, but a negative value in the y-direction.

Is there a way to apply source terms to individual components of vectors???

I want to do something along the lines of:
+ fvm::Sp(S,U.component(0)) - fvm::Sp(S,U.component(1))

Any ideas??

thanks,
Ken

eugene · February 10, 2009, 13:53

As far as I know its not currently possible. The diagonal coefficient always has to be a scalar, you are trying to make it a vector.

kdarc · February 10, 2009, 15:17

One way to do this may be to simply specify Ux and Uy as scalars, and solve them separately with their own different source terms. The U vector can simply be recomposed afterwards for other standard operations.

Any obvious problems (that aren't occurring to me!!) with doing this???

Thanks.

eugene · February 11, 2009, 07:52

Well the Laplacian operator that you need for the pressure equation, is not defined for a vector. I am not sure how you would get around this problem. Perhaps use a sparse tensor?

Other than this I don't see anything fundamental, but I'm sure there will be a lot of operational issues that will complicate the problem.

I have a vague recollection that someone else wanted to do this kind of thing a few years back. Have you tried searching the forum?

February 10, 2009, 12:00	Hi. I have a momentum source	#1
kdarc Member Ken Darcovich Join Date: Mar 2009 Location: Ottawa, Canada Posts: 31 Rep Power: 17	Hi. I have a momentum source I need to use (in 2d) that operates with a positive value in the x-direction, but a negative value in the y-direction. Is there a way to apply source terms to individual components of vectors??? I want to do something along the lines of: + fvm::Sp(S,U.component(0)) - fvm::Sp(S,U.component(1)) Any ideas?? thanks, Ken

February 10, 2009, 13:53	As far as I know its not curre	#2
eugene Senior Member Eugene de Villiers Join Date: Mar 2009 Posts: 725 Rep Power: 21	As far as I know its not currently possible. The diagonal coefficient always has to be a scalar, you are trying to make it a vector.

February 10, 2009, 15:17	One way to do this may be to s	#3
kdarc Member Ken Darcovich Join Date: Mar 2009 Location: Ottawa, Canada Posts: 31 Rep Power: 17	One way to do this may be to simply specify Ux and Uy as scalars, and solve them separately with their own different source terms. The U vector can simply be recomposed afterwards for other standard operations. Any obvious problems (that aren't occurring to me!!) with doing this??? Thanks.

February 11, 2009, 07:52	Well the Laplacian operator th	#4
eugene Senior Member Eugene de Villiers Join Date: Mar 2009 Posts: 725 Rep Power: 21	Well the Laplacian operator that you need for the pressure equation, is not defined for a vector. I am not sure how you would get around this problem. Perhaps use a sparse tensor? Other than this I don't see anything fundamental, but I'm sure there will be a lot of operational issues that will complicate the problem. I have a vague recollection that someone else wanted to do this kind of thing a few years back. Have you tried searching the forum?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
momentum source term	zwdi	FLUENT	14	June 27, 2017 15:40
UDF for source term in momentum equation	Enrico	FLUENT	9	May 30, 2014 11:34
How can I import a momentum source term from a dat	May Huang	FLUENT	1	September 5, 2011 05:50
Adding a momentum source term	segersson	OpenFOAM Running, Solving & CFD	5	March 2, 2006 23:06
PDE in udf for defining momentum source term	shekhar	FLUENT	0	April 12, 2004 01:43