Next Article in Journal
Antimicrobial Potential of Tara Hydroalcoholic Extract (Caesalpinia spinosa) against Streptococcus Associated with Strep Throat
Previous Article in Journal
An Investigation on Optimized Performance of Voluteless Centrifugal Fans by a Class and Shape Transformation Function
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Hybrid Approach for Modeling and Optimisation of Phosphoric Acid Production through the Integration of AspenTech, SciLab Unit Operation, Artificial Neural Networks and Genetic Algorithm

1
Faculty of Technology Novi Sad, University of Novi Sad, Bulevar cara Lazara 1, 21000 Novi Sad, Serbia
2
Institute of General and Physical Chemistry, University of Belgrade, Studentski Trg 12-16, 11000 Belgrade, Serbia
3
Laboratory for Thermal Engineering and Energy, Institute of Nuclear Sciences “Vinča”, University of Belgrade, P.O. Box 522, 11001 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Processes 2023, 11(6), 1753; https://doi.org/10.3390/pr11061753
Submission received: 4 May 2023 / Revised: 29 May 2023 / Accepted: 6 June 2023 / Published: 8 June 2023

Abstract

:
The purpose of the study was to identify and predict the optimized parameters for phosphoric acid production. This involved modeling the crystal reactor, UCEGO filter (as a detailed model of the filter is not available in Aspen Plus or other simulation software), and acid separator using Sci-Lab to develop Cape-Open models. The simulation was conducted using Aspen Plus and involved analyzing 10 different phosphates with varying qualities and fractions of P2O5 and other minerals. After a successful simulation, a sensitivity analysis was conducted by varying parameters such as capacity, filter speed, vacuum, particle size, water temperature for washing the filtration cake, flow of recycled acid and strong acid from the separator below the filter, flow of slurry to reactor 1, temperature in reactors, and flow of H2SO4, resulting in nearly one million combinations. To create an algorithm for predicting process parameters and the maximal extent of recovering H3PO4 from slurry, ANN models were developed with a determination coefficient of 99%. Multi-objective optimization was then performed using a genetic algorithm to find the most suitable parameters that would lead to a higher reaction degree (96–97%) and quantity of separated H3PO4 and lower losses of gypsum. The results indicated that it is possible to predict the influence of process parameters on the quality of produced acid and minimize losses during production. The developed model was confirmed to be viable when compared to results found in the literature.

1. Introduction

The current global consumption of H3PO4 is on the rise, while the amount of basic raw material required for its production decreases each year. The world’s reserves of phosphate rocks are limited, with a total of 65,000 million tons [1]. The increased demand for phosphoric acid in the food, mineral waste, and pharmaceutical industries is driving market growth. It is estimated that the global phosphoric acid market will reach USD 47,718.13 million by 2028 [2]. Based on these trends, the depletion of phosphate rock reserves will continue. Therefore, it is necessary to take action and implement cost-saving measures to prolong the increase in consumption for as long as possible.
Although there are alternative ways to produce phosphoric acid [1], they are cost prohibitive, especially for developing countries. Consequently, the most efficient approach would be to increase the utilization of P2O5 from phosphate rock and minimize losses of acid in gypsum to slow down the declining trend of phosphate.
The morphological composition and quality of phosphate ore vary depending on its geographical origin. Even neighboring mines can have different compositions, which presents a challenge in maximizing the utilization of its main component, P2O5. Various production methods for phosphoric acid (such as Hemihydrate, Dihydrate, dry process, etc.) can be found in the literature [3,4,5,6,7,8,9,10,11,12,13]. However, the most commonly used process for the production of phosphoric acid is the Dihydrate process [8]. Different compositions of phosphate ore can lead to varying results in UCEGO filtration and the reaction in the reactor. The filter is sensitive to factors such as crystal size and shape, viscosity, and density, which can affect the filtration and reaction processes. Additionally, the presence of organic components can result in the formation of a large amount of foam, which can clog the filter mesh and hinder filtration, resulting in a greater loss of phosphoric acid in the gypsum cake. The separation processes of H3PO4 from gypsum using UCEGO filters have received little attention in the literature. While books [14,15] provide detailed descriptions of the filter used in phosphoric acid production, research papers typically focus on the reaction, without emphasizing the importance of the filtration [4,6,9,10,11,12]. It should be mentioned that population modeling and optimization of phosphoric acid production are presented in [13], where multi-objective optimization (MO) was conducted with high precision (relative difference 0.01–4.5%).
The results of the investigations discussed in this study enable the attainment of a high yield of H3PO4 in the reactor, taking into account various parameters and conditions in the filtration section (UCEGO filter and separator) as empirical assumptions and predictions. It is important to note that the processes are simulated statically, while the filter and reactor conditions, as well as the composition of each material flow, are dynamic. This study takes into consideration the fact that phosphate ore is present in many varieties of composition and that small variations in composition can lead to different process conditions necessary to obtain optimal results for the reaction and filtration. Therefore, the main aims of this study were to investigate the large number of different qualities of phosphate ore and to identify key process parameters that have a significant impact on phosphoric acid production on an industrial scale (both reactors and UCEGO filters). In addition, to establish connections between these parameters, as well as investigate their relationship with different qualities of phosphate ore, an algorithm was developed to help the broader scientific and engineering community predict the final outcome and quality of the produced acid for any type and quality of phosphate ore, as well as minimize losses of gypsum. This will enable the best process parameters to be found in order to increase the lifespan of phosphate ore and, consequently, lower the energy consumption of the process.

2. Materials and Methods

2.1. Process Simulation of Phosphoric Acid Production

To obtain the model and simulation, the software tools Matlab (version 2021b), SciLab (version 6.1.1), and Aspen Plus (version V12) were utilized. Each model was developed as a rigorous model, where all operating conditions and balances were taken into account to match the industrial conditions. The capacity of the simulation plant was 13,000 P2O5 per month or 156,000 tons of P2O5 per year.
To ensure the feasibility of this technology and obtain all balances, a complete simulation was performed. For this work, only the reaction and filtration sections were considered. In the production of phosphoric acid, it should be recognized that the process can be divided into two distinct parts. The first part involves the production of diluted phosphoric acid, which plays a pivotal role in determining the quality of the resulting acid and the losses encountered in the gypsum. The second part comprises the vacuum evaporation process, which aims to concentrate the acid by removing water according to market requirements and desired quality for mineral waste production. Significant consequences arise when the quality of the diluted acid falls below satisfactory limits during the initial stage of acid production. Such a scenario necessitates greater energy consumption during vacuum evaporation and may result in higher gypsum traces within the acid. These gypsum traces can subsequently contribute to the formation of deposits in heat exchangers, ultimately causing a reduction in the capacity of the concentration line. To address these critical factors, this research deliberately focuses on tracking the first part of production, encompassing the reaction and filtration sections. This choice allowed monitoring of the quality of the acid and its losses in gypsum, as these aspects significantly impact the overall process efficiency. Although there are expected differences between simulations and real production, this type of process simulator can be used to obtain reliable results for process operation, thanks to the available methods for modeling, the large component database, and the available thermodynamic packages. Initially, the components that constitute the phosphate composition were separated from the Aspen Plus database as the first step towards creating this simulation. Then, models for the reactor, filter, and separator were developed by combining SciLab unit operations with Aspen Plus. Finally, the operating conditions were chosen to imitate real-life production. To the best of the researchers’ knowledge, existing techniques often overlook the intricate interplay and coordination among the various sections of the UCEGO filter. They tend to focus on individual sections in isolation, without recognizing the potential cross-effects and interferences that can occur. This manuscript seeks to bridge this gap by analyzing and tracking changes in each part of the filter and establishing connections between them. The conventional approach of considering the filter as a single entity is not sufficient to address the challenges posed by inconsistent filtration quality across different sections.

2.1.1. Experimental Determination of Specific Resistance

To calculate and model the filter and obtain a reliable filter-specific resistance for the process simulation, a laboratory experiment was conducted. Through a series of carefully conducted experiments and referring to the literature [15], we established a correlation between the quantity of filtered acid and the time required for its filtration over the woven material. Determining the gradient and intercept of this correlation enabled us to derive the coefficients used in Equations (1) and (2), which facilitated calculation of the permeability and specific resistance, as demonstrated in the literature [15]. The woven filter used in the experiment was produced by Arena, Valjevo, Serbia.
To obtain a pulp sample, a measuring bottle with a volume of 1 L is used, and the sample is taken directly from the pipeline that supplies the mixture distributor. The bottle’s content is poured through a Bihner funnel, where a woven piece is placed, through which filtration takes place (manufacturer: Arena). The time necessary for complete liquid filtration from the pulp was measured using the same vacuum as in the process (Table 1).
The base equations used for calculating permeability and specific resistance are given in Equations (1) and (2), respectively, [15].
gradient = α av c μ l A f 2 Δ P f
intercept = μ l R A f Δ P f
These equations and experimental values are implemented in the model of the UCEGO filter, Negotin, Serbia. The composition and condition of the filter feed change frequently, and as such, it was necessary to create code that follows and implements these new parameters and their influences on the filtration of phosphoric acid from gypsum cake.

2.1.2. Obtaining Data for Simulation

Only a few articles published in the literature have focused on the composition of phosphate ore, its classification based on the bone phosphate of lime (BPL), and the variation in CaO, P2O5, and the mineral composition. Becker [16] addressed this gap and provided an opportunity to utilize the data from his book in simulations to predict and validate the data with real-life production. The primary challenge lies in the fact that varying compositions of a single component in the ore can result in different processes.
In this research, the different quality and composition of phosphate ores as selected from [16] are presented in Table 2.
Conversely, numerous articles provide information on and insight into the effects of specific impurities on the reaction and filtration processes [17,18,19,20,21]. Chaabouni’s research [21] outlines how the presence of Al2O3 and Fe2O3 affects filtration, as demonstrated by Equation (3).
A l 2 O 3 + F e 2 O 3 P 2 O 5 = f r o m   0.08   t o   0.1
Low values of Al2O3 aluminum promote the growth of gypsum and, consequently, the formation of uniform crystals, thus improving the filtration rate; however, high ratios of aluminum/fluorine increase the viscosity of the acid and its density, which has a detrimental effect on the filtration [21].
The influence of components SiO2, Al2O3, and F [21] on filterability can be predicted by Equation (4).
R C F = % F 19 % A l 2 O 3 17 + % S i O 2 15
The filterability of phosphoric muds increases with a decrease in the ratio of fluorine complexation (RCF) within the range of 0.42 ≤ RCF ≤ 1.378. This was noted for both raw phosphate 1 (RCF = 0.571) and raw phosphate 2 (RCF = 0.868), where the mass content of acid as P2O5 increases and corrosion decreases. However, if the ratio falls below 0.42, the filterability decreases due to the presence of undissolved silica salts and poor crystallization [21].

2.1.3. Description of the Chosen Filter and Separator Model

The UCEGO model is the most commonly used filter in phosphoric acid plants. Its characteristics are listed in Table 3 [22]. Depending on the plant’s capacity, filters of various sizes are produced. For this research, type UCEGO No. 9 was selected to model the essential data required to optimize the process (Table 3). Additionally, values such as filter area, number of segments, and diffusivity remain constant. Conversely, variables such as pressure, filter speed, number of segments allocated to each section (controlled by barriers and gondolas), process water flow, process water temperature, and the material flow of strong acid to storage tank vary.
Only one separator, labeled as position 13.01.02 (Figure 1), provides separation of the filtrate and suctioned air through a filter cake for each filtration sector. It consists of two borders, with the first border shorter than the second one, since the liquid level in that section is controlled by the pump speed. When the pump speed is reduced, the level rises and overflows in the second sector (recycled acid), which controls the density and fraction of the H3PO4 that returns to the reactor. Thus, the model (model pertains to the B23) was designed to simulate this process and has the ability to change and control the flow of strong acid and transfer excess to the return acid sector. This parameter is critical in phosphoric acid production, since the quality and quantity of recycled acid control the state in the reactor. If there is an insufficient amount of return acid, it can lead to reaction blockage and a rapid rise in density, which can make filtration impossible. Typically, it takes several hours up to a shift to regulate the process when this occurs.

2.1.4. Artificial Neural Network

In this study, Artificial neural network (ANN) models were developed to establish a correlation between input parameters and their impact on the production process. The use of an ANN model provides distinct advantages over traditional statistical and regression models. Its advanced capabilities lie in its adaptivity and ability to capture nonlinearity. Given the complexity of parameter changes in the production of phosphoric acid, tracking these changes during regular production becomes challenging. It is noteworthy that altering a single parameter may affect some but not all outlet parameters, while changing multiple input parameters simultaneously may influence other output parameters without affecting the previously modified parameter. To establish a model capable of learning and predicting these intricate relationships, it was deemed necessary to employ an ANN model. By leveraging the adaptive and nonlinear properties, we could effectively discover the connections between these changes. The percentage of samples used for training, validation, and testing was 70%, 15%, 15% respectively, the type of fitting neural network was fitnet, and the type of training algorithm was Levenberg–Marquardt.
The first model was designed to study the filtration process, which included all filtration sections, i.e., strong acid, recycled acid, weak acid, and gypsum disposal sections. This model aimed to identify the combined influence of various parameters on the outlet. The second model focused on monitoring the changes in reactor parameters and their impact on the outlet value. These models were developed separately to address the specific requirements and characteristics of each section. To achieve the best results, the number of neurons in the ANN model was varied from 2 to 30. The ANN model consists of input layers, hidden layers, and output layers.
By employing advanced artificial neural network techniques and employing a multi-objective optimization genetic algorithm, this manuscript aims to optimize the filtration performance across all sections simultaneously, thereby minimizing loss and improving the overall production efficiency. The novelty of this approach lies in its comprehensive nature and the integration of advanced artificial intelligence techniques to address the challenges associated with filtration in phosphoric acid production. To the best of authors’ knowledge, no previous studies have specifically focused on developing a model for the UCEGO filter that accounts for the individual sections and their coordination.

2.1.5. Multi-Objective Optimization

The procedure and steps involved in performing multi-objective optimization is presented in Figure 2. The procedure begins with modeling and simulation of phosphoric acid production, where the static change of parameters was achieved (step one). Step two involves choosing the range of input variables that are going to be used in Aspen Plus by performing sensitivity analysis in order to determine outputs. Using the parameters from the sensitivity analysis, we sought to find correlation between them by implementing Matlab. After creating the ANN model, it was possible to perform code generation in the MATLAB Matrix-Only function, which was later used in the genetic algorithm.
The obtained mathematical function is then used in the final step, the multi-objective analysis, employing the selected GA (genetic algorithm) module gamultiobj. The gamultiobj function in MATLAB utilizes a variation of the traditional GA known as the non-dominated sorting genetic algorithm II (NSGA-II). NSGA-II is an elitist algorithm that uses a combination of genetic operators such as selection, crossover, and mutation to evolve a population of candidate solutions. The result is a set of Pareto plots that displays the optimal parameters for phosphoric acid production. The path from modeling to simulation, sensitivity analysis, and ANN, up to MO, is used to solve the next tasks:
  • Limiting loses of phosphoric acid in gypsum;
  • Increasing filtration rate;
  • Achieving high quantity of strong phosphoric acid;
  • Increasing fraction of H3PO4 in strong acid.

Procedure

Several objectives were targeted in this multi-objective analysis. First, the reaction and filtration sections were observed independently. As the UCEGO filter consists of several zones, each zone was included, along with its effect on outlets, in the MO analysis.
As mentioned earlier, it is essential to consider the filtration and reaction section as a whole because an alteration in one affects the other. Hence, it was logical to proceed to the second step, which involves tracking the variations of reactor parameters and their impact on filtration.
The final steps involved monitoring the impact of parameter variations in both the reactor and the filter on the quality of the product and the filtration process. Thus, new MO analysis was conducted by taking into consideration the effect of all input changes on output values. This path is presented in Figure 3.

3. Results and Discussion

3.1. Model Validation

A static simulation of phosphoric acid production, which includes only the reaction and filtration sections, was performed using Aspen Plus. The flow-sheet diagram obtained from Aspen Plus is shown in Figure 4, and a list of the most important block explanations is presented in Table 4.
The phosphoric acid production section can be divided in two parts, the reaction and filtration sections. To perform sensitivity analysis, a separate stream from the SciLab model that represents the final solution or calculated value was created and combined with the Aspen model to maintain material balance. This was necessary, since Aspen cannot recognize the SciLab model and use it directly. For instance, the outlet stream from the UCEGO model, washing section 2 (VF2, UW2, VW2, Y2, F2, MLS2, MLD2, TW2, θP2), represents the calculated parameters, and a “dummy” inlet stream was added to the model to ensure mass balance.
An algorithmic scheme that provides a comprehensive visualization of the model’s architecture is included as Supplementary Materials to enhance readers’ understanding of the design. This illustrates the different components and their connections, providing a clearer representation of the model’s structure.

3.2. Confirmation and Validation of Simulation

The model developed for this simulation was compared to the data from the literature [16] to determine the accuracy. The phosphoric acid production process was simulated under the same conditions used in practice, including ore capacity, H2SO4 flow, return acid flow, temperature, return slurry flow, filter speed, water flow for washing the filter cake, vacuum on the filter, temperature of the process water, and other factors.
Table 5 presents a comparison between the simulation results and the literature data. The table indicates an acceptable agreement between the simulation and conditions in real production. Sensitivity analysis was used to identify the optimal solution and the necessary parameters to achieve the results that match those in the literature.

3.3. Selection of Input Variables for ANN Model

The production of phosphoric acid is a challenging process as it is influenced by numerous parameters, such as the capacity of phosphate, temperature, filter speed, the vacuum on the filter, recycled acid flow, H2SO4 flow, strong acid flow, return slurry flow, water for the filtration flow, particle size, steam for heating the water, and water temperature. These parameters were selected as input variables for the sensitivity analysis carried out in Aspen Plus, which generated the results. The value of input variables are given in Table 6 and Table 7.
All inputs were changed as specified in Table 6 and Table 7. Depending on the phosphate ore capacity and type of phosphate, H2SO4 flow increased by 1–3%, according to Becker [16]. For values such as filtration water flow and strong acid flow, parameters were chosen depending on phosphate ore capacity and the BPL of ore.

3.3.1. Effect of Pressure on Filtration

The pressure effect on filtration is illustrated in Figure 5. Figure 5a shows that as the vacuum pressure increases, the fractional separation of soluble components decreases. This occurs due to a decrease in specific resistance and permeability of the solid phase, leading to insufficient separation. Following the strong acid sector, the superficial velocity (u) in the return and weak acid sectors decreases, since it is directly proportional to pressure, causing a further decrease in liquid pore velocity (v). The decrease in these parameters results in a reduction of Reynolds and Smith numbers, Dl/D, Dn, and W, which in turn leads to a decrease in fractional efficiency (F). In this work, F1 represents the first washing section, and F2 the second washing section. The permeability of soluble components needs to be minimized. As shown in Figure 5, increasing the vacuum pressure results in a decrease in the fractional separation of soluble components due to the decrease in specific resistance and permeability of the solid phase. After the strong acid sector, in the return acid and weak acid sector, the superficial velocity is decreased due to its direct proportionality to pressure, leading to further decreases in liquid pore velocity. Applying a higher vacuum pressure will lead to the most soluble components, along with gypsum, ending up in the gypsum warehouse (Figure 5), but this also represents a loss of economy.

3.3.2. Effect of Filter Speed on Filtration

The influence of filter speed on the quantity of separated liquid, cake thickness, and throughput of the liquid and solute phase during the filtration process is shown in Figure 6 and Figure 7. With a decrease in the filter speed, higher separation of liquid Vf from pulp is achieved, reducing liquid and solid throughput, Ml and Ms, respectively.
Decreasing the liquid phase’s throughput leads to a higher separation of H3PO4 in strong acid (material flow M-5/1, Figure 4). This separation step is crucial for the second wash (material flow M-7/1, Figure 4) since all quantities left unused from the return acid and weak acid sector end up in gypsum (material flow M-10, Figure 4) and represent losses in production. Compensation for the concentration of H3PO4 in the return acid (material flow M-6/1, Figure 4) involves reducing the strong acid flow to the storage tank. As a result, some of the strong acid that is not directed to the storage tank overflows into the separator section of the return acid.

3.3.3. Process Water Temperature Effect on Return Acid and Weak Acid Sector

The diagram in Figure 8 illustrates the effect of the process temperature on the fractional separation of soluble components. As the temperature of the process water (material flow M-3, Figure 4) increases, the viscosity and density of the acid decrease. This results in an increase in surface acceleration, which is inversely proportional to viscosity and density. As a result, the superficial velocity increases, leading to an increase in liquid pore velocity. Consequently, the liquid pore velocity is directly proportional to the relationship of Re/Sc, which causes an increase in Dn and results in the change in W and F.
Higher temperature of the process water (material flow M-3, Figure 4) benefits the filtration process. It increases fraction recovery of soluble components, as the liquid quantity separates from the gypsum cake (Figure 8). As a result, the process loss is minimized by decreasing the throughput of soluble components.

3.3.4. Influence of Filter Feed Temperature

Figure 9 illustrates how the temperature of the feed flow (material flow Slurry, Figure 7) affects the cumulative volume of filtrate. Increasing the temperature of the feed flow, as described in the previous chapter, reduces the viscosity and density of acid, resulting in a higher cumulative volume of filtrate. The temperature was varied at four points: the lowest temperature of 65 °C, which represents the temperature achieved when the process starts after a technical day or equipment washing, and the highest temperature of 80 °C, which is maintained during the production process.

3.4. ANN Development

Complete Filtration and Reaction Section

The first step was to develop an ANN topology for the complete filtration section by changing the filtration section parameters and separately changing the reaction section parameters. After conducting sensitivity analysis in Aspen Tech simulation software, the Matlab toolbox was utilized. The goal of the evaluation was to find the mean square error and correlation coefficient (R) to determine the best performance of the ANN. As R approaches 1, the gap between predicted and actual data decreases, and as R approaches 0, it indicates a mismatch between the two sets of data. Figure 10 and Figure 11 illustrate the ANN topology based on changes in the filtration and reaction section parameters. It was necessary to manipulate the number of hidden neurons to find the best correlation.
For the filtration section, the inlet parameters of the ANN topology include temperature, vacuum, filter speed, particle size, steam flow, strong acid flow, and phosphate ore capacity. These parameters were iterated and used to determine the outlet parameter values, which include cake thickness, strong acid flow, MSOLID, MLIQUID, MSOLUB, H3PO4 in strong acid, gypsum in strong acid, % H3PO4 in strong acid, and SO4 in strong acid.
The recycled acid sector is located after the strong acid section and is considered one of the most crucial parts of the process, since its quality can impact the conditions in the reactor and affect the filtration process, resulting in lower quality of strong acid and higher losses of P2O5 in gypsum. The outlet monitored parameters include the extracted recycled acid, the concentration of H3PO4 in the extracted acid, the recycled acid sent to R1, the concentration of H3PO4 in the acid sent to R1, F1, Y1, MLD1, and MLS1.
The section that follows the recycled acid sector is the weak acid sector, and its significance in the filtration process cannot be overlooked. If the process quality in this section is satisfactory, it minimizes the loss of phosphoric acid in gypsum and improves the quality of the recycled acid and, consequently, the reactor’s condition. The outlet parameters that must be monitored include F2, Y2, MLD2, MLS2, and the temperature of the process water.
The last section under the vacuum is the gypsum section; the quality of the gypsum is an indicator of how the reaction and filtration processes have performed. If a large amount of phosphoric acid remains in the gypsum, it indicates that something went wrong during the reaction and filtration processes, and prompt action is required to reduce losses. The monitored outlet values are the H3PO4 loss in gypsum (ton/h) and the fraction of H3PO4 removed with gypsum from the filter.
The degree of correlation achieved for the filtration and reaction sections is 0.94655 and 0.97667, respectively, using 29 hidden neurons, as shown in Figure 12.
Disconnection on the parameter line is caused by the parameters between (210–240 ton/h) not being taken in consideration. These parameters represents recycled acid flow. As we used phosphate ore capacity between 70 ton/h and 90 ton/h, the recycled acid flow sent to the reactor must consider roughly three times this capacity, or at least 240 ton/h. The maximum capacity of recycled acid flow sent to the reactor is constrained by the pump capacity, which is 300 ton/h, and by the fact that the more recycled acid that is sent to the reactor, the lower the average resident time in reactor and, consequently, the lower the conversion degree. In a phosphoric acid plant, the ideal resident time is 3 h to 9 h, depending on the volume of the reactor and the capacity of the plant. As it is possible to produce more recycled acid by lowering the output of strong acid to the warehouse and increasing the quantity of water for washing the filtration cake, with pump construction constraints, this would lead to increased levels of fluid in the separator, causing the three sections to be mixed. This mixing will decrease the concentration of acid sent back to the reactor, lowering the density of pulp and leading to a low conversion degree of P2O5 from phosphate ore. After acid is mixed in the separator, filtration on first and second washing sections is disturbed, leading to poor extraction of H3PO4 from pulp and high losses of acid in the gypsum. Figure 12 shows correlation between the results of the sensitivity analysis and the parameters predicted with ANN.

3.5. MO (Multi-Objective Analysis)

3.5.1. MO (Multi-Objective Analysis Using GA) Strong Acid Sector

In Figure 13a, the correlation between the extracted H3PO4 and gypsum in the strong acid material flow (F-5-1, Figure 4) is shown. The main objective of the process is to achieve the highest quantity and fraction of extracted H3PO4 in the material flow (F-5-1, Figure 4). To achieve this goal, it is necessary to increase MSOLID and decrease MLIQUID. However, finding the best compromise of many parameters, such as filtration speed, cake thickness, density of infeed material flow, quantity of infeed, concentration of solids in infeed, viscosity, and vacuum pressure on the filter, is necessary. All these parameters are interdependent and will change during the filtration process, highlighting the complexity of achieving the desired results during the production process.
As shown in Figure 13, increasing the fraction of H3PO4 in the strong acid leads to an increase in the losses of acid in gypsum. To achieve a high mass fraction of phosphoric acid in the strong acid, it is necessary to increase the density and mass fraction of H3PO4 that is returned to the reactor as recycled acid. This requires reducing the quantity of strong acid that is sent back to the warehouse, so that all excess acid that is not sent can move to the recycled acid sector, increasing its density and H3PO4 fraction. However, when the recycled acid with these characteristics is sent back to the reactor, it can cause blockage of the reaction, requiring more H3PO4 and leading to the creation of smaller gypsum crystals, which impairs filtration. This also increases the quantity of H3PO4 in the slurry from the reactor to the filter. Therefore, a faster filter rotation speed is required, which reduces the time for filtration in each section, ultimately leading to an increase in the quantity of H3PO4 in gypsum. However, this may cause the filter to reach its capacity and fail to extract the required quantity of H3PO4 from the slurry before it reaches the final section (gypsum section).
In Table 8, case scenarios are presented regarding the relation between H3PO4 extracted from the filter and losses of phosphoric acid in gypsum, achieved by using the optimtool in Matlab and the gamultiobj (multiobjective optimization using genetic algorithm) solver. It is clear how important the difference between the extracted acid and the acid sent to warehouse is in the quantity of H3PO4 lost in gypsum. When there is a significant difference, all acid not used will be held in the separator and poured from the strong acid sector to the recycled acid sector. If the mass fraction of P2O5 (main component of H3PO4) increases to above 19%, it will lead to larger losses of gypsum, since not all P2O5 will manage to react, nor will it be possible to extract it from the filter (considering the construction limitations of the filter). In order to lower the speed filtration and to allow time for each section to perform at an optimal level, the thickness of the gypsum cake must be 4–5 cm. Lowering the rotation speed of UCEGO can lead to an increase in cake thickness, resulting in even worse filtration conditions; this is mostly influenced by the capacity of the phosphate ore. Lowering the capacity of the plant can lead to better filtration; however, economic considerations often lead to a demand for increased production.
In Table 9, some case scenarios for the reaction section are presented regarding mass fraction of H3PO4 in strong acid and losses of phosphoric acid in gypsum, as determined using the optimtool in Matlab and the gamultiobj (multiobjective optimization using genetic algorithm) solver. Besides the maintenance of temperature, slight changes in the flow of recycled acid, H2SO4, and water for filtration, as well as the difference between extracted and used strong acid, can significantly impact the quality of the produced acid and its losses.
First, the flow of recycled acid in the reactor has a significant impact on the reaction and, subsequently, the filtration condition. If the recycled acid has a sufficient concentration of P2O5 and flow is between 2.7–3 times the capacity of ore, reaction and filtration conditions will be of such that the filtration is optimized. If there is insufficient recycled acid to be sent to the reactor, the density of the pulp will dramatically increase, leading to worsening of filtration conditions and consequently leading to an even smaller extraction of acid from the recycled acid section, finally causing reaction blockage.
Second, the flow of H2SO4 can either lead to satisfactory reaction utilization or to its blocking. The excess SO42+ in strong acid must be 2.2–2.6%. If there is too little acid added to the reactor, then it will lead to poor reaction efficiency and excess of non-reacted P2O5 and CaO. On the other hand, an excess of H2SO4 added to the reactor will cause gypsum to surround P2O5 particles, preventing conversion to H3PO4 and leading to losses of this element in the gypsum.

3.5.2. Multi-Objective Analysis Using GA: Recycled Acid Sector

In Figure 14, a correlation is shown between liquid throughput and fractional solute recovery (Y1) and between solute throughput and fractional solute recovery in the recycled acid sector’s first washing sector. The main objective of the process is to achieve the lowest liquid and solute throughput and the highest fractional solute recovery, which will result in the highest possible amount of liquid and soluble phase extracted from the filter cake in the first washing sector. This will increase the concentration of the recycled material flow (F-6-1, Figure 4) and, consequently, the concentration of phosphoric acid in the reactor. Depending on the quality of the phosphate, different compositions and quantities of recycled acid need to be sent back to the reactor. The quantity of recycled acid necessary to be returned to reactor can be calculated in accordance with Becker [16]. All equations are presented in the Supplementary Materials.
In Table 10 and Table 11, some case scenarios for the reaction and filtration sections are presented regarding the relation between throughput of liquid and fractional solute, and Y1 and fractional solute recovery, determined using optimtool.
First, the influence of temperature is absolute; by increasing temperature, we achieve higher yield in the reactor, which consequently leads to better filtration conditions. The temperature should be over 80 °C as this will lead to a loss of one molecule of water in gypsum, which should be avoided. At above 82 °C, gypsum will lose both molecules of water.
Second, vacuum and filter speed are the main parameters on the UCEGO filter. It is mandatory to adjust these such that thickness of gypsum cake is between 4–5 cm. If speed is too low in comparation to phosphate ore capacity, the filter cake can become almost non-filtrable, which leads to high losses.
Third, particles can have different shapes and sizes, depending on the reaction conditions and recycled acid sent back to reactor. Small crystals such as the needle type will lead to blockage of pores on the woven filter, causing bad filtration and increasing losses of acid in gypsum. On the other hand, bigger plate-like crystal will benefit the process of filtration, increasing the mass fraction of H3PO4 in strong acid and leading to smaller losses in gypsum.
Forth, steam is used to preheat the water that is used for washing the gypsum cake. This water will first be sent to the second washing section, and after filtration, move to the first washing section; see Figure 8. For the filtration section, a lower bounds search space of [70 65 45 0.45 4 × 10−5 3 60] was established, while the corresponding upper bounds were [90 80 60 0.6 7 × 10−5 5 110]. These bounds ensure that the algorithm’s search is confined within a reasonable and meaningful range.
Similarly, for the reaction section, the lower bounds search space was set as [70 65 240 55 60 8 70], while the upper bounds were [90 80 270 78 100 9.5 100]. These bounds were carefully chosen to maintain the feasibility and practicality of the algorithm’s exploration in the reaction domain.

3.5.3. Multi-Objective Analysis Using GA: Weak Acid Sector

Figure 15a shows the correlation between liquid throughput–fractional solute recovery, and Figure 15b shows the Y2–fractional solute recovery in the weak acid sector and second washing sector. The main goal of process is to achieve the lowest throughput of liquid and highest fractional solute recovery. By achieving this, in the second washing sector, the lowest possible amount of liquid and soluble phase will be left in the filter cake as gypsum waste from the filter. If this is not achieved, a higher quantity of phosphoric acid will remain trapped in the filtration cake before the last section, which removes the gypsum from the filter. This represents losses in the production process; depending on capacity, such losses can be significantly large.
In Table 12 and Table 13, case scenarios for the reaction and filtration sections are presented regarding the relation between the throughput of liquid and fractional solute, and Y2 and fractional solute recovery, for the second washing section, achieved using the same optimization tools as in previous chapters.
The study demonstrates the feasibility of predicting the influence of process parameters on the quality of produced acid and minimizing losses during production. The viability of the developed model was confirmed through comparison with results reported in the literature.

4. Conclusions

The developed model and simulation have the potential to forecast process parameters and conditions in reactors and UCEGO filters. This presents an opportunity for plants to reduce energy consumption and, more importantly, enhance phosphoric acid extraction while decreasing acid losses in gypsum. By achieving this, the life span of phosphate ore could potentially be extended.
Working conditions of the UCEGO filter and reactor were determined using ANN and multi-objective analysis. Datasets were obtained using compositions of different phosphate ores from literature. The aim was to take in consideration different quality levels of phosphate ore, with different values of BPL, in order to create an achievable method that takes into consideration all relevant advantages and disadvantages.
Using the collected data and process parameters, approximately one million combinations were generated by modifying all feasible parameters during the production of phosphoric acid. This can assist engineers and operators in forecasting the optimal course of action and maintaining it throughout production. Subsequently, ANN models were established, which exhibited a high determination value of over 0.99.
Changing process parameters in this simulation proved to be adequate to decrease losses of H3PO4 in gypsum and increase its concentration in the produced acid. Constraints such as temperature in the reactor, mass flow of recycled acid from the filter to the reactor, H2SO4 mass flow, water for filtration, strong acid flow, vacuum on the filter, filter speed, and temperature of the filtration water were necessary in order to achieve optimal operation for the filter and reactor, a conversion degree of 96–97%, a fractional solute recovery over 0.96 for both washing sections, and the lowest liquid throughput.
The results achieved can be compared with those in the literature and in phosphoric acid plants, showing that simulation programs such as AspenTech, Matlab, SciLab, and SciLab Unit operation can be used to predict and optimize industrial production.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pr11061753/s1.

Author Contributions

Conceptualization, M.P. (Marko Pavlović) and P.K.; methodology, M.P. (Marko Pavlović); software, M.P. (Marko Pavlović); validation, M.P. (Marko Pavlović); formal analysis, M.P. (Marko Pavlović); investigation, M.P. (Marko Pavlović); resources, M.P. (Marko Pavlović); data curation, M.P. (Marko Pavlović), J.L., L.P., M.P. (Milada Pezo), O.B. and P.K.; writing—original draft preparation, M.P. (Marko Pavlović), J.L., L.P., M.P. (Milada Pezo), O.B. and P.K.; writing—review and editing, M.P. (Marko Pavlović), J.L., L.P., M.P. (Milada Pezo), O.B. and P.K.; visualization, M.P. (Marko Pavlović); supervision, J.L. and P.K.; project administration, J.L. and P.K.; funding acquisition, M.P. (Marko Pavlović), J.L., L.P., M.P. (Milada Pezo), O.B. and P.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to acknowledge the Ministry of Education, Science and Technological Development of the Republic of Serbia for their financial support, Project No. 451-03-47/2023-01/200134.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

AFiltration area, m2
AfFiltration area of woven filter in Buhner funnel, m2
ANNArtificial Neural Network
BMass fraction of gypsum per unit mass of slurry in reactor
BeckerConversion degree in a reactor based on calculation from Becker
cSolid concentration in feed, kg m−3
DMolecular diffusivity of solute, m2 s−1
DaSpecific gravity of phosphoric acid
DI1Impeller diameter for Reactor 1, m
DI2Impeller diameter for Reactor 2, m
DI3Impeller diameter for Reactor 2, m
DT1Reactor tank 1 diameter, m
DT2Reactor tank 2 diameter, m
DT3Reactor tank 2 diameter, m
dPVacuum pressure, kPa
F1Fraction of solute removed from a filter cake in washing section 1
F1Fraction of solute removed from a filter cake in washing section 2
GDegree of supersaturation in Reactor 1
G2Degree of supersaturation in Reactor 1
G3Degree of supersaturation in Reactor 3
GdiGypsum cake production, gypsum cake/ton rock
IImpurities apart from the P2O5 component in phosphate rock
KLPhysical mass transfer coefficient of H+ in solution (m/h)
LfThe thickness of the filter cake, m
MLD1Throughput of liquids in filter cake, ton h−1
MLD1Throughput of liquids in filter cake, ton h−1
MLS1Throughput of solute in filter cake, ton h−1
MLs2Throughput of solute in filter cake, ton h−1
MliquidThroughput of liquids in filter cake, ton h−1
MOMulti-objective
MsolidThroughput of solids in filter cake, ton h−1
MsolubThroughput of solute in filter cake, ton h−1
RResistance of fluid flow through a filter cloth, m−1
Re1Reynolds number in Reactor 1
Re2Reynolds number in Reactor 2
Re3Reynolds number in Reactor 3
SSupersaturation of calcium sulfate in Reactor 1
S2Supersaturation of calcium sulfate in Reactor 2
S3Supersaturation of calcium sulfate in Reactor 3
TTemperature in Reactor 1
T2Temperature in Reactor 2
T3Temperature in Reactor 2
TavgAverage or mean residence time in Reactor 1
Tavg2Average or mean residence time in Reactor 2
Tavg3Average or mean residence time in Reactor 3
TDDeliquoring time, s
TD1Deliquoring time for first section, s
TD2Deliquoring time for second section, s
TRComplete dissolution time of a single particle of radius R, h
TR1Complete dissolution time of a single particle of radius R in R1, h
TR2Complete dissolution time of a single particle of radius R in R2, h
TR3Complete dissolution time of a single particle of radius R in R3, h
Tw1Time of filtration for first washing section, s
Tw2Time of filtration for second washing section, s
Uw1Superficial velocity of a fluid first washing section, m3 m2 s−1
Uw2Superficial velocity of a fluid second washing section, m3 m2 s−1
VVolume of reactor 1, m3
V2Volume of reactor 2, m3
V3Volume of reactor 3, m3
VfFiltrate volume—strong acid sector, m3
Vf2Filtrate volume—second washing sector, m3
Vf1Filtrate volume—first washing sector, m3
VISKOZViscosity of reaction solution in Reactor 1, kg m h−1
Viskoz2Viscosity of reaction solution in Reactor 2, kg m h−1
Viskoz3Viscosity of reaction solution in Reactor 3, kg m h−1
Vw1Pore velocity of a fluid, first washing section, m s−1
Vw2Pore velocity of a fluid, second washing section, m s−1
VsiVolume of slurry produced in the reactor (slurry/ton rock)
X1Conversion degree in a Reactor 1
X2Conversion degree in a Reactor 2
X3Conversion degree in a Reactor 3
XavParticle diameter, μm
Greek Letter
αavSpecific resistance of a filter cake, m kg−1
ΔPfVacuum pressure, kPa
θP1Dewatering time—first Deliquoring section, s
θP2Dewatering time—second Deliquoring time, s
μlViscosity of the liquid in a feed, Pa*s
υeCrystal growth rate in Reactor 1, kg h−1 m2
υe2Crystal growth rate in Reactor 2, kg h−1 m2
υe3Crystal growth rate in Reactor 3, kg h−1 m2
υlLinear crystal growth rate in Reactor 1, m h−1
υl2Linear crystal growth rate in Reactor 2, m h−1
υl3Linear crystal growth rate in Reactor 3, m h−1
υm1Dissolution rate of phosphate rock per unit of particle surface in R1, kg h−1 m2
υm2Dissolution rate of phosphate rock per unit of particle surface in R2, kg h−1 m2
υm3Dissolution rate of phosphate rock per unit of particle surface in R3, kg h−1 m2
Φm1Mineral particle shape factor in Reactor 1
Φm2Mineral particle shape factor in Reactor 2
Φm3Mineral particle shape factor in Reactor 3
φd1Number of filtration segments dedicated to first Deliquoring sector
φd2Number of filtration segments dedicated to second Deliquoring sector
φfNumber of filtration segments dedicated to pre-sector
φf2Number of filtration segments dedicated to strong acid sector
φTTotal number of UCEGO filter segments
φw1Number of filtration segments dedicated to first washing sector
φw2Number of filtration segments dedicated to second washing sector
ωRotation speed of filter, rpm
ω1Reactor 1 impeller speed, rpm
ω2Reactor 2 impeller speed, rpm
ω3Reactor 3 impeller speed, rpm

References

  1. Salas, B.V.; Wiener, M.S.; Martinez, J.R.S. Phosphoric Acid Industry—Problems and Solutions; InTech: London, UK, 2017. [Google Scholar]
  2. Global Phosphoric Acid Market—Industry Trends and Forecast to 2028. Available online: https://www.databridgemarketresearch.com/reports/global-phosphoric-acid-market (accessed on 8 August 2022).
  3. Mathias, P.M.; Mendez, M. Simulation of phosphoric acid production by the dihydrate process. In Proceedings of the 22nd Clearwater Convention on Phosphate Fertilizer & Sulfuric Acid Technology Sheraton Sand Key Resort, Clearwater Beach, FL, USA, 22–23 May 1998. [Google Scholar]
  4. Kannamma, B.; Kannadasan, T.; Prabhakaran, D. Analysis and Simulation of Dihydrate Process for the Production of Phosphoric Acid (Reactor Section). Am. J. Eng. Res. 2013, 2, 1–8. [Google Scholar]
  5. Messnaoui, B.; Bouhaouss, A.; Derja, A. A Steady State Model for Describing the Phosphate Lattice Loss. Procedia Eng. 2012, 46, 134–142. [Google Scholar] [CrossRef] [Green Version]
  6. Mathias, P.M.; Chen, C.-C.; Walters, M. Modeling the Complex Chemical Reactions and Mass Transfer in a Phosphoric Acid Reactor. In Proceedings of the Third Joint China/USA Chemical Engineering Conference, Beijing, China, 25–28 September 2000. [Google Scholar]
  7. Papadopoulos, A.I.; Theodosiadis, K.; Seferlis, P. Modeling, Design and Optimisation of Industrial Phosphoric Acid Production Processes. In Proceedings of the 10th Conference on Process Integration, Modelling and Optimisation for Energy Saving and Pollution Reduction, Thessaloniki, Greece, 2007; Available online: https://folk.ntnu.no/skoge/prost/proceedings/icheap8-pres07/pres07webpapers/28%20Papadopoulos.pdf (accessed on 8 August 2022).
  8. Bouchkira, I.; Latifi, A.M.; Khamar, L.; Benjelloun, S. Modeling and Multi-Objective Optimization of the Digestion Tank of an Industrial Process for Manufacturing Phosphoric Acid by Wet Process. Comput. Chem. Eng. 2022, 156, 107536. [Google Scholar] [CrossRef]
  9. Abu-Eishah, S.I.; Abu-Jabal, N.M. Parametric Study on the Production of Phosphoric Acid by the Dihydrate Process. J. Chem. Eng. 2001, 81, 231–250. [Google Scholar] [CrossRef]
  10. Peng, Y.; Zhu, Z.; Braatz, R.D.; Myerson, A.S. Gypsum Crystallization during Phosphoric Acid Production: Modeling and Experiments Using the Mixed-Solvent-Electrolyte Thermodynamic Model. Ind. Eng. Chem. Res. 2015, 54, 7914–7924. [Google Scholar] [CrossRef]
  11. Yeo, Y.K.; Cho, Y.S.; Moon, B.K. Simulation of the Dihydrate Process for the Production of Phosphoric Acid. Proc. Korean Soc. Autom. Control. 1988, 1, 875–878. [Google Scholar] [CrossRef]
  12. Grema, A.S.; Iyodo, H.M.; Grema, A.S.; Imam, Y.Y.; Mohammed, H.I. Modeling and Simulation of Hemihydrate Phosphoric Acid Plant. Arid Zone J. Eng. Technol. Environ. 2018, 14, 169–171. [Google Scholar]
  13. Papadopoulos, A.I.; Seferlis, P. Generic Modelling, Design and Optimization of Industrial Phosphoric Acid Production Processes. Chem. Eng.Process. Process Intensif. 2009, 48, 493–506. [Google Scholar] [CrossRef]
  14. Tarleton, E.S.; Wakeman, R.J. Solid/Liquid Separation; Elsevier: Amsterdam, The Netherlands, 2007; ISBN 9781856174213. [Google Scholar]
  15. Richard, J.; Wakeman, E.S. Tarleton Filtration Equipment Selection, Modelling and Process Simulation, 1st ed.; Elsevier Advanced Technology: Amsterdam, The Netherlands, 1999. [Google Scholar]
  16. Becker, P. Phosphates and Phosphoric Acid: Raw Materials, Technology, and Economics of the Wet Process. Revised and Expanded; Marcel Dekker, Inc.: New York, NY, USA, 1989; Volume 6. [Google Scholar]
  17. Kybartiene, N.; Valancius, Z.; Leskeviciene, V.; Urbonas, L. Influence of the composition of phosphate rock on the amount of water-insoluble phosphate impurities in semi-hydrate phosphogypsum. Ceramics-Silikáty 2015, 29, 29–36. [Google Scholar]
  18. Pyagai, I.; Zubkova, O.; Babykin, R.; Toropchina, M.; Fediuk, R. Influence of Impurities on the Process of Obtaining Calcium Carbonate during the Processing of Phosphogypsum. Materials 2022, 15, 4335. [Google Scholar] [CrossRef] [PubMed]
  19. Fernandes, N.J.; Galvão, M.A.; Araujo, L.R.; Ataíde, C.H.; Barrozo, M.A.S. Effect of the Impurities on the Phosphoric Acid Process. In Materials Science Forum; Trans Tech Publications Ltd.: Stafa-Zurich, Switzerland, 2012; Volume 727, pp. 386–391. [Google Scholar]
  20. Zarenezhad, B. Effect of Impurities in Phosphoric Acid on the Granulometery of the Produced DAP(Di-Ammonium Phosphate) in Petrochemical Plants. Korean J. Chem. Eng. 2002, 19, 928–931. [Google Scholar] [CrossRef]
  21. Chaabouni, A.; Chtara, C.; Nzihou, C.A.; El Feki, H. Study the Nature and the Effects of the Impurities of Phosphate Rock in the Plants of Production of Phosphoric Acid. Adv. J. Chem. 2014, 7, 1296–1299. [Google Scholar] [CrossRef]
  22. Little, D.; Berberat, C. Ucego Filters in Phosphoric Acid Production; Peninsular Florida Section of the American Institute of Chemical Engineering: Tallahassee, FL, USA, 1981. [Google Scholar]
Figure 1. Separator.
Figure 1. Separator.
Processes 11 01753 g001
Figure 2. Procedure and steps involved in multi-objective optimization.
Figure 2. Procedure and steps involved in multi-objective optimization.
Processes 11 01753 g002
Figure 3. Combined ANN for filtration section.
Figure 3. Combined ANN for filtration section.
Processes 11 01753 g003
Figure 4. Flow sheet simulation of phosphoric acid production.
Figure 4. Flow sheet simulation of phosphoric acid production.
Processes 11 01753 g004
Figure 5. Influence of pressure on (a) fractional separation of soluble components and (b) permeability of liquid components.
Figure 5. Influence of pressure on (a) fractional separation of soluble components and (b) permeability of liquid components.
Processes 11 01753 g005
Figure 6. Influence of filtration speed on (a) amount of separated liquid and (b) liquid-phase throughput.
Figure 6. Influence of filtration speed on (a) amount of separated liquid and (b) liquid-phase throughput.
Processes 11 01753 g006
Figure 7. Influence of filtration speed on (a) amount of separated soluble and (b) cake thickness.
Figure 7. Influence of filtration speed on (a) amount of separated soluble and (b) cake thickness.
Processes 11 01753 g007
Figure 8. Influence of process water temperature on (a) fractional separation of soluble components and (b) amount of separated liquid from filtration cake.
Figure 8. Influence of process water temperature on (a) fractional separation of soluble components and (b) amount of separated liquid from filtration cake.
Processes 11 01753 g008
Figure 9. Influence of inlet temperature on amount of separated liquid from filtration cake.
Figure 9. Influence of inlet temperature on amount of separated liquid from filtration cake.
Processes 11 01753 g009
Figure 10. ANN model for all sectors—filtration section parameters.
Figure 10. ANN model for all sectors—filtration section parameters.
Processes 11 01753 g010
Figure 11. ANN for all sectors—reaction section parameters.
Figure 11. ANN for all sectors—reaction section parameters.
Processes 11 01753 g011
Figure 12. ANN results for all sectors: (a) filtration section parameters and (b) reaction section parameters.
Figure 12. ANN results for all sectors: (a) filtration section parameters and (b) reaction section parameters.
Processes 11 01753 g012
Figure 13. (a) MO correlation between extracted H3PO4 and gypsum in strong acid material flow achieved by changing filtration parameters; (b) MO correlation between fraction of H3PO4 in strong acid and losses of acid in gypsum.
Figure 13. (a) MO correlation between extracted H3PO4 and gypsum in strong acid material flow achieved by changing filtration parameters; (b) MO correlation between fraction of H3PO4 in strong acid and losses of acid in gypsum.
Processes 11 01753 g013
Figure 14. (a) MO correlation between throughput of liquid and fractional solute recovery achieved after changing filtration and reaction parameters, (b) MO correlation between Y1 and fractional solute recovery achieved after changing filtration and reaction parameters.
Figure 14. (a) MO correlation between throughput of liquid and fractional solute recovery achieved after changing filtration and reaction parameters, (b) MO correlation between Y1 and fractional solute recovery achieved after changing filtration and reaction parameters.
Processes 11 01753 g014
Figure 15. (a) MO correlation between throughput of liquid and fractional solute recovery achieved after changing filtration and reaction parameters; (b) MO correlation between Y2 and fractional solute recovery achieved after changing filtration and reaction parameters.
Figure 15. (a) MO correlation between throughput of liquid and fractional solute recovery achieved after changing filtration and reaction parameters; (b) MO correlation between Y2 and fractional solute recovery achieved after changing filtration and reaction parameters.
Processes 11 01753 g015
Table 1. Experimental values of pulp filterability.
Table 1. Experimental values of pulp filterability.
Time, s1440.852.963.284117.4151.3196
Volume, L0.050460.14390.185050.219560.285040.395750.50120.63709
Table 2. Phosphate compositions used in simulation [16].
Table 2. Phosphate compositions used in simulation [16].
Composition/
Phosphate
Morocco—Khourigba
65–66 BPL
Russia Kola
85 BPL
Morocco—Khourigba
70–71 BPL
Israel—Negev
72 BPL
H2O1.311.30.3
P2O530.5539.13233.2
SO31.8301.882
F3.72.893.93.7
SiO22.11.22.11
CO27.1264.9
Cl0.0200.020.04
CaO50.250.551.253
Al2O30.40.850.30.1
Fe2O30.20.450.20.05
MgO0.90.120.50.2
Na2O0.90.520.90.7
K2O0.070.130.070.01
C(org)0.2100.220
Composition/
Phosphate
Algeria—Dyobel Onk 65 BPLBrazil
76–78 BPL
Algeria—Dyobel Onk 75–77 BPLUS—North Carolina 59 BPL
H2O///1.3
P2O529.635.834.628
SO32.401.62.5
F3.71.5643.4
SiO22.41.462.214.4
CO274.881.24.1
Cl0.020.003/0.004
CaO48.552.954.244.1
Al2O30.40.350.50.4
Fe2O30.40.250.40.6
MgO1.40.790.80.5
Na2O1.40.260.60.9
K2O0.150.150.010.1
C(org)0.3001.6
Composition/
Phosphate
US—North Carolina
72 BPL
Finland
80–81 BPL
H2O0.28
P2O532.836.8
SO330
F3.92.6
SiO21.62.7
CO22.44.3
Cl0.0180.006
CaO52.551.6
Al2O30.40.3
Fe2O30.650.6
MgO0.631.07
Na2O10.27
K2O0.120.26
C(org)0.150
Table 3. UCEGO filter characteristics.
Table 3. UCEGO filter characteristics.
Filter Characteristic
A/m2φTφfφf2φw1φd1φw2φd2
109361810292
Working Conditions
ΔPf/kPaxav/μmD/m2 s−1ω/rpm
50201 × 10−90.6
Table 4. Model used in Aspen Plus.
Table 4. Model used in Aspen Plus.
Models/StreamSimulation IDModel Usage
Reactor 1Reactor 1SciLab model of reactor
Reactor 2Reactor 2SciLab model of reactor
Reactor 3Reactor 3SciLab model of reactor
UCEGO filterUCEGOUCEGO filter: for separation of phosphoric acid from gypsum
Pump13-01-03Strong acid pump: pump for sending collected strong acid in separator to warehouse
Pump13-01-04Return/recycled acid pump: pump for sending collected recycled acid from separator back to reactor
Pump13-01-05Weak acid pump: pump for sending collected weak acid from separator to recycled acid sector through filter for washing gypsum cake
Absorber 1B54For removing fluoride from gas phase produced in reactor
Absorber 2B55For removing fluoride from gas phase produced in reactor
PumpB58Pump for sending separated liquid phase from Absorber B54, which is used as washing media
Vacuum coolerB66Vacuum cooling is used to cool down part of produced mixture from the reactor, in order to control the temperature in the reactor
PumpB62Pump sends water for washing weak acid section through filter
SciLab modelB60Model control quantity of water sent to the pump B62 and, afterward, to the filter section
H2SO4 tankB24Storage acid tank
H2SO4 pumpB29, B31Pumps for sending H2SO4 to reactor. Depending on necessary capacity, they can work parallel or separately
SeparatorB23SciLab model of separator, which is used to store and separate different types of H3PO4 acid removed from UCEGO filter
PumpB27Pump sends weak acid as washing medium for recycled acid section through filter
Recycled acid correctionS242This stream is used in combination with the model of separator B23, as replacement for the frequency regulator of pump 13-01-04, to control the quantity of recycled acid that will be sent back to reactor
Strong acid correctionKOREKThis stream is used in combination with the model of separator B23, as replacement of the frequency regulator of pump 13-01-03, to control the quantity of strong acid that is sent to the storage tank. Excess of acid will overflow to the recycled acid sector in the separator
Table 5. Comparison between developed model and the literature.
Table 5. Comparison between developed model and the literature.
Morocco—Khouribga 65–66 BPLRussia Kolav 85 BPLMorocco-Khouribga 70–71 BPL
ComponentModelLiteratureComponentModelLiteratureComponentModelLiterature
P2O529.964/P2O530.128P2O529.98430.2
Izrael—Negev 72 BPLAlgeria—Dyobel Onk 65 BPLBrazil 76–79
ComponentModelLiteratureComponentModelLiteratureComponentModelLiterature
P2O529.15729P2O529.91529.1P2O528.45628
Algeria—Djebel Onk 75–77 BPLUS—North Carolina 59 BPLUS—North Carolina 72 BPL
ComponentModelLiteratureComponentModelLiteratureComponentModelLiterature
P2O527.65927P2O526.954/P2O528.63727.6
Finland 80–81 BPL
ComponentModelLiterature
P2O527.65927
Table 6. Input values for ANN model for phosphoric acid production–filtration section.
Table 6. Input values for ANN model for phosphoric acid production–filtration section.
InputsValueExtent/Increment
Phosphate ore capacity (ton/h)70–905
Temperature (°C)65–805
Filter speed (rpm)0.45–0.60.03
Vacuum on filter (mbar)450–60030
Strong acid flow (ton/h)60–110*
Particle size (µm)4–71
Steam flow (ton/h)3–50.5
* Values chosen depending on phosphate ore capacity and the BPL of ore.
Table 7. Input value for ANN model for phosphoric acid production–reaction section.
Table 7. Input value for ANN model for phosphoric acid production–reaction section.
InputsValueExtent/Increment
Phosphate ore capacity (ton/h)70–905
Temperature (°C)65–805
Recycled acid flow (ton/h)240–29010
H2SO4 flow (ton/h)55–78*
Strong acid flow (ton/h)60–100*
Return slurry flow (ton/h)8–9.50.3
Filtration water flow (ton/h)70–120*
* Values chosen depending on phosphate ore capacity and the BPL of ore.
Table 8. MO GA optimized parameters between produced strong acid and losses of phosphoric acid in gypsum.
Table 8. MO GA optimized parameters between produced strong acid and losses of phosphoric acid in gypsum.
Filtration Parameter
MinMax
Phosphate Ore Capacity (ton/h)74.1304874.7
Temperature (°C)78.31278.07
Vacuum (kPa)5151.21
Filter speed (rpm)0.540.41
Particle size (μm)5.69 ×10−55.4 × 10−5
Steam flow (ton/h)4.6874.682
Strong acid flow (ton/h)90.491.08
H3PO4 (ton/h)124.19863110.862
H3PO4 losses in gypsum (ton/h)4.50040.055
Table 9. MO GA optimized parameters between mass fraction of H3PO4 in produced strong acid and losses of phosphoric acid in gypsum.
Table 9. MO GA optimized parameters between mass fraction of H3PO4 in produced strong acid and losses of phosphoric acid in gypsum.
Reaction Parameters
MinMax
Phosphate Ore Capacity (ton/h)86.48988.554
Temperature (°C)73.976.26
Recycled acid flow (ton/h)227.14250.565
H2SO4 mass flow (ton/h)69.12865.5
Strong acid flow (ton/h)85.389.6
Slurry return (ton/h)8.0611.99
Water for filtration (ton/h)76.0382.65
H3PO4 mass fraction in strong acid0.4640.385
H3PO4 losses in gypsum (ton/h)3.0780.672
F10.9280.943
F20.9320.936
SO4 fraction in strong acid0.01650.026
Table 10. MO GA optimized parameters between throughput of liquid and fractional solute for first washing section, including best- and worst-case scenarios for filtration and reaction parameters.
Table 10. MO GA optimized parameters between throughput of liquid and fractional solute for first washing section, including best- and worst-case scenarios for filtration and reaction parameters.
Filter Parameters Reaction Parameters
MinMax MinMax
Phosphate ore capacity (ton/h)70.0099589.992949Phosphate ore capacity (ton/h)70.026789.76
Temperature (°C)79.04889570.058808Temperature (°C)79.62975.84
Vacuum (kPa)54.98439745.002289Recycled acid flow (ton/h)240.896254.99
Filter speed (rpm)0.4070620.402549H2SO4 mass flow (ton/h)60.174.54
Particle size (μm)5.59 × 10−54.05 × 10−5Strong acid flow (ton/h)70.02387.5
Steam flow (ton/h)4.5943344.0367736Slurry return (ton/h)11.9811.7
Strong acid flow (ton/h)94.97276261.679838Water for filtration (ton/h)83.3973
F10.96252540.8673016F10.9680.913
MLD1 (ton/h)36.57453673.55024MLD1 (ton/h)27.0160.111
Table 11. MO GA optimized parameters between Y1 and fractional solute recovery for first washing section, including best- and worst-case scenarios for filtration and reaction parameters.
Table 11. MO GA optimized parameters between Y1 and fractional solute recovery for first washing section, including best- and worst-case scenarios for filtration and reaction parameters.
Filter Parameters Reaction Parameters
MinMax MinMax
Phosphate ore capacity (ton/h)70.0085689.982Phosphate ore capacity (ton/h)7089.76
Temperature (°C)79.471.36Temperature (°C)79.00475.84
Vacuum (kPa)54.945.057Recycled acid flow (ton/h)246.93255
Filter speed (rpm)0.5480.402H2SO4 mass flow (ton/h)6074.54
Particle size (μm)5.6 × 10−54.02 × 10−5Strong acid flow (ton/h)70.287.5
Steam flow (ton/h)4.4684.521Slurry return (ton/h)1211.7
Strong acid flow (ton/h)72.965.1Water for filtration (ton/h)7073
F10.96760.889F10.95870.913
Y10.01460.1871Y10.01340.112
Table 12. MO GA optimized parameters between throughput of liquid and fractional solute for second washing section, including best- and worst-case scenarios for filtration and reaction parameters.
Table 12. MO GA optimized parameters between throughput of liquid and fractional solute for second washing section, including best- and worst-case scenarios for filtration and reaction parameters.
Filter Parameters Reaction Parameters
MinMax MinMax
Phosphate ore capacity (ton/h)70.0689.97166Phosphate ore capacity (ton/h)83.0188.6
Temperature (°C)78.406974.60292Temperature (°C)79.87876.9
Vacuum (kPa)54.92245.40905Recycled acid flow (ton/h)265.2265.3
Filter speed (rpm)0.43740.416919H2SO4 mass flow (ton/h)60.2373.7
Particle size (μm)5.85 × 10−54.05 × 10−5Strong acid flow (ton/h)60.2374
Steam flow (ton/h)4.9334.118Slurry return (ton/h)11.89.05
Strong acid flow (ton/h)70.274.257Water for filtration (ton/h)77.3684.46
F20.95870.86873F20.95880.9083
MLD2 (ton/h)13.229133.43067MLD2 (ton/h)12.17630.436
Table 13. MO GA optimized parameters between Y2 and fractional solute recovery for second washing section, including best- and worst-case scenarios for filtration and reaction parameters.
Table 13. MO GA optimized parameters between Y2 and fractional solute recovery for second washing section, including best- and worst-case scenarios for filtration and reaction parameters.
Filter Parameters Reaction Parameters
MinMax MinMax
Phosphate ore capacity (ton/h)70.82788.827Phosphate ore capacity (ton/h)83.7388
Temperature (°C)77.348870.58Temperature (°C)77.4670.44
Vacuum (kPa)53.908945.68Recycled acid flow (ton/h)276.1261.6
Filter speed (rpm)0.513080.4352H2SO4 mass flow (ton/h)60.10774.07
Particle size (μm)6 × 10−54.06 × 10−5Strong acid flow (ton/h)95.573.6
Steam flow (ton/h)4.7944.275Slurry return (ton/h)11.410.5
Strong acid flow (ton/h)88.460.5Water for filtration (ton/h)77.2779.05
F20.96670.8751F20.9580.90763
Y20.015120.202Y20.02640.1308
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pavlović, M.; Lubura, J.; Pezo, L.; Pezo, M.; Bera, O.; Kojić, P. A Novel Hybrid Approach for Modeling and Optimisation of Phosphoric Acid Production through the Integration of AspenTech, SciLab Unit Operation, Artificial Neural Networks and Genetic Algorithm. Processes 2023, 11, 1753. https://doi.org/10.3390/pr11061753

AMA Style

Pavlović M, Lubura J, Pezo L, Pezo M, Bera O, Kojić P. A Novel Hybrid Approach for Modeling and Optimisation of Phosphoric Acid Production through the Integration of AspenTech, SciLab Unit Operation, Artificial Neural Networks and Genetic Algorithm. Processes. 2023; 11(6):1753. https://doi.org/10.3390/pr11061753

Chicago/Turabian Style

Pavlović, Marko, Jelena Lubura, Lato Pezo, Milada Pezo, Oskar Bera, and Predrag Kojić. 2023. "A Novel Hybrid Approach for Modeling and Optimisation of Phosphoric Acid Production through the Integration of AspenTech, SciLab Unit Operation, Artificial Neural Networks and Genetic Algorithm" Processes 11, no. 6: 1753. https://doi.org/10.3390/pr11061753

APA Style

Pavlović, M., Lubura, J., Pezo, L., Pezo, M., Bera, O., & Kojić, P. (2023). A Novel Hybrid Approach for Modeling and Optimisation of Phosphoric Acid Production through the Integration of AspenTech, SciLab Unit Operation, Artificial Neural Networks and Genetic Algorithm. Processes, 11(6), 1753. https://doi.org/10.3390/pr11061753

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