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:
and
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
and total sodium chloride (
c) in solution (a
1–a
4, 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 CO
2 introduced into the solution and is limited by the amount of CaCO
3 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 CO
2 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]:
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, ≈ 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, ≥ 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 . It was assumed that 90% of ammonia would be absorbed; therefore, . 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:
The driving force of the mass exchange process can be expressed in terms of difference:
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:
From the material balance, it is known that:
From the comparison of the above quantities we obtain:
Considering
V from the material balance equation, we obtain:
Transforming the equation and integrating the limits from 0 to
F and from
y2 to
y1, we obtain:
Therefore, we can obtain:
where:
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:
Eventually, combining the above equations with the basic differential equation, we obtain:
Integrating the above relationship from 0 to
Hw and from
y2 to
y1, we obtain:
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
, where X
up ; hence, X
down = 0.2217. By drawing a graph of
, 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:
The number of mass transfer units is 0.2217.
We assume that the height of mass penetration for the liquid phase is:
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 NH
3 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 NH
3·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: , , and .
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 NH
3 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:
Within the studied range of temperature changes, ammonia brine chemistry, and inert volume flow rate, the desorption intensity of NH
3 was approximated by the equation:
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 CO
2 and the co-crystallization of NH
4HCO
3 or NaCl. The absorption of ammonia in aqueous solution can be represented by the following equation:
The absorption of carbon dioxide in an aqueous ammonia solution can be presented as follows:
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:
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 CO
2 contained in the reaction solution. Increasing the concentration of NH
2COO
− ions caused an increase in supersaturation in the crystallization area of NaHCO
3. 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 NH
3 and CO
2 in the gas fed into the SAB. The fact that the amount of CO
2 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 CO
2 and NH
3 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 CO
2 and NH
3 concentrations in the exhaust gases decreased. Additional ammoniated brine dosing was alternated with periods of no dosing. The CO
2 concentration exceeded 25%
v/
v without additional ammoniated brine and was on average 3.85% lower with ammoniated brine. Likewise, the concentration of NH
3 was 0.63% lower on average. Lowering the concentration of CO
2 and NH
3 in the waste gas has a very large ecological benefit in terms of total NaHCO
3 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 NaHCO
3.
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.