CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

1D COAL pyrolysis - Temperature JUMP

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 13, 2019, 04:06
Default 1D COAL pyrolysis - Temperature JUMP
  #1
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 11
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Dear Foamers,
I have created my own solver (cokeOvenFOAM) based on fireFOAM for 1D pyrolysis [considering dry coal]. As a next step, I want to create a SUBMODEL considering moisture evaporation [wet coal] and couple with my pyrolysis solver.

In attachment (DRY COAL) --> Temperature graph (cokeOvenFOAM - my prolysis solver)
In attachment (WET COAL) --> Temperature graph (Expecting result)

The difference between two graphs --> temperature jump at the bottom, which is due to moisture evaporation. I couldn't able to figure out as how to introduce or couple the moisture evaporation (phase change) in 1D pyrolysis solver.

Kindly someone share their ideas to proceed on.
Thank you
Attached Images
File Type: jpg Query.JPG (97.8 KB, 16 views)
Kummi is offline   Reply With Quote

Old   July 15, 2019, 00:12
Default
  #2
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 11
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Can someone help me out here please?
Kummi is offline   Reply With Quote

Old   July 15, 2019, 18:25
Default
  #3
Senior Member
 
Michael Alletto
Join Date: Jun 2018
Location: Bremen
Posts: 615
Rep Power: 15
mAlletto will become famous soon enough
What happens I guess is when the coal is heated up, is that partly of the heat is needed for the phase change of the water: -> you need to consider this in the energy equation of the pyrolisis of the model maybe as additional source term.

The water vapor is leaving your pyrolis region -> this has to be considered in the mass convervation equation also as source term.

What you additionlly need I guess is how to calculate the rate of evaporation: Here you have to look in the literature for a relation which gives you the evaporation rate inside of liquids inside a porous media as function of the tempartur and maybe pressure.

This is the way I would proceed to solve the problem you described.
mAlletto is offline   Reply With Quote

Old   July 17, 2019, 01:52
Default Phase transition phenomenon ==> with Only TEqn.H
  #4
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 11
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Dear Michael Alletto,
Thank you for your response. The pyrolysis model in openFOAM is constructed based on solidChemistry model using Arrhenius equation.
So firstly to make it simple, I constructed ONLY phase change as separate solver in OF211 to understand flow physics [Attachment ~ Solver explains the movement of interface from right to left for non-isothermal phase change {Stefan 1D Problem - moving boundary is met} ]

In this case, alpha (moisture content) as a function of T ~ empirical relation ~ as mentioned in TEqn.H in attachment
Quote:
alpha = 0.1*Foam::erf(4.0*(T-Tboil)/(Tv-Tl))+scalar(0.5);
But in my case, alpha is not function of T, instead the position of the interface (boiling plane) is fixed with the condition as follows:
Quote:
Initial moisture content alpha = 0.1 (10%)
@T=100deg (moiture content_alpha=0) (near the heated wall)
@T=100deg >> H2O(l) = H2O (v)
~ And I am confused here - dont know how to implement the conditional statement based on temperature for phase change phenomenon

ATTACHMENT FIGURE: CLEAR info
REF - Clear details here in this MANUSCRIPT ~ https://sci-hub.tw/https://doi.org/10.1016/0016-2361(83)90225-9
*And you are right, additional source term is added in the energy equation for the phase change. In my case, boiling and condensation source terms are not added in main energy equation, instead only convection source term is added in energy equation [details in MANUSCRIPT]
*For an initial case, the movement of water vapor is not considered. So, no mass conservation Equ. The current problem is based on only energy equation [TEqu].
Attached Images
File Type: jpg Model with conditionss.jpg (134.9 KB, 9 views)
File Type: jpg Result.jpg (107.9 KB, 7 views)
Attached Files
File Type: gz myboilingFoam_solver.tar.gz (112.7 KB, 7 views)
Kummi is offline   Reply With Quote

Old   July 17, 2019, 02:24
Default
  #5
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 11
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Based on the conditional statement of temperature, evaporation rate (mass balance) and heat balance are calculated as,
Quote:
Mass balance (r) = rho * alpha * (dx/dt) [kg/m2.s]
Heat balance (-K dT/dx) = (r) * latent heat [W/m2]
Heat balance of boiling is expressed in unit of W/m2 ==> where it couldnot be added as source term in energy equation which carries the unit W/m3 ~ And I am confused here as how to implement it as a source term in energy equation?
KINDLY SHARE YOUR IDEAS PLEASE.... THANK YOU !!
Kummi is offline   Reply With Quote

Old   July 17, 2019, 16:05
Default
  #6
Senior Member
 
Michael Alletto
Join Date: Jun 2018
Location: Bremen
Posts: 615
Rep Power: 15
mAlletto will become famous soon enough
Do you know this paper: https://reader.elsevier.com/reader/s...8B2BAEDFDB8CE0




They propose a simple relation for the source term in the energy equation to describe the heat source term due to evaporation/condensation.


This is also helpful https://indico.cern.ch/event/578092/...7/C3OrK-05.pdf


They solve two additional equations transport equations for the liquid and vapour content.
mAlletto is offline   Reply With Quote

Old   July 19, 2019, 10:12
Default Movement of boiling plane (interface) with CONDITION
  #7
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 11
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Dear Michael Alletto,
Thank you so much for your response. I have gone through these solvers [interThermalPhaseChangeFOAM, interCondensingEvaporatingFOAM] before. However, I checked once again as per your suggestion.

In interThermalPhaseChangeFOAM ==> EEqn.H
Quote:
//Use limited alpha1:
const volScalarField limAlpha1( min(max(alpha1, scalar(0)), scalar(1)) );
--> Here alpha varies between the values 0 and 1
--> alpha1 depends upon thermal phase change model which follows phase change heat, phase change volume and production
But my case looks to be very simple: REF - Clear details here in this MANUSCRIPT ~ https://sci-hub.tw/https://doi.org/10.1016/0016-2361(83)90225-9

I HAVE 2 UNCLEARED POINTS:

(i) BOILING: In my case, w (moisture content) doesn't vary with expression. Moisture content varies with higher (-) and lower (+) temperature sides. How to model Moisture content depending upon high and low-temperature side in openFOAM methodology. (ATTACHMENT 1)
Quote:
Position of boiling plane --> INTERFACE CONDITION
At T=100deg --> (moisture content -> high temp side) w(-) =0 (near the heated wall)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------
(moisture content -> low temp side) w(+) =calculated ==> and substituted in Mass balance (r) = rho * w(+) * (dx/dt)
(ii) The heat balance carries the unit of (ATTACHMENT 1) of W/m2 - So how the heat can be taken into account in source term of ENERGY Equ. because ENERGY Equ. carries the unit W/m3? (I dont know to implement it in OpenFOAM)
ENERGY EQU: rho*Cp*(delT/delt) = del/delx[K*delT/delx ] + r_gas*c_gas*delT/delx + rho*delQ/delT (W/m3)

~ Once BOILING is completed, I will take into account of condensation.
~ At T>100deg -- convection source term is added in energy equation (ATTACHMENT 2)

Could you please give me some direction about it. Thank you once again !!
Attached Images
File Type: jpg ATTACHMENT_1_Boiling.jpg (110.9 KB, 16 views)
File Type: png ATTACHMENT_2_Convection.png (153.7 KB, 12 views)
File Type: png Phase Changes.png (79.0 KB, 9 views)
Kummi is offline   Reply With Quote

Old   July 22, 2019, 11:47
Default
  #8
Member
 
Zhiheng Wang
Join Date: Mar 2016
Posts: 72
Rep Power: 10
Zhiheng Wang is on a distinguished road
Hi Kumaresh.
I got your message the problem look similar to my problem statement.
please send the email as I have mentioned in message
Zhiheng Wang is offline   Reply With Quote

Old   July 22, 2019, 21:53
Default
  #9
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 346
Rep Power: 11
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Hello Zhiheng Wang,
Thank you for your response.
I will mail you shortly.

^^
Kummi 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
UDF to store and use the coal particle temperature history ZJUXJK Fluent UDF and Scheme Programming 3 March 21, 2017 11:45
heat transfer coefficient value in floefd jason kid FloEFD, FloWorks & FloTHERM 10 October 16, 2016 17:55
Implementation temperature jump boundary condition in buoyantBoussinesqSimpleFoam hoseinhd OpenFOAM Programming & Development 1 December 14, 2015 03:28
unexpected constant Temperature on a clip surface Sungki OpenFOAM Running, Solving & CFD 0 August 4, 2015 05:50
Free stream temperature saharesobh FLUENT 1 October 11, 2012 21:52


All times are GMT -4. The time now is 05:39.