Adaptive process control for achieving consistent particles' states in atmospheric plasma spray process

The coatings produced by an atmospheric plasma spray process (APSP) must be of uniform quality. However, the complexity of the process and the random introduction of noise variables such as fluctuations in the powder injection rate and the arc voltage make it difficult to control the coating quality that has been shown to depend upon mean values of powder particles’ temperature and speed, collectively called mean particles’ states (MPSs), just before they impact the substrate. Here, we use a science-based methodology to develop a stable and adaptive controller for achieving consistent MPSs and thereby decrease the manufacturing cost. We first identify inputs into the APSP that significantly affect the MPSs and then formulate a relationship between these two quantities. When the MPSs deviate from their desired values, the adaptive controller is shown to successfully adjust the input parameters to correct them. The performance of the controller is tested via numerical experiments using the software, LAVA-P, that has been shown to well simulate the APSP.


Introduction
The atmospheric plasma spray process (APSP) used to produce a variety of coatings such as low porosity, thermal barrier, wear and corrosion resistant, and functionally graded [1] is schematically shown in Fig. 1. Coatings of thickness less than a few micrometers are categorized as thin, and others as thick. Besides an APSP, coating techniques include electroplating, chemical treatments, hot dipping, chemical vapor deposition, physical vapor deposition (e.g., evaporation, puttering, ion plating), and pulsed laser deposition. Primary characteristics of coating processes are summarized in Table 2.1 of Fauchais et al. 's [2] chapter. The methodology of finding significant input parameters and controlling them to produce a good quality coating described in this paper is applicable to the above-mentioned processes and to other technologies like hot forging and machining. Details are provided below only for an APSP.
An APSP involves plasma formation within a gun nozzle, particle injection into the plasma jet flowing out of the nozzle, particle heating and acceleration during their motion in the plasma, and coating formation on a substrate. Usually, a mixture of argon (Ar, ~ 30-60 slm (standard liter per minute)) and hydrogen (H 2 , ~ 0-15 slm) injected into the gun nozzle is ignited by a high-intensity arc between tips of a cathode and a water-cooled anode. The arc intensity depends on the current (~ 300-600 A), the voltage (~ 30-70 V), the distance between the anode and the cathode, their geometries, and the gun efficiency.
The resulting plasma jet emanating from the nozzle outlet has the maximum temperature of ~ 14,000 K, and peak velocity of ~ 2.5 km/s at its centerline, and steep temperature and velocity gradients in the radial direction (~ 10 4 K/ mm and 5 × 10 5 /s) [3,4]. In general, micron-sized powder particles (~ 10-100 µm) and a carrier gas (typically Ar) are injected vertically downward through the powder port to reach the hot core of the plasma jet. These particles are heated, partially or fully melted, accelerated towards the substrate situated at ~ 10-15 cm from the nozzle, and are deposited on it to form a coating. Complex interactions between the plasma and the powder particles significantly influence coating properties and repeatability, and hence the process cost [5][6][7].
A small disturbance in process parameters can significantly alter mean particles' temperature and axial velocity (collectively called mean particles' states, MPSs) and thus the coating quality [8]. These disturbances can be caused by nozzle wear, powder port (injector) wear, pulsing of powder due to leaks, worn parts in the powder feeder, and powder dampness. The powder port's wear can be caused by angular-shaped particles interacting with its inner surface. Fluctuations in particles' mass flow rate can occur due to the powder clogging and/or changes in the carrier gas flow rate. Whereas the particle dwell time in the plasma is ~ 1-4 ms, the arc root fluctuations have time scales of the order of 10 -1 -10 -2 ms [7,9], time constants of the powder feeder instabilities are in the range of 1-200 s, and the electrode/injector erosion occurs in minutes. Hence, mitigating their effects on MPSs in a few ms by suitably varying significant input parameters using a real-time feedback controller can maintain consistency in the coating quality.
The development of a controller requires how input parameters (e.g., current, gas flow rates) affect MPSs.
Choudhury et al. [10] and Zhang et al. [11] used an artificial neural network (ANN) to establish such relationships. Even though an ANN can reveal complex nonlinear interactions amongst dependent and independent variables, it is a 'black box' that is prone to overfitting and is computationally expensive [12]. Some of these limitations can be overcome by using statistical design techniques such as screening, full factorial, and the response surface methodology. One can then focus on delineating effects of the main process parameters and cross-interactions amongst them on the output parameters; e.g., see Mawdsley [13] for the 2-level factorial design based multiple regression, and Li et al. [14] for the uniform design of experiments to develop regression equations. These investigations did not consider noise variables.
Srinivasan et al. [15] explored controlling the MPSs by using a monitoring system with a feedback controller to vary the primary gas flow rate and the current. Kanta et al. [16] used two ANNs-one to predict the process input parameters and the other to predict in-flight MPSs, and a fuzzy logic (FL) controller to adjust the former for keeping MPSs within their prescribed bounds. Planche et al. [17] developed an automatic system using an ANN to predict MPSs for the desired coating properties, another ANN to predict the needed input process parameters, and two FL controllers to correct errors in predictions from the ANNs. Finke et al. [7] developed real-time diagnostics and control by monitoring the MPSs and adjusting the input parameters with a proportional-integral-derivative controller that employs a manually selected constant gain factor. These works did not consider effects of random noise variables.
The motivation, then, is to develop a stable and adaptive controller for an APSP to achieve the desired MPSs. The adaptive controller presented here has an adjustable tuning mechanism to cope up with effects of unknown Fig. 1 Schematics of an APSP [3]. Gases injected into the gas gun pass over the arc to form a plasma that exits the gun. Powder particles through the powder port are injected into the plasma and travel with it to the substrate. The particles are heated during their trajectory in the plasma. The plasma exiting the gun outlet and trajectories of particles injected into it are numerically simulated to find particles' temperature and velocity as they move through the plasma disturbances and has been designed using a model reference adaptive control (MRAC) scheme based on the linearized dynamics of the MPSs.

Development of mathematical model of the problem
A major challenge in developing a science-based model of a physical problem is to capture physics relevant to all important phenomena such as the plasma jet exiting from the nozzle outlet, chemical interactions amongst the jet constituents, the drag force experienced by powder particles injected into the plasma, heat transfer between particles and the plasma, and plasma's interaction with the air surrounding it. The mathematical model of the APSP developed in this paper is based on the following assumptions. (1) The plasma jet exiting the nozzle is axisymmetric, is in local thermodynamic equilibrium, and is composed of a mixture of chemically reacting constituents with temperature-dependent thermodynamic and transport properties.
(2) Chemical reactions among different species including ionization, dissociation, and recombination are considered, and the turbulence is simulated by using the standard k− ∈ model. (3) The arc voltage fluctuations in the plasma nozzle are neglected because of difficulties in expressing them in terms of gases injected into the nozzle, the anode and the cathode materials, and the distance between them. (4) Powder particles are rigid spheres, randomly vary in diameter, do not interact with each other, exchange heat with the plasma, and can melt due to temperature rise, and the internal convection within a molten particle has a negligible effect on the heat transfer. (5) A powder particle is regarded as a point mass for simulating its heating, melting, and solidification [6]. The volume fraction of the melted material is computed by considering the latent heat of melting. (6) Effects of gravity and buoyancy on the plasma jet and particles' trajectories are negligible as compared to that of the viscous drag forces [18,19]. (7) The effect of the carrier gas flowing through the powder port at ~ 5 slm on the plasma jet is negligible [6].

Justification of assumptions
Assumption 5 considerably simplifies the analysis of heating of particles, has been shown by other investigators to give fairly accurate results that are close to the experimental observations, and saves considerable computational resources. Assumption 6 is reasonable because of large differences between particles' and the plasma velocities, and drag forces are dominant. For assumption 7, we note that even though the injection of the carrier gas and the powder particles adversely affects axisymmetry of the plasma jet near the powder injection port, however, the effect is local and the plasma jet is nearly axisymmetric away from the powder port. The injection velocity (speed and direction) of particles varying in diameters from 5 to 125 μm are randomly assigned within their pre-specified ranges to achieve the desired powder mass injection rate.

Governing equations and solution techniques
For brevity, we refer the reader to Sects. 3, 4 and 5 of Shang et al. [3] for details of the governing equations for the plasma modeled as a mixture of different constituents, and particles' trajectories and heating. The software LAVA-P numerically solves these equations by the finite volume method over the simulation region exhibited in Fig. 1 and using the following boundary conditions (BCs) [20] for the temperature, the velocity, and the turbulence variables at the nozzle exit plane.
Here, T f is the plasma temperature; u and v, respectively, are the radial and the axial velocity components; R i and R o , respectively, are the inner and the outer radii of the nozzle exit; T 0 and V 0 , respectively, are the plasma temperature and the axial velocity at the nozzle exit; T w is the torch wall temperature; and k and , respectively, are the turbulence kinetic energy and its dissipation rate. Coefficients n T and n v indicate how fast the temperature and the axial velocity decrease in the radial direction; the x-axis is in the radial direction, ( u∕ x) max equals the maximum value of the radial velocity gradient; the parameter κ determines the shape of the profile, and δ 0.1 equals the jet width at the location where v = 0.1V 0 . The torch wall is assumed to be rigid where no-slip condition holds [20]. T ∞ represents the ambient temperature. Open boundaries and downstream boundaries are considered as either inward or outward flow conditions [20]. On the axis of symmetry, x = 0, BCs are k(x, 0) = 0 (x, 0) = 0 u = 0 and u y + v x = 0. Initially, the simulation region is filled with quiescent air at room temperature.
Following [20], we choose n T and n v , and iteratively find T 0 and V 0 so that the total mass flow rate and the energy output from the nozzle exit differ by less than 1%, respectively, from their input values into the nozzle.
For the specified powder mass injection rate, the software calculates the number of particles for a given average injection speed and assigns them a random velocity when exiting the powder port.

Parameter values
Values assigned to parameters are: R i = 0.39 mm, R o = 15 mm, the gun efficiency = 0.7, T ∞ = 300 K, P ∞ = 85.5 kPa, T w = 700 K, κ = 0.015. Values of material parameters used for the zirconia (Zr) powder particles are listed in [6]. Figure s1 of the supporting information represents the distribution of particles' diameters used in simulations, a schematic of their trajectories in the plasma, and the location of the observation window used to measure particles' temperature and velocity.

Measure of quality of computed results
The deviation of the computed results from those measured experimentally is quantified by using the following L 2 -norm, e 2 , of the difference between them.
where represents either the temperature or a velocity component. This provides the average of deviations between the computed and the experimental values of rather than the maximum difference between them at any location. A local spike in the difference over a small distance has the same effect on e 2 as a small difference over a large distance.

Finding optimal input parameters for the MPSs
Three steps involved in finding values of the process input parameters for achieving desired MPSs are screening, finding response functions relating significant input parameters to the MPSs, and using a genetic algorithm (GA) to solve nonlinear algebraic response function equations for the input parameters when the MPSs are prescribed. The screening analysis ranks the influence of each process parameter on the MPSs to identify which ones should be included in developing response functions for the MPSs. (2) To investigate statistical significance of process parameters and interactions between any two of them, a complete polynomial of at least degree 2 is needed. Here, we use a complete polynomial of degree 2, employ the analysis of variance (ANOVA) to statistically quantify the significance of various terms in the polynomial, construct an objective function in terms of the two response functions (one for the mean particles' temperature and the other for the velocity), and use an GA to find values of the significant process parameters from the response functions for desired MPSs.

Screening of process parameters
We use Morris' [21] global one-factor-at-a-time (OFAT) screening method to rank their effects on the MPSs as: (1) negligible, (2) linear and additive, and (3) nonlinear and/ or involved in interactions with other process parameters [21]. As explained in [21,22], it requires calculating elementary effects (EEs) of a randomly chosen value of one process parameter at a time on the MPSs. We assume that each MPS is a scalar-valued function of 7 process parameters u i (i = 1, … , 7) that have been normalized to vary between 0 and 1. Each , the EE of the i th parameter at a randomly chosen base point u * = u 1 , … u 7 is defined by The value of EE i u * equals the sensitivity of the MPS ( denoted by y in Eq.3) to the ith parameter at the base point u * . The mean, i , and the standard deviation, i , of the EE i of the ith process parameter are computed using Here, u (j) (j = 1, … , r) are r randomly chosen base points, and EE (j) i represents the EE of the i th parameter at the base point u (j) . Values of i and i identify process parameters that have negligible (low i ), linear (high i and low i ), and nonlinear and/or interactions (high i and high i ) effects [21,23].

Identification of significant input parameters using analysis of variance (ANOVA)
The results of ANOVA are meaningful only if the input values used in the design of experiments are independent of each other [24], i.e., the correlation coefficient, cr , between the input variables ξ 1 and ξ 2 defined by is very small. Here, superscript (i) refers to the ith observation out of n observations, and a bar superimposed on a variable indicates the mean value of n observations. The cr has the range [-1, 1], its magnitude indicate the strength and sign of the direction of the linear dependence between the two variables [24]. The -ve sign means that variables 1 and 2 have opposite effects.
To investigate effects of the main and the cross-interactions between process parameters on MPSs, we use ANOVA by considering a complete polynomial degree of 2 given by Eq. (6) to fit the computed data to the MPSs.
In Eq. (6), y model represents the normalized fit to the MPSs, u = P, Q, I, d y , PQ, PI, Pd y , QI, Qd y , Id y , P 2 , Q 2 , I 2 , d 2 y and the vector b represents the sensitivity coefficients that quantify the effect of an element of u on the output [24]. To delineate the importance of the parameter u i 0 in the array u , we define a new model that omits the term u i 0 from the sum in Eq. (6). The sum of squares (SS) associated with u i 0 is given by Here, N equals the number of samples and y the mean of the vector y . The mean squares are defined as where DF i 0 is the degree of freedom of the parameter u i0 . The value of MS i 0 is used to test the null hypothesis by calculating the F-ratio defined as where MS E is the mean squares associated with the residual error. If F i 0 is less than the standard value of its F -distribution, then the parameter u i 0 is significant. The significance of terms in Eq. (6) is calculated by using the null hypothesis with 0.5% rejection, i.e., if the p-value (Probability > F) for b i is less than 0.005, then u i will be significant [25] in evaluating y model .

Finding the response functions, i.e., expressions for MPSs in terms of significant process parameters
The response functions or equivalently expressions for the MPSs in terms of significant process parameters are found by performing numerical experiments by using the full factorial approach. For the number of significant processing parameters = a and b levels for each parameter, b a simulations are conducted. The coefficients of the response functions in Eq. (6) are estimated by using the linear regression analysis, and finding R 2 where The value of R 2 close to 1 means that the model accurately captures the variability of the MPSs with the change in process parameters.
We note that Batra and Taetragool [26] used the response surface methodology and the regression analysis to relate MPSs to the significant input parameters.

Solution of nonlinear algebraic equations using the genetic algorithm
The GA explained in [27] is used to solve, within a prescribed error, the nonlinear Eq. (6)for the significant input parameters for desired values of the MPSs. We use the 'ga' toolbox in MATLAB with default values of parameters to minimize the Error defined by where the subscript 'des' represents the desired value of the variable, and P , Q and I , respectively, denote the Ar flow rate, the H 2 flow rate, and the current. We accept the solution that reduces the normalized error to less than 10 -6 .

Development of adaptive controller
The objective is to develop an adaptive controller for the APSP to maintain MPSs within the prescribed range when the average powder injection velocity and the arc voltage change. The design criterion is to attain the desired MPSs in less than 50 ms for 20slm ≤ P ≤ 60slm, 0 ≤ Q ≤ 20slm, and 300A ≤ I ≤ 600A.
Developing the adaptive controller involves the following five steps: a mathematical model, the system identification, controller design, implementation, and testing. These are briefly explained below.

Approximate mathematical model for the MPSs
We find below a relation between the inputs and the outputs by linearizing the nonlinear dynamics of particles' states around a known steady-state solution called an equilibrium point. Recall that a particle's temperature is found by assuming it as a rigid sphere [3], and its motion in the plasma is affected by only the drag force exerted on it by the plasma. The particle's acceleration, governed by Newton's second law of motion, equals the drag force divided by its mass. That is, where U p = u p v p w p T is the particle velocity vector, U f the plasma velocity, means mass density of the plasma, p the particle mass density, r p the particle radius, U rel the velocity of the particle relative to that of the plasma, U rel = U f − U p , and C D the drag coefficient. The heating of a particle is governed by where V p is the particle volume, T f the plasma temperature, T p the particle temperature, c p the specific heat capacity, A p the surface area of the particle, and h the heat transfer coefficient. Using A p = 4 r 2 and V p = 4 3 r 3 (r equals the particle radius) in Eqs. 11 and 12, the particle acceleration in the axial (y-) direction and its temperature are given by where v p (t) and v f (t) are, respectively, components of U p (t) and U f (t) in the y-direction. Equations 13 and 14 are firstorder nonlinear ordinary differential equations with v f and T f as excitations.
To develop an approximate model for the MPSs as a function of the Ar and the H 2 flow rates and the current I, we assume the following affine relations: Substitution from Eq. 15 into Eqs. (13) and (14) gives Here, | | U s rel | | equals its value in the steady-state configuration around which the equations are linearized. Neglecting the bias terms c 0,v and c 0,T in Eq. (16) results in the following multi-input multi-output (MIMO; 3 inputs and 2 outputs) state-space (SS) model for the MPSs.
Values of constants a v , a T , b 11 , ..., b 23 in Eq. 17 depend on the equilibrium point, value of | | U s rel | | at the equilibrium point, and are estimated using the input-output data.
We will verify in Results section that Eq. 17 well represents deviations from the equilibrium point of the MPSs.

System identification
For estimating constants in Eq. 17, we use the following sinusoidal variations in the input variable u = P, Q or I: where u b is the base value, u a1 , u a2 , u a3 , u a4 , and u a5 are amplitudes of perturbations; 1 , 2 , 3 , 4 , and 5 are frequencies, and t is the time in ms. According to [28], the input variations given by Eq. 18 are sufficiently rich as they contain enough frequencies of excitation. Values of the 11 variables u b , u a1 , u a2 , … , u a5 , 1 , 2 , … 5 are chosen using the Latin Hypercube sampling approach based on the normal distributions with the mean ( ) and the standard deviation ( ). The corresponding MPSs are calculated in the observation window at a sampling time of 0.01 ms. The values of parameters in Eq. 17 are estimated by using the simulated (with LAVA-P) input-output data and the MATLAB toolbox 'ident' lucidly explained by Ljung [29] the prediction error method (PEM) framework, and the MAT-LAB function 'ssest' with the order of the system = 1, input delay = 0, form = canonical, feed through = 0, disturbance model = none, initial states = 0, and the weighting scheme for singular value decomposition = the canonical variable algorithm. The entire input-MPS data of each sample are first processed by subtracting from it the mean value. Next, the output data are smoothened by using a moving average of 100 data points to reduce noise in it. The resulting input-output data are divided into two equal parts; the first part is used for estimation of parameter values and the second part for validation. The quality of model prediction is measured by the comparison index, 'Fit, ' defined as where ŷ is the predicted output for the input y , y is the mean of the output, and ‖⋅‖ indicates the L 2 norm (or the length) of a vector.

Controller design
The linearized MIMO SS model, Eq. 17, for the MPSs around an equilibrium point is rewritten as Here, y(t) = {v(t), T (t)} T is the output vec tor, u(t) = {P(t), Q(t), I(t)} T the input vector, and matrices 'A' and 'B' depend on the equilibrium point (i.e., the operating point). We assume that matrix A is Herwitz and the plant is controllable, i.e., for any desired output y des , there exists a sequence of inputs u(0), u(1), … , u(k) such that y(k) = y des within an acceptable tolerance. The model reference adaptive control (MRAC) algorithm involves two loops-the inner loop consisting of a feedback and the outer loop of an adaptive mechanism for the controller as shown in Fig. 2. Based on Eq. 20, we choose the following reference model for updating the controller gains Using MPSs at the instant of a perturbation initiation as the initial conditions, the outputs v m and T m of the reference model asymptotically converge, respectively, to the desired outputs v des and T des with the rate of convergence depending on the parameters a mv and a mT . In order for the MPSs v(t) and T (t) from LAVA-P to asymptotically track the reference outputs v m (t) and T m (t) , we consider the following control law: The objective is to design an adaptive law that guarantees (a) the stability of the control system (controller + process dynamics), and (b) the MPS y(t) tracks y m (t) within a small bound for all finite inputs u(t) . The theory for designing the adaptive controller is given in the supporting information. As explained there, for positive definite matrices Λ and P, we find gain matrices K (t) and L(t) by using the following adaptive law:

Implementation
By assuming that the LAVA-P represents an actual plant, we discuss implementation of the adaptive controller. We assume that the MPSs reach their desired values in 50 ms, and set for the reference model. For the tracking error to be asymptotically stable, there should exist a positive definite and symmetric matrix R such that the time derivative of the Lyapunov function (refer to the supporting material), V = −e T Re ≤ 0 [28]. A different choice of R affects the transient response but not the boundedness and the asymptotic convergence of the error. We choose where u sol = P sol Q sol I sol T is the optimal solution obtained from the GA for the desired output The next step is to choose an appropriate adaptive gain matrix Λ for the APSP. The theorem presented in the supporting material only requires that Λ be positive definite. Lacking any systematic way to find Λ, the following form is found by trial and error: where 1 , 2 , and 3 are tuning parameters for P, Q, and I, respectively. By tuning these parameters, one can adjust the MPSs in real-time remembering that a large value of can make the MRAC unstable. Figure 3 shows the proposed intelligent adaptive process control scheme for the APSP to get consistent MPSs. At the start of the process, values of P, Q, and I are predicted using the response function, the desired values of MPSs, and the GA. These are used to find initial estimates of the controller gains from Eq. 24 which are then used to calculate inputs from Eq. 22 to be used in LAVA-P. During the spray process, the MPSs are measured in the observation window from LAVA-P at every sampling time, and the inputs are adaptively varied by the MRAC that updates the gains to minimize the tracking error using Eqs. 21 and 23 and achieve desired MPSs. The controller stops the process (or the plant) if any of the inputs reach their limiting values.

Testing
The proposed adaptive control scheme is tested with software LAVA-P for disturbances introduced either in the average injection velocity or in the arc voltage. The computed and the experimental distributions of the steady-state plasma velocity and temperature along the jet axis (x = 0) and on the radial line y = 2 cm are compared in Fig. 4. With time t = 0, when the plasma jet exits the nozzle, the plasma was found to reach a steady state at 4 ms with at most 1% error (pointwise) in the axial velocity and the temperature in the computational domain.

Particles' axial velocity and temperature
Powder particles are injected one at a time into the plasma through the powder port located at the point (x = 8 mm, y = 6 mm) at t = 5 ms, and their trajectories till t = 10 ms are recorded. The computed average values of particles' axial velocity and temperature in the 1 cm wide window (9.5 ≤ y ≤ 10.5 cm, shown in Fig. s1 of the supporting information) at t = 9.1, 9.2 …, 9.9, 10 ms are compared in Fig. 5   Fig. 4 Comparison between the simulated and the experimental results of Williamson et al. [30]: a plasma axial velocity (e 2 = 14.82%) and plasma temperature (e 2 = 5.59%) distributions along the jet (y-) axis, b plasma axial velocity (e 2 = 5.24%) and temperature (e 2 = 11.9%) versus the distance from the jet axis on the radial line y = 2 cm. The decaying trends along and perpendicular to the jet axis of the computed and the experimental temperatures and the axial velocities agree qualitatively and quantitatively with their experimental values of Smith et al. [31]. Smith et al. [31] used a two-color pyrometer, IPP2000, and a spray pattern trajectory (SPT) sensor to measure, respectively, the average temperature and the axial velocity of particles flowing through the observation window. They emphasized that the SPT sensor was focused on a narrow zone near the plasma jet axis but did not list coordinates of the sensor location. To compare the computed results with the experimental findings of Smith et al., they are divided into three groups based on their diameters in the ranges 30-42 µm, 56-71 µm, and 85-100 µm. The variations with the axial distance y from the nozzle of the average axial velocity and the average temperature of particles at time t = 10 ms are shown in Fig. 5a, b. The computed particle characteristics of small particles are close to those measured experimentally. However, computed results for the medium and the large size particles differ from the corresponding experimental findings. At some axial locations, the difference in the mean axial velocities of particles in each of the two groups is about 50%. The small-sized particles have higher axial velocity because they are easily accelerated by the plasma jet. The results depicted in Fig. 5b reveal that the computed average temperature of particles agrees well with that measured experimentally. The temperatures of the small-size particles are higher than that of large particles due to their in-flight motion near the hotter jet axis. The smaller size particles have melted by the time they go out of the observation window and arrive at the substrate. By observing the plot (not included here) between temperature vs size of injected particles, we note that the temperature of the mediumsize particles (45-65 μm) stays close to the melting point (2,950 K) during their in-flight path. Many of the large particles are only partially melted or even not melted at all. Recalling that the heating model (cf. Equation 12) regards a particle as a point mass, a large particle requires more time to be heated to the melting temperature. Larger particles have low axial velocity resulting in more residence time in the plasma jet which helps in their being heated up. Also, some large particles travel farther away from the jet axis where the plasma temperature is not very high and hence are either partially melted or not melted at all depending upon their locations in the plasma jet.

Input parameters that significantly influence the MPSs
Keeping the applied voltage, the average powder injection speed and the particle size fixed at 50 V, 10 m/s, and 30-100 μm, respectively, we present results of screening the seven process parameters, namely the current, the Ar flow rate, the H 2 flow rate, the powder feed rate, the spray distance, and the injector location (y, x) along and perpendicular to the jet axis. Their minimum and maximum values considered are: current (A) 300, 600; H 2 flow rate (slm) 2, 12; Ar flow rate (slm) 30, 60; powder feed rate (g/min) 10, 40; stand-off distance (cm) 7.5, 12.5; y-location of the injector (cm) 0.4, 1.0, and x-location of the injector (cm) 0.6, 1.0.

Screening using elementary effects
We choose r = 8 in Eq. (4) and investigate the influence of the seven parameters for p = 12 and p = 20 resulting in 8 × (7 + 1) = 64 simulations needed for this study. The 130 μm are injected into the plasma at t = 5 ms to achieve a mass rate of 20 g/min. The particle size distribution and the location of the observation window are exhibited in Fig. s1 of the supplement material mean and the standard deviations of EEs of the MPSs are depicted in Fig. 6. The results associated with p = 12 and p = 20 are qualitatively similar. The stand-off distance, the powder feed rate, and the x-location of the injector are the least significant parameters since their mean EEs have low values. The y-location of the injector, the current, and the Ar flow rate have nonlinear effects and/or interactions among them as signified by high values of the mean and the standard deviation for them. The H 2 flow rate has the least effect on the mean axial velocity but a strong effect on the mean temperature. Accordingly, the Ar flow rate, the current, the y-location of the injector, and the H 2 flow rate are considered as significant process parameters for developing the response functions.

Screening using ANOVA
The The results of the ANOVA test for the MPSs are given in Table 1. Based on the null hypothesis, we find that all process parameters in Eq. 6 for both MPSs are significant since their p-values are almost zero. Even though the p-value for the interaction between Q and d y is greater than 0.005, those for Q and d y are less than 0.005. From results of the ANOVA and the EEs, and for illustrating the development of the controller, we henceforth take the Ar flow rate, the current, the y-location of the injector, and the H 2 flow rate as significant process parameters for developing the response functions. Values of the probability are indicated in bold font to easily decide which parameters are significant.

Response functions for the mean particles' states
We now present results for the response functions for the MPSs from the numerical experiments designed by using the full factorial approach for the 4 significant process parameters and 5 levels for each parameter that require Values of cr (computed using the definition in Eq. 5) providing the correlation among all parameters (process parameters and MPSs) are listed in Table 2. The correlations between the process parameters are indeed zero; hence, they are independent of each other. Thus, the corresponding MPSs are statistically significant. The correlation coefficients suggest that the dependence of the MPSs upon the Ar flow rate is strong. The dependence of the mean axial velocity upon the current is minimal. The correlation between the mean temperature and the Ar flow rate is negative, which implies that there exists an inverse dependence between them. There is also inverse dependence between the MPSs and the injector location along the y-axis. The correlation coefficient between the mean Table 2 Correlation coefficients between the process parameters and the MPSs (Symbols P, Q, I, and d y are defined in Fig. 6)  axial velocity and the mean temperature is -0.07, which implies that the mean axial velocity does not depend upon the mean temperature.

Response functions
The coefficients of the quadratic expression in Eq. 6 are estimated by using the linear regression analysis on the data from the full factorial approach (625 simulations) resulting in the following response functions for the MPSs.
Here, subscript 'n' represents the normalized parameter with values between 0 and 1. The R 2 values of 1 and 0.98, respectively, for the mean axial velocity and the (27) v n, mod el = 0.248 + 0.476P n + 0.105Q n + 0.308I n − 0.276d y,n − 0.040P n Q n + 0.047P n I n − 0.034P n d y,n +0.091Q n I n − 0.013I n d y,n − 0.078P 2 n − 0.061Q 2 n − 0.097I 2 n + 0.048d 2 y,n ;R 2 = 1 (28) T n, mod el = 0.651 − 0.323P n + 0.331Q n + 0.298I n − 0.127d y,n + 0.122P n Q n + 0.152P n I n − 0.106P n d y,n − 0.108Q n I n −0.017Q n d y,n + 0.139I n d y,n − 0.044P 2 n − 0.162Q 2 n − 0.126I 2 n − 0.044d 2 y,n ;R 2 = 0.98 mean temperature imply that the postulated expression in Eq. 6 is a good representation of the numerical data. The validity of the response functions in Eqs. 27 and 28 is tested for a new set of process parameters, and the corresponding results are listed in Table s1 of the supporting information. For these cases, the percentage error is less than 2%.

Values of process parameters for desired MPSs obtained by solving Eqs. 27 and 28
The solutions for arbitrarily selected values for the desired MPSs computed by minimizing the error defined in Eq. 10 using the GA are summarized in Table 3. For a given gas gun, the powder inject port location is fixed. Hence, d y cannot be changed and we work with three significant input process parameters, namely the current, the Ar flow rate, and the H 2 flow rate. Recalling that the GA solution depends upon the seed number, results for the first seed number are provided in Table 3 and those for the second set of seed numbers in the supporting information. Results from the two sets of seed numbers are similar to each other. It is clear that the GA values of the three process parameters provide MPSs that differ by less than 4.5% from their desired values that occurs for the desired mean axial velocity = 100 m/s and the desired mean temperature = 3,100 K. We emphasize that for accurately predicting values of the three process parameters, the desired values of the MPSs should be in the range used to generate the response functions given by Eqs. 27 and 28. The desired mean particles' velocity of 3,100 m/s is outside the range employed to develop Eqs. 27 and 28.

Validation of three-input two-output model represented by Eq. 17
The 10 randomly generated functions u in Eq. 18 for exciting the Ar flow rate are listed in Table 4, and the corresponding variations of inputs are shown in Fig. 8a Fig. 8b,c are smoothened by taking the moving average of 100 trailing points, we get solid curves displayed in Fig. 8d. Similar analyses are conducted for randomly generated vectors for exciting the H 2 flow rate and the current. Their plots as well as those of results for the MPSs are provided in the supporting information.
Predictions from models for the MPSs agree well with those found using LAVA-P software with an average fit of 85% for training and 82% for validation for the mean axial velocity, and 72% for training and 66% for validation for the mean temperature.
The parameters for the MPSs estimated using results for the 50 samples, depicted as box plots in Fig. 9, are distributed close to each other with a few outliers enclosed in red circles. Thus, Eq. 17 satisfactorily represents the dynamics of the MPSs for variations in the input variables.

Results for the adaptive controller
Based on the knowledge gained from preliminary results, we select 1 = 2, 2 = 2 and 3 = 5 in Eq. 26 and investigate the performance of the proposed control scheme for the problem presented in subsection 3.5. The objective of the MRAC is to attain the desired mean particles' axial velocity = 90 m/s and temperature = 2,850 K within 50 ms of the introduction of a disturbance. Results depicted in Fig. 10 confirm that indeed the objective is achieved, and the proposed controller minimizes deviations in the MPSs caused by disturbances. Fig. 11 The closed-loop control responses of the MRAC for a increasing desired mean particles' axial velocity but decreasing mean temperature, and b increasing desired mean particles' axial velocity and mean particles' temperature disturbances. Of course, in practice disturbances are not limited to those stipulated here. Even though the controller will not work for all disturbances, one can use the methodology to design controllers for a pre-specified class of disturbances. The real-time performance of the proposed controller will depend upon how quickly a motor can change the process input parameters, and the time lag of the plant to respond to these changes.

Conclusions
This paper proposes an approach to design an adaptive controller for an atmospheric plasma spray process (APSP) to get consistency in mean particles' states (the axial velocity and the temperature, MPSs) before they impact the substrate. Throughout this work, the physical tests are replaced by numerical experiments that use the software, LAVA-P, whose predictions are first shown to agree well with the experimentally measured values of particles' velocity and temperature. From the screening analysis, the argon (Ar) and the hydrogen (H 2 ) flow rates, the current, and the y-location of the injector are identified as significant input parameters for the MPSs. The response functions (i.e., relations between the MPSs and the significant input parameters) are first established as nonlinear algebraic equations between the four significant input parameters and the MPSs. For a given gas gun, the powder injector location is generally fixed. For desired values of the MPSs, the trial input values of Ar and H 2 flow rates and the current are predicted by solving these nonlinear algebraic equations using a genetic algorithm (GA) included in MATLAB. Using equations for the evolution of particle's axial velocity and temperature linearized around steady-state values of these variables, an adaptive controller is designed to modify the trial input values of the three process parameters to obtain desired values of the MPSs in the presence of a noise variable (either a disturbance in the powder injection velocity or in the arc voltage). The performance of the controller has been successfully tested for two example problems.
Additional research is needed to design a robust controller for mitigating effects of all noise variables and either warn the operator or shut down the process when MPSs cannot be controlled. The methodology presented here is applicable to other coating processes and manufacturing processes such as hot forging and machining.