In this paper, we focus on determining the safe operational domain of a coupled actuator–valve configuration. The so-called “smart valves” system has increasingly been used in critical applications and missions including municipal piping networks, oil and gas fields, petrochemical plants, and more importantly, the U.S. Navy ships. A comprehensive dynamic analysis is hence needed to be carried out for capturing dangerous behaviors observed repeatedly in practice. Using some powerful tools of nonlinear dynamic analysis including Lyapunov exponents and Poincaré map, a comprehensive stability map is provided in order to determine the safe operational domain of the network in addition to characterizing the responses obtained. Coupled chaotic and hyperchaotic dynamics of two coupled solenoid-actuated butterfly valves are captured by running the network for some critical values through interconnected flow loads affected by the coupled actuators' variables. The significant effect of an unstable configuration of the valve–actuator on another set is thoroughly investigated to discuss the expected stability issues of a remote set due to others and vice versa.

## Introduction

Multidisciplinary electromechanical–fluid systems have been widely used in many megascale networks. Municipal piping systems, oil and gas fields, petrochemical plants, and more critically, the U.S. Navy are the immediate ones which need to utilize a reliable, safe, and efficient coupled flow distribution network.

Future smart cities would inevitably need an autonomous flow control network in order to help improve the safety of such a critical system and also to decrease the incremental costs of operation and maintenance. Malfunctions of the flow network have occurred repeatedly resulting in the flow interruption of small towns/districts. Although significant cost and energy would be needed to be spent in order to restore the whole system in addition to human resources required to be recruited. Economical and even social impact of these malfunctions can be expected to be dramatic and consequently, a fully automated flow network is needed to be designed and operated.

The same issues exist for oil and gas fields and petrochemical plants. The industry of oil and gas is one of the most sensitive elements of the global economy and plays important roles even in global politics. The flow control network is the essential part of these fields and is therefore required to be safely designed to minimize the flow interruptions leading to much higher oil/gas production.

The U.S. National Defense and Homeland Security is undoubtedly a highly important priority which requires to be addressed and investigated carefully. The U.S. Navy broadly employs the network of coupled electromechanical valve sets for cooling purposes, mainly for chilled water systems. The proper performance of the network is remarkably effective for other critical units of radar, sonar, defense systems, etc.

We have carried out broad analytical and experimental studies from nonlinear modeling to design optimization of both an isolated and interconnected symmetric butterfly valves driven by solenoid actuators [1–8]. The multidisciplinary couplings, including electromagnetics and fluid mechanics, had to be thoroughly regarded in the modeling phase in order to yield an accurate nonlinear model of such a complex system. A third-order nondimensional dynamic model of the single set was derived to be used in nonlinear dynamic analysis [3] and optimal design [4].

The dynamic analysis expectedly yielded practically observed crisis and transient chaotic dynamics of a single actuated valve for some critical physical parameters. A comprehensive stability map was also presented as an efficient tool to determine the safe domain of operation which in turn could serve for identifying the lower and upper bounds of the design optimization efforts. The design optimization was then carried out [4] to select the optimal actuation unit's parameters coupled with the mechanical and fluid parts in order to significantly reduce the amount of energy consumption (upward of 40%).

Note that the applications addressed earlier contain thousands of actuated valves in which a high level of dynamic coupling has been repeatedly observed in practice. These dynamic couplings among different sets need to be captured through analytical studies. We have developed [5] a novel nonlinear model for two sets of solenoid-actuated butterfly valves operating in series. The closing/opening valves were modeled as changing resistors and the flow between them as a constant one. A sixth-order nonlinear coupled model revealed the high dynamic sensitivity of each element of a set, the valve and the actuator, to another one and vice versa. The power spectrum was used in confirming the same frequency response of a neighbor set due to the external periodic noise applied on another set of the valve and actuator.

In further studies, we optimized the design of coupled actuation units of two sets operating in series [8] subject to a sudden contraction. The pipe contraction imposed an additional resistance to be modeled, and therefore, the coupled dynamic equations derived in Ref. [5] had to be slightly modified. We represent the modeling process here for completeness. We have surprisingly established an interesting coupling between currents of the actuation units through the interconnected flow loads, including hydrodynamic and bearing torques, which affect both the valves' dynamics.

Important nonlinear phenomena in electromechanical systems have also received considerable attention. Sun et al. [9] studied the hyperchaotic behavior of the newly presented simplified Lorenz system by using a sinusoidal parameter variation and hyperchaos control of the forced system via feedback. Banerjee et al. [10] investigated the synchronization of chaos and hyperchaos in first-order time-delayed systems that are coupled using the nonlinear time-delay excitatory coupling by assigning two characteristic time delays: the system delay that is same for both the systems and the coupling delay associated with the coupling path. Many efforts for analyzing hyperchaotic dynamics have been reported in Refs. [11–26].

We have addressed the modeling process subject to the pipe contraction in Ref. [8] and represent here for completeness. The contribution of this work is the inclusion of interconnected electromechanical–fluid nonlinearities between two actuator–valve configurations to thoroughly analyze the effects of an unstable set of the valve–actuator on another one. Through this comprehensive analysis, chaotic and hyperchaotic dynamics of two coupled configurations are captured by exposing the network to some critical values. The responses are then characterized using some powerful tools including Lyapunov exponents and Poincaré map.

## Mathematical Modeling

Figure 1(a) shows a pair of symmetric butterfly valves driven by solenoid actuators through rack and pinion arrangements. The rack and pinion mechanism provides a kinematic constraint which connects the dynamics of the valve and actuator. Applying DC voltages, as being used in the Navy ships for chilled water systems, the motive forces give translational motions to the actuators' moving parts (plungers), and subsequently, the valves rotate to desirable angles. Note that a return spring has been a common practice among industries to open the valves.

Interconnected modeling of such a multiphysics system undoubtedly needs some simplifying assumptions to neglect useless and tremendously time-consuming numerical calculations. The magnetic force resulted from the magnetic field needs an extremely short period of time to reach its maximum value. This period is the so-called “diffusion time” and has an inverse relationship with the amount of current used. Note that using current of 4 A would yield a negligible diffusion time of $\tau d\u224820\u2009ms$ [1] with respect to the nominal operation time of 40 s. We have to also assume dominant laminar flow for both the coupled valves. Note that developing an analytical model is a necessity to carry out the dynamic analysis and optimization which would lead us to make such a commonly used assumption and also to avoid the numerical difficulties involved with a turbulent regime. However, a crucial question needs to be carefully answered with respect to the validity of such an assumption. Using the values of pipe diameter and flow mean velocity listed in Table 1, one can easily distinguish the existence of the turbulent regime which invalidates the assumption we have made. From another aspect, the analytical formulas derived for the flow loads, including the hydrodynamic and bearing torques, have been developed based on the assumption of laminar flow [27,28]. To address the issues discussed above, we have carried out experimental work to measure the sum of the hydrodynamic and bearing torques as the most affecting loads on the valves' and subsequently actuators' dynamics [8]. The experiment yielded the total torque (Fig. 2) for the inlet velocity of $v\u22482.7\u2009m/s$ and valve diameter of $Dv=2\u2009in.$ validating the laminar flow assumption [29]. The flow torques have shown highly important roles for the dynamics of an isolated solenoid-actuated butterfly valve and we hence expect to observe such effects for the interconnected sets [5].

Ρ | 1000 kg/m^{3} | v | 0.1 m/s |

$J1,2$ | $0.104\xd710\u22121$ kg m^{2} | N_{2} | 3000 |

N_{1} | 3000 | $C11,22$ | 1.56 × 10^{6}^{ }H^{−1} |

$gm1,m2$ | 0.1 m | $V1,2$ | 24 V |

$Dv1$ | 0.2032 m | $Dv2$ | 0.127 m |

$Ds1,s2$ | 0.01 m | P_{out} | 2 kPa |

$k1,2$ | 1000 N·m^{−}^{1} | $C21,22$ | 6.32 × 10^{8}^{ }H^{−1} |

L_{1} | 2 m | L_{2} | 1 m |

μ_{f} | 0.018 kg m^{−1} s^{−1} | $R1,2$ | 6 Ω |

$r1,2$ | 0.05 m | θ | 90 deg |

P_{in} | 256 kPa |

Ρ | 1000 kg/m^{3} | v | 0.1 m/s |

$J1,2$ | $0.104\xd710\u22121$ kg m^{2} | N_{2} | 3000 |

N_{1} | 3000 | $C11,22$ | 1.56 × 10^{6}^{ }H^{−1} |

$gm1,m2$ | 0.1 m | $V1,2$ | 24 V |

$Dv1$ | 0.2032 m | $Dv2$ | 0.127 m |

$Ds1,s2$ | 0.01 m | P_{out} | 2 kPa |

$k1,2$ | 1000 N·m^{−}^{1} | $C21,22$ | 6.32 × 10^{8}^{ }H^{−1} |

L_{1} | 2 m | L_{2} | 1 m |

μ_{f} | 0.018 kg m^{−1} s^{−1} | $R1,2$ | 6 Ω |

$r1,2$ | 0.05 m | θ | 90 deg |

P_{in} | 256 kPa |

*q*is the volumetric flow rate,

_{v}*μ*indicates the fluid dynamic viscosity, $Dv1$ and $Dv2$ are the valves' diameters,

_{f}*L*

_{1}and

*L*

_{2}stand for the pipe lengths before and after contraction, $RL1$ and $RL2$ indicate the constant resistances, and $Pcon1$ and $Pcon2$ are the flow pressures before and after contraction.

*K*

_{con}is calculated as follows:

*R*)” and “coefficient (

*c*),” as important parameters of the regulating valves, are nonlinear functions of the valve rotation angle to be stated [32] as follows:

_{v}*α*indicates the valve rotation angle,

*ρ*is the density of the media, and

*v*stands for the flow velocity. Rewriting Eq. (8) in the standard form gives

*T*) and bearing (

_{h}*T*) torques [27,28] have expectedly shown the high sensitivity to the pressure drop obtained via Eq. (9) leading us to rewrite them as follows:

_{b}where *D _{s}* stands for the stem diameter of the valve,

*μ*indicates the friction coefficient of the bearing area, $Ci=(\pi /8)\mu Dvi2Ds$, and

*T*and

_{ci}*C*are the hydrodynamic torque and the sum of upper and lower contraction coefficients, respectively, depending on the valve rotation angle [1].

_{cci}*c*and

_{vi}*R*in order to model the system analytically. For our case study of $Dv1$ = 8 in. and $Dv2$ = 5 in., the valves' coefficients and resistances are developed as follows:

_{ni}*P*

_{1}and

*P*

_{2}terms are derived by combining Eqs. (6) and (18) as follows:

*P*

_{1}and

*P*

_{2}to $Rn1,\u2009Rn2,\u2009RL1,\u2009RL2$, and

*R*

_{con}are distinguishable through Eqs. (19) and (20), as observed in the practice. Any slight dynamic changes of the upstream set of the valve–actuator would be expected to be observed for the downstream one. The hydrodynamic and bearing torques' dependencies on all the resistances are reformulated as follows:

*f*is a nonlinear function of the changing

_{i}*T*,

_{ci}*C*, and the valve rotation angle. To carry out a systematic dynamic analysis, the following functions are fitted to the $Dvi3fi$ of each valve [5,29]:

_{cci}where $a1=0.4249,\u2009a1\u2032=0.1022,\u2009b1=\u221218.52,\u2009b1\u2032=\u221217.0795,\u2009c1=\u22127.823\xd710\u22124,\u2009c1\u2032=\u22122\xd710\u22124,\u2009d1=\u22121.084$, and $d1\u2032=\u22121.0973$.

where *b _{d}* indicates the equivalent torsional damping,

*K*is the equivalent torsional stiffness,

_{t}*V*stands for the supply voltage,

*x*is the plunger displacement,

*r*indicates the radius of the pinion,

*C*

_{1}and

*C*

_{2}are the reluctances of the magnetic path without air gap and that of the air gap, respectively,

*F*is the motive force,

_{m}*N*stands for the number of coils,

*i*indicates the applied current,

*g*is the nominal airgap,

_{m}*J*indicates the polar moment of inertia of the valve's disk, and

*R*is the electrical resistance of coil.

## Linear Dynamic Analysis

*b*) and the friction coefficient of bearing area (

_{di}*μ*) of both the sets to evaluate their effects on the stability/instability of the coupled sets

_{i}where $j22=\u221296bd1\u2212177.13\mu 1$ and $j55=\u221296.15bd2\u2212323.71\mu 2$.

*b*and

_{di}*μ*:

_{i}Using the numerical approach, the coupled sets' normalized eigenvalues are presented in Fig. 3 revealing an interesting stability map by assuming that *b _{di}*'s and

*μ*'s change equally for a critical range of $10\u22128\u2264bdi=\mu i\u22643\xd710\u22121$.

_{i}Figure 3 shows the instability and stability of the coupled sets for the ranges of $10\u22128\u2264bdi=\mu i\u22649\xd710\u22122$ and $10\u22121\u2264bdi=\mu i\u22643\xd710\u22121$ by presenting positive and negative real parts of eigenvalues, respectively. Such a stability map would help us select some critical values in order to capture practically observed chaotic and hyperchaotic dynamics.

## Results

Chaotic motions are known to be sensitive to even slight changes of initial conditions and hence examining different initial conditions potentially serve to characterize the responses obtained. The first set, Initial_{1}, would be a realistic option for the so-called “modulating” valves to regulate/reroute flow for many applications addressed earlier. The second set, Initial_{2}, is chosen to be close enough to the system's physically feasible equilibrium point and also to avoid numerical singularities.

Figures 4(a) and 4(b) show the phase portraits of the coupled valves for two sets of the initial conditions. Figures 4(a) and 4(b) reveal coupled chaotic and hyperchaotic dynamics. The hyperchaotic attractors by having two or more positive lyapunov exponents [33] are also known to be sensitive to initial conditions, and subsequently, orbits initiated from two close points move expectedly away from each other until the separation reaches the size of attractor. Some power tools of the nonlinear dynamic analysis would potentially help us characterize the responses obtained, in particular, Lyapunov exponents and Poincaré maps shown in Figs. 5(a)–7(b). A chaotic attractor presents one positive Lyapunov exponent [33], as shown in Fig. 5(a) ($L6=+0.1535$), indicating the chaotic motions of the interconnected valves–actuators configuration. Figure 5(b) presents two positive Lyapunov exponents not only for the sudden pipe contraction ($\theta =90\u2009deg$) but also for a broad spectrum of the approach angles ($35\u2009deg\u2264\theta \u226485\u2009deg$), which can be served as a proof of the network's hyperchaotic dynamics. Poincaré map is another power tool to distinguish among periodic, quasi-periodic, and chaotic responses. Note that for an *n*-dimensional system ($n\u22653$), this tool may not yield a clear nature of the response to determine whether the motion is chaotic or two-period quasi-periodic [33]. Although the Lyapunov exponents, as discussed earlier, would firmly confirm the chaotic and hyperchaotic motions of the interconnected actuated valves along with irregular Poincaré maps (almost different sets of 635 points) as shown in Figs. 6(a)–7(b) for each set of the upstream and downstream valves and for two sets of the initial conditions.

Figures 8(a) and 8(b) show the total flow loads, the sum of both the hydrodynamic and bearing torques versus the motive forces for two initial conditions. The squared areas remarkably magnify the differences between the chaotic and hyperchaotic responses of the coupled sets by presenting relatively larger attractor sizes for the hyperchaotic ones.

It is also of a great interest to present the transition from chaos to hyperchaos, in particular for an interconnected system. Figures 9, 10(a), and 10(b) reveal such a transition. Figure 9 shows a broad spectrum of Lyapunov exponents versus the equivalent *b*_{eqi}'s and *μ _{i}*'s (as discussed in Sec. 3), which explicitly separates the chaotic and hyperchaotic domains. The chaotic domain, for the range of $2\xd710\u22124\u2264bdi=\mu i\u226410\u22123$ by presenting one positive Lyapunov exponent, is transmitted to the hyperchaotic one for the range of $10\u22127\u2264bdi=\mu i<2\xd710\u22124$, which presents two positive Lyapunov exponents. As discussed earlier, the attraction domain of a hyperchaotic response is relatively larger than that of the chaotic one, as can be observed in Figs. 10(a) and 10(b). Figures 10(a) and 10(b) present bifurcation diagrams versus the equivalent

*b*

_{eqi}'s and

*μ*'s for the upstream and downstream valves, respectively. The larger domains of attractions of the hyperchaotic responses, for both the valves, are distinguishable with respect to the chaotic ones. The hyperchaotic and chaotic domains of the downstream set are expectedly larger than those of the upstream set.

_{i}One of the crucial issues which needs to be investigated is the effects of dangerous behavior of a valve–actuator set on another one. Figures 11(a) and 11(b) present this interesting situation in which the upstream set is assumed to be chaotic by exposing to the critical values of $bd1=\mu 1=10\u22127$, and the downstream set is operated safely using $bd2=\mu 2=10\u22121$, for the second set of initial condition.

The hyperchaotic motion of the upstream valve is again expected to be observed, but another smaller chaotic attractor, with a positive Lyapunov exponent of $L1=+0.013$, surprisingly emerges (Fig. 11(a)) which is significantly different from the previous case of Fig. 4(b). Its Poincaré map shown in Fig. 12(a) confirms the chaotic dynamics of the upstream set containing irregular points but expectedly different from the map presented in Fig. 7(a). It is of great interest to observe a very weak chaotic motion of the downstream valve as shown in Fig. 11(b), whereas a stable response was logically expected to be seen. Figure 12(b) is the Poincaré map of the downstream set revealing irregular points but too close to its equilibrium point. It is fairly straightforward to conclude that the chaotic motion of the upstream set is transmitted to the downstream one through the media trapped between them and subsequently affects its dynamics. Increasing the chaotic attractor domain of a set would accordingly magnify the domains of neighbor ones which would gradually cause failure of the whole network of thousands of valves–actuators. Such failures have to be avoided to reduce the considerable cost needed to restore the flow line.

## Conclusions

This paper represented an interconnected nonlinear model of two actuators and valves subject to the sudden contraction. These dependencies among different components were formalized to yield a sixth-order dynamic model of the whole system. We established a stability map to yield a clear picture of stability/instability of the coupled network for some critical values of the equivalent viscous damping and the friction coefficient of the bearing area.

The coupled chaotic and hyperchaotic dynamics were captured and discussed. Some powerful tools of the nonlinear dynamic analysis were then employed, including Lyapunov exponents and Poincaré map, to characterize the responses obtained. We presented the expected larger hyperchaotic attractor domains in comparison with the chaotic ones. One and two positive Lyapunov exponents were shown to confirm the chaotic and hyperchaotic dynamics of the coupled actuated valves, respectively. The irregular Poincaré maps were also presented to support both the chaotic and hyperchaotic dynamics along with the positive Lyapunov exponents.

The upstream valve–actuator set was intentionally operated with the same initial condition and critical values of the hyperchaotic dynamics to evaluate its effects on a stable downstream set. The dynamics of the upstream set was surprisingly different by revealing a chaotic attractor demonstrated by its Poincaré map and a positive Lyapunov exponent. The downstream set was also affected by the chaotic dynamics of the upstream one by showing the irregular Poincaré map but too close to its equilibrium point.

We are currently focusing our efforts on developing a comprehensive model for *n* valves and actuators to be operated optimally and safely in series.

## Acknowledgment

The experimental work of this research was supported by the Office of Naval Research Grant No. N00014/2008/1/0435. We appreciate this grant and the advice and direction provided by Mr. Anthony Seman III, the program manager.