Next Article in Journal
Nanoplasmonic Sensor Based on Surface Plasmon-Coupled Emission: Review
Previous Article in Journal
An Accuracy-Efficiency-Power Consumption Hybrid Optimization Method for CNC Milling Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transient Numerical Simulation of Regenerative Systems with Waste Gas Recirculation Strategies in Glass Production Plant

by
Carlo Cravero
* and
Alessandro Spoladore
Department DIME, Università di Genova, 16145 Genova, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(7), 1496; https://doi.org/10.3390/app9071496
Submission received: 22 February 2019 / Revised: 26 March 2019 / Accepted: 27 March 2019 / Published: 10 April 2019
(This article belongs to the Section Mechanical Engineering)

Abstract

:

Featured Application

The numerical approach is applied to the simulation of regenerative chambers for glass production plants without or with gas recirculation strategies to reduce the NOx emissions.

Abstract

The glass production industry has one of the highest energy consumption rates and environmental emission impacts with respect to the existing industrial sectors. Glass furnaces nowadays are conceived with regenerative systems to take advantage of the residual heat from the combustion exhausts in order to increase the thermal efficiency of the system. The exhaust gases are also used in innovative systems to reduce the NOx emissions in specifically designed gas recirculation systems tailored to the glass furnace. In this paper, a numerical model for the simulation of the regenerative chambers in both steady and unsteady conditions is presented. The option of gas recirculation is also included. Special attention has been focused on the radiative heat transfer from the recirculated gases using a gas emissivity model previously developed by the authors.

1. Introduction

Regenerative heat exchanger for heat recovery in glass production processes is a type of heat exchanger where the thermal energy is intermittently stored by a solid medium usually made by a matrix of high thermal capacity and high surface to volume ratio refractories bricks. In the charging phase, a hot flow releases heat to the medium; in the discharge phase, a cold flow receives heat from the medium. In 1816, Robert Stirling invented the first thermal regenerator as a part of his well known engine [1]. Since the Industrial Revolution, thermal regenerators have been used for heat recovery from waste gases discharged by industrial processes. In 1857, Sir C.W. Siemens developed a regenerative furnace [2] where a continuous heat recovery was accomplished by a couple of regenerator chambers (RC) made of refractory material. The Siemens regenerator furnace layout is still largely used in industrial sectors where high temperature waste heat is produced, as in the case of glass industry. Figure 1 shows the layout of a typical end-port (EP) glass furnace, a regenerative plant commonly used in the container glass production industry.
The raw material is heated and molten in a 1 to 1.5 m deep rectangular shaped tank. The glass pool is mainly heated by a flame in the combustion chamber (CC). Natural gas or heavy oil are used as fuel. The combustion process typically takes over the 90% of the energy input. Additional heat is frequently introduced by means of electrodes immersed into the glass with the purpose of increasing the convective flows into the glass bath. The raw material is introduced in the so-called “dog house” placed in the first part of the glass pool. The glass melts in the central part of the tank, and then it flows into a throat to be distributed to the channels that feed the forming machines. The glass bath is heated up to temperatures that can reach 1550 °C, and the exhaust gas leaves the CC at very high temperatures, over 1450 °C, corresponding to about 80–85% of the total energy consumption of the furnace. For this reason, the heat recovery of the waste gases is essential. Usually, heat regenerators are preferred over heat recuperators due to their higher conversion efficiency that significantly compensates the higher installation costs. In fact, a well-designed heat regenerator can easily recover 60% of the heat available from the exhaust gas, while a recuperator, due to the metallic materials involved, can recover up to 40% only.
The EP furnace has the typical regenerative furnace layout for a medium size plant. These two ports connect the CC to a couple of RCs—two large brick structures that contain perforated or cruciform brick stacks made of refractory material. The stack geometry produces internal ducts shaped to optimize the fluid to solid heat transfer.
The plant design is symmetrical. The combustion is fed alternatively from the preheated air by the left or the right RC. In the first phase, the air is preheated, flowing from the bottom to the top of the left RC, and it flows in the CC by the left port. The air is mixed with the fuel injected by the burners located under the port, creating a flame that covers about two thirds of the tank length. The exhaust gases come across the right ports in the right RC. Exhaust gases flow from the top to the bottom and release heat to the solid medium. A continuous combustion process and heat recovery is guaranteed by switching the cold and the hot flows into the two RCs.
Since 1930, there has been a significant scientific effort in the development of numerical models for the design of heat regenerators [3]. Hausen [4] published a mathematical model called the “ideal regenerator model” that uses differential equations to rigorously describe the heat transfer process in thermal regenerators. To overcome the mathematical difficulty of the rigorous model by Hausen, Rummel [5], in 1931, presented a simplified model called the “pseudo-recuperator model”, where the thermal regenerator at cyclic equilibrium is modeled as a conventional recuperator where fluid and solid temperatures are the time-mean temperatures in the real regenerator phases. Only 30 years later, Lambertson [6] and Willmott [7] published a solution of the “ideal regenerator model” differential equation system. Wilmott used the model to predict the regenerator performance under variable mass flow conditions [8] and to study the effects of the gas storage on the stack during the flow inversions [9]. In recent years, Sardeshpande et al. [10] and Tapasa et al. [11] developed simulation tools to investigate the energy performance of a regenerative glass furnace. Sardeshpande et al. [12] developed a model based on the Rummel’s hypothesis that is useful to analyze the regenerator performance under different operating conditions. The numerical solution of the Hausen’s model is currently no longer a limit. In [13], Sadrameli et al. developed a numerical model of a glass furnace regenerator based on the Hausen assumptions. Nowadays, heat regenerators can be efficiently simulated with state-of-the-art computational fluid dynamic (CFD) 3D techniques. Reboussin et al. [14] developed a CFD analysis to study the heat transfer in a glass furnace regenerator during the cold and the hot phases. In [15,16], the authors described a similar CFD procedure where a specific model for the heat transfer and pressure drop phenomena into the stacking is detailed. A lower order numerical model for the regenerator—under steady operation—was set up and validated using the results from high fidelity CFD [17]. CFD models have also been developed for the Strategic Waste Gas Recirculation (WGR) system [18]. It is a thermal regenerator add-on module developed within the PrimeGlass project framework [19]. The WGR system, developed to reduce NOx emissions, is a ventilation system that feeds, at its base, the RC operating in the cold phase with a fraction of the waste gas flow rate extracted from the base of the RC operating in the hot phase. A preliminary numerical model developed for steady analysis [17] has been extended with a model for gas thermal radiation [20] in order to predict the additional thermal effects in the air phase in the presence of CO-CO2-H2O molecules (due to the exhaust gases using a WGR system). In this paper, a lower order numerical approach for the unsteady analysis of the thermal regenerator, with or without the WGR system, is presented.
The model is firstly applied to predict the thermal performance of regenerators with different geometrical dimensions to demonstrate its applicability for design purposes. The model is then used to simulate the thermal effects of the WGR system.
The present lower order model is developed for the design and development of the regenerative system to set up and optimize the plant with the WGR solution. The transient analysis option gives an increased opportunity to setup the WGR system and to optimize the switch time between hot and cold phases. In fact, the heat transfer process is altered by the presence of the exhaust gases into the combustion air due to the increase of radiation. The present model with transient analysis is a first step in the development of a digital twin for the entire furnace system.

2. Simulation Model

The numerical transient model of the thermal regenerator is implemented using the software Matlab with the Simulink tool. As in the real system, the numerical model is composed of two RCs sub-models.
Each RC model is modeled with a solid domain representing the refractory and with a fluid domain that models the cold or the hot fluid phase. In order to reproduce the temperature profile, the solid and the fluid domains are discretized with N finite cells, as shown in Figure 2. For each cell, a transient energy balance equation is added to describe the heat and mass transfer between the adjacent cells using the following hypothesis:
  • No heat transfer through the casing with the external environment,
  • No seepage or flux leakage from the casing,
  • Fluid thermal properties function of temperature only,
  • Incompressible flow.
Figure 3 shows a sketch of the energy interaction between adjacent cells. The energy balance of the generic ith fluid cell can be evaluated by Equation (1): the internal energy time variation (left side of the equation) is equal to the sum of the contributions of the heat exchanged with the adjacent solid cell and the net heat transported by the mass flow from the inlet to the outlet of the fluid cell.
The energy balance of the generic ith solid cell can be evaluated by Equation (2): the internal energy time variation (left side of the equation) is equal to the sum of the contributions of the heat exchanged with the adjacent fluid cell and upper and lower solid cells.
For   i = 1 , , N m f , i c v f , i d T ˙ f , i d t = m ˙ f ( c p f , i 1 T f , i 1 c p f , i T f , i ) + ( U s f , i A s f , i ) ( T s , i T f , i )
m s , i c v s , i d T ˙ s , i d t = ( U s s , i 1 A s s , i 1 ) ( T s , i 1 T s , i ) + ( U s s , i + 1 A s s , i + 1 ) ( T s , i + 1 T s , i ) + ( U s f , i A s f , i ) ( T f , i T s , i )
where the subscript f and s refer to the fluid and solid domain, respectively, m ˙ s is the cell mass, m ˙ f is the mass flow, c v and c p are respectively the specific heat flow at constant volume and pressure. The parameters A s f , i and A s s , i are respectively the heat transfer surface at the solid/fluid interface and at the solid/solid interface of the ith layer. U s f , i and U s s , i are respectively the thermal transmittances between the solid and the fluid cells and between the two solid cells, defined in the following Equations (3) and (4).
For   i = 1 , , N U s f , i = 1 1 h t o t , i + t h k s [ W m 2 · K ]
U s s , i = k s H N [ W m 2 · K ]
where h t o t is the total heat transfer coefficient; it includes the convective and radiant effects between the gas flow and the checkers. t h is half of the refractories thickness, k s is its thermal conduction coefficient, and H is the total checkers height.
The equation system described by Equations (1) and (2) can be re-arranged as in Equations (5) and (6).
d d t [ T f ] = [ a f ] [ T f ] + [ b f ] [ T s ] + [ c f ] T f , i n i
d d t [ T s ] = [ a s ] [ T s ] + [ b s ] [ T f ]
where [ T f ] and [ T s ] are the fluid and solid temperature arrays (Equations (7) and (8)), and T f , i n i is the fluid inlet temperature.
[ T f ] = [ T f , 1 , T f , 2 , , T f , i , , T f , N 1 , T f , N ]
[ T s ] = [ T s , 1 , T s , 2 , , T s , i , , T s , N 1 , T s , N ]
In the above equations [ a f ] and [ b f ] are squared matrices of order N with non-zero value and:
For   i = 1 , , N a f , ( i , i ) = ( m ˙ f c p f , i + U s f , i A s f , i ) m f , i c v f , i
b f , ( i , i ) = U s f , i A s f , i m f , i c v f , i
For   i = 2 , , N a f , ( i , i 1 ) = m ˙ f c p f , i m f , i c v f , i
[ c f ] is an array with order N with a non-zero value at the first element:
c f , ( 1 ) = m ˙ f c p f , i m f , i c v f , i
The terms [ a s ] and [ b s ] are squared matrices with order N with non-zero value in the following cells:
For   i = 1 , , N a s , ( i , i ) = ( U s s , i 1 A s s , i 1 + U s f , i A s f , i + U s s , i + 1 A s s , i + 1 ) m s , i c v s , i
b s , ( i , i ) = U s f , i A s f , i m s , i c v s , i
For   i = 2 , , N a s , ( i , i 1 ) = U s s , i 1 A s s , i 1 m s , i c v s , i
For   i = 1 , , N 1 a s , ( i , i + 1 ) = U s s , i + 1 A s s , i + 1 m s , i c v s , i

2.1. Fluid Properties

The model is set up to simulate a flow mixture made by a combination of oxygen, nitrogen, carbon dioxide, and water vapor. For each regeneration phase, the mass concentration w j of each j t h chemical element is given as reported in Table 1. The mixture specific heat at constant pressure c p f , is calculated as a weighted sum of the specific heats, as in Equation 17. The specific heat of each j t h element is calculated using a 4th degree polynomial function with two temperature ranges. The mixtures are considered ideal gases, i.e., the specific gas constant of the mixtures R f is a function of the chemical mass fraction and the molar mass M M j , as in Equation (18). The fluids density ρ f can then be obtained using the equation of state for an ideal gas, as in Equation (19).
The fluid specific heat at constant volume c v f can be evaluated using the Mayer relation, as in Equation (20).
c p f ( T f ) = w j   c p j ( T f ) [ J kg · K ]
R f = 8.314 w j M M j [ J kg · K ]
ρ f ( T f ) = p t o t R f T f [ kg m 3 ]
c v f ( T f ) = c p f ( T f ) R f   [ J kg · K ]
The thermal properties of the solid domain, i.e., the refractories, are kept constant with temperature and along the regenerator height. It is assumed there exists a thermal conduction coefficient k s equal to 5 W/m2·K, a specific heat c v s equal to 1200 J/kg·K, and a density of 3500 kg/m3.

2.2. Geometrical Parameters

The following geometrical data for the regenerator are given:
  • The total RC volume V t o t and the fluid-volume fraction F f .
  • The total heat transfer surface between the fluid and the solid A s f , t o t .
  • The hydraulic diameter D c h of the RC channels.
  • The refractories average thickness t h .
With these assumptions, the model geometrical parameters can be set:
A s f , i = A s f , t o t N [ m 2 ]
V s , i = ( 1 F f ) V t o t N [ m 3 ]
V f , i = F f V t o t N [ m 3 ]
Finally, the mass of each cell is calculated from Equations (24) and (25).
m f , i = ρ f ( T f ) V f , i [ kg ]
m s , i = ρ s V s , i [ kg ]

2.3. Heat Transfer Coefficients

The heat transfer characterization plays a crucial role in all heat exchanger modeling. The heat transfer in the RC is a complex phenomenon—the cold phase is characterized by a mixed forced-natural convection regime, whereas in the hot phase, the heat transfer is governed by the radiant heat emitted by the exhaust gases. In [15,16,17], the authors simulated the heat transfer in a single channel of an RC for a conventional glass plant, evaluating the total heat transfer coefficient across the RC height for cold or hot fluid phases. In [20], the attention was focused on a tool for the evaluation of the radiant gas emissivity crucial to study the WGR system strategies. The present lower order numerical model uses the values obtained in [17] for the heat transfer coefficients for the cold and the hot phases. In the case of the WGR system, the tool described in [20] is used for the evaluation of the radiant heat transfer coefficients to be implemented in the cold phase.

2.4. Simulation Model Setup

A real glass furnace regenerator is made by a couple of twin chambers, where the cold and hot flows are switched to maintain the thermal performance of the system. The alternate regenerator system is numerically modeled as a rotary regenerator with a non-continuous rotation.
The regenerator model couples two RC sub-models. The first RC, called Cold RC, simulates the cold phase fluid. The second RC, called Hot RC, simulates the hot phase fluid. The phase inversion between the RCs is modeled by switching, with a given timing, the temperature profile across the solid domain between the Cold and the Hot RC. Initial temperature profiles are set for the cold and the hot flows and for the solid domain. The mass flow rates, inlet temperature, and chemical composition are set to be constant in time. The model running in unsteady mode reaches the equilibrium after the simulation of about 100 to 200 phase inversions depending on the initial temperature conditions specified.

3. Applications to a Conventional Glass Thermal Regenerator

The regenerator model is applied to analyze the system performance with different size and operative conditions. As reported in [21], the performance of a periodic thermal storage can be evaluated using three non-dimensional parameters:
(1) The thermal efficiency η T , i.e., the ratio between the air outlet and the exhaust gas inlet temperatures (Equation (26)):
η T = T a i r , o u t T e x h , i n
(2) The storage effectiveness R , i.e., the ratio between the energy absorbed by the air in a cycle and the energy available from the exhaust gas flow (Equation (27)):
R = m ˙ a i r c p a i r τ a i r ( T ˜ a i r , o u t T a i r , i n ) m ˙ e x c p e x τ e x ( T ˜ e x , o u t T e x , i n )
(3) The capacitance utilization C U , i.e., the ratio between the energy absorbed by the air in a cycle and the maximum energy that can be stored in the solid matrix (Equation (28)):
C U = m ˙ a i r c p a i r τ a i r ( T ˜ a i r , o u t T a i r , i n ) m s c v s ( T e x h , i n T a i r , i n )
In glass furnace regenerative systems, η T can commonly approach values close to one. This is possible due to the mass flow and specific heat imbalance between the air and the waste gas phases. From Equation (27), it is clear that if η T is equal to one, R reaches its maximum value (Equation (29)). With this assumption, a specific effectiveness ε can be defined (Equation (30)) as the ratio between the actual storage effectiveness R and R m a x .
R m a x = m ˙ a i r c p a i r m ˙ e x h c p e x h
ε = R R m a x
The analyses are performed with the following reference case: a thermal regenerator where each RC has a total volume of 100 m3, F f equal to 0.7, and a total heat transfer area equal to 2900 m2. The regenerator is set up to work with a reverse time τ equal to 20 min.
The regenerator is fed with an exhaust mass flow rate equal to 4.88 kg/s at a temperature of 1135 °C and with an air mass flow rate equal to 3.77 kg/s at a temperature of 145 °C. The chemical composition of the air and exhaust gases are set as reported in Table 1. A constant linear profile of the heat exchange coefficients is set for each gas phase, as reported in [17]. Figure 4 shows the time-averaged temperature profile of the two flows and of the refractories. Figure 5 shows the same variables on the top and the bottom layers of the RC over a regeneration cycle. From the above data, the η ˜ T is 0.883 with an oscillation of about 0.02 over the cycle, and the R value is 0.482, which means ε equal to 0.834.
The unsteady regenerator model is then applied to compare three regenerative chambers of varying size, as detailed in Table 2. The heat transfer surface and the volume are gradually reduced. The inversion time equal to 20 min is kept. Figure 6 shows the value of the regenerative performance coefficients for each case. The η ˜ T increases with the RC size. The increase rate is more pronounced for changes in small RC. The R shows a trend similar to η ˜ T , and the C U value decreases by increasing the RC volume. In fact, referring to Equation (28), the increase of the thermal capacity of the storage m s c v s is more effective than the increase of the regenerator energy performance.
The model is then applied to the set of configurations from Table 2 by increasing the inversion time τ : 20 - 30 - 40 min and the results are shown in Figure 7 and Figure 8. As before, with fixed RC geometry, the time-averaged outlet temperatures of the fluid T ˜ a i r , o u t and T ˜ e x h , o u t are almost constant for the different τ values. The values of η ˜ T , R , and ε do not change with τ .

4. Model Application to Thermal Regenerator with WGR System

The unsteady regenerator model is applied to the analysis of the WGR system. The WGR system is a low NOx strategy where a portion of exhaust gas is recirculated at the base of the RC in the air phase. The experimental tests on pilot furnaces performed in the PrimeGlass project [19] have confirmed that the WGR system does have an impact on the thermal balance of the regenerative chambers due to the increase of the mass flow rate and the increase in heat exchange caused by the introduction of radiant chemical species into the combustion air.
In previous studies [16,18], the WGR system was studied using CFD models in order to optimize the distribution of the waste gas in the stacker and the port to CC of the Cold RC. The CFD domain included the complete 3D geometry of the Cold RC [17], and the WGR ports were added at the bottom of the chamber. Multiblock structured meshes were generated, and the CFD solutions were obtained using the Ansys CFX CFD platform. The CFD model solved the Reynolds Averaged Navier-Stokes (RANS) governing equations (continuity, momentum, energy) for the gas mixture using the second order upwind numerical scheme with the shear stress tensor (SST) turbulence closure added to model the effects of turbulence [15,16,17,18].
The WGR system is investigated using CFD by changing the waste gas mass flow rate and the layout of the inlet ports. The following typical post-processing of the CFD solution is performed to understand the WGR system performance. The numerical CFD results are post-processed in order to track the waste gas flow inside the RC by detecting of the O2 concentration. The resulting O2 patterns in the case of waste gas/air mass flow ratio of 20% are reported in Figure 9, Figure 10, and Figure 11. In Figure 9, the isosurface at 16% of oxygen concentration is highlighted in order to show the waste gas tracked across the Cold chamber. Figure 10 shows the oxygen concentration pattern in the stacker cross section at the stacking top. In Figure 11, the previous control surface of Figure 10 is discretized into twelve rectangular sections where the O2 and the mass flow distributions are plotted. From these analyses, it is clear that the waste gas flow fills only the first part of the chamber. The above experience gained from the application of 3D CFD has been used to inspire the WGR effects modeling into the present lower order code. The Cold RC is cut across its symmetry axis into three geometrical identical sub-sections in order to simulate the thermal effects from the WGR system. At each subsection, the mass flow, the fluid temperature, and the chemicals mass fraction are specified as boundary conditions (see Table 3). The above values can be obtained from a reference CFD solution (as in the present example) or taken from a range based on the experience gained using CFD on similar configurations.
The profiles of the heat exchange coefficient for the Cold RC are obtained by adding to the reference linearized profiles [17] (valid for a standard regenerative system) the radiant effect introduced by the gas emission model presented in [20].
Figure 12 shows the heat transfer coefficient profile for each Cold RC sub-section. The first section is not fed with the waste gas (no radiant effect added), while the second and third sections have significant waste gas fraction; this explains the important heat transfer coefficient increase due to the radiant effect.
Figure 13 shows the mean temperature profile for each Cold RC sub-section. It can be seen that the fluid temperature profile is closer to the solid thermal profile in the S-2 and S-3 sections, i.e., where there is a waste gas fraction, than it is in the S-1 section.
In order to understand the WGR thermal effects on the exchanged thermal power P t h , each regenerator subsection is calculated using Equation (31).
P t h = m ˙ f ( c p f , i n T f , i n c p f , o u t T f , o u t )
The specific thermal power (the thermal power referred to one square meter of heat transfer surface) is compared to the numerical results obtained without the WGR system (see Table 4).
It can be observed that, as far as the fluid exit temperature is concerned, the first section has a low temperature, about 983 °C, due to the low heat transfer. It is interesting to note that the second section gives the maximum outlet temperature of about 1052 °C, although it has a smaller heat transfer compared to the third one. This effect is due to the distribution of the mass flow rate and to the different inlet temperature of the flowing mixture.
The specific heat power increases from 0.97 kW/m2 in the first regeneration section to 1.87 kW/m2 in the third section. The average value in the case of the WGR with 20% of waste gas fraction is 1.29 kW/m2 compared to 1.06 kW/m2 with the system without the WGR. The total heat flux to the air with the WGR is about 3.75 MW with an increase of about 17% with respect to a standard regenerator.
The bulk temperature of the preheated air is 1018 °C with the WGR system; an increase of about 50 °C is obtained with respect to the standard regeneration system. The above effects confirm that the WGR system is effective not only for the NOx reduction process (as confirmed by the experimental tests [19]) but also for enhancing the thermal efficiency of the thermal regenerative process.

5. Conclusions

A one-dimensional transient model for the analysis of regeneration systems for glass furnaces is presented. This lower order model that allows a quick simulation of the regenerative set performance uses heat transfer coefficient profiles obtained from the application of specific CFD models. The heat transfer coefficients are locally re-computed in the code, taking into account the local composition of the gas according to the WGR strategy. A gas emissivity routine previously developed by the authors is used for the above heat transfer model. The model has shown its ability to simulate the effects of regenerator geometry or inversion time change. It can therefore be routinely used for the optimization of the refractory material distribution along the axis of the regenerative chambers during the system design. The model, extended to take into account the WGR effects, is applied to compare the regenerative system performance change with respect to the standard design. The analyses show an increase of the thermal performance of the regenerator when the WGR system is activated (air outlet temperature increase of about 50 °C in the example). The model proved to be a useful tool for understanding and quantifying the additional thermal effect introduced with the addition of radiant gases into the chamber when operating in the air phase according to the WGR strategy.

Author Contributions

Data curation, Alessadro Spoladore; Funding acquisition: C.C.; Methodology, C.C.; Software, A.S.; Supervision, C.C.

Acknowledgments

This model and the CFD applications presented have been developed within the PRIMEGLASS European Project (LIFE12 ENV/IT/001020) framework by the research team from the University of Genova.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Variables
aMatrix [a] coefficients
AArea [m2]
bMatrix [b] coefficients
cMatrix [c] coefficients
CpSpecific heat capacity at const. pressure [J/kg·K]
CvSpecific heat capacity at const. volume [J/kg·K]
CUCapacitance utilization
DDiameter [m]
F f Fluid fraction
HRegenerator height [m]
hConvective heat transfer coefficient [W/m2·K]
mMass [kg]
Mass flow rate [kg/s]
MMMolar mass [g/mol]
NNumber of cells
RStorage effectiveness
R f Specific gas constant [J/kg·K]
TTemperature [°C]
ththickness [m]
UThermal transmittance [W/m2·K]
VVolume [m3]
w j Mass concentration
ΕSpecific effectiveness
η T Thermal efficiency
ΡDensity [kg/m3]
ΤRegenerator inversion time [s]
Subscripts
airAir
chChannel
exhExhaust gas
fFluid
inInlet
maxMaximum
outOutlet
sSolid
sfSolid-fluid
ssSolid-solid
Abbreviations
CCCombustion chamber
CFDComputational fluid dynamic
EPEnd Port
RCRegenerative chamber
WGRWaste Gas Recirculation

References

  1. Stirling, R. Stirling Air Engine and the Heat Regenerator. Patent No. 4081, 1816. [Google Scholar]
  2. Siemens, C.W. XXXIII.—On the regenerative gas furnace as applied to the manufacture of cast steel. J. Chem. Soc. 1868, 21, 279–310. [Google Scholar] [CrossRef]
  3. Willmott, A.J. The development of thermal regenerator theory: 1931–the present. J. Inst. Energy 1993, 54. [Google Scholar]
  4. Hausen, H. Über die Theorie des Wärmeaustausches in Regeneratoren. The Theory of Heat Exchange in Regenerators. ZAMM-J. Appl. Math. Mech. 1929, 9, 173–200. [Google Scholar] [CrossRef]
  5. Rummel, K. The calculation of the thermal characteristics of regenerators. J. Inst. Fuel 1931, 3, 160–174. [Google Scholar]
  6. Lambertson, T.J. Performance Factors of a Periodic-Flow Heat Exchanger. Ph.D. Dissertation, Naval Postgraduate School, Monterey, CA, USA, 1957. [Google Scholar]
  7. Willmott, A.J. Digital computer simulation of a thermal regenerator. Int. J. Heat Mass Transf. 1964, 7, 1291–1302. [Google Scholar] [CrossRef]
  8. Willmott, A.J. Simulation of a thermal regenerator under conditions of variable mass flow. Int. J. Heat Mass Transf. 1968, 11, 1105–1116. [Google Scholar] [CrossRef]
  9. Willmott, A.J.; Hinchcliffe, C. The effect of gas heat storage upon the performance of the thermal regenerator. Int. J. Heat Mass Transf. 1976, 19, 821–826. [Google Scholar] [CrossRef]
  10. Sardeshpande, V.; Gaintonde, U.N.; Banerjee, R. Model based energy benchmarking for glass furnace. Energy Convers. Manag. 2007, 48, 2718–2738. [Google Scholar] [CrossRef]
  11. Tapasa, K.; Jitwatcharakomol, T. Thermodynamic calculation of exploited heat used in glass melting furnace. Procedia Eng. 2012, 32, 969–975. [Google Scholar] [CrossRef]
  12. Sardeshpande, V.; Anthony, R.; Gaintonde, U.N.; Banerjee, R. Performance analysis for glass furnace regenerator. J. Appl. Energy 2011, 88, 4451–4458. [Google Scholar] [CrossRef]
  13. Sadrameli, S.M.; Zarrinehkafsh, M.T. Simulation of fixed bed regenerative heat exchangers for flue gas heat recovery. J. Appl. Therm. Eng. 2004, 24, 373–382. [Google Scholar]
  14. Reboussin, Y.; Fourmigue, J.F.; Marty, P.; Citti, O. A numerical approach for the study of glass furnace regenerators. J. Appl. Therm. Eng. 2005, 25, 2299–2320. [Google Scholar] [CrossRef]
  15. Cravero, C.; Marsano, D. Numerical simulation of regenerative chambers for glass production plants with a non-equilibrium heat transfer model. In Proceedings of the Conference on Energy, Enviroment, Ecosystem and Sustainable Development, Roma, Italy, 27–29 January 2017. [Google Scholar]
  16. Basso, D.; Cravero, C.; Reverberi, A.P.; Fabiano, B. CFD analysis of regenerative chambers for energy efficiency improvement in glass production plants. Energies 2015, 8, 8945–8961. [Google Scholar] [CrossRef]
  17. Cravero, C.; Marsano, D.; Spoladore, A. Numerical strategies for fluid-dynamic and heat transfer simulation for regenerative chambers in glass production plants. NAUN Int. J. Math. Models Methods Appl. Sci. 2017, 11, 82–87. [Google Scholar]
  18. Basso, D.; Briasco, G.; Carretta, M.; Cravero, C.; Mola, A. CFD simulation of regenerative chambers in glass industry to support the design process for thermal efficiency improvement. In Proceedings of the International CAE Conference, Pacengo del Garda, Italy, 27–28 October 2014. [Google Scholar]
  19. Available online: www.primeglass.it (accessed on 1 November 2017).
  20. Cogliandro, S.; Cravero, C.; Marini, M.; Spoladore, A. Simulation strategies for regenerative chambers in glass production plants with strategic exhaust gas recirculation system. IJHT Int. J. Heat Technol. 2017, 35, S449–S455. [Google Scholar] [CrossRef]
  21. Romie, F.E. Periodic thermal storage: The regenerator. ASME J. Heat Transf. 1979, 101, 726–731. [Google Scholar] [CrossRef]
Figure 1. End-port (EP) glass furnaces design.
Figure 1. End-port (EP) glass furnaces design.
Applsci 09 01496 g001
Figure 2. Sketch of a single chamber domain discretization.
Figure 2. Sketch of a single chamber domain discretization.
Applsci 09 01496 g002
Figure 3. Heat flow for the ith layer.
Figure 3. Heat flow for the ith layer.
Applsci 09 01496 g003
Figure 4. Average temperature profile across the regenerative chamber (RC) height (reference case).
Figure 4. Average temperature profile across the regenerative chamber (RC) height (reference case).
Applsci 09 01496 g004
Figure 5. Temperatures variation with time at the top and bottom of the RC (reference case).
Figure 5. Temperatures variation with time at the top and bottom of the RC (reference case).
Applsci 09 01496 g005
Figure 6. Thermal efficiency, storage effectiveness, and capacitance utilization averaged (τ = 20′).
Figure 6. Thermal efficiency, storage effectiveness, and capacitance utilization averaged (τ = 20′).
Applsci 09 01496 g006
Figure 7. Thermal efficiency variation for different RC sizing and inversion time.
Figure 7. Thermal efficiency variation for different RC sizing and inversion time.
Applsci 09 01496 g007
Figure 8. Capacitance of utilization variation for different RC sizing and inversion time.
Figure 8. Capacitance of utilization variation for different RC sizing and inversion time.
Applsci 09 01496 g008
Figure 9. Waste gas tracking in the Cold phase chamber (20% average waste gas fraction) [20].
Figure 9. Waste gas tracking in the Cold phase chamber (20% average waste gas fraction) [20].
Applsci 09 01496 g009
Figure 10. Oxygen mass fraction profile in the Cold phase chamber stacking (20% average waste gas fraction) [20].
Figure 10. Oxygen mass fraction profile in the Cold phase chamber stacking (20% average waste gas fraction) [20].
Applsci 09 01496 g010
Figure 11. Oxygen mass fraction and waste gas mass fraction distribution in the Cold phase chamber stacking (20% average waste gas fraction) [20].
Figure 11. Oxygen mass fraction and waste gas mass fraction distribution in the Cold phase chamber stacking (20% average waste gas fraction) [20].
Applsci 09 01496 g011
Figure 12. Total and radiant heat transfer coefficients for each Cold RC sub-section.
Figure 12. Total and radiant heat transfer coefficients for each Cold RC sub-section.
Applsci 09 01496 g012
Figure 13. Average temperature profiles across the Cold RC with the WGR system application.
Figure 13. Average temperature profiles across the Cold RC with the WGR system application.
Applsci 09 01496 g013
Table 1. Air and exhaust gases chemicals mass concentration example.
Table 1. Air and exhaust gases chemicals mass concentration example.
AirWaste Gas
w j [ kg j kg tot ] N20.7910.687
O20.2090.023
CO200.119
H2O00.171
Table 2. Geometrical parameters of simulated cases.
Table 2. Geometrical parameters of simulated cases.
Case ACase BCase C
A s f , t o t 29001930970m2
V t o t 10066.633.3m3
F f 0.70.70.7
Table 3. Parameters values of the Cold RC with the Waste Gas Recirculation (WGR) system.
Table 3. Parameters values of the Cold RC with the Waste Gas Recirculation (WGR) system.
TotS-1S-2S-3
Volume [m3]10033.333.333.3
Heat transfer surface [m2]2900966.7966.7966.7
Mass flow [kg/s]4.521.131.132.26
Waste gas mass fraction20%0%10%35%
Inlet temperature [C]205147186.5243.4
Table 4. Model output of the Cold RC with the WGR system.
Table 4. Model output of the Cold RC with the WGR system.
WGR offWGR on - Average WG Fraction 20%
-TotS-1S-2S-3
Mass flow [kg/s]4.524.521.131.132.26
WG mass fraction [%]02001035
Inlet temperature [C]150205147186243
Outlet temperature [C]968101898310521019
Thermal power P t h [MW]3.093.740.9450.9911.81
Specific Thermal power [kW/m2]1.0651.2900.9781.0251.872

Share and Cite

MDPI and ACS Style

Cravero, C.; Spoladore, A. Transient Numerical Simulation of Regenerative Systems with Waste Gas Recirculation Strategies in Glass Production Plant. Appl. Sci. 2019, 9, 1496. https://doi.org/10.3390/app9071496

AMA Style

Cravero C, Spoladore A. Transient Numerical Simulation of Regenerative Systems with Waste Gas Recirculation Strategies in Glass Production Plant. Applied Sciences. 2019; 9(7):1496. https://doi.org/10.3390/app9071496

Chicago/Turabian Style

Cravero, Carlo, and Alessandro Spoladore. 2019. "Transient Numerical Simulation of Regenerative Systems with Waste Gas Recirculation Strategies in Glass Production Plant" Applied Sciences 9, no. 7: 1496. https://doi.org/10.3390/app9071496

APA Style

Cravero, C., & Spoladore, A. (2019). Transient Numerical Simulation of Regenerative Systems with Waste Gas Recirculation Strategies in Glass Production Plant. Applied Sciences, 9(7), 1496. https://doi.org/10.3390/app9071496

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop