Source: Xiaofeng Liu, Jose Roberto Moreto, and Jaime Dorado, Department of Aerospace Engineering, San Diego State University, San Diego, California
A boundary layer is a thin flow region immediately adjacent to the surface of a solid body immersed in flow field. In this region, viscous effects, such as the viscous shear stress, dominate, and the flow is retarded due to the influence of friction between the fluid and the solid surface. Outside of the boundary layer, the flow is inviscid, i.e., there is no dissipative effects due to friction, thermal conduction or mass diffusion.
The boundary layer concept was introduced by Ludwig Prandtl in 1904, which enables significant simplification to the Navier-Stokes (NS) equation for the treatment of flow over a solid body. Inside the boundary layer, the NS equation is reduced to the boundary layer equation, while outside of the boundary layer, the flow can be described by the Euler equation, which is a simplified version of the NS equation.
Figure 1. Boundary layer development over a flat plate.
The simplest case for boundary layer development occurs on a flat plate at zero angle of incidence. When considering boundary layer development on a flat plate, the velocity outside of the boundary layer is constant so that the pressure gradient along the wall is considered to be zero.
The boundary layer, which naturally develops on a solid body surface, typically undergoes the following stages: first, the laminar boundary layer state; second, the transition state, and third, the turbulent boundary layer state. Each state has its own law(s) describing the flow structure of the boundary layer.
The research into the development and structure of boundary layer is of great importance for both theoretical study and practical applications. For example, boundary layer theory is the basis for calculating the skin friction drag on ships, aircraft and the blades of turbomachines. Skin friction drag is created on the body surface within the boundary layer and is due to the viscous shear stress exerted on the surface through fluid particles in direct contact with it. The skin friction is proportional to the fluid viscosity and the local velocity gradient on the surface in the surface normal direction. Skin friction drag is present on the entire surface, thus it becomes significant over large areas, such as an airplane wing. Additionally, turbulent fluid flow creates more skin friction drag. The macro-turbulent fluid motion enhances the momentum transfer within the boundary layer by bringing fluid particles with high momentum down to the surface.
This demonstration focuses on the turbulent boundary layer over a flat plate, in which the flow is irregular, such as in mixing or eddying, and the fluctuations are superimposed on the mean flow. Thus, the velocity at any point in a turbulent boundary layer is a function of time. In this demonstration, constant temperature hot-wire anemometry, or CTA, will be used to conduct a boundary layer survey. Then, the Clauser chart method will be used to calculate the skin friction coefficient in a turbulent boundary layer.
A turbulent flow is one in which irregular fluctuations, such as mixing or eddying motions, are superimposed on the mean flow. The velocity at any point in a turbulent boundary layer is a function of time. The fluctuations can occur in any direction of the flow field, and they affect macroscopic lumps of fluid. Thus, whereas momentum transport occurs on a microscopic (or molecular) scale in a laminar boundary layer, it occurs on a macroscopic scale in a turbulent boundary layer. The size of these macroscopic lumps determines the scale of the turbulence.
The effects caused by the fluctuations are as if the viscosity were increased. As a result, the shear forces at the wall and the skin-friction component of the drag are much larger when the boundary layer is turbulent. However, since a turbulent boundary layer can negotiate an adverse pressure gradient for a longer distance, boundary layer separation may be delayed or even avoided altogether.
When describing a turbulent flow, it is convenient to express the local velocity components as the sum of a mean motion plus a fluctuating motion:
where is the time-averaged value of the u component of velocity, and is the velocity of the fluctuation. The time-averaged value at a given point in space is calculated as:
The integration interval, Δt, should be much larger than any significant period of the fluctuation velocity, , so as to converge to a mean velocity value. Thus, by definition, the converged mean value is independent of time, i.e., .
For a boundary layer on a flat plate, the external velocity is a constant. Therefore the pressure gradient term is zero. Even with this simplification, there is no exact solution for a turbulent boundary layer. However, through extensive experimental and analytical investigations on the boundary layer, the flow structure and the empirically-determined relationships describing the profile of the tangential component of the mean velocity have been established.
Very near the wall, the viscous shear dominates. To a first order, the velocity profile is linear; that is, is proportional to y. Thus, the wall shear stress can be expressed as:
where is called the skin-friction velocity and is defined as:
where τw is the skin friction, i.e., the wall shear stress. The skin friction is usually expressed in terms of the skin friction coefficient, Cf, which is defined as:
With these definitions, it is clear that for the laminar sublayer, the following relation is valid:
In the laminar sublayer, the velocity is so small that viscous forces dominate and there is no turbulence. The edge of the laminar sublayer corresponds to a y+ of 5 to 10.
In 1933, Prandtl deduced that the mean velocity in the inner region of the boundary layer must depend on the wall shear stress, i.e., skin friction, the fluid physical properties, and the distance, y, from the wall. The velocity in the inner region is thus described by the log law-of-the-wall:
In 1930, Von Kármán deduced that in the outer region of the turbulent boundary layer, the mean velocity, , is reduced below the free-stream value, , in a manner that is independent of the viscosity but is dependent on the wall shear stress and the distance, y, over which its effect has diffused. The velocity in the outer region is given by:
which is known as Law of the Wake. In this equation, is the boundary layer thickness, and is the skin-friction velocity, which is defined as:
For incompressible flow past a flat plate, the constants are defined as follows:
A suitable technique to measure turbulent boundary layer properties is by hot-wire anemometry, which is based on two principles related to the cooling effect of flow on a heated wire. The first principle is based on the heat transfer of a flow over a surface. When a fluid flows over a hot surface, the convective heat coefficient changes, which in turn affects the heat exchange rate on that surface, and consequently, may further affect the surface temperature.
The second principle is Joule's law, which states that the heat dissipation from an electrical conductor is proportional to the electric potential that is applied to the electrical conductor, squared as shown in the following equation:
where is the heat flux, I is the electric current through a conductor, R is the electrical resistance of the conductor, and U is the electric potential. One may use these two principles to correlate the velocity of fluid flow surrounding a heated metallic wire probe by measuring the electric potential applied to the probe terminals. The applied potential can be used to maintain a constant current through the wire, which is constant current anemometry or CCA, or a constant temperature on the wire, which is constant temperature anemometry or CTA.
In this demonstration, we use constant temperature anemometry (CTA) to conduct a turbulent boundary layer survey. CTA is a widely-used conventional flow diagnostic technique that has a high frequency response and can measure the small scales of turbulence without large interferences. The CTA technique utilizes a very thin metallic wire (≈ 5 μm, usually made of platinum or tungsten), which is connected to an arm of a Wheatstone bridge (Figure 2). The wire is heated to a constant temperature by applying an electrical current. Any cooling is caused by fluid flow around the wire. The Wheatstone bridge controls the electric potential applied to the wire in response to the flow velocity changes so that the heated wire resistance, and therefore the wire temperature, is held constant. The electric potential change of the Wheatstone bridge defines the signal output of the CTA.
Thus, the change in the bridge potential is a function of the heat transfer coefficient, where the heat transfer coefficient is a function of the velocity. We can obtain an empirical correlation between the airspeed and the electric potential of the bridge by calibrating the hot-wire apparatus experimentally. This involves fitting the experimental data using known heat transfer relationships.
Figure 2. TSI constant temperature anemometer model 1750. (a) Anemometer and cable connectors. (b) Electric circuit diagram, in which Rs represents the hot-wire probe.
Once air speed is calculated using CTA, we can deduce the skin friction coefficient, Cf, on the flat plate. Unfortunately, direct measurement of the skin friction drag is not available, therefore, indirect methods are used to determine its value. The Clauser Chart method is one such method. In the Clauser chart method, the measured value of the skin-friction coefficient, Cf, is determined by comparing the measured boundary layer velocity profile with a family of curves derived from the log law-of-the-wall with prescribed skin-friction coefficient values. The curve that best overlaps with the log law portion of the measured velocity profile on the semi-log plots gives the value of the measured skin friction coefficient.
1. Hot-wire system dynamic response determination
The purpose of this procedure is to understand how fast the anemometer system can respond to flow signal changes. This capability is gauged by measuring the frequency response when the signal turns on and off by applying a square wave.
2. Hot-wire calibration
The purpose of this procedure is to establish the correlation between the airspeed and the electric potential of the Wheatstone bridge. This allows the flow velocity to be measured.
Figure 3. Schematic for width of signal overshoot, τ as observed on an oscilloscope during a square wave test.
3. Boundary layer survey
A boundary layer is a thin flow region immediately adjacent to the surface of a solid body in a flow field. The region of flow outside of the boundary layer, called the free stream region has a constant velocity. However, within the boundary layer there is a velocity gradient due to friction at the surface. The boundary layer typically undergoes several stages.
First the laminar boundary state, followed by the transition state and finally, the turbulent boundary layer state, which involves irregular flow and fluctuations, like mixing or eddying. The boundary layer is the basis for the calculation of skin friction drag on aircraft.
Skin friction drag is created within the boundary layer and is due to the viscous shear stress exerted on the surface. Skin friction drag is proportional to fluid dynamic viscosity, mu, and the local stream wise velocity shear strain rate, which is the gradient of the streamwise velocity in the normal direction. So it becomes significant for large areas, such as an airplane wing. Additionally, skin friction drag is higher in turbulent flow, since the fluid particles interact with the surface at high momentum.
One way to measure turbulent boundary layer properties is using hot wire anemometry, which is based on two principles related to the cooling effect of flow on a heated wire. According to the first principle, when a fluid flows over a hot surface, the convective heat coefficient changes, which results in changes in the surface temperature.
The second principle is Joule's law, which states that an electrical conductors heat dissipation, Q, is proportional to the square of the electric current, I, applied to the conductor. We can use the two principles to determine the velocity of fluid flow surrounding a heated metallic wire probe, by measuring the electrical potential E, that has to be applied to maintain a constant temperature of the wire.
A commonly used hot wire technique is Constant Temperature Anemometry or CTA. CTA consists of a very thin metallic wire, called the probe, which is connected to the arm of a Wheatstone bridge. The Wheatstone bridge controls the electrical potential and adjusts it as needed in order to maintain a constant temperature across the wire. Any cooling is caused by fluid flow around the wire. Thus, the change in the potential is a function of the heat transfer coefficient and by extension is a function of velocity.
In this experiment, we will demonstrate the use of a Constant Temperature Anemometry setup to measure the turbulent boundary layer over a flat plate.
First, we will learn how the Constant Temperature Anemometer, or CTA, system responds to flow signal changes using a wind tunnel. To begin, secure the hot wire probe of the CTA system inside of the wind tunnel using a support shaft.
Then, set up a DC power supply, signal generator, and oscilloscope. The components are connected as shown. To begin, turn on the hot wire power supply, the signal generator and the oscilloscope. Set the signal generator to supply a square wave input to the Wheatstone bridge with a 150 mV amplitude and a 10 kHz frequency.
Observe the output signal in the oscilloscope to make sure that the frequency and amplitude are correct. Now close the test section, plug in the serial cable, turn on the wind tunnel and set the wind speed to 40 mph. Once the airflow stabilizes, measure the width of the signal overshoot, tau, observed on the oscilloscope. Use the measured value of tau to calculate the cut-off frequency for the hot wire system using this equation. Then turn off the wind tunnel airflow.
Next we will establish the correlation between wind speed and the electrical potential of the Wheatstone bridge. To begin, raise the CTA probe in the vertical direction so that it is in the free stream region. Start the wind tunnel control software and then start the virtual instrument software. Set the sampling rate to 10 kHz and the number of samples to 100,000.
Now, with the wind tunnel airspeed set to 0 mph, record the voltage on the Wheatstone bridge. Then, increase the wind tunnel airspeed at increments of 3 mph up to 15 mph, measuring the voltage at each increment. Be sure to allow the air flow to stabilize before recording the voltage measurement.
Next, increase the wind tunnel air speed up to 60 mph in 5-mph increments, measuring the voltage at each increment. When all measurements are complete, reduce the airspeed to 30 mph and then turn off the wind tunnel airflow.
Using the same setup as before, lower the CTA probe slowly until it touches the test section floor, which will act as the flat plate. Set the airflow to 40 mph. Keep the sampling frequency at 10 kilohertz and the number of samples at 100,000. Record the voltage at the lowest vertical setting, which is next to the flat plate and in the boundary layer.
Now, move the probe vertically in steps of 0.05 mm up to a height of 0. 5 mm, recording the voltage reading at each position. Then, increase the probe height in increments of 0.1 mm up to a height of 1. 5 mm. And then in steps of 0.25 mm up to a final height of 4 mm, while recording the voltage at each increment.
When all of the measurements have been made, reduce the wind speed to 20 mph and then turn off the airflow. Then shut off the power supply, signal generator, and oscilloscope.
The first step in analyzing the data is to use the data acquired during the calibration step of the experiment, to determine the correlation between the hot wire voltage and air speed. There are several different methods to do this, which involve fitting the data to known heat transfer relationships, and it's covered in detail in the appendix of this video.
Once the mathematical relationship has been determined, use the voltage measurements to calculate velocity at each vertical height. After adjusting the nominal height to account for any artifacts from an overbent probe, plot the velocity profile u(y), which can then be used to determine the boundary layer displacement thickness.
This value represents the distance that the plate would have to be moved vertically in order to obtain the same flow rate that occurs between the surface and the fluid. We can also calculate the momentum thickness, defined as shown, which is the distance the plate would have to be moved vertically in order to have the same momentum that exists between the fluid and plate.
From these two parameters, we can calculate the shape factor, H. The shape factor is used to determine the nature of the flow, where a shape factor of about 1.3 indicates fully turbulent flow and about 2.6 for laminar flow. Between these values is transitional flow. In the case of this experiment, the shape factor was calculated as 1.9, indicating transitional flow.
In summary, we learned about boundary layer flow development, and then used a Constant Temperature Anemometry setup to analyze the turbulent boundary layer over a flat plate and observe low behavior.
The CTA was calibrated in Section 2 of the protocol by measuring the voltage of the hot wire at different air speeds. This data was then used to determine the mathematical relationship between the measured variable, voltage, and the indirect variable, air speed. There are many approaches to fitting the experimental data to mathematical relationships for velocity, several of which are covered in the appendix. After the mathematical relationship is determined, velocity is easily calculated from the voltage in further experiments with the CTA.
In section 3 of the protocol, the air speed was measured using the CTA at different vertical positions in the wind tunnel. This represented different distances, y, from the flat plate. From the measured instantaneous flow velocity at each point, the average boundary layer velocity profile can be obtained. The velocity profile, u(y), can be used to determine the vertical distance that the plate would have to be moved perpendicular to itself for an inviscid flow to obtain the same flow rate that occurs between the surface and fluid, called the boundary layer displacement thickness, *. This is defined as:
where is the free stream velocity. The momentum thickness, θ, or the distance the plate would have to be moved in the direction parallel to itself in order to have the same momentum that exists in between the fluid and itself, is defined as:
Then, the shape factor, H, which can be used to determine the nature of the flow, is defined as:
where a shape factor of 1.3 indicates fully turbulent flow, a shape factor of 2.6 indicates laminar flow, and any value in between represents transition or turbulent yet not fully developed flow.
For the turbulent boundary layer case, several properties can be further examined. The skin friction can be determined using the Clauser chart method (see Figure 4). The Clauser chart method can be used to obtain the skin friction coefficient, Cf, from the measured velocity, u(y). From log law-of-the-wall, we have the following:
where κ ≈ 0.40 ~ 0.41 and B=5.0 to 5.5. Practically, κ=0.4 and B=5.5. From the definition, the skin friction coefficient is given by:
where q is the dynamic pressure of the free stream and τw is the shear stress at the wall. The log law-of-the-wall can then be expressed as (See Appendix):
where, .
Given a series of Cf values, a family of curves can be generated for vs. Ry. Several values of Ry ranging from 100 to 100,000 and Cf values ranging from 0.001 to 0.006 should be used to plot the curves in a log-linear format. This forms the Clauser chart, which can be used to determine the skin friction coefficient, Cf, as shown in Figure 4. By comparing the measured boundary layer velocity profile with the family of curves that are based on the log law-of-the-wall with the prescribed skin-friction coefficient values, the curve that best overlaps with the log law portion of the measured velocity profile gives the value of the measured skin friction coefficient.
Figure 4: Clauser Chart.
This result can be compared to the result obtained using the integral equation method. Also, the velocity fluctuation profile can be obtained and the experimental result can be compared against the log law-of-the-wall. See the Appendix for more information.
The demonstration shows how to use constant temperature anemometry, a powerful tool used to study turbulent flow over a surface, which in this specific case was a flat plate. This method is simpler and less expensive than other methods, such as PIV, PTV, and LDV, and it provides a high temporal resolution. The application of hot-wire anemometry to a turbulent boundary layer provides a cost effective and hands-on approach to demonstrate the behavior of turbulent flows.
Constant temperature anemometry has numerous applications. This technique can be used to survey both turbulent and laminar flows. Hot-wire anemometry can be used to study the wake flows of an airfoil or an airplane model, thus providing information such as the drag of the airfoil and the level of wake turbulence, which provides valuable information for aircraft design.
Hot-wire anemometry can also be used in environmental fluid dynamics investigations, such as to study plume flows, which are responsible for the mass and momentum transport and mixing of a variety of processes found in the Earth`s atmosphere.
A variation to hot-wire anemometry is hot-film anemometry, which is typically used in liquid flows that require robust and reliable performance. For example, monitoring of the air flow at the air intake duct of an automobile engine is often performed by a sensor made of hot film.
The application of hotwire anemometry is not restricted to the mechanical engineering realm. CTA can also be used for example in biomedical applications to measure respiration rate.
Materials List
Name | Company | Catalog Number | コメント |
Equipment | |||
Instructional Subsonic Wind Tunnel | Jetstream | The dimensions of the test section of the wind tunnel are as follows: 5.25" (width) x 5.25"(height) x 16" (length). The wind tunnel should be able to attain air speeds of 0 – 80 mph. | |
The Wall | The wall of the test section is made of glass. | ||
CTA model 1750 | TSI Corp. | ||
Hot-wire probe | TSI Corp | TSI 1218-T1.5 | Tungsten-platinum coated, standard boundary layer probe. The diameter of the probe is 3.81 μm. The length of the sensing area of the wire is 1.27 mm. |
A/D Board | National Instruments | NI USB 6003 | Maximum sampling rate of 100 kHz with 16-bit resolution |
Traverse System | Newport | Newport 370-RC Rack-And-Pinion Rod Clamp & 75 Damped Optical Support Rod Assembly | |
Pitot tube | The dynamic pressure of the free stream will be sensed by a tiny Pitot tube installed at the beginning region of the test section. The resolution of the Pitot tube is 0.1 mph. | ||
Software | LabView software will be used for data acquisition. | ||
Power Supply | Heath | 2718 | Heath 2718 Tri-Power Supply with 15V DC output is used to power the hot-wire anemometer. |
Oscilloscope | Tektronix | 2232 | |
Signal Generator | Agilent | 33110A |
A boundary layer is a thin flow region immediately adjacent to the surface of a solid body in a flow field. The region of flow outside of the boundary layer, called the free stream region has a constant velocity. However, within the boundary layer there is a velocity gradient due to friction at the surface. The boundary layer typically undergoes several stages.
First the laminar boundary state, followed by the transition state and finally, the turbulent boundary layer state, which involves irregular flow and fluctuations, like mixing or eddying. The boundary layer is the basis for the calculation of skin friction drag on aircraft.
Skin friction drag is created within the boundary layer and is due to the viscous shear stress exerted on the surface. Skin friction drag is proportional to fluid dynamic viscosity, mu, and the local stream wise velocity shear strain rate, which is the gradient of the streamwise velocity in the normal direction. So it becomes significant for large areas, such as an airplane wing. Additionally, skin friction drag is higher in turbulent flow, since the fluid particles interact with the surface at high momentum.
One way to measure turbulent boundary layer properties is using hot wire anemometry, which is based on two principles related to the cooling effect of flow on a heated wire. According to the first principle, when a fluid flows over a hot surface, the convective heat coefficient changes, which results in changes in the surface temperature.
The second principle is Joule’s law, which states that an electrical conductors heat dissipation, Q, is proportional to the square of the electric current, I, applied to the conductor. We can use the two principles to determine the velocity of fluid flow surrounding a heated metallic wire probe, by measuring the electrical potential E, that has to be applied to maintain a constant temperature of the wire.
A commonly used hot wire technique is Constant Temperature Anemometry or CTA. CTA consists of a very thin metallic wire, called the probe, which is connected to the arm of a Wheatstone bridge. The Wheatstone bridge controls the electrical potential and adjusts it as needed in order to maintain a constant temperature across the wire. Any cooling is caused by fluid flow around the wire. Thus, the change in the potential is a function of the heat transfer coefficient and by extension is a function of velocity.
In this experiment, we will demonstrate the use of a Constant Temperature Anemometry setup to measure the turbulent boundary layer over a flat plate.
First, we will learn how the Constant Temperature Anemometer, or CTA, system responds to flow signal changes using a wind tunnel. To begin, secure the hot wire probe of the CTA system inside of the wind tunnel using a support shaft.
Then, set up a DC power supply, signal generator, and oscilloscope. The components are connected as shown. To begin, turn on the hot wire power supply, the signal generator and the oscilloscope. Set the signal generator to supply a square wave input to the Wheatstone bridge with a 150 mV amplitude and a 10 kHz frequency.
Observe the output signal in the oscilloscope to make sure that the frequency and amplitude are correct. Now close the test section, plug in the serial cable, turn on the wind tunnel and set the wind speed to 40 mph. Once the airflow stabilizes, measure the width of the signal overshoot, tau, observed on the oscilloscope. Use the measured value of tau to calculate the cut-off frequency for the hot wire system using this equation. Then turn off the wind tunnel airflow.
Next we will establish the correlation between wind speed and the electrical potential of the Wheatstone bridge. To begin, raise the CTA probe in the vertical direction so that it is in the free stream region. Start the wind tunnel control software and then start the virtual instrument software. Set the sampling rate to 10 kHz and the number of samples to 100,000.
Now, with the wind tunnel airspeed set to 0 mph, record the voltage on the Wheatstone bridge. Then, increase the wind tunnel airspeed at increments of 3 mph up to 15 mph, measuring the voltage at each increment. Be sure to allow the air flow to stabilize before recording the voltage measurement.
Next, increase the wind tunnel air speed up to 60 mph in 5-mph increments, measuring the voltage at each increment. When all measurements are complete, reduce the airspeed to 30 mph and then turn off the wind tunnel airflow.
Using the same setup as before, lower the CTA probe slowly until it touches the test section floor, which will act as the flat plate. Set the airflow to 40 mph. Keep the sampling frequency at 10 kilohertz and the number of samples at 100,000. Record the voltage at the lowest vertical setting, which is next to the flat plate and in the boundary layer.
Now, move the probe vertically in steps of 0.05 mm up to a height of 0. 5 mm, recording the voltage reading at each position. Then, increase the probe height in increments of 0.1 mm up to a height of 1. 5 mm. And then in steps of 0.25 mm up to a final height of 4 mm, while recording the voltage at each increment.
When all of the measurements have been made, reduce the wind speed to 20 mph and then turn off the airflow. Then shut off the power supply, signal generator, and oscilloscope.
The first step in analyzing the data is to use the data acquired during the calibration step of the experiment, to determine the correlation between the hot wire voltage and air speed. There are several different methods to do this, which involve fitting the data to known heat transfer relationships, and it’s covered in detail in the appendix of this video.
Once the mathematical relationship has been determined, use the voltage measurements to calculate velocity at each vertical height. After adjusting the nominal height to account for any artifacts from an overbent probe, plot the velocity profile u(y), which can then be used to determine the boundary layer displacement thickness.
This value represents the distance that the plate would have to be moved vertically in order to obtain the same flow rate that occurs between the surface and the fluid. We can also calculate the momentum thickness, defined as shown, which is the distance the plate would have to be moved vertically in order to have the same momentum that exists between the fluid and plate.
From these two parameters, we can calculate the shape factor, H. The shape factor is used to determine the nature of the flow, where a shape factor of about 1.3 indicates fully turbulent flow and about 2.6 for laminar flow. Between these values is transitional flow. In the case of this experiment, the shape factor was calculated as 1.9, indicating transitional flow.
In summary, we learned about boundary layer flow development, and then used a Constant Temperature Anemometry setup to analyze the turbulent boundary layer over a flat plate and observe low behavior.