[jpo]

Hello,

I tried to add a z-momentum source term to model constant gravity of 9.81 m/s^2. The UDF is below.
Gravity in operating conditions panel was turned off.

Then I removed the UDF and turned gravity on in operating conditions panel and set "operating density" to 0.


The UDF failed to model 9.81m/s^2 gravity. Would anyone know why?


UDF file>>>>>>
#include "udf.h"
DEFINE_SOURCE(gravity_source,c,t,dS,eqn)
{
/* VAR */
double xc[ND_ND];
double gravity;
double source;


/* BEGIN */
C_CENTROID(xc,c,t);
gravity = 9.81;
source = gravity*C_R(c,t); /*rho g */
return source;
}

[ComputerGuy]

JPO,

First: think about posting udf-related questions in the UDF forum (http://www.cfd-online.com/Forums/fluent-udf/)

Next: how is it failing? What are the error messages?

ComputerGuy

Originally Posted by jpo

Hello,

I tried to add a z-momentum source term to model constant gravity of 9.81 m/s^2. The UDF is below.
Gravity in operating conditions panel was turned off.

Then I removed the UDF and turned gravity on in operating conditions panel and set "operating density" to 0.


The UDF failed to model 9.81m/s^2 gravity. Would anyone know why?


UDF file>>>>>>
#include "udf.h"
DEFINE_SOURCE(gravity_source,c,t,dS,eqn)
{
/* VAR */
double xc[ND_ND];
double gravity;
double source;


/* BEGIN */
C_CENTROID(xc,c,t);
gravity = 9.81;
source = gravity*C_R(c,t); /*rho g */
return source;
}

[jpo]

ComputerGuy, Thank you for suggesting the UDF forum. I was not aware of it.

How it fails? - Residuals oscillate and do not converge... The same case without UDF (gravity is defined in the oper. conditions panel) converges well.