CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

Source term for EVAPORATION in Energy Equ. - technical difficulty ?

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

LinkBack Thread Tools Search this Thread Display Modes
Old   September 8, 2019, 03:43
Question Source term for EVAPORATION in Energy Equ. - technical difficulty ?
Senior Member
Join Date: Oct 2016
Posts: 349
Rep Power: 11
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Hello Foamers,
My problem is a kind of 1D pyrolysis (Wet Coal) modelling, which is focussed mainly on heat loss - thus only Energy (heat) equation is included in my work.
ENERGY EQU: rho*Cp*(dT/dt) = del/delx[K*dT/dx ] + rho*dQ/dT + r*cp*dT/dx (W/m3)
Unsteady conduction = Diffusion term + SOURCE TERM 1 + SOURCE TERM 2
My query is all about the SOURCE TERM 2. Concerning it, my query is all about.

SOURCE TERM 1 (Pyrolysis) = heat loss due to pyrolysis solved by Arrhenius-like degradation chemistry (SOLVED this source term based on FireFOAM 1D pyrolysis model code- NO ISSUES HERE)
SOURCE TERM 2 (Evaporation) = phase change - moisture to vapor (moisture embedded in the wet coal).
CONDITION for SOURCE TERM 2: T= 100 ==> alpha (moisture content) =0

Based on it, mass and heat balances are calculated at the boiling plane (surface) as,
Mass balance (r) = rho * alpha * Velocity [kg/m2.s]
Heat balance (-K dT/dx) = (r) * latent heat [W/m2]
Moisture Evaporation is calculated at the surface (Boiling Plane). The calculated mass balance (r) is multiplied with cp and dT/dx as SOURCE TERM 2 gives the unit of W/m3.

REF - Clear details here in this MANUSCRIPT ~ (ATKINSON, B., & MERRICK, D. (1983). Mathematical models of the thermal decomposition of coal4. Heat transfer and temperature profiles in a coke-oven charge. Fuel, 62(5), 553–561)
How the source terms can be included in such problems? I have contacted certain people, few replied as they are not familier with pyrolysis, others asked me to look into interFoam solvers. But the approach of interFoam seems to different with separate transport equation for volume fraction, which is not the same in my case. I'm trying to figure out and learn as how the problem can be attacked based on algorithm and importantly how such problems can be approached step by step in OpenFOAM//
How the SOURCE TERM 2 can be technically approached and included in OpenFOAM?
PS ~ Based on real physics, Pyrolysis (source term 1) and Evaporation of moisture (source term 2) occur simultaneously inside the chamber.

Please share your ideas, it will be highly helpful.
Thank you.
Kummi is offline   Reply With Quote

Old   September 9, 2019, 09:32
Senior Member
Join Date: Oct 2016
Posts: 349
Rep Power: 11
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Hello Foamers,
As explained above, I've been trying to setup the SOURCE TERM 2 for the pyrolysis solver. I have two ideas to discuss here: (not sure which will go really well for my case)

1) Coupling SOURCE TERM 2 inside the pyrolysis solver under "reactingOneDimCoeffs". Doing this, I guess both pyrolysis and moisture evaporation will occur simultaneously.
active true;
pyrolysisModel reactingOneDim;
regionName pyrolysisRegion;
filmCoupled true;
gasHSource true;
QrHSource false;
radFluxName Qr;
moveMesh false; // true;
minimumDelta 1e-6;
useChemistrySolvers true;
2) Introducing fvOptions for the source term. My source term is similar to "solidificationMeltingSource" under fvOption, which I need to modify it as per my need.
But in pyrolysis solver, "solidificationMeltingSource" can be a good choice? If I introduce fvOption, I have to edit my main source code pertaining to solidificaionMeltingSource to get it to work right ?
type solidificationMeltingSource;
active on;
selectionMode all;
Tmelt 273;
L 334000; // enthalpy of fusion for water [J/kg]
thermoMode thermo; // retrieve thermo properties from themo model
beta 50e-6; // thermal expansion coeff [1/K]
rhoRef 1000;
Does anyone know how this SOURCE TERM 2 will be coded in pyrolysis solver?
Any ideas about above methods?
Kindly give some suggestions. Thank you

Last edited by Kummi; September 9, 2019 at 09:37. Reason: Inclusion of SOURCE TERM
Kummi is offline   Reply With Quote


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
[Other] Tabulated thermophysicalProperties library chriss85 OpenFOAM Community Contributions 62 October 2, 2022 03:50
Source Term due to evaporation in energy transport equation styleworker OpenFOAM Programming & Development 3 September 7, 2022 03:09
polynomial BC srv537 OpenFOAM Pre-Processing 4 December 3, 2016 09:07
DxFoam reader update hjasak OpenFOAM Post-Processing 69 April 24, 2008 01:24
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 11, 2000 03:43

All times are GMT -4. The time now is 06:11.