confusion on how to write udf for density of a mixture solution in multiphase problem

Hirakjyoti · September 8, 2017, 09:57

Hi I am new to CFD and this forum. I am simulating a problem where an aqueous solution of LiCl wih mass fraction of licl 0.3(secondary phase) is used to dehumidify moist air (primary phase) and i need to write udf for density of the solution which is a function of mass fraction and temperature. I have gone through the udf tutorials but still I am struggling with this. The equation is attached below. Please help me in writing the udf.
where rho(0,1,2,3)are constants.

Hirakjyoti · September 11, 2017, 01:21

I have tried to write the udf, but i m not sure whether it is correct or not.I need help.

#include"udf.h"
DEFINE_PROPERTY(Variable_density_rho,c,t,C_Y3)
{
Thread*c_thread;
cell_t,c;
real C_Y3(c,t,3);
return C_R(c,t)=1000*((1+pow((C_Y3/(1-C_Y3)),1)*0.540966)+(-.30372*
(pow((C_Y3/(1-C_Y3)),2))+(.100791*(pow(C_Y3/(1-C_Y3)),3))));
}

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
UDF for coal combustion in multiphase flow	niisha24	Fluent UDF and Scheme Programming	13	September 19, 2022 04:01
ATTN ALL: SOLUTON TO UDF COMPILE PROBLEM	Rizwan	Fluent UDF and Scheme Programming	40	March 18, 2018 06:05
Problem with DPM simulation with particles injection and EXECUTE_AT_THE_END UDF.	Ari	Fluent UDF and Scheme Programming	4	May 31, 2016 08:51
The problem of Stiff Chemistry Solver in Multiphase model	shenzhou1987	FLUENT	1	January 11, 2016 02:52
fluent udf problem: write specific data for every iteration in a file.	nnvoro	Fluent UDF and Scheme Programming	1	May 27, 2013 15:26

September 11, 2017, 01:21		#2
Hirakjyoti New Member Hirakjyoti Barman Join Date: Sep 2017 Posts: 7 Rep Power: 8	I have tried to write the udf, but i m not sure whether it is correct or not.I need help. #include"udf.h" DEFINE_PROPERTY(Variable_density_rho,c,t,C_Y3) { Threadc_thread; cell_t,c; real C_Y3(c,t,3); return C_R(c,t)=1000((1+pow((C_Y3/(1-C_Y3)),1)0.540966)+(-.30372 (pow((C_Y3/(1-C_Y3)),2))+(.100791*(pow(C_Y3/(1-C_Y3)),3)))); }