International Journal of Control Science and Engineering
p-ISSN: 2168-4952 e-ISSN: 2168-4960
2022; 12(1): 1-25
doi:10.5923/j.control.20221201.01
Received: Jun. 2, 2022; Accepted: Jul. 6, 2022; Published: Aug. 23, 2022

Vincent A. Akpan1, Ioakeim K. Samaras2, George D. Hassapis3
1Department of Biomedical Technology, The Federal University of Technology, Akure, Ondo State, Nigeria
2Intracom Telecom, Software Development Center, 19.7 Km Markopoulou Avenue, Peania, Greece
3Department of Electrical and Computer Engineering, Aristotle University of Thessaloniki, Thessaloniki, Greece
Correspondence to: Vincent A. Akpan, Department of Biomedical Technology, The Federal University of Technology, Akure, Ondo State, Nigeria.
| Email: | ![]() |
Copyright © 2022 The Author(s). Published by Scientific & Academic Publishing.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/

Industrial control networks play significant roles in industrial distributed networked control systems since it enables all the system components to be interconnected as well as monitor and control the physical equipment in industrial environments [2]. The integration of control and communication in networked control systems (NCSs) has made the design and analysis of NCSs a great theoretical challenge for conventional control theory [11]. A major trend in modern industrial and commercial control systems is to integrate computing, communication and control into different levels of machine/factory operations and information processes [12]. NCSs provide a natural platform for distributed learning control. However, it seems that, apart from the remote-tuning of PID controllers, there is no strong research activity to combine NCS study with adaptive control, learning control, network communications and so on [12]. The transmission time for the data packets introduces network-induced delays to NCSs, which are well known to degrade the performance of the control systems [11]. This paper presents a NCS implemented over on a specific service-oriented-architecture (SOA) computer network based on device profile web services (DPWS) for industrial control applications with emphasis on reduced network-induced delay. The performance of the proposed NCS based on four-level hierarchical structure is compared with other networks by considering its real-time implementation for online neural network-based model identification and a nonlinear model adaptive model predictive control (NAMPC) of a fluidized bed furnace reactor (FBFR). Even though SOA connections offer flexibility and scalability advantages, their large communication overhead makes it difficult to satisfy the real-time requirements of the control algorithm when it is implemented over traditional Ethernet networks. However and contrary to [11], in the proposed NCS over the SOA computer network based on DPWS, every component conforms to a SOA technology while the exchange of messages follows a new format technique which significantly reduces transmission delays and overheads without prohibiting the SOA high level interfacing. Furthermore these components are interconnected with each other by utilizing the switched Ethernet architecture which further reduces the transmission delays. Simulation results for the FBFR pilot plant model identification and control have shown that the proposed computer network allows the satisfaction of the real-time constraints of the considered model-based predictive control of the FBFR pilot plant. The aforementioned results render the proposed computer network suitable for advanced control of industrial processes with time constants similar to those of the FBFR process.
Keywords: Device profile for web services (DPWS), Network control system (NCS), Neural Networks, Nonlinear Adaptive Model Predictive Control (NAMPC), Nonlinear Model Identification, Service Oriented Architecture (SOA), Switched Ethernet
Cite this paper: Vincent A. Akpan, Ioakeim K. Samaras, George D. Hassapis, Implementation of Distributed Network Control System over a Service-Oriented-Architecture Computer Network Based on Device Profile for Web Services for Industrial Control Applications, International Journal of Control Science and Engineering, Vol. 12 No. 1, 2022, pp. 1-25. doi: 10.5923/j.control.20221201.01.
![]() | Figure 1. The proposed four-level network control system (NCS) architecture based on hierarchical structure |
|
|
![]() | Figure 2. DPWS WS-Eventing notification message transformation according to the IASFT specification |
sensors and
actuators transmit data simultaneously. This delay can be defined as follows:![]() | (1) |
is the processing transmission delay at the sensors and actuators,
is the transmission delay, i.e. the delay in queue plus the delay in the network, and
is the frame reception delay at the control system. When a TCP connection is established, its delay must also be taken into account. This connection is established by the exchange of a CONNECTION-REQUEST message and a CONNECTION-ACCEPTED segment as it is documented in [31]. Again the worst case TCP connection establishment delay is observed when all the devices that are located at the device level request such a connection simultaneously. Therefore the worst case CONNECTION-REQUEST delay from the device level to control system is: ![]() | (2) |
is the TCP processing transmission delay at the sensors and actuators while
is TCP processing reception delay at the control system. These delays correspond to the flow of data from the transport layer to the PHY layer and vice versa.
is the transmission delay of the TCP request segment transmitted from the device level to the control system. The worst case CONNECTION-ACCEPTED delay which is denoted as
is the delay experienced when the last CONNECTION-ACCEPTED segment is sent by the control system to device level and is the same with
. Now, the worst case transmission delay that a TCP data segment experiences when it is transmitted from the device level to the control system is defined by using Equations (1) and (2) as: ![]() | (3) |
be the overall processing delay that the TCP data segment experiences when it is sent from the device level to the control system and
be the overall transmission delay that a TCP data segment experiences for the same path. Then Equation (3) is formed as:![]() | (4) |
where
is the computational time that the identification and control algorithms need for computing the control input signals to the process,
is the processing transmission delay at the control system and
is the frame reception delay at the device level and
is the transmission delay a frame sent from the control system to the device level. When a TCP connection is established, it must be taken into account the TCP connection establishment delay too. Following the same way with the one presented previously, the worst case transmission delay that a TCP data segment experiences sent from the control system to the device level is defined as:![]() | (5) |
is the transmission delay of the TCP request segment transmitted from control system to the device level. Let
be the overall processing delay a TCP data segment experiences due to transmission by the control system to the device level and let
be the overall transmission delay a TCP data segment experiences for the same path [2,11,13,14]. Then Equation (5) is formed as:![]() | (6) |
![]() | (7) |
as this is the bounded delay that it can offer as long as there are no overflow events in the switches.![]() | Figure 3. Structure of a SOAP message |
![]() | Figure 4. Simplified diagram of the steam deactivation unit (SDU) of the FCC pilot plant |
![]() | Figure 5. Schematic of the vertical cross-section of the cylindrical fluidized bed furnace reactor (FBFR) of the SDU |
, the FBFR process can be represented by as a p-input q-output nonlinear discrete-time system with disturbance term
by the following nonlinear autoregressive moving average with exogenous inputs (NARMAX) model [32–34]:![]() | (8) |
is a nonlinear function of its arguments, and
are the past input vector,
are the past output vector,
is the current output,
and
are the number of past inputs and outputs respectively that define the order of the system, and
is time delay. The predictor form of Equation (8) based on the information up to time
can be expressed as [32–34]:![]() | (9) |
is the regression (state) vector, and
is an unknown parameter vector which must be selected such that
,
is the error between Equations (8) and (9) defined as:![]() | (10) |
![]() | (11) |
is the sampling time of the system outputs. Then, the minimization of Equation (10) can be stated as follows:![]() | (12) |
is formulated as a mean square error (MSE) type cost function which can be stated as:![]() | (13) |
as an argument in
is to account for the desired model
dependency on
. Thus, given an initial small random value of
,
and Equation (11), the model identification problem reduces to the minimization of Equation (12) to obtain
.The minimization of Equation (12) is approached by considering
as a neural network model. The complete NN model identification scheme based on the teacher-forcing method is illustrated in Figure 6 [34]. According to this scheme, the validated mathematical model of the FBFR process is placed in parallel with its NNARMAX model (in the dashed box) where the TDL (tapped delay line memory) are used to store temporal NN input information. The architecture of the “Neural Network Model” of Figure 6 is a dynamic feedforward NN (DFNN) shown in Figure 7. The inputs to the DFNN model of Figure 7 are
,
and
which are concatenated into
as shown in Figure 7. The output of the NN model of Figure 6 in terms of the network parameters of Figure 7 is given as:![]() | (14) |
and
are the number of hidden neurons and number of regressors respectively; i is the number of outputs,
and
are the hidden and output weights respectively;
and
are the hidden and output biases;
is a linear activation function for the output layer and
is an hyperbolic tangent activation function for the hidden layer defined here as:![]() | Figure 6. The FBFR NNARMAX model identification scheme using ARLS training algorithm |
![]() | Figure 7. Architecture of the dynamic feedforward NN (DFNN) model |
![]() | (15) |
is a collection of all network weights and biases in Equation (14) in term of the matrices
and
. Equation (14) is here referred to as NN NARMAX (NNARMAX) model predictor for simplicity.The
in Equation (8) is usually unknown but can be estimated as a covariance noise matrix using
with an iterative algorithm described in [32–34]. Thus, using, Equation (13) becomes:![]() | (16) |
is a penalty norm and also removes ill-conditioning, where
is an identity matrix,
and
are the weight decay values for the input-to-hidden and hidden-to-output layers respectively. Note that both
and
are adjusted simultaneously with
during network training and are used to update
iteratively [32–34].Several methods have been proposed in literatures for the minimization of Equation (10) [32–34]. The formulation of Equations (16) from (10) follows from the work of Akpan and co-workers where efficient training algorithms have been developed and validated [32–34]. The training algorithm adopted in this work is an online adaptive recursive least squares (ARLS) algorithm due to its proven efficiency [32–34].
and
are the desired reference signal, prediction error, control input, system output,
step-delay prediction model output,
step-ahead predicted output, noise/input disturbances and
step-delay operator respectively and
is the number of samples based on the new measurement data sample.![]() | Figure 8. The NN-based NAMPC scheme for the FBFR process with the NNARMAX neural network model |
![]() | (17) |
and
are the desired reference and the filtered reference signals respectively; Am and Bm are the denominator and numerator polynomials of the filter. In this way, the NAMPC is deigned, in part, based on the filter tracking error capability; where Am and Bm serves as tuning parameters used to improve the robustness and internal stability of the NAMPC algorithm respectively.The NAMPC strategy is based on a receding horizon principle illustrated in Figure 9 and depends on an explicit model of the system. The NAMPC strategy of Figure 9 is summarized as follows:1). At the current sampling time k, the NNARMAX model predictor uses the past m-inputs, n-outputs and the current system information to identify the nonlinear discrete-time NNARMAX model of the FBFR process.2). Assuming that at Assuming that the identified NNARMAX model is stable, proper and deterministic; then the NAMPC algorithm uses the linearized model parameters of the identified nonlinear NNARMAX model to accurately predict the current system output
at that same sample time instant k.3). At time
, the NAMPC algorithm calculates a sequence of control inputs
consisting of the current
and future inputs
. The current input
is held constant after
control moves; where
is the maximum control horizon. The input
is calculated such that a set of
approaches the desired reference signal in an optimal manner over a specified prediction horizon
where
and
are the minimum and maximum prediction horizons respectively.4). The predicted values are used to calculate the control signals by minimizing an objective function of the form:![]() | (18) |
![]() | Figure 9. The NAMPC strategy |
![]() | (19) |
![]() | (20) |

where
is the change in control signal;
and
are two weighting matrices penalizing changes on
and
in Equation (18). Although a sequence of
moves is calculated at each sampling instant, only the first control move
is actually implemented and applied to control the process. The remaining control signals are not applied because at the next sampling instant
a new output
is known based on new measurements. The NAMPC strategy enters a new optimization loop while the remaining control signals
are used to initialize the optimizer. This is indeed the receding horizon principle inherent in MPC strategy.The NAMPC algorithm based on the full-Newton optimization used in this work is taken from the [33,34], where it has been demonstrated to be suitable for the control of the FBFR process [33]. Hence, the effort in this work is directed towards the online closed-loop implementation of the model identification and NAMPC schemes on the proposed NCS based on SOA technology via DPWS clients and servers to reduce the computation time at each sampling instant as well as reduced wiring and connections over long distance between the control computer and the industrial nonlinear FBFR process.![]() | (21) |
and
are the proportional, integral and derivative gains respectively,
is the sampling time and
is the error term defined as the difference between the process
and desired reference
given as![]() | (22) |
The controlled outputs of the FBFR are the temperatures of the six sections of the FBFR namely: reactor’s interior
, interior reactor wall
, air gap between the reactor and the heater
, heater
, insulator
, and outer reactor metal wall
; and is expressed as
In the work carried out in [35], a well-tuned PID and MPC controllers were developed to operate the FBFR process in the range of 3.76 and 3.66 kW respectively out of the total 5.04 kW heat energy available for the process and a sampling time (T) of 2 minutes was considered for 22 hours operating cycles.This means that 1320 data samples can be obtained from the process in one minute. In order to develop a neural network to accurately model the FBFR process in the present study, the heat supplied (Q = 5.04 kW) is varied by + 20% in order to cover the entire operating range of the pilot plant, during both initial heat-up and deactivation, and to account for the possible uncertainties in the plant model outside the operating region.Thus, the ± 20% lower and upper values of Q are 3.528 kW and 6.552 kW respectively. Using these values of Q, the complete validated mathematical model of the FBFR process is simulated in open-loop with a sampling step of 1 minute to obtain 1320 input-output data pairs for the NN training while 300 input-output validation data is obtained from the FBFR pilot plant under normal operating condition.The input vector
to the NNARMAX consist of the current state input regression vector
and the regression vector of the six output state derivatives
. The desired outputs of the NNARMAX model are the predicted values of Th and Tri given as
. In this case, the change in the values of Q affects the FBFR outputs and can be viewed as external disturbances
(see Figure 1).The training data is scaled to zero mean and unit variance using their mean values and standard deviations according to the following equations:![]() | (23) |
and
are the mean and standard deviation of the input and output training data pair; and
and
are the scaled inputs and outputs respectively.The scaling based on Equation (23) is to prevent signals of large magnitudes from dominating the identified model [34]. The network is trained for 20 epochs with the following parameters selected as: j=6, l=2, i=2, m=2, n=2 and
. The four design parameters for the ARLS algorithm are selected to be: α=0.8, β=1.2,
and π=0.98 resulting in γ=0.0204 which gives initial values for
and
equal to 0.8384 and 1.0788 respectively.After the network training, the joint weights are rescaled afterwards according to the following expression![]() | (24) |
![]() | Figure 10. Convergence of the ARLS training algorithm for 20 iterations for the FBFR process |
![]() | Figure 11. Comparison of 5-step ahead output predictions by the trained network (red --*) with the original unscaled training data (blue -) for (a) and (b) ![]() |
defined [34,35] is:![]() | (25) |
. The optimal tuning parameters obtained for the PID are and NAMPC controllers are given in Table 4. Next, the FBFR NNARMAX model is used to simulate the PID and NAMPC controllers in open-loop.
|
|
![]() | Figure 13. Online closed-loop implementation of the NNARMAX model identification and adaptive control of the FBFR using the proposed NCS based on SOA computer network via DPWS clients and servers |
in a first-in first-out fashion as discussed in sub-section 4.1 and they also constitute the
and
data that are the inputs to the NNARMAX identification scheme of Figure 7. Since
has been selected in sub-section 5.1.1, thus the NNARMAX model input vector
consists of the current input and output states of the FBFR process at each time sample. The current inputs
as well as the new outputs
produced by the FBFR process due to changes in Q are delivered to the proposed identification and control scheme over six networks: a DPWS-based Ethernet network (first network), a DPWS-based Ethernet network which uses the EXI format (second network), a DPWS-based Ethernet network which uses the IASFT format (third network), a DPWS-based switched Ethernet network (fourth network), a DPWS-based switched Ethernet network which uses the EXI format (fifth network) and a DPWS-based switched Ethernet network which uses the IASFT format (sixth network). The first three networks reside in the Ethernet networks while the last three reside in the switched Ethernet networks used in this study. The sixth network constitutes the proposed computer network. All of them consist of six sensors corresponding to six state output derivatives, two actuators for the current FBFR process inputs and one component associated with the proposed identification and control scheme. When Ethernet network is used the aforementioned components are interconnected with each other through an Ethernet bus while when the switched Ethernet is used three Ethernet switches are utilized for interconnections complying by this way with the architecture depicted in Figure 1. The performance of all networks is studied during eventing and control level interactions due to the reasons explained in Section 2.In all networks the sensors and actuators are DPWS servers and transmit data to the proposed identification and control scheme by using the eventing level interaction while the control system is a DPWS client and communicates with the actuators by utilizing the control level interaction. All the aforementioned transactions are accomplished at each sample time. During these interactions HTTP is used and so TCP connections are established. Therefore Equation (7) can be used for calculating the worst case overall control loop delay when the switched Ethernet is used. The same equation can be used for determining the worst case overall control loop delay when the Ethernet network is used, as long as
and
are the overall processing delays that TCP data segments experience in the Ethernet network.A simulation study has been made using the network simulator (NS)–2 version 2.34 [37] for determining the
and
in Equation (7). NS-2 supports transmissions and receptions of packets on different wires and nodes regenerate the information and only forwards it to the port on which the destination is attached. So the switched Ethernet architecture is supported. Furthermore NS-2 implements IEEE 802.3 specification and therefore the Ethernet network can be used. The adoption of the simulator is based on the fact that the worst case transmission delay requires the simultaneous transmission of data as discussed in Section 2. As soon as the simulator provides better synchronization between nodes, more accurate results could be obtained than from a real network. The simulation was made for 120 sampling periods by using half and full duplex links as a medium for connecting the components in Ethernet networks and in the switched Ethernet networks respectively. The link capacity was set to 10Mbps, and the propagation delay to 0.1 μs. In the simulator, the node that corresponds to the proposed identification and control algorithms was programmed to transmit data to the actuators as soon as it receives the current state input vector and the vector of the six output state derivatives. Moreover, a constant bit rate (CBR) application was developed that produces the eventing and control messages with 1 minute rate. This application has been developed over TCP/IP technology. Lastly, in all networks ten more nodes were used as traffic generators. These generators transmit 256 bytes of data every 1 ms in order to add additional traffic.In Table 5 are listed the data volume of the exchanged messages that were produced from the FBFR control application from the different formats. Next a comparison is made between the six networks in order to verify the efficacy of the proposed NCS based on SOA computer network.
|
due to the non-deterministic characteristic of the Ethernet network. As it can also be seen in Figure 14, in the first network the transmission delay some times exceeds the sampling period of the FBFR process.![]() | Figure 14. delay between the FBFR process and the control system obtained by NS-2 when the Ethernet networks are used |
is observed to be 77.21 seconds.
which is the DPWS protocol stack response time and is defined to be approximately 10 ms [37]. Also
and is observed to be approximately 300 μs while the average
is approximately 2.7 seconds as evident in Figure 15. All the aforementioned delays were computed using an Intel® Core™ 2 CPU running at 2.66GHz. Therefore, the worst case overall control loop delay in the first network is calculated using Equation (7) to be equal with 79.95 seconds. So the DPWS-based Ethernet network cannot fulfill the real time characteristics of the FBFR process as evident in the online step response simulation result of Figure 16. The poor performance of the NAMPC in tracking the desired reference is due to the transmission delay introduced by the network. As depicted most especially in Figure 16(a), the NAMPC sometimes tracks and sometimes does not track the desired reference according to the transmission delay introduced by the network which is sometimes below or above the sampling time of the FBFR process (see Figure 14).![]() | Figure 15. Computation time for the online FBFR model identification and control at each time sample |
![]() | Figure 16. Online identification and control of the FBFR process implemented over the first network: (a) Th and (b) Tri predictions with their respective control signals (c) HRP and (d) DWP |
is observed to be 37 seconds. Here
is the DPWS protocol stack response time plus the execution time for encoding the DPWS notification message to the EXI notification message. This delay was observed to be 22 ms. Moreover
is the DPWS protocol stack response time plus the execution time for decoding the EXI notification message to the DPWS notification message. This delay was observed to be 18 ms.
is the DPWS protocol stack response time plus the execution time for encoding the DPWS control message to the EXI control message and is 13 ms while
is the DPWS protocol stack response time plus the execution time for decoding the EXI control message to the DPWS control message and is 22 ms.
and
seconds. All the aforementioned delays were computed using an Intel® Core™ 2 CPU running at 2.66GHz. Therefore, the worst case overall control loop delay in the second network is calculated using Equation (7) to be equal with 39.81 seconds.In the third network, the maximum
is observed to be 48.9 seconds. Here
is the DPWS protocol stack response time plus the execution time of the XSLT processor for producing the IASFT notification message from the DPWS notification message. This delay was observed to be 13 ms while in Figure 2 it has been shown how the IASFT mechanism transforms the DPWS notification message from the FBFR application to the equivalent IASFT notification message. Moreover
is the DPWS protocol stack response time plus of the XSLT processor for producing the DPWS notification message from the IASFT notification message. This delay was observed to be 5 ms.
is the DPWS protocol stack response time plus the execution time of the XSLT processor for producing the IASFT control message from the DPWS control message and is 14 ms while
is the DPWS protocol stack response time plus the execution time of the XSLT processor for producing the DPWS control message from the IASFT control message and is 5 ms.
and
= 2.7 seconds. All the aforementioned delays were computed using an Intel® Core™ 2 CPU running at 2.66GHz. Therefore, the worst case overall control loop delay in the third network is calculated using Equation (7) to be equal with 51.67 seconds.In the second and third network the worst case overall control loop delay is below the sampling period of the FBFR process. Therefore these two networks fulfill the real time requirement of the FBFR process as shown in the online step response simulation result of Figure 17 where the NAMPC tracks the desired reference signal at each sampling instant. When the EXI format is used the DPWS-based Ethernet network introduces the minimum worst case overall control loop delay between the Ethernet networks. However the interoperability feature is lost as binary-based encoding is utilized. On the other hand, the IASFT renders the DPWS-based Ethernet network suitable for controlling the FBFR process while XML and DPWS standards are not distorted.
as well as the data volumes of the exchanged messages are the same with the ones computed for the first, second and third networks respectively.The simulation results obtained using the NS-2 are shown in Figure 18. In the fourth, fifth and sixth network the
was observed to be 0.023 seconds, 0.012 seconds and 0.014 seconds respectively at each sampling time (except from the first one). Therefore the worst case overall control loop delays for the fourth, fifth and sixth networks are calculated using Equation (7) to be equal with 2.77, 2.82 and 2.78 seconds respectively.![]() | Figure 18. delay between the FBFR process and the control system obtained by NS-2 when the switched Ethernet networks are used |
![]() | Figure 19. Round trip delays between DPWS clients (on Tri and Th controllers) and sensors (a) Tri, (b) Tirw, (c) Tbrwh, (d) Th, (e) Tins and (f) Tormw |