CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

error in adsorber simulation on adsorption chiller

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 30, 2018, 02:47
Default error in adsorber simulation on adsorption chiller
  #1
New Member
 
Andre Kurniawan
Join Date: Sep 2017
Posts: 6
Rep Power: 8
andrk is on a distinguished road
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;
}
Thee error message is The fl process could not be started.


Hope you notice the error, thank you so much
andrk is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 15:59.