body force

Goutam · February 8, 2012, 11:39

fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
- fvm::laplacian(nu, U)

-beta*(T - T0)*g
;

fvScalarMatrix TEqn(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(DT, T)
;

I add body force with the momentum equation, when I run my problem,

I am getting the value of p, T and U. only problem is, no variation is found

for U. I consider g as (9.71 0 0). moving top wall temp is 1, cold bottom

wall temp is 0. others walls temp is zeroGradient. its a simple 2D case.

Goutam · February 8, 2012, 11:40

fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
- fvm::laplacian(nu, U)

-beta*(T - T0)*g
);

fvScalarMatrix TEqn(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(DT, T)
);

I add body force with the momentum equation, when I run my problem,
I am getting the value of p, T and U. only problem is, no variation is found
for U. I consider g as (9.71 0 0). moving top wall temp is 1, cold bottom
wall temp is 0. others walls temp is zeroGradient. its a simple 2D case.

February 8, 2012, 11:39	body force	#1
Goutam Senior Member Goutam Saha Join Date: Dec 2011 Location: UK Posts: 131 Rep Power: 14	fvVectorMatrix UEqn ( fvm::ddt(U) + fvm::div(phi, U) - fvm::laplacian(nu, U) -beta(T - T0)g ; fvScalarMatrix TEqn( fvm::ddt(T) + fvm::div(phi, T) - fvm::laplacian(DT, T) ; I add body force with the momentum equation, when I run my problem, I am getting the value of p, T and U. only problem is, no variation is found for U. I consider g as (9.71 0 0). moving top wall temp is 1, cold bottom wall temp is 0. others walls temp is zeroGradient. its a simple 2D case.

February 8, 2012, 11:40	body force	#2
Goutam Senior Member Goutam Saha Join Date: Dec 2011 Location: UK Posts: 131 Rep Power: 14	fvVectorMatrix UEqn ( fvm::ddt(U) + fvm::div(phi, U) - fvm::laplacian(nu, U) -beta(T - T0)g ); fvScalarMatrix TEqn( fvm::ddt(T) + fvm::div(phi, T) - fvm::laplacian(DT, T) ); I add body force with the momentum equation, when I run my problem, I am getting the value of p, T and U. only problem is, no variation is found for U. I consider g as (9.71 0 0). moving top wall temp is 1, cold bottom wall temp is 0. others walls temp is zeroGradient. its a simple 2D case.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
How to add a source term (body force) to icoFoam?	chris Liu	OpenFOAM Programming & Development	4	October 27, 2014 16:38
Force can not converge	colopolo	CFX	13	October 4, 2011 23:03
Body force at the cell face	Souviktor	Fluent UDF and Scheme Programming	0	March 31, 2009 09:54
how to include body force in cfx	selvam R.P	CFX	4	November 25, 2005 05:01
UDF for lift force on a bluff body	sawa	FLUENT	2	April 11, 2005 04:06