Home Evaluation of the operational performance of Karbala waste water treatment plant under variable flow using GPS-X model
Article Open Access

Evaluation of the operational performance of Karbala waste water treatment plant under variable flow using GPS-X model

  • Noor Abdulhussein Khudhair EMAIL logo , Basim K. Nile and Jabbar Hammoud Al-Baidani
Published/Copyright: March 4, 2024
Become an author with De Gruyter Brill

Abstract

As Karbala City is a religious tourism destination, millions of tourists visit the city annually, so there is a high fluctuation of flow in the wastewater between the plants during the year. The aim of this study is to evaluate the performance of the wastewater treatment plant (WWTP) in Karbala, Iraq, in removing pollutants for different flow rates in five scenarios using the GPS-X model. The most important phase in modeling, which greatly impacts simulation accuracy, is characterizing the influent composition to meet the mass balance. As a result, the influent wastewater was initially described and thoroughly examined. The model has been calibrated, followed by the collected data’s validation. The sensitivity of different stoichiometric and kinetic factors in the GPS-X was examined and screened to calibrate the model. To demonstrate the consistency between the simulated and measured data, the route mean square error was used in this instance. The result showed that the Karbala WWTP complies with Iraqi environmental regulations for water discharged to surface water or water for other uses and has an appropriate efficiency of wastewater treatment even if the flow entering the plant reaches the peak flow rate of 180,000/day. An improvement in orthophosphate removal efficiency was observed as the flow rate increased because of the contact time in the anaerobic basins, which gradually decreased as the flow increased until it reached the appropriate time for PO 4 3 removal. The outcomes of the present study provide an impression to the operators of the treatment plant of the impact of fluctuating flow on the treatment plant. The developed model can also be used for future studies.

1 Introduction

The volume of wastewater entering the treatment facilities has increased due to urbanization and industrial growth. Wastewater contains substances that are harmful to both people and the environment. As a result, dumping untreated wastewater in natural water bodies creates serious sustainability issues, like endangering aquatic and terrestrial life and raising the expense of treating the polluted water [1]. Consistency with discharge characteristic criteria, harm inhibition on the environment and human health, feasible water reuse, and recovery concerns connected to energy and materials (e.g., nutrients) make proper wastewater management a significant global task [2]. To adhere to environmental regulations, enterprises that produce pollutants must effectively treat their wastewater before releasing it to maintain the long-term viability of their water resources [3]. Activated sludge-based systems are currently the most popular choice since they are the most cost-effective and efficient of the several biological wastewater treatment methods used today. Severe weather might bring either dry or rainy spells, and the population expansion will create new difficulties for operating a wastewater treatment facility. The wastewater treatment process (WWTP) can be simulated under various intake circumstances to prepare for these issues. In harsh climatic circumstances, the feed stream flow rates may exceed the intended maximum values, and the compositions may fluctuate greatly [4].

It is crucial to employ simulation to make the most of the model’s predictive skills, which may be used to quickly and effectively weed out the optimal design option and cut down on the time and expense of laboratory testing. Modeling is the process of simplification of actual representation. Mathematical operations and equations used time-dependent variables and parameters to define the model [5]. Recent times have seen a rise in dynamic modeling and simulation in wastewater treatment. Many models are created to enhance activated sludge. Several software programs, such as SIMBA, GPS-X, AQUASIM, BioWin, STOAT, FOR, and WEST that promote dynamic modeling of wastewater treatment facilities use these models. Tools for simulation and modeling are used to evaluate process procedures, optimize designs, and analyze costs [6]. The calibration and validation steps are crucial for these models to give accurate results. To precisely manage and calibrate the plant’s operation, the stoichiometry of the reactions inside the reactor controls the materials produced and consumed. To investigate the microbial transformation process in removing organic matter and nutrients, kinetics and stoichiometric measurements in the titration and verification procedures were considered necessary [7].

The GPS-X model has been used in a wide range of research. Faris et al. [5] investigated the medium volume, dissolve oxygen content, and sludge return to the moving bed biofilm reactor by using sensitivity analysis. The results demonstrated that the nitrification process is impacted when the dissolved oxygen content exceeds 3 mg/L in the ammonia-rich side streams generated from the rejected water of the anaerobic digester. Alwardy et al. [8] evaluated Al-muamirah WWTP in Hilla City, the results revealed that the facility has an acceptable level of wastewater treatment efficiency and produces water that complies with Iraqi environmental standards for water discharged to surface water or water for other uses. Kobeyeb et al. [9] Presented a study to design and model an on-site greywater treatment system for a hotel building in Los Angeles, California. Many greywater treatment plant options were rigorously examined and modeled using the GPS-X software. The MBR plant was the best option for considering location and standards and was advised for use in hotel structures. Hammed et al. [10] used the GPS-X model to investigate how to lower nutrient concentrations in a full-scale plant. The findings demonstrated that the lack of rbCOD in significant quantities to aid in reducing nitrates rendered having a proportion of internal recycle (IR) of 3% worthless. Cao et al. [7] introduced an innovative way using GPS-X and response surface methodology to improve the removal of total nitrogen (TN) in WWTPs. 61 parameters’ sensitivities were checked and examined. The findings demonstrated that the denitrification rate was significantly impacted by the DO concentration that diffused into various biological compartments. SRT and TN elimination go hand in hand. Key parameter relevance and optimization orders were examined.

Due to the city of Karbala’s popularity as a destination for religious travelers, there is a significant variation in wastewater flow into the plant. The present study aims to study the effect of variation in flow on the operational efficiency of the Karbala WWTP by using five scenarios with different flow rates (minimum flow to peak flow) in GPS-X software to model the plant. Additionally, the optimal operation of the WWTP was suggested for each scenario of different flow rates depending on the model results after calibration and validation of the model.

2 Materials and methods

2.1 Karbala WWTP

The WWTP in the city of Karbala serves 2.5 million people. The conventional activated sludge method is vital to Karbala’s WWTPs. Four sewage treatment facilities utilizing the conventional activated sludge system of type A2/O are used in an integrated wastewater treatment project in Karbala City. Each facility has a daily discharge capacity of 100,000 m3. This plant’s geographical coordinates are 32.525590° North and 44.074909° East. Figure 1 shows a satellite.00000000000e image of the Karbala WWTP.

Figure 1 
                  Google Earth satellite view of the Karbala WWTP.
Figure 1

Google Earth satellite view of the Karbala WWTP.

The plant involves five steps of treatment. The first step is preliminary (physical) treatment, which includes four units, which are coarse and fine screens used for removing coarse and fine solid materials. A grit and oil removal chamber is used to remove sand and oil by a physical machine. The next step is primary treatment, carried out by four primary sedimentation basins. After that, there is the secondary treatment, which involves four tanks. There are two anaerobic reactors with a volume of 8,736 m3 for both of them, two anoxic reactors with a volume of 14,112 m3, eight aeration reactor tanks with a volume of 54,054 m3, and eight sedimentation tanks with a surface area of 6,432 m2. The next step is tertiary treatment, which includes chemical disinfection through the chlorination tank with a surface area of 3,000 m3. The last stage is sludge treatment, which involves two gravity thickeners with a surface area of 400 m2, three mechanical thickeners with a surface area of 60 m2, four anaerobic digesters, and 60 cells of drying beds with a surface area of 50,000 m2.

2.2 Collection of samples

The purpose of the present study was to study the effect of variable flow on the operational performance of Karbala WWTP. Influent sewage has been analyzed for its properties to measure contaminants’ concentrations prior to treatment. Moreover, concentrations of effluent contaminants following treatment have been evaluated and compared with the Iraqi standard limits. Four tests for each pollutant were conducted monthly for 1 year, from January 2022 to December 2022. According to the Standard Procedures for the Examination of Water and Wastewater, the tests were conducted in the laboratory of the Karbala Sewage Department (APHA 2017). The biological oxygen demand (BOD5), chemical oxygen demand (COD), total suspended solids (TSS), ammonia-nitrogen ( NH 3 + N ), orthophosphate ( PO 4 3 ), nitrate ( NO 3 -N ), and nitrite ( NO 2 -N ) are measured. Table 1 displays the concentrations of contaminants in the influent and effluent wastewater.

Table 1

Concentrations of influent and effluent sewage of Karbala WWTP

Parameter Influent, mg/l Effluent, mg/l Iraqi standard
COD 250 23 100
BOD5 115 4 40
TSS 140 11 60
NO3 0 16 8
NO2 0 0.2
NH 3 + 28 0.5 10
PO4 3− 4.5 2.5 3
H2S 30 0.3 3
DO 0 3.5
SO4 850 860 600

Note: all parameter units are in (mg/l).

2.3 Modeling of Karbala WWTP using GPS-X model

Mathematical modeling has emerged as an essential instrument for sustainable wastewater management, particularly in simulating intricate procedures inherent in activated sludge processes. Biological systems are very complicated. Temperature, flow fluctuation, the quantity of the wastewater pollutant, and the operational circumstances are just a few of the many variables significantly affecting treatment efficiency. The modeling technology company Hydromantis’ GPS-X program, version 8.0.1 (education license), was utilized in the present study. GPS-X is a multipurpose modeling platform designed to simulate WWTPs. This allows for the interactive and dynamic examination of the complicated relationships between the plant’s many unit processes [10]. It was constructed with integrated biological WWTP and many other processes that involve physical and chemical reactions. Data were obtained from the WWTP of Karbala and fed into the GPS-X program. A comprehensive carbon, nitrogen, and phosphorus pH (mantis2lib) library was selected. The dimensions of each unit and the input concentration were entered into the model. The assumptions used to create the model in this study are: the pH is stable and close to neutral, and the biological process takes place at a temperature of 20°C (Figure 2).

Figure 2 
                  Schematic diagram of Karbala WWTP.
Figure 2

Schematic diagram of Karbala WWTP.

To establish the best match between simulation results and actual plant effluent, the model was calibrated by adjusting stoichiometric, kinetic, and other important parameters to fit the simulation. After that, mathematical validation was used to check the model, and the following equations illustrate the results for root mean square error (RMSE) and correlation coefficient (R). The reasonable limits of these statistical criteria are 1 ≥ R > 0.8 and 0 ≤ RMSE < 1.5.

(1) R = ( C o C o ¯ ) ( C p C p ¯ ) ¯ σ c o σ c p ,

(2) RMSE = ( C o C p ) 2 ¯ C o ¯ C p ¯ ,

where C o is the actual data, C p is the modeled data, C o is the average of actual data, C p is the average of modeled data, and σ is the standard deviation over the dataset. Many simulation runs were carried out under different flow conditions to learn how varying flows may affect the plant’s efficiency. The steps listed above, illustrated in Figure 3, are used to simulate the modeling of the plant [11].

Figure 3 
                  Flow chart of the present study.
Figure 3

Flow chart of the present study.

3 Results and discussion

The average monthly flow entering the WWTP is 60,000 m3/day, under the design value of the plant, which has a value of 100,000 m3/day. There is little variation in flow during the year; the flow reaches its peak value during 2 months of the year because of the visitors coming to the city to revive religious rituals. The average monthly influent concentrations of TSS, COD, BOD5, NH4 +, and PO 4 3 are 140, 250, 115, 28, and 4.5, and there is a variation in the parameters mentioned earlier during the year. This fluctuation is due to seasonal variations in wastewater concentrations and industrial wastewater from industrial regions that enter the Karbala wastewater treatment system through the municipal wastewater system.

3.1 Modeling of Karbala WWTP by using GPS-X program

The Karbala WWTP was modeled using GPS-X by utilizing five different scenarios.

3.2 Model calibration

In order to evaluate the effectiveness of the model, it is recommended to perform calibration procedures utilizing organic or default data. The model was verified using effluent COD, effluent TSS, BOD, orthophosphate, ammonia nitrogen, nitrite, and nitrate as indicators. This is compared with the calibration of the Albarrakiya WWTP by GPS-X software done by ref. [12], which uses the data of COD, BOD, and TSS to calibrate the model. Calibration is necessary due to the system’s limited understanding of real physical, chemical, and biological processes and the desire to balance model simplicity with model accuracy. The model ran with the default data inside; hence, the outcomes were wildly inaccurate [13]. The model was calibrated with data on the average monthly pollutant concentration for 1 year in 2022. Compared with the study by Alwardy et al. [8], which uses data from 2 months to calibrate the Al-muamirah WWTP model by GPS-X extended from 1/2021 to 2/2021. Several critical parameters that influenced the alteration of outcomes were adjusted. Characterizing the influent wastewater is seen as the most important phase in the modeling process, requiring rigorous investigation. The laboratory data pertain to the influent and effluent data for an influent flow of 60,000 m3/day. This represents the second scenario, which involves adjustments to achieve minimal differences between the simulated and observed pollutants in the effluent wastewater. The model used for influent characterization is the carbon, nitrogen, phosphorus, and pH (mantis3lib) library. In this regard, the GPS-X already contains the default settings for the COD fractions; however, these default values have been modified to ensure better model calibration [14]. In Table 2, the significant influent raw wastewater portions are mentioned.

Table 2

Default and calibrated values of influent stoichiometry composition

Influent stoichiometry composition
Parameter Symbol Default value Calibrated value
Soluble inert fraction of the total COD frsi 0.05 (raw) 0.08 (primary) 0.065
Readily biodegradable fraction of the total COD frss 0.2 (raw) 0.32 (primary) 0.09
Particulate inert fraction of the total COD frxi 0.13 (raw) 0.12 (primary) 0.168

In this instance, the effluent TSS readings were found to be more than 30 mg/L, which was far from the result of the actual influent flow of 60,000 m3/day. As a result, the GPS-X model’s sensitive parameters with regard to the secondary sedimentation basin were changed. Some delicate parameters have been adjusted, including the maximum Vesilind settling velocity, which has gone from 410 to 981.95 m/day; the maximum settling velocity, which has been adjusted from 274 to 356 m/day; and the feeding point from the bottom, which has been changed from the default 1 m to the real 4.3 m. Along with other operational characteristics of a different dimension, these three parameters were the most sensitive.

The effluent COD reading was higher than the result of the plant settling correction factor for xii, and the settling correction factor for xbai was adjusted from their default values 1, 1 to 0.07, 0.26. Acetate fraction of total COD was changed from the default value of 0 to 0.102 to adjust the value of PO4 3−, whose reading was 5; after the calibration, the reading was 3.122. The default nitrate values were shown to be substantially lower than the actual quantities, but the adjustment resulted in fitting the projected nitrate concentration with reality after modifying some sensitive parameters on the kinetics of ammonia and nitrates, including (ammonium fraction of soluble TKN) and (Aerobic heterotrophic yield on soluble substrate). The effluent stoichiometry composition is shown in Table 3.

Table 3

Effluent stoichiometry composition

Effluent stoichiometry composition
Calibrated fraction Default Calibrated
Maximum settling velocity 274 365
Maximum Vesilind settling velocity 410 981.95
Settling correction factor for xii 1 0.07
Settling correction factor for xbai 1 0.26
Acetate fraction of total COD 0 0.102
Aerobic heterotrophic yield on soluble substrate 0.6666 0.36963

Figure 4 shows the technique to handle the calibration process. After calibration, the model is ready to be applied to the scenarios.

Figure 4 
                  Illustration diagram of calibration process.
Figure 4

Illustration diagram of calibration process.

3.3 Statistical analysis

The next task after the calibration of the model for a whole year in 2022 is to examine the calibration process using statistical analysis. Figure 5 shows the actual and simulated concentrations of each pollutant. The statistical equations used to examine the calibration process are the correlation coefficient (R) and RMSE.

Figure 5 
                  The calibration of the actual and the simulated values.
Figure 5

The calibration of the actual and the simulated values.

The correlation coefficient for each actual and simulated parameter was above 0.8, which means that there is a positive linear relationship so that as one variable increases, the other variable also increases proportionally. The RMSE for each parameter was close to zero, and that indicates that the model has more accurate predictions and matches the data well. Table 4 shows the correlation and the RMSE for each parameter of actual and simulated values for average monthly data for a period of 12 months during the year 2022.

Table 4

Value of R and RMSE of calibration process

Parameter R-value RMSE
TSS 0.92 0.011
BOD5 0.87 0.082
COD 0.85 0.021
NH 3 + -N 0.89 0.022
NO2–N 0.83 0.027
NO3–N 0.86 0.138
PO 4 3 0.81 0.011

3.4 Sensitivity analysis

In order to assess the potential effects of variation in flow on the output variables of BOD, COD, TSS, PO 4 3 , NH 3 + -N , NO 3 -N , and NO 2 -N of effluent, a basic sensitivity analysis was carried out as part of this study. This research is advantageous because it improves model prediction and reduces the number of model parameters that require calibration as per EPA standards [15].

Figures 6 and 7 demonstrate how the influent discharge affects the plant’s production. The hydraulic and organic loads on the plant’s reactors increase with the increase in the discharge. This has a negative effect on the treatment processes and increases the concentrations of BOD, COD, and TSS in the treated effluent wastewater. The BOD, COD, and TSS concentrations did not exceed the Iraqi standards even while the Karbala wastewater treatment facilities reached 180,000 m3/day, the plant’s peak flow. Reducing the levels of pollutants in the treated wastewater might be possible if the mass balance could be immediately altered in response to the quantity and quality realities. It has been demonstrated that phosphate removal efficiency increases as discharge increases. For the removal of phosphorus during rbCOD fermentation, acetate is required. For rbCOD fermentation, retention times of 0.25–1.0 h are sufficient [16].

Figure 6 
                  Effect of flow rate on effluent COD, BOD, and TSS concentrations.
Figure 6

Effect of flow rate on effluent COD, BOD, and TSS concentrations.

Figure 7 
                  Effect of flow rate on effluent concentrations of ortho-phosphate, ammonia nitrogen, nitrite, and nitrate.
Figure 7

Effect of flow rate on effluent concentrations of ortho-phosphate, ammonia nitrogen, nitrite, and nitrate.

3.5 Applied scenarios

Upon completion of the calibration processes of the model, the model is ready to be applied in the required scenario and estimate the plant’s response to the variable flow. Five scenarios were used that depended on variable flow. In the first scenario, the number of primary sedimentation basins was decreased from 4 to 2, which resulted in a reduction in energy consumption by 50% as well as maintenance and operating expenses. The number of secondary sedimentation basins was reduced from 8 to 3, which helped reduce energy consumption and maintenance costs by about 38%. In the second scenario, the number of primary and secondary clarifiers was reduced to half, which reduced 50% of the energy and cost of operation. While in the third scenario, the number of secondary clarifiers was reduced from 8 to 6. All the primary and secondary clarifiers were used in the fourth and fifth scenarios.

Table 5 demonstrates that the COD and BOD removal efficiencies for each scenario’s examined under different flows were over 90%, thus demonstrating that the organic matter was removed highly effectively. Microorganisms, particularly heterotrophic bacteria, assisted in removing organic materials once the appropriate dissolved oxygen concentrations were present and there was enough mixing inside the reactor to cause these compounds to break down and turn into fixed substances [17]. When the flow increases, COD and BOD readings increase; this increase could result from a shorter detention time period. The time wastewater is allowed to remain in the treatment system is known as the detention time. Shorter detention times occur when the influent flow is increased because the wastewater is treated more quickly. As a result, the microorganisms in charge of breaking down organic waste have less time to work, which raises the COD and BOD levels in the effluent, as indicated in Figures 8 and 9.

Table 5

Effect of fluctuation of flow on effluent COD and BOD concentrations

Flow rate, m3/day 40,000 60,000 100,000 140,000 180,000
COD mg/L 24.36 26.97 31.22 35.65 40.86
COD, removal efficiency 90% 89% 87% 85% 83%
BOD mg/L 3.807 4.155 5.434 6.748 7.829
BOD, removal efficiency 96% 96% 95% 94% 93%
Figure 8 
                  The variation in the effluent concentrations of COD with various flow rates.
Figure 8

The variation in the effluent concentrations of COD with various flow rates.

Figure 9 
                  The variation in the effluent concentrations of BOD with various flow rates.
Figure 9

The variation in the effluent concentrations of BOD with various flow rates.

The increasing flow rate through the settling tanks leads to decreased settling efficiency. Table 6 demonstrates that TSS values also increase with the increase in flow rate. As a result, the amount of suspended particles in the effluent may increase, increasing the TSS values, which can be attributed to shorter retention times. Higher TSS levels in the effluent might result from shorter contact times between the wastewater and the treatment processes, which can prevent suspended particles from settling. For all scenarios, the secondary sedimentation basins achieved a removal efficiency of greater than 90% of the suspended solids. The removal of nutrients was impacted by the very low concentration of rbCOD in the wastewater of the city of Karbala [18]. Figure 10 shows the increase in TSS concentrations with the increase in the plant flow rates.

Table 6

Effect of fluctuation of flow on effluent TSS concentration

Flow rate, m3/day 40,000 60,000 100,000 140,000 180,000
TSS, mg/L 7.9 11.6 17.23 22.05 27.37
TSS removal efficiency 94% 91% 87% 84% 80%
Figure 10 
                  The variation in the effluent concentrations of TSS with the variation in flow rates.
Figure 10

The variation in the effluent concentrations of TSS with the variation in flow rates.

The NH 3 + -N values change for each scenario, as shown in Figure 11. The value of NH 3 + -N reaches 3.835 when the flow is 180,000 m3/day because the HRT decreases due to the increase in flow rate. With less time available for treatment, the nitrification process, which converts ammonium ( NH 3 + -N ) to nitrate ( NO 3 -N ), may be incomplete. As a result, higher ammonium levels can be observed in the effluent. Another reason for increasing the concentrations of NH 3 + -N is that when the influent flow increases, the contact time between the wastewater and the air or aeration mechanism tends to be reduced. This can limit the stripping of ammonia, leading to higher NH 3 + -N concentrations in the effluent. As the effluent concentration of ammonia nitrogen increases with the increase in the flow, the removal efficiency decreases. It was 98% when the flow was at its minimum value, which reduced the efficiency to 86% when it increased to 180,000 m3/day, as shown in Table 7. The nitrification process was functioning effectively at the Karbala WWTP, with ammonia being oxidized by over 95% for minimum and actual flow. The denitrification process at the Karbala WWTP is experiencing suboptimal performance due to an inadequate supply of appropriate organic substrate.

Figure 11 
                  The variation in effluent concentrations of 
                        
                           
                           
                              
                                 
                                    NH
                                 
                                 
                                    3
                                 
                                 
                                    +
                                 
                              
                              -N
                           
                           {\text{NH}}_{3}^{+}\text{-N}
                        
                      with different values of flow rates.
Figure 11

The variation in effluent concentrations of NH 3 + -N with different values of flow rates.

Table 7

Effect of fluctuation of flow on effluent NH 3 + -N concentrations

Flow rate, m3/day 40,000 60,000 100,000 140,000 180,000
NH 3 + -N , mg/L 0.329 0.41 0.7273 1.491 3.835
NH 3 + -N removal efficiency 98% 98% 97% 94% 86%

The potential impact of nitrogen and phosphorous removal on the high levels of dissolved oxygen discharged through IR and return activated sludge flow rates should also be considered [19]. The primary source of nutrition for phosphorous bacteria is rbCOD. Additionally, it serves as a source of feeding and functions as an electron donor in nitrate reduction. The competition for removing phosphates and nitrates by this compound may result in an uneven distribution of these nutrients. This parameter is highly responsive to the depletion of nutrients [20]. Table 8 shows the removal efficiency of phosphorous of the Karbala WWTP is considered low for minimum and actual flow, 31 and 29%, and this is in agreement with the study by Hammed et al. [10]. An increase in rbCOD concentration correlates with improved efficiency in phosphate removal and enhanced nitrification and denitrification processes. Due to the low concentrations present, it is necessary to incorporate an external carbon source to enhance the nutrient removal processes at the Karbala WWTP. The introduction or increase in rbCOD levels has the potential to impact the BOD and COD concentrations present in the plant’s effluent.

Table 8

Effect of fluctuation of flow rates on effluent PO 4 3 concentrations

Flowrate m3/day 40,000 60,000 100,000 140,000 180,000
PO 4 3 mg/L 3.03 3.122 1.543 0.932 0.66
PO 4 3 removal efficiency 31% 29% 67% 81% 86%

It was noticed that when the discharge increases, the nitrates and nitrites are only slightly influenced. Figures 12 and 13 show that, unlike NO 2 -N , the concentration of NO 3 -N decreases when the flow increases (Figure 14).

Figure 12 
                  The variation in the effluent concentrations of orthophosphate with different values of flow rates.
Figure 12

The variation in the effluent concentrations of orthophosphate with different values of flow rates.

Figure 13 
                  The variation in the effluent concentrations of nitrate with different values of flow rates.
Figure 13

The variation in the effluent concentrations of nitrate with different values of flow rates.

Figure 14 
                  The variation in the effluent concentrations of nitrite with different values flow rates.
Figure 14

The variation in the effluent concentrations of nitrite with different values flow rates.

4 Conclusion

The GPS-X program was successfully used in the present study to model the largest WWTP in Karbala, Iraq. The default settings of some model parameters, such as the kinetic and stoichiometric parameters and the influent wastewater characterization, have been calibrated. The developed model can be used for future studies because the model’s accuracy was acceptable, and the simulation result was close to the actual parameter concentrations of the plant. The results of the five scenarios with different flows show that the removal efficiency of the plant remains acceptable according to Iraqi standards when the influent flow is the minimum and actual flow for all pollutants except PO4 3−, which has a removal efficiency of 31% for minimum flow and 29% for basic flow. The removal efficiency is acceptable even if the influent flow reaches 180,000 m3/day, which is above the peak flow of the plant. In the future, it is suggested to decrease the number of primary clarifiers to two and secondary clarifiers to three when the influent flow is 40,000 m3/day or less and use two primary clarifiers and four secondary clarifiers when the influent flow to the treatment plant reaches 60,000 m3/day in summer or dry days of the year.

  1. Funding information: We declare that the manuscript was done depending on the personal effort of the author, and there is no funding effort from any side or organization.

  2. Conflict of interest: The authors state no conflict of interest.

  3. Data availability statement: Most datasets generated and analyzed in this study are in this submitted manuscript. The other datasets are available on reasonable request from the corresponding author with the attached information.

References

[1] Aghalari Z, Dahms HU, Sillanpää M, Sosa-Hernandez JE, Parra-Saldívar R. Effectiveness of wastewater treatment systems in removing microbial agents: a systematic review. Globalization Health. 2020;16:1–11.10.1186/s12992-020-0546-ySearch in Google Scholar PubMed PubMed Central

[2] Jafarinejad S, Jiang SC. Current technologies and future directions for treating petroleum refineries and petrochemical plants (PRPP) wastewaters. J Environ Chem Eng. 2019;7(5):103326.10.1016/j.jece.2019.103326Search in Google Scholar

[3] Ismail MS, Yahya MD, Auta M, Obayomi KS. Adsorption behavior of heavy metal ions by hybrid inulin-TEOS for water treatment. Civ Eng J. 2022;8(9):1787–98.10.28991/CEJ-2022-08-09-03Search in Google Scholar

[4] Lahdenperä E, Koiranen T. Effects of variable feeds in wastewater plants by means of process simulation. 11th ICARD, IMWA, MWD Conference – Risk to Opportunity; 2018.Search in Google Scholar

[5] Faris AM, Zwain HM, Hosseinzadeh M, Majdi HS, Siadatmousavi SM. Start-up and operation of novel EN-MBBR system for sidestreams treatment and sensitivity analysis modeling using GPS-X simulation. Alex Eng J. 2022;61(12):10805–18.10.1016/j.aej.2022.04.026Search in Google Scholar

[6] Hernandez-Sancho F, Sala-Garrido R. Cost modelling in waste water treatment processes: an empirical analysis for Spain. Dangerous pollutants (xenobiotics) in urban water cycle. Netherlands: Springer; 2008.Search in Google Scholar

[7] Cao J, Yang E, Xu C, Zhang T, Xu R, Fu B, et al. Model-based strategy for nitrogen removal enhancement in full-scale wastewater treatment plants by GPS-X integrated with response surface methodology. Sci Total Environ. 2021;769:144851.10.1016/j.scitotenv.2020.144851Search in Google Scholar PubMed

[8] Alwardy A, Alquzweeni SS, Al-Saadi RJ. Modelling and simulation of Al-muamirah wastewater treatment plant by GPS-X software. Kerbala J Eng Sci. 2021;1(2):113–30.10.1088/1755-1315/877/1/012027Search in Google Scholar

[9] Kobeyev S, Tokbolat S, Nazipov F, Satyanaga A. Design and modeling of an on-site greywater treatment system for a hotel building. Int J Build Pathol Adapt. 2023;41(1):201–24.10.1108/IJBPA-08-2021-0109Search in Google Scholar

[10] Hammed AN, Nile BK, Al-Baidani JH. Improving the performance of the Karbala wastewater treatment plant by using GPS-X Model. International Symposium on Multidisciplinary Studies and Innovative Technologies (ISMSIT); vol. 2022; 2022. IEEE.10.1109/ISMSIT56059.2022.9932706Search in Google Scholar

[11] Zwain HM, Nile BK, Faris AM, Vakili M, Dahlan I. Modelling of hydrogen sulfide fate and emissions in extended aeration sewage treatment plant using TOXCHEM simulations. Sci Rep. 2020;10(1):22209.10.1038/s41598-020-79395-8Search in Google Scholar PubMed PubMed Central

[12] Awad A, Jafar R, Ali K. Modeling Banias refinery treatment plant using the GPS-X. Tishreen Univ J Eng Sci Ser. 2019;41:6.Search in Google Scholar

[13] Mu’azu ND, Alagha O, Anil I. Systematic modeling of municipal wastewater activated sludge process and treatment plant capacity analysis using GPS-X. Sustainability. 2020;12(19):8182.10.3390/su12198182Search in Google Scholar

[14] Sadri Moghaddam S, Pirali MR. Modeling and calibration of a full-scale wastewater treatment plant using GPS-X model (A case study of Tehran). Numer Methods Civ Eng. 2021;5(4):67–76.10.52547/nmce.5.4.67Search in Google Scholar

[15] EPA. Environmental protection agency (EPA). Wastewater Technology Fact Sheet: Oxidation Ditches. Washington, DC, USA: US EPA; 2000.Search in Google Scholar

[16] Metcalf EAO, Mohammad B, Gregory B, Franklin LP, William S, H David Tchobanoglous, et al. Wastewater engineering: treatment and resource recovery. McGraw Hill Education; 2014.Search in Google Scholar

[17] Bankston E, Wang Q, Higgins BT. Algae support populations of heterotrophic, nitrifying, and phosphate-accumulating bacteria in the treatment of poultry litter anaerobic digestate. Chem Eng J. 2020;398:125550.10.1016/j.cej.2020.125550Search in Google Scholar

[18] Khursheed A, Gaur RZ, Sharma MK, Tyagi VK, Khan AA, Kazmi AA. Dependence of enhanced biological nitrogen removal on carbon to nitrogen and rbCOD to sbCOD ratios during sewage treatment in sequencing batch reactor. J Clean Prod. 2018;171:1244–54.10.1016/j.jclepro.2017.10.055Search in Google Scholar

[19] Stewart RD, Bashar R, Amstadt C, Uribe-Santos GA, McMahon KD, Seib M, et al. Pilot-scale comparison of biological nutrient removal (BNR) using intermittent and continuous ammonia-based low dissolved oxygen aeration control systems. Water Sci Technol. 2022;85(2):578–90.10.2166/wst.2021.630Search in Google Scholar PubMed

[20] Figdore BA, Stensel HD, Winkler M-KH. Bioaugmentation of sidestream nitrifying-denitrifying phosphorus-accumulating granules in a low-SRT activated sludge system at low temperature. Water Res. 2018;135:241–50.10.1016/j.watres.2018.02.035Search in Google Scholar PubMed

Received: 2023-09-21
Revised: 2023-11-04
Accepted: 2023-11-08
Published Online: 2024-03-04

© 2024 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Methodology of automated quality management
  3. Influence of vibratory conveyor design parameters on the trough motion and the self-synchronization of inertial vibrators
  4. Application of finite element method in industrial design, example of an electric motorcycle design project
  5. Correlative evaluation of the corrosion resilience and passivation properties of zinc and aluminum alloys in neutral chloride and acid-chloride solutions
  6. Will COVID “encourage” B2B and data exchange engineering in logistic firms?
  7. Influence of unsupported sleepers on flange climb derailment of two freight wagons
  8. A hybrid detection algorithm for 5G OTFS waveform for 64 and 256 QAM with Rayleigh and Rician channels
  9. Effect of short heat treatment on mechanical properties and shape memory properties of Cu–Al–Ni shape memory alloy
  10. Exploring the potential of ammonia and hydrogen as alternative fuels for transportation
  11. Impact of insulation on energy consumption and CO2 emissions in high-rise commercial buildings at various climate zones
  12. Advanced autopilot design with extremum-seeking control for aircraft control
  13. Adaptive multidimensional trust-based recommendation model for peer to peer applications
  14. Effects of CFRP sheets on the flexural behavior of high-strength concrete beam
  15. Enhancing urban sustainability through industrial synergy: A multidisciplinary framework for integrating sustainable industrial practices within urban settings – The case of Hamadan industrial city
  16. Advanced vibrant controller results of an energetic framework structure
  17. Application of the Taguchi method and RSM for process parameter optimization in AWSJ machining of CFRP composite-based orthopedic implants
  18. Improved correlation of soil modulus with SPT N values
  19. Technologies for high-temperature batch annealing of grain-oriented electrical steel: An overview
  20. Assessing the need for the adoption of digitalization in Indian small and medium enterprises
  21. A non-ideal hybridization issue for vertical TFET-based dielectric-modulated biosensor
  22. Optimizing data retrieval for enhanced data integrity verification in cloud environments
  23. Performance analysis of nonlinear crosstalk of WDM systems using modulation schemes criteria
  24. Nonlinear finite-element analysis of RC beams with various opening near supports
  25. Thermal analysis of Fe3O4–Cu/water over a cone: a fractional Maxwell model
  26. Radial–axial runner blade design using the coordinate slice technique
  27. Theoretical and experimental comparison between straight and curved continuous box girders
  28. Effect of the reinforcement ratio on the mechanical behaviour of textile-reinforced concrete composite: Experiment and numerical modeling
  29. Experimental and numerical investigation on composite beam–column joint connection behavior using different types of connection schemes
  30. Enhanced performance and robustness in anti-lock brake systems using barrier function-based integral sliding mode control
  31. Evaluation of the creep strength of samples produced by fused deposition modeling
  32. A combined feedforward-feedback controller design for nonlinear systems
  33. Effect of adjacent structures on footing settlement for different multi-building arrangements
  34. Analyzing the impact of curved tracks on wheel flange thickness reduction in railway systems
  35. Review Articles
  36. Mechanical and smart properties of cement nanocomposites containing nanomaterials: A brief review
  37. Applications of nanotechnology and nanoproduction techniques
  38. Relationship between indoor environmental quality and guests’ comfort and satisfaction at green hotels: A comprehensive review
  39. Communication
  40. Techniques to mitigate the admission of radon inside buildings
  41. Erratum
  42. Erratum to “Effect of short heat treatment on mechanical properties and shape memory properties of Cu–Al–Ni shape memory alloy”
  43. Special Issue: AESMT-3 - Part II
  44. Integrated fuzzy logic and multicriteria decision model methods for selecting suitable sites for wastewater treatment plant: A case study in the center of Basrah, Iraq
  45. Physical and mechanical response of porous metals composites with nano-natural additives
  46. Special Issue: AESMT-4 - Part II
  47. New recycling method of lubricant oil and the effect on the viscosity and viscous shear as an environmentally friendly
  48. Identify the effect of Fe2O3 nanoparticles on mechanical and microstructural characteristics of aluminum matrix composite produced by powder metallurgy technique
  49. Static behavior of piled raft foundation in clay
  50. Ultra-low-power CMOS ring oscillator with minimum power consumption of 2.9 pW using low-voltage biasing technique
  51. Using ANN for well type identifying and increasing production from Sa’di formation of Halfaya oil field – Iraq
  52. Optimizing the performance of concrete tiles using nano-papyrus and carbon fibers
  53. Special Issue: AESMT-5 - Part II
  54. Comparative the effect of distribution transformer coil shape on electromagnetic forces and their distribution using the FEM
  55. The complex of Weyl module in free characteristic in the event of a partition (7,5,3)
  56. Restrained captive domination number
  57. Experimental study of improving hot mix asphalt reinforced with carbon fibers
  58. Asphalt binder modified with recycled tyre rubber
  59. Thermal performance of radiant floor cooling with phase change material for energy-efficient buildings
  60. Surveying the prediction of risks in cryptocurrency investments using recurrent neural networks
  61. A deep reinforcement learning framework to modify LQR for an active vibration control applied to 2D building models
  62. Evaluation of mechanically stabilized earth retaining walls for different soil–structure interaction methods: A review
  63. Assessment of heat transfer in a triangular duct with different configurations of ribs using computational fluid dynamics
  64. Sulfate removal from wastewater by using waste material as an adsorbent
  65. Experimental investigation on strengthening lap joints subjected to bending in glulam timber beams using CFRP sheets
  66. A study of the vibrations of a rotor bearing suspended by a hybrid spring system of shape memory alloys
  67. Stability analysis of Hub dam under rapid drawdown
  68. Developing ANFIS-FMEA model for assessment and prioritization of potential trouble factors in Iraqi building projects
  69. Numerical and experimental comparison study of piled raft foundation
  70. Effect of asphalt modified with waste engine oil on the durability properties of hot asphalt mixtures with reclaimed asphalt pavement
  71. Hydraulic model for flood inundation in Diyala River Basin using HEC-RAS, PMP, and neural network
  72. Numerical study on discharge capacity of piano key side weir with various ratios of the crest length to the width
  73. The optimal allocation of thyristor-controlled series compensators for enhancement HVAC transmission lines Iraqi super grid by using seeker optimization algorithm
  74. Numerical and experimental study of the impact on aerodynamic characteristics of the NACA0012 airfoil
  75. Effect of nano-TiO2 on physical and rheological properties of asphalt cement
  76. Performance evolution of novel palm leaf powder used for enhancing hot mix asphalt
  77. Performance analysis, evaluation, and improvement of selected unsignalized intersection using SIDRA software – Case study
  78. Flexural behavior of RC beams externally reinforced with CFRP composites using various strategies
  79. Influence of fiber types on the properties of the artificial cold-bonded lightweight aggregates
  80. Experimental investigation of RC beams strengthened with externally bonded BFRP composites
  81. Generalized RKM methods for solving fifth-order quasi-linear fractional partial differential equation
  82. An experimental and numerical study investigating sediment transport position in the bed of sewer pipes in Karbala
  83. Role of individual component failure in the performance of a 1-out-of-3 cold standby system: A Markov model approach
  84. Implementation for the cases (5, 4) and (5, 4)/(2, 0)
  85. Center group actions and related concepts
  86. Experimental investigation of the effect of horizontal construction joints on the behavior of deep beams
  87. Deletion of a vertex in even sum domination
  88. Deep learning techniques in concrete powder mix designing
  89. Effect of loading type in concrete deep beam with strut reinforcement
  90. Studying the effect of using CFRP warping on strength of husk rice concrete columns
  91. Parametric analysis of the influence of climatic factors on the formation of traditional buildings in the city of Al Najaf
  92. Suitability location for landfill using a fuzzy-GIS model: A case study in Hillah, Iraq
  93. Hybrid approach for cost estimation of sustainable building projects using artificial neural networks
  94. Assessment of indirect tensile stress and tensile–strength ratio and creep compliance in HMA mixes with micro-silica and PMB
  95. Density functional theory to study stopping power of proton in water, lung, bladder, and intestine
  96. A review of single flow, flow boiling, and coating microchannel studies
  97. Effect of GFRP bar length on the flexural behavior of hybrid concrete beams strengthened with NSM bars
  98. Exploring the impact of parameters on flow boiling heat transfer in microchannels and coated microtubes: A comprehensive review
  99. Crumb rubber modification for enhanced rutting resistance in asphalt mixtures
  100. Special Issue: AESMT-6
  101. Design of a new sorting colors system based on PLC, TIA portal, and factory I/O programs
  102. Forecasting empirical formula for suspended sediment load prediction at upstream of Al-Kufa barrage, Kufa City, Iraq
  103. Optimization and characterization of sustainable geopolymer mortars based on palygorskite clay, water glass, and sodium hydroxide
  104. Sediment transport modelling upstream of Al Kufa Barrage
  105. Study of energy loss, range, and stopping time for proton in germanium and copper materials
  106. Effect of internal and external recycle ratios on the nutrient removal efficiency of anaerobic/anoxic/oxic (VIP) wastewater treatment plant
  107. Enhancing structural behaviour of polypropylene fibre concrete columns longitudinally reinforced with fibreglass bars
  108. Sustainable road paving: Enhancing concrete paver blocks with zeolite-enhanced cement
  109. Evaluation of the operational performance of Karbala waste water treatment plant under variable flow using GPS-X model
  110. Design and simulation of photonic crystal fiber for highly sensitive chemical sensing applications
  111. Optimization and design of a new column sequencing for crude oil distillation at Basrah refinery
  112. Inductive 3D numerical modelling of the tibia bone using MRI to examine von Mises stress and overall deformation
  113. An image encryption method based on modified elliptic curve Diffie-Hellman key exchange protocol and Hill Cipher
  114. Experimental investigation of generating superheated steam using a parabolic dish with a cylindrical cavity receiver: A case study
  115. Effect of surface roughness on the interface behavior of clayey soils
  116. Investigated of the optical properties for SiO2 by using Lorentz model
  117. Measurements of induced vibrations due to steel pipe pile driving in Al-Fao soil: Effect of partial end closure
  118. Experimental and numerical studies of ballistic resistance of hybrid sandwich composite body armor
  119. Evaluation of clay layer presence on shallow foundation settlement in dry sand under an earthquake
  120. Optimal design of mechanical performances of asphalt mixtures comprising nano-clay additives
  121. Advancing seismic performance: Isolators, TMDs, and multi-level strategies in reinforced concrete buildings
  122. Predicted evaporation in Basrah using artificial neural networks
  123. Energy management system for a small town to enhance quality of life
  124. Numerical study on entropy minimization in pipes with helical airfoil and CuO nanoparticle integration
  125. Equations and methodologies of inlet drainage system discharge coefficients: A review
  126. Thermal buckling analysis for hybrid and composite laminated plate by using new displacement function
  127. Investigation into the mechanical and thermal properties of lightweight mortar using commercial beads or recycled expanded polystyrene
  128. Experimental and theoretical analysis of single-jet column and concrete column using double-jet grouting technique applied at Al-Rashdia site
  129. The impact of incorporating waste materials on the mechanical and physical characteristics of tile adhesive materials
  130. Seismic resilience: Innovations in structural engineering for earthquake-prone areas
  131. Automatic human identification using fingerprint images based on Gabor filter and SIFT features fusion
  132. Performance of GRKM-method for solving classes of ordinary and partial differential equations of sixth-orders
  133. Visible light-boosted photodegradation activity of Ag–AgVO3/Zn0.5Mn0.5Fe2O4 supported heterojunctions for effective degradation of organic contaminates
  134. Production of sustainable concrete with treated cement kiln dust and iron slag waste aggregate
  135. Key effects on the structural behavior of fiber-reinforced lightweight concrete-ribbed slabs: A review
  136. A comparative analysis of the energy dissipation efficiency of various piano key weir types
  137. Special Issue: Transport 2022 - Part II
  138. Variability in road surface temperature in urban road network – A case study making use of mobile measurements
  139. Special Issue: BCEE5-2023
  140. Evaluation of reclaimed asphalt mixtures rejuvenated with waste engine oil to resist rutting deformation
  141. Assessment of potential resistance to moisture damage and fatigue cracks of asphalt mixture modified with ground granulated blast furnace slag
  142. Investigating seismic response in adjacent structures: A study on the impact of buildings’ orientation and distance considering soil–structure interaction
  143. Improvement of porosity of mortar using polyethylene glycol pre-polymer-impregnated mortar
  144. Three-dimensional analysis of steel beam-column bolted connections
  145. Assessment of agricultural drought in Iraq employing Landsat and MODIS imagery
  146. Performance evaluation of grouted porous asphalt concrete
  147. Optimization of local modified metakaolin-based geopolymer concrete by Taguchi method
  148. Effect of waste tire products on some characteristics of roller-compacted concrete
  149. Studying the lateral displacement of retaining wall supporting sandy soil under dynamic loads
  150. Seismic performance evaluation of concrete buttress dram (Dynamic linear analysis)
  151. Behavior of soil reinforced with micropiles
  152. Possibility of production high strength lightweight concrete containing organic waste aggregate and recycled steel fibers
  153. An investigation of self-sensing and mechanical properties of smart engineered cementitious composites reinforced with functional materials
  154. Forecasting changes in precipitation and temperatures of a regional watershed in Northern Iraq using LARS-WG model
  155. Experimental investigation of dynamic soil properties for modeling energy-absorbing layers
  156. Numerical investigation of the effect of longitudinal steel reinforcement ratio on the ductility of concrete beams
  157. An experimental study on the tensile properties of reinforced asphalt pavement
  158. Self-sensing behavior of hot asphalt mixture with steel fiber-based additive
  159. Behavior of ultra-high-performance concrete deep beams reinforced by basalt fibers
  160. Optimizing asphalt binder performance with various PET types
  161. Investigation of the hydraulic characteristics and homogeneity of the microstructure of the air voids in the sustainable rigid pavement
  162. Enhanced biogas production from municipal solid waste via digestion with cow manure: A case study
  163. Special Issue: AESMT-7 - Part I
  164. Preparation and investigation of cobalt nanoparticles by laser ablation: Structure, linear, and nonlinear optical properties
  165. Seismic analysis of RC building with plan irregularity in Baghdad/Iraq to obtain the optimal behavior
  166. The effect of urban environment on large-scale path loss model’s main parameters for mmWave 5G mobile network in Iraq
  167. Formatting a questionnaire for the quality control of river bank roads
  168. Vibration suppression of smart composite beam using model predictive controller
  169. Machine learning-based compressive strength estimation in nanomaterial-modified lightweight concrete
  170. In-depth analysis of critical factors affecting Iraqi construction projects performance
  171. Behavior of container berth structure under the influence of environmental and operational loads
  172. Energy absorption and impact response of ballistic resistance laminate
  173. Effect of water-absorbent polymer balls in internal curing on punching shear behavior of bubble slabs
  174. Effect of surface roughness on interface shear strength parameters of sandy soils
  175. Evaluating the interaction for embedded H-steel section in normal concrete under monotonic and repeated loads
  176. Estimation of the settlement of pile head using ANN and multivariate linear regression based on the results of load transfer method
  177. Enhancing communication: Deep learning for Arabic sign language translation
  178. A review of recent studies of both heat pipe and evaporative cooling in passive heat recovery
  179. Effect of nano-silica on the mechanical properties of LWC
  180. An experimental study of some mechanical properties and absorption for polymer-modified cement mortar modified with superplasticizer
  181. Digital beamforming enhancement with LSTM-based deep learning for millimeter wave transmission
  182. Developing an efficient planning process for heritage buildings maintenance in Iraq
  183. Design and optimization of two-stage controller for three-phase multi-converter/multi-machine electric vehicle
  184. Evaluation of microstructure and mechanical properties of Al1050/Al2O3/Gr composite processed by forming operation ECAP
  185. Calculations of mass stopping power and range of protons in organic compounds (CH3OH, CH2O, and CO2) at energy range of 0.01–1,000 MeV
  186. Investigation of in vitro behavior of composite coating hydroxyapatite-nano silver on 316L stainless steel substrate by electrophoretic technic for biomedical tools
  187. A review: Enhancing tribological properties of journal bearings composite materials
  188. Improvements in the randomness and security of digital currency using the photon sponge hash function through Maiorana–McFarland S-box replacement
  189. Design a new scheme for image security using a deep learning technique of hierarchical parameters
  190. Special Issue: ICES 2023
  191. Comparative geotechnical analysis for ultimate bearing capacity of precast concrete piles using cone resistance measurements
  192. Visualizing sustainable rainwater harvesting: A case study of Karbala Province
  193. Geogrid reinforcement for improving bearing capacity and stability of square foundations
  194. Evaluation of the effluent concentrations of Karbala wastewater treatment plant using reliability analysis
  195. Adsorbent made with inexpensive, local resources
  196. Effect of drain pipes on seepage and slope stability through a zoned earth dam
  197. Sediment accumulation in an 8 inch sewer pipe for a sample of various particles obtained from the streets of Karbala city, Iraq
  198. Special Issue: IETAS 2024 - Part I
  199. Analyzing the impact of transfer learning on explanation accuracy in deep learning-based ECG recognition systems
  200. Effect of scale factor on the dynamic response of frame foundations
  201. Improving multi-object detection and tracking with deep learning, DeepSORT, and frame cancellation techniques
  202. The impact of using prestressed CFRP bars on the development of flexural strength
  203. Assessment of surface hardness and impact strength of denture base resins reinforced with silver–titanium dioxide and silver–zirconium dioxide nanoparticles: In vitro study
  204. A data augmentation approach to enhance breast cancer detection using generative adversarial and artificial neural networks
  205. Modification of the 5D Lorenz chaotic map with fuzzy numbers for video encryption in cloud computing
  206. Special Issue: 51st KKBN - Part I
  207. Evaluation of static bending caused damage of glass-fiber composite structure using terahertz inspection
Downloaded on 18.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/eng-2022-0558/html
Scroll to top button