**Open Access Library Journal**

Vol.04 No.05(2017), Article ID:75878,7 pages

10.4236/oalib.1103546

An Absorption-Adsorption Apparatus for Gases Purification from SO_{2} in Power Plants

Christo Boyadjiev, Boyan Boyadjiev^{ }

Institute of Chemical Engineering, Bulgarian Academy of Sciences, Sofia, Bulgaria

Copyright © 2017 by authors and Open Access Library Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: March 21, 2017; Accepted: April 30, 2017; Published: May 3, 2017

ABSTRACT

Different companies propose methods and apparatuses for waste gases purification from SО_{2} using two-phase absorbent (CaCO_{3} suspension).The basic problem of the carbonate absorbents is that its chemical reaction with SO_{2} leads to CO_{2} emission, which leads to similar ecological problems. The use of synthetic anionites as adsorbents for gas purification from SO_{2} provides possibilities for adsorbent regeneration. In a proposed method, the waste gas purification is realized in two steps-physical absorption of SO_{2} with water in a counter-current column and chemical adsorption of
${\text{HSO}}_{3}^{-}$
from the water solution by synthetic anionite particles. The adsorbent regeneration is made with NH_{4}OH solution. The obtained (NH_{4})_{2}SO_{3} (NH_{4}HSO_{3}) is used (after reaction with HNO_{3}) for production of concentrated SO_{2} (gas) and NH_{4}NO_{3} (solution). The proposed method permits to be used the absorption columns, where are used the CaCO_{3} suspensions. In this paper, the efficiency of the waste gases purification from SО_{2} is increased in an absorption-adsorption apparatus, where the absorption is realized in co-current flows and the adsorption takes place in the flexible adsorbent. For this is proposed a new absorption-adsorption column apparatus with bubbling plates. A mathematical model of the absorption-adsorption process is presented, too.

**Subject Areas:**

Chemical Engineering & Technology, Industrial Engineering

**Keywords:**

Absorption, Adsorption, Gas Purification, SO_{2}, Power Plant, Convection- Diffusion Model, Average Concentration Model

1. Introduction

Different companies (Babcock & Wilcox Power Generation Group, Inc., Alstom Power Italy, Idreco-Insigma-Consortium) propose methods and apparatuses for waste gases purification from SО_{2} using two-phase absorbent (CaCO_{3} suspension) [1] . The basic problem of the carbonate absorbents is that its chemical reaction with SO_{2} lead to CO_{2} emission (every molecule of SO_{2} absorbed from the air is equivalent to a molecule of CO_{2} emitted in the air), because the ecological problems (greenhouse effects) of SO_{2} and CO_{2} are similar. The large quantity of by-products is a problem, too. Another drawback of these methods is the impossibility for regeneration of the absorbents.

The theoretical analysis [2] - [8] of the method and apparatus for waste gases purification from SО_{2} using two-phase absorbent (CaCO_{3} suspension) shows the possibility to use an absorption-adsorption method. The use of synthetic anionites (basic anion-exchange resins-R-OH form of Amberlite, Duolite, Kastel, Varion, Wofatit) as adsorbents [9] [10] [11] for gas purification from SO_{2} provides possibilities for adsorbent regeneration. In the proposed method [12] [13] the waste gas purification is realized in two steps-physical absorption of SO_{2} with water and chemical adsorption of
${\text{HSO}}_{3}^{-}$
from the water solution by synthetic anionite particles in a fixed bed adsorber. The adsorbent regeneration is made with NH_{4}OH solution. The obtained (NH_{4})_{2}SO_{3} (NH_{4}HSO_{3}) is used (after reaction with HNO_{3}) for production of concentrated SO_{2} (gas) and NH_{4}NO_{3} (solution). The proposed patent [12] makes it possible to create a waste-free technology for waste gases purification from sulfur dioxide by means of regenerable absorbent and adsorbent. The proposed method [13] permits to be used for the absorption columns, where are used the CaCO_{3} suspensions.

The efficiency of the waste gases purification from SО_{2} can be increased if the absorption is realized in co-current flows and the adsorption takes place in the flexible adsorbent. For this is proposed a new absorption-adsorption column apparatus with bubbling plates.

2. An Absorption-Adsorption Apparatus

In the proposed absorption-adsorption method for waste gases purification from sulfur dioxide [12] [13] , the absorption is realized in the counter-current absorber (where practical gas velocity does not exceed 5 m∙s^{−1}) and the adsorption is carried out in a fixed bed adsorber. The efficiency of the process can be increased if the absorption is realized in co-current flows and the adsorption takes place in the flexible adsorbent. For this it can use a new absorption- adsorption apparatus with bubbling plates. The bubbling of the gas at a plate through a layer of aqueous suspension of synthetic anionite allows an increasing of the gas velocity (reduction of the diameter of the absorption column), elimination of the adsorption column and carrying out the adsorption in a flexible adsorbent. The movement of the gas between the plates leads to mixing in the gas phase, which increases the absorption rate because the absorption of sulfur dioxide in water is limited by mass transfer in the gas phase.

The new absorption-adsorption apparatus with bubble plates [14] is shown in Figure 1. The gas enters tangentially into the column 1 through the inlet 2,

Figure 1. Absorption-adsorption apparatus.

passes through the distribution pipes 5, concentric bubble caps 6 of the plates 4 and exits the column through the outlet 3. The aqueous suspension of synthetic anionite enters the column 1 via the valve 11 and the pipes 7 and creates of plates 4 layer with a certain thickness. After saturation of the adsorbent with sulfur dioxide, the aqueous suspension is output from the column through the pipes 8 and valves 11 and enters in the system 9 for the regeneration of the adsorbent. The suspension of the regenerated adsorbent is removed from 9 by pump 10 is returned to the plates 4 in the column 1.

The operation of the absorption-adsorption apparatus is the following cyclic scheme:

1. Supply all the plates with the necessary amount (volume) of aqueous suspension of synthetic anionite.

2. Start the absorption-adsorption process and monitor the SO_{2} concentration of the exit 3.

3. When the increasing the SO_{2} concentration at the outlet of the gas 3 exceeds permissible limits, the suspension of the first (bottom) plate is transferred in the regeneration system 9 and the plate is loaded with a new (regenerated) suspension.

4. When the increasing the SO_{2} concentration at the outlet of the gas 3 exceeds permissible limits, the suspension of the second (next) plate is transferred in the regeneration system 9 and the plate is loaded with a new (regenerated) suspension.

5. The procedures are repeated until reaching the top plate, then starts again from the first plate.

3. Absorption-Adsorption Process Modeling

Tangentially supplying of the gas in the column [15] reduces the radial non-un- iformity of the velocity at the column cross-sectional area and the gas velocity is a constant
$\left({u}_{1}=const.\right)$
, practically. Under these conditions, the concentration of SO_{2} in the gas phase is changed only in the height of the column (
${c}_{11}={c}_{11}\left(t,x\right)$
), which increases the rate of mass transfer rate. On the other hand the gas bubbling creates an ideal mixing regime in the liquid phase and the concentration of SO_{2} in the liquid phase is
${c}_{12}={c}_{12}\left(t\right)$
. The concentration of SO_{2} in the solid phase (capillaries) is
${c}_{13}$
and the concentration of active sites in the adsorbent is
${c}_{23}$
.

The interphase mass transfer rate of SO_{2} from the gas to the liquid is
${k}_{0}\left({c}_{11}-\chi {c}_{12}\right)$
, while the liquid to the adsorbent is
${k}_{1}\left({c}_{12}-{c}_{13}\right)$
. The chemical reaction rate of the SO_{2} with the adsorbent is
$k\text{\hspace{0.17em}}{c}_{13}{c}_{23}$
.

The modeling of a non-stationary (as a result of the adsorbent saturation) absorption-adsorption process on the plate number
$n\text{\hspace{0.17em}}\left(n=1,\cdots ,N\right)$
, for gas purification from SO_{2} [14] , uses a combination of the physical absorption and chemical adsorption models [1] :

$\begin{array}{l}\frac{\partial {c}_{11}^{\left(n\right)}}{\partial t}+{u}_{1}\frac{\partial {c}_{11}^{\left(n\right)}}{\partial z}={D}_{1}\frac{{\partial}^{2}{c}_{11}^{\left(n\right)}}{\partial \text{\hspace{0.17em}}{z}^{2}}-{k}_{0}\left({c}_{11}^{\left(n\right)}-\chi {c}_{12}^{\left(n\right)}\right);\\ \frac{\text{d}{c}_{13}^{\left(n\right)}}{\text{d}t}={k}_{1}\left({c}_{12}^{\left(n\right)}-{c}_{13}^{\left(n\right)}\right)-k{c}_{13}^{\left(n\right)}{c}_{23}^{\left(n\right)};\text{\hspace{1em}}\frac{\text{d}{c}_{23}^{\left(n\right)}}{\text{d}t}=-k{c}_{13}^{\left(n\right)}{c}_{23}^{\left(n\right)};\\ t=0,\text{\hspace{1em}}{c}_{11}^{\left(n\right)}\equiv {c}_{11}^{0},\text{\hspace{1em}}{c}_{13}^{\left(n\right)}\equiv 0,\text{\hspace{1em}}{c}_{23}^{\left(n\right)}\equiv {c}_{23}^{0};\\ z=0,\text{}{c}_{11}^{\left(n\right)}\left(0\right)\equiv {c}_{11}^{\left(n-1\right)}\left(l\right),\text{\hspace{1em}}0\equiv {\left(\frac{\partial {c}_{11}^{\left(n\right)}}{\partial z}\right)}_{z=0},\text{\hspace{1em}}{c}_{11}^{\left(1\right)}\left(0\right)\equiv {c}_{11}^{0};\text{\hspace{1em}}n=1,\cdots ,N.\end{array}$ (1)

In (1)
$z=0$
is the entrance to the gas input of the plate number
$n\text{\hspace{0.17em}}\left(n=1,\cdots ,N\right)$
,
$l$
is the distance between the plates,
${D}_{1}$
is the diffusivity of SO_{2} in the gas phase.

The concentration of SO_{2} in the water of each plate is determined by the amount of absorbed SO_{2} (
${W}_{1}$
) and its distribution between liquid phase (
${W}_{2}$
), the solid phase (capillaries) (
${W}_{3}^{1}$
) and adsorbed on the surface of the capillaries is (
${W}_{3}^{2}$
):

${W}_{1}={V}_{1}\left[{c}_{11}^{\left(n\right)}\left(0\right)-{c}_{11}^{\left(n\right)}\left(l\right)\right],\text{\hspace{0.17em}}\text{\hspace{0.17em}}{W}_{2}={V}_{2}{c}_{12}^{\left(n\right)},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{W}_{3}^{1}={V}_{3}{c}_{13}^{\left(n\right)},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{W}_{3}^{2}={V}_{3}\left({c}_{23}^{0}-{c}_{23}^{\left(n\right)}\right),$ (2)

i.e.

${c}_{12}^{\left(n\right)}=\frac{{V}_{1}\left[{c}_{11}^{\left(n\right)}\left(0\right)-{c}_{11}^{\left(n\right)}\left(l\right)\right]-{V}_{3}{c}_{13}^{\left(n\right)}-{V}_{3}\left({c}_{23}^{0}-{c}_{23}^{\left(n\right)}\right)}{{V}_{2}},$ (3)

where ${V}_{1},{V}_{2},{V}_{3}$ are the gas, liquid and solid phase volumes on the plates.

4. Generalized Analysis

The use of dimensionless (generalized) variables [16] allows to make a qualitative analysis of the models (1), (3), where as characteristic scales are used the inlet and initial concentrations, the characteristic time ${t}_{0}$ (saturation time of the adsorbent) and the distance between the plates $l$ :

$T=\frac{t}{{t}_{0}},\text{\hspace{1em}}Z=\frac{z}{l},\text{\hspace{1em}}{C}_{11}=\frac{{c}_{11}}{{c}_{11}^{0}},\text{\hspace{1em}}{C}_{12}=\frac{{c}_{12}\chi}{{c}_{11}^{0}},\text{\hspace{1em}}{C}_{13}=\frac{{c}_{13}\chi}{{c}_{11}^{0}},\text{\hspace{1em}}{C}_{23}=\frac{{c}_{23}}{{c}_{23}^{0}}.$ (4)

When (4) is put into (1), (3), the model in generalized variables takes the form:

$\begin{array}{l}\gamma \frac{\partial {C}_{11}^{\left(n\right)}}{\partial T}+\frac{\partial {C}_{11}^{\left(n\right)}}{\partial Z}={\mathrm{Pe}}^{-1}\frac{{\text{d}}^{2}{C}_{11}^{\left(n\right)}}{\text{d}{Z}^{2}}-{K}_{0}\left({C}_{11}^{\left(n\right)}-{C}_{12}^{\left(n\right)}\right);\\ \frac{\text{d}{C}_{13}^{\left(n\right)}}{\text{d}T}={K}_{1}\left({C}_{12}^{\left(n\right)}-{C}_{13}^{\left(n\right)}\right)-K{C}_{13}^{\left(n\right)}{C}_{23}^{\left(n\right)};\text{\hspace{1em}}\frac{\text{d}{C}_{23}^{\left(n\right)}}{\text{d}T}=-K{\alpha}^{-1}{C}_{13}^{\left(n\right)}{C}_{23}^{\left(n\right)};\\ T=0,\text{\hspace{1em}}{C}_{11}^{\left(n\right)}\equiv 1,\text{\hspace{1em}}{C}_{13}^{\left(n\right)}\equiv 0,\text{\hspace{1em}}{C}_{23}^{\left(n\right)}\equiv 1.\\ Z=0,\text{}{C}_{11}^{\left(n\right)}\left(0\right)\equiv {C}_{11}^{\left(n-1\right)}\left(1\right),\text{\hspace{1em}}0\equiv {\left(\frac{\text{d}{C}_{11}^{\left(n\right)}}{\text{d}Z}\right)}_{Z=0},\text{\hspace{1em}}{C}_{11}^{\left(1\right)}\left(0\right)\equiv 1;\text{\hspace{1em}}n=1,\cdots ,N.\end{array}$ (5)

${C}_{12}^{\left(n\right)}=\frac{\chi {V}_{1}\left[{C}_{11}^{\left(n\right)}\left(0\right)-{C}_{11}^{\left(n\right)}\left(l\right)\right]-{V}_{3}{C}_{13}^{\left(n\right)}-\alpha {V}_{3}\left(1-{C}_{23}^{\left(n\right)}\right)}{{V}_{2}},\text{\hspace{1em}}n=1,\cdots ,N.$ (6)

The following parameters are used in (5), (6):

$\mathrm{Pe}=\frac{{u}_{1}l}{{D}_{1}},\text{\hspace{1em}}K=k{t}_{0}{c}_{23}^{0},\text{\hspace{1em}}{K}_{0}=\frac{{k}_{0}l}{{u}_{1}},\text{\hspace{1em}}{K}_{1}={k}_{1}{t}_{0},\text{\hspace{1em}}\alpha =\frac{\chi {c}_{23}^{0}}{{c}_{11}^{0}},\text{\hspace{1em}}\gamma =\frac{l}{{t}_{0}{u}_{1}}.$ (7)

Practically $0=\gamma <{10}^{-2},\text{\hspace{0.17em}}0={\mathrm{Pe}}^{-1}<{10}^{-2}$ and as a result

$\begin{array}{l}\frac{\text{d}{C}_{11}^{\left(n\right)}}{\text{d}Z}=-{K}_{0}\left({C}_{11}^{\left(n\right)}-{C}_{12}^{\left(n\right)}\right);\\ Z=0,\text{}{C}_{11}^{\left(n\right)}\left(0\right)\equiv {C}_{11}^{\left(n-1\right)}\left(1\right),\text{\hspace{1em}}{C}_{11}^{\left(1\right)}\left(0\right)\equiv 1;\text{\hspace{1em}}n=1,\cdots ,N.\end{array}$ (8)

$\begin{array}{l}\frac{\text{d}{C}_{13}^{\left(n\right)}}{\text{d}T}={K}_{1}\left({C}_{12}^{\left(n\right)}-{C}_{13}^{\left(n\right)}\right)-K{C}_{13}^{\left(n\right)}{C}_{23}^{\left(n\right)};\text{\hspace{1em}}\frac{\text{d}{C}_{23}^{\left(n\right)}}{\text{d}T}=-K\alpha {C}_{13}^{\left(n\right)}{C}_{23}^{\left(n\right)};\\ T=0,\text{\hspace{1em}}{C}_{13}^{\left(n\right)}\equiv 0,\text{\hspace{1em}}{C}_{23}^{\left(n\right)}\equiv 1.\end{array}$ (9)

The solution of the equations of the model (6), (8), (9) uses a two-stage algorithm. In the first stage must be solved the equations for $n=1$ . In the second stage must be applied consistently for every plate the algorithm for $n=1$ .

Algorithm of the solution

1. Put $n=1$ .

2. Put ${C}_{12}^{\left(1\right)}={X}_{i}=0.1i,\text{\hspace{0.17em}}i=1,\cdots ,10.$

3. The solution of leads to

${C}_{11}^{\left(1i\right)}={C}_{11}^{\left(1\right)}\left(Z,{X}_{i}\right),\text{\hspace{0.17em}}{C}_{11}^{\left(1\right)}\left(0,{X}_{i}\right),\text{\hspace{0.17em}}{C}_{11}^{\left(1\right)}\left(1,{X}_{i}\right),\text{\hspace{0.17em}}i=1,\cdots ,10.$

4. The solution of (9) leads to ${C}_{13}^{\left(1i\right)}={C}_{13}^{\left(1\right)}\left(T,{X}_{i}\right),\text{\hspace{0.17em}}{C}_{23}^{\left(1i\right)}={C}_{23}^{\left(1\right)}\left(T,{X}_{i}\right),\text{\hspace{0.17em}}i=1,\cdots ,10.$

5. The solutions in 3 and 4 must be introduced in (6) and as a result is obtained ${C}_{12}^{\left(1i\right)}={C}_{12}^{\left(1\right)}\left(T,{X}_{i}\right),\text{\hspace{0.17em}}i=1,\cdots ,10.$

6. Put $T={T}_{i}=0.1j,\text{\hspace{0.17em}}j=1,\cdots ,10$ in ${C}_{12}^{\left(1i\right)}={C}_{12}^{\left(1\right)}\left(T,{X}_{i}\right),\text{\hspace{0.17em}}i=1,\cdots ,10$ an as a result is obtained ${\overline{C}}_{12}^{\left(1i\right)}={C}_{12}^{\left(1\right)}\left({T}_{j},{X}_{i}\right),\text{\hspace{0.17em}}i=1,\cdots ,10,\text{\hspace{0.17em}}j=1,\cdots ,10.$

7. A polynomial approximation ${P}_{12}^{\left(1j\right)}\left({T}_{j},X\right),\text{\hspace{0.17em}}j=1,\cdots ,10$ of ${\overline{C}}_{12}^{\left(1i\right)}={C}_{12}^{\left(1\right)}\left({T}_{j},{X}_{i}\right),$ $i=1,\cdots ,10,\text{\hspace{0.17em}}j=1,\cdots ,10$ with respect to $X$ must be obtained.

8. The solutions of the equations $X={P}_{12}^{\left(1j\right)}\left({T}_{j},X\right),\text{\hspace{0.17em}}j=1,\cdots ,10$ with respect to $X$ permit to be obtained ${X}_{j},\text{\hspace{0.17em}}j=1,\cdots ,10$ and the obtained solutions must be denoted as ${X}_{j}={C}_{12}^{\left(1\right)}\left({T}_{j}\right),\text{\hspace{0.17em}}j=1,\cdots ,10.$

9. A polynomial approximation of ${C}_{12}^{\left(1\right)}\left({T}_{j}\right),\text{\hspace{0.17em}}j=1,\cdots ,10$ with respect to $T$ permit to be obtained ${C}_{12}^{\left(1\right)}\left(T\right)={C}_{12}^{\left(1\right)}\left({T}_{j}\right),\text{\hspace{0.17em}}j=1,\cdots ,10.$

10. The introducing of ${C}_{12}^{\left(1\right)}\left(T\right)$ in (8) and (9) permits to be obtained its solutions for $n=1$ .

11. The obtained solution of (8) ${C}_{11}^{\left(1\right)}={C}_{11}^{\left(1\right)}\left(T,Z\right)$ in 10 permits to be obtained ${C}_{11}^{\left(1\right)}={C}_{11}^{\left(1\right)}\left(T,1\right)$ and as a result to be used the algorithm 1-11 for consistent solutions of the equations set (6), (8) for $n=2,\cdots ,N.$

5. Parameters Identification

The parameters in the model (6), (8), (9), which are subject to experimental determination are
$K,\text{\hspace{0.17em}}{K}_{0},{K}_{1}$
. They may be obtained from experimental data of the SO_{2} concentration at the gas outlet from the first plate
${C}_{11\mathrm{exp}}^{\left(1j\right)}={C}_{11}^{\left(1\right)}\left({T}_{j},1\right),\text{\hspace{0.17em}}{T}_{i}=0.1j,\text{\hspace{0.17em}}j=1,\cdots ,10$
, where
$T=1$
is the time for the fully saturation of the adsorbent on the first plate. For this purpose must be minimized the function of the least squares with respect to
$K,\text{\hspace{0.17em}}{K}_{0},{K}_{1}$
:

$F\left(K,{K}_{0},{K}_{1}\right)={\displaystyle \sum _{j=1}^{10}{\left[{C}_{11}^{\left(1\right)}\left({T}_{j},1\right)-{\overline{C}}_{11\mathrm{exp}}^{\left(1j\right)}\right]}^{2}}.$ (10)

6. Conclusion

The proposed utility model [14] uses an absorption-adsorption column and gives the possibility to create a waste-free technology for waste gases purification from sulfur dioxide by means of regenerable adsorbent, where the absorbent regeneration system is similar to the regeneration system in the patent [12] [13] . The efficiency of the processes is increased in an absorption-adsorption apparatus, where the absorption is realized in co-current flows and the adsorption takes place in the flexible adsorbent. That’s why a new absorption-adsorption column apparatus with bubbling plates is offered. A mathematical model of the absorption-adsorption process in a plate column is presented too.

Cite this paper

Boyadjiev, C. and Boyadjiev, B. (2017) An Absorption- Adsorption Apparatus for Gases Purifica- tion from SO_{2} in Power Plants. Open Access Library Journal, 4: e3546. https://doi.org/10.4236/oalib.1103546

References

- 1. Boyadjiev, Chr., Doichinova, M., Boyadjiev, B. and Popova-Krumova, P. (2016) Modeling of Column Apparatus Processes. Springer-Verlag, Berlin Heidelberg.
- 2. Boyadjiev, Chr. (2011) Mechanism of Gas Absorption with Two-Phase Absorbents. International Journal of Heat and Mass Transfer, 54, 3004-3008.

https://doi.org/10.1016/j.ijheatmasstransfer.2011.02.050 - 3. Boyadjiev, Chr. (2011) On the SO2 Problem in Power Engineering. Energy Forum, Bulgaria, 114-125.
- 4. Boyadjiev, Chr. (2012) On the SO2 Problem in Power Engineering. Proceedings, Asia-Pacific Power and Energy Engineering Conference (APPEEC 2012), 1.
- 5. Boyadjiev, Chr. (2012) On the SO2 Problem of Solid Fuel Combustion. VIII All-Russian Conference with International Participation “Solid Fuel Combustion”, Novosibirsk, 13-16 November 2012.
- 6. Boyadjiev, Chr., Doichinova. M. and Popova, P. (2011) On the SO2 Problem in Power Engineering. 1. Gas Absorption. 15th Workshop on Transport Phenomena in Two-Phase Flow. Sunny Beach Resort, Bulgaria, 94-103.
- 7. Boyadjiev, Chr., Popova, P. and Doichinova, M. (2011) On the SO2 Problem in Power Engineering. 2. Two-Phase Absorbents. 15th Workshop on Transport Phenomena in Two-Phase Flow, Bulgaria, 104-115.
- 8. Boyadjiev, Chr., Doichinova, M. and Popova, P. (2012) On the SO2 Problem in Power Engineering. Transactions of Academenergo, No. 2, 44-65.
- 9. Pantofchieva, L. and Boyadjiev, Chr. (1995) Adsorption of Sulphur Dioxide by Synthetic Anion Exchangers. Bulgarian Chemical Communications, 28, 780-794.
- 10. Boyadjiev, Chr., Pantofchieva, L. and Hristov, J. (2000) Sulphur Dioxide Adsorption in a Fixed Bed of a Synthetic Anionite. Theoretical Foundations of Chemical Engineering, 34, 141-144.

https://doi.org/10.1007/BF02757831 - 11. Hristov, J., Boyadjiev, Chr. and Pantofchieva, L. (2000) Sulphur Dioxide Adsorption in a Magnetically Stabilized Bed of a Synthetic Anionite. Theoretical Foundations of Chemical Engineering, 34, 439-443.

https://doi.org/10.1007/BF02827387 - 12. Boyadjiev, Chr., Boyadjiev, B., Doichinova, M. and Popova-Krumova, P. (2014) Method and Apparatus for Gas Cleaning from Sulfur Dioxide. Bulgarian Patent Application No. 111398.
- 13. Boyadjiev, Chr. and Boyadjiev, B. (2017) An Absorption-Adsorption Method for Gases Purification from SO2 in Power Plants. Recent Innovations in Chemical Engineering. (In Press)
- 14. Boyadjiev, Chr. and Boyadjiev, B. (2015) Absorption-Adsorption Apparatus for Gas Cleaning from Sulfur Dioxide. Bulgarian Utility Model BG 2196 U1.
- 15. Boyadjiev, Chr. and Boyadjiev, B. (2013) Column Reactor for Chemical Processes. Bulgarian Utility Model BG 1776 U1.
- 16. Boyadjiev, Chr. (2010) Theoretical Chemical Engineering. Modeling and Simulation. Springer-Verlag, Berlin Heidelberg.