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

Coal gasification (coalChemistryFoam)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 6, 2010, 14:24
Default Coal gasification (coalChemistryFoam)
  #1
Member
 
N. A.
Join Date: May 2010
Posts: 64
Rep Power: 7
N. A. is on a distinguished road
Hi Folks,

I am trying to use coalChemistryFoam to simulate coal gasification. I am using coal particles mixed with steam coming out of inlet of a gasifier. However I get some instability when solving for chemistry indicating that the temperature is out of bounds when slving for chemistry (see error below). My questions are:

1. has anyone used coalChemistryFoam to simulate coal-gasification.
2. if yes, can you share the sample input files and set-up files.

Thanks,
Nir

-------------------------
Error
-------------------------
Solving chemistry
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = -172.4576252
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = -172.4576252
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = -168.8897125
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = -168.8897125
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.06252264858, Final residual = 2.128743217e-07, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.05006504269, Final residual = 1.698478568e-07, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.04849053656, Final residual = 1.647714921e-07, No Iterations 3
DILUPBiCG: Solving for O2, Initial residual = 9.278860932e-06, Final residual = 5.681887612e-08, No Iterations 1
DILUPBiCG: Solving for CH4, Initial residual = 1.523074818e-06, Final residual = 1.138947373e-08, No Iterations 1
DILUPBiCG: Solving for H2, Initial residual = 1.53011947e-06, Final residual = 1.144182535e-08, No Iterations 1
DILUPBiCG: Solving for CO2, Initial residual = 1.523069754e-06, Final residual = 1.138943648e-08, No Iterations 1
DILUPBiCG: Solving for H2O, Initial residual = 1.242080567e-06, Final residual = 7.387563002e-09, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 7.844066942e-07, Final residual = 7.844066942e-07, No Iterations 0
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = -165.0947803
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = -165.0947803
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = -162.1038572
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = -162.1038572
DICPCG: Solving for G, Initial residual = 7.414788851e-06, Final residual = 7.414788851e-06, No Iterations 0
DILUPBiCG: Solving for p, Initial residual = 0.06213669809, Final residual = 2.115563004e-11, No Iterations 4
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.5111303e-14, global = -1.98763565e-14, cumulative = -6.750422694e-11
DILUPBiCG: Solving for p, Initial residual = 7.208921875e-09, Final residual = 7.208921875e-09, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.36121857e-12, global = -1.98763636e-14, cumulative = -6.75241033e-11
T gas min/max = 200, 450
DILUPBiCG: Solving for k, Initial residual = 0.962630372, Final residual = 3.331572521e-09, No Iterations 5
ExecutionTime = 2.91 s ClockTime = 3 s
N. A. is offline   Reply With Quote

Old   July 7, 2010, 02:42
Default
  #2
Member
 
Leosding
Join Date: Mar 2009
Posts: 51
Rep Power: 8
leosding is on a distinguished road
na,
maybe you can refer the post "Patches for OpenFOAM 1.7 on MacOS X"

maybe a bug, or the problems of configurations.

Leo
leosding is offline   Reply With Quote

Old   June 13, 2011, 23:20
Default
  #3
Member
 
桂莹
Join Date: Apr 2011
Posts: 36
Rep Power: 6
yingkun is on a distinguished road
I have the same problem,have you solved it?
yingkun is offline   Reply With Quote

Old   June 14, 2011, 03:49
Default
  #4
Member
 
Gregor Olenik
Join Date: Jun 2009
Location: http://greole.github.io/
Posts: 78
Rep Power: 8
gregor is on a distinguished road
Hi,

an easy way would be using openFOAM 1.7.1 (which solves equations for hs instead of h). At least for me it avoided some janaF errors.

gregor
gregor is offline   Reply With Quote

Old   June 14, 2011, 04:43
Default
  #5
Member
 
桂莹
Join Date: Apr 2011
Posts: 36
Rep Power: 6
yingkun is on a distinguished road
hi,
I do use OF1.7.1, and it shows the problem
yingkun is offline   Reply With Quote

Old   June 14, 2011, 04:58
Default
  #6
Member
 
Gregor Olenik
Join Date: Jun 2009
Location: http://greole.github.io/
Posts: 78
Rep Power: 8
gregor is on a distinguished road
ok sorry the inital poster did not but it is an old post anyway.
janaF errors can have numerous reasons (boundary conditons, models etc.) and without posting any details about your case there is no way to find a solution.
gregor is offline   Reply With Quote

Old   June 14, 2011, 05:28
Default
  #7
Member
 
桂莹
Join Date: Apr 2011
Posts: 36
Rep Power: 6
yingkun is on a distinguished road
Quote:
Originally Posted by gregor View Post
ok sorry the inital poster did not but it is an old post anyway.
janaF errors can have numerous reasons (boundary conditons, models etc.) and without posting any details about your case there is no way to find a solution.
hi,
the model is just as the tutorial case simplifiedSiwek, I just change the mesh and boundary conditions and relational files in constant, my case has two inlets(inlet1 jets coal and primary air at 300K,inlet2 jets preheated air at 1500K):
Screenshot.jpg

the BC like follows:
BC.tar.gz
or can you give me your QQ or sth (mine :549844146)and then I will show you more detail informations
yingkun is offline   Reply With Quote

Old   June 15, 2011, 03:35
Default
  #8
Member
 
Gregor Olenik
Join Date: Jun 2009
Location: http://greole.github.io/
Posts: 78
Rep Power: 8
gregor is on a distinguished road
Hi,

first of all i have experienced difficulties when using "mixed" boundary conditions for different species. You have a zeroGradient condition for H2,CH4,CO at the inlets and fixedValue for the other species. So use fixedValue = 0.0 for H2 as well.
What time scheme are you using ? I had trouble when using CN with a high explicit ratio . And you could try to switch of radiation,particles etc. just to eliminate possibilities.


gregor
gregor is offline   Reply With Quote

Old   June 15, 2011, 04:06
Default
  #9
Member
 
桂莹
Join Date: Apr 2011
Posts: 36
Rep Power: 6
yingkun is on a distinguished road
Quote:
Originally Posted by gregor View Post
Hi,

first of all i have experienced difficulties when using "mixed" boundary conditions for different species. You have a zeroGradient condition for H2,CH4,CO at the inlets and fixedValue for the other species. So use fixedValue = 0.0 for H2 as well.
What time scheme are you using ? I had trouble when using CN with a high explicit ratio . And you could try to switch of radiation,particles etc. just to eliminate possibilities.


gregor
yes,I do change the radiation model,time step and something, but do not work.The paticles--I don't really know where and how to change,the only correlative files are the constant\coalcloud1properties and coalcloud1positions which I don't know how to set
controlDict.txt
yingkun is offline   Reply With Quote

Old   June 15, 2011, 05:01
Default
  #10
Member
 
Gregor Olenik
Join Date: Jun 2009
Location: http://greole.github.io/
Posts: 78
Rep Power: 8
gregor is on a distinguished road
you can deactivate the particles in constant\coalcloud1properties by setting active to false. but try to change your boundary conditions first.

gregor
gregor is offline   Reply With Quote

Old   August 4, 2011, 05:28
Default yingkun
  #11
New Member
 
Naveed Iqbal
Join Date: Oct 2009
Location: Germany
Posts: 19
Rep Power: 7
niqbal is on a distinguished road
Hi,
I enter in the thread as I am also simulating the similar case, but with some variations, has you solved your problem or not ? I have changed this case to a cylindrical domain, where I have modified the coalCloud by increasing the number of particles and and some other details, it is working fine, but then I changed the boundray condition (inlet U) to some flow rate option, it is still working but it is taking too long. If the thread is still active we can share some of our experience.
Cheers,
Naveed
niqbal is offline   Reply With Quote

Old   August 8, 2011, 22:10
Default
  #12
Member
 
桂莹
Join Date: Apr 2011
Posts: 36
Rep Power: 6
yingkun is on a distinguished road
hi
I came back home those days,that why I reply so late.
I change the geometry and boundary conditions,and eject coal horizontally from the left face,now I meet some problems,it just can calculates a short time,may be we can discuss it.
yingkun is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Nonpremixed Combustion Model for Coal Gasification C.Y.Ler FLUENT 2 April 29, 2014 08:51
Coal gasification DYANA FLUENT 20 April 17, 2013 10:19
COal GAsification: defining Coal coalgas FLUENT 0 May 1, 2009 01:10
Coal Gasification Macy FLUENT 1 September 24, 2007 13:48
Coal gasification HELP C.Y.Ler FLUENT 3 January 3, 2007 18:43


All times are GMT -4. The time now is 02:03.