Open access peer-reviewed chapter - ONLINE FIRST

# Simulation and Optimization of an Integrated Process Flow Sheet for Cement Production

By Oluwafemi M. Fadayini, Adekunle A. Obisanya, Gloria O. Ajiboye, Clement Madu, Tajudeen O. Ipaye, Taiwo O. Rabiu, Shola J. Ajayi and Joseph T. Akintola

Submitted: September 17th 2020Reviewed: November 27th 2020Published: March 9th 2021

DOI: 10.5772/intechopen.95269

## Abstract

In this study the process flow diagram for the cement production was simulated using Aspen HYSYS 8.8 software to achieve high energy optimization and optimum cement flow rate by varying the flow rate of calcium oxide and silica in the clinker feed. Central composite Design (C.C.D) of Response Surface Methodology was used to design the ten experiments for the simulation using Design Expert 10.0.3. Energy efficiency optimization is also carried out using Aspen Energy Analyser. The optimum cement flow rate is found from the contour plot and 3D surface plot to be 47.239 tonnes/day at CaO flow rate of 152.346 tonnes/day and the SiO2 flow rate of 56.8241 tonnes/day. The R2 value of 0.9356 determined from the statistical analysis shows a good significance of the model. The overall utilities in terms of energy are found to be optimised by 81.4% from 6.511 x 107 kcal/h actual value of 1.211 x 107 kcal/h with 297.4 tonnes/day the carbon emission savings.

### Keywords

• central composite design
• optimisation
• response surface methodology
• cement production
• design expert

## 1. Introduction

Optimization is a mathematical technique used to find the best solution to objective function (s) by maximising the desired variable and minimising the undesired variables under some set of constraints with the sole aim of improving performance and cost [23]. The optimisation technique in cement and concrete studies can be broadly classified as a meta-heuristic approach and statistical experimental design methods [24]. The meta-heuristic approach is an iterative method that intelligently exploits search space at learning strategies. It includes Genetic Algorithm (GAs), Particle swarm optimization (PSO), Harmony Search (HS), Ant Colony Optimization (ACO), Charged System Search (CSS), Big Bang-Big Crunch (BB-BC), Artificial Bee Colony algorithm (ABC), spherical interpolation of the objective function, Colliding bodies optimization (CBO), Vibrating Particles System (VPS), simulated annealing, krill herd (KH), Whale Optimization Algorithm (WOA), hybrid Harmony Search, force method and genetic algorithm, mine and improved mine blast algorithms [24, 25, 26] which are modelled from natural and social behaviours as well as physics laws.

Statistical experimental design methods are widely used to obtain desired optimise solution for a set of constraints [27]. Response Surface Methodology (RSM) is a statistical optimisation technique used to model and analyse a process to determine the effect of independent multivariable on the process response and to evaluate the relations between these variables [28]. RSM is based on understanding the topography of the response surface to determine the most appropriate response region [29]. RSM experimental design can be categorised into Box–Behnken Design (BBD), Central Composite Design (CCD), Dohlert design, Mixture response and three-level factorial design [30, 31, 32]. The BBD is created from 3 level factorial design [32] and gives quadratic response model with three minimum number of factors requiring three levels of factors (upper, centre, lower) for each factor and specific positioning of design points [33, 34]. The CCD is developed from the 2 factorial design and gives the quadratic response model with five levels for each factor. Hence it is more robust and insensitive to missing data or experimental runs [34].

In recent years, Response Surface Methodology has been applied to optimise several chemical processes such as extraction [29], adsorption [23], pharmaceutical wastewater treatment [28], leaching [35]. Studies on cement production optimisation have been carried out on clinker simulation using AspenTech [36], cement raw materials blending using a general nonlinear time-varying model [37], cement grinding using population balance model [6], clinker chemistry and kiln energy efficiency using metaheuristic optimization techniques [38], numerical and computational fluid dynamics study of cement calciner [16]. RSM has been efficient and accurate in studies on cement and concrete technology [39, 40, 41, 42, 43]. This study focused on the simulation of an integrated wet cement process flow sheet using Aspen HYSYS and optimisation of the cement production rate at minimum raw material feed using CCD of response surface methodology.

## 2. Methodology

### 2.1 Cement production simulation

Aspen HYSYS 8.8 was used for the steady-state simulation of the integrated process flow sheet for the cement production [44]. Within the simulation environment, topological optimization (proper arrangement of equipment) was done to enable very high energy savings or optimization. A pure component such as water, CO2 and air are added as conventional components, while non-conventional components are added as hypothetical components to the HYSYS environment based on their physical properties (molecular weight & density). Figure 1 shows the block diagram of the production of cement while the HYSYS process flow diagram for the cement production simulation is shown in Figure A1 in the appendix. Limestone is decomposed in the first reactor to give off CO2 as gas, while the produced CaO is the feed to the Section A reactor to react with silicate to form Calcium disilicate. The produced Calcium disilicate reacts further with the unreacted CaO in reactor B to produce Calcium trisilicate. Calcium oxide (CaO) further reacts with Aluminium oxide to produce tricalcium aluminate, another constituent of cement, while the final product component is produced in section C, where CaO reacts with Aluminate and Ferric oxide to produce tetracalciumaluminoferrite. These separate components produced at a different section of the simulated Kiln are mixed to achieve a matrix compound of the cement product, having over 70% of CaO.

### 2.2 Aspen Hysys simulation

Aspen Hysys was used for the steady-state simulation of the integrated process flow sheet for the cement production. Within the simulation environment, topological optimization (proper arrangement of equipment) was done to enable very high energy savings or optimization. A pure component such as water, CO2 and air are added as conventional components, while non-conventional components are added as hypothetical components to the HYSYS environment based on their physical properties (molecular weight and density). Based on the process description, the different reactions taking place in each simulated reactor, as presented in the flowchart are:

CaCO3CaO+CO2Limestone decompositionE1
2CaO+SiO2Ca2SiO4SectionAE2
CaO+Ca2SiO4Ca3SiO5SectionBE3
3CaO+Al2O3Ca3Al2O6SectionCE4
4CaO+Al2O3+Fe2O3Ca4Al2Fe2O10SectionDE5

The various products in the various sections of the process reactors are; Tricalcium silicate (Ca2SiO4) which is responsible for early strength and the initial set of the cement; Dicalcium silicate (Ca3SiO5) which increases the strength as it age; Tricalcium aluminate (Ca3Al2O6) which contributes to the concrete strength development in the first few days but least desirable due to its reactiveness with sulphate containing soils and water; Tetracaliumaluminoferrite (Ca4Al2Fe2O10) which reduces clinkering temperature. The equipment design parameters employed in this work are provided in Table 1 .

Simple separatorDelta PStream fractions
0.000Solids in vapour0.0100
Solids in liquid0.0100
Liquid in bottoms0.0100
ExchangerDelta P (bar)Delta T (°C)Duty (kcal/h)
Cooler0.0006552.691 x0 107
Component SplitterTop Temp. (°C)Bottom Temp. (°C)Top Pressure (bar)Bottom Pressure (bar)Duty
Splitter 13030112.054 x 105
Splitter 212521252111.825 x 104
Splitter 312521252111.923 x 104
ReactorsDelta PVessel volume (m3)Liquid level (%)Liquid volume (m3)Duty (kcal/h)
CaCO3 decomp reactor0.00000.000050.0025.002.6764 x 107
Section A reactor0.000050.0050.0025.007.6209 x 106
Section B reactor0.000050.0050.0025.002.6131 x 106
Section C reactor0.000050.0050.0025.004.5634 x 105
Section D reactor0.000050.0050.0025.007.4721 x 105

### Table 1.

Equipment design parameter.

The flow rate of the major raw materials for the production of cement in the clinkering reactor as depicted by Eqns. (6–12) are carefully chosen based on the standard provided by Winter N. B. [45]. The Chemical parameters based on the oxide composition are very useful in describing clinker characteristics. The following parameters are widely used.

1. Lime Saturation Factor (LSF):is the measure of the ratio of alite to belite in the clinker. It is estimated by the ratio of CaO to the sum of other three main oxides SiO2, Fe2O3 and Al2O3. The equation is given by:

LSF=CaO2.8SiO2+1.2Al2O3+0.65Fe2O3E6
LSF=1902.860+1.215+0.6510=0.98E7

LSF values in clinkers range between 0.92–0.98. The LSF value of 0.98 falls within an acceptable range.

1. Silica Ratio (SR):This is also known as Silica Modulus. The expression of SR is given as:

SR=SiO2Al2O3+Fe2O3E8

Based on the experimental design for the simulated cement production process. The flow rate of Al2O3 and Fe2O3 are 15 tonnes/day and 10 tonnes/day respectively. The low level and high levels of SiO2 are found to be 50 tonnes/day and 60 tonnes/day respectively. Hence, the SR values are the high and low value of the SiO2 flow rate are calculated as follows:

SR=5015+10=2.0E9
SR=6015+10=2.4E10

A high silicate ratio means that more calcium silicates are present in the clinker and less aluminate and ferrite. SR is typical, between 2.0 and 3.0. The SR values of 2.0 and 2.4 fall within an acceptable range of 2.0 and 3.0.

1. Aluminate Ratio (AR):This is the ratio of aluminate and ferrite phases in the clinker. AR value ranges between 1–4 in Portland clinkers. The flow rate of Al2O3 and Fe2O3 used in the process simulation are 15 tonnes/day and 10 tonnes/day respectively. The equation governing the AR of the oxide is given by

AR=Al2O3Fe2O3E11
AR=1510=1.5E12

The mass flow and corresponding clinker quality parameters are presented in Table 2 .

CaOSiO2Al2O3Fe2O3LSFSRAR
High level1906015100.982.41.5
Low level1355015100.822.01.5

### Table 2.

Raw material mass flow and clinker quality parameter.

### 2.3 Multivariate design of experiment

The central composite design of response surface methodology was used to analyse the effect of CaO and SiO2 on cement production rate. The total number of experimental runs (N) required for nindependent variables and ncnumber of replica centre points is given by Eq. 13

N=2n+2n+ncE13

Design Expert 10.0.3 software was used to generate the experimental design from the ten experimental runs to study the combined effect of two variables on the response. For two variables factor in the experiments; four factorial points (2n), four axial points (2n) and two replicates at the central points (nc) at distance α = 1.414 from the centre were used for the CCD design. A polynomial empirical model was developed from the ten experimental runs to correlate the response with the independent variables. The mathematical expression can be expressed as:

Y=βo+i=1nβiXi+i=1nβiiXi2+1ijnβijXiXj+εE14

Where Y= Predicted response, βo= constant coefficient, βi= Linear coefficient, βii= Quadratic, βij= Interaction coefficients and ε = model random error, nthe number of variable factors, Xiand Xjare the coded values of the variable parameters [35].

The response generated function distance from the centre α=2n/4. The codes are calculated as a function of the range the factors as shown in Table 3 .

CodeMathematical relationship
Lower axial point- αXmin
Lower level- 1Xmax+Xmin/2αXmaxXmin/2
Centre point0Xmax+Xmin/2
Upper level1Xmax+Xmin/2+αXmaxXmin/2
Upper axial pointαXmax

### Table 3.

Relationship between the variable values and their assigned codes.

The central composite experimental design for the synthesis of cement via simulation is depicted in Table 4 . The mass flow rate of CaO and SiO2 measured in tonnes/day are the independent variables or predictors which are studied for their effect on the response variable (cement flow rate) at a constant Al2O3 and Fe2O3 flow rates.

FactorNameUnitsLowHigh
AFlow rate of CaOtonnes/day135190123.609201.391
BFlow rate of SiO2tonnes/day506047.928962.0711

### Table 4.

Design of Experiment using central composite (C.C.D) Design of Response Surface Methodology (R.S.M).

### 2.4 Model fitting and statistical analysis

The interaction between the variables and the response data as well as the statistical parameters were analysed graphically by analysis of variance (ANOVA) in the Design-Expert software. Regression analysis, significance, F-test, surface and contour plots of the response were also generated from the software. A probability value of 95% confidence level was used to evaluate the significance of the model terms and coefficients.

## 3. Results and discussion

### 3.1 Simulation and optimisation of cement flow

All simulations were done in duplicate and the experimental design were generated by the Central Composite Design (C.C.D) of the Design-Expert Software, which resulted in a total of 10 experimental (simulation) runs and the results of the experiments (simulations) are shown in Table 5 .

Run No.FactorsResponse
Coded levelsActual valuesThe flow rate of cement
ABABHYSYS SimulationC.C.D Model
11.0001.0001906046.645.6674
2−1.000−1.0001355030.7431.9851
30.0000.000162.55546.0446.04
40.000−1.414162.547.928930.5731.271
51.000−1.0001905044.342.1551
6−1.4140.000123.6095546.0443.4867
70.0000.000162.55546.0446.04
8−1.0001.0001356046.649.0574
91.4140.000201.3915546.0448.2808
100.0001.414162.562.071146.8445.8265

### Table 5.

Simulation and predicted results from central composite design (C.C.D).

The change in mean response in cement flow per unit increase in variable occurs when other predictors area kept constant and is estimated by the coefficient of estimation and is presented Table 6 .

CoefficientStandard95% CI95% CI
FactorEstimatedfErrorLowHighVIF
Intercept46.0410.8644.0048.08
A-Flowrate of CaO1.7010.680.0803.311.00
B-Flow rate of Silica5.1510.683.536.761.00
AB−3.3910.97−5.67−1.111.00
A2−0.07810.73−1.811.651.02
B2−3.7510.73−5.48−2.011.02

### Table 6.

Coefficient estimation for cement flow rate in terms of coded factors.

The empirical quadratic equation for the optimal cement product rate as a function of CaO and SiO2 mass flow in coded form as derived from Table 7 was obtained according to the CCD and is given in Eq. 15

SourceSum of squaresdfMeanF valuep-value (Prob > F)
Model379.61575.9220.350.0005
A22.98122.986.160.0421
B211.861211.8656.780.0001
AB45.97145.9712.320.0099
A20.04210.0420.0110.9180
B297.60197.6026.160.0014
Residual26.1273.73
Lack of Fit26.1238.71
Pure Error0.00040.000
Cor. Total405.7312

### Table 7.

ANOVA results for the statistical model for cement flow rate.

C=46.04+1.70A+5.15B3.39AB0.078A23.75B2E15

The test of the significance and the adequacy of the model and its coefficients lack fitness which was based on F-value or P-value at 95% confidence level was tested from analysis of variance (ANOVA) and the result is presented in Table 7 . The result shows that the model at an F value of 20.35 and a very low P-value of 0.0005 indicates that the statistical regression model was significant. The result also shows that A, B, ABand B2are significant terms.

The regression statistical analysis is summarised in Table 8 . The R squared value of 0.9356 is in good agreement with the adjusted R-square value of 0.8896, showing a good fit of the model, as the closer the R squared value to 1.00, the more significant the model. The adequate precision of 13.55 indicates the low noise level and a strong signal for optimisation. Hence, this indicates the two predictors (flow rate of CaO and flow rate of SiO2) could predict the flow rate of cement, thus the model equation, contour plot and 3D surface plot could be used to predict the response (flow rate of the cement).

ParameterValues
R-Squared0.9356
Pred R-Squared0.5422
Mean43.69
Std. Dev.1.93
C.V. %4.42

### Table 8.

Statistical information for the statistical model for cement flow rate.

### 3.2 Surface plots

The contour plot which shows the possible relationship between the CaO, SiO2 and cement product mass flow is presented in Figure 2 . The darker red regions indicate higher C (response) values. Here, the optimum flow rate of cement is found from the isolines to be 47.748 tonnes/day at a flow rate of 152.346 tonnes/day of CaO and 56.8241 tonnes/day SiO2.

The three-dimensional (3D) response surface plots obtained from the model equation using Design Expert 10.03 is depicted in Figure 3 . This depicts the effect of the flow rate of CaO, the flow rate of SiO2 on the flow rate of cement. The flow rate of cement was observed to increase with an increase in the flow rate of CaO. Conversely, increasing the flow rate of SiO2 did not increase the flow rate of cement. Hence, the major predictor in cement production in the clinkering section is the flow rate of CaO.

The optimization plot of the cement output is shown in Figure 4 . The optimum cement flow rate of 47.239 tonnes/day is found to be at CaO flow rate of 152.346 tonnes/day and SiO2 flow rate of 56.8241 tonnes/day.

### 3.3 Energy optimization

Aspen Energy Analyser was used to determine the percentage of energy savings based on converged steady-state simulation of the process flow sheet in Figure 1 . The total energy savings as a function of process utilities and carbon emissions are present in Figure 5 .

The overall utilities in terms of energy are found to be optimised from the actual value of 6.511 x 107 kcal/h to 1.211 x 107 kcal/h and indicating available energy savings of 5.3 x 107 kcal/h, with overall energy savings of 81.40% which also correspond to 297.4 tonnes/day carbon emission reduction.

## 4. Conclusion

Process flow diagram for the cement production was simulated to achieve high energy optimization and optimum cement flow rate by minimising the flow rate of the feed (CaO and SiO2). Central composite Design (C.C.D) of Response Surface Methodology used to design the experiment for the simulation using Design Expert 10.0.3. The optimum cement flow rate is found from surface and contour plots to be 47.239 tonnes/day at CaO flow rate of 152.346 tonnes/day and SiO2 flow rate of 56.8241 tonnes/day. The R squared value of 0.9356 determined from the statistical analysis shows a very high significance of the model. Energy efficiency optimization is also carried out using Aspen Energy Analyser. The overall utilities in terms of energy are found to be optimised by 81.4 % from 6.511 x 107 kcal/hactual value to 1.211 x 107 kcal/h with 297.4 tonnes/day the carbon emission savings.

Further work could be performed on fault identification and diagnosis of the process plant. Incorporated with an automated plant to guarantee the safety of workers, reduce environmental problems and increase yield to sustain production improvement.

## 5. Recommendations

This research work sought to recommend the following concerns in which cement production could be improved:

1. Research and development (R&D) in the cement production, processing and utilisation should be encouraged. This will play a vital role in the construction industry, operation and maintenance of efficient road network and effective transportation system.

2. Automation of cement and kiln sections of the cement production is recommended

3. Optimization of the cement production can be tailored into the fabrication of high tech cement equipment and parts.

4. Optimization of the limestone crusher to quantify the amount of crushed limestone is needed.

5. Looking into future the results obtained in this research will open up several important possibilities in the cement production at optimum conditions. This will have a multiplier effect on infrastructural amenities development.

## Conflict of interest

There is no conflict of interest associated with this work.

chapter PDF

## More

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Oluwafemi M. Fadayini, Adekunle A. Obisanya, Gloria O. Ajiboye, Clement Madu, Tajudeen O. Ipaye, Taiwo O. Rabiu, Shola J. Ajayi and Joseph T. Akintola (March 9th 2021). Simulation and Optimization of an Integrated Process Flow Sheet for Cement Production [Online First], IntechOpen, DOI: 10.5772/intechopen.95269. Available from: