andrk |
March 30, 2018 02:47 |
error in adsorber simulation on adsorption chiller
Hi everyone, would you like to help me to check whats wrong with my step or anything?
This is what I have done...
fluent launcher -> 2D, double presition, parallel 2 processor
time= transient
multiphase = mixture, 2 phases, slip velocity, implicit volume fraction, dispersed interface model, implicit body force
energy = on
viscous = k-epsilon (2 eq)
user defined -> function hooks adjust = 4 selected
I create new material named silica-gel, change density, cp, and thermal conductivity
Phase-1 = water-vapor
Phase-2 = silica-gel (diameter = 0.25 mm)
zone name = silicagel, porous zone (porosity 0.8),enable 1 energy source term (udf sumber_energi)
initialize from inlet
time step 1s for 100 (max 5 iterations/time step)
The geometry is like this...
https://imagizer.imageshack.com/v2/xq90/924/Zr9sef.png
And this is the UDF code...
Code:
#include "udf.h"
#define NO 71.6
#define ALPHA 3080
#define BETA 18.9
#define R 8.31429
/* PERSAMAAN ADSORPSI ISOTERMAL */
DEFINE_ADJUST(adsorpsi_isotermal,domain)
{
Thread *t;
cell_t c;
real E,A,P;
thread_loop_c(t,domain) /* loop over cell threads in domain d */
{
begin_c_loop(c,t)
{
real temp = C_T(c,t); /* temp as temperature cell field variable macro */
P = C_P(c,t) + RP_Get_Real("operating-pressure");
E = ALPHA + BETA * temp;
A = R * temp * log(1.47E9/P); /* A = RT ln (Po/P) */
C_UDSI(c,t,0) = NO * exp(-A*A/(E*E));
}
end_c_loop(c,t)
}
}
/* PERSAMAAN KAPASITAS ADSORPSI */
DEFINE_ADJUST(kapasitas_adsorpsi, domain)
{
Thread *t;
cell_t c;
real dt;
thread_loop_c(t,domain)
{
begin_c_loop(c,t)
{
dt = CURRENT_TIMESTEP;
C_UDSI(c,t,1)=C_UDSI_M1(c,t,1)+0.15*dt*(C_UDSI_M1(c,t,0)-C_UDSI_M1(c,t,1));
}
end_c_loop(c,t)
}
}
/* PERSAMAAN KONSERVASI MASSA */
DEFINE_ADJUST(konservasi_massa, domain)
{
Thread*t; cell_t c;
thread_loop_c(t,domain)
{
begin_c_loop(c,t)
{
C_UDSI(c,t,2)=-0.532224*0.15*(C_UDSI(c,t,0)-C_UDSI(c,t,1));
}
end_c_loop(c,t)
}
}
/* PERSAMAAN KONSERVASI MASSA SEBELUMNYA DIJADIKAN SOURCE TERM */
DEFINE_SOURCE(src_massa,c,t,dS,eqn)
{
dS[eqn]=0.0;
return C_UDSI_M1(c,t,2);
}
/* PERSAMAAN PANAS ADSORPSI */
DEFINE_ADJUST(panas_adsorpsi,domain)
{
Thread*t;
cell_t c;
real p_operating, p_static;
thread_loop_c(t,domain)
{
begin_c_loop(c,t)
{
C_UDSI(c,t,3)=ALPHA*sqrt(log(NO/C_UDSI_M1(c,t,0)));
}
end_c_loop(c,t)
}
}
/* PERSAMAAN ENERGI */
DEFINE_SOURCE(sumber_energi,c,t,dS,eqn)
{
real sumber_energi;
sumber_energi = -C_UDSI_M1(c,t,2)*0.496056352e3*ALPHA*sqrt(log(NO/C_UDSI_M1(c,t,0)));
dS[eqn] = 0.0;
return sumber_energi;
}
Thee error message is The fl process could not be started.
Hope you notice the error, thank you so much
|