body force

Goutam · February 8, 2012, 10: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, 10: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, 10:39	body force	#1
Goutam Member Goutam Saha Join Date: Dec 2011 Location: UK Posts: 91 Rep Power: 3	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, 10:40	body force	#2
Goutam Member Goutam Saha Join Date: Dec 2011 Location: UK Posts: 91 Rep Power: 3	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.

Thread Tools
Show Printable Version Email this Page
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

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