1. Introduction
Sodium bicarbonate (NaHCO
3) is an inorganic, white, water-soluble salt that has a vast number of applications. It is used in the manufacture of cleaning products and plastic foams, in paper bleaching, as an additive in the food industry, as a powder for fire extinguishers, in animal feed, and in soft water treatment. In medicine, it is widely used because it is a weak base, and in cooking it is used as yeast thanks to the release of carbon dioxide [
1,
2]. This salt can be found naturally but is currently mostly produced industrially.
Sodium bicarbonate is a salt widely used in different industry sectors; this salt can be found naturally, but it is mostly obtained through industrial manufacturing. Because of its high demand in society, sodium bicarbonate’s market size was valued at USD 3387.3 million in 2023. The sodium bicarbonate industry is projected to grow from USD 3521.1 million in 2024 to USD 4800.4 million by 2032 [
3]. Two big names are responsible for how its production takes place; Nicolas Leblanc and Ernest Solvay were creators of the routes with the same name in the industry.
The origin of producing sodium bicarbonate is linked to that of sodium carbonate. The production of sodium carbonate involves evaporating extracts obtained from the ashes of various plants, including Spanish Barilla [
4], the coastal plant known as barrilheira [
5], and several marine plants such as Chenopodium, Salicornia, and Salsola, among others [
6]. Until the end of the 18th century, these plant sources were the primary providers of sodium carbonate [
7]. The French government, in 1775, seeking to avoid foreign dependence on the import of natural sodium carbonate, established a prize for those who proposed a satisfactory process for obtaining Na
2CO
3 from NaCl [
8,
9,
10]. The French chemist Nicholas Leblanc (1742–1806) competed with the process known as the Leblanc Route, which is based on the conversion of NaCl into Na
2SO
4 by treatment with H
2SO
4. A solid mixture of Na
2CO
3 and CaS is then formed by heating the Na
2SO
4 with coal and CaCO
3, and the soda ash is extracted with water [
9]. Leblanc patented the synthesis of soda ash, but after his death, his patent was confiscated, causing his process to be adopted in many countries and continuously improved [
10,
11].
Because of the limitations of the Leblanc process, the Solvay process was created by Belgian industrial chemist Ernest Solvay in 1860. This process was rapidly disseminated throughout Europe, and the Leblanc process was almost completely replaced by this new process [
10,
12]. The main reagents of the Solvay process are CO
2, NH
3, and NaCl, which are responsible for the production of carbonate compounds. Ammonia (NH
3) acts as an intermediary and is recycled at the end of the process, resulting in minimal losses in relation to the final product [
13].
From 1890 onward, most soda ash and sodium bicarbonate were produced by the Solvay process, a scenario that continues to this day [
14].
Other raw materials have been studied as alternatives to produce sodium bicarbonate; one of the most promising is sodium sulfate. The reaction of sodium sulfate with carbon dioxide and ammonia has been studied as a route to produce sodium bicarbonate, with the intention of replacing traditional routes [
15,
16].
In this route, sodium sulfate participates in a double-exchange reaction, providing sodium (Na
+) and sulfate (SO
42−) ions. Then, carbon dioxide is used to precipitate the sodium bicarbonate, forming reactive bicarbonate ions, and a positive ion is fed into the reaction medium to obtain the solubility potential [
17]. Thus, this route promotes mitigation of CO
2 emissions and produces ammonium sulfate as a by-product, which is a fertilizer of great commercial interest. However, baking soda is very susceptible to salt contamination by this route.
4. Results and Discussion
The simulation effectively demonstrated the generation of the sodium bicarbonate; the reaction occurred in the crystallizer, and the stream coming out of this model was separated in a centrifuge. The mass fractions of the “C-LIQ” stream, corresponding to the mother water stream, are presented in
Table 3.
With the results of the mother water stream, we started the separation studies.
Figure 5 illustrates the result of sensitivity analysis in a single flash.
After the results of the sensitivity analysis, we found that from these conditions, only two salts were precipitated including salt 5, ammonium sulfate, and salt 7, sodium sulfate; therefore, the other salts are not shown in
Figure 5. This first result is consistent with the literature, where Kresnyak, Halldorson, and Hantke [
16] observed that with a change in the temperature of the mother water, the formation of sodium sulfate and ammonium sulfate occurred.
From this analysis, we observed that there was no possibility of precipitating sodium bicarbonate. Thus, there was no concern about exceeding the decomposition temperatures of sodium bicarbonate found, taking into account only other factors in the choice of equipment operating conditions, such as the removal of gaseous components or water recovery.
According to
Figure 6, the removal of gases was hampered at vacuum pressures. At ambient pressure, we observed that there was a greater withdrawal of CO
2 working towards the temperature condition of 90 °C. At higher temperatures the condition of supersaturation occurred, causing the precipitation of salts. In this way, temperatures close to 95 °C should not be exceeded.
To check if there is a possibility of greater carbon dioxide withdrawal, a second flash operating at a pressure of 1 bar was inserted immediately after the first. The first flash operated at a temperature where there was recovery of gaseous components; however, there was no precipitation, so there was no problem in the operation of the second stage. Thus, the operating conditions for the first flash were used according to the first sensitivity analysis. The results of the analysis are shown in
Figure 6.
For temperatures below 55 °C, the data presented an error status in the program, possibly because separation by distillation did not occur in both flashes. Because the temperature of the second flash was lower than that of the first, the separation occurred only in the first flash. Thus, these data were removed before constructing the graphs represented in
Figure 6. At approximately 65 °C, the sodium sulfate began to precipitate, so the working temperature was set to 60 °C. Even without analyzing the cost of implementing this equipment, it is possible to see that it is not favorable to insert a second flash in sequence with the first because there was an increase of only 3.29% in the removal of carbon dioxide.
When comparing the flashes in recycling and in sequence, we observed that there was a minimal increase in CO
2 removal, and there was no sense in using a second flash. To verify the possibility of greater CO
2 withdrawal, the addition of a second flash was analyzed again; however, it worked in recycling with the first flash. The result of the sensitivity analysis is shown in
Figure 7.
When a second flash was added in recycling to the first, we observed that before the precipitation of any salt occurred, it was possible to remove 684.82 kg/h of CO
2, representing an increase of 3.28% compared with the removal of a single flash. To complement this study, we analyzed whether it would make sense to add a second flash from an economic perspective. Thus,
Table 4 compares the cost of the equipment, installation, quantity of utility, and its price, as well as the diameters and volumes of the equipment
According to
Table 4, the addition of a second flash in both cases is practically the same as the first one when considered in relation to its size, i.e., the cost of adding a second flash will not compensate for the amount of CO
2 removed. In addition, when comparing the burst and recycled flashes, there is a large decrease in utility usage in the second arrangement, but it is not worth it for the amount of CO
2 removed. Based on the sensitivity analyses for the flash distillation step, it can be observed that the increase in the number of flashes will cause an increase in CO
2 removal, but their addition is not economically viable. Thus, only one flash is necessary to remove the gas compounds.
For the evaporator, the study started with an analysis performed by varying the steam flow rate and the evaporator pressure. As with the flash, there was no formation of salts other than salt 5 and salt 7. According to the result, it was seen that for larger amounts of water withdrawn, the temperature of the evaporator must be very high, causing the formed salts to degrade. Thus, we decided to work in a vacuum, with a pressure of 0.1 bar, so that the equipment could work properly. This analysis was used to determine the pressure of the evaporator, and the resulting data are not shown here. A second sensitivity analysis was performed in order to verify the required steam flow. The results are shown in
Figure 8.
According to
Figure 8, it is possible to remove approximately 24% of water from the mother water, using approximately 1360 kg/h of steam. If the temperature of the evaporator increases, the precipitation of salt 7 will begin to occur, which is not expected in this equipment. From these results, it is possible to estimate the pressure drops for multi-effect evaporators.
Figure 9 and
Figure 10 show the result of the sensitivity analysis for the evaporator of two effects and three effects, respectively.
For a two-effect evaporator, it is clear that before the precipitation of any salt occurs, it is possible to remove approximately 19% of the water using 940 kg/h of saturated steam. For a three-effect evaporator, the required inlet steam flow rate is 835 kg/h, since after this flow, the formation of sodium sulfate occurs. In this condition, it is possible to remove around 17% of the water from the solution.
A comparison of the evaporators is shown in
Table 5An important answer obtained by the sensitivity analyses for the evaporators is that there are no underlying amounts of sodium bicarbonate to be precipitated; there is only the possibility of forming sodium sulfate and ammonium sulfate. Thus, these results show that the quantities of raw material were optimized adequately by Yoshi et al. [
21].
The sensibility analysis shows that sodium sulfate is in higher concentration in the solution and that it is not possible to precipitate ammonium sulfate without contamination by sodium sulfate. As such, precipitating sodium sulfate may not be advantageous, as it is an essential process reagent and would need to be dissolved again to be reused. Even so, a precipitation study was carried out to determine if it was possible to extract ammonium sulfate with high purity. This investigation was important to determine whether it was possible to withdraw large amounts of water, as well as to ensure that by-products could be managed efficiently and sustainably.
Before starting the precipitation itself, it was necessary to understand how the salts behave in the solution they were in. The solubility diagrams for the pressures of 0.1, 0.5, and 1 bar are presented in
Figure 11.
It can be seen that the greatest difference between the solubility curves occurs in the pressure curves of 0.1 bar and ambient pressures. For the 0.1 bar curves, the solubilities of sodium sulfate and ammonium sulfate do not reach 100 °C, but the solubility of ammonium sulfate reaches 57.5 °C and sodium sulfate reaches 50 °C. This is probably related to the thermal decomposition of these salts, so the software cannot compute the solubilities because the material is in the gas phase and no longer in the solid phase.
In relation to sodium sulfate, which was the first salt to be precipitated, it is observed that at all pressures, it presents lower solubility and that its curve only occurs above 30 °C. In addition, its curve is related to a dissolution phenomenon that occurs exothermically; that is, solubility decreases with increasing temperature, so it is preferable to work at higher temperatures to facilitate the precipitation of this salt.
It is seen that the solubility curve has a similar behavior between the pressures of 0.1 and 0.5 bar and possibly at 1 bar as well. However, the curve for 1 bar should exceed 100 °C to see the whole behavior.
In general, we found that below 25 °C, the solubility of Glauber’s salt is lower than that of ammonium sulfate, and above 30 °C, the solubility of sodium sulfate is lower than that of ammonium sulfate. Therefore, it is ideal to work between values close to 25 and 30 °C to avoid contamination from another salt when aiming to precipitate ammonium sulfate.
It is inferred that the ideal process would be the removal of the ammonium sulfate first. However, this analysis was performed, and the crystallization of any salt was not obtained; so, to remove the sodium sulfate, an evaporator–crystallizer was used, increasing the temperature, and to precipitate the sodium sulfate, a crystallizer was used, cooling the solution.
To find the conditions of the crystallizing evaporator, and the equipment responsible for precipitating sodium sulfate, a sensitivity analysis was performed. The results of this analysis are shown in
Figure 12.
As expected, there is no precipitation of other salts in addition to sodium sulfate and ammonium sulfate salts. As the purpose of this equipment is now crystallization in addition to the removal of water from the solution, the best operating condition is a steam flow at the inlet that supplies enough energy to crystallize the largest amount of sodium sulfate without precipitation of ammonium sulfate. According to
Figure 12, in approximately 6900 kg/h of steam, ammonium sulfate begins to precipitate, so a flow rate of 6890 kg/h was chosen. In this flow, there is a precipitation of 1278.89 kg/h of sodium sulfate. Unlike a common evaporator, where you should not exceed the saturation temperature of the solution, in the evaporator with crystallization, this temperature is exceeded. Consequently, the precipitation of salts will occur and more water can be evaporated. In this case, the vapor removal goes from approximately 24% of the common evaporator to close to 86%. The composition of the steam stream is shown in
Table 6.
Table 6 shows that the vapor stream contains three components that are raw materials for the process, thus indicating the possibility of reintroducing this stream back into the process.
The result of the sensitivity analysis for the crystallizer are shown in
Figure 13.
It can be seen that by lowering the temperature, there is a considerable increase in the crystallization of ammonium sulfate, without crystallization of sodium sulfate; however, below 24 °C, there is also crystallization of Glauber’s salt. As Glauber’s salt is not economically attractive, and to increase the purity of ammonium sulfate, the most suitable crystallizer temperature is 24 °C.
Table 7 and
Table 8 show the global composition of the “SULSO” and “SULAMO” currents, representing sodium sulfate and ammonium sulfate composition, respectively.
It is observed that sodium sulfate presents little contamination, having a 0.9434 mass fraction. To determine if the crystallization of sodium sulfate as a by-product of the route is attractive, it is necessary to carry out an economic study, considering the added equipment and the cost for the crystallization of sodium sulfate, followed by its dissolution for entry into the process.
An important finding of this study is that the use of evaporation equipment is not advantageous without considering the crystallization of other salts. The first point is that a limited amount of vapor is obtained, which means that even with process recycling, the impact would be minimal compared with the cost of adding this new equipment. Furthermore, attempts to increase evaporation effects to improve the process results in facilitating salt precipitation rather than increasing evaporation. In contrast, when evaporation occurs alongside salt precipitation, it is possible to use these salts as raw materials or by-products, in addition to allowing the removal of a large amount of vapor, even in a single effect. Thus, evaporation combined with crystallization presents a promising opportunity for optimizing this route.