Next Article in Journal
Optimization of Cogeneration Power-Desalination Plants
Previous Article in Journal
Fault Diagnosis of Induction Motor Using D-Q Simplified Model and Parity Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Ammonia Concentration on Solvay Soda Process Parameters and Associated Environmental and Energetic Effects

1
Department of Chemical Technology, Faculty of Chemistry, Nicolaus Copernicus University in Toruń, 7 Gagarin Street, 87-100 Toruń, Poland
2
MCMP Sp. z o.o., 5 Świerkowa Street, 86-300 Grudziądz, Poland
3
Faculty of Civil Engineering and Environmental Sciences, Bialystok University of Technology, 45A Wiejska Street, 15-351 Białystok, Poland
4
CIECH R&D Sp. z o.o., 62 Wspólna Street, 00-684 Warszawa, Poland
*
Author to whom correspondence should be addressed.
Energies 2022, 15(22), 8370; https://doi.org/10.3390/en15228370
Submission received: 3 October 2022 / Revised: 3 November 2022 / Accepted: 6 November 2022 / Published: 9 November 2022
(This article belongs to the Section B: Energy and Environment)

Abstract

:
Modifying the absorption process in soda production by the Solvay method requires performing many calculations and determining a new equilibrium process. An increase in ammonia concentration in the reaction solution causes kinetic changes in equilibrium. Changes to the Solvay soda production technology were determined using chemical and instrumental analysis methods. A modification of the process in the form of SAB was introduced. Information allowing the design of an additional absorber and its location in the network of technological devices was presented in the form of parameters using typical chemical engineering assumptions. Spectroscopic and electrochemical techniques were used for this purpose. The increase in total alkalinity due to the addition of ammonia to 135 mmol·20 cm−3 resulted in cooling savings of about 152.4 MJ·Mg−1 of soda. The ammonia desorption rate and process energy parameters were determined for the new system. The temperature requirements for the carbonation column were defined, and in particular, a technique was developed to minimize the cooling of the lower part of the reactor, which reduces the consumption of process energy. Emissions of CO2 were reduced from 11.70 to 7.85% and NH3 from 5.52 to 4.89% in exhaust gases from the carbonation column.

Graphical Abstract

1. Introduction

The long-term development of the soda industry worldwide has led to a number of modifications to the Solvay process. The main one is the dual process, which allows manufacturing units to produce ammonium chloride (which can be used as a fertilizer) in almost equal amounts as the product. Several plants in the world operate in this mode; most of them are located in China. The Akzo dry lime process, which uses lime instead of lime milk to recover ammonia, is currently considered the most advanced technology. There are technologies in which soda is produced as a by-product [1]. Energy, water consumption, and environmental pollution are a difficult problem in the world. Due to the rapid expansion of technology, the demand for water and energy is increasing, which also enhances the dependence on the availability of these sources [2,3]. The Solvay soda production process is very energy-intensive [4] and is characterized by the large production of solid, liquid, and gaseous waste [5,6,7]. In many countries worldwide, especially in developing countries, severe water shortages have been caused by population growth, industrialization, global warming, and climate change. Water is the most important resource on Earth, and is essential to human life. It is estimated that approximately 1.8 billion people will be without access to clean drinking water by 2025. The proposed solution in the Solvay soda industry allows for large savings in the consumption of cooling water and energy [8]. There are several energy-consuming processes in the Solvay soda technology: distillation, carbonation (cooling), calcination, and calcination of limestone.
The main issue is the efficiency and energy of the carbonation process of ammoniated brine, especially based on seasonality. In summer, cooling the column is a big challenge. Considering the increasingly higher temperatures and availability of cooling water, this is one of the main economic factors in the operation of a carbonation column. To increase the efficiency of the process when there is less cooling of the column, studies on increasing the concentration of ammonia in the carbonation process were carried out [9].
One of the diffusion operations in the soda production process using the Solvay technique is the absorption of ammonia. Permeation, or mass exchange, is the basis of many unit operations, including distillation, rectification, drying, hydration, absorption, adsorption, crystallization, dissolution, etc. [10]. There is a gap in the research in terms of determining the influence of the ammonia concentration on the efficiency of the soda process under the real conditions of the Solvay process.
Absorption is the uptake of a component or components of a gaseous mixture by a liquid, which is the absorbent. The gas, due to diffusion, penetrates into the liquid through the interfacial surface to form a solution. A gaseous mixture can contain components that are soluble and components that are practically insoluble in the liquid. The former are called active components, the latter neutral or inert [11,12]. This is one of the main processes of Solvay soda technology.
The statics of the absorption process, i.e., the equilibrium between the liquid and the gas phase, depend on the composition of the phases, temperature and pressure, since the two-phase system under consideration has three degrees of freedom according to Gibbs’s phase rule [10]. Therefore, at a given temperature and pressure, a fixed composition of one phase corresponds to a well-defined composition of the other phase. Absorption will occur when the content of the component in the gas phase is greater than the corresponding equilibrium state.
Empirical formulas are known in the literature that allow for the approximate determination of the selected parameters of the Solvay process, but their technological usefulness has not been confirmed. Previous research has shown that ensuring the maximum value of the degree of transformation of the sodium ion WNa requires the particularly careful operation of crude brine purification, ammonization, and carbonation.
From the analysis of dependencies, we can formulate:
[ N H 3 ,   b o u n d ] = 68.0 + ( c 95.0 ) + 0.56 ( ( o + e + a + m ) 92.0 ) + 0.57 ( R 190.0 ) 0.14 ( T 25.0 )
and
W N a = a 1 · R + a 2 · [ N H 3 ,   t o t a l ] + a 3 · [ C l ] t o t a l a 4 · T
It is apparent from the analysis that the figure in question is directly proportional to the value of the degree of carbonation (R) and the concentration of ammonia ( o + e + a + m ) and total sodium chloride (c) in solution (a1–a4, numerical coeff.). An increase in the temperature of ammoniated brine at the exit from the carbonation column reduces the efficiency of the carbonation process [9].
Previous industrial experience indicates that maintaining the optimal level of process efficiency primarily depends on the magnitude of ammonia desorption in the carbonation apparatus, which is associated with temperature–hydrodynamic conditions and the distribution of the partial pressure of NH3 over the carbonated ammonia brine solution.
When considering the equation determining the efficiency of the soda process, it should be noted that only one quantity can be changed, taking into account the other parameters. The concentration of chloride ions is the maximum concentration that occurs in saturated brine. The degree of carbonation depends on the amount of CO2 introduced into the solution and is limited by the amount of CaCO3 decomposition. A change in the concentration of total ammonia remains [13].
The kinetics of absorption, that is, the rate of mass exchange, are determined by how far the system is from equilibrium. Of course, the statics, as well as the kinetics, further depend on the properties of the absorbing liquid, the absorbed component, and the inert gas. The kinetics of the absorption process are also affected by the phase contact conditions, which are related to the design of the absorption device [14].
Absorption is accompanied by a thermal effect caused by a change in state of aggregation [15]. Isothermal absorption is possible when the amount of liquid supplied is large compared to the amount of gas (low solubility of the gas) or when the heat of absorption is properly dissipated from the absorption apparatus; otherwise, the absorption will proceed with a change in the temperature of the system. In many cases, it is difficult to delineate a strict boundary between physical absorption and absorption associated with a chemical reaction [16]. As part of ongoing research at the Inowrocław production plant, it was decided to introduce a new “small” absorber (SAB) to increase the concentration of ammonia in the ammonia brine, pre-carbonated. The optimal operation of the additional absorption process translates into higher utilization of CO2 and the increased efficiency of the soda process.
First, the input and output parameters of the media for SAB were determined. Second, the number of mass transfer units was calculated, with one mass transfer unit corresponding to a section of the apparatus length over which the change in working concentrations of absorbate equals the average driving force over that section. The absorption process is an exothermic process [14,15].
The height of the mass transfer unit is calculated by the following formula, based on [10]:
h y = V K v · F 0
The goal of the work carried out was to obtain a stable flux of ammoniated brine, characterized by a direct titer of about 130 mmol∙20 cm−3 and determine the chloride ion concentration and CO2 before and after the absorption process.

2. SAB Calculation and Parameters in Modified Process

The main change in the process is the addition of ammonia absorber. This is an essential part of the modification to achieve the stated objectives of this paper. The input parameters of the ammonia brine were as follows: the absorber load was assumed at a level corresponding to 20–30% of the current load of the carbonation column, Vmax = 15–20 m3∙h−1, C NH 3 ≈ 100 mmol∙20 cm−3 ≈ 5.00 mol·dm−3, CNaCl ≈ 88 mmol∙20 cm−3 ≈ 4.40 mol·dm−3, T ≈ 313 K. The output parameters of ammonia brine with higher ammonia concentration were as follows: Vmax = 15–20 m3∙h−1, C NH 3 ≥ 130 mmol∙20 cm−3 ≈ 6.50 mol·dm−3, CNaCl ≈ 83 mmol∙20 cm−3 ≈ 4.15 mol·dm−3, T ≈ 333K.
A general diagram of the proposed absorber to increase the concentration of ammonia in ammonia brine, pre-carbonated, is shown in Figure 1.
An analysis of the adsorbed gas was performed using FT-IR spectroscopy (Vertex 70V FT-IR spectrometer, Bruker Optics, Ettlingen, Germany). The content of 1 Nm3 of the analyzed gas was 0.6482 Nm3 NH3, 0.1804 Nm3 CO2, and 0.1714 Nm3 H2O. The ammonia content of the input gas (bottom of the absorber) was Y d o w n = 1 . 8425   mol   NH 3 inerts   mol . It was assumed that 90% of ammonia would be absorbed; therefore, Y u p   90 = 0 . 18425   mol   NH 3 inerts   mol . The flow rate of the carbonated brine was set at 15–20 m3·h−1, that is, 16,935.0–22,580.0 kg·h−1. The density of the ammonia brine (with values of 100.76 mmol∙20 cm−3 NH3 and 89.93 mmol∙20 cm−3 Cl) was 1.1290 kg·dm−1.
Identifying the solution included determining direct alkalinity, total alkalinity (n), the amount of CO2 (d), and the content of chloride ions (Cl) (c) and ammonium chloride (NH4Cl) (b).
The method involves distilling off the ammonia and absorbing it in a standard solution of sulfuric(VI) acid. Excess sulfuric(VI) acid is titrated with a standard solution of sodium hydroxide in the presence of methyl orange. Ammonia contained in liquids in the form of carbonate and bicarbonate salts passes into the gas phase as a result of thermal decomposition of these compounds; this is called free ammonia.
Bound ammonia in the form of NH4Cl decomposes with a strong alkali such as sodium hydroxide.

2.1. Direct Alkalinity (n)

To determine alkalinity, the sample was titrated with 0.5 M sulfuric(VI) acid in the presence of methyl orange. The titration was carried out until the color of the solution changed from yellow to orange-yellow.

2.2. Amount of CO2 (d)

The determination was performed on a Scheibler apparatus. The carbonate and bicarbonate ions contained in the sample were decomposed using hydrochloric acid, with the release of free CO2, which displaces the liquid from the gas burette tube so that the amount of gas released can be determined.

2.3. Content of Cl (c)

The determination of chloride ions was performed using a Metrohm Ti-Touch 916 compact potentiometric titrator. To perform the measurement, the sample was neutralized, then diluted and acidified with sulfuric(VI) acid. The sample was then titrated with 0.03333 M silver nitrate solution in the presence of a silver electrode, using the appropriate instrument software.

2.4. Content of NH4Cl (b)

The potentiometric titrator was used to determine the amount of ammonium chloride in the sample. To carry out the measurement, the sample was boiled in lye; this process involved adding the sample to a specific amount of sodium hydroxide and then leaving it to stand on a hot plate to boil off the ammonia. The remaining excess alkali was titrated with 1 M sulfuric(VI) acid in the presence of a glass electrode, using the appropriate program of the instrument.

2.5. Content of NH2COONH4 (e)

The presence of carbamate ions was confirmed by 13C NMR (Avance 300NMR spectrometer, Bruker, Germany), and the ion content was analyzed using IR spectrometry at a wavelength of 687 cm−1.
All analyses were performed using the standards and methodologies presented in [17].
For the absorption process, based on Whitman–Lewis boundary layer theory and Fick’s first law [18], the mass exchange equation can be formulated as:
G = K · F · Δ y
The driving force of the mass exchange process can be expressed in terms of difference:
Δ y = ( y y )
As the absorption process proceeds, the driving force of the mass exchange process changes along the absorption apparatus (i.e., the distance of the system from equilibrium changes). For a sufficiently small mass exchange area dF, the mass exchange equation can be written as:
d G = K · d F · ( y y )
From the material balance, it is known that:
G = V · ( y 1 y 2 )   or   d G = V d y
From the comparison of the above quantities we obtain:
d G = K · d F · ( y y ) = V d y
Considering V from the material balance equation, we obtain:
K · d F · ( y y ) = G y 1 y 2 d y
Transforming the equation and integrating the limits from 0 to F and from y2 to y1, we obtain:
0 F d F = G y 1 y 2 y 2 y 1 d y K ( y y )
Therefore, we can obtain:
G = F · y 1 y 2 y 2 y 1 d y K · ( y y ) = K · F · Δ y s r
where:
Δ y s r = y 1 y 2 y 2 y 1 d y y y
Most often, however, strict determination of the size of the surface area of the separation of phases is not possible; therefore, in the formulas used for practical calculations, this size is replaced with a readily available volume size. Assuming that the magnitude of the phase separation area is proportional to the volume, we can calculate:
F = α · F 0 · H w   or   d F = α · F 0 · d H w
Eventually, combining the above equations with the basic differential equation, we obtain:
d H w = V d y K v · F 0 · ( y y )
Integrating the above relationship from 0 to Hw and from y2 to y1, we obtain:
H w = V K v · F 0 y 2 y 1 d y y y   or ,   in   short ,   H w = h y m y
As an example, it is useful to present the process of calculating the absorber for water. The solubility of ammonia in water at 313 K at a pressure of 1 atm assuming an isothermal process is shown in Table 1.
We can graph ( Y down Y up ) ( X down X up ) , where Xup  = 1.8425 0.18425 ( 6.5 5.0 ) = 1.65825 1.5 = 1.1055 × 1.2 = 1.3266 ; hence, Xdown = 0.2217. By drawing a graph of mol   NH 3 mol   of   inert = f ( mol   NH 3 mol   H 2 O ) , the number of mass transfer units can be determined graphically, as shown in Figure 2.
The column fill height is the product of the number and height of mass transfer units:
H w = h y m y
The number of mass transfer units is 0.2217.
We assume that the height of mass penetration for the liquid phase is:
h y = L k x a S
From here, the parameters of the absorber can be easily calculated.

3. Experimental Results

In order to effectively increase the concentration of ammonia in the pre-carbonated brine, it is necessary to determine the effect of the changes on the solubility of NaCl in the pre-carbonated brine so as not to cause NaCl precipitation.
In order to obtain information on the solubility of NaCl, equilibrium solubility tests were conducted on sodium chloride at 303 K in ammonia solution with equilibrium concentrations of 118 and 106.5 mmol∙20 cm−3 (initial concentrations 136.5 and 121.6 mmol∙20 cm−3, respectively), depending on the degree of carbonation of the solution (R). The results of the experiments are shown in Figure 3 and Figure 4.
The graphs in Figure 3 and Figure 4 show that the solubility of NaCl decreases with increasing NH3 in solution, and simultaneously increases as the degree of carbonation of the solution rises. At the same time, the trend of increasing solubility with increasing R is more noticeable. That is, considering the analyzed system, with an increase in the degree of carbonation, the solubility of NaCl increases, and this increase is more characteristic when R is in the range of 42–46%.
The traditional Solvay process absorption system implemented at CIECH Soda Polska S.A.’s Inowrocław soda plants is shown in Figure 5.
The process of making ammonia brine containing about 4.8 mol NH3·dm−3 in a saturating apparatus can be carried out with an ammonia concentration of up to 6.5 mol·dm−3. In view of the above, ammonia gases from the distillation process should be directed as a gas phase to a small countercurrent absorber with a cooler (SAB). Pre-carbonated ammonia brine is the absorbing liquid. The media are contacted in countercurrent. The cooled liquid after absorption is dispensed into the selected area of the carbonation column. A diagram of how the task is carried out is shown in Figure 6.
The choice of dosing sites for enriched ammonia brine to Cl is a matter of research, which will be published in a subsequent paper. To properly characterize the procedure, the process of ammonia desorption from the system must also be analyzed. Ammonia desorption measurements were carried out to identify the ammonia loss in the sodium process for different ammonia concentrations in the starting solution, in an open system using our apparatus, shown in Figure 7.
The apparatus made it possible to perform the tests at a strictly assumed volume rate of inert qin, flowing from a buffer tank through a thermostated reactor containing a well-defined volume of ammonia brine of known chemical composition. The pressure and gas flow control block allowed full control of hydrodynamic parameters during the testing of the ammonia desorption process under intensive barbotage conditions in the reactor. The scrubber system, with a standard solution of sulfuric acid(VI) and sodium lye, provided measurements of the amount of ammonia and carbon dioxide in the effluent gas.
In order to determine the effect of the concentration of ammonia in the solution and the degree of carbonation on the intensity of the process under consideration, six solutions of ammonia brine with various chemical compositions were prepared: C NH 3 = 4.404 6.258 mol · dm 3 , C NaCl = 4.147 4.267 mol · dm 3 , and C CO 2 = 0.000 1.966 mol · dm 3 .
An investigation of the ammonia desorption process at 303, 313, and 323 K was performed, analyzing the change in concentrations of ammonia, carbon dioxide, and sodium chloride in the liquid phase at specific time intervals.
Table 2 and Table 3 list the measurement results of two series of tests of the kinetics of the ammonia desorption process with the extreme concentration of NH3 in the brine.
In the following, we summarize the chemical composition of the ammonia brine entering the reactor, the change in concentration of total ammonia in solution, the temperature of the process, the value of the volumetric flow rate of the inert gas, and the time of ammonia desorption.
The effect of individual process parameters on the studied ammonia desorption was evaluated by function analysis:
K N H 3 d e s = f ( T ,   R , c N H 3 c O , q i n )
Within the studied range of temperature changes, ammonia brine chemistry, and inert volume flow rate, the desorption intensity of NH3 was approximated by the equation:
N N H 3 d e s = K N H 3 d e s · c N H 3 c
The experimental data of the desorption process confirm the assumptions of the kinetics of this process. The higher the temperature, the greater the desorption. Because the stream of desorbed ammonia is not linear, experimental studies of the desorption rate were conducted. The desorption tests confirmed that in order for additional ammonia to be bound, it must be dosed to the zone of the carbonation column, where the total CO2 concentration is as high as possible. As a result, the desorption time will be the longest and the concentration of carbamate ions the highest. Research on this process was carried out at the Faculty of Chemistry of the Nicolaus Copernicus University in Toruń.

4. Solution Implementation and Energy Changes

In order to implement the ammonization process of ammonia brine and pre-carbonated ammonia brine using the SAB plant, a number of technical, technological, and analytical experiments were carried out. In the course of the work, a well-defined flux of ammoniated brine was obtained with direct alkalinity of 130–140 mmol∙20 cm−3 (i.e., 6.50–7.00 mol∙dm−3).
In the course of the work, favorable conditions were created to carry out the research on a higher concentration of ammonia in the brine, and thus a higher concentration of CO2 and the co-crystallization of NH4HCO3 or NaCl. The absorption of ammonia in aqueous solution can be represented by the following equation:
NH3(g) + H2O(c) ⇌ NH3∙H2O(r) − 35.9 kJ
The absorption of carbon dioxide in an aqueous ammonia solution can be presented as follows:
2NH(r) + CO2(g) + H2O(c) ⇌ (NH4)2CO3(r) − 91.7 kJ
Ammonia is a gas that is easily soluble in water. In this case, the absorption of the gas in water is affected only by the resistance of the gas boundary layer. The rate of this process can be expressed by the formula:
N = Kg S (Pg − Pr) (kmol·h−1·m−2)
The expression (Pg − Pr) is referred to as the driving force of the absorption process, or absorption potential. This potential depends on the partial pressure difference of the absorbate (Pg) in the feed gas and its equilibrium pressure over the liquid (Pr). These values are expressed in mm Hg. In the above equation, Kg expresses the absorption constant in the gaseous boundary layer, the unit of which is kmol·m−2·h−1·mmHg−1. For the product of Kg S, an approximation is assumed, which leads to a constant value during the absorption process at each stage under certain conditions. Given this assumption, it can be concluded that the rate of absorption is affected only by the value of the potential. The matter gets a bit more complicated if we also take into account the influence of CO2. As can be seen from the above reaction equations, the absorption of both components into the liquid causes an exoenergetic reaction, which results in the release of heat and ultimately an increase in the temperature of the solution. This adversely affects the equilibrium pressure of carbon dioxide, by which its absorption efficiency decreases. An increase in temperature significantly accelerates the rate of reaction and reduces the viscosity of the solution, causing favorable conditions for ammonia diffusion, thereby facilitating the absorption of CO2 according to the chemical reaction. To avoid this, it was assumed that the absorption process would take place at an elevated temperature not exceeding 343 K.
The basic factors that guided the design, construction, and operation of the SAB absorber were as follows:
(a)
Method of conducting the process, continuous or semi-continuous;
(b)
Mass exchange area;
(c)
Liquid retention time;
(d)
Determination of the phase that controls the speed of the process;
(e)
Exoenergetic process: heat removal, use of cooling;
(f)
Probability of corrosion of apparatus; and
(g)
Rate of gas and liquid flow entering the apparatus.
The implementation of the task was carried out according to the layout shown in Figure 6 and Figure 8.
The basic scheme of the apparatus was rearranged in the following manner. Gas after RG-CC (FIC 4080.4) was introduced into the countercurrent SAB absorber, with Pall fill and cellular fill at the top of the apparatus, in countercurrent to the brine solution introduced through two stubs: onto the cellular fill (FIC 4080.9) and onto the Pall fill (FIC 4080.13), with the flux directed onto the Pall fill connected to the return flux (FIC 4080.15). The waste gas was directed to the corresponding absorber.
The ammoniated brine was directed to the cooler (A1.4082.1), from which the stream was divided into two parts: the turnback (FIC 4080.15), which is part of the flux directed to the Pall fill, and the flux of ammoniated brine (FIC 5104.15), which is directed in the start-up and test stage to the absorber tank, and in the carbonation stage to the corresponding CL location. The parameters of the SAB were as follows: The column was installed at +9900, designed as an apparatus consisting of three segments and a gas discharge head. The diameter and height of the apparatus were D = 1000 mm and H = 8000 mm. In the lower segment, there was a brine tank with a volume of V = 0.78 m3, above which is a pre-absorption zone on a bed of Pall rings with a contact area of 240 m2. The filling in the lower segment had the following characteristics: Pall rings, 50 × 50 × 2, 316 L, V = 1.80 m3. In the middle segment, flow was designed through two layers of cellular fill. Absorption of about 65% of the ammonia contained in the gas after passing through the Pall rings was assumed on the cellular fill. The characteristics of the fill in the middle segment were as follows: process cell fill of 2 × 3 layers, 316 L, φ 600. The apparatus was designed so that it would also be possible to remove the cell fill and replace it with another bed of Pall rings.
In the upper segment of the column was a designed demister, comprising a separation cell fill of three layers, 316 L, φ 600. A tray was placed under the column, from which all possible leaks from the stubs of the column could be drained by gravity.
The implementation of the SAB operation was carried out in such a way that the gas was started first and the brine flow was second. Once the flows were established, the appropriate amount of cooled brine per turn was determined and the process parameters were observed.
The research objective of the absorption process was to determine the motion parameters of the SAB absorber, in particular the relationship of total alkalinity = f (V brine V1, V brine V2) (Figure 9 and Figure 10), and the parameters under production conditions.
The analysis of the data, presented in Figure 9, allowed us to determine the optimal volume ratio for SAB liquids in order to obtain the highest possible alkalinity. The green–red area is optimal. For the selected area in Figure 9, the volume of the recycle stream was determined. From Figure 10, it can be concluded that in the range of 124–135 mmol∙20 cm−3, it increased linearly. The volume of the recycle stream (except from the recycle stream temperature) may be a parameter for controlling the direct alkalinity of the additional ammoniated brine.
The carbamate ion concentration is responsible for the increase in efficiency of the carbonation process. As seen in the results of the data presented in Table 4, the concentration of carbamate ions increased significantly during the addition of ammonification, while the parameters of the brine were maintained. The higher alkalinity due to the higher concentration of ammonia caused a greater fixation of the CO2 contained in the reaction solution. Increasing the concentration of NH2COO ions caused an increase in supersaturation in the crystallization area of NaHCO3. This effect is very favorable for the efficiency of the soda process.
Process energy changes are best presented in the form of a Sankey diagram [20,21]. The energy effect of the carbonation process is presented in Figure 11.
In the experiment, after increasing the concentration of ammonia in the carbonation process, there was an increase in temperature at the exit from the column. The efficiency of the soda process increased slightly, due to the suction effect on NaHCO3. This also has a beneficial effect on the process in summer conditions, when it is difficult to cool the carbonation column.
From the data presented in Table 4, it can be seen that the composition of the ammoniated brine strictly depends on the composition of the absorbed gas, and in particular on the amount of NH3 and CO2 in the gas fed into the SAB. The fact that the amount of CO2 in the gas after RG-CC translates directly into the amount of carbamate ions in the solution after SAB should receive special attention.
For the carbonation process, ammoniated brine solutions containing the highest possible concentration of NH2COO ions are preferred, which will cause a direct increase in supersaturation in CL.

5. Environmental Aspects

Minimized cooling water consumption resulting from the reduced cooling by 152.4 MJ∙Mg−1 was the main effect of the change in process regime. An additional ecological effect is presented in Figure 12. The concentrations of CO2 and NH3 in the exhaust gas were measured by infrared absorption using an LDS-6 laser process gas analyzer (Siemens, Germany). The diagram shows that during dosing of additional ammoniated brine, the CO2 and NH3 concentrations in the exhaust gases decreased. Additional ammoniated brine dosing was alternated with periods of no dosing. The CO2 concentration exceeded 25% v/v without additional ammoniated brine and was on average 3.85% lower with ammoniated brine. Likewise, the concentration of NH3 was 0.63% lower on average. Lowering the concentration of CO2 and NH3 in the waste gas has a very large ecological benefit in terms of total NaHCO3 production.
As a result of extensive analytical and industrial research at CIECH Soda Polska S.A.’s Inowrocław soda plants using ammoniated brine, the soda yield of the process was increased using the developed SAB device. Acknowledging that the soda industry is troublesome in terms of waste [1,16,18], this reduction in the amount of substrate at higher yields also minimizes the amount of process by-product. An additional environmental aspect, due to the thermodynamics of the process, is the increase in temperature at the outlet of the carbonation column, with unchanged crystal quality. This results in minimized cooling water consumption, which minimizes the cost of producing NaHCO3.

6. Conclusions

The study shows that the rate of ammonia desorption increases exponentially with increasing temperature and inert volumetric flow rate, while it decreases with increasing carbonation of the system. The value of the kinetic constant of the ammonia desorption process at an inert volume flow rate of 8.0 dm3·h−1 varied in the range 2–24 × 103 mmol·dm−3·s−1, while at a volumetric flow rate of 12.0 dm3·h−1 the range was 4–27 × 103 mmol·dm−3·s−1. The above relationships provide the basis for determining the optimal range of process parameters for the initial carbonation of ammonia brine using Solvay technology. Considering the hydrodynamic data of the SAB absorber, after conducting a series of experiments, the results of which are presented above, it can be assumed that the model parameters of the system are as follows: Average flux of pre-carbonated brine directed to the Pall fill (V2) is about 8 m3∙h−1; average flux of pre-carbonated brine directed to the cell fill (V1) is about 12 m3∙h−1; average turning flux is about 10 m3∙h−1; and increase in total alkalinity is about 36 mmol∙20 cm−3. The increase in overall CO2 in ammoniated brine can be described by the relation: y = 0.86x − 29.2 where y is overall CO2 and x is general alkalinity; the average gas flow for the above streams is about 2800 Nm3∙h−1. The MAB operating temperature is 334–340 K, and the mass exchange flux of ammonia between gas phase and liquid is 61.2 kg NH3∙h−1∙m−3; thus, this is the value of the transfer coefficient of ammonia under fixed hydrodynamic conditions. The temperature of the bottom of the carbonation column may be 307 K without additional cooling. This saves energy due to the need for less cooling water.

Author Contributions

Conceptualization, M.C. and Ł.K.; methodology, M.C. and B.K.; software, Ł.K. and S.Ł; validation, B.K., D.Ż. and M.S.; formal analysis, M.C. and B.K.; data curation, K.S.; writing—original draft preparation, M.C.; writing—review and editing, U.K.; visualization, S.Ł.; supervision, U.K.; project administration, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was carried out within a development project by CIECH R&D in cooperation with Nicolaus Copernicus University in Toruń (subcontractor) on a world innovative technology for the carbonation of ammonia brine enabling an increase in the soda yield of the soda ash production process (POIR.01.02.00-0031/2016) co-financed within Operational Programme Smart Growth 2014–2020, Measure 1.2. Sectoral R&D programs. Special thanks go to the Bialystok University of Technology Faculty of Construction and Environmental Sciences for facilitating the calculations.

Data Availability Statement

Data sharing is not applicable to this article.

Acknowledgments

Special thanks go to the Bialystok University of Technology Faculty of Construction and Environmental Sciences for facilitating the calculations.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

PLMlime mixer (pre-limer)
RH-DSdistiller heater
DSdistiller
RH-CDdiaphragm filter liquor heater
RG-RHammonia regeneration node gas cooler
DTdistiller liquid expander (detander)
CCammonia distillation column from condensates
MABsmall absorber (enrichment of ammonia brine into ammonia)
LCLcarbonation gas scrubber
ABM Iabsorber upper part
ABM IIabsorber lower part
SAammonia brine
SB–CL (SB/CL)scrubber of column gases
LPV–FLRfilter gas scrubber to vacuum pump
LPV-ABMabsorber gas scrubber to vacuum pump
RS–ABbrine tank after ammonia absorption
CLsettling (production) carbonation column
SAB small absorber (solution being implemented)
Vflow rate of the gaseous phase
Koverall mass exchange coefficient
Kv = K∙∆ is the size of phase separation area per unit volume
F0cross-sectional area of column
Xup, Xdownammonia content of input and output liquid from SAB
Yup, Ydownammonia content of input and output gas from SAB
Lup, Ldownvolume of liquid flow per unit time, input and output from SAB
Gup, Gdownvolume of gas flow per unit time, input and output from SAB
Gamount of substance passing from one phase to another
Fsize of phase separation area
ydriving force of mass exchange process, i.e., difference between actual concentration of a substance in a given phase and its achievable concentration at equilibrium
yactual concentration of substance in a specific phase
y*concentration of substance in same phase at equilibrium
y1absorbate concentration in gas phase before absorption
y2absorbate concentration in gas phase after absorption
ysrmean driving force
Hwcolumn fill height
Lliquid flow rate (kg∙s−1)
k x mass penetration coefficient for liquid film (kg∙m−2∙s−1)
aarea of filling unit (m2∙m−3)
Scross-sectional area of column (m2)
qinspecific flow of inertia
K N H 3 d e s   ammonia desorption rate constant (mmol·dm−3·s−1)
Rdegree of carbonation of solution
C N H 3 c O initial concentration of ammonia in working liquid
ctNH3concentration of ammonia in liquid after time t
Pgpartial pressure of absorbate
Prequilibrium pressure of adsorbate over liquid
Kgabsorption constant in gaseous boundary layer (kmol·m−2·h−1·mmHg−1)
Vmaxmaximum flow per unit time
CNaClsodium chloride concentration
Ttemperature
ttime
a[NH3] total ammonia concentration, compliant with ISO 21877
c[Cl] chloride concentration, compliant with ISO 5314:1981
q[Na+] sodium concentration, compliant with BN-90 9567-29
m[NH3(aq)] free ammonia concentration physically dissolved in solution
n[alkalinity] total concentration of alkaline compounds, compliant with ISO 9963-1
o[NH4+] modification of ISO 5314:1981
d[CO2] total concentration of carbon dioxide in solution, compliant with Standard Method 4500-CO2 Carbon Dioxide (2017)
w[CO32−] + [HCO3] + [CO2(r)] total concentration of carbonates and free carbon dioxide
e[NH2COO] carbamate ion concentration, NMR and precipitation methods
f[CO32−] carbonate ion concentration, compliant with EN ISO 10693:2014
g[HCO3] hydrocarbonate concentration, compliant with EN ISO 10693:2014
Iionic force
WNasodium yield of process

References

  1. Cichosz, M.; Igliński, B.; Buczkowski, R.; Rzymyszkiewicz, P. Rozwój i stan przemysłu sodowego na świecie. Przemysł Chem. 2018, 3, 319–325. [Google Scholar] [CrossRef]
  2. Kang, Z.; Jia, X.; Zhang, Y.; Kang, X.; Ge, M.; Liu, D.; Wang, C.; He, Z. A Review on Application of Biochar in the Removal of Pharmaceutical Pollutants through Adsorption and Persulfate-Based AOPs. Sustainability 2022, 14, 10128. [Google Scholar] [CrossRef]
  3. Zahmatkesh, S.; Amesho, K.T.T.; Sillanpaa, M.; Wang, C. Integration of renewable energy in wastewater treatment during COVID-19 pandemic: Challenges, opportunities, and progressive research trends. Clean. Chem. Eng. 2022, 3, 100036. [Google Scholar] [CrossRef]
  4. Liu, Z. National carbon emissions from the industry process: Production of glass, soda ash, ammonia, calcium carbide and alumina. Appl. Energy 2016, 166, 239–244. [Google Scholar] [CrossRef] [Green Version]
  5. Li, C.; Liang, Y.; Jiang, L.; Zhang, C.; Wang, Q. Characteristics of ammonia-soda residue and its reuse in magnesium oxychloride cement pastes. Constr. Build. Mater. 2021, 300, 123981. [Google Scholar] [CrossRef]
  6. Wang, X.-b.; Yan, X.; Li, X.-y. Environmental risk for application of ammonia-soda white mud in soils in China. J. Integr. Agric. 2020, 19, 601–611. [Google Scholar] [CrossRef]
  7. Cormos, A.-M.; Cormos, C.-C.; Agachi, P.Ş. Making soda ash manufacture more sustainable. A modeling study using ASPEN Plus. In Computer Aided Chemical Engineering; Pleşu, V., Agachi, P.Ş., Eds.; Elsevier: Amsterdam, The Netherlands, 2007; Volume 24, pp. 551–556. [Google Scholar]
  8. Zahmatkesh, S.; Klemeš, J.J.; Bokhari, A.; Rezakhani, Y.; Wang, C.; Sillanpaa, M.; Amesho, K.T.T.; Ahmed, W.S. Reducing chemical oxygen demand from low strength wastewater: A novel application of fuzzy logic based simulation in MATLAB. Comput. Chem. Eng. 2022, 166, 107944. [Google Scholar] [CrossRef]
  9. Hou, T.-P. Manufacture of Soda; Reinhold Publishing Corporation: New York, NY, USA, 1942. [Google Scholar]
  10. Zarzycki, R.C.A.; Coulson, J.M. Absorption Fundamentals and Applications; Pergamon Press: Oxford, UK; New York, NY, USA; Seoul, Korea; Tokyo, Japan, 1993; p. 647. [Google Scholar]
  11. Chu, F.; Yang, L.; Du, X.; Yang, Y. Mass transfer and energy consumption for CO2 absorption by ammonia solution in bubble column. Appl. Energy 2017, 190, 1068–1080. [Google Scholar] [CrossRef]
  12. Brian, P.L.T.; Hurley, J.F.; Hasseltine, E.H. Penetration theory for gas absorption accompanied by a second order chemical reaction. AIChE J. 1961, 7, 226–231. [Google Scholar] [CrossRef]
  13. Fedotieff, P.P. Der Ammoniaksodaprozess vom Standpunkte der Phasenlehre. Z. Phys. Chem. 1904, 49U, 162–188. [Google Scholar] [CrossRef]
  14. Qin, F.; Wang, S.; Hartono, A.; Svendsen, H.F.; Chen, C. Kinetics of CO2 absorption in aqueous ammonia solution. International J. Greenh. Gas Control 2010, 4, 729–738. [Google Scholar] [CrossRef]
  15. Mann, R.; Moyes, H. Exothermic gas absorption with chemical reaction. AIChE J. 1977, 23, 17–23. [Google Scholar] [CrossRef]
  16. Fortescue, G.E.; Pearson, J.R.A. On gas absorption into a turbulent liquid. Chem. Eng. Sci. 1967, 22, 1163–1176. [Google Scholar] [CrossRef]
  17. Cichosz, M.; Kiełkowska, U.; Skowron, K.; Kiedzik, Ł.; Łazarski, S.; Szkudlarek, M.; Kowalska, B.; Żurawski, D. Changes in Synthetic Soda Ash Production and Its Consequences for the Environment. Materials 2022, 15, 4828. [Google Scholar] [CrossRef] [PubMed]
  18. Onda, K.; Sada, E.; Kobayashi, T.; Fujine, M. Gas absorption accompanied by complex chemical reactions—III Parallel chemical reactions. Chem. Eng. Sci. 1970, 25, 1023–1031. [Google Scholar] [CrossRef]
  19. Skowron, K.; Sobczak, W.; Kiedzik, L.; Zelazny, R.; Buczkowski, R.; Wachowiak, M.; Cichosz, M. Method for Brine Carbonization with Additional Amount of Ammonia in a Process of Soda Production and System of Devices for Carrying Out This Method. European Patent EP3312143A1, 19 October 2017. [Google Scholar]
  20. de-Córdoba, G.F.; Molinari, B. Sankey diagrams for macroeconomics: A teaching complement bridging undergraduate and graduate Macro. Heliyon 2022, 8, e10717. [Google Scholar] [CrossRef] [PubMed]
  21. Glover, R.E.; Al-Haboubi, M.; Petticrew, M.P.; Eastmure, E.; Peacock, S.J.; Mays, N. Sankey diagrams can clarify ‘evidence attrition’: A systematic review and meta-analysis of the effectiveness of rapid diagnostic tests for antimicrobial resistance. J. Clin. Epidemiol. 2022, 144, 173–184. [Google Scholar] [CrossRef]
Figure 1. Diagram of absorber and process flux.
Figure 1. Diagram of absorber and process flux.
Energies 15 08370 g001
Figure 2. Amount of mass transfer for research absorption system.
Figure 2. Amount of mass transfer for research absorption system.
Energies 15 08370 g002
Figure 3. Solubility of NaCl in ammonia solution depending on degree of carbonation.
Figure 3. Solubility of NaCl in ammonia solution depending on degree of carbonation.
Energies 15 08370 g003
Figure 4. Solubility of NaCl decreases with increasing NH3 in solution, and simultaneously increases as degree of carbonation of solution rises.
Figure 4. Solubility of NaCl decreases with increasing NH3 in solution, and simultaneously increases as degree of carbonation of solution rises.
Energies 15 08370 g004
Figure 5. Traditional Solvay process absorption system implemented at CIECH Soda Polska S.A.’s Inowrocław soda plants (based on [19]).
Figure 5. Traditional Solvay process absorption system implemented at CIECH Soda Polska S.A.’s Inowrocław soda plants (based on [19]).
Energies 15 08370 g005
Figure 6. Modified Solvay process absorption system implemented at CIECH Soda Polska S.A.’s Inowrocław soda plants (based on [19]).
Figure 6. Modified Solvay process absorption system implemented at CIECH Soda Polska S.A.’s Inowrocław soda plants (based on [19]).
Energies 15 08370 g006
Figure 7. Ammonia desorption open system.
Figure 7. Ammonia desorption open system.
Energies 15 08370 g007
Figure 8. Coupled system to conduct the ammonization process (SAB).
Figure 8. Coupled system to conduct the ammonization process (SAB).
Energies 15 08370 g008
Figure 9. Total alkalinity dependent on V1 and V2 streams.
Figure 9. Total alkalinity dependent on V1 and V2 streams.
Energies 15 08370 g009
Figure 10. Total alkalinity dependent on Vrecycle stream.
Figure 10. Total alkalinity dependent on Vrecycle stream.
Energies 15 08370 g010
Figure 11. Energetic effects of soda in carbonation process, in MJ∙Mg−1.
Figure 11. Energetic effects of soda in carbonation process, in MJ∙Mg−1.
Energies 15 08370 g011
Figure 12. Concentrations of CO2 and NH3 in carbonation column exhaust gases.
Figure 12. Concentrations of CO2 and NH3 in carbonation column exhaust gases.
Energies 15 08370 g012
Table 1. Solubility of ammonia in water at 313 K at pressure of 1 atm.
Table 1. Solubility of ammonia in water at 313 K at pressure of 1 atm.
X = mol   NH 3 mol   H 2 O Y = mol   NH 3 mol   of   inert
0.26472.3628
0.21171.0822
0.15880.5605
0.10580.2816
0.07940.1865
0.05300.1112
0.04230.0869
0.03180.0629
0.02640.0520
0.02110.0410
0.01690.0327
0.01270.0246
0.01050.0207
Table 2. Effect of temperature, gas flow rate, and chemical composition of ammonia brine on ammonia desorption rate (experimental data).
Table 2. Effect of temperature, gas flow rate, and chemical composition of ammonia brine on ammonia desorption rate (experimental data).
Brine Chemical Composition of
coNH3
coNaCl
coCO2
Ammonia Concentration in ctNH3
Solution
Temperature T (K)
303313323
Gas flow Rate (qin dm3 h−1)
8.012.08.012.08.012.0
Desorption Time
(mol·dm−3)/(mval/20 cm3)t (s)
12345678
4.404/88.08
4.267/85.34
0.000/0.00
4.404/88.08
4.286/85.72
4.167/83.34
4.048/80.96
3.930/78.60
0
665
1210
1765
2320
0
370
740
1090
1490
0
380
765
1095
1460
0
190
385
605
800
0
175
350
515
685
0
100
200
305
415
4.394/87.88
4.247/84.94
0.953/38.12
4.394/87.88
4.276/85.52
4.157/83.14
4.038/80.76
3.920/78.40
0
1100
2190
3355
4565
0
680
1330
2035
2705
0
610
1165
1740
2355
0
525
1025
1640
2215
0
425
920
1330
1775
0
200
385
560
772
4.384/87.68
4.217/84.34
1.375/55.00
4.384/87.68
4.266/85.32
4.147/82.94
4.028/80.56
3.910/78.20
0
1270
2470
3640
4880
0
890
1770
2650
3675
0
940
1870
2790
3800
0
570
1125
1715
2270
0
502
989
1559
2101
0
360
658
1163
1541
Table 3. Effect of temperature, gas flow rate, and chemical composition of additionally ammoniated brine on ammonia desorption rate (experimental data).
Table 3. Effect of temperature, gas flow rate, and chemical composition of additionally ammoniated brine on ammonia desorption rate (experimental data).
Brine Chemical Composition of
coNH3
coNaCl
coCO2
Ammonia Concentration in ctNH3
Solution
Temperature T (K)
303313323
Gas Flow Rate (qin dm3 h−1)
8.012.08.012.08.012.0
Desorption Time
(mol·dm−3)/(mval/20 cm3)t (s)
12345678
6.258/125.2
4.148/82.96
0.000/0.00
6.258/125.2
6.140/122.8
6.021/120.4
5.902/118.0
5.784/115.7
0
460
980
1350
1720
0
275
550
780
1080
0
222
455
690
910
0
140
275
415
550
0
120
250
370
500
0
85
170
250
340
6.175/123.5
4.147/82.94
1.265/50.6
6.175/123.5
6.057/121.1
5.938/118.8
5.819/116.4
5.701/114.0
0
520
1030
1565
2095
0
345
685
1030
1385
0
370
715
1090
1460
0
245
495
730
980
0
205
420
615
815
0
145
295
440
590
6.141/122.8
4.209/84.18
1.966/78.64
6.141/122.8
6.041/120.8
5.903/118.1
5.785/115.7
5.666/113.3
0
665
1335
1990
2665
0
475
940
1340
1795
0
500
990
1485
1995
0
370
755
1105
1470
0
330
670
995
1315
0
240
490
725
965
Table 4. Results of analysis of and calculations for ammoniated brine.
Table 4. Results of analysis of and calculations for ammoniated brine.
Analyte
(mval∙20 cm−3)
Brine after SABBrine before SAB
12345
a ≡ [NH3]130.6132.8135.0129.1126.698.9
c ≡ [Cl]83.181.180.981.280.887.6
q ≡ [Na+]83.181.180.981.280.887.6
m ≡ [NH3(r)]54.455.381.9103.476.979.7
n ≡ [alkalinity]130.6132.8135.0129.1126.698.9
o ≡ [NH4+]49.550.336.119.334.312.9
pH9.619.629.9310.309.929.40
d ≡ [CO2]78.379.754.025.850.619.8
w ≡ [CO32−] + [HCO3] + [CO2(r)]25.125.220.113.419.77.3
e ≡ [NH2COO]53.254.533.912.430.912.5
f ≡ [CO32−]20.720.818.312.917.96.0
g ≡ [HCO3]4.44.41.90.561.91.3
k ≡ [acidity]78.3379.754.025.850.619.8
I (ionic force)197.2198.4168.9129.5165.995.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cichosz, M.; Kiełkowska, U.; Łazarski, S.; Kiedzik, Ł.; Szkudlarek, M.; Skowron, K.; Kowalska, B.; Żurawski, D. Influence of Ammonia Concentration on Solvay Soda Process Parameters and Associated Environmental and Energetic Effects. Energies 2022, 15, 8370. https://doi.org/10.3390/en15228370

AMA Style

Cichosz M, Kiełkowska U, Łazarski S, Kiedzik Ł, Szkudlarek M, Skowron K, Kowalska B, Żurawski D. Influence of Ammonia Concentration on Solvay Soda Process Parameters and Associated Environmental and Energetic Effects. Energies. 2022; 15(22):8370. https://doi.org/10.3390/en15228370

Chicago/Turabian Style

Cichosz, Marcin, Urszula Kiełkowska, Sławomir Łazarski, Łukasz Kiedzik, Marian Szkudlarek, Kazimierz Skowron, Beata Kowalska, and Damian Żurawski. 2022. "Influence of Ammonia Concentration on Solvay Soda Process Parameters and Associated Environmental and Energetic Effects" Energies 15, no. 22: 8370. https://doi.org/10.3390/en15228370

APA Style

Cichosz, M., Kiełkowska, U., Łazarski, S., Kiedzik, Ł., Szkudlarek, M., Skowron, K., Kowalska, B., & Żurawski, D. (2022). Influence of Ammonia Concentration on Solvay Soda Process Parameters and Associated Environmental and Energetic Effects. Energies, 15(22), 8370. https://doi.org/10.3390/en15228370

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