rev. 11/89

Note: At this time NCAR operates only one aircraft, a C-130Q Hercules (tail number N130AR). NCAR operated the three aircraft mentioned in this Bulletin when it was written in 1989. The basic equipment and instrumentation, equations and principles discusssed here are still pertinent and being used at our facility.

**INTRODUCTION**NCAR presently operates three aircraft equipped with inertial navigation systems (INS) and a gust probe or a radome with pressure ports for air motion measurements. These aircraft (King Air, Sabreliner, and Electra) can be instrumented to measure both mean horizontal wind and turbulent fluctuations of the three air velocity components. The components are obtained by subtracting the velocity of the airplane, measured by the INS, from the velocity of the air, measured by the radome or gust probe. Since the flow measurements are in the aircraft frame of reference, these measurements are first converted to the earth's frame of reference by using the three airplane attitude angles measured by the INS, and are then subtracted from the airplane velocity components, which are obtained from integrated accelerometer outputs (corrected for Coriolis accelerations) in an earth-based frame of reference. Most applications of air-motion sensing systems involve either

(1) measurement of the mean horizontal wind, which is useful mainly for synoptic or mesoscale meteorological experiments, or

(2) measurement of high-frequency fluctuations of all three velocity components, which are used mainly to calculate turbulence statistics such as variances, and turbulence fluxes of heat, moisture, momentum, and turbulence energy.**AIR VELOCITY EQUATIONS**In this section we derive the equations for calculating the air velocity components from the INS and the flow-sensing probes. The derivation is similar to Axford's (1968). However, the terminology and coordinate systems are modified to conform to standard aeronautical engineering usage (Etkin, 1959) for the aircraft coordinate system and to standard inertial navigation (Kayton and Fried, 1969) and meteorological usage for the air velocity in the local earth coordinate system.

The velocity of the air with respect to the earth,

, is obtained by adding the velocity of the aircraft with respect to the earth,*V*= i*u*+ j*v*+ k*w*, and the velocity of the air with respect to the aircraft,*v*_{p}. If all the measurements are made at the same point, we have*v*_{a}The components of

are measured with sensors mounted on the aircraft, either on a boom forward of the aircraft to reduce the effects of the distortion of the airflow by the aircraft or alternatively, from pressure ports on a radome to take advantage of the pressure differences produced by asymmetric flow about the aircraft's nose. On the NCAR aircraft the components of*v*_{a}are obtained from the integrated outputs of accelerometers mounted on the stabilized platform of an INS. Alternatively,*v*_{p}could be obtained from radiation transmitting and receiving devices, such as Doppler radars or radar altimeters, or from radio or radar navigation techniques.*v*_{p}If the aircraft velocity is obtained from integrated accelerometer outputs, the velocity is measured in an inertial frame of reference. In an earth-based coordinate system, we have

where

is the measured acceleration

and are the angular velocities of the earth and platform, respectively; and

is the gravitational acceleration (which includes the centripetal acceleration).Since the accelerometers may not be collocated with the air velocity sensors, it may be necessary to add the term to (2.2), where is the angular acceleration of the aircraft and

is the distance from the accelerometers to the air velocity sensors. Thus, (2.1) can be combined with the integral of (2.2) and rewritten as*R*=*iX*+*jY*+*kZ*In the general case, the measured components of

and*v*_{a}may not be in the local earth coordinate system. The airflow sensors, for example, are usually rigidly connected to the aircraft fuselage. Therefore, it is necessary to know the angles between the coordinates of the airflow sensors and the local earth coordinates. The measured velocity components can then be rotated by means of the appropriate angular transformation to the local earth axes so that the rotation angles of the platform gimbals are equivalent to the aircraft attitude angles. From the inside platform gimbal outward, i.e., from the stabililzed accelerometer and gyro cluster outward to the aircraft frame, the order of the rotations is as follows: first, a rotation , about the*v*_{p}-axis of the earth coordinate system; second, a rotation , about the*z*-axis, which is the y-axis rotated in the horizontal plane by the angle ; and third, a rotation , about the*y'*-axis, which is the*x'*-axis in the aircraft coordinate system. The aircraft axes are shown in Fig. 2.1; a right-handed angular rotation is positive. The symbols and axes have been defined to correspond to standard aerodynamic conventions (Etkin, 1959).*x***Fig. 2.1. Top: Coordinate systems used in deriving equations for calculating the air velocity components. Bottom: Airplane attitude angles and axes used in equations for calculating the air velocity components.**In matrix notation, the transformation of the air velocity from aircraft (

,*x'*,*y'*) to platform (*z'*,*x*,*y*) coordinates {same as (b) coordinates in Fig. 2.1} is given by*z*where

The terms of must also be transformed to the local earth coordinate system.

is more than an order of magnitude larger than*X*and*Y*on all the NCAR aircraft installations, so that terms involving*Z*and*Y*can be neglected. Setting*Z*, we have*X*=*L*Similarly for we have

Thus, in vector notation,

The angle of attack, , is the angle of the airstream with respect to the aircraft in the aircraft's vertical plane, with positive in the downward direction; the angle of sideslip, , is the angle of the airstream with respect to the aircraft in the aircraft's horizontal plane, with clockwise (looking from above) rotation positive, as shown in Fig. 2.1. The air velocity components with respect to the aircraft are calculated by first correcting the measured true airspeed for angle of attack and sideslip sensitivities (See Section 4.); the magnitude of the corrected true airspeed is defined as

. We then rotate the airflow vector to the airplane coordinate system. If we first rotate through the angle to the horizontal plane, then rotate to the longitudinal aircraft axis, the components of air velocity in the aircraft frame of reference are*U*_{a}along the*-U*_{a}D^{-1}-axis,*x'*, along the*-U*_{a}D^{-1}tan-axis, and*y'*along the*-U*_{a}D^{-1}tan-axis, where*z'**D = (1 + tan*^{2}+ tan^{2})^{1/2}[This differs from the expressions given by Lenschow (1972). For small flow angles, however, the magnitudes of the differences are negligible. We are grateful to A. Weinheimer and J. Leise, who each independently pointed out the error and corrected the previous derivation.]

So far, the equations have been derived using standard aeronautical engineering terminology. In converting to a meteorological and inertial navigation frame of reference, we make the following changes: first, the sign of vertical velocity is changed so that

is positive in an upward direction; second, , which now corresponds to the true heading, is measured from the*w*-axis, or from north, and the*y*-axis now points east, as shown in Fig. 3. As a result,*z***(new)**=**(old) - 90°**,, and*v*(new) = -*v*(old). The final equations used for calculating the air velocity with respect to the earth are:*w*(new) = -*w*(old)These are the exact equations that are used in the standard NCAR data processing routines for calculating the air velocity. For purposes of illustration, however, it can be shown that for approximate straight and level flight (i.e., when the pilot or autopilot attempts to keep the aircraft level, but air velocity fluctuations still cause perturbations in the aircraft velocity and attitude angles) many of the terms in (2.4) are negligible.

First we assume that terms involving the separation distance,

, are small. This is generally true if the separation distance is less than about 10 m and the aircraft is not undergoing a pilot-induced pitching maneuver. Next we make the following small-angle approximations: for the horizontal components,*L*and*u*, we assume that the cosines of , and are unity; that terms that involve the products of sines of two of the above angles are negligible; and that*v***tan (, ) sin (, )**. For the vertical component, we make the same assumptions for the angles , and . Equations (2.10) then reduce to a form commonly used for approximate calculations of the three velocity components,*w*We can estimate the required accuracy of attitude and airstream incidence angle measurements from (2.11). In order to obtain the mean horizontal wind, the angle ( + ) must be measured accurately. For a mean wind accuracy of 0.5 m/s and an aircraft speed of 100 m/s, this angle must be measured with an absolute accuracy of 0.005 radians, or 0.06°. A heading angle of this accuracy is difficult to measure with a standard aircraft magnetic compass. An INS, however, can measure true heading to well within this accuracy.

A reasonable figure for short-term velocity accuracy necessary to estimate turbulence fluctuations is 0.1 m/s. At an airplane speed of 100 m/s, the required angular accuracy for the first terms on the right side of (2.11) is 0.001 radians, or 0.060°. This can be a difficult requirement to meet. Fortunately, an INS has even more stringent accuracy requirements than this for accurate navigation, so attitude angles of this short-term accuracy can be measured.

**DESCRIPTION OF INSTRUMENTATION**__Gust Probe Sensors__The three air velocity components with respect to the airplane are obtained from measurements of the two airflow angles and the true airspeed. Figures 3.1 and 3.2 show the two types of sensors that are used on NCAR aircraft to measure the airflow angles; the differential pressure probe and the instrumented radome sense the pressure difference across a symmetric set of ports at varying how angles. On the differential pressure probe, airspeed is obtained from a differential pressure measurement across the dynamic and static ports of the probe.

**Fig. 3.1a. The NCAR Sabreliner aircraft equipped with both the differential pressure probe (at the end of the boom) and the radome with pressure ports. Both are used to obtain flow angles by pressure difference measurements.****Fig. 3.1b. Schematic representation of the differential pressure probe manufactured by Rosemount Engineering, Inc., Model 858. Pressures at the ports are labeled on the figure.**- Radome with fluted deiced area
- Differential pressure ports
- Nose weather radar
- Sensor mount pods
- Pressure transducers

**Fig. 3.2. Schematic views of radome differential pressure gust probe system.**__Differential Pressure Probe and Radome Technique__An inevitable consequence of aircraft flight is distortion of the how in the vicinity of the aircraft. Therefore, any probe mounted on the aircraft is sensing in a perturbed flow field (Wyngaard et al., 1985). This means that for air velocity measurements using a boom-mounted gust probe, corrections must be made for the effects of this perturbation. These corrections are determined from the flight maneuvers discussed in Section 4. Alternatively, in the case of the radome technique, the pressure field resulting from the flow perturbations induced by the aircraft itself is used to obtain the flow angles. As discussed by Brown et al. (1983), for a hemispherically-shaped radome, this pressure difference is a nearly linear function of the flow angle.

The differential-pressure flow-angle sensor is constructed in the shape of a hemisphere capping a cylindrical tube (Fig. 3.1b). The flow around it leads to the pressure distribution shown in Fig. 3.3 when the tube axis is parallel to the flow. If the flow is at an angle with respect to the tube axis, the pressure distribution is no longer symmetric with respect to the tube axis. Thus, for a pair of ports symmetrically located above and below the tube axis, as in Fig. 3.3, an angular rotation causes an increased pressure in one port and a reduced pressure in the other. The net effect is a differential pressure between the top and bottom ports which is a function of flow angle and indicated airspeed. By placing a second pair of ports at right angles to the first, simultaneous measurements of both airflow angles with respect to local flow are obtained. Further details of the differential-pressure measuring technique are discussed by Brown et al. (1983).

**Fig. 3.3. Top: pressure distribution normalized by the dynamic pressure on the surface of a hemisphere-capped cylinder embedded in a uniform flow parallel to the axis of the cylinder. Bottom: A schematic cross-sectional view of the cylinder shows the arrangement of pressure sensing ports for sensing airflow angles, and Pitot and static pressures.**RAF currently has available an electrically-deiced differential-pressure probe (Rosemount Model No. 858AJ) for use on the Sabreliner aircraft. This probe can be used for limited measurements in clouds and in icing conditions. However, long flight legs in clouds or precipitation can lead to water accumulation in the pressure lines, which can cause errors in the pressure measurements. The probe has symmetric sets of holes drilled in the 2.8 cm-diameter hemisphere at an angle of 45° and a forward-looking port 0.14 cm in diameter. The probe also has a ring of static ports surrounding the cylindrical part of the probe 8.9 cm from the probe tip. Pressure differences are sensed by transducers mounted at the base of the nose boom connected to the ports by pressure lines. The inside diameter of the tubing inside the probe is 0.36 cm; the inside diameter of the lines connecting the probe tubing to the pressure transducer is 0.43 cm. The length of these pressure lines limits the instrument's frequency response. In the Sabreliner, for example, the tubing length is about 5 m. The time responses of both the dynamic pressure and the flow angle measurements on the Sabreliner with the differential pressure probe are discussed in Appendix A

The radome technique is similar, in principle, to the differential-pressure technique. The major difference is that the nose of the airplane itself is used as a probe instead of a separate flow angle sensor. Therefore, a nose boom is not needed for this technique. The technique cannot be calibrated in commonly-available wind tunnels, except possibly by using scaled models of the aircraft. In situ calibration techniques are used, and the system must be calibrated for each aircraft and whenever the nose configuration of the aircraft is significantly modified.

The radome system is capable of higher-frequency response than the differential pressure probe, since the transducers can be mounted closer to the pressure ports. The frequency response of the radome system used on the Sabreliner is > 10 Hz. Details of the radome technique are discussed by Brown, et al. (1983). All three NCAR aircraft can be deployed using the radome technique. The Sabreliner radome is unheated and therefore is limited to operation in clear air and non-precipitating cloud. The King Air and Electra radomes are deiced and contain provisions for trapping ingested water; therefore, they can be used in all weather. Wind tunnel, flight test and theoretical analyses indicate that, for both the radome and the differential pressure probe technique, the flow angles are related to the pressure difference across the respective angle-sensing pressure ports. For the attack angle, , this relation is

where

is the Pitot-static pressure difference, and is a bias in the measured flow angle. The sensitivity coefficients can be estimated from measurements of at known values of*q*, as well as or . For the Rosemount 858AJ,*q***=**0.079/° for Mach number,*C*< 0.52. This coefficient has been determined theoretically, as well as from wind tunnel calibrations conducted in the Cornell Aeronautical Laboratory (Andersen, 1955), and by the RAF in the U.S. Air Force Academy Subsonic Wind Tunnel. For*M*> 0.52, data obtained from the manufacturer indicate that*M*monotonically decreases to 0.056/° at*C*= 0.9. A fourth-order polynomial is used to approximate*M*for 0.52 <=*C*<= 0.90.*M*The sensitivity coefficients for the radome differential pressure system have been determined from theoretical studies, flow model calculations (for the Electra), wind tunnel scale model studies (for the KingAir) and from the flight maneuvers discussed in Section 4, and are generally in close agreement; however, the flight-derived coeflicients are used in practice for calculating the three air velocity components. These values are listed in Table 1, along with the tube lengths between the pressure ports and the transducers, and the angles and , which are the angles between the ports and the longitudinal airplane axis in the vertical and horizontal planes, respectively.

Table 1

Sensitivity coefficients, pressure port attitude angles ( for the attack angle and for the sideslip angle) and tube lengths for the radome differential pressure systems on all three NCAR aircraft. In all cases the tube diameter is 4.3 mm.Aircraft /° /° Tube length (m) King Air

(N312D)0.08207 33° 0.07448 33° 0.7 Sabreliner

(N307D)0.0570 33° 0.0410 33° 1.8 Electra

(N308D)0.07155

= 0.409533° 0.06577

= 0.037529° 3.5 The differential pressure transducers are calibrated on the airplane by applying known values of pressure difference across the pairs of probe ports.

__Pitot and Static Pressure Probes__All NCAR airplanes have dynamic or ram and static pressure ports for measuring dynamic (total - static) pressure

which is used to calculate airspeed. Rosemount Models 850L-2 and 852J Pitot-static tubes are mounted on the King Air and Electra wingtips; furthermore, the Rosemount Model 858AJ differential pressure probe (Fig. 3.1) already contains the necessary static and dynamic pressure ports. Dynamic pressure can also be measured at the front of the radomes, although they are more sensitive to airflow angles than the Pitot-static tubes. These pressure ports are connected to pressure transducers, which are used to measure both the static pressure*q*, as well as*P*_{s}. As discussed in Appendix A, the 858AJ probe has longer pressure lines which limit its frequency response for pressure fluctuations to < 1 Hz on the Sabreliner.*q*The true airspeed is calculated from

,*q*and the measured temperature,*P*_{s}, by application of Bernoulli's equation for the conservation of energy in a compressible fluid. (See Appendix B for details.) The Mach number,*T*_{r}(which is the ratio of the true airspeed to the speed of sound), is calculated from the relation*M*where =

/*C*_{p}is the ratio of the specific heat of air at constant pressure to the specific heat at constant volume. Then the true airspeed is obtained from the relation*C*_{v}where

(*r*r -*T*_{}s)/(*T*_{}-*T*_{t}) is the recovery factor of the thermometer,*T*_{s}is the gas constant (= 287.04 J/kg/K for dry air),*R*is the ambient air temperature, and*T*_{s}is the total air temperature which is the temperature attained by an adiabatic reduction of the air velocity to zero (with respect to the airplane).*T*_{t}

__INS Measurements__The inertial navigation system (INS) consists of an orthogonal triad of accelerometers mounted on a gyrostabilized gimbaled platform which is kept level by continuous computation and application of the required torquing rates. These rates are determined by integrating the horizontal accelerometer outputs and applying corrections based on the earth's rate of rotation and the geographic coordinates of the platform. The geographic coordinates are obtained by a second integration of the accelerometer outputs, again with appropriate corrections. The two corrected horizontal velocity components are available from the INS computer. The vertical velocity is obtained from a counter that accumulates pulses from the vertical accelerometer. Thus the INS can provide measurements of the three airplane velocity components and, by recording the output of angular resolvers mounted on the platform gimbals, the three attitude angles of the airplane with respect to the earth.

The INS platform used by NCAR (Litton Systems, Inc. LTN-51) is not torqued about its vertical axis while operating in the navigation mode, and therefore the heading or azimuth of the platform remains constant in inertial space. In this mechanization, which is called a free-azimuth, or azimuth-wander system, platform heading in the earth-based coordinate system changes with time and with east-west displacement of the aircraft. Its initial heading is determined by the orientation of the INS during alignment; early in the alignment sequence the platform is caged so that the platform axes are aligned to the INS case. True heading, therefore, is a calculated variable, obtained initially by a computer gyrocompassing of the system. Thereafter, during navigation the INS calculates and continuously updates the true heading. In the LTN-51, the angle between true north and the

-component of the platform (along which is measured) is called the azimuth-wander angle, , as shown in Fig. 3.4. The airplane velocity components in an earth-based (geographic) coordinate system are then calculated from the rotation*x***Fig. 3.4. Angular orientation of an azimuth-wander INS platform with respect to geographic coordinates.**The INS computer solves and integrates (2.2) in the horizontal plane of the local earth coordinate system; thus the instantaneous velocity and position are computed and can be displayed, along with related quantities, such as ground speed and drift angle, and range and bearing to selected targets or destinations. A block diagram of an inertial navigation system is shown in Fig. 3.5.

**Figure 3.5. Block diagram of an inertial navigation system. The dashed rectangle represents the gimbaled stable element which is kept level with respect to local earth, and R**_{p}is the position vector.The inertial system's computer also calculates the angular velocity required to keep the platform level. Torque is then applied (by electromagnetic coupling) to the appropriate gyro axes perpendicular to the spin axes of the gyro. The angular velocity of the gyro rotor along an axis mutually perpendicular to both the gyro spin axis and torque input axis is directly proportional to the applied torque. The axes of the stable platform are slaved to the gyro output axes by rotating the appropriate gimbals so that the angular deviation between the two sets of axes is negligible. Thus, the rotation angles of the gimbals are identical to the attitude angles of the aircraft in the local earth coordinate system. These angles are measured with dual-speed resolvers which have both single and eight-pole outputs. The outputs of each of the resolvers are thus approximately proportional to the sine and cosine of the attitude angle (for the single-pole output) and eight times the attitude angle (for the eight-pole output). Using the pitch angle as an example,

**/n = tan**, where^{-1}(E_{1}/E_{2})is the number of poles and*n*and*E*_{1}are the sine and cosine outputs, respectively. In practice, the single-pole output is used to determine the octant of the eight-pole output, and the eight-pole output is used for the attitude angle measurement.*E*_{2}A complete error analysis of an inertial navigation system is complicated. However, several salient features can be easily demonstrated by considering a simplified two-dimensional case. As an example, we investigate the error propagation of an accelerometer bias error or an angular misalignment of the platform and an accelerometer or gyro drift error, assuming that the Coriolis terms of (2.2) can be neglected. The accelerometers along the horizontal axes sense an acceleration approximately equal to the angular deviation from the local geopotential surface , times

, which is the gravitational acceleration at some reference level,*g*_{0}, from the center of the earth. The inertial system's computer integrates this acceleration twice, changes its computed position a distance, ; and rotates the platform through an angle*R*_{0}**/**. Thus, the acceleration sensed by the inertial system is*R*_{0}where is the acceleration error.

Along the vertical axis, the gravitational acceleration as a function of the vertical displacement, , away from some reference level and the latitude can be obtained from the equation

According to Kayton and Fried (1969), this equation is accurate to within 0.02 cm/s

^{2}at sea level. The resulting equations for horizontal and vertical position errors (neglecting the latitudinal variation of) are*g*where = = 1.24 × 10

^{-3}rad/s (equivalent to a period of 84.4 min, known as the Schuler period). The major effect of the neglected Coriolis terms involvingand*u*_{p}, which are typically about 10% of , is to divide the horizontal error into two sinusoidal components with angular velocities of approximately ± (1/2)( + sin ).*v*_{p}Thus, the two sinusoidal components added together result in a modulated error with a carrier frequency of and a modulation frequency of (1/2)( + sin ).

As an example, we assume

The solutions of (3.8) for the velocity error are

As can be seen, the vertical velocity error diverges rapidly. To compare the two errors, we assume

= 2 × 10*a*_{0}^{-2}cm/s^{2}and= 2 × 10*a*_{1}^{-5}cm/s^{3}. Substituting these values into (3.9), we obtainThese variables are plotted in Fig. 3.6.

**Fig. 3.6. Velocity errors in an inertial navigation system. (The error-generating function is described in text.) Curve A-vertical velocity, curve B-horizontal velocity.**The position errors, i.e., the integrals of (3.10) are

These variables are plotted in Fig. 3.7.

**Fig. 3.7. Position errors in an inertial navigation system. (The error-generating function is described in text.) Curve A-inertial altitude, curve B-horizontal position.**Other errors in position and velocity exist as well. Generally, the INS errors can be approximated by a monotonically-increasing, time-dependent error and a periodic error at the Schuler frequency. These errors can change drastically from one flight to another on the same system since unpredictable errors occur in INS gyros and accelerometers, and the initial alignment, which varies from flight to flight, plays an important role in subsequent INS performance. The alignment sequence, which normallv takes 15 to 20 minutes, includes initializing the INS coordinates, leveling the platform, and gyrocompassing to determine true north. Once the INS is in its navigation mode, any interruption can cause it to lose its navigation capability for the duration of the flight. Much more detailed analyses of inertial systems and their errors are presented by Broxmeyer (1964), Kayton and Fried (1969) and Britting (1971).

__Baro-inertial Loop__In the vertical, no closed loop exists for inertial measurements. However, pressure altitude can be used to limit errors introduced by errors in vertical acceleration. (An offset, or bias in vertical acceleration, when integrated twice, results in an error in altitude that increases with time squared.) At NCAR, the pressure altitude is combined with the vertical acceleration in a third-order loop (see Fig. 3.8) to provide an absolute long-term reference for altitude and vertical velocity. The advantage of a third-order baro-inertial loop (Blanchard, 1971) over a second-order loop system is the elimination of velocity and position error due to bias in the measured acceleration. For the NCAR aircraft (Lenschow, et al., 1978a), a 60-second time constant is used in the mechanization of this loop. This value is based on a compromise between minimizing the effect of high-frequency noise from the pressure altitude measurement (e.g., data system resolution, transducer noise, and pressure line "ringing"); and (1) minimizing the recovery time of the loop for extraneous disturbances (i.e., glitches) and (2) improving the long-term stability of the vertical accelerometer output.

**Fig. 3.8. Block diagram of a third-order baro-inertial loop for calculating altitude and vertical airplane velocity from a pressure altimeter and INS vertical airplane acceleration. Feedback loops involving C**_{1}, C_{2}, and C_{3}are set for a time constant of 60 s on NCAR aircraft. The term involving 2/*g*is a correction for the change of gravity with altitude (R is the distance to the center of the earth).*R*

**IN-FLIGHT CALIBRATION MANEUVERS**Regardless of where instruments are located and how carefully they are calibrated, errors are likely to be present in the measured variables. Ground tests are not useful for calculating velocity-related errors. Wind tunnel tests are difficult and prohibitively expensive for exact simulation of flight conditions. Therefore, in-flight calibrations play an important role in estimating errors and correcting aircraft measurements.

Because of upwash ahead of the aircraft, the airflow angles (attack and sideslip) and airspeed measured at the tip of a nose boom may differ considerably from the actual values that would be measured far away from the aircraft. The upwash affects not only the sensitivity, but also the zero offset of angle measurements, which, therefore, must also be determined from in-flight calibrations.

Maneuvers used for this purpose involve changes in aircraft speed and attitude angles. The following list summarizes several maneuvers used on NCAR aircraft equipped with an INS and the information that can be obtained from them; examples of these maneuvers are shown in Fig. 4.1:

**Fig. 4.1a. Examples of reverse heading and airspeed maneuvers used to check the quality of air velocity measurements. The lateral and longitudinal velocity components are measured with respect to the aircraft; therefore, the measured wind should change sign, but not amplitude, after the 180° turn, if the wind field remains constant and is measured without error. An error in airspeed will result in a difference in the amplitude of only the longitudinal component before and after the turn, while an error in the sideslip angle will similarly affect only the lateral component, which simplifies correction procedures. The airspeed maneuver modulates the attack and pitch angles; if pitch angle is measured accurately, the error in attack angle a can be determined by comparing the vertical wind component with respect to the airplane (***U*sin ) with the vertical wind component with respect to the earth (_{a}*w*). In this example, there is little correlation between the two, so the fluctuations in*w*are presumed to be due to turbulence rather than an inaccurate measurement of . The airspeed maneuver can also be used to estimate airspeed-dependent errors in other variables and the temperature recovery factor.**Fig. 4.1b. Examples of pitch and yaw maneuvers. The pitch maneuver is used as an overall check on the accuracy of the***w*measurement; in this example, there is little modulation of*w*during the pitching maneuver, which implies that fluctuations in*w*are measured accurately. Similarly, the yaw maneuver is used as an overall check on the lateral (with respect to the aircraft) component; again there is little modulation of*u*and*v*(in geographic coordinates) by the yawing maneuver.**REVERSE HEADING MANEUVER**Fly at constant altitude and heading (usually in smooth air above the boundary layer) for several minutes. Then turn 180° by first turning 90° in one direction, then 270° in the other direction at a constant rate so the airplane flies through the same volume of air on its return track. This maneuver modulates the airspeed and sideslip angle errors, since they are measured in the airplane coordinate system. The INS errors are not modulated, however, since they are measured in an inertial frame of reference. If the wind along the flight track is assumed to stay constant during this maneuver, differences in the two wind components between the two headings can be used to independently estimate errors in both airspeed and sideslip angle.

**SPEED VARIATION MANEUVER**Fly at constant altitude and heading, and smoothly vary the aircraft speed from close to stall to close to maximum cruise speed. Since the lifting force on the aircraft is directly proportional to the attack angle and the square of the airspeed, modulating airspeed also modulates attack angle. For level flight,

in (2.11); if*w*= 0_{p}, is small,*w***=**. Thus, can be calibrated in flight by this technique, provided that is measured accurately. The attitude angle transducers, in contrast to airflow angle sensors, can be accurately calibrated in the laboratory.If airspeed is measured incorrectly, temperature may also be affected. Temperature recovery factors (B.8) can also be measured or corrected with this maneuver, since airspeed variations modulate the measured temperature because of dynamic heating effects. Any other measurements affected by either airspeed or attack angle variations are also modulated by this maneuver.

**PITCH MANEUVER**Vary the aircraft elevator angle while holding the heading constant to obtain a sinusoidal pitching motion with a period of 10 to 20 s and a maximum rate of ascent and descent of 2.5 to 4 m/s. This maneuver modulates vertical airplane velocity,

; airspeed,*w*_{p}; and, to a lesser extent, attack angle, . If any of these variables have significant errors, a periodic error in vertical air velocity should be evident. Since the terms do not have the same phase angle, in practice it is often possible to determine which of the variables is in error simply by determining the phase of the error in*U*_{a}and comparing it with the phase of*w*, ,*w*_{p}and . This maneuver also can be used to detect dynamic errors in static pressure or rate-of-climb instruments by comparing their outputs with the integrated lNS vertical acceleration.*U*_{a}**YAW MANEUVER**Vary the aircraft rudder angle while holding the roll and altitude constant to obtain a sinusoidal skidding or sideslip motion with a period of about 10 s and a maximum amplitude of about 2° sideslip. This maneuver modulates heading , horizontal airplane velocity components

and*u*_{p}, and . As with the pitch maneuver, errors in any of these variables cause a periodic variation in the horizontal wind velocity.*v*_{p}On the NCAR aircraft, the system performance is judged to be satisfactory if the vertical air velocity error is less than 10% of the vertical airplane velocity for the pitch maneuver, and the lateral air velocity component is less than 10% of

sin for the yaw maneuver.*U*_{a}We estimate that short-term (i.e., not including long-term INS drift) velocity errors can be reduced to < 0.3 m/s by in-flight calibrations, which agrees well with the analysis by Nicholls (1983).

__Use of Aircraft Response Characteristics__An alternative to measuring all the components of

and the three attitude angle in order to obtain*V*_{a}is to measure the response of the aircraft to fluctuations in the air velocity and use the aircraft transfer function (i.e., the ratio of the airplane response to the input forcing function) to help estimate the air velocity. This process, however, has inherent limitations. First, the airplane transfer function is not known exactly; furthermore, it depends upon the airplane's mass, which is not constant with time because of fuel consumption. Second, the aircraft cannot respond to high-frequency variations in the air velocity. Third, the airplane responds indistinguishably to both horizontal and vertical gusts. Nevertheless, for some applications (e.g., Lenschow, 1976), useful results can be obtained by this technique. This approach was used by Bunker (1955) for turbulence flux measurements, and more recently by Kyle et al. (1976) to estimate updraft profiles in thunderstorms.*V*An intermediate approach, which is considerably more accurate, is to measure the airspeed and pitch angle of the aircraft and use its aerodynamic characteristics to calculate angle of attack, from which the vertical air velocity can be obtained. This approach is discussed by Kelly and Lenschow (1978). Here again, the accuracy depends upon aircraft flight characteristics and response time.

**LIMITATIONS OF AIR MOTION MEASUREMENTS**As with all physical measuring systems, there are limitations, both in the accuracy of the measurement and in the bandwidth (range of wavelengths) over which air velocity can be measured. Furthermore, the limitations are somewhat different for different velocity components. The horizontal velocity error, for example, has two components: a cyclic component with a period of 84.4 min and a nonperiodic component that usually increases with time. A typical navigational accuracy of the INS is approximately 0.5 m/s. Because of the periodic component, the velocity error for a position error rate of 0.5 m/s is usually larger than 0.5 m/s/h (or 1 kt/h). Exact numbers cannot be specified, since sources of error are multifarious. Some of these are the following:

- The length of the initial alignment time.
- The alignment environment (movement inside the airplane can degrade the alignment by disturbing the platform).
- The flight pattern and turbulence during the flight (and perhaps even the condition of the runway).
- Incorrect biases on the gyroscopes.
- Random changes in instrument performance.

In the vertical direction, the airplane velocity error is limited by using the measured pressure altitude as a long-term reference. Since this is a relatively accurate and stable reference, the vertical velocity of the airplane can be measured to better than 0.1 m/s if horizontal pressure gradients are not unusually large. The wind component along the longitudinal axis of the airplane is obtained from the difference between the airplane true airspeed and the INS-measured longitudinal airplane velocity. Since the airplane speed is usually several times the wind speed, this involves a small difference of two large numbers (the faster the airplane flies, the larger the numbers). Because of the upstream effects of the air flowing around the airplane and the sensor-mounting structures, the Pitot-static and static pressure measurements used in the true airspeed calculation must be empirically corrected. This is done by:

- tower fly-bys,
- flying in opposite directions through a volume of air (reverse heading maneuver) with the wind speed assumed to be constant and removing any observed difference in wind speed along the flight path by adjusting the pressure calibration, and
- the trailing cone method. (The trailing cone method {Brown, 1985 } employs a 30 m long section of tubing 0.95 cm in diameter trailing behind the aircraft with an aerodynamically shaped cone attached to the end for stabilization. The static pressure is obtained from four rings of pressure ports 2.1 m ahead of the cone. This method can be used at all normal research speeds and altitudes.)
- The heading may not be sampled frequently enough during rapid turns to provide sufficient resolution for instantaneous wind messurements.
- The air-sensing probes are not as accurate and flow distortion effects may become more important at the larger airflow angles encountered during rapid turns.

Because of the small angles involved in measuring the airflow angles (a change of 0.1° at an airspeed of 100 m/s is equivalent to a change of 17 cm/s in the lateral or vertical air velocity) and the upstream effects of flow around the airplane, the absolute calibration of the airflow angle measurement is (as discussed in Section 3.1.1) dependent upon flight tests rather than laboratory calibrations. An additional problem is the dependency of the measurement on the shape of the sensor. For example, if a new port or radome is used, or if a port becomes dented or scratched, the bias of the sensor may change significantly.

The vertical wind component measurement has an additional constraint. Flight tests cannot determine a mean absolute vertical air velocity. One approach is to assume that the vertical velocity is zero over a particular segment of a flight and reference the velocity measurement to this. It should be kept in mind, however, that this is not an absolute mean vertical velocity and furthermore that the vertical air velocity averaged over a path of longer than a few kilometers may be too small to be measurable.

Some of the drift in the vertical air velocity measurement may, under certain conditions, be reduced by using the airplane itself as a probe. This amounts to calculating the attack angle from the lift characteristics of the airplane. This can only be done by flying with constant airplane power without turns and keeping track of the airplane's weight. But even this technique does not provide an absolute mean vertical air velocity.

Because of these limitations, a "quick look" at the data processed immediately after a flight may not reveal offsets in both the lateral and vertical air velocity components. The lateral velocity offset is removed by comparing wind measurements along reciprocal headings. The vertical velocity component offset can only be approximately removed by assuming that the vertical velocity component is zero along some portion of the flight track. A further problem that should be kept in mind is the degradation in the accuracy of the air velocity measurement during turns or other abrupt changes in aircraft attitude. Although the equations used to calculate the air velocity are exact in the sense that no significant terms are omitted and small angle approximations are not used, system errors are accentuated in turns for the following reasons:

Thus far we have discussed mainly factors that affect the long-wavelength measurement of air velocity. The primary limitations at short wavelengths are vibration and frequency response limitations of the pressure lines. For the Rosemount probe, which is mounted on a nose boom protruding several meters in front of the airplane, nose-boom vibration is a limitation. The location of the air-sensing probe is a compromise between measuring in a region of minimal perturbation of the airflow by the airplane and maintaining a rigid coupling between the air-sensing probe and the INS.

This discussion is summarized in Table 2. Keep in mind that these numbers represent our present capabilities. In some cases, we expect improvements or some relaxation of limits in the future. An example is the improvement in horizontal wind measurement that is possible by updating the INS with external position or velocity information (e.g., LORAN C or other radio navigation). These potential improvements await further testing and software development.

Table 2 Present Specifications of Air Velocity Component Measurement Accuracies | |||||
---|---|---|---|---|---|

Airplane | Normal True Airspeed Range (m/sec) | Upper Limit Frequency Response (Hz) | Accuracy of Wind Components | ||

(Relative) Short-Term ( u,v,w)(< 10 min) (m/s) |
(Absolute) Long-Term (horizontal, u,v)(m/s, t in hours) | ||||

Electra |
90-125 | 10 | ±0.1 | ±(1.0 + 0.5t) | |

Sabreliner |
140-220 | 10 | ±0.1 | ±(1.0 + 0.5t) | |

King Air |
70-100 | 10 | ±0.1 | ±(1.0 + 0.5t) | |

t = time in hours from start of flight. |

The level of accuracy summarized in Table 2 is useful for a wide variety
of research objectives but is not sufficient for all proposed applications.
One class of observations for which this accuracy is not sufficient is
measurement of the integrals of velocity components around closed paths.
Divergence
and vorticity
are examples of this type of measurement. Ching
(1975)
has used this integral approach to estimate the drag coefficient
(**C _{D}**) from the vorticity equation, with the result

where

We note that by Stokes' theorem, . Thus all terms in (5.1) involve integrals of velocity components (or products of velocities) around a closed path.is the height of the boundary layer,h

is the Coriolis parameter.f

The errors in divergence or vorticity are due either to errors in the air velocity measured in the airplane coordinate system or in the velocity and attitude angles measured in an earth-based coordinate system with an INS. Constant velocity and heading errors normal to the flight path cancel out for vorticity measurements, and constant velocity errors parallel to the flight path cancel out for divergence measurements. Similarly, constant INS velocity errors integrated around a closed path cancel out for divergence and vorticity measurements. However, the denominator of (5.1) also depends upon the mean wind speed, which does not cancel out for a closed-path integration of a constant error in INS velocity.

Ching
(1975)
has estimated the terms in (5.1) from Barbados Oceanographic and
Meteorological Experiment (BOMEX) ship and rawinsonde data for a
station separation,
,
of about 5 × 10^{5} m. Typical values that he found for
these terms (excluding the time derivative term) in this maritime tropical
environment are:

If the airplane is assume to fly around a closed path in the shape of a
regular polygon, with ** n** sides of length
,
the enclosed area is
.
The ratio of the circumference,

The error in divergence or vorticity due to an error in the velocity component normal or parallel to the flight path, , integrated around the closed flight path, is

If the error is assumed to be a constant plus a linear function of time,

where ** U_{a}**, is the airplane velocity.

If the airplane is flown around the same closed path successively in opposite directions, the contribution of the mean velocity error to the divergence of vorticity should vanish. However, the term in the denominator of (5.1) depends upon the mean wind speed, and thus the mean velocity error contribution does not vanish. Therefore, for divergence and vorticity,

and for the denominator of (5.1),

For the BOMEX case, using the error given in the first paragraph,
** n** = 4 and

In order to determine the effect of long tube lengths on the time response
of pressure measurements, flight tests were conducted on a Queen Air flying
both the Rosemount probe and, for reference, a fixed-vane probe which has a
much faster response time (on the order of 0.01 s) for both Pitot pressure
and attack angle. None of the variables was altered. The power spectral
densities of the attack angle pressure difference, using tube lengths
identical to the Sabreliner's (about 5 m), are plotted in Fig. Al. There
is little obvious indication of attenuation at high frequencies or of a
large resonance peak. Because of this insensitivity, we will not consider
attenuation characteristics, but only phase angle differences. Figure A2
shows that the phase angle of the Rosemount Pitot pressure is
-6** f**°,
where

**Figure Al. Spectral density of the Rosemount differential pressure probe
attack angle output (uncalibrated units). Line of -5/3 slope is drawn for
reference. The tube length was about 5 m, and ambient static pressure was
about 790 mb.**

**Fig. A2. Phase angle of fixed-vane-probe Pitot pressure versus the
Rosemount differential pressure probe Pitot pressure. Negative phase angle
means that the Rosemount lags the fixed vane Pitot. Dashed line is a phase
lag of
= -6f°.**

**Fig. A3. Phase angle of attack angle differential pressure versus fixed
vane. Positive phase angle means that the differential pressure probe lags
the fixed vane. Dashed line is a phase lag of
= -9f°.**

We can compare these results with the theoretical results obtained by Iberall (1950) for the attenuation and phase lag of a pressure measurement at one end of a tube with an oscillating pressure applied at the other end. The theory takes into account compressibility and the volume of air contained in the instrument sensing chamber. However, it does not apply exactly to the Sabreliner, system since the tubes are terminated at ports which have a smaller diameter than the tubes. We can make one simplification; the instrument volume is negligible compared to the tube volume.

For the Pitot pressure tube with an assumed diameter of 0.36 cm and
a ratio of instrument volume to tube volume of
** V_{i}**/

For the attack angle differential pressure, if we use a tube diameter of 0.14 cm, the theoretical phase angle is about 90° at 10 Hz which agrees well with the observed phase lag.

Obviously, the effect of a varying diameter in the pressure lines is
significant. Since Iberall
(1950)
does not cover such complex geometry, we can only use his results for
qualitative comparisons with the observed results. However, we can apply
the theory quantitatively to estimate what happens when we vary the geometry
or the fluid properties. In particular, the Sabreliner typically flies
much higher than the altitude of these test results. The time constant of
the pressure lines is inversely proportional to the static pressure. At
12,200 m (40,000 ft), the atmospheric pressure is 188 mb. We add to
this 80 mb for the Sabreliner Pitot pressure. The ratio of the flight
test pressure to this pressure is 3.06. Thus, at this level the phase lags
for Sabreliner measurements of true airspeed and incidence flow angles are
-18** f**° and -27

D. H. Lenschow, P. Spyers-Duran, and C. Friehe

The true airspeed of the aircraft, ** U_{a}**, is
calculated from measurements of the Pitot-static pressure difference,

For an adiabatic process,

where
** C_{p}**/

or equivalently, using ** T_{s}**,

where

and ** R** =

Table B1 Correction factors ( COR) for specific heat of air at constant
pressure , and ratios of moist air
to dry air Mach numbers (C_{p}) and true airspeeds
(M) for several values of dewpoint temperature
(U_{a}). The values of T_{d} are
saturated mixing ratio at a static pressure of 1,000 mb.r |
||||||
---|---|---|---|---|---|---|

T_{d}(°C) |
r(× l0 ^{3}) |
e(mb) |
COR |
M (moist air)M (dry air) |
U_{a} (moist air)U_{a}
(dry air) |
% change in U_{a} |

0 | 3.84 | 6.11 | 1.0084 | 1.00071 | 1.0017 | 0.17 |

10 | 7.76 | 12.27 | 1.0160 | 1.00144 | 1.0033 | 0.33 |

20 | 14.95 | 23.27 | 1.0300 | 1.00274 | 1.0067 | 0.67 |

30 | 27.69 | 42.43 | 1.0542 | 1.00476 | 1.012 | 1.2 |

In practice, neither the total nor the ambient temperature is measured.
Instead, the measured temperature, ** T_{r}** lies between

is determined for each thermometer on the airplane. This factor is assumed to be constant and independent of humidity. Substituting (B.8) into (B.2),

and finally, substituting (B.7) and (B.9) into (B.6),

The values of
and ** C_{p}** are not completely independent of humidity.
The water vapor dependence of

where ** r** is water vapor mixing ratio (g H

The effect of water vapor on is given by (List, 1971)

where ** e** is the partial pressure of water vapor and

Andersen, N.Y., 1955: An Estimate of the Probable Behavior of
the Cornell Aeronautical Laboratory Yaw-Pitch Pitot-Static Probe in the
Mach Number Range of 0.3 to 2.0. *Report No. IH-933-P-1*. Cornell
Aeronautical Laboratory, Inc., 33 pp.

Axford, D.N., 1968: On the accuracy of wind measurements using an
inertial platform in an aircraft, and an example of a measurement of the
vertical mesostructure of the atmosphere. *J. Appl. Meteorol,*
**7**, 645-666.

Blanchard, R.L., 1971: A new algorithm for computing inertial
altitude and vertical velocity. *IEEE Trans. Aerosp., Electron.
Syst.*, **AES-7**, 1143-1146.

Britting, K.R., 1971: *Inertial Navigation Systems Analysis*.
John Wiley & Sons, Inc., New York, NY, 249 pp.

Brown, E.N., C.A. Friehe, and D.H. Lenschow, 1983: The use of pressure
fluctuations on the nose of an aircraft for measuring air motion.
*J. Clim. Appl. Meteor.*, **22**, 171-180.

Brown, E.N., 1985: Calibration of the NCAR Sabreliner Research Static
Pressure Source with a Trailing Cone Assembly. *Technical Note
NCAR/TN-253+EDD*. NCAR, Boulder, Colo., 16 pp.

Broxmeyer, C., 1964: *Inertial Navigation Systems*. McGraw-Hill Book
Co., New York, N.Y., 254 pp.

Bunker, A.F., 1955: Turbulence and shearing stresses measured over
the North Atlantic by an airplane acceleration technique. *J.
Meteorol.*, **12**, 445-455.

Ching, J.K.S., 1975: Determining the drag coefficient from
vorticity, momentum, and mass budget analysis. *J. Atmos. Sci.*,
**32**, 1898-1908.

Etkin, B., 1959: *Dynamics of Flight*. John Wiley & Sons, Inc.,
New York, N.Y., 519 pp.

Gultepe, I., A.J. Heymsfield and D.H. Lenschow, 1989: A comparison
of vertical velocity in cirrus obtained from aircraft and lidar divergence
measurements during FIRE. *J. Atmos. and Oceanic Tech.*, in press.

Johnson, H.P., D.H. Lenschow, and K. Danninger, 1978: A new fixed
vane for air motion sensing. In *Preprints, Fourth Symposium on
Meteorological Observations and Instrumentation*. 10-14 April, Denver,
Colo., American Meteorological Society, Boston, Mass., 467-470.

Kayton, M., and W.R. Fried, 1969: *Avionics Navigation Systems*.
John Wiley & Sons, Inc., New York, N.Y., 666 pp.

Kelly, T.J., and D.H. Lenschow, 1978: Thunderstorm updraft measurements
from aircraft. In *Preprints, Fourth Symposium on Meteorological
Observations and Instrumentation*, 10-14 April, Denver, Colo., American
Meteorological Society, Boston, Mass., 474-478.

Lenschow, D.H., 1971: Vanes for sensing incidence angles of the air
from aircraft. *J. Appl. Meteorol.*, **10**, 1339-1343.

Lenschow, D.H., 1972: The Measurement of Air Velocity and Temperature
Using the NCAR Buffalo Aircraft Measuring System. *Technical Note
TN/STR-74*, NCAR, Boulder, Colo., 39 pp.

Lenschow, D.H., 1976: Estimating updraft velocity from an airplane
response. *Mon. Weather Rev.*, **104**, 618-627.

Lenschow, D.H., C.A. Cullian, R.B. Friesen, and E.N. Brown, 1978a:
Status of air motion measurements on NCAR aircraft. In *Preprints,
Fourth Symposium on Meteorological Observations and Instrumentation*,
10-14 April, Denver, Colo., American Meteorological Society, Boston, Mass.,
433-438.

List, R.L., 1971: *Smithsonian Meteorological Tables*. Smithsonian
Institution, Washington, D.C., 527 pp.

Nicholls, S., 1983: *An Observational Study of the Mid-latitude,
Marine Atmospheric Boundary Layer*. Ph.D. dissertation, University of
Southampton, U.K., 307 pp.

Raymond, D.J., and M.H. Wilkening, 1982: Flow and mixing in New
Mexico mountain cumuli. *J. Atmos. Sci.*, **38**, 2211-2228.

Wyngaard, J.C., L. Rockwell and C.A. Friehe, 1985: Errors in the
measurement of turbulence upstream of an axisymmetric body. *J. Atmos.
Ocean. Tech.*, **2**, 605-614.

For Further Reading:

Lenschow, D.H., Ed., 1986: *Probing the Atmospheric Boundary Layer:
Aircraft Measurements in the Boundary Layer*. American Meteorological
Society, Boston, Mass., 39-55.

Last update: Mon Nov 12 18:35:01 MST 2001