1. Introduction
Chitosan is a linear N-acetyl polysaccharide derived from chitin. It is the second most abundant biopolymer in nature and can be produced from the exoskeletons of crustaceans such as shrimp, prawns, and lobsters, among others [
1]. Chitosan finds extensive applications in agriculture [
2,
3], food [
4,
5], cosmetics, and pharmaceuticals [
6], due to its hydrophilic nature [
7], biocompatibility [
6], biodegradability [
8], antibacterial properties [
9], remarkable affinity for biomacromolecules [
10], and non-toxicity [
11]. Moreover, it can be chemically modified to enhance its physical and chemical properties [
12]. With its amino functional groups (-NH
2), chitosan has been widely used as an adsorbent for water treatment, removing dyes, organic pollutants, and heavy metals [
13].
The increasing pollution of water bodies has driven the development of bioadsorbents with an enhanced mechanical strength, adsorption capacity, selectivity, and easy recovery. One of the modifying agents used to improve the recovery of chitosan-based bioadsorbents from treated waters is magnetite nanoparticles (Fe
3O
4), owing to their mixed-valence ferromagnetic nature with divalent and trivalent cations [
14]. When particles are scaled down to the nanometer range, they become superparamagnetic, preventing self-aggregation [
15]. The modification of biopolymers, such as chitosan, lignin, cellulose, and poly-and (γ-glutamic acid), among others, with magnetite nanoparticles has been studied for the adsorption of contaminants such as methylene blue, eriochrome black T, and heavy metals. Specifically, magnetically modified chitosan offers advantages such as a strong metal chelation capacity due to the presence of amino and hydroxyl groups in the chitosan chain [
16].
Likewise, the modification of biopolymers such as chitosan and alginate with thiourea improves the adsorption capacity of the biopolymers by increasing the number of active sites. This is because the thiol group forms stable complexes with metals such as Hg, Pt, Ag, Au, and to a lesser extent, Cd and Zn [
17]. Thiourea-crosslinked chitosan microbeads exhibited an enhanced adsorption capacity for platinum due to their increased number of amino groups compared to unmodified chitosan, going from 80.6 mg/g to 129.9 mg/g using modified and unmodified chitosan, respectively [
17]. Additionally, chitosan–magnetite nanocomposites have been successfully applied in the adsorption of toxic metals such as Pb(II) and Ni(II), and the beads were easily recovered through the application of a magnetic field [
18].
Despite all the advantages of thiourea- and magnetically functionalized chitosan adsorbents, their production involves the use of a large amount of water, making it necessary to evaluate the use of optimization methodologies in this process. Thus, the design of an industrial-scale process using mass integration techniques can reduce the generation of wastewater and the consumption of freshwater [
19,
20]. Mass integration provides a holistic approach to the design, modernization, and operation of processes that emphasizes the unity of the process itself [
21]. By considering the interaction between process units, resources, flows, and objectives, process integration enables a fundamental understanding of the overall knowledge of the process, systematically determines achievable performance targets, and facilitates systematic decision making to achieve these targets [
22]. Recently, mass integration was implemented to minimize freshwater consumption in a yeast company through the optimization of the water utilization network using pinch methodology. This allowed for the recycling of 173,480.4 t/year of water within the system, demonstrating that 27% of freshwater usage can be saved [
23].
Different bioadsorbent production processes have been evaluated from an environmental perspective using the Waste Reduction Algorithm (WAR), which takes into account the Potential Environmental Impacts (PEI) associated with chemical processes [
24]. The evaluation of magnetite nanoparticle production using the co-precipitation method and the WAR algorithm revealed that the process does not generate significant negative impacts that could affect the ecosystem. It transforms feed streams with a high PEI into final products with a lower PEI, leading to the conclusion that large-scale magnetite production is environmentally applicable and respectful [
25]. Regarding the synthesis of titanium dioxide, magnetite, and/or thiourea-modified chitosan microbeads, it was found that the presence of monovalent alcohols, such as ethanol (used in washing stages) and propanol (formed as an undesired product during hydrolysis in nanoparticle formation), increases the environmental effects related to the Photochemical Oxidation Potential (PCOP) category. However, overall, the process showed a moderately good performance [
26].
This work aims to assess the environmental evaluation of the industrial-scale production of bioadsorbents from chitosan modified with iron nanoparticles and functionalized with thiourea as a chelating agent (BMCMN) combined with mass integration, using computer-aided engineering tools. This approach shows, for the first time, the use of computer-assisted process engineering (CAPE) strategies and the mass integration of processes to enhance the technical performance of a production topology for thiourea-modified chitosan bioadsorbents with magnetite, as well as the quantification of their generated environmental impacts, both upstream and under atmospheric and toxicological categories.
2. Materials and Methods
The methodology applied in this study consisted of the following stages: (i) mass integration, (ii) the process simulation of the integrated case, and (iii) an environmental evaluation using the Waste Reduction Algorithm (WAR). These stages were carried out to assess the proposal for the industrial production of bioadsorbents modified with chelating agents and magnetic nanoparticles (BMCMN). To develop the second stage, the Aspen Plus® software was used.
2.1. Process Production Modeling
Initially, the block diagram was designed based on the production process of BMCMN and the simulations were based on experimental data and previous information published by the research group [
27]. For the scaling of the project, the availability of the raw materials was taken into account, resulting in a design of a plant with a processing capacity of 1746.2 tons of chitosan per year [
28]. Objective flows were established, taking into account the experimentally reported flows that served as a basis for the process simulation. The block diagram for the designed topology is shown in
Figure 1.
For the synthesis of magnetic nanoparticles using the co-precipitation method of FeCl
3·6H
2O and FeCl
2·4H
2O, which were initially fed into the process in a molar ratio of 2:1 [
29], they were fed with mass flow rates of 219.298 kg/h and 79.169 kg/h, respectively, and diluted in 4642.51 kg/h of water. The resulting solutions were mixed together via mechanical stirring (MIX-1), then heated to 80 °C [
30], and sent to the reaction vessel, where 3505.58 kg/h of a 3M NaOH solution was added, raising the pH of the medium to approximately 12, allowing for the formation and precipitation of magnetic nanoparticles with a 98% conversion rate [
31]. The resulting slurry was cooled to 28 °C and sent to magnetic separation, where 725.915 kg/h of nanoparticles were collected and separated from 7722.65 kg/h of the resting solution (wastewater). They were washed in three stages with 21,777.81 kg/h of water (wash 1) and 1842.89 kg/h of ethanol (wash 2 and 3) [
31], and then dried in an oven at 105 °C and 1 bar [
32]. It has been reported that the magnetic nanoparticles obtained using this method have a size distribution in the range of 10–500 nm, with a high purity [
33].
For the synthesis of BMCMN, 180.71 kg/h of chitosan was initially dissolved in 9071.79 kg/h of a solution of acetic acid at 28 °C and 1 bar of pressure (Mixing 2). Subsequently, the solubilized chitosan was mixed with 90.36 kg/h of thiourea (Mixing 3), and the resulting stream (functionalized chitosan) was mixed with the synthesized magnetic nanoparticles (Mixing 4) for the formation of BMCMN under the same standard conditions (28 °C and 1 bar). The synthesized bioadsorbent was washed with water twice (wash 4 and 5), then dried at 105 °C to remove the excess of water, and cooled to 28 °C for storage.
For the study, the process was divided into two categories, as shown in
Table 1:
2.2. Simulation Procedure
The production of BMCMN was simulated using the commercial package Aspen Plus v12 (Bedford, MAX, USA), since the processes involve substances in the solid, liquid, and gas states; further, it has shown satisfactory results in the modeling and simulation of these type of processes, capturing the behavior [
35,
36]. The requirements for the consumables, utilities, and energy needs were calculated for the production of 1746.2 tons of chitosan per year. To accomplish this, the substances involved in the processes were initially selected from the simulator’s chemical database. Subsequently, the thermodynamic model was chosen based on the algorithm described by Carlson [
37]. Considering the aforementioned, the process conditions, and the chemical substances involved, a selection of appropriate models for predicting the thermodynamic properties was recommended. Among these models, ELECNRTL, a modified version of the Non-Random Two Liquids (NRTL) model for electrolytes, was suggested. This choice was motivated by the presence of both weak electrolytes, such as acetic acid, and strong electrolytes, such as NaOH and NaCl, in the system. Next, the conditions of the streams, including their mass, pressure, and temperature, were inputted. The necessary equipment for the different stages of the processes was selected, and the results were verified to ensure consistency with experimental data.
2.3. Mass Integration of Water Streams
Mass integration refers to a systematic methodology that provides a fundamental understanding of the overall mass flow within a process, enabling the identification of performance targets and the optimization of species generation and routing throughout the process [
22]. A pinch analysis for water is a convenient tool for determining the minimum amount of water to be used, considering the recirculation of wastewater streams back into the process. In essence, this analysis takes into account the water requirements, both in terms of quantity and quality, for each process in the system. Quality is represented by the concentration of critical contaminants in the stream and the contaminant load that can enter each stage [
38]. These values can be established using algebraic approximations [
39].
In the pinch analysis methodology, two important terms are used: the streams of residual water that can be utilized for reuse in the process are referred to as “sources”, and the processes that can receive the reused water are known as “sinks”. Firstly, the sinks were organized in ascending order based on their maximum allowable composition. Then, the sources were organized in ascending order according to their contaminant composition. The loads of the sinks and sources were calculated, taking into account Equations (2) and (3).
Subsequently, the accumulated loads were calculated by summing the individual loads of the sinks and sources, for each interval k, as shown in Equations (4)–(6). The accumulated loads were organized in ascending order, and load-interval diagrams (LID) were developed. Cascade diagrams were constructed, and residual streams were calculated for each interval. The most negative residual was set as the minimum target for freshwater consumption. Finally, the cascade diagram was reviewed, adding the maximum residual to the first interval and setting the residual interval to zero as the critical point for material recycling. The minimum waste discharge objective was determined as the residual at the outlet of the last interval [
39].
For each interval k, the load and flow rates of the sources and sinks were calculated using Equations (4)–(6), respectively [
21].
where ΔM
k is the interval load, ΔW
k is the source flow per interval, and ΔG
k is the sink flow per interval
On the other hand, the construction of the cascade diagrams was carried out using the data obtained from the LID.
Figure 2a illustrates the flow balance for a given interval k, and the compilation of different intervals generated cascade diagrams, as shown in
Figure 2b,c. The most negative value of the residual (
) represents the objective for minimizing the consumption of freshwater to eliminate the existing infeasibilities. To achieve this, a fresh resource flow equal to
was added at the input of interval 1 in the revised cascade diagram. The residual obtained in the last interval represents the objective for minimizing the generation of wastewater [
39]. The described procedure provides the necessary information for the design of substance distribution networks.
The mass integration of water streams in the synthesis process of BMCMN was carried out using the base case [
40] as a starting point. The objective function was defined as the minimization of the freshwater requirements and wastewater effluents for the two stages considered in
Section 2.1. Based on the flow and composition data obtained in the simulation, the sources and sinks to be considered for the integration were selected, taking into account the incoming freshwater streams and the residual streams predominantly composed of water. The critical contaminant in the effluents was identified and the addition of treatment units was established based on it. Subsequently, load-interval diagrams (LID) and cascade diagrams were generated, which were corrected if they had any negative residuals, allowing for the achievement of the minimum objectives for freshwater and wastewater. Finally, water distribution networks were designed for each studied subprocess, indicating which effluents should be recirculated to the different processes and in what quantity, to meet the minimum objectives established in the pinch analysis.
2.4. Environmental Assessment
The Waste Reduction Algorithm (WAR) (developed by the National Risk Management Research Laboratory of the U.S. Environmental Protection Agency) was used for the environmental impact assessment [
41]. The potential environmental impact (PEI) of a given quantity of material and energy is a conservation relationship regarding the potential environmental impact based on the input and output impact flows of a process. The WAR–GUI software was used for calculating the PEI generated, considering eight environmental impact categories, four of which were related to atmospheric effects, while the remaining ones referred to toxicological effects. Within the atmospheric impact indices, there were two global categories: global warming potential (GWP) and ozone depletion potential (ODP), and two regional categories: photochemical oxidation potential (PCOP) and acidification potential (AP). On the other hand, the toxicological impact indices could be classified into human and ecological categories. The human indices included human toxicity potential by ingestion (HTPI) and that by inhalation or dermal exposure (HTPE), while the ecological indices corresponded to aquatic toxicity potential (ATP) and terrestrial toxicity potential (TTP) [
24]. In this work, Total Potential Environmental Impacts were used for the comparison of the processes, because they involved all the categories and were calculated by multiplying the mass flow rate by its chemical potential.
The first index, PEI
OUT, refers to the potential effects emitted to the environment external to the process and is useful for determining if a plant can produce a desired product while generating minimal environmental impact. The other index, PEIGEN, refers to the impacts generated by the system, and this indicator allows us to determine if the process was internally environmentally efficient. These PEIs can be calculated per unit of time using Equations (7) and (8).
Additionally, the potential environmental impacts per unit mass of the products can be expressed using Equations (7) and (8).
where
and
are the input rates of the PEI to the chemical process,
and
are the input and output rates of the PEI for the energy generation of the process,
and
are the PEIs related to the residual energy;
is the mass flow rate of the stream
j,
is the mass fraction of the component
k in the stream
j, and
is the PEI of the
k specie [
24].
For the environmental assessment of the process, missing chemical species were initially added to the software database along with their environmental parameters. Then, the chemical species involved in the process that were already present in the database were selected. If they were not available, their CAS number, chemical formula, exact chemical formula, molecular weight, and toxicological data from the Material Safety Data Sheet (MSDS) were provided. With the simulation results from Aspen Plus, the compositions of each stream were identified, thus quantifying the generated waste and products. Once all the necessary information was provided to the software, the amounts of generated and emitted global potential impacts, both per unit mass and per unit time, for each category of the waste reduction algorithm, were obtained. The calculations also took into account the effect of the type of fuel on the environmental impacts of the processes, which were represented in bar graphs and data tables.
4. Conclusions
The environmental assessment using the Waste Reduction Algorithm (WAR) for a mass-integrated approach to the industrial production of bioadsorbents modified with chelating agents and magnetic nanoparticles (BMCMN) was carried out to determine its Potential Environmental Impacts (PEI). From the simulation, it was found that the mass-integrated topology achieved a product flow of 2029.89 t/year, with a processing capacity of 5503.36 t/year. The application of mass integration through the design of water distribution networks based on a pinch analysis enabled a recirculation of 130.6 t/h, reducing the consumption of freshwater by 49.5. This allowed for a decrease of 40.9 in wastewater compared to the base case.
The processes evaluated in this research demonstrated a good environmental performance, as evidenced by the negative total generated impacts. The optimization of the process in terms of water usage did not have significant effects on the environmental and process safety indicators. From a toxicological and atmospheric perspective, it was found that the HTPI and TTP categories presented the highest PEI indices, which was attributed to the handling of substances such as FeCl₂·4H₂O, FeCl₃·6H₂O, and thiourea. In this regard, ethanol had effects on the atmospheric categories, with the PCOP category contributing approximately 70% to the environmental impacts from these sources. The process generated fewer environmental impacts than similar processes, such as the production of bioadsorbents modified with TiO2 nanoparticles and magnetite-TiO2 (1870 PEI/h vs. 4000 PEI/h vs. 4970 PEI/h, respectively). In terms of process safety, the associated indices remained constant or varied in a way that kept the index unchanged.