Next Article in Journal
Formulation and Characterization of Double Emulsions W/O/W Stabilized by Two Natural Polymers with Two Manufacturing Processes (Comparative Study)
Next Article in Special Issue
Robust Fault Detection in Monitoring Chemical Processes Using Multi-Scale PCA with KD Approach
Previous Article in Journal
Investigation on UV Degradation and Mechanism of 6:2 Fluorotelomer Sulfonamide Alkyl Betaine, Based on Model Compound Perfluorooctanoic Acid
Previous Article in Special Issue
Numerical Simulation of a Valorisation-Oriented Hybrid Process for the Bio-Oil-Related Separation of Acetol and Acetic Acid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Industrial Control System for Cement Sulfates Content Using a Feedforward and Feedback Mechanism

by
Dimitris Tsamatsoulis
Heidelberg Materials Hellas S.A., Heidelberg Materials Group, 17th Km Nat. Rd. Athens–Korinthos, 19300 Aspropyrgos, Greece
ChemEngineering 2024, 8(2), 33; https://doi.org/10.3390/chemengineering8020033
Submission received: 23 November 2023 / Revised: 13 February 2024 / Accepted: 1 March 2024 / Published: 7 March 2024
(This article belongs to the Special Issue Feature Papers in Chemical Engineering)

Abstract

:
This study examines the design and long-term implementation of a feedforward and feedback (FF–FB) mechanism in a control system for cement sulfates applied to all types of cement produced in two mills at a production facility. We compared the results with those of a previous controller (SC) that operated in the same unit. The Shewhart charts of the annual SO3 mean values and the nonparametric Mann–Whitney test demonstrate that, for the FF–FB controller, the mean values more effectively approach the SO3 target than the older controller in two out of the three cement types. The s-charts for the annual standard deviation of all cement types and mills indicate that the ratio of the central lines of FF–FB to SC ranges from 0.39 to 0.59, representing a significant improvement. The application of the error propagation technique validates and explains these improvements. The effectiveness of the installed system is due to two main factors. The feedforward (FF) component tracks the set point of SO3 when the mill begins grinding a different type of cement, while the feedback (FB) component effectively attenuates the fluctuations in the sulfates of the raw materials.

1. Introduction

There is widespread agreement that the sulfate (SO3) content of cement is a crucial quality parameter because it affects the compressive strength, setting time, and long-term performance of cement. For this reason, cement standards [1,2] stipulate that cement must contain clinker and calcium sulfate (Cs) and define a high SO3 limit for each product type. Gypsum is the primary form of calcium sulfate and is accurately fed into the cement mill (CM) during the grinding process.
One of the essential mineral phases of clinker is the tricalcium aluminate (3CaO ∙ Al2O3 or C3A), which reacts very fast with water (H). Gypsum addition retards the fast hydration of C3A by generating ettringite (C3A.3Cs.32H) [3], according to Equation (1).
C 3 A + 3 C s + 32 H C 3 A .3 C s .32 H
In concrete production, transfer, and placement, the formation of ettringite prevents flash setting caused by rapid C3A hydration. Conversely, an excess of gypsum leads to harmful expansion and a decrease in the strength of concrete and mortar [3,4,5]. Therefore, there is an optimal value for the sulfates. To the best of the author’s knowledge, numerous researchers have investigated the ideal SO3 level over the past 80 years due to its significance, researching the effects of sulfates on several important cement properties. Lerch [6], a pioneer in cement research, conducted the first in-depth study on sulfate optimization. Several researchers [7,8,9,10,11,12,13] have investigated and identified the modifications in the hydration rate of clinker mineral phases after adding different amounts of sulfates. These modifications affect fundamental cement properties, such as water demand for normal consistency, setting times, compressive strength at various ages, heat of hydration, and hydration degree.
The cement industry has prioritized reducing its carbon footprint in recent years by decreasing CO2 emissions in clinker production and clinker consumption per ton of product. The latter is achievable by incorporating supplementary cementitious materials (SCM) into the cement composition. Optimizing SO3 can reduce the incorporation of clinker into cement while maintaining or improving product performance. The best sulfate content in systems containing clinker and one or more SCMs has been extensively investigated [4,14,15,16,17,18,19,20]. The optimal position depends on the property that needs optimization. Niemuth [4] studied the impact of incorporating fly ash into Portland cement on the optimum sulfate content. At various SO3 levels, he presented experimental data on strength development and heat release during early hydration. Adu-Amankwah et al. [14] investigated the effects of sulfate additions on the hydration and performance of ternary slag–limestone composite cement through porosity and strength measurements. Han et al. (2015) examined the influence of gypsum on the characteristics of composite binders containing slag and iron tailing powder using a range of measurement techniques. Yamashita et al. [16] studied the not negligible impact of limestone powder on the optimal SO3 for Portland cement with varying Al2O3 content, using compressive strength as a criterion. Liu et al. [17] examined the effect of gypsum content on cementitious mixtures containing limestone, fly ash, and slag by studying various properties, including initial and final setting time, paste fluidity, water demand, and strength. Fiscan [18] studied the optimal sulfates in cement–slag blends using calorimetry and early strength results. Tsamatsoulis et al. [19] attempted to determine the SO3 optimum of Portland, Portland composite, and pozzolanic cement types by implementing a unified approach and shallow artificial neural networks. Andrade Neto et al. [20] compiled laboratory techniques for estimating the optimal sulfate content and described the benefits and drawbacks of each method.
Simply knowing the optimal level of sulfates for each type of cement and setting it as a target in daily production is insufficient for cement manufacturing. The measured SO3 levels should closely align with this target with minimal variance. Therefore, continuous regulation of gypsum, particularly with a controller, is essential to achieve this goal. To the best of the author’s knowledge, it is hard to find a description of cement sulfate controllers installed in milling systems in the literature. In the author’s experience, most cement plants use manual step-change rules for SO3 regulation. In a previous study [21], we developed simulations to compare the results of manual regulation with those of a controller comprising both feedback (FB) and feedforward (FF) parts. The FB component attenuates process disturbances, whereas the FF tracks changes at the set point (SP). Combining these two independent regulators has several practical applications. Ko et al. [22] analyzed an FF–FB regulator for an electro-hydraulic valve system utilizing a proportional control valve. The FB component was a proportional-integral-derivative (PID) controller. Wang et al. [23] designed a composite control model containing FF and FB controllers for optical fiber alignment using a piezoelectric actuator. Araque et al. [24] implemented the same technique by combining the two control types for temperature uniformity control. The authors demonstrated that incorporating a model-based feedforward loop improves the tracking of reference signals.
This study analyzes the design and implementation of an FF–FB system to control the SO3 content in the cement mill outlet by adjusting the percentage of gypsum in the CM inlet. We applied this control technique to two CMs of the Halyps plant for the cement types (CEM) produced. The simulation presented in [21] used the same cement mills. The main novelty of this study is the design and long-term implementation of such a system in cement manufacturing, as it is difficult to find a description of such controllers installed in milling systems in the literature. The structure of the paper is as follows. Section 2 provides a brief description of the grinding process, the types of CEM used, and raw materials analyses. Subsequently, we present the design of the FF–FB control system and its digital implementation in the quality control of sulfates during cement production. The author developed all the software in C# 9.0. Additionally, we briefly describe the rules previously applied to adjust SO3 in the same milling facilities. We conclude Section 2 by comparing the two control techniques using a process simulator. Section 3 analyzes the long-term results of the controller by comparing them with the results of previous SO3 adjustments applied to the same installations. We conducted the assessment based on industrial data from 19 consecutive years, covering the period from 2005 to 2023. Finally, Section 4 summarizes the primary findings of this industrially applied research.

2. Process Description and Control Technique

2.1. Process Description and Materials Analysis

Cement plants typically grind cement in closed milling systems. Figure 1 shows a simplified flowchart of a grinding circuit, including all essential installations. We used the same configuration in [25] in a study of optimization of the process control of cement milling. The weight feeders feed the raw materials to either the ball mill or the separator (fly ash). The recycling elevator directs the output from the mill to the dynamic separator. The fine stream from the classifier constitutes the final product, whereas the coarse material returns to the CM for further grinding. The critical parameters related to quality include (i) cement fineness, (ii) separator speed, (iii) ratio of the coarse material-flow rate to the mill-feed rate, (iv) recycling elevator power, and (v) air-flow rate through the mill and pipes.
The plant quality department regulates the sulfates by sampling the cement in the mill outlet, measuring its SO3 content, and adjusting the gypsum proportion in the CM feed. The control was applied to five CEM types produced according to EN 197-1:2011 [1] and is shown in Table 1. We presented the same Table in [19], where we optimized the sulfate content of the same CEM types. The range of the products is broad, covering Portland (I, II) and pozzolanic (IV) types as well as all three strength classes (32.5, 42.5, and 52.5). Table 2 presents the long-term statistics of the SO3 content in the raw materials. The lab conducted analyses of Portland CEM types on Oxford Instruments (Oxfordshire, UK) LAB X 3000 (2005–2015) and Hitachi (Tokyo, Japan) X-Supreme 8000 (2016–2023) XRF analyzers, while Malvern-Panalytical (Almelo, The Netherlands) Axios-Cement and Zetium carried out the analyses of the raw materials and pozzolanic CEM types.
Both clinker and fly ash contain significant amounts of SO3, with noticeable variations. The coefficient of variation %CV (=Std. Dev./Aver. × 100) lies within the range of 35.5% to 44.8%. The above causes two types of disturbance. (a) When the CEM type changes, the running composition can lead to a low-frequency step disturbance in the flow rate of sulfates due to differing percentages of clinker and (or) fly ash. Therefore, it is likely necessary to adjust the gypsum to achieve the current target. (b) The variation in SO3 content of the two mentioned materials causes disturbances during operation with the same CEM type, necessitating attenuation by adjusting the proportion of gypsum.

2.2. Controller Design

The presence of two distinct types of disturbances in the process variable provides the benefit of employing two controllers, acting separately on the control variable. Our earlier study introduced a dual regulator [21], comprising feedforward and feedback modules (FF–FB). Figure 2 depicts the block diagram of the transfer functions and signals related to the sulfate control.
The signal SP represents the SO3 target of the current CEM, and DSP is the signal for the SO3 target in case the mill starts to grind another CEM type. DSP is a low-frequency disturbance in the control loop that takes nonzero values only when there is a change in the cement type. XFB and XFF are signals expressing gypsum percentages derived from the FB and FF controllers. Dp is the SO3 disturbance inserted into the process through the ingredient feeders (clinker, gypsum, and fly ash) due to the variance in the raw-materials composition. Signal Y represents the sulfate content of the product exiting the closed grinding circuit. Si expresses the SO3 percentage of the cement after sampling and measurement. Figure 2 shows the following transfer functions. GP refers to the gypsum mixing within the milling system. GM is a time delay function for sampling and sulfate measurement. GFB and GFF indicate the FB and FF controllers, respectively. If DSP = 0, then AFB = 1 and AFF = 0. On the contrary, if DSP ≠ 0, then AFB = 0 and AFF = 1.
Equations (2) and (3) provide, in the Laplace domain, the open loop transfer function GOL and the transfer function from SP and disturbances to the output Y. GOL represents the function of the system where the SO3 output is not fed back to control the gypsum percentage.
G O L = G P · G M · A F B · G F B + A F F · G F F  
Y = G P · A F B · G F B + A F F · G F F 1 + G O L · S P + G P · G F F 1 + G O L · D S P + G P 1 + G O L · D P
Equation (4) expresses the GM function. The average sampling and measuring time, TM, which is a pure delay, is 0.25 h. As shown in reference [21], the transfer function GP between the gypsum percentage in the CM feed and the %SO3 in the final product can be modeled using first-order dynamics with time delay (FOTD). Equation (5) describes the model, where TD is the delay time, T0 is the time constant, and KV is the gain.
G M = e T M · s  
G P = K V · e T D · s 1 + T 0 · s  
According to [21], the milling circuit’s dynamic parameters are KV = 0.4, TD = 0.133 h, and T0 = 0.233 h. The gain meaning is the increase in %SO3 for a 1% increase in gypsum dosage. We conclude that the gain value is near the respective value computed from Table 2 (43.67/100 = 0.44). The meaning of TD and T0 is that after TD, a step increase in gypsum will affect SO3 in the CM outlet. After T0, the SO3 change is 63% of the total increase. Equation (6) provides the GP transient response in the time domain after a ΔG step change in %gypsum.
S O 3 t S O 3 0 G = K V 1 e t T D / T 0  
SO3(0) and SO3(t) are the %SO3 values at the beginning of the step increase and at time t, respectively. In the steady state, the maximum %SO3 increase is ΔSO3 = KV·ΔG. Equation (7) computes the fraction a(t) of ΔSO3 at time t.
a t = 1 e t T D / T 0  
Assuming that the system is near the steady state when the fraction α(t) reaches the value of 0.98, the required time calculated from Equation (7) is TTP = 1.04 h. After this transient period, the system is in equilibrium with respect to step changes in gypsum feeding. For spot sampling, the delay time between the next and previous feedback controller outlets is the sum of TTP and TM, which is equal to Ts,Min = 1.29 h. Ts,Min is the minimum sampling period to avoid transient phenomena.
A simple integral controller (I) appropriately regulates the feedback control loop [21] with gain ki. Similarly, a proportional controller of gain KFF attenuates the low-frequency disturbances of the feedforward loop. Equation (8) provides the respective transfer functions.
G F B = k i s   ;   G F F = K F F  

2.3. Digital Implementation

The feedback element of the controller calculates the gypsum setting of the CM feeder after each SO3 measurement of an instantaneous (spot) sample, which is performed at regular time intervals. In contrast, the feedforward element acts only when the CM starts producing a different CEM type. Consequently, the controller operates in discrete time intervals characterized by the sampling period Ts. Equation (9) computes the error ei between the SO3 set point, SSP, and the SO3 of the sample taken at time i, Si.
e i = S S P S i  
The discrete implementation of the feedback integral controller utilizes the backward form [26] to calculate the gypsum percentage Gi at time i by adding the control action to the gypsum content Gi−1 of time i − 1. Equation (10) expresses this function. The sampling period is Ts = 2 h, and the optimal gain ki is 0.8, as found in an earlier simulation study [21]. The controller is unconstrained. The set of Equations (11)–(15) implements the feedforward proportional controller.
G i = G i 1 + k i · e i · T s  
d S P = C l C E M , N C l C E M , P · S C l + A s h C E M , N A s h C E M , P · S A s h 100  
G P r e v = G i   ;   S P = S i   ;   S N = S P + d S P   ;   S i = S N  
K F F = 100 S G   ;   D G = K F F · S S P S i  
I f   A B S D G > M a r g :   i f   D G < 0   : D G = M a r g   e l s e i f   D G > 0   D G = M a r g  
G i = G P r e v + D G  
SP is the SO3 measured at time i, GPrev is the feedback controller output at time i, dSP is the disturbance due to the CEM type change, and SN is the SO3 content considering the disturbance. ClCEM,P, ClCEM,N, AshCEM,P, and AshCEM,N are the average clinker and fly-ash contents of the previous and current CEM types, respectively. SCl, SAsh, and SG are the mean SO3 contents of clinker, fly ash, and gypsum, as shown in Table 2. Equation (13) calculates the unconstrained output of the controller, DG. However, the optimal feedforward controller is constrained, as proven in [21]. The conditions in (14) implement the constraints for the maximum absolute change of the DG. Equation (16) calculates the value of the margin Marg, which depends on the previous and current CEM types.
M a r g = A B S S S P , N S S P , P d S P · K F F + M 0  
SSP,P and SSP,N are the SO3 targets of the previous and current CEM, respectively, and M0 is an additional margin of gypsum. In our application, M0 = 0.5.

2.4. Comparisons Using a Process Simulator

For the past eleven years, the FF–FB controller designed in Section 2.2 has been operating in CM5 and CM6 of the Halyps plant to regulate the sulfates for all the CEM types produced. Before using this regulator, the plant employed the step rules (SC) controller, as mentioned in [21] and shown in Equation (17). This regulator consists of a dead band of 0.4 for SO3 and provides step changes in the gypsum feed with a gain of 0.5 or a multiple thereof. The description of Equations (11)–(15) provides the physical meaning of the parameters of Equation (17).
e i   ϵ   0.2   ,   0.2   D G = 0 e i   ϵ   0.2   ,   0.6   D G = 0.5 e i   ϵ   0.6   ,   1.0   D G = 1.0 e i > 1.0   D G = 1.5 e i   ϵ   0.6   , 0.2   D G = 0.5 e i   ϵ   1.0   , 0.6   D G = 1.0 e i < 1.0   D G = 1.5 G i = G i 1 + D G  
Process simulators allow comparisons between SC and FF–FB controllers. The simulator runs the CM for 600 h. The cement type changes every 20 h between CEM II B-M (P-L) 32.5 and CEM II A-L 42.5. A disturbance occurs in clinker SO3 every 10 h, i.e., two disturbances appear every 20 h. Clinker SO3 never changes when the CEM type changes. The simulator creates the magnitude of each disturbance using a predefined mean and standard deviation of clinker SO3, a randomly generated probability, and the inverse normal distribution. The simulator also calculates a low variation in gypsum sulfate every 2 h. Table 3 lists the parameters of the developed simulation. Figure 3 shows the SO3 results for the SC and FF–FB controllers for the CEM II B-M (P-L) 32.5. The simulator used the same disturbances for both control techniques. Compared to SC, the FF–FB results have less dispersion around the target.
A quantitative comparison of the two controllers’ performances is feasible by implementing the simulator multiple times. The simulator used the same disturbances for each run. Table 4 presents the average statistical results after 100 implementations.
We used three criteria to evaluate the closeness to SSP: (a) the average standard deviation, (b) the percentage of the population out of the interval [0.95·SSP, 1.05·SSP], and (c) the same statistic as (b) but for the first SO3 values when the CEM type changes. Criterion (c) assesses the efficiency in set-point tracking, whereas criteria (a) and (b) evaluate the degree of disturbance rejection. The three statistics’ ratios are consistently smaller than one, indicating that the FF–FB system outperforms the SC in disturbance attenuation and set-point tracking.

3. Long-Term Results and Analysis

3.1. Shewhart Control Charts and Nonparametric Analysis

The Shewhart control charts [27] (pp. 8–9) are suitable for comparing the results of the two control techniques in the long term. These charts require data in rational subgroups taken at approximately regular intervals during the process. In this study, each subgroup contains all the SO3 results of samples taken during one year per CM and CEM type. We performed the comparison by generating mean ( X ¯ ) and standard deviation (s) charts. We separated the results into two groups: (a) the period of SO3 adjustment using Equation (16) (2005–2012) and (b) the period of FF–FB controller application (2013–2023). The central line (CL) of the X ¯ -chart is a prespecified process parameter equal to the SO3 target per CEM type. Because the number of samples ni varies annually, we used the pooled standard deviation sPk [28] (p. 93) to determine the central line of the s-chart and the control limits, as shown in Equation (18).
s P k = i = 1 M n i k 1 · s i k 2 i = 1 M k n i k M k 1 2   i = 1 , 2 M k   k = 1 :   S C   a n d   k = 2 : F F / F B  
where Mk is the number of years in the selected period and sik is the SO3 standard deviation of year i and period k. The annual number of samples is sufficient to calculate statistics when nik ≥ 20.
There is a maximum size of 25 samples in each subgroup in Table 2 of ISO 7870-2 [27] (p. 9), which includes the factors to compute the lower and upper control limits, LCL and U C L . This table is not suitable because all nik exceed this value. Reference [29] provides the general formulae to determine the control limits LCL and U C L , applicable to any number N of samples and given by Equations (19)–(21).
X ¯ c h a r t :   L C L = C L 3 · s P k c 4 N   U C L = C L + 3 · s P k c 4 N  
s c h a r t :   L C L = s P k 3 · s P k · 1 c 4 2 c 4   U C L = s P k + 3 · s P k · 1 c 4 2 c 4  
c 4 = 2 N 1 · N 2 1 ! N 1 2 1 !  
where the coefficient c4 uses the noninteger factorial, determined by the Gamma function and its properties: Γ(x + 1) = x·Γ(x) and Γ(½) = π1/2.
Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 illustrate the two kinds of control charts for the CEM types produced in both periods under examination. For each period, the number N needed for determining the upper and lower control limits is the average of the populations of yearly samples.
The mean charts depict the closeness of the annual mean SO3 to the target. In all CEM types and both CMs, applying the FF–FB controller leads these two variables into proximity. In the case of the SC regulator, the average sulfates and their target are close for the Portland types but not for the pozzolanic cement. The comparison of the differences between the target, STik, and the realized SO3, Savik, requires a statistical test. Equation (22) computes the absolute value dik of this difference, and Figure 9a,b shows this function for two CEM types.
d i k = S T i k S a v i k   i = 1 .. M k   a n d   k = 1 , 2  
Due to the asymmetry of all distributions, statistical tests based on the normal distribution are not applicable. Therefore, estimating the difference of means ([30], pp. 18–19) is inapplicable, and the test shall be nonparametric. Mann and Whitney [31] developed such a statistic, which the relative literature [32,33,34] continuously refers to, finding implementation in several research fields [35,36,37,38]. This test concerns the sets D1 and D2 with populations M1 and M2 described in (23).
D 1 = d 11 , d 21 , d M 1 1   D 1 = d 12 , d 22 , d M 2 2  
The method considers the set D = D1D2 and finds the rank of each element within the union D after sorting them in increasing order. R1 and R2 are the sum of the ranks for observations one and two. Equation (24) provides the test statistic U , which shall be compared with the critical values U c r .
U 1 = M 1 M 2 + M 1 M 1 + 1 2 R 1     U 2 = M 1 M 2 + M 2 M 2 + 1 2 R 2   U          = min U 1 ,   U 2
The null hypothesis H0 is that there is no tendency for the ranks of D1 to be significantly higher than that of D2 occurring when U > U c r . The alternative hypothesis HA is that the ranks of D1 are systematically higher than that of D2 occurring when U U c r . The hypotheses show that the test is one tail, and reference [39] provides the critical values for probabilities a = 0.01 and a = 0.05. Table 5 presents the results of the Mann–Whitney test for sets D1 and D2, leading to the following conclusions. The two controllers have equivalent performance concerning the CEM II B-M (P-L) 32.5. On the contrary, in CEM II A-L 42.5, the annual mean values better approximate the target using the FF–FB controller than applying the SC, with a probability of 95%. This improvement is most noticeable at CEM B (P-W) 32.5, where the values of the D1 set are higher than those of D2 with a probability of 99%.
The s-charts shown in Figure 4b, Figure 5b, Figure 6b, Figure 7b and Figure 8b illustrate a pronounced decrease in the annual standard deviation after the FF–FB controller started operating. The U C L of FF–FB is always lower than the LCL of SC, indicating that the former controller attenuates disturbances caused by sulfate variability in the raw materials better than the latter. Table 6 shows the pooled standard deviations for the five CEM types produced and their ratios when a CEM covers both periods (SC and FF–FB).
The severe and systematic improvement of SO3 stability using the FF–FB controller is apparent in this table. The ratio of sP2sP1 ranges from 0.51 to 0.59 for Portland CEM types, where the regulating action is the attenuation of clinker variability in SO3. The improvement is better in pozzolanic cement CEM IV B (P-W) 32.5, where SO3 disturbances originate from fly ash and clinker. The two pozzolanic types show similar standard deviations. Despite its high clinker content, CEM I 52.5 shows a lower pooled standard deviation during FF–FB operation than the respective values of all CEM types during SC operation.

3.2. Assessing Controllers’ Quality by Combining Standard Uncertainties

Cement SO3 is the sum of the sulfates in the clinker, fly ash, and gypsum. Consequently, accounting for error propagation, the variance and covariance of input variables affect cement SO3 variability. We assume that between an output y and input variables, there is a functional relationship y = f(x1, x2, … xN). Then, Equation (25) provides the square of the combined standard uncertainty u c 2 y as a function of the uncertainties u(xi) [28] (pp. 18–23).
u c 2 y = i = 1 N f x i 2 · u 2 x i + 2 · i = 1 N 1 j = i + 1 N r i , j · f x i · f x j · u x i · u x j  
where r(i, j) is the correlation coefficient between xi and xj. r(i, j) = 0 for uncorrelated variables, and it is positive or negative for positively or negatively correlated variables.
Equation (26) illustrates the sulfate mass balance of cement, where CL, G, and FA denote the fractions of clinker, gypsum, and fly ash in the cement composition, SCL, SG, and SFA are the sulfates coming from the raw materials, and SO3,CEM, SO3,CL, SO3,G, and SO3,FA are the SO3 percentages in cement and in the three materials.
S O 3 , C E M = S C L + S G + S F A   S C L = C L · S O 3 , C L   S G = G · S O 3 , G   S F A = F A · S O 3 , F A  
The pairs of variables (CL, SO3,CL), (G, SO3,G), and (FA, SO3,FA) are uncorrelated. Therefore, Equations (27)–(29) give the uncertainties uS,CL, uS,G, and uS,FA of SO3 of each component within the CEM composition.
u S , C L 2 = C L 2 · u S O 3 , C L 2 + u C L 2 · S O 3 , C L 2  
u S , G 2 = G 2 · u S O 3 , G 2 + u G 2 · S O 3 , G 2  
u S , F A 2 = F A 2 · u S O 3 , F A 2 + u F A 2 · S O 3 , F A 2  
where uCL, uG, uFA are the uncertainties of the fractions CL, G, and FA in the CM feeders, SO3,CL, SO3,G, SO3,FA are the sulfates of CL, G, and FA, and uSO3,CL, uSO3,G, uSO3,FA are the respective uncertainties. If a controller regulates cement SO3 by changing the gypsum, the variables SCL and SG are negatively correlated. The same occurs for SFA and SG. Equation (30) specifies this case for the uncertainty of cement sulfates uS,CEM.
u S , C E M 2 = u S ,   C L 2 + u S , G 2 + u S , F A 2 2 r C L , G · u S , C L · u S , G 2 r F A , G · u S , F A · u S , G  
where the correlation coefficients of CL and G, r(CL, G), and of FA and G, r(FA, G) belong to the interval [0, 1]. The model of Equations (26)–(30) requires the estimation of correlation coefficients based on actual quality data of each CEM type and raw material. Table 7 demonstrates the type of data used for all the variables.
In the case of cement containing up to four components (clinker, gypsum, limestone, and pozzolan), we determined the daily fractions of clinker and gypsum in the composition by using SO3, loss on ignition, and insoluble residue [40] of cement and raw materials and solving the respective linear system. In the case of cement with fly ash (CEM IV B (P-V) 32.5), we also utilized the oxides CaO, SiO2, and Al2O3 to estimate the fly-ash content using the generalized reduced gradient nonlinear regression technique. We applied the same nonlinear method to calculate the correlation coefficients. A strong relationship exists between these coefficients and the ability of a controller to reject or attenuate disturbances. The larger the correlation coefficient between gypsum and clinker or fly ash, the more robustly the regulator adjusts the gypsum percentage to compensate for clinker or fly-ash sulfate disturbances or changes in cement composition. Table 8 shows the regression analysis of the standard deviation results after each controller was applied. We provide Table S1 with all the lab data used in the error propagation analysis for the two applied control techniques and all CEM types.
The correlation coefficient between clinker and gypsum, r(CL, G), using the FF–FB is significantly higher than that using the SC controller. The above perfectly explains the considerable improvement in the annual standard deviation in Portland cement, shown in Table 8. Applying the SC controller, the r(FA, G) is around null, becoming significant using the FF–FB. The conclusion is that FA and G are uncorrelated, and the SC cannot compensate for SO3 disturbances in fly ash, which the FF–FB controller satisfactorily achieves, resulting in a noticeable drop in the pooled standard deviation for the pozzolanic cement. Figure 10a,b compares actual and calculated values, yact and ycalc, whereas Equations (31) and (32) provide these two variables. This figure clearly explains the significant increase in R2 in the FF–FB case compared with the R2 of SC.
y a c t = u s , C E M 2 u s ,   C L 2 u s , G 2 u s , F A 2  
y c a l c = 2 r C L , G · u s , C L · u s , G 2 r F A , G · u s , F A · u s , G  

4. Conclusions

This study analyzes the design of a feedforward and feedback mechanism and its implementation in an industrial control system for cement sulfates. The controller considers all the fundamental aspects and particularities of the grinding process and quality requirements. (a) Variability of the raw materials SO3; (b) CM dynamics; (c) sampling period and measuring delays; (d) cement composition and feeders’ accuracy; and (e) grinding of various CEM types with different sulfate targets. The results of the FF–FB controller, long-term applied in two CMs of the Halyps plant for all CEM types, have been compared with those of the previously applied controller, which used step rules to adjust gypsum. The main conclusions of this study are as follows.
(1)
The Shewhart X ¯ -charts of the annual SO3 mean values and the nonparametric Mann–Whitney statistical test prove that using the FF–FB controller, the mean values approach better the SO3 target than the SC controller in two out of the three CEM types produced continuously for eighteen years;
(2)
FF–FB is better than SC in target approximation with a probability of 95% (a = 0.05) in CEM II A-L 42.5. The two controllers do not show distinguishable performance for the same test level a in CEM II B-M 32.5. This resulted from the second CEM type reduced clinker content and the consequent milder variance of clinker SO3 within the cement composition. In contrast, the ability of FF–FB to regulate gypsum is better than SC so the SO3 values are closer to the target and appear in CEM II A-L 42.5. Compared with CEM II B-M 32.5, this cement has a higher clinker content, which causes a higher variation in SO3 within the composition. The enhanced performance of FF–FB is more distinct in the pozzolanic cement, where clinker and fly ash are the two independent sources of sulfate disturbances because the test rejects the null hypothesis of equivalence with a probability of 99%;
(3)
The Shewhart s-charts of the annual standard deviation per CEM type and CM show that the FF–FB controller performs substantially better than the SC. The U C L of the former is always lower than the LCL of the latter. The ratio of the central lines of FF–FB to SC ranges from 0.51 to 0.59 for the Portland CEM types. This ratio is further reduced to 0.39 in pozzolanic cement CEM IV B (P-W) 32.5, where SO3 disturbances originate from fly ash and clinker;
(4)
Our analysis illustrates that the error propagation method is appropriate for comparing controller performance. If a controller regulates the cement sulfates by changing the gypsum, the SO3 contained in the gypsum and clinker are negatively correlated. The same occurs for SO3 in gypsum and fly ash. The larger the absolute value of correlation coefficients, the more robustly the controller regulates the gypsum content to compensate for clinker or fly-ash sulfate disturbances or changes in cement composition. The coefficients r(CL, G) and r(FA, G) are 0.876 and 0.006, respectively, using the SC controller. In the FF–FB case, the values are essentially higher, 0.962 and 0.647, respectively. The above clearly explains the higher performance of the feedforward–feedback system compared with SC.
To the best of the author’s knowledge, it is hard to find a description of cement sulfate controllers installed in milling systems in the literature. The technical novelty of this research is the design and long-term industrial implementation of such a system comprising a feedforward and feedback component. The FF component tracks the set-point changes of SO3 when the cement mill starts to grind another CEM type, and the FB controller effectively attenuates the variations of the raw materials’ sulfates.
Cement factories today use various alternative fuels to reduce their carbon footprint per clinker ton. Their highly changeable mix composition and sulfur level increase the SO3 variance of clinker, making it indispensable to implement an optimized controller to regulate sulfates around the target. An optimal SO3 target provides the maximum compressive strength [19] and permits clinker reduction in cement composition, further contributing to the decrease in CO2 per ton of product. Consequently, the actions to optimize and regulate SO3 are interconnected, resulting in a positive environmental impact.
We selected spot cement sampling and a sampling period such that transient dynamic phenomena are negligible. Further development of the research on optimal controllers regulating SO3 in the CM outlet involves designing a control system where the sampling is continuous, providing an average sample for each Ts period. This type of controller must account for the transient phenomena that occur during the mean sample preparation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/chemengineering8020033/s1, Table S1. data for the calculation of correlation coefficients.

Funding

This research received no external funding.

Data Availability Statement

The data and results presented in this paper are available upon request from the authors.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. EN 197-1:2011; Cement—Part 1: Composition, Specifications and Conformity Criteria for common Cements. CEN/TC 51; CEN: Brussels, Belgium, 2011; pp. 10–15.
  2. C150/C150M-22; Standard Specification for Portland Cement. ASTM International: West Conshohocken, PA, USA, 2022.
  3. Mindess, S.; Young, J.F.; Darwin, D. Concrete, 2nd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2003; Volume 31, pp. 60–62. [Google Scholar]
  4. Niemuth, M. Effect of Fly Ash on the Optimum Sulfate of Portland Cement. Ph.D. Dissertation, Purdue University, West Lafayette, IN, USA, December 2012; p. 68. Available online: https://www.researchgate.net/publication/266077758_Effect_of_Fly_Ash_on_Optimum_Sulfate_of_Portland_Cement (accessed on 23 November 2023).
  5. Evans, K.A. The Optimum Sulphate Content in Portland Cement. p. 3. Available online: https://tspace.library.utoronto.ca/bitstream/1807/11621/1/MQ29389.pdf (accessed on 7 April 2023).
  6. Lerch, W. The Influence of Gypsum on the Hydration and Properties of Portland Cement Pastes. In Proceedings of the American Society for Testing Materials; American Concrete Institute (ACI): Farmington Hills, MI, USA, 1946; Volume 46, pp. 1252–1291. [Google Scholar]
  7. Bentur, A. Effect of Gypsum on the Hydration and Strength of C3S Pastes. J. Am. Ceram. Soc. 1976, 59, 210–213. [Google Scholar] [CrossRef]
  8. Soroka, I.; Abayneh, M. Effect of gypsum on properties and internal structure of PC paste. Cem. Concr. Res. 1986, 16, 495–504. [Google Scholar] [CrossRef]
  9. Sersale, R.; Cioffi, R.; Frigione, G.; Zenone, F. Relationship between gypsum content, porosity and strength in cement. I. Effect of SO3 on the physical microstructure of Portland cement mortars. Cem. Concr. Res. 1991, 21, 120–126. [Google Scholar] [CrossRef]
  10. Gunay, S.A.A. Influence of Aluminates Phases Hydration in Presence of Calcium Sulfate on Silicates Phases Hydration: Consequences on Cement Optimum Sulfate. Ph.D. Thesis, University of Bourgogne, Bourgogne, France, 2012. Available online: https://theses.hal.science/tel-00767768 (accessed on 7 April 2023).
  11. Zunino, F.; Scrivener, K. The influence of sulfate addition on hydration kinetics and C-S-H morphology of C3S and C3S/C3A systems. Cem. Concr. Res. 2022, 160, 106930. [Google Scholar] [CrossRef]
  12. Andrade Neto, J.S.; de Matos, P.R.; De la Torre, A.G.; Campos, C.E.M.; Torres, S.M.; Monteiro, P.J.M.; Kirchheim, A.P. Hydration and interactions between pure and doped C3S and C3A in the presence of different calcium sulfates. Cem. Concr. Res. 2022, 159, 106893. [Google Scholar] [CrossRef]
  13. Mohammed, S.; Safiullah, O. Optimization of the SO3 content of an Algerian Portland cement: Study on the effect of various amounts of gypsum on cement properties. Constr. Build. Mater. 2018, 164, 262–370. [Google Scholar] [CrossRef]
  14. Adu-Amankwah, S.; Black, L.; Skocek, J.; Ben Haha, M.; Zajac, M. Effect of sulfate additions on hydration and performance of ternary slag-limestone composite cements. Constr. Build. Mater. 2018, 164, 451–462. [Google Scholar] [CrossRef]
  15. Han, F.; Zhou, Y.; Zhang, Z. Effect of gypsum on the properties of composite binder containing high- volume slag and iron tailing powder. Constr. Build. Mater. 2020, 252, 119023. [Google Scholar] [CrossRef]
  16. Yamashita, H.; Yamada, K.; Hirao, H.; Hoshino, S. Influence of Limestone Powder on the Optimum Gypsum Content for Portland Cement with Different Alumina Content. Available online: https://www.researchgate.net/publication/285554157_Influence_of_Limestone_Powder_on_the_Optimum_Gypsum_Content_for_Portland_Cement_with_Different_Alumina_Content (accessed on 23 November 2023).
  17. Liu, F.; Lan, M.Z. Effects of Gypsum on Cementitious Systems with Different Mineral Mixtures. Available online: https://www.researchgate.net/publication/269647645_Effects_of_Gypsum_on_Cementitious_Systems_with_Different_Mineral_Mixtures (accessed on 23 November 2023).
  18. Fincan, M. Sulfate Optimization in the Cement-Slag Blended System Based on Calorimetry and Strength Studies. Ph.D. Thesis, University of South Florida, Tampa, FL, USA, 2021; p. 2. Available online: https://digitalcommons.usf.edu/cgi/viewcontent.cgi?article=9967&context=etd (accessed on 23 November 2023).
  19. Tsamatsoulis, D.C.; Korologos, C.A.; Tsiftsoglou, D.V. Optimizing the Sulfates Content of Cement Using Neural Networks and Uncertainty Analysis. ChemEngineering 2023, 7, 58. [Google Scholar] [CrossRef]
  20. Andrade Neto, J.S.; De la Torre, A.G.; Kirchheim, A.P. Effects of sulfates on the hydration of Portland cement—A review. Constr. Build. Mater. 2021, 279, 122428. [Google Scholar] [CrossRef]
  21. Tsamatsoulis, D. Simulation of Cement Grinding Process for Optimal Control of SO3 Content. Chem. Biochem. Eng. Q. 2014, 28, 13–25. Available online: http://silverstripe.fkit.hr/cabeq/past-issues/article/27 (accessed on 23 November 2023).
  22. Ko, Y.R.; Kim, T.H. Feedforward Plus Feedback Control of an Electro-Hydraulic Valve System Using a Proportional Control Valve. Actuators 2020, 9, 45. [Google Scholar] [CrossRef]
  23. Wang, X.; Li, J.; Lu, X. Design and Control of a Trapezoidal Piezoelectric Bimorph Actuator for Optical Fiber Alignment. Materials 2023, 16, 5811. [Google Scholar] [CrossRef] [PubMed]
  24. Araque, J.G.; Angel, L.; Viola, J.; Chen, Y. Design and Implementation of a Recursive Feedforward-Based Virtual Reference Feedback Tuning (VRFT) Controller for Temperature Uniformity Control Applications. Machines 2023, 11, 975. [Google Scholar] [CrossRef]
  25. Tsamatsoulis, D.C. Optimizing the Control System of Cement Milling: Process Modeling and Controller Tuning Based on Loop Shaping Procedures and Process Simulations. Braz. J. Chem. Eng. 2014, 31, 155. [Google Scholar] [CrossRef]
  26. Astrom, K.; Hagglund, T. Advanced PID Control; Instrumentation, Systems and Automatic Society: Research Triangle Park, NJ, USA, 2006; pp. 414–418. [Google Scholar]
  27. ISO 7870-2:2013; Control Charts—Part 2: Shewhart Control Charts. ISO/TC 69; ISO: Geneva, Switzerland, 2013; pp. 8–9.
  28. Joint Committee for Guides in Metrology/Working Group 1 (JCGM/WG 1) Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement. pp. 18–23. Available online: https://www.bipm.org/documents/20126/2071204/JCGM_100_2008_E.pdf/cb0ef43f-baa5-11cf-3f85-4dcd86f77bd6 (accessed on 23 November 2023).
  29. NIST/SEMATECH Engineering Statistics Handbook. Available online: https://www.itl.nist.gov/div898/handbook/pmc/section3/pmc321.htm (accessed on 9 November 2023).
  30. ISO 2854:1976; Statistical Interpretation of Data—Techniques of Estimation and Tests Relating to Means and Variances. ISO/TC 69; ISO: Geneva, Switzerland, 1976; pp. 18–19.
  31. Mann, H.B.; Whitney, D.R. On a Test of Whether one of Two Random Variables is Stochastically Larger than the Other. Ann. Math. Statist. 1947, 18, 50–60. [Google Scholar] [CrossRef]
  32. Jar, J.H. Biostatistical Analyis, 5th ed.; Pearson Education. Inc.: Upper Saddle River, NJ, USA, 2010; pp. 163–172. [Google Scholar]
  33. Deshpande, J.V.; Naik-Nimbalkar, U.; Dewan, I. Nonparametric Statistics; World Scientific Publishing Co. Pte. Ltd.: Hackensack, NJ, USA, 2017; pp. 114–143. Available online: https://lccn.loc.gov/2017029415 (accessed on 23 November 2023).
  34. Corder, G.W.; Foreman, D.I. Nonparametric Statistics: A Step-by-Step Approach, 2nd ed.; Wiley & Sons, Inc.: Hoboken, NJ, USA, 2014; pp. 69–80. [Google Scholar]
  35. Nachar, N. The Mann-Whitney U: A Test for Assessing Whether Two Independent Samples Come from the Same Distribution. Tutor. Quant. Methods Psychol. 2008, 4, 13–20. Available online: https://www.researchgate.net/publication/49619432_The_Mann-Whitney_U_A_Test_for_Assessing_Whether_Two_Independent_Samples_Come_from_the_Same_Distribution (accessed on 23 November 2023). [CrossRef]
  36. García-Marín, A.P.; Estévez, J.; Morbidelli, R.; Saltalippi, C.; Ayuso-Muñoz, J.L.; Flammini, A. Assessing Inhomogeneities in Extreme Annual Rainfall Data Series by Multifractal Approach. Water 2020, 12, 1030. [Google Scholar] [CrossRef]
  37. Rubarth, K.; Sattler, P.; Zimmermann, H.G.; Konietschke, F. Estimation and Testing of Wilcoxon–Mann–Whitney Effects in Factorial Clustered Data Designs. Symmetry 2022, 14, 244. [Google Scholar] [CrossRef]
  38. Wahyudi, R.D.; Singgih, M.L.; Suef, M. Investigation of Product–Service System Components as Control Points for Value Creation and Development Process. Sustainability 2022, 14, 16216. [Google Scholar] [CrossRef]
  39. Real Statistics Using Excel. Mann–Whitney Table. Available online: https://real-statistics.com/statistics-tables/mann-whitney-table/ (accessed on 23 November 2023).
  40. EN 196-2:2013; Methods of Testing Cement—Part 2: Chemical Analysis of Cement. CEN/TC 51; CEN Management Centre: Brussels, Belgium, 2013.
Figure 1. Flowchart of a closed grinding circuit.
Figure 1. Flowchart of a closed grinding circuit.
Chemengineering 08 00033 g001
Figure 2. Block diagram of sulfates control.
Figure 2. Block diagram of sulfates control.
Chemengineering 08 00033 g002
Figure 3. Comparison of SC and FF–FB systems for CEM II B-M (P-L) 32.5.
Figure 3. Comparison of SC and FF–FB systems for CEM II B-M (P-L) 32.5.
Chemengineering 08 00033 g003
Figure 4. Control charts of CEM II B-M (P-L) 32.5 produced in CM6: (a) X ¯ -chart and (b) s-chart.
Figure 4. Control charts of CEM II B-M (P-L) 32.5 produced in CM6: (a) X ¯ -chart and (b) s-chart.
Chemengineering 08 00033 g004
Figure 5. Control charts of CEM II A-L 42.5 produced in CM6: (a) X ¯ -chart and (b) s-chart.
Figure 5. Control charts of CEM II A-L 42.5 produced in CM6: (a) X ¯ -chart and (b) s-chart.
Chemengineering 08 00033 g005
Figure 6. Control charts of CEM II B-M (P-L) 32.5 produced in CM5: (a) X ¯ -chart and (b) s-chart.
Figure 6. Control charts of CEM II B-M (P-L) 32.5 produced in CM5: (a) X ¯ -chart and (b) s-chart.
Chemengineering 08 00033 g006
Figure 7. Control charts of CEM II A-L 42.5 produced in CM5: (a) X ¯ -chart and (b) s-chart.
Figure 7. Control charts of CEM II A-L 42.5 produced in CM5: (a) X ¯ -chart and (b) s-chart.
Chemengineering 08 00033 g007
Figure 8. Control charts of CEM IV B (P-W) 32.5 produced in CM5: (a) X ¯ -chart and (b) s-chart.
Figure 8. Control charts of CEM IV B (P-W) 32.5 produced in CM5: (a) X ¯ -chart and (b) s-chart.
Chemengineering 08 00033 g008
Figure 9. Frequency distributions of the yearly SO3 mean values for (a) CEM II B-M (P-L) 32.5 and (b) CEM II A-L 42.5.
Figure 9. Frequency distributions of the yearly SO3 mean values for (a) CEM II B-M (P-L) 32.5 and (b) CEM II A-L 42.5.
Chemengineering 08 00033 g009
Figure 10. yact and ycalc values for (a) SC and (b) FF–FB implementation.
Figure 10. yact and ycalc values for (a) SC and (b) FF–FB implementation.
Chemengineering 08 00033 g010
Table 1. CEM types.
Table 1. CEM types.
CEMConstituent (%) 128-Day Strength Limits (MPa)
ClinkerLimestone
(L)
Pozzolan
(P)
Fly Ash
(W)
MinorLowHigh
CEM I 52.5 N95–100 0–552.5
CEM II A-L 42.5 N80–94<-- 6–20 --> 0–542.562.5
CEM II B-M (P-L) 32.5 N65–79<-------- 21–35 --------> 0–532.552.5
CEM IV B(P) 32.5 N-SR45–64 <- 36–55 -> 0–532.552.5
CEM IV B(P-W) 32.5 N45–64 <-------- 36–55 ------->0–532.552.5
1 Gypsum is not included in the composition but added according to SO3 target.
Table 2. SO3 of raw materials.
Table 2. SO3 of raw materials.
ClinkerLimestonePozzolanFly AshGypsum
Count493103347719
Average0.930.020.02.4942.78
Std. Dev.0.330.020.01.432.88
%CV 35.5 44.86.7
Table 3. Simulator parameters.
Table 3. Simulator parameters.
CEM II B-M (P-L) 32.5CEM II A-L 42.5
SO3 target (%)2.53
Clinker content (%)6580
SO3 low limit (%) (=0.95 of SO3 target)2.3752.85
SO3 high limit (%) (=1.05 of SO3 target)2.6253.15
Initial gypsum (%)4.0
Initial SO3 (%)2.80
Sampling period (h)2.0
Clinker mean SO3 (%)0.93
Clinker SO3 standard deviation (%)0.20
Gypsum mean SO3 (%)42.78
Gypsum SO3 standard deviation (%)1.0
Table 4. Statistical results of simulator application.
Table 4. Statistical results of simulator application.
CEMII B-M (P-L) 32.5II A-L 42.5
SCFF–FBRatioSCFF–FBRatio
Standard deviation, sSC, sFF/FB (%)0.1340.109 0.1460.125
sFF/FB/sSC 0.812 0.854
(%) of population out of [0.95·SSP, 1.05·SSP] PSC, PFF/FB32.115.4 27.314.6
PFF/FB/PSC 0.481 0.534
(%) of population out of [0.95·SSP, 1.05·SSP] in CEM type change, CSC, CFF/FB38.19.6 36.55.6
CFF/FB/CSC 0.251 0.154
Table 5. Mann–Whitney test for the sets D1 and D2.
Table 5. Mann–Whitney test for the sets D1 and D2.
CEM II B-M (P-L) 32.5CEM II A-L 42.5 CEM IV B (P-W) 32.5
M116156
M2211410
U 140623
U c r for a = 0.0192518
U c r for a = 0.051136614
Table 6. Pooled standard deviations.
Table 6. Pooled standard deviations.
CEMII B-M (P-L) 32.5II A-L 42.5IV B (P-W) 32.5I 52.5IV B (P) 32.5
CM5 sP10.1450.1420.241
CM5 sP20.0860.0830.0950.1260.096
CM5 sP2/sP10.590.580.39
CM6 sP10.1510.134
CM6 sP20.0780.072
CM6 sP2/sP10.510.53
Table 7. Variables of error propagation model and type of quality data.
Table 7. Variables of error propagation model and type of quality data.
VariableType of Quality Data
uS,CEMAnnual standard deviation of the SO3 daily data
CL, G, FAAnnual average of clinker, gypsum, and fly-ash fractions in CEM composition calculated from daily data chemical analysis
uCL, uG, uFAAnnual standard deviation of clinker, gypsum, and fly-ash fractions in CEM composition calculated from daily data chemical analysis
SO3,CL, uSO3,CLAnnual average and standard deviation of clinker SO3 calculated from daily data
SO3,G, uSO3,G,
SO3,FA, uSO3,FA
Annual average and standard deviation of gypsum and fly-ash SO3 calculated from the samples taken in one year
Table 8. Regression results analysis of the standard deviation results.
Table 8. Regression results analysis of the standard deviation results.
ControllerCountr(CL, G) r(FA, G)Standard ErrorStandard DeviationR2
SC370.8760.0060.013590.035420.853
FF–FB500.9620.6470.008330.039880.956
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tsamatsoulis, D. An Industrial Control System for Cement Sulfates Content Using a Feedforward and Feedback Mechanism. ChemEngineering 2024, 8, 33. https://doi.org/10.3390/chemengineering8020033

AMA Style

Tsamatsoulis D. An Industrial Control System for Cement Sulfates Content Using a Feedforward and Feedback Mechanism. ChemEngineering. 2024; 8(2):33. https://doi.org/10.3390/chemengineering8020033

Chicago/Turabian Style

Tsamatsoulis, Dimitris. 2024. "An Industrial Control System for Cement Sulfates Content Using a Feedforward and Feedback Mechanism" ChemEngineering 8, no. 2: 33. https://doi.org/10.3390/chemengineering8020033

APA Style

Tsamatsoulis, D. (2024). An Industrial Control System for Cement Sulfates Content Using a Feedforward and Feedback Mechanism. ChemEngineering, 8(2), 33. https://doi.org/10.3390/chemengineering8020033

Article Metrics

Back to TopTop