2.9 Feedback and Feedforward Control

Size: px
Start display at page:

Download "2.9 Feedback and Feedforward Control"

Transcription

1 2.9 Feedback and Feedforward Control M. F. HORDESKI (985) B. G. LIPTÁK (995) F. G. SHINSKEY (970, 2005) Feedback control is the action of oving a anipulated variable in response to a deviation or error e between the controlled variable c and its set point r in such a way as to reduce, or if possible, to eliinate the error. Feedback control cannot anticipate and thereby prevent errors, because it can only initiate its corrective action after an error has already developed. Because of dynaic lags and delays in the response of the controlled variable to anipulation, soe tie will elapse before the correction will take effect. During this interval, the deviation will tend to grow, before eventually diinishing. Feedback control cannot therefore achieve perfect results; its effectiveness is liited by the responsiveness of the process to anipulation. By contrast, feedforward correction can be initiated as soon as any change is detected in a load variable as it enters the process; if the feedforward odel is accurate and the load dynaics are favorable, the upset caused by the load change is canceled before it affects the controlled variable. Because feedforward odels and sensors are both iperfect, feedforward loops are usually corrected by feedback triing. FEEDBACK CONTROL Set, r Controller Load, q Valve or flowloop Controller var., c FIG. 2.9a Load changes can enter through different gain and dynaics fro the controller output. The purpose of any process control syste is to aintain the controlled variable at a desired value, the set point, in the face of disturbances. The control syste regulates the process by balancing the variable load(s) with equivalent changes in one or ore anipulated variables. For the controlled variable to reain stationary, the controlled process ust be in a balanced state. Regulation through feedback control is achieved by acting on the change in the controlled variable that was induced by a change in load. Deviations in the controlled variable are converted into changes in the anipulated variable and sent back to the process to restore the balance. Figure 2.9a shows the backward flow of inforation fro the output of the process back to its anipulated input. The load q is a flow of ass or energy entering or leaving a process, which ust be balanced by an equal flow leaving or entering. It ay have ore than one coponent for exaple in a teperature loop, both the flow and teperature of a strea are coponents of its energy deand, but they ay be balanced by a single anipulated variable such as stea flow. The steadystate gain K and dynaic gain vector g in the paths of the anipulated and load variables ay differ, and therefore have been given different subscripts and q. Liitations of Feedback Control Feedback, by its nature, is incapable of correcting a deviation in the controlled variable at the tie of detection. In any process, a finite delay exists between a changing of the anipulated variable and the effect of the change on the controlled variable. Perfect control is not even theoretically attainable because a deviation in the controlled variable ust appear before any corrective action can begin. In addition, the value of the anipulated variable needed to balance the load ust be sought by trial and error, with the feedback controller observing the effect of its output on the controlled variable. Best-Possible Feedback Control An absolute liitation to the perforance of a feedback control loop is iposed by any deadtie in the loop. Figure 2.9b describes the results of a step-load change applied to a process whose dynaics consist of deadtie and a single lag in both the load path g q and the anipulated-variable path g. The tie scale is given in units of deadtie τ d in the path of the anipulated variable. To siplify the illustration, the deadties in both paths are identical (this is not essential any K q g q Process K g 73

2 74 Control Theory c r e b E b Manual These values of e b and E b are the best that can be obtained for any feedback controller on a process whose dynaics consist of deadtie and a single lag. A real controller will yield a soewhat larger peak deviation and an area at least twice as great, depending on its odes and their tuning. If the process contains two or ore lags, the value of E b reains the sae as estiated above but is ore difficult to approach with a real controller. q Tie, t/t d FIG. 2.9b The best-possible load response for a process having deadtie and a single lag. deadtie in the load path, or none, will produce the sae response, siply shifting its location in tie). Also for siplification, the steady-state gains in both paths are ade equal. At tie, a step decrease in load q enters the process. After the deadtie in the load path expires, at tie zero, the controlled variable c responds, beginning to rise along an exponential curve deterined by the gain and lag in the load path. If the controller were left in Manual, it would continue to a new steady state. In this exaple, K q is 2.0, leaving the final value of c in Manual having changed twice as uch as load q; the tie constant τ q of the load lag in this exaple is 2.0τ d. Also shown is the possible return trajectory of c to set point r if the anipulated variable were to step as soon as any deviation is detected, i.e., at tie zero, by an aount = q( K / K )( + / ε τ τ ) q d q 2.9() where ε is the exponential operator This ove turns out to be the best that is possible by a feedback controller, as it causes the deviation to decay to zero during the next deadtie. The peak deviation reached during this excursion is e b d q = qk ( / ε τ τ ) q 2.9(2) At the tie the peak is reached, the controller output ust be stepped to atch the new steady-state load. The leading edge of the curve, i.e., up to tie.0, is deterined copletely by the load step and the gain and dynaics in its path; the trailing edge is deterined by the controller and its tuning. The leading and trailing edges of the best response are copleentary, so that the area that they enclose is d q E = e τ = qk τ ε τ / ( τ ) b b d q d 2.9(3) Integrated Error The actual effectiveness of feedback control depends on the dynaic gain of the controller, which is a function of its control odes and their tuning. Although a high controller gain is desirable, the dynaic gain of the closed loop at its period of oscillation ust be less than unity if the loop is to reain stable. In effect, then, the dynaic gain of the process dictates the allowable dynaic gain of the controller. For each process and controller, optiu settings exist that iniize soe objective function such as Integrated Absolute Error (IAE). However, ost controllers are not tuned optially, for various reasons, such as process nonlinearities. Therefore, controller perforance needs to be stated in ters of the actual settings being used. As an exaple, the integrated error produced by a load change to a process under ideal Proportional- Integral-Derivative (PID) control will be evaluated based only on its ode settings. The controller output at tie t is related to the deviation e between the controlled variable and its set point by 00 P e t edt D de = + + C I + t0 dt 2.9(4) where P, I, and D are the percent Proportional Band, Integral tie and Derivative tie, respectively, and C is the output of the controller at tie, t 0 when it was first placed in Autoatic. Let t be a steady state, so that e = 0 and its derivative is also zero. Then a load change arises, causing the controller to change its output to return the deviation to zero. When a new steady state is reached, the controller will have an output 2 : t 00 P e 2 edt D de 2 = + + C I t dt (5) where e 2 and its derivative are again zero. The difference in output values between the two steady states is t 00 2 = 2 = edt PI t Solving for the integrated error: t2 E = e dt = PI t (6) 2.9(7)

3 2.9 Feedback and Feedforward Control 75 For any sustained load change, there ust be a corresponding change in the controller output, and in fact the current load is often reported as the steady-state controller output. The wider the proportional band and the longer the integral tie, the greater will be the integrated error per unit load change. (While the derivative setting does not appear in the integrated-error function, its use allows a lower integral tie than is optiu for a PI controller.) Therefore, loops with controllers having large PI products are candidates for feedforward control. FEEDFORWARD CONTROL Feedforward provides a ore direct solution to the control proble than finding the correct value of the anipulated variable by trial and error, as occurs in feedback control. In the feedforward syste, the ajor coponents of load are entered into a odel to calculate the value of the anipulated variable required to aintain control at the set point. Figure 2.9c shows how inforation flows forward fro the load to the anipulated variable input of the process. The set point is used to give the syste a coand. (If the controlled variable were used in the calculation instead of the set point, a positive feedback loop would be fored.) A syste, rather than a single control device, is norally used for feedforward loops because it is not always convenient to provide the coputing functions required by the forward loop with a single device or function. Instead, the feedforward syste consists of several devices if ipleented in hardware or several blocks of software if ipleented digitally. The function of these blocks is to ipleent a atheatical odel of the process. Load Balancing A dynaic balance is required to keep the control variable at set point. It is achieved by solving the process aterial- and/or energy-balance equations continuously. When a change in load is sensed, the anipulated variable is autoatically adjusted to the correct value at a rate that keeps the process continually in balance. Although it is theoretically possible to achieve such Set, r Feedfwd syste Load, q Valve or Flowloop K g FIG. 2.9c Feedforward calculates the value of the anipulated variable based on current load and set point. K q g q Process c perfect control, in practice the syste cannot be ade to duplicate the process balance exactly. The aterial- and energy-balance equations are not usually difficult to write for a given process. Variations in nonstationary paraeters such as heat-transfer and ass-transfer coefficients do not ordinarily affect the perforance of a feedforward syste. The load coponents are usually inflow or outflow when level or pressure is to be controlled, feed flow and feed coposition where product coposition is to be controlled, and feed flow and teperature where product teperature is to be controlled. Flow is the priary coponent of load in alost every application because it can change widely and rapidly. Coposition and teperature are less likely to exhibit such wide excursions, and their rate of change is usually liited by upstrea capacity. In ost feedforward systes these secondary load coponents are left out; their effects are accoodated by the feedback controller. The output of a feedforward syste should be an accurately controlled flow rate, if possible. This controlled flow rate cannot usually be obtained by anipulating the control valve directly, since valve characteristics are nonlinear and inconsistent, and the delivered flow is subject to such external influences as upstrea and downstrea pressure variations. Therefore, ost feedforward systes depend on soe easureent and feedback control of flow to obtain an accurate anipulation of the flow rate. Only when the range of the anipulated variable exceeds that which is available in floweters, or the required speed of response exceeds that of a flow loop, should one consider having the valves positioned directly, and in such cases care ust be taken to obtain a reproducible response. Steady-State Model The first step in designing a feedforward control syste is to for a steady-state atheatical odel of the process. The odel equations are solved for the anipulated variable, which is to be the output of the syste. Then the set point is substituted for the controlled variable in the odel. This process will be deonstrated by using the exaple of teperature control in a heat exchanger (Figure 2.9d). 2 A liquid flowing at rate W is to be heated fro teperature T to controlled teperature T 2 by stea flowing at rate W s. The energy balance, excluding any losses, is WC ( T T ) = W λ p 2 2.9(8) Coefficient C p is the heat capacity of the liquid and λ is the latent heat given up by the stea in condensing. Solving for W s yields W WK ( T T ) s = s 2set 2.9(9)

4 76 Control Theory X W FY T FY T 2set T W s FIC Contr. var. Set Offset fro 0% error Accurate static odel Liquid FT TT TT HIC T 2set Heat Exchanger FIG. 2.9d Stea flow is calculated here to satisfy the steady-state heat balance. where K = C p /λ, set point T 2set replaces T 2, and W and T are load coponents to which the control syste ust respond. Note that this energy balance is iperfect, because it disregards inor loads such as the heat losses to the environent and any variations in the enthalpy of the stea. Another factor is that even the two ajor loads cannot be easured perfectly, as there are no error-free sensors of flow and teperature. The net result is a odel that is approxiately an energy balance around this process. An ipleentation of Equation 2.9(9) is given in Figure 2.9d. A control station introduces T 2set into the difference block [ ] where input T is subtracted. The gain of this block is adjusted to provide the constant K. Its output is then ultiplied by the liquid flow signal W to produce the required set point for stea flow. During startup, if the actual value of the controlled variable does not equal the set point, an adjustent is ade to gain K. Then, after aking this adjustent, if the controlled variable does not return to the set point following a change in load, the presence of an error in the syste is indicated. The error ay be in one of the calculations, a sensor, or soe other factor affecting the heat balance that was not included in the syste. Heat losses and variations in stea pressure are two possible sources of error in the odel. Since any error can cause a proportional offset, the designer ust weigh the various sources of error and copensate for the largest or ost changeable coponents where practical. For exaple, if stea pressure variations were a source of error, the stea floweter could be pressure copensated. Because of all these errors, feedforward control as shown in Figure 2.9d is seldo used without feedback tri. Later in this section, a feedback controller will be added to autoatically correct the loop for the errors in the odel. T FT Stea Load Manip. var. Tie FIG. 2.9e The dynaic response of the uncopensated feedforward syste when the controlled variable responds faster to changes in load than to changes in anipulated flow. control systes. Two typical dynaic responses of the controlled variable are shown in Figure 2.9e; the offset in the dashed curve results fro a 0% error in the static calculations. If there is no static error, the teperature returns to set point eventually. However, the dynaic balance is issing. The decreased stea flow does not result in an instantaneous decrease in heat-transfer rate, because the teperature difference across the heat-transfer surface ust first be decreased, and this requires a decrease in shell pressure. The lower shell pressure at lower loads eans that the shell contains less stea as the load decreases. Since the static feedforward syste does not overtly lower the stea inventory in the shell, it is lowered through a lagging heat-transfer rate, resulting in a transient rise in exit teperature on a falling load; conversely, a transient teperature drop follows a rising load. This transient error reveals a dynaic ibalance in the syste on a drop in load, the exchanger teporarily needs less heat than the stea flow controller is allowing. In order to correct the transient teperature error, a further cutback in the flow of energy ust be applied to the exchanger beyond what is required for the steady-state balance. A ore general explanation is provided through the use of Figure 2.9f. This figure shows the load and the anipulated Feedforward dynaics + t s + t q s + t s Dynaic Model With the feedforward odel as ipleented in Figure 2.9d, a step decrease in liquid flow results in a siultaneous and proportional step decrease in stea flow. Transient errors following a change in load are to be expected in feedforward Load q + t q s Process FIG. 2.9f A feedforward dynaic odel for a general process. c

5 2.9 Feedback and Feedforward Control 77 variable as entering the process at different points, where they encounter different dynaic eleents. In the heatexchanger exaple, the liquid enters the tubes while the stea enters the shell. The heat capacities of the two locations are different. As a result, the controlled variable (the liquid teperature) responds ore rapidly to a change in liquid flow than to a change in stea flow. Thus, the lag tie of the load input is less than that of the anipulated variable. The objective of the feedforward syste is to balance the anipulated variable against the load, providing the forward loop with copensating dynaic eleents. Neglecting steadystate considerations, Figure 2.9f shows what the dynaic odel ust contain to balance a process having parallel firstorder lags. The lag tie τ q of the load input ust be duplicated in the forward loop, and the lag tie τ in the anipulated input of the process ust be cancelled. Thus the forward loop should contain a lag divided by a lag. Since the inverse of a lag is a lead, the dynaic copensating function is a leadlag. The lead tie constant should be equal to τ, and the lag tie constant should equal τ q. In the case of the heat exchanger, the fact that τ is longer than τ q causes the teperature rise on a load decrease. This is the direction in which the load change would drive the process in the absence of control. In transfer function for, the response of a lead-lag unit is s Gs ()= + τ + τ s 2.9(0) where τ is the lead tie constant, τ 2 is the lag tie constant, and s is the Laplace operator. The frequency response of a feedforward loop is usually not significant, since forward loops cannot oscillate. The ost severe test for a forward loop is a step change in load. The response of the lead-lag unit output y to a step in its input x is τ τ2 yt () = xt () + ε τ 2.9() Its axiu dynaic gain (at t = 0) is the lead/lag ratio τ /τ 2. The response curve then decays exponentially fro this axiu at the rate of the lag tie constant τ 2, with 63.2% recovery reached at t = τ 2 as shown in Figure 2.9g. The figure shows soe deadtie copensation added to the lead-lag function however, it is not required for the heat exchanger. When properly adjusted, the dynaic copensation afforded by the lead-lag unit can produce the controlled variable response having the for shown in Figure 2.9h. This is a signature curve of a feedforward-controlled process since ost processes do not consist of siple first-order lags, a first-order lead-lag unit cannot produce a perfect dynaic balance. Yet a ajor reduction in peak deviation is achieved over the uncopensated response, and the area can be equally distributed on both sides of the set point. 2 2 t/ τ2 t t2 Load t2 FIG. 2.9g The step response of a lead-lag unit with deadtie copensation. The lead-lag unit is applied to the flow signal as [ f(t)] in Figure 2.9i. Its settings are particular to that load input, and therefore it ust be placed where it can odify that signal and no other. No dynaic copensator is included for T, as it is not expected to change quickly. Tuning the Lead-Lag Unit First, a load change should be introduced without dynaic copensation by setting the lead and lag tie constants at equal values or zero to observe the direction of the error. If the resulting response is in the direction that the load would drive it, the lead tie should exceed the lag tie; if not, the lag tie should be greater. Next, easure the tie required for the controlled variable to reach its axiu or iniu value. This tie should be introduced into the lead-lag unit as the saller tie constant. Thus, if the lead doinates, this would be the lag setting. For exaple, in Figure 2.9e the teperature responds in the direction the load would drive it, and therefore the lead tie setting should exceed the lag tie setting. The tie required for the teperature in Figure 2.9e to rise to its peak should be set equal to the saller tie constant, in this case equal to the lag tie setting. Set the greater tie constant at twice Contr. var. Manip. var. Load Manip. var. Deadtie 63% resp. Set Tie Tie FIG. 2.9h The typical step response of a dynaically copensated feedforward syste.

6 78 Control Theory Liquid f(t) FY FT TT T TT FIG. 2.9i Feedforward feedback control of the heat exchanger, with dynaic copensation. this value and repeat the load change. If the peak deviation is reduced, but the curve is still not equally distributed around the set point, increase the greater tie constant and repeat the load change. The area (integrated error) of the response curve will be equally distributed about the set point if the difference between the lead and lag settings is correct. Once this equalization is obtained, both settings should be increased or decreased with their difference constant until a iniu error aplitude is achieved. When the controller is properly tuned, the stea energy equivalent of the transient (area between teperature and its set point) in Figure 2.9e should atch the extra stea introduced (area above input step change) in Figure 2.9g. If there is a significantly greater deadtie in the load path than in the anipulated-variable path, deadtie copensation ay be added, which will reduce the peak deviation. However, if the deadtie in the anipulated-variable path is the longer of the two, exact copensation is ipossible, and this is the case in the heat exchanger. Still, careful adjustent of the lead-lag copensator can result in zero integrated error. Adding a Feedback Loop X W FY FY T TIC Heat Exchanger FIC Any offset resulting fro steady-state errors can be eliinated by adding feedback. This can be done by replacing the feedforward set point with a controller, as shown in Figure 2.9i. The feedback controller adjusts the set point of the feedforward syste in cascade, while the feedforward syste adjusts the set point of the anipulated flow controller in cascade. The feedback controller should have the sae control odes as it would without feedforward control, but the settings need not be as tight. The feedback controller reacts to a disturbance by creating another disturbance in the opposite direction one-half cycle later. However, the feedforward syste positions the anipulated variable so that the error in the controlled variable disappears. If acted upon by a tightly W s T FT Stea set feedback controller, the correct position calculated by feedforward will be altered, producing another disturbance that prolongs the settling tie of the syste. The feedback controller ust have the integral ode to eliinate any steady-state offset that ight be caused by errors in the sensors, the odel, or the calculations. However, the integrated error it sustains following a load change is arkedly reduced copared to what it would be following a load change without feedforward. Now its output changes only an aount equal to the change in the error in the feedforward calculation. For sall load changes, this change in error would approach zero, leading to an integrated error of zero, with or without dynaic copensation. In the presence of feedback, but without dynaic copensation, the transient in Figure 2.9e would be balanced by a following transient on the other side of set point, until the area is equalized. This usually prolongs the response without reducing the peak deviation. The dynaic copensator therefore needs to be tuned with the feedback controller in anual, to approach zero integrated error before adding feedback. Linear and Bilinear Feedforward Systes The relationship shown in the energy balance of Equation 2.9(8) is bilinear: stea flow is related to both liquid flow and to teperature rise in a linear anner, but as their product rather than their su or difference. This distinction is crucial to the successful functioning of feedforward control because a bilinear syste has variable gains. The feedforward gain the ratio of stea flow to liquid flow varies directly with the required teperature rise. This variable gain is accoodated in the syste by the ultiplier [ ]. If inlet teperature were to vary so slowly that it did not require feedforward copensation or if it were not easured, then feedforward fro flow alone would be required. Yet the required feedforward gain the ratio W s /W would vary directly with T 2set T. If either T 2set or T were to change, the feedforward gain should change with it. For a bilinear process, then, a linear feedforward syste having a constant gain is a isfit it is accurate at only one cobination of secondary load and set point. To be sure, feedback control can tri any errors in the feedforward calculation and return the deviation to zero. But if the error is in the gain ter, the feedback controller should adjust the feedforward gain, through a ultiplier otherwise the gain will be incorrect for all subsequent flow changes, resulting in only partial feedforward correction or even overcorrection. Eliinating inlet teperature as a variable in Figure 2.9i also eliinates the difference block [ ]. The controller output then goes directly to the ultiplier, and in reducing the teperature deviation to zero, will assue a value representing K(T 2set T ). If either T 2set or T were then to vary, the feedback controller would respond to the resulting deviation by changing its output until the deviation returns to zero. In so doing, it has changed the gain of the ultiplier to the new value of K(T 2set T ).

7 2.9 Feedback and Feedforward Control 79 Control of coposition is also a bilinear process. The required ratio of two flows entering a blender, for exaple, is a function of their individual copositions and that of the blend. Siilarly, the ratio of product flow to feed rate or stea flow to feed rate in a distillation colun or evaporator also varies with feed and product copositions. In general, teperature and coposition loops are bilinear, and their flow ratios should always be adjusted through a ultiplier. By contrast, pressure and level loops are linear, and their feedforward gain can be constant. An exaple of the latter is three-eleent dru-level control, where each unit of stea reoved fro the dru ust be replaced by an equal unit of feedwater the feedforward gain is constant at.0. This distinction is crucial because any advanced ultivariable control systes are based on a linear atrix with constant coefficients. The plant is tested to develop the coefficients that relate the variables, and the control atrix is built fro these results. When the plant operating conditions ove away fro the original test conditions, the fixed coefficients in the control atrix ay no longer represent the process relationships accurately, degrading the perforance of the feedforward loops. This ay require frequent retesting and recalibration of the atrix. Soe ultivariable systes contain ultiple atrices, each corresponding to its own set of operating conditions and switched into service when those conditions develop. Self-tuning feedforward control 3 is also available, applied either in a linear (additive) or bilinear (ultiplicative) anner, as configured anually. The paraeters that are tuned adaptively are the steady-state gain and a lag copensator. The steady-state gain is adjusted to iniize integrated error following a change in easured load, and the lag is adjusted to iniize integral-square error. Load changes ust be sharp, clean steps, with settling tie allowed between the, for tuning to proceed effectively. But once tuned, rando disturbances can be accoodated. However, self-tuning is only recoended where precise odeling and flow control are unavailable. Perforance The use of feedback in a feedforward syste does not detract fro the perforance iproveent that was gained by feedforward control. Without feedforward, the feedback controller was required to change its output to follow all changes in load. With feedforward, the feedback controller ust only change its output by an aount equal to what the feedforward syste fails to correct. If a feedforward syste applied to a heat exchanger could control the stea flow to within 2% of that required by the load, the feedback controller would only be required to adjust its output to copensate for 2% of a load change, rather than the full aount. This reduction of in Equation 2.9(7) by 50/ results in the reduction of E by the sae ratio. Reductions by 0/ in errors resulting fro load changes are relatively coon, and iproveents of 00/ have been achieved in soe systes. Teperature PID control Dynaic feedforward Tie, s FIG. 2.9j Coparison of feedforward control of a heat exchanger against optially tuned PID feedback control. Figure 2.9j illustrates the control perforance of a steaheated exchanger that has experienced a step decrease in process flow. The static feedforward odel reduces the peak deviation substantially but extends settling tie and does not iprove on the integrated error. Dynaic copensation is seen to be essential in axiizing perforance in this exaple. It is not perfect because the deadtie in the path of the anipulated variable is soewhat longer than that in the load path, and that lost tie cannot be ade up. Still, the lead-lag copensator has been tuned to eliinate any integrated error. The feedforward syste is ore costly and requires ore engineering effort than a feedback syste does, so prior to design and installation, the control iproveent it brings ust be deterined to be worthwhile. Most feedforward systes have been applied to processes that are very sensitive to disturbances and slow to respond to corrective action and to product streas that are relatively high in value. Distillation coluns of 50 trays or ore have been the principal systes controlled with this technology. Boilers, ultipleeffect evaporators, direct-fired heaters, waste neutralization plants, solids dryers, copressors, and other hard-to-control processes have also benefited fro feedforward control. Variable Paraeters Static feedforward Heat exchangers are characterized by gain and dynaics that vary with flow. The settings of both the PID controller and the feedforward copensator that produced the results shown in Figure 2.9j were optiized for the final flow and would not be optiu for any other flow. Differentiation of the heat-balance equation, Equation 2.9(8), shows the process gain to vary inversely with process flow W: dt 2 = λ dw WC s 2.9(2) If PID control alone is applied, this gain variation can be copensated through the use of an equal-percentage stea p 240

8 80 Control Theory valve, whose gain varies directly with stea flow. However, when the PID controller is cobined with bilinear feedforward, its output passes through the ultiplier shown in Figure 2.9i, where it is ultiplied by process flow. This operation keeps the PID loop gain constant, while anipulating stea flow linearly. Heat-exchanger dynaics deadtie and lags also vary inversely with process flow, as is typical of once-through processes, i.e., where no recirculation takes place. This proble is less easily solved. The PID controller ust have its integral tie set relative to the slowest response (lowest expected flow) and its derivative tie set relative to its fastest response (highest expected flow), or ideally, prograed as a function of easured flow. Otherwise, instability ay result at extree flow rates. Ideally, lead and lag settings of the feedforward dynaic copensator should also be prograed as a function of easured flow. However, the penalty for not doing so is not severe. The lead-lag ratio is not subject to change because the lags on both sides of the process change in the sae proportion with flow. The priary purpose of the dynaic copensator is to iniize the peak deviation following the load change, and this is accoplished by the dynaic gain of the copensator, which is its lead-lag ratio. Any subsequent variation in integrated error is eliinated by the PID controller. Although paraeter variations can also be accoodated by a self-tuning feedforward copensator and feedback controller, these ethods are less accurate in arriving at optiu settings and are always late, having tuned for the last flow condition and not for the next one. References. Shinskey, F. G., Process Control Systes, 4th ed., New York: McGraw- Hill, 996, pp Shinskey, F. G., Feedforward Control Applied, ISA Journal, Noveber Hansen, P. D., Badavas, P. C., and Kinney, T. B., New Adaptive Multivariable Controller, presented at the 3rd Intl. Control Eng. Conference, Chicago, 994. Bibliography Butterworths, W. R., Advanced Process Control, New York: McGraw-Hill, 98. Gore, F., Add Copensation for Feedforward Control, Instruents and Control Systes, March 979. Paler, R., Nonlinear Feedforward Can Reduce Servo Settling Tie, Control Engineering, March 978. Seborg, D., et al., Process Dynaics and Control, New York: Wiley, 989. Shinskey, F.G., Reflections on CPC-III, Cheical Process Control-CPC III, T. J. McAvoy, et al. (Eds.), Asterda: Elserier Sci. Pub., 986. Shinskey, F. G., Soothing Out Copressor Control, Cheical Engineering, February 999.

CHAPTER 19: Single-Loop IMC Control

CHAPTER 19: Single-Loop IMC Control When I coplete this chapter, I want to be able to do the following. Recognize that other feedback algoriths are possible Understand the IMC structure and how it provides the essential control features

More information

Model Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon

Model Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential

More information

Experimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis

Experimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna

More information

Reducing Vibration and Providing Robustness with Multi-Input Shapers

Reducing Vibration and Providing Robustness with Multi-Input Shapers 29 Aerican Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -2, 29 WeA6.4 Reducing Vibration and Providing Robustness with Multi-Input Shapers Joshua Vaughan and Willia Singhose Abstract

More information

COS 424: Interacting with Data. Written Exercises

COS 424: Interacting with Data. Written Exercises COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well

More information

A method to determine relative stroke detection efficiencies from multiplicity distributions

A method to determine relative stroke detection efficiencies from multiplicity distributions A ethod to deterine relative stroke detection eiciencies ro ultiplicity distributions Schulz W. and Cuins K. 2. Austrian Lightning Detection and Inoration Syste (ALDIS), Kahlenberger Str.2A, 90 Vienna,

More information

Smith Predictor Based-Sliding Mode Controller for Integrating Process with Elevated Deadtime

Smith Predictor Based-Sliding Mode Controller for Integrating Process with Elevated Deadtime Sith Predictor Based-Sliding Mode Controller for Integrating Process with Elevated Deadtie Oscar Caacho, a, * Francisco De la Cruz b a Postgrado en Autoatización e Instruentación. Grupo en Nuevas Estrategias

More information

Analyzing Simulation Results

Analyzing Simulation Results Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient

More information

Q5 We know that a mass at the end of a spring when displaced will perform simple m harmonic oscillations with a period given by T = 2!

Q5 We know that a mass at the end of a spring when displaced will perform simple m harmonic oscillations with a period given by T = 2! Chapter 4.1 Q1 n oscillation is any otion in which the displaceent of a particle fro a fixed point keeps changing direction and there is a periodicity in the otion i.e. the otion repeats in soe way. In

More information

Ch 12: Variations on Backpropagation

Ch 12: Variations on Backpropagation Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith

More information

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving

More information

MSEC MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL SOLUTION FOR MAINTENANCE AND PERFORMANCE

MSEC MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL SOLUTION FOR MAINTENANCE AND PERFORMANCE Proceeding of the ASME 9 International Manufacturing Science and Engineering Conference MSEC9 October 4-7, 9, West Lafayette, Indiana, USA MSEC9-8466 MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL

More information

An Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control

An Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control An Extension to the Tactical Planning Model for a Job Shop: Continuous-Tie Control Chee Chong. Teo, Rohit Bhatnagar, and Stephen C. Graves Singapore-MIT Alliance, Nanyang Technological Univ., and Massachusetts

More information

Data-Driven Imaging in Anisotropic Media

Data-Driven Imaging in Anisotropic Media 18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands

More information

High-Speed Smooth Cyclical Motion of a Hydraulic Cylinder

High-Speed Smooth Cyclical Motion of a Hydraulic Cylinder 846 1 th International Conference on Fleible Autoation and Intelligent Manufacturing 00, Dresden, Gerany High-Speed Sooth Cyclical Motion of a Hydraulic Cylinder Nebojsa I. Jaksic Departent of Industrial

More information

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co

More information

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators Suppleentary Inforation for Design of Bending Multi-Layer Electroactive Polyer Actuators Bavani Balakrisnan, Alek Nacev, and Elisabeth Sela University of Maryland, College Park, Maryland 074 1 Analytical

More information

Support Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization

Support Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering

More information

Torsion Experiment. Encoder #3 ( 3 ) Third encoder/disk for Model 205a only. Figure 1: ECP Torsion Experiment

Torsion Experiment. Encoder #3 ( 3 ) Third encoder/disk for Model 205a only. Figure 1: ECP Torsion Experiment Torsion Experient Introduction For the Torsion lab, there are two required experients to perfor and one extra credit assignent at the end. In experient 1, the syste paraeters need to be identified so that

More information

Topic 5a Introduction to Curve Fitting & Linear Regression

Topic 5a Introduction to Curve Fitting & Linear Regression /7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline

More information

An Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period

An Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction

More information

Figure 1: Equivalent electric (RC) circuit of a neurons membrane

Figure 1: Equivalent electric (RC) circuit of a neurons membrane Exercise: Leaky integrate and fire odel of neural spike generation This exercise investigates a siplified odel of how neurons spike in response to current inputs, one of the ost fundaental properties of

More information

Lecture #8-3 Oscillations, Simple Harmonic Motion

Lecture #8-3 Oscillations, Simple Harmonic Motion Lecture #8-3 Oscillations Siple Haronic Motion So far we have considered two basic types of otion: translation and rotation. But these are not the only two types of otion we can observe in every day life.

More information

Block designs and statistics

Block designs and statistics Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent

More information

REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION

REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION ISSN 139 14X INFORMATION TECHNOLOGY AND CONTROL, 008, Vol.37, No.3 REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION Riantas Barauskas, Vidantas Riavičius Departent of Syste Analysis, Kaunas

More information

Pseudo Dynamic Model Reference Nonlinear Control of a Continuous Bioreactor with Input Multiplicities

Pseudo Dynamic Model Reference Nonlinear Control of a Continuous Bioreactor with Input Multiplicities seudo Dynaic Model Reference Nonlinear Control of a Continuous Bioreactor with Input Multiplicities G.rabhaker Reddy*,V. Raesh Kuar and B.Spandana Abstract In the present work, the odel reference nonlinear

More information

Ocean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers

Ocean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers Ocean 40 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers 1. Hydrostatic Balance a) Set all of the levels on one of the coluns to the lowest possible density.

More information

Polygonal Designs: Existence and Construction

Polygonal Designs: Existence and Construction Polygonal Designs: Existence and Construction John Hegean Departent of Matheatics, Stanford University, Stanford, CA 9405 Jeff Langford Departent of Matheatics, Drake University, Des Moines, IA 5011 G

More information

Physically Based Modeling CS Notes Spring 1997 Particle Collision and Contact

Physically Based Modeling CS Notes Spring 1997 Particle Collision and Contact Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which

More information

General Properties of Radiation Detectors Supplements

General Properties of Radiation Detectors Supplements Phys. 649: Nuclear Techniques Physics Departent Yarouk University Chapter 4: General Properties of Radiation Detectors Suppleents Dr. Nidal M. Ershaidat Overview Phys. 649: Nuclear Techniques Physics Departent

More information

ma x = -bv x + F rod.

ma x = -bv x + F rod. Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous

More information

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay

More information

A Simple Regression Problem

A Simple Regression Problem A Siple Regression Proble R. M. Castro March 23, 2 In this brief note a siple regression proble will be introduced, illustrating clearly the bias-variance tradeoff. Let Y i f(x i ) + W i, i,..., n, where

More information

Interactive Markov Models of Evolutionary Algorithms

Interactive Markov Models of Evolutionary Algorithms Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary

More information

Ph 20.3 Numerical Solution of Ordinary Differential Equations

Ph 20.3 Numerical Solution of Ordinary Differential Equations Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing

More information

ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER

ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER IEPC 003-0034 ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER A. Bober, M. Guelan Asher Space Research Institute, Technion-Israel Institute of Technology, 3000 Haifa, Israel

More information

All Excuses must be taken to 233 Loomis before 4:15, Monday, April 30.

All Excuses must be taken to 233 Loomis before 4:15, Monday, April 30. Miscellaneous Notes he end is near don t get behind. All Excuses ust be taken to 233 Loois before 4:15, Monday, April 30. he PHYS 213 final exa ties are * 8-10 AM, Monday, May 7 * 8-10 AM, uesday, May

More information

Time-of-flight Identification of Ions in CESR and ERL

Time-of-flight Identification of Ions in CESR and ERL Tie-of-flight Identification of Ions in CESR and ERL Eric Edwards Departent of Physics, University of Alabaa, Tuscaloosa, AL, 35486 (Dated: August 8, 2008) The accuulation of ion densities in the bea pipe

More information

CONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. IX Uncertainty Models For Robustness Analysis - A. Garulli, A. Tesi and A. Vicino

CONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. IX Uncertainty Models For Robustness Analysis - A. Garulli, A. Tesi and A. Vicino UNCERTAINTY MODELS FOR ROBUSTNESS ANALYSIS A. Garulli Dipartiento di Ingegneria dell Inforazione, Università di Siena, Italy A. Tesi Dipartiento di Sistei e Inforatica, Università di Firenze, Italy A.

More information

Force and dynamics with a spring, analytic approach

Force and dynamics with a spring, analytic approach Force and dynaics with a spring, analytic approach It ay strie you as strange that the first force we will discuss will be that of a spring. It is not one of the four Universal forces and we don t use

More information

Feature Extraction Techniques

Feature Extraction Techniques Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that

More information

A Better Algorithm For an Ancient Scheduling Problem. David R. Karger Steven J. Phillips Eric Torng. Department of Computer Science

A Better Algorithm For an Ancient Scheduling Problem. David R. Karger Steven J. Phillips Eric Torng. Department of Computer Science A Better Algorith For an Ancient Scheduling Proble David R. Karger Steven J. Phillips Eric Torng Departent of Coputer Science Stanford University Stanford, CA 9435-4 Abstract One of the oldest and siplest

More information

Reading from Young & Freedman: For this topic, read the introduction to chapter 25 and sections 25.1 to 25.3 & 25.6.

Reading from Young & Freedman: For this topic, read the introduction to chapter 25 and sections 25.1 to 25.3 & 25.6. PHY10 Electricity Topic 6 (Lectures 9 & 10) Electric Current and Resistance n this topic, we will cover: 1) Current in a conductor ) Resistivity 3) Resistance 4) Oh s Law 5) The Drude Model of conduction

More information

Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi.

Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi. Seisic Analysis of Structures by K Dutta, Civil Departent, II Delhi, New Delhi. Module 5: Response Spectru Method of Analysis Exercise Probles : 5.8. or the stick odel of a building shear frae shown in

More information

26 Impulse and Momentum

26 Impulse and Momentum 6 Ipulse and Moentu First, a Few More Words on Work and Energy, for Coparison Purposes Iagine a gigantic air hockey table with a whole bunch of pucks of various asses, none of which experiences any friction

More information

16.512, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 30: Dynamics of Turbopump Systems: The Shuttle Engine

16.512, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 30: Dynamics of Turbopump Systems: The Shuttle Engine 6.5, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 30: Dynaics of Turbopup Systes: The Shuttle Engine Dynaics of the Space Shuttle Main Engine Oxidizer Pressurization Subsystes Selected Sub-Model

More information

International Journal of Scientific & Engineering Research, Volume 4, Issue 9, September ISSN

International Journal of Scientific & Engineering Research, Volume 4, Issue 9, September ISSN International Journal of Scientific & Engineering Research, Volue 4, Issue 9, Septeber-3 44 ISSN 9-558 he unscented Kalan Filter for the Estiation the States of he Boiler-urbin Model Halieh Noorohaadi,

More information

Some Perspective. Forces and Newton s Laws

Some Perspective. Forces and Newton s Laws Soe Perspective The language of Kineatics provides us with an efficient ethod for describing the otion of aterial objects, and we ll continue to ake refineents to it as we introduce additional types of

More information

EMPIRICAL COMPLEXITY ANALYSIS OF A MILP-APPROACH FOR OPTIMIZATION OF HYBRID SYSTEMS

EMPIRICAL COMPLEXITY ANALYSIS OF A MILP-APPROACH FOR OPTIMIZATION OF HYBRID SYSTEMS EMPIRICAL COMPLEXITY ANALYSIS OF A MILP-APPROACH FOR OPTIMIZATION OF HYBRID SYSTEMS Jochen Till, Sebastian Engell, Sebastian Panek, and Olaf Stursberg Process Control Lab (CT-AST), University of Dortund,

More information

SUPERIOR-ORDER CURVATURE-CORRECTED PROGRAMMABLE VOLTAGE REFERENCES

SUPERIOR-ORDER CURVATURE-CORRECTED PROGRAMMABLE VOLTAGE REFERENCES SUPEIO-ODE CUATUE-COECTED POGAMMABLE OLTAGE EFEENCES Cosin Popa e-ail: cosin@golanapubro Faculty of Electronics and Telecounications, University Politehnica of Bucharest, B dul Iuliu Maniu 1-3, Bucuresti,

More information

Ensemble Based on Data Envelopment Analysis

Ensemble Based on Data Envelopment Analysis Enseble Based on Data Envelopent Analysis So Young Sohn & Hong Choi Departent of Coputer Science & Industrial Systes Engineering, Yonsei University, Seoul, Korea Tel) 82-2-223-404, Fax) 82-2- 364-7807

More information

Analysis and Implementation of a Hardware-in-the-Loop Simulation

Analysis and Implementation of a Hardware-in-the-Loop Simulation DINAME 27 - Proceedings of the XVII International Syposiu on Dynaic Probles of Mechanics A. T. Fleury, D. A. Rade, P. R. G. Kurka (Editors), ABCM, São Sebastião, SP, Brail, March 5-, 27 Analysis and Ipleentation

More information

Optical Properties of Plasmas of High-Z Elements

Optical Properties of Plasmas of High-Z Elements Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer

More information

ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics

ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents

More information

Distillation. The Continuous Column. Learning Outcomes. Recap - VLE for Meth H 2 O. Gavin Duffy School of Electrical Engineering DIT Kevin Street

Distillation. The Continuous Column. Learning Outcomes. Recap - VLE for Meth H 2 O. Gavin Duffy School of Electrical Engineering DIT Kevin Street Distillation The Continuous Colun Gavin Duffy School of Electrical Engineering DIT Kevin Street Learning Outcoes After this lecture you should be able to.. Describe how continuous distillation works List

More information

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either

More information

Analysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning

Analysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning Analysis of Ipulsive Natural Phenoena through Finite Difference Methods A MATLAB Coputational Project-Based Learning Nicholas Kuia, Christopher Chariah, Mechatronics Engineering, Vaughn College of Aeronautics

More information

Chapter 18. Process Control

Chapter 18. Process Control Chapter 18. Process Control 18. INTRODUCTION The atheatical odels described in previous chapters serve adequately for successful configuration and ipleentation of unit processes as well as distributed

More information

DESIGN, SIMULATION AND IMPLEMENTATION OF DECENTRALIZED PI CONTROLLERS FOR A MULTIVARIABLE FEED BLENDING SYSTEM AT THE FALCONBRIDGE SUDBURY SMELTER

DESIGN, SIMULATION AND IMPLEMENTATION OF DECENTRALIZED PI CONTROLLERS FOR A MULTIVARIABLE FEED BLENDING SYSTEM AT THE FALCONBRIDGE SUDBURY SMELTER DESIGN, SIMULATION AND IMPLEMENTATION OF DECENTRALIZED PI CONTROLLERS FOR A MULTIVARIABLE FEED BLENDING SYSTEM AT THE FALCONBRIDGE SUDBURY SMELTER B. Boulet McGill University Departent of Electrical and

More information

Non-Parametric Non-Line-of-Sight Identification 1

Non-Parametric Non-Line-of-Sight Identification 1 Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,

More information

Using a De-Convolution Window for Operating Modal Analysis

Using a De-Convolution Window for Operating Modal Analysis Using a De-Convolution Window for Operating Modal Analysis Brian Schwarz Vibrant Technology, Inc. Scotts Valley, CA Mark Richardson Vibrant Technology, Inc. Scotts Valley, CA Abstract Operating Modal Analysis

More information

III.H Zeroth Order Hydrodynamics

III.H Zeroth Order Hydrodynamics III.H Zeroth Order Hydrodynaics As a first approxiation, we shall assue that in local equilibriu, the density f 1 at each point in space can be represented as in eq.iii.56, i.e. f 0 1 p, q, t = n q, t

More information

Architectural issues in the control of tall multiple-input multiple-output plants

Architectural issues in the control of tall multiple-input multiple-output plants Preprints of the 19th World Congress The International Federation of Autoatic Control Architectural issues in the control of tall ultiple-input ultiple-output plants Pedro Albertos Mario E. Salgado Departent

More information

Self-Erecting Inverted Pendulum: Swing up and Stabilization Control

Self-Erecting Inverted Pendulum: Swing up and Stabilization Control Self-Erecting Inverted Pendulu: Swing up and Stabilization Control S. McGilvray (Winner of the IEEE Life Meber ward for best paper fro the Central Canada Council for the IEEE Student Paper Contest, ) Contents

More information

Inspection; structural health monitoring; reliability; Bayesian analysis; updating; decision analysis; value of information

Inspection; structural health monitoring; reliability; Bayesian analysis; updating; decision analysis; value of information Cite as: Straub D. (2014). Value of inforation analysis with structural reliability ethods. Structural Safety, 49: 75-86. Value of Inforation Analysis with Structural Reliability Methods Daniel Straub

More information

Multi-Scale/Multi-Resolution: Wavelet Transform

Multi-Scale/Multi-Resolution: Wavelet Transform Multi-Scale/Multi-Resolution: Wavelet Transfor Proble with Fourier Fourier analysis -- breaks down a signal into constituent sinusoids of different frequencies. A serious drawback in transforing to the

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering 2.010: Systems Modeling and Dynamics III. Final Examination Review Problems

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering 2.010: Systems Modeling and Dynamics III. Final Examination Review Problems ASSACHUSETTS INSTITUTE OF TECHNOLOGY Departent of echanical Engineering 2.010: Systes odeling and Dynaics III Final Eaination Review Probles Fall 2000 Good Luck And have a great winter break! page 1 Proble

More information

PHY 171. Lecture 14. (February 16, 2012)

PHY 171. Lecture 14. (February 16, 2012) PHY 171 Lecture 14 (February 16, 212) In the last lecture, we looked at a quantitative connection between acroscopic and icroscopic quantities by deriving an expression for pressure based on the assuptions

More information

Hyperbolic Horn Helical Mass Spectrometer (3HMS) James G. Hagerman Hagerman Technology LLC & Pacific Environmental Technologies April 2005

Hyperbolic Horn Helical Mass Spectrometer (3HMS) James G. Hagerman Hagerman Technology LLC & Pacific Environmental Technologies April 2005 Hyperbolic Horn Helical Mass Spectroeter (3HMS) Jaes G Hageran Hageran Technology LLC & Pacific Environental Technologies April 5 ABSTRACT This paper describes a new type of ass filter based on the REFIMS

More information

OPTIMIZATION OF SPECIFIC FACTORS TO PRODUCE SPECIAL ALLOYS

OPTIMIZATION OF SPECIFIC FACTORS TO PRODUCE SPECIAL ALLOYS 5 th International Conference Coputational Mechanics and Virtual Engineering COMEC 2013 24-25 October 2013, Braşov, Roania OPTIMIZATION OF SPECIFIC FACTORS TO PRODUCE SPECIAL ALLOYS I. Milosan 1 1 Transilvania

More information

Fault Tree Modeling for Redundant Multi-Functional Digital Systems

Fault Tree Modeling for Redundant Multi-Functional Digital Systems International Journal of Perforability Engineering, Vol. 3, No. 3, July, 2007, pp. 329-336 RAMS Consultants Printed in India Fault Tree Modeling for Redundant Multi-Functional Digital Systes HYUN GOOK

More information

DETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION

DETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION DETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION Masaki WAKUI 1 and Jun IYAMA and Tsuyoshi KOYAMA 3 ABSTRACT This paper shows a criteria to detect

More information

Birthday Paradox Calculations and Approximation

Birthday Paradox Calculations and Approximation Birthday Paradox Calculations and Approxiation Joshua E. Hill InfoGard Laboratories -March- v. Birthday Proble In the birthday proble, we have a group of n randoly selected people. If we assue that birthdays

More information

Optimal nonlinear Bayesian experimental design: an application to amplitude versus offset experiments

Optimal nonlinear Bayesian experimental design: an application to amplitude versus offset experiments Geophys. J. Int. (23) 155, 411 421 Optial nonlinear Bayesian experiental design: an application to aplitude versus offset experients Jojanneke van den Berg, 1, Andrew Curtis 2,3 and Jeannot Trapert 1 1

More information

Adapting the Pheromone Evaporation Rate in Dynamic Routing Problems

Adapting the Pheromone Evaporation Rate in Dynamic Routing Problems Adapting the Pheroone Evaporation Rate in Dynaic Routing Probles Michalis Mavrovouniotis and Shengxiang Yang School of Coputer Science and Inforatics, De Montfort University The Gateway, Leicester LE1

More information

Neural Dynamic Optimization for Control Systems Part III: Applications

Neural Dynamic Optimization for Control Systems Part III: Applications 502 IEEE TRANSACTIONS ON SYSTEMS, MAN, AND CYBERNETICS PART B: CYBERNETICS, VOL. 31, NO. 4, AUGUST 2001 Neural Dynaic Optiization for Control Systes Part III: Applications Chang-Yun Seong, Meber, IEEE,

More information

When Short Runs Beat Long Runs

When Short Runs Beat Long Runs When Short Runs Beat Long Runs Sean Luke George Mason University http://www.cs.gu.edu/ sean/ Abstract What will yield the best results: doing one run n generations long or doing runs n/ generations long

More information

Randomized Accuracy-Aware Program Transformations For Efficient Approximate Computations

Randomized Accuracy-Aware Program Transformations For Efficient Approximate Computations Randoized Accuracy-Aware Progra Transforations For Efficient Approxiate Coputations Zeyuan Allen Zhu Sasa Misailovic Jonathan A. Kelner Martin Rinard MIT CSAIL zeyuan@csail.it.edu isailo@it.edu kelner@it.edu

More information

Numerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term

Numerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term Nuerical Studies of a Nonlinear Heat Equation with Square Root Reaction Ter Ron Bucire, 1 Karl McMurtry, 1 Ronald E. Micens 2 1 Matheatics Departent, Occidental College, Los Angeles, California 90041 2

More information

AVOIDING PITFALLS IN MEASUREMENT UNCERTAINTY ANALYSIS

AVOIDING PITFALLS IN MEASUREMENT UNCERTAINTY ANALYSIS VOIDING ITFLLS IN ESREENT NERTINTY NLYSIS Benny R. Sith Inchwor Solutions Santa Rosa, Suary: itfalls, both subtle and obvious, await the new or casual practitioner of easureent uncertainty analysis. This

More information

COULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY?

COULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY? Kragujevac J. Sci. 3 (8) 3-44. UDC 53.35 3 COULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY? Nebojša Danilović, Milan Kovačević and Vukota Babović Institute of Physics, Faculty of Science,

More information

Pattern Recognition and Machine Learning. Artificial Neural networks

Pattern Recognition and Machine Learning. Artificial Neural networks Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lessons 7 20 Dec 2017 Outline Artificial Neural networks Notation...2 Introduction...3 Key Equations... 3 Artificial

More information

OSCILLATIONS AND WAVES

OSCILLATIONS AND WAVES OSCILLATIONS AND WAVES OSCILLATION IS AN EXAMPLE OF PERIODIC MOTION No stories this tie, we are going to get straight to the topic. We say that an event is Periodic in nature when it repeats itself in

More information

Use of PSO in Parameter Estimation of Robot Dynamics; Part One: No Need for Parameterization

Use of PSO in Parameter Estimation of Robot Dynamics; Part One: No Need for Parameterization Use of PSO in Paraeter Estiation of Robot Dynaics; Part One: No Need for Paraeterization Hossein Jahandideh, Mehrzad Navar Abstract Offline procedures for estiating paraeters of robot dynaics are practically

More information

INTELLECTUAL DATA ANALYSIS IN AIRCRAFT DESIGN

INTELLECTUAL DATA ANALYSIS IN AIRCRAFT DESIGN INTELLECTUAL DATA ANALYSIS IN AIRCRAFT DESIGN V.A. Koarov 1, S.A. Piyavskiy 2 1 Saara National Research University, Saara, Russia 2 Saara State Architectural University, Saara, Russia Abstract. This article

More information

Risk & Safety in Engineering. Dr. Jochen Köhler

Risk & Safety in Engineering. Dr. Jochen Köhler Risk & afety in Engineering Dr. Jochen Köhler Contents of Today's Lecture Introduction to Classical Reliability Theory tructural Reliability The fundaental case afety argin Introduction to Classical Reliability

More information

Analysis of ground vibration transmission in high precision equipment by Frequency Based Substructuring

Analysis of ground vibration transmission in high precision equipment by Frequency Based Substructuring Analysis of ground vibration transission in high precision equipent by Frequency Based Substructuring G. van Schothorst 1, M.A. Boogaard 2, G.W. van der Poel 1, D.J. Rixen 2 1 Philips Innovation Services,

More information

Chapter 2 General Properties of Radiation Detectors

Chapter 2 General Properties of Radiation Detectors Med Phys 4RA3, 4RB3/6R3 Radioisotopes and Radiation Methodology -1 Chapter General Properties of Radiation Detectors Ionizing radiation is ost coonly detected by the charge created when radiation interacts

More information

Example A1: Preparation of a Calibration Standard

Example A1: Preparation of a Calibration Standard Suary Goal A calibration standard is prepared fro a high purity etal (cadiu) with a concentration of ca.1000 g l -1. Measureent procedure The surface of the high purity etal is cleaned to reove any etal-oxide

More information

Fast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials

Fast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials Fast Montgoery-like Square Root Coputation over GF( ) for All Trinoials Yin Li a, Yu Zhang a, a Departent of Coputer Science and Technology, Xinyang Noral University, Henan, P.R.China Abstract This letter

More information

Energy and Momentum: The Ballistic Pendulum

Energy and Momentum: The Ballistic Pendulum Physics Departent Handout -10 Energy and Moentu: The Ballistic Pendulu The ballistic pendulu, first described in the id-eighteenth century, applies principles of echanics to the proble of easuring the

More information

Kinematics and dynamics, a computational approach

Kinematics and dynamics, a computational approach Kineatics and dynaics, a coputational approach We begin the discussion of nuerical approaches to echanics with the definition for the velocity r r ( t t) r ( t) v( t) li li or r( t t) r( t) v( t) t for

More information

Principles of Optimal Control Spring 2008

Principles of Optimal Control Spring 2008 MIT OpenCourseWare http://ocw.it.edu 16.323 Principles of Optial Control Spring 2008 For inforation about citing these aterials or our Ters of Use, visit: http://ocw.it.edu/ters. 16.323 Lecture 10 Singular

More information

Sensorless Control of Induction Motor Drive Using SVPWM - MRAS Speed Observer

Sensorless Control of Induction Motor Drive Using SVPWM - MRAS Speed Observer Journal of Eerging Trends in Engineering and Applied Sciences (JETEAS) 2 (3): 509-513 Journal Scholarlink of Eerging Research Trends Institute in Engineering Journals, 2011 and Applied (ISSN: 2141-7016)

More information

An earlier article in this column considered the problem

An earlier article in this column considered the problem --- CALC CORNER Estiating nternal Air Cooling Teperature Reduction in a Closed Box Utilizing Theroelectrically Enhanced Heat Rejection Previously published in February, 2013 Bob Sions BM Retired The following

More information

Probability Distributions

Probability Distributions Probability Distributions In Chapter, we ephasized the central role played by probability theory in the solution of pattern recognition probles. We turn now to an exploration of soe particular exaples

More information

NBN Algorithm Introduction Computational Fundamentals. Bogdan M. Wilamoswki Auburn University. Hao Yu Auburn University

NBN Algorithm Introduction Computational Fundamentals. Bogdan M. Wilamoswki Auburn University. Hao Yu Auburn University NBN Algorith Bogdan M. Wilaoswki Auburn University Hao Yu Auburn University Nicholas Cotton Auburn University. Introduction. -. Coputational Fundaentals - Definition of Basic Concepts in Neural Network

More information

Uniaxial compressive stress strain model for clay brick masonry

Uniaxial compressive stress strain model for clay brick masonry Uniaxial copressive stress strain odel for clay brick asonry Heant B. Kaushik, Durgesh C. Rai* and Sudhir K. Jain Departent of Civil Engineering, Indian Institute of Technology Kanpur, Kanpur 208 016,

More information

Lower Bounds for Quantized Matrix Completion

Lower Bounds for Quantized Matrix Completion Lower Bounds for Quantized Matrix Copletion Mary Wootters and Yaniv Plan Departent of Matheatics University of Michigan Ann Arbor, MI Eail: wootters, yplan}@uich.edu Mark A. Davenport School of Elec. &

More information

Ufuk Demirci* and Feza Kerestecioglu**

Ufuk Demirci* and Feza Kerestecioglu** 1 INDIRECT ADAPTIVE CONTROL OF MISSILES Ufuk Deirci* and Feza Kerestecioglu** *Turkish Navy Guided Missile Test Station, Beykoz, Istanbul, TURKEY **Departent of Electrical and Electronics Engineering,

More information