CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

A simple problem with mfix

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 10, 2009, 20:58
Default A simple problem with mfix
  #1
New Member
 
ZRR
Join Date: Apr 2009
Posts: 9
Rep Power: 16
joshzrr is on a distinguished road
Hello:
First of all, thank you for your time.
I used the MFIX (the stable one) make a simple test, just like the Fig in the attchments (geometry of vessel and array of numerical grids). It is a 2-D problem, there are two mass_inflow BC, one PO BC, for the wall NSW is adopted. One of the MI BC MIA is in the IC, the black area is define as NSW too.
But with the same install, TESTA cannot converge, but the TESTB can. I output the profile of velocity, both TESTA and TESTB, the value is different, it clear that TESTA is greater than TESTB, but with a similar tendency .
I do not know how to solve it, hope you can help me.
Looking forward to your reply!
Best regards!

PS: The input files MFIX. dat for TESTA is as follows:
! It is made available in MFIX documents
! Run-control section
RUN_NAME = 'des_fb1'
DESCRIPTION = 'fluid bed with single central jet'
RUN_TYPE ='NEW' ! 'RESTART_1'!
UNITS = 'cgs'
TIME = 0.0 !start time
TSTOP = 10.0
DT = 5E-4 !time step
DT_MAX = 1E-3
ENERGY_EQ = .FALSE. !do not solve energy eq
SPECIES_EQ = .FALSE. .FALSE. !do not solve species eq
MOMENTUM_X_EQ(1) = .FALSE.
MOMENTUM_Y_EQ(1) = .FALSE.
MOMENTUM_Z_EQ(1) = .FALSE.
! DT_FAC = 1.0
MAX_NIT = 250
TOL_RESID = 1.E-3
! Geometry Section
COORDINATES = 'cartesian'
XLENGTH = 3. !X length
IMAX = 3 !cells in i direction
YLENGTH = 5. !height
JMAX = 5 !cells in j direction
NO_K = .TRUE. !2D, no k direction
DZ(1) = 0.4d0 !Define DZ to be particle
!for 2D DEM cases
GRAVITY = 980
! Gas-phase Section
MU_g0 = 1.8E-4 !constant gas viscosity
RO_g0 = 1.205E-3 !constant gas density
! Solids-phase Section
RO_s = 2.7 !solids density
D_p0 = 0.4 !particle diameter
e = 0.9 !restitution coefficient
Phi = 30.0 !angle of internal friction
EP_star = 0.42 !void fraction at minimum
! fluidization
! Initial Conditions Section
! 1. bed
IC_X_w(1) = 0.0 !domain
IC_X_e(1) = 3.0 ! 0 < x < 15, 0 < y < 30
IC_Y_s(1) = 0.0
IC_Y_n(1) = 5.0
! initial values in the region
IC_EP_g(1) = 1.00 ! void fraction
IC_U_g(1) = 0.0 ! x-dir gas velocity
IC_V_g(1) = 0 ! y-dir gas velocity
IC_U_s(1,1) = 0.0
IC_V_s(1,1) = 0.0
IC_T_G(1) =273.
IC_T_S(1,1)=273
! Boundary Conditions Section
! 1. Central jet
BC_X_w(1) = 1.0 ! central jet
BC_X_e(1) = 2.0 ! 7 < x < 8, y = 0
BC_Y_s(1) = 0.0
BC_Y_n(1) = 0.0
BC_TYPE(1) = 'MI' !specified mass inflow
BC_EP_g(1) = 1.0
BC_U_g(1) = 0.0
BC_V_g(1) = 0 !inlet jet vel
BC_P_g(1) = 0.0
! 1. Secondary flow 1
BC_X_w(2) = 0.0
BC_X_e(2) = 1.0
BC_Y_s(2) = 0.0
BC_Y_n(2) = 0.0
BC_TYPE(2) = 'MI' !specified mass inflow
BC_EP_g(2) = 1.0
BC_U_g(2) = 0.0
BC_V_g(2) = 0.0
BC_P_g(2) = 0.0
! 1. Secondary flow 2
BC_X_w(3) = 2.0
BC_X_e(3) = 3.0
BC_Y_s(3) = 2.0
BC_Y_n(3) = 2.0
BC_TYPE(3) = 'MI' !specified mass inflow
BC_EP_g(3) = 1.0
BC_U_g(3) = 0.0
BC_V_g(3) = 0.0
BC_P_g(3) = 0.0
! 1. Secondary flow 2
BC_X_w(5) = 2.0
BC_X_e(5) = 3.0
BC_Y_s(5) = 0.0
BC_Y_n(5) = 2.0
BC_TYPE(5) = 'NSW' !specified mass inflow
! 2. Exit
BC_X_w(4) = 0.0 ! top exit
BC_X_e(4) = 3.0 ! 0 < x < 15, y = 90
BC_Y_s(4) = 5.0
BC_Y_n(4) = 5.0
BC_TYPE(4) = 'PO' !specified pressure outflow
BC_P_g(4) = 0.0
!
! DES INPUT
PARTICLES = 1 ! Number of particles
MN = 14 ! Maximum number of neighbors allowed per particle
KN = 800000 ! Normal inter-particle collision spring constant
KT = 800000 ! Tangential " " " "
ETA_DES_N = 18 ! Normal inter-particle damping coefficient
ETA_DES_T = 18 ! Tangential " " "
MEW = 0.2 ! Inter-particle friction coeffienct
KN_W = 1200000 ! Normal particle-wall collision spring constant
KT_W = 1200000 ! Tangential " " " "
ETA_N_W = 22 ! Normal particle-wall damping coefficient
ETA_T_W = 22 ! Tangential " " "
MEW_W = 0.1 ! Particle-wall friction coefficient
DES_NEIGHBOR_SEARCH = 2 ! Particle neighbor search method; 1=n-square; 2=quadtree; 3=octree
NEIGHBOR_SEARCH_N = 1
DTSOLID_FACTOR = 0.1D0
P_TIME = 10.0
! DES Logicals
DISCRETE_ELEMENT = .TRUE.
WALLDTSPLIT = .TRUE.
DES_CONTINUUM_COUPLED = .TRUE.
TSUJI_DRAG = .TRUE.
! End DES Input
! Output Control
!
OUT_DT = 0.05 !write text file BUB01.OUT
! every 0.1s
RES_DT = 0.01 !write binary restart file
! BUB01.RES every 0.01 s
NLOG = 25 !write logfile BUB01.LOG
!every 25 time steps
FULL_LOG = .TRUE. !display residuals on screen
!SPX_DT values determine how often SPx files are written. Here BUB01.SP1, which
!contains void fraction (EP_g), is written every 0.01 s, BUB01.SP2, which contains
! gas and solids pressure (P_g, P_star), is written every 0.1 s, and so forth.
!
! EP_g P_g U_g U_s ROP_s T_g X_g
! P_star V_g V_s T_s X_s Theta Scalar
! W_g W_s
SPX_DT = 0.01 0.01 0.01 0.1 100. 100. 100. 100.0 100. 100. 100.
! The decomposition in I, J, and K directions for a Distributed Memory Parallel machine
NODESI = 1 NODESJ = 1 NODESK = 1
! Sweep Direction
LEQ_SWEEP(1) = 'ISIS'
LEQ_SWEEP(2) = 'ISIS'
LEQ_SWEEP(3) = 'ISIS'
LEQ_SWEEP(4) = 'ISIS'
LEQ_SWEEP(5) = 'ISIS'
LEQ_SWEEP(6) = 'ISIS'
LEQ_SWEEP(7) = 'ISIS'
LEQ_SWEEP(8) = 'ISIS'
LEQ_SWEEP(9) = 'ISIS'
Attached Images
File Type: jpg TEST.jpg (32.2 KB, 20 views)
joshzrr is offline   Reply With Quote

Old   December 14, 2009, 20:37
Default
  #2
New Member
 
ZRR
Join Date: Apr 2009
Posts: 9
Rep Power: 16
joshzrr is on a distinguished road
Is anyting i can do?
joshzrr is offline   Reply With Quote

Old   October 31, 2013, 05:11
Default
  #3
Member
 
mohsen
Join Date: Sep 2013
Posts: 42
Rep Power: 12
mohsen0488 is on a distinguished road
Dear joshzrr
i downloaded the MFIX files from MFIX site, and i read the readme.pdf but i don't know how to use the MFIX code in windows 7 64 bit .Fortran software is installed on my computer.
can you help me please to use MFIX code ?

Thanks a lot in advance.
mohsen0488 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Problem with simple if-comparison sega OpenFOAM Programming & Development 3 September 8, 2009 08:21
Cyclic BC problem simple geo sripplinger OpenFOAM Running, Solving & CFD 5 June 2, 2008 05:25
a simple posproc. problem, help needed Sasha FLUENT 0 February 6, 2007 17:23
Simple Meshing Problem with Gambit Mario FLUENT 4 April 18, 2003 10:52
? Simple CFD problem Davy Main CFD Forum 0 August 21, 1998 07:19


All times are GMT -4. The time now is 09:19.