|
[Sponsors] |
error in adsorber simulation on adsorption chiller |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 30, 2018, 02:47 |
error in adsorber simulation on adsorption chiller
|
#1 |
New Member
Andre Kurniawan
Join Date: Sep 2017
Posts: 6
Rep Power: 9 |
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... 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; } Hope you notice the error, thank you so much |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
adsorption simulation by fluent! | owen | FLUENT | 2 | April 16, 2021 11:09 |
Surface Source - Fixed Temperature? | robtheslob | FloEFD, FloWorks & FloTHERM | 18 | May 12, 2017 02:28 |
Simulation FPEs - turbulence for transient and steady-state? | DaveR | OpenFOAM Running, Solving & CFD | 5 | March 5, 2017 15:06 |
FSI TWO-WAY SIMULATION | Smagmon | CFX | 1 | March 6, 2009 13:24 |
3-D Contaminant Dispersal Simulation | Apple L S Chan | Main CFD Forum | 1 | December 23, 1998 10:06 |