Measurement and Prediction of Discharge Coefficients in Highly Compressible Pulsating Flows to Improve EGR Flow Estimation and Modeling of Engine Flows

An assumption of constant discharge coefficient (Cd) is often made when modeling highly compressible pulsating engine flows through valves or other restrictions. Similarly, orifices and flow nozzles used for real-time EGR flow estimation are often calibrated at a few steady-state points with one single constant Cd that minimizes the error over the selected points. This assumption is based on near constant Cd observed at high Reynolds number for steady flow. It has been shown in this work that this assumption is not reasonable for pulsating flow, particularly at large amplitudes and low flow rates. The discharge coefficient of a square-edged orifice placed in the exhaust stream of a diesel engine produced Cd’s varying between 0.60-0.90 for the resulting critical/near-critical flows. A novel pulsating flow measurement apparatus that allowed independent variation of pressure, flow rate and frequency and allowed reproducible measurements independent of transducer characteristics, produced Cd’s in the range of 0.25-0.60 with a similar square-edge orifice. The variation in Cd was characterized by two dimensionless variables that normalized the standard deviation of the pulsating signal with dynamic pressure and average differential pressure drop respectively. The results raise important questions that can potentially initiate fundamental work to fill the gap in the literature that exists for highly compressible pulsating flows.

this possibility was to measure the discharge coefficient (C d ) of an orifice located in pulsating flow. A preliminary study investigating the C d of a square-edged orifice installed in the exhaust stream of a diesel engine found that C d kept increasing with Reynold's numbers exceeding 10,000 for steady pulsating flow . Steady-flow C d ′ s are asymptotically constant beyond that Re. The current work was initiated as a result of these observations.
The five engine operating conditions in the preliminary work were expanded to 105 steady state data points for the current work. The engine pulsation parameters were correlated, e.g., pulsation amplitude increased with flow rate. To separate causative factors, a pulsating flow apparatus was constructed in which each factor could be independently varied. This allowed a systematic investigation of dimensionless groups that might be correlated with flow resistance.
A brief theoretical/historical background of steady-flow C d and its direct application to pulsating flow is provided below. This is followed by a review of the literature addressing C d measurement in pulsating flow. The literature on pulsating flow resistance has historically focused on sources of error when applying steady-flow knowledge to unsteady flow situations. Based on the literature and current work, the quasi-static assumption appears to be the biggest source of error. This assumption has also been questioned in related areas such as pulsatile pipe flow and the modeling of turbocharger performance; related literature has been discussed.
The discharge coefficient of an orifice or nozzle is defined as the ratio of the actual to theoretical flow rate: The theoretical flow rateṁ theoretical for steady compressible flow through an orifice or nozzle is given by Fluid Meters (1971a): Subscripts "1" and "2" refer to upstream and downstream conditions, respectively, p is the pressure, ρ is density, k is the specific heat ratio for the gas and A is the flow area. Equation (2) can be derived by assuming isentropic expansion through the restriction, and applying conservation laws of mass and energy.
The discharge coefficient C d accounts for friction, heat transfer, boundary layer thickness and the approximation that properties measured downstream of the orifice or nozzle exit represent the conditions at the exit. The theoretical flow rate can also be calculated by incorporating a theoretical or empirical "compressibility factor" Y in the incompressible orifice equation derived by applying the Bernoulli's equation across an orifice (Fluid Meters, 1971a): The theoretical compressibility factor can be obtained by rewriting Equation (2) in the form of Equation (3). The ASME research committee on flow meters first documented empirical equations and tables for C d for a variety of orifices, nozzles and nozzle-venturi devices in 1924, and continued updating results till it produced the sixth edition in 1971 (Fluid Meters, 1971a). An empirical equation for C d for a square-edged orifice in steady flow with flange taps for pressure measurement was used in the current work (Fluid Meters, 1971b). The report recommends damping out any pulsations prior to compressible flow measurement, based on the work of Sparks (1966), who found no relationship between the average differential pressure and the average velocity of pulsating air-flow through an orifice. The report recommends using surge bottles and pulsation dampeners according to the design recommendations of Chilton and Handley (1952). However, this is not always a practical proposition for internal combustion engines. Exhaust Gas Recirculation (EGR) measurement is often based on differential pressure measurement of pulsating exhaust gases. The current industrial practice is to assume that the discharge coefficient of the orifice/flow nozzle is constant at the high Reynolds numbers encountered in the EGR loop, and to experimentally determine C d in a test cell by running the engine at select steady-state operating conditions. This assumption is based on the data in the ASME report originating in the work by Bean et al. (1929), as well as work on viscous liquids by Tuve and Sprenkle (1933). The data was collected for steady-state, no-pulsating flows. Similarly, it is common practice to model engine flows, e.g., through intake and exhaust valves, using C d values determined from non-pulsating data. For example, the commonly used engine modeling software, "GT-Power, " uses discharge coefficients based on the work of Lichtarowicz et al. (1965), for modeling pulsating flow through all flow restrictions such as valves and orifices, using a quasi-static assumption. The current work demonstrates that the quasi-static assumption is incompatible with the range of pulsating C d 's measured for many flow conditions. For example many pulsating measurements at high pressure and low flow rates were in the range 0.25-0.40, well outside the range of 0.60-0.62 if quasi-static assumptions were applied to these flow conditions for the particular orifice used.
The existing literature on discharge coefficients in pulsating flow does not offer a solution to the practical problems mentioned above; rather it allows a better understanding of the problem. It discusses four main aspects, mainly in the context of the incompressible flow using Equation (3). a. Square Root Error: The error associated with averaging the pulsating differential pressure signal rather than averaging the square root of p.
b. Measurement Errors due to acoustic effects arising in the pressure lines and fittings leading to the differential pressure transducer. Slow response transducers and aliasing could compound this problem.
Square root errors, acoustic effects and aliasing are wellunderstood, e.g., Nakamura et al. (2005) and González et al. (2016) have resolved these errors when using Pitot tubes to measure flow rate at the exhaust tailpipe. c. Error due to neglecting the temporal inertia term: The unsteady Bernoulli's equation across an orifice can be integrated (along the direction of flow) to obtain the following equation, valid at any instant, for incompressible, inviscid flow (Doblhoff-Dier et al., 2011): where C 1 and C 2 are constants calculated from the geometry of the orifice and measurement setup. The second term on the right-hand-side is the temporal inertia term that accounts for the convective acceleration of the flow along the streamline, as a result of increasing pressure of the pressure pulse. The term is zero for non-pulsating flow, in which case Equation (4) reduces to Equation (3), with C 1 encompassing both C d and Y to accommodate viscous, compressible flows. For pulsating flow, since part of the pressure differential is used to accelerate the flow, neglecting the inertial term can result in overestimation of C d . In their 1991 review paper, Gajan et al. (1992) have reported that the inertial error increases with frequency (Strouhal number) and amplitude, and cannot be ignored; errors up to 25% have been reported (see Figure 3 of their paper). Doblhoff-Dier et al. have examined this question in a more recent 2011 paper (Doblhoff-Dier et al., 2011) and found that the inertial term can only be ignored if the ratio of the pulsating velocity component to the product of frequency and orifice diameter exceeds 10. Most of the data presented in this work exceeded that ratio, ranging between 1 and 200. Note that Equation (4) was derived from incompressible assumptions and a simple inertial term cannot be derived for compressible flows. Doblhoff-Dier et al. were able to ignore compressibility effects for their experimental setup, but since the density change across the orifice was quite significant for the engine data presented in this work (up to 100%), the inertial term of Equation (4) might not represent most engine flows. d. Quasi-Static Assumption: All of the above assume that the C d of pulsating flow is identical to that of steady state flow (quasistatic assumption), and accurate metering of pulsating flow can occur if there exists no square-root and measurement errors and inertia of the accelerating flow is accounted for. However, as pointed out by Doblhoff-Dier et al., many researchers (Gajan et al., 1992;Mottram, 1992;ISO, 1998) believe that pulsating C d (and therefore C 1 in Equation 4) is different from steady state C d , owing to differences in boundary layer thickness and flow structures in the separation zone. Gajan et al. (1992) have said that "the biggest gap in our knowledge of pulsation effects on flow through an orifice is probably the variation of the discharge coefficient with various pulsation and acoustic parameters in a given flow system." They have described work performed at CERT/ONERA (France) where it was found that the plateau shaped velocity profile at the exit of the orifice disappeared faster in the free stream direction for pulsating flow relative to steady flow, i.e., diffusion of momentum was increased by flow pulsation. Variations in C d were observed to increase with amplitude of pulsation, but could not be correlated to any particular variable, e.g., Strouhal number. Similar variations in pulsating C d have been reported by Reis and Hanriot (2017), very recently for low Reynolds number pulsating liquid flow. Deviation of C d from steady state values in accelerating and de-accelerating flows have been studied as early as 1951 (Daily et al., 1955).
In addition, substantial literature on pulsating flow through pipes (not orifices) exists, and the effect of oscillations on flow structure, boundary layer and shear-stress is a sub-field in itself, see for example relatively recent highly cited work by He and Jackson (2009), or review papers by Gundogdu and Carpinlioglu (1999a), Gundogdu and Carpinlioglu (1999b), Çarpinlioglu and Gündogdu (2001), Iguchi and Ohmi (1982), Ohmi et al. (1982), or Ahn and Ibrahim (1992).
These investigations are inconsistent with the quasi-static assumption and so are the results presented in this work. According to the quasi-static assumption, for the particular orifice installation and range of Reynolds numbers tested, all C d 's should be within a narrow range between 0.60 and 0.62, owing to asymptotic behavior at high Reynolds number. But most of the pulsating measurements were outside that range, and high pressure-low flow conditions produced C d 's in the range of 0.25-0.40. No combination of instantaneous quasi-static C d 's exceeding 0.60 can result in and average C d below 0.40. The inertial term is too small to explain this discrepancy. For all the data points in this work, Inertial force was a small fraction of the total pressure force driving the flow (<0.2%). For pulsating flow conditions producing low C d 's, most system-level engine simulation packages would overestimate flow rates even if they simulated the wave dynamics accurately, since they generally use the quasi-steady assumption to choose C d at the predicted instantaneous Reynolds number using steady flow results. For example, GT-Power uses the steady flow results of Lichtarowicz et al. (1965) to assign default (but adjustable) values of discharge coefficients for orifices of various lengths. If the low C d 's for pulsating flow are a result of differences in boundary layer and/or flow structures in the separation zone, they cannot be described by system level models. This is where experimental data such as that generated by this work can be useful. The same argument can be made for on-board EGR flow estimation.
The errors listed above are also encountered when predicting turbocharger performance with maps generated by steady flow (non-pulsating) test benches. But unlike discharge coefficients for pulsating flow, interest in unsteady turbocharger performance has continued to this day, and some insight can be drawn from this body of work. Automotive turbochargers experience pulsating flow outside the laboratory but their performance has historically been characterized by steady flow maps, using the quasi-steady assumption. However, several researchers have found significant deviations in the radial-inflow turbine characteristics under pulsating flow conditions and corresponding quasi-steady predictions (Dale and Watson, 1986). These researchers include Dale and Watson (1986) and Baines et al. (1994) at the Imperial college, Winterbone et al. (1990Winterbone et al. ( , 1991 at the University of Manchester, Kosuge et al. (1976) at Keio University in Japan and Capobianco et al. (1989) and Capobianco and Gambarotta (1992) at the University of Genoa in Italy. Most of these investigations have focused on different turbocharger turbines, e.g., single/double entry, mixed flow etc., rather than the compressor presumably because the pulsations are of much larger magnitude at turbine inlet. Kosuge et al. (1976) have reported increasing deviations with pulse frequency and amplitude (both were correlated and the effect of each factor was not studied separately). More recently, Copeland et al. (2011) have found a significant decrease in measured turbine efficiency for high amplitude-low frequency pulsating flow through a double-entry circumferentially divided turbocharger, relative to quasi-steady predictions. They have also summarized the insights gained from the previous works cited above: The emptying and filling of mass and the wave dynamics occurring within the turbine volute/passage determines the flow dynamics for pulsating flow. Further, the pressure fluctuates faster than the mass flow rate. Quasi-steady predictions would neglect these effects. Consequently, Liu and Copeland have recently developed a new method for turbine mapping that uses experimental data generated from a pulse rig (Liu and Copeland, 2018). They have also reported negative turbine power during the trough period of the pulsation. Marelli et al. (2014) have performed similar investigations on turbocharger compressors.
The disagreement between quasi-steady and dynamic behavior for both turbochargers and orifices probably stem from different phenomena. Boundary layer and flow separation differences between steady and pulsating flow might be more important for orifices. The effect of volume capacitance and wave dynamics on mass flow rate might be more important for turbochargers. But it is interesting to note that turbines have been continuously investigated for non-quasi-steady behavior since the 70's but not orifices. In fact, no work could be found in the literature that investigated highly compressible pulsating flows that occur in engines. None of the pulsating orifice works cited above had density changes of the magnitude encountered in engine flows. For the current work, the ratio of upstream to downstream density for the engine data ranged between 1.08 and 1.91. On the other hand, many researchers have investigated highly compressible (critical) but non-pulsating flows (Grace and Lapple, 1951;Jackson, 1963;Kastner et al., 1964;Deckker and Chang, 1965;Rohde et al., 1969;Brain and Reid, 1975;Ward-Smith, 1979).
While there exists a gap in literature, the ubiquitous use of variable geometry turbochargers and high pressure EGR systems have created a need for accurate modeling and estimation of pulsating flows across a wide range of flow conditions. The current work is not aimed at understanding pulsating C d in engine flows at a fundamental level, but rather to measure and characterize it across a wide range of frequencies, flow rates and pressure, including highly compressible pulsating flows and critical/near-critical pulsating flows that have never been investigated till now. Instead of trying to understand the observed deviation from steady flow C d in terms convective acceleration and non-quasi-static behavior, attempts have been made to characterize C d in terms of dimensionless parameters, such that the results are useful for engine modeling and control. In addition, the data can be used by other researchers for fundamental studies that are required to answer some of the questions originating from the results of this work.

DESCRIPTION OF APPARATUS AND TEST PLAN
A schematic of the engine apparatus is shown in Figure 1, top. An orifice was located in the exhaust stream of a naturally aspirated four-stroke diesel engine. High frequency pressure and temperature measurements were made at and across the orifice. The mass flow rate through the orifice was equal to the reference intake mass flow rate measured with a laminar flow element, plus the fuel flow rate measured with a rotary fuel flow meter. These measurements were used to calculate C d using Equation (2), for steady-state data points spanning the torque curve, ranging between 40 and 100% throttle position (increment 10%) and 1,000-2,900 RPM (increment 500 RPM). Additionally, a finer RPM sweep with 50 RPM increment was undertaken at 70% throttle position. Data from 51 unique speed-throttle combinations was obtained, and repeated several times. Since the dynamometer used throttle rather than torque control, every throttle repeat produced slightly different fueling. Including the repeats and after discarding data that failed quality checks, 105 unique steady-state points with different mass flow rates were obtained in all. Plots of the operating conditions are available in Figure A1 in Supplementary Material. A detailed description of the engine apparatus used can be found in previous work .
The engine data (presented later) was not ideal for investigating dimensionless groups correlated to pulsating C d . Amplitude, frequency and flow rate were correlated and could not be independently varied. Orifice pressure could not be increased without excess smoke/engine stalling. The pulsations produced by the four cylinders were of unequal magnitude as opposed to consistent pulsations produced by the apparatus, compare An apparatus was therefore constructed to obtain consistent pulsations while varying frequency, flow velocity and system pressure independently. These variations also produced differently shaped p pulses; the pulsating patterns of nine of the ten repeat points selected for uncertainty analysis are shown in Figure 1 to illustrate this. The pulsating patterns of all of the 189 data points used for this work is available in Figure A2 in Supplementary Material.
A sketch of the pulsating apparatus is shown in Figure 2. Pulsating air is generated by a reciprocating compressor and flows to the test section through a bleed valve (V1) to control mass flow rate and therefore flow velocity. Downstream of the orifice, the orifice pressure valve (V2) controls the system pressure. V1 and V2 are intended to control velocity and pressure independently; in practice it is a multiple inputmultiple output (MIMO) system in that system pressure FIGURE 1 | The schematic on the top shows the apparatus used to generate the engine data, details available in Brahma et al. (2014). Bottom plots: Comparison of engine differential pressure pulses (bottom-right) with those generated by the pulsating apparatus (bottom-left). Each engine cycle produced four pulses that were different from each other. Possible reasons are differences in combustion and different interference patterns dependent on the geometry of the exhaust manifold. The pulsating apparatus used a single cylinder compressor to produce consistent pulsations. is weakly dependent on V1 and flow velocity is weakly related to the position of V2. After V2, the pulsating flow is cooled to room temperature in a heat exchanger before being sent to a surge tank to remove pulsations prior to the reference Pitot-static tube measurement, made with steady air at room temperature.
The Pitot tube measurement is calibrated with non-pulsating steady flow shop air, which replaces the compressor during shop air calibration. Each component of the apparatus and its justification and evolution is discussed below.
A two horsepower "Makita" brand single cylinder compressor designed for industrial use was used to generate pulsating flow, shown in Figure 3E. It was connected to the main test section with 9.5 mm outer diameter (3/8th inch) copper pipe of length 1,168 mm. In order to vary pulsation frequency, the compressor's single-phase AC motor was replaced by a 3-phase AC motor whose speed was controlled by a variable frequency drive (VFD). A specially machined adapter made this possible. The VFD was used to vary the compressor speed between 1,000 and 3,000 RPM during testing, i.e., between 20 and 50 Hz. For a 4-cylinder engine this frequency range corresponds to 600-1,500 engine RPM while the engine data ranged between 1,000 and 3,000 RPM. The partial overlap of frequency was due to safety reasons. The compressor was built to run at 1,800 RPM only, and running it above 3,000 RPM produced significant vibration and thermal stresses. Pulsation frequency by itself was uncorrelated to C d (see Figure A3 in Supplementary Material); the main difference between the engine and apparatus data resulted from differences in flow velocities. It will be shown later (by Figure 12) that both datasets can be reconciled with two non-dimensional groupings. Shows the shell-and-tube heat exchanger used to cool the gas prior to the pitot-static reference measurement (D). The compressor used to generate pulsating flow is shown by (E). (F) Shows the surge tank and pressure regulation equipment used to provide steady flow for the shop air calibration process. Figure 3A shows the main test section. Two 14 mm ID (5/8th inch OD) copper pipes, 559 mm long, were soldered on either side of a hollow brass cylinder bored to the same internal diameter. The brass cylinder was split into two 38 mm sections that were bolted together with the aluminum orifice plate sandwiched in between. The orifice plate was therefore installed between approximately 597 mm straight pipe sections on both sides, corresponding 37 pipe diameters (ASME recommends at least 6 pipe diameters for this kind of installation (Fluid Meters, 1971c). The orifice diameter was 4.2 mm, so the ratio of orifice to pipe ID was 0.30 (d/D or β = 0.30), chosen to be close to the engine data (β = 0.26). Lower β's allow a better comparison between pulsating and non-pulsating flow since non-pulsating C d shows smaller variation and becomes constant at a lower Re d . The other dimensions of the system were selected by matching compressor performance (3.3 CFM at 621 kPa gage pressure) with measurement capabilities, with the selected β.
The concentric square-edged orifice was 3.2 mm thick (1/8th) and machined out of Aluminum according to ASME standards (MFC-3M standards). Flange taps for differential pressure measurement were located 25.4 mm (1 inch) from the orifice plate, and machined onto the brass cylinder. The orifice pressure measurement was teed off the upstream pressure tap. A taps was drilled 12.7 mm (½ inch) upstream of the orifice. An additional temperature surface measurement was made just upstream of the brass cylinder; an external thermocouple was attached with thermocouple tape and wrapped with insulating material (covering part of the brass cylinder in Figure 3A). The flow temperature measurement was used to calculate density; the surface temperature measurement was used for corroboration, to detect bad data.
Two copper pipes and fittings of combined length 222 mm were threaded into the flange taps to cool the gas by natural convection prior to the p measurement, to maintain the transducer at room temperature and prevent thermal drift. The transducers were laboratory grade variable reluctance transducers with interchangeable diaphragms to adjust range: Validyne Model DP 15 for orifice p (selected range 138 kPa, 0.5% FS) and upstream orifice pressure p 1 (selected range 552 kPa, 0.5% FS), Validyne Model DP103 for Pitot-Static p measurement (selected range 5.5 kPa, 0.25% FS). Vibrations were mitigated by the four threaded supports as well as hanging weights, seen in Figure 3A. Nylon tubing was used to isolate the transducers from test section vibrations. The total length of nylon tubing and cooling pipe between the p sensor and orifice was identical for upstream and downstream location.
Measurement of the pulsating pressure signal can be influenced by the frequency response of the sensor, data acquisition frequency, the natural frequency of the vibrating column of air existing between the diaphragm and the flow, and the vibration characteristics of the combined system. The frequency response of the Validyne transducer exceeded 10 kHz and data was acquired at 10,000 Hz. The natural frequency of the trapped column of air, including transducer cavity volume was calculated to be 68 Hz (see Supplementary Materials for calculations), which was higher than the highest signal frequency of 50 Hz, but not ideal since some amplification can be expected in the vicinity of the natural frequency. For systems with significant damping however, some distortion in the opposite direction (clipping of the peaks) can be expected due to the compressibility of the air trapped in the transducer cavity (∼0.07 cm 3 ), tubing (∼14.20 cm 3 ), and fittings (∼16.51 cm 3 including the convective cooling pipes). Figures A4, A5 in Supplementary Material show that this was the primary problem with the pulsating signal measurement. Ideally the transducer would have been mounted at the point of measurement without the cooling pipes to avoid these issues. However, this was not feasible due to the high temperatures resulting in thermal drift of the transducer and the high vibrations at that location. Attempts were made to cool the flow prior to the orifice but this resulted in condensation at some operating conditions, causing damage to the equipment. The transducer used for the engine data was located at the point of measurement and cooled by blasting shop air on it. This was not a desirable or practical solution for the pulsating apparatus due to the proximity of people as well as the reference pitot measurement near the transducer. Since the primary measurement required to calculate C d is the average rather than pulsating pressure signal, a decision was made to focus on an accurate average signal measurement and isolate the sensor from temperature and vibration. Damping tubes were fabricated (by stuffing a nylon tube with cotton, packing it tightly, with set screws on both ends to retain the compression) to obtain an accurate physically averaged measurement. Commercial snubbers were rejected because they could not eliminate the pulsations completely. Given the importance of transducer location, it is reasonable to believe that pulsating C d 's estimated from physical rather than electronic averaging is the only way to achieve agreement between different investigators. The pulsating signal was still recorded and used for characterizing the C d (by normalizing the standard deviation of the pulsating signal to create dimensionless groups) but this characterization was empirical and had no bearing on the accuracy of the fundamental C d measurements.
The dual capability of measuring damped and pulsating signal was achieved by using a high pressure valve that would bypass the damping tube when opened. Closed valve position allowed only damped pressure measurement, while open position allowed pulsed measurement. The system is shown by Figure 3B.
The damping tube and valve arrangement allowed physically averaged p andp 1 measurements, while also recording the pulsating p and p 1 at 10,000 Hz. The downstream pressure was calculated byp 2 =p 1 − p. Both measurements were made for every data point acquired in this work. As expected, differences between the average of the pulsating signal and the (physically averaged) damped signal were observed. This disagreement nearly disappeared during a limited number of experiments with the transducer located very close to the flow (and managing temperature by turning the compressor off intermittently). This confirmed the accuracy of the damped measurement. It also provided direct evidence that the trapped air column was clipping the pulsating signal (see Figure A5 in Supplementary Material).
A temperature measurement was necessary to at the Pitotstatic tube to calculate density for the reference mass flow rate calculation. This measurement needed to be non-intrusive and also account for the thermal inertia of the mass within the surge tank, i.e., the stabilization time prior to acquiring a new data point would need to be long enough to ensure that the temperature measurement was not influenced by the previous data point. Both these problems were eliminated by cooling the air to room temperature prior to entering the surge tank. A shelland-tube heat exchanger was used for this purpose, and is shown by Figure 3C. A transparent PVC pipe was used as the shell which enclosed the 15.9 mm OD (5/8th inch) copper pipe. Compression fittings (Swagelok brand) were used with tee fittings (see inset) to create a water-tight seal. Cold city water flowed through the shell continuously at a high flow rate exceeding 25 gallons per minute, and was drained without recirculation. The "U" shaped 15.9 mm OD (5/8th inch) copper pipe and surge tank downstream of the heat exchanger, both at room temperature, ensured that pitot exit was at room temperature. This was confirmed by a non-intrusive surface temperature measurement prior to exit to the Pitot tube shown by Figure 3D. Room temperature Pitot air eliminated the need of an intrusive thermocouple near the exit that could potentially interfere with the delicate Pitot-static measurement (about 305 mm H 2 O p at maximum flow rate). The heat exchanger was initially located upstream of the test section but this resulted in condensation and two-phase flow through the orifice. Shifting the heat exchanger downstream of the test section eliminated the condensation problem because cooling occurred at low pressure. Figure 3D also shows the mechanism (left) that allowed the Pitot-static tube to move in all three axes. An adapter was made to center the Pitot tube relative to the exit pipe at a distance of 3 diameters. The effective flow area, which, when multiplied by the velocity and density would equal the mass flow rate through the apparatus was determined by flowing non-pulsating steady shop air through the system. Orifice p and orifice temperature measurements were used to calculate the mass-flow rate during this steady flow process, henceforth referred to as "shop-air calibration." Shop-air calibration, i.e., calibrating the effective area of the exiting jet, utilized the fact that the C d for a squareedged orifice in steady flow is well-known. ASME recommended mathematical functions were used to determine C d (Fluid Meters, 1971b), which ranged between 0.60 and 0.62 depending on the Reynolds number. The shop air calibration process was therefore the reverse of regular data acquisition: instead of using the pitotstatic measurement to determine orifice C d in pulsating flow, the known C d of steady flow was used to determine the effective area to be used with the pitot-tube velocity and density. It was discovered that the effective area was not constant, but a function of Reynolds number. The ratio of the effective area to the area of the exit pipe (effective area ratio) ranged between 0.74 and 0.84 for Re < 10,000 and between 0.82 and 0.84 thereafter, with virtually no variation after Re > 25,000 (calibration curve included in Figure A6 in Supplementary Material). A calibration curve between the voltage output of the pitot transducer and the effective area ratio was used to determine the mass flow rate during pulsating data acquisition. Note that the pulsation made no difference to the Pitot tube measurement because it was downstream of the 14 gallon surge tank. The effective area ratio during the shop-air calibration process and the pulsating data acquisition can therefore be expected to be identical for the same Re (or the same Pitot Voltage, since other properties were at nearly identical room temperature and pressure).
Shop air at 965 kPa gage pressure was available, and although free of pulsations, its pressure tended to change gradually based on its usage by other machines in the building. A target of at most 0.05% difference over a 10 second time period, corresponding to each shop-air calibration data point, was set. To meet this target, another bigger 30 gallon surge tank, shown by Figure 3F, was used with two pressure regulators (upstream and downstream). The pressure was first stepped down to 689 kPa gage pressure upstream of the surge tank and then to 414 kPa gage pressure downstream of surge tank. The supply pressure was therefore maintained below the threshold above which building supply fluctuated.
The density at the Pitot-static tube and the orifice was calculated by treating the fluid as a mixture of dry air and water vapor. Humidity was measured with a sling psychrometer. Minute-resolution barometric data from a barometer installed on the roof of the building (adjusted for height of the instrument) was used to determine atmospheric pressure. Barometric data was corroborated by local weather station data points. Every data point during pulsating flow data acquisition as well as during the shop-air calibration process utilized the barometric pressure, ambient temperature and humidity measurements.
In order to make the system leak-proof, the pipe exit at the Pitot-static tube was plugged, and with bleed valve V1 closed, the apparatus was pressurized to 552 kPa, the maximum pressure during pulsating data acquisition. The pressure drop was measured as a function of time to calculate leakage rate. Leaks that were large enough to be detected with soap water were systematically eliminated by either tightening or changing the fittings or by re-taping the pipe threads with Teflon tape. At least four layers of Teflon tape were used for every fitting. The entire apparatus was redone at one point to include fewer fittings and change the quality and quantity of Teflon tape. The final hardware used to generate the data reported in this work had a leakage rate of 0.26 grams per hour, a number too small to affect the data, and too small to be detected by the soap-bubble method.
All the data was collected with National Instrument DAQs at a 1,000 Hz. LABVIEW was used to obtain 10 s worth of data at every data point (defined uniquely by compressor RPM, flow rate and valve V2 position, test plan described next), i.e., 10,000 instances of orifice pressure/temperature, orifice P, Pitot-Static P/temperature, and compressor outlet temperature over a 10 second time-period were collected for every data point. More than 1,500 pulsating data points were collected while the apparatus hardware evolved to be leak-proof, free of twophase flow (no condensation), allowed cooling to enable room temperature Pitot-static measurements, allowed a physically averaged as well as pulsating orifice P measurement, and allowed repeatable shop-air calibration with <0.05% variation in system pressure. The 189 data points presented in this work were acquired only after the hardware had evolved to its final form over a period of about 2 years.
The 189 point dataset (1,890,000 instances) spanned the entire range of the apparatus, covering variation in pulsation frequency (20-50 Hz), flow rate (1-10 kg/h corresponding to 5-30 m/s average orifice velocity), and system pressure (1-6 bar). This variation was achieved by changing the positions of the bleed valve V1 and orifice pressure valve V2 at different frequencies (see Figures 2A, 3 for location of valves). Frequency was varied in steps of 5 Hz from 20 Hz to 50 Hz. For each frequency, orifice pressure was set at three different levels with V2 with the bleed valve V1 closed. Then for every frequency and V2 setting, the flow rate was changed by adjusting the bleed valve V1, starting from fully closed to fully open. An average of 9 points was collected for each of the 3 V2 settings at each of the 7 frequencies, resulting in 7 × 3 × 9 = 189 points. The test plan is shown in Table 1. In order to record data at steady state, the apparatus was allowed to stabilize for at least 10 min when changing to a new frequency, 3 min after a change in V2 and 1 min after a change in V1. After the damped P was measured, both high pressure valves were opened to allow pulsation measurement, and 25 s were allowed to elapse before the pulsating signal was recorded.

Pulsating Apparatus Results
The results presented below were found to be repeatable within ±0.01 when measurements were made on 10 different occasions with different barometric pressure, temperature and humidity. Orifice density was calculated by treating the gas as an ideal mixture of dry air and water vapor. The specific heat ratio was calculated as a function of temperature using ideal gas tables for interpolation. The reference mass flow rate was calculated from the velocity calculated from the Pitot-static p measurement downstream of the surge tank. The density of the air at the Pitot-static tube was calculated using barometric pressure and measured temperature, which was maintained at virtually room temperature by the heat exchanger. The effective area of the flow at the Pitot tube was determined for every data point by using the calibration curve generated by the Shop Air Calibration process. Figure 4 (left) shows the discharge coefficient calculated for the 189 data points using the physically averaged p 1 and p 2 (using p) measurement with Equations (1) and (2). A large variation in C d from 0.25 up to 0.60 is observed. Figure 4 (right) shows C d plotted against the Strouhal number, commonly used to characterize pulsating flows. The Strouhal St number is defined as: where f is the frequency (Hz), d is orifice diameter and V is average pipe velocity at the orifice. As with Reynolds number, even though a general trend is noticeable, a large variation in C d values at any given Strouhal number is observed. The lack of correlation with Strouhal number has also been reported by CERT researchers (Gajan et al., 1992). Dimensional analysis was performed to find dimensionless groupings to explain the large variation of C d 's. No single dimensionless parameter could reduce the variation to below the general uncertainty level of ±0.01 (uncertainty analysis is presented later), hence pairs of dimensionless groups were investigated. The following pair of variables produced the best results, C d = C d (η , ξ ) where: Here, σ p is the standard deviation of the pulsating p (nondamped measurement), p is the physically averaged pressure using the damped tube, ρ is the density at the orifice and V is the average pipe velocity at the orifice. Both η and ξ are ξ are non-dimensional measures of pulsation, normalized by the upstream dynamic pressure V 2 and the average pressure drop across the orifice p, respectively. These two latter variables (the two denominators) were almost perfectly correlated for steady flow at high Reynolds numbers (since ρV 2 α YC d p from Equation 3), but not for pulsating flow; a plot is available in Figure A7 in Supplementary Material. Since both p and 1 2 ρV 2 are dominant terms of an integral momentum or energy balance across the orifice, and since they are not strongly correlated for pulsating flow, both are required to characterize pulsation. As either η or ξ is increased, flow resistance increased, i.e., C d decreased. This is shown by Figures 5A,B. The highest η and ξ were obtained at higher system pressures (larger pulsations achieved by closing V2) and low flow rates (lower velocity and pressure drop achieved by opening bleed valve V1). For similar velocity and pressure drop, higher system pressure resulted in larger pulsations.
The variation in C d at any given η could be explained by ξ . This is illustrated by Figure 5C which shows a quadratic surface C d = C d η , ξ , η 2 , ξ 2 , ηξ generated by linear regression with all the data plotted on the response surface (left subplot). In general, all 189 points were located very close to the quadratic surface.
The corresponding contour plot is shown by Figure 5D. Only those parts of the surface in the vicinity of data points have been shown; the rest of the response surface is extrapolation and not shown. The data occupies only a fraction of the two-dimensional η -ξ space because both variables are correlated.
The accuracy of the η -ξ surface was assessed by 10-fold crossvalidation. The widely used recommendations made by Weiss and Kulikowski (1991) for estimating the predictive capability of an empirical model built with between 100 and 1,000 data points, were followed. Nine-tenths of the randomized data was used to build the regression model which was then evaluated on the remaining one-tenth, and the process was repeated ten times, with different on-tenth splits chosen each time to cover all the data. Figure 6 shows predictions made with the regression model ( Figure 6A) and k-nearest neighbor (kNN) algorithm ( Figure 6B).
The kNN algorithm simply chooses k training data points that are nearest to the test point to be predicted (k = 1 for this work). It is independent of the shape of the response surface; hence it was chosen to compare parameter pairs that would best predict C d . The η -ξ pair performed better than all other possible pairs,  as shown by Table 2. All parameters are listed in ascending order of RMS error resulting from the 10-fold cross-validation process. More than two parameters did not produce any significant improvement in predicting C d . For regression, the quadratic surface is only an approximation. It performed better than a linear or pure quadratic surface (no interaction terms) bit it is not known if the "true" surface is quadratic. The regression errors could possibly be further reduced by a different surface, e.g., exponential surface, but such an exploration is not feasible with this data for reasons described later.
Two interesting observations were made. If the numerical average of pulsating signal was used to calculate C d (instead of the physical average), the predictions improved. Second, the pulsating signal used for calculating σ p was acquired at 1,000 Hz, but the results did not suffer significantly until the data acquisition frequency was reduced to 2 Hz, much below the pulsating frequency of 20-50 Hz. This suggests that a data acquisition frequency allowing accurate standard deviation estimation is sufficient. Why σ p is necessary, and whether it is sufficient to characterize complex pulsating patterns is not known.
The results presented so far are interesting starting points for fundamental investigations but not useful from a practical standpoint. This is because the determination of η requires knowledge of the average pipe velocity, which is unknown in most practical situations. In fact the usual purpose of using a FIGURE 7 | Left: Engine data showed a large variation in C d and suggests that Re d and C d are unrelated (though correlated) or that other additional variables determine C d . The lines are best fit lines at each RPM. The non-pulsating C d was calculated by using a mathematical function from the (Fluid Meters, 1971b) ASME report (Brahma, 2011). Right: The variation of C d is not explained by the Strouhal number either. discharge co-efficient is to determine the mass flow rate, which in turn is required to calculate velocity. An alternative is to use the velocity resulting from an isentropic expansion assumption: The p o ≈ p 1 approximation is justified for this dataset because none of the 189 data points had more than a 0.025% difference between p 1 p 2 and p o p 2 (and the dynamic pressure term did not exceed 0.05% of p 1 for any data point). Since the velocity calculated from Equation (8) corresponds to imaginary isentropic conditions, reasonable predictions could no longer be made with only η isentropic and ξ isentropic . It was found that both kNN and regression generated quadratic surface required the use of all dimensionless parameters i.e., C d = C d (η isentropic , ξ isentropic , Re isentropic , Ma isentropic , St isentropic ), unlike the situation with measured average velocity. Corresponding 10-fold cross-validation results for regression and the kNN algorithm are shown in Figures 6C,D, respectively. Both algorithms used the logarithmic values of the five dimensionless parameters, to reduce correlation between variables and avoid rank deficient matrices. As expected, the results are worse than those shown by Figures 6A,B where the measured average velocity was used to calculate η . Still, a practical estimate of C d within ±0.05 can be useful for a practical application if C d ranges between 0.25 and 0.60.

Engine Results
All the observations described above also held true for the engine data. Figure 7 shows the variation of C d with Reynolds number (left) and Strouhal number (right). The range of Reynolds number, about 40,000-160,000 is higher than the apparatus (5,000-40,000) due to the difference in the flow rates between the four cylinder engine and the single cylinder compressor. A wide range of C d 's ranging between 0.60 and 0.90 is observed. The higher values compared to the apparatus (0.25-0.60) might be due to the engine flows being in the critical or near-critical range. Similar numbers have been reported in the literature for steady (non-pulsating) critical flows (Grace and Lapple, 1951;Jackson, 1963;Kastner et al., 1964;Deckker and Chang, 1965;Rohde et al., 1969;Brain and Reid, 1975;Ward-Smith, 1979) in the 0.80-1.00 range. As mentioned before, no literature on pulsating critical flow could be found.
For steady flow, according to the equation/tables in 1971 ASME report (Fluid Meters, 1971b), the C d for the same diameter ratio β = d/D = 0.26 was virtually constant at about 0.60 for the same Re d range (40,000-160,000).
Although the discharge coefficient increased with Re d , large variations in C d at any given Re d can be seen, suggesting they are unrelated (but correlated), or other additional variables determine C d . Similar to the apparatus, C d could not be satisfactorily correlated to any one dimensionless group, but was could be described by η and ξ . Out of the 105 engine data points, 62 were critical, based on peak pressure ratio across the orifice. For these data points, ξ was normalized using p 1 rather than p as the denominator in Equation (7). show the surface and contour plots, respectively, for critical data points (non-critical points could not be plotted on this surface since ξ is defined differently). The surface and contours are qualitatively similar to the pulsating apparatus results shown by Figures 5C,D but more linear.
The accuracy of regression predictions is illustrated by Figures 8C,D. Both subplots show the results of 10-fold crossvalidation, using measured average velocity (left) and using isentropic velocity (right), corresponding to the regression subplots of Figures 6A,C, respectively. Therefore, the left subplot shows results for a quadratic η -ξ surface, while the right subplot corresponds to a quadratic multivariate surface, i.e., C d = C d (η peak , ξ peak , Re peak , Ma peak , St peak ) for predictions without knowledge of average velocity. Such a multivariate surface can be potentially useful for estimating EGR flow or modeling engine valves, where the average velocity is a-priori unknown. Predictive accuracy was within ±0.025; better than the pulsating data. Results from the 62 critical data points and 43 non-critical data points were combined for Figures 8C,D. Interestingly, the critical data when plotted by itself had significantly lower error (see Figure A8 in Supplementary Material). No clear criteria for classification of pulsating flow as critical or not could be found in literature. The peak pressure ratio of the pulsating signal was used to estimate Mach number at the orifice assuming isentropic expansion. This criterion and the implications of shifting it are discussed in a related work (Brahma, 2019) that has explored the possibility that pulsating flow might become critical earlier than steady flow due to acoustic effects.

Repeatability and Random Error
Ten data points spanning the range of C d values were chosen to determine repeatability of the data. These repeat points are overlaid over the 189 data points in Figure 9. The left subplot shows how the repeat points are uniformly distributed across the range of C d 's as well as the range of the primary dimensionless group η. Since C d and η manly depended on the p and σ p measurement, respectively, the repeat points were chosen to be concentrated around the lower ranges of the p-σ p space, where the repeatability was expected to be worse. This is shown by the right subplot.
Acquisition of these ten data points was repeated on five different days spanning a month. The orifice temperature, pressure, differential pressure and Pitot-static measurement were maintained at identical values for each of the ten points. The five repeats were then processed using two separate Shop Air Calibrations (for calibrating Pitot-static area factor), performed on two different days. This resulted in ten different repeat datasets, shown by Table 3.
Note that changes in ambient temperature, barometric pressure and humidity over the month-long repeat data acquisition were accounted for. Overall, the data was highly repeatable with the average standard deviation of nearly 0.005 in the value of C d . The maximum absolute deviation from mean for any of the 10 repeat points was about 0.01 while the highest percent standard deviation was 1.8% (for the smallest C d ).
The random uncertainty Px around the mean C d for each of the ten repeat points, i.e., for each row of Table 3, was calculated FIGURE 9 | The repeat points chosen for uncertainty analysis overlaid over the 189 data points. They were chosen to span the range of C d values, yet represent the lower range of the p measurement, where uncertainty is higher. Note that the relationship between C d and η appears linear on a logarithmic x-axis, as opposed to the non-linear relationship shown in Figure 5A. at a 95% level using the t-distribution according to standard practice (Wheeler and Ganji, 2009): Here, S x is the sample standard deviation, M is the sample size (M = 10) and S x √ M is the estimated standard deviation of the population.

Systematic Error and Total Uncertainty
Systematic error Bx was calculated for each of the 10, repeat points using: Here, w i 's are the different measurements, i.e., orifice pressure, differential pressure, orifice temperature, pitot pressure/temperature, while the E w i 's are the systematic errors. The derivatives were computed numerically by the MATLAB code used to process the raw data to final C d values, including code used for shop air calibration. Combined systematic errors were 0.5% of full scale for pressure measurements (Validyne Transducers) and 2.2 • C for temperature measurements (Omega type "k" thermocouples). Since the P term was overwhelmingly dominant, a more accurate determination of systematic errors for that measurement was undertaken. Non-linearity and hysteresis for each of the 10 P measurements were determined using a NIST traceable reference pressure gage (plot available in Figure A9 in Supplementary Material). The error of the reference pressure gage (0.25% full scale) was also accounted for. Therefore, Equation (10) had three separate terms for the p measurement and four separate terms for each of the other FIGURE 10 | Random Error, Systematic Error and Net Uncertainty for the ten repeat points.
four measurements. The total systematic error and the total uncertainty for all ten measurements are plotted and tabulated in Figure 10.
Based on this analysis it can be conservatively stated that total uncertainty is ±0.01 at any data point for which C d > 0.30, otherwise it is ±0.05. The average over all ten repeat points was about 0.01. The root cause for the higher uncertainties at lower C d is the low magnitude of p. It can be seen from Figure 10 that low C d 's occur at higher η and ξ . The denominators of both these variables become smaller at low p's. In theory, this uncertainty could be reduced by a factor of 2 by using a transducer with lower full scale value, but this would require two transducers of different ranges for the data presented here.

DISCUSSION
Overall, the RMS errors were generally comparable to the net uncertainty of about ±0.01 when measured average velocity was used to predict C d in terms of η and ξ (Figures 6A,B,  8C), the only exception being regression predictions of Figure 6A.Without prior knowledge of the average velocity upstream of the orifice (using dimensionless variables based on isentropic velocity), RMS errors (Figures 6C,D, 8D) were about ±0.02. As noted previously that the η -ξ surface might not be quadratic at all. An empirical exploration of the "correct" functional dependence could not be undertaken because any mathematical function that approximates the ratio of the two variables, i.e., P 0.5ρV 2 , will "fit" the data very well, since that ratio is the approximate definition of C d . For example, although a linear relationship between C d and log(η) was observed, any response surface with logarithmic variables produced near-perfect predictions, because the regression coefficients approximated the ratio P 0.5ρV 2 . The quadratic model form was chosen because such confounding was not possible, but the "true" mathematical form, if it exists, is unknown. Similarly, the role of the dimensionless variables Ma, Re and St is unclear. Each of them, used in addition to η and ξ , resulted in only a small improvement, and did not appear vital. But all five variables produced a significant improvement when used together.
Based on the combined data generated from both systems, four inferences can be drawn with reasonable certainty: FIGURE 11 | Combined data from engine and pulsating apparatus illustrate the greater range of C d relative to non-pulsating flow. a. Pulsating flow discharge coefficients are much more variable than non-pulsating cases. The current practice of using a constant C d at high Re d probably introduces significant error. This is evident from Figure 11 which shows data from both the pulsating apparatus and the engine on a single plot. For this combined data, C d ranges from 0.20 to 0.90 as Re d increases from 5,000 to 160,000.
b. The quasi-static assumption is not a good approximation: The quasi-static assumption predicts a C d ranging between 0.60 and 0.62. No combination of quasi-static values can produce most of the pulsating C d measurements shown in Figure 11. This suggests that the pulsating boundary layer and separation zone might be fundamentally different from steady flow phenomena. Hence a key assumption made by on-board EGR control strategies and system level engine simulation tools is not accurate and significant error is possible even if the pulsating wave dynamics is described accurately.
c. Pulsating critical/near critical C d 's are much higher than pulsating non-critical C d 's for about the same diameter ratio β, as evident from Figure 11, reasons unknown. A separate but related work examines the criteria for classifying pulsating flow as critical or not (Brahma, 2019), a question unanswered by literature because of the absence of work on pulsating critical flows.
d. The variables η and ξ are important variables for characterizing pulsating flow resistance. Figures 12A,B show C d to be strongly correlated to both η and ξ for the combined dataset. As mentioned previously, the engine C d 's are higher (and slopes are different) possibly due to critical/near-critical flow. Engine η and ξ are smaller relative to the apparatus data because of higher flow rates, velocities and pressure drops. Since there was significant uncertainty in the measurement of the pulsating signal, the characterization is possibly better than shown.
The original motivation of this work was to determine the root cause of inaccuracy in EGR flow and Volumetric Efficiency estimation during the turbocharger lag period, responsible for smoke spikes (Bean et al., 1929;Tuve and Sprenkle, 1933;Chilton and Handley, 1952). Based on this work it can be speculated that high pressure and low flow conditions encountered in the intake manifold and EGR system during the turbocharger lag period can possibly increase both η and ξ , thus increasing the flow resistance through engine valves, EGR valve and the airhandling plumbing, resulting in inaccurate ECM estimates of fuel-Oxygen ratio responsible for smoke spikes in electronically controlled engines.

CONCLUSIONS
The wide variation of discharge coefficients for a square-edged orifice in highly compressible pulsating flow was investigated and characterized by dimensionless groups. The work was motivated by inaccurate estimation of EGR flow and Volumetric Efficiency during the turbocharger lag period, based on constant discharge coefficient (C d ) assumption at high Reynolds numbers (Re d ).
Results from an orifice placed in the exhaust stream of a fourcylinder diesel engine showed that C d varied between 0.60 and 0.90 for Re d ranging between 40,000 to 160,000 for critical/nearcritical flows, typical of engines. Although correlated with C d , Re d did not explain the large variation in C d . A pulsating flow apparatus was constructed to further investigate pulsating flow C d under more controlled conditions with the ability to independently vary frequency, flow rate and system pressure with a wide range of consistent pulsating patterns. Based on experiments, physical rather than numerical filtering was performed to obtain a mean p signal, resulting in repeatable and reproducible C d 's independent of the transducer, transducer location and plumbing. Again a wide variation was observed, with C d ranging between 0.25 and 0.60 for Re d ranging from about 5,000-40,000.
For the particular orifice tested, the quasi-static assumption predicts a C d ranging between 0.60 and 0.62. Since the majority of pulsating C d measurements were well outside that range, and given that the inertial effects were estimated to be small, the quasi-static assumption does not hold true for the data. This suggests that the pulsating boundary layer and separation zone might be fundamentally different from steady flow phenomena. A key assumption made by one-dimensional engine simulation tools might not be accurate. Accurate prediction of the wave dynamics through the engine system might not sufficient to determine the flow resistance of the system if quasi-steady assumptions are used.
The variation of C d for both datasets was best explained by a pair of dimensionless variables η and ξ , defined as the standard deviation of the pulsating differential pressure signal across the orifice, normalized by the dynamic pressure and mean differential pressure, respectively. These variables were able to predict C d with errors approaching experimental uncertainty, for both engine and pulsating apparatus data. The response surface for both datasets was qualitatively similar. For practical situations where the average velocity used to calculate η is unknown, reasonable predictions were made by using the isentropic velocity to normalize σ p , although more than two dimensionless variables were required. There is room for improvement in the measurement of σ p , by locating the transducer closer to the orifice and minimizing the volume of connecting tubing.
The flow resistance of cylinder valves, EGR valve, flow measurement orifices/flow-nozzles (if any) and engine plumbing might depend on variables similar to η and ξ . ECM estimates of EGR fraction and volumetric efficiency for emission control could benefit from a variable flow resistance assumption.
Fundamental studies involving the visualization and modeling of the pulsating boundary layer through an orifice are required to better understand the dimensionless variables η and ξ , their functional relationship with C d , the role of St, Re and Ma, and the precise contribution of convective acceleration. The data from this work might be useful for this purpose.

AUTHOR CONTRIBUTIONS
IB designed experiments, collected data and performed analysis.

ACKNOWLEDGMENTS
The author would like to thank Bucknell machinists Aaron Clark and Daniel Johnson for design and fabrication of the pulsating apparatus, undergraduate student Kevin VanDelden for apparatus design and initial data collection, and Bucknell laboratory director Hugh Weber for assembling the Variable Frequency Drive and assistance and guidance with fittings and measurement. The author is grateful to the mechanical engineering department at Bucknell University for providing the resources required for this work.