28 Time-Dependent Flows
Introduction
In fluid dynamics and aerodynamics, a time-dependent or unsteady flow refers to the conditions where the properties (such as velocity, pressure, temperature, and density) at a particular point in space will change with time. This behavior contrasts with steady flow, where the fluid properties at a given point remain constant over time. Unsteady flows occur, and indeed are expected, in many real-world situations, such as when a fluid flow is subjected to changes in operating conditions, the initiation of transient events, or when time-varying external events produce non-steady forces on flows or flow systems.
In aerospace engineering, unsteady flow effects occur frequently. For example, they are found and must be considered to be able to predict the fluid behavior in the fields of aeroelasticity and flutter, rotating machinery, including jet engines and helicopter rotors, mixing processes, hydraulic and pneumatic control systems, and combustion processes in air-breathing and rocket engines. Analyzing and predicting unsteady flow can be more challenging than steady flow because it requires careful consideration of all relevant time-dependent factors. Closed-form solutions to most unsteady flow problems are generally only realizable in cases with specific assumptions. Therefore, experimental techniques and computational fluid dynamics (CFD) are often used to study unsteady flows in various engineering applications.
Learning Objectives
- Understand how to differentiate between steady and unsteady flows and recognize situations in which unsteady flow properties are relevant to problem-solving.
- Know about reduced frequency and Strouhal number and why they are essential for characterizing unsteady flows and the degree of unsteadiness.
- Be able to apply the conservation principles of fluid dynamics to solving some exemplar unsteady flow problems.
Classification of Time-Dependent Flows
Steady flows are flows with properties that do not change as a function of time. A steady-state flow refers to the condition where the macroscopic flow properties, such as the velocity and pressure at a point, do not change with respect to time, as shown in the figure below on the left. Mathematically, for steady flows, then
(1)
where is any property such as pressure, temperature, velocity, density, etc. However, in a time-dependent flow, also known as an unsteady flow, the flow properties at a point will change in time, as shown in the figure below on the right. In this case, all of the unsteady terms in any equations used to describe the flow must be retained and solved, i.e.,
(2)
Unsteady flow phenomena are encountered in many engineering applications. Examples include the flows in turbomachinery and other internal combustion engines, helicopter aerodynamics, wind turbines, and numerous problems in aeroacoustics where creating time-varying aerodynamic loads produces sound (noise). A turbulent flow is an unsteady flow, by definition. However, a turbulent flow can still be statistically steady. This definition means that the average flow velocity and other quantities are constant with respect to time, and all the statistically varying properties, such as the fluctuating velocity component, are constant in time.
The figure below shows the difference between a statistically steady turbulent flow and a statistically unsteady flow. Turbulence enhances the mixing of fluids, promoting the transport of momentum, heat, and other properties within the fluid. The irregular motion of fluid particles makes it difficult to predict the behavior of the flow over time precisely; this randomness, or mathematically stochastic behavior, is a fundamental characteristic of turbulence. However, statistically a flow property can be decomposed into a mean or average part, , and a mean fluctuating part, , i.e.,
(3)
This latter process is called a Reynolds decomposition, which is a helpful concept in modeling turbulent flows.
One reason it is helpful to distinguish between steady and unsteady flows is that the former is often more tractable to understand and predict. Eliminating time from the solution of the equations that govern fluid flow problems usually results in a significant simplification of the governing equations as well as the mathematical and/or numerical techniques needed to solve these equations. Characterizing steady, quasi-steady, or unsteady flows is often done using parameters such as reduced frequency or Strouhal number.
Reduced Frequency
One nondimensional parameter used to help categorize whether a flow is steady or unsteady is the reduced frequency, given the symbol . The reduced frequency is defined, in general, as
(4)
where is a characteristic physical frequency of the unsteady flow (in units of radians per second), is a characteristic length scale, and is a reference flow velocity. Base units should always be used to obtain the correct numerical value for the reduced frequency. For a wing or airfoil, such as one oscillating in angle of attack or an oscillatory vertical (heaving) motion as shown in the figure below, the reduced frequency is often defined in terms of its semi-chord, i.e., , and the free-stream velocity, , i.e., in this case, it is defined as
(5)
The aerodynamics of oscillating airfoil sections and wings must be addressed in the field of aeroelasticity and flutter.
For = 0, the flow is steady, and all the usual aerodynamic results apply regarding the relationships between the aerodynamic quantities and the angle of attack. For 0 0.05, the flow can be considered quasi-steady; that is, unsteady effects are generally minor, and for some problems, they may be neglected completely. However, such reduced frequency bounds may be somewhat subjective.
Flows with characteristic reduced frequencies of 0.05 and above are usually considered unsteady, so all unsteady terms must be retained in the governing flow equations. Such problems are more challenging to determine for two reasons. First, the local flow properties also depend on the previous time, i.e., what has happened in the time history of the lift and other aerodynamic forces. Second, the effects of compressibility may be significant even if the free-stream flow velocity is low. This means that the time scales of the unsteady motion become of the order of propagating pressure disturbances in the flow, which occur at the speed of sound.
Strouhal Number
The Strouhal number, given the symbol , is another dimensionless number used in fluid dynamics to characterize the behavior of oscillating or vibrating objects within a fluid flow. It is beneficial for analyzing phenomena such as vortex shedding and sound generation. The Strouhal number is defined as the ratio of the oscillation frequency to the product of the object’s characteristic length (or diameter) of the object and the reference or free-stream velocity of the fluid, e.g., . The equation is given by
(6)
Notice that the physical frequency (in units of per second or Hz) is used to determine the Strouhal number rather than the circular frequency in the reduced frequency. The Strouhal number is named after the Czech physicist and engineer Vincenc Strouhal, who significantly contributed to the study of oscillating flows in the late 19th and early 20th centuries.
The Strouhal number is often used to analyze various fluid flow phenomena, such as the shedding of vortices behind a cylinder or another type of bluff (non-streamlined) body. It helps characterize the vortex shedding frequency relative to the fluid flow velocity and the object’s size. For different types of flows and body shapes, there are often characteristic Strouhal number ranges corresponding to specific flow regimes and behaviors, i.e., the behavior locks into a particular shedding frequency.
Check Your Understanding #1 – Calculating the value of the Strouhal number
A circular cylinder with a diameter of 0.1 m is immersed in a fluid flow with a velocity of 1.0 m/s. The shedding frequency of vortex shedding behind the cylinder is measured to be 2 Hz. What is the Strouhal number in this case?
Show solution/hide solution.
The Strouhal number is defined by
Using the numerical values gives
This value of 0.2 is often associated with the onset of vortex shedding behind a cylindrical object in a fluid flow. For example, suppose a particular combination of flow speeds and length scales results in a Strouhal number of 0.2. In that case, vortex shedding can be expected, and vortex-induced vibrations on a structure are possible. However, it is essential to note that specific Strouhal number values may vary for different flow conditions and other body shapes.
Time-Dependent Fluid Flows
Time-dependent flows are often more challenging to deal with in fluid dynamics, but the same principles of conservation of mass, momentum, and energy will apply. Studying exemplars is an excellent way to learn how to solve unsteady-flow problems. In this regard, several classic exemplars in fluid dynamics can be used to lay down the solution principles when considering the additional dimension of time.
Torricelli’s Law
Torricelli’s law relates the velocity of a liquid flowing out of an orifice to the height of the liquid above the level of the orifice, as shown in the figure below. The derivation of this law assumes ideal conditions where the liquid is quasi-steady, incompressible, and inviscid. While the problem is fundamentally unsteady in that the fluid properties at different points in the flow change with time as the liquid drains out of the tank because their time rates of change are small, then the flow can be assumed quasi-steady. Therefore, the applicability of the Bernoulli equation can be assumed as it describes the instantaneous balance of energy between different points in the flow.
The application of Bernoulli’s equation between levels 1 and 2 gives
(7)
where is the pressure at the surface, and is the pressure at the discharge from the orifice. The fluid is a liquid, so = constant. It is assumed that there are no losses from the discharge through the orifice. Assuming (i.e., both are equal to atmospheric pressure, ), and , the Bernoulli equation simplifies to
(8)
Solving for the discharge or “jet” velocity, , gives
(9)
This result shows that the discharge velocity is proportional to the square root of the difference in hydrostatic height or “head” between the free surface at and the orifice at .
To find the horizontal distance, , reached by a jet from any one of the orifices requires the time, , it takes for the jet to reach the ground, i.e.,
(10)
where is the height of the orifice above the ground. Therefore, the horizontal distance reached by the jet, , is given by
(11)
This latter result shows that the distance traveled by the fluid jet depends on the hydrostatic head of fluid above the orifice, , and the height above the ground, .
Time for Liquid to Drain from a Tank
Now consider the problem of a tank of circular cross-section, , that lets a liquid leave through a drain valve at the bottom, as shown below. The discharge velocity, , of the flow from the drain varies with the height, , of the fluid level above the drain according to the relationship . Notice that the height, , will decrease with time as the liquid leaves the tank.
The application of the continuity equation gives
(12)
where is the area of liquid discharge from the drain. The volume of liquid in the tank for any height is
(13)
where is the cross-sectional area of the tank. It is given that
(14)
and this flow rate depends on the instantaneous height of the liquid in the tank, , which, as previously derived, is called Torricelli’s theorem. Therefore, using the conservation of mass gives
(15)
Separating the variables and integrating them gives
(16)
where the limits of integration are: At , then , and when , the tank is considered empty. Performing the integration gives
(17)
In terms of the diameter of the tank, , and the outlet diameter, , of the drain, then
(18)
Therefore, the time required to empty the tank will be
(19)
where is the initial height of the liquid in the tank when the drain is first opened.
Filling an Air Tank
Consider now a rigid tank of volume with air pumped at a constant mass flow rate, as shown in the figure below. In this problem, the effect of compressibility must be considered. Still, it can be assumed that the process is slow enough that it is isothermal, and the heat generated by compressing the air is dissipated through radiation by the tank.
This is also an unsteady flow problem because a mass of air is pumped into a fixed volume, and for mass conservation, the air density must increase in time. The initial density and pressure are and , respectively. The general form of the continuity equation is
(20)
In this case
(21)
and so
(22)
Assume a uniform mixing of air so the density inside the volume is uniform. Therefore,
(23)
Integrating using the separation of variables gives
(24)
so that
(25)
which shows that the density of the air increases linearly with time. The corresponding pressure can be obtained from the equation of state, i.e., . If the process is assumed to be isothermal, then , and so
(26)
Hydraulic Shock
Hydraulic shock, also known as “water hammer” or “hydraulic hammer,” is a time-dependent flow phenomenon that occurs in fluid flow systems, such as hydraulic and fuel systems, when there is a sudden change in the flow velocity. This abrupt change can result in strong pressure waves propagating through the fluid, causing transitory surges in pressure inside the system, as shown in the schematic below. Water hammer is common in fluid flow systems where valves and regulators are rapidly opened or closed. The rapid, time-dependent changes in pressure can generate loud banging or knocking sounds, often audible throughout the piping system.
The fundamental cause of water hammer is the fluid’s inertia and the need for energy conservation when a sudden stoppage of the flow occurs. When a fluid is in motion and experiences a sudden change in velocity, its kinetic energy must be either absorbed or dissipated. If a valve is closed suddenly, the fluid comes to an abrupt stop, all kinetic energy is converted into an increased pressure, and the flow stagnates in the pipe.
This process then creates a pressure wave that travels back upstream through the pipe; this pressure wave travels through the fluid at the speed of sound in the medial, , known as the wave speed. When the pressure wave reaches another point where the fluid is stagnant, i.e., the reservoir in this case, the higher pressure head starts the wave downstream again, with continuing oscillations until the energy of the pressure wave is dissipated through frictional effects. The increase in pressure causes a “hammer” effect, i.e., a thumping effect, which can be felt and/or heard by an observer.
The resulting pressure pulses can be several times higher than the normal operating pressure. The change in pressure from the water hammer depends on how quickly the value is closed, the type of fluid, the length of the pipe, its elastic properties, and how the pipe is mounted, e.g., supported at its ends and/or clamped along its length. Prolonged or severe water hammer can lead to structural damage to the pipes, their fittings, and other system components. The phenomenon can cause fatigue, leaks, and pipe bursts and accelerate the wear and tear on valves and pumps, increasing maintenance requirements. Preventing or mitigating water hammer effects involves engineering solutions and using various devices, such as surge tanks, air chambers, and water hammer arrestors, which are designed to absorb the excess pressure and prevent it from causing damage to the system.
Consider the analysis of this problem. A liquid stored in a tank flows steadily through a pipe of length , as shown in the figure below. At the time , the valve at the downstream end is quickly closed, producing the classic pressure pulse of a water hammer. Water hammer transients are typically axisymmetric because the axial mass, momentum, and energy changes are significantly more significant than their radial counterparts.
Nikolay Joukowsky laid the foundation of the water hammer theory[1] where the pressure amplitude, , in the liquid in the pipe is related to the change in flow velocity, , using
(27)
where is the wave propagation speed (assumed to be constant) and is the density of the fluid; this equation is basically what is referred to today as piston theory. The negative sign in Eq. 27 describes a water hammer wave moving downstream, while the positive sign describes the wave moving upstream. Because pressure head is often used in the field of hydraulics, then Eq. 27 can also be written as
(28)
There are three cases of interest when it comes to predicting water hammer effects, all of which involve the wave propagation speed:
- Gradual closure of the valve.
- Sudden closure of the valve.
- Closure of the valve when allowing for the elasticity of the pipe, which will affect the wave propagation speed.
The transit time, also known as reflection time, taken for the pressure wave to propagate from the valve to the tank and then back to the valve will be
(29)
If the time it takes to close the valve is , then if is longer than the reflection time, the valve closure is considered gradual. If is shorter than time, the reflection of the valve closure is referred to as being sudden.
If , which is a gradual closure of the valve, then the increase in pressure from the water hammer is given by
(30)
is the average flow velocity inside the pipe before the valve closes. The equivalent pressure head of the water hammer is
(31)
If , which is a sudden closure of the valve, then the increase in pressure from the water hammer is given by
(32)
and the equivalent pressure head is
(33)
where represents the average flow velocity of the fluid inside the pipe before the valve is closed. This velocity is critical in determining the magnitude of the pressure increase caused by the water hammer effect. It reflects the kinetic energy of the moving fluid that is converted into pressure when the flow is suddenly stopped.
Notice that in both cases, the propagation speed of the pressure wave, , is needed. If the pipe is rigid, then
(34)
where is called the bulk modulus of the liquid, for which values for various liquids are available. This latter equation is often referred to as the Newton-Laplace equation.
If the pipe is elastic, for which most will be to a lesser or greater degree, then
(35)
where is called the effective bulk modulus of the liquid in the pipe. The value of can be obtained using
(36)
where is the modulus of elasticity of the pipe, is the pipe diameter, and is the wall thickness of the pipe. The value of depends on exactly how the pipe is mounted and anchored; usually . The modulus of elasticity of the pipe is a material property for which values are also widely available.
Finally, consider the valve’s sudden closure when the pipe’s elasticity is considered. In this case
(37)
so the elasticity of the system reduces the pressure effects from water hammer.
Check Your Understanding #2 – Determining hydraulic shock pressure in a pipeline
A hydraulic fluid with a density is flowing through a titanium pipe in an aircraft’s hydraulic system. The pipe has a diameter of and a wall thickness of . The bulk modulus of the hydraulic fluid is , and the modulus of elasticity of titanium is . The hydraulic fluid has an average velocity before a valve is suddenly closed in the system. Assume . Calculate the pressure increase in the pipe from the hydraulic show (water hammer) effect.
Show solution/hide solution.
To calculate the pressure increase, we first need the wave speed , which depends on the effective bulk modulus , i.e.,
where is calculated using
Substituting the given values into the equation for gives
The wave speed, , can now be calculated, i.e.,
The pressure increase from the water hammer effect (sudden valve closure) is given by
Substituting the known values gives the change in pressure as
Why is the speed of sound so much higher in a liquid than in air?
The speed of sound in a medium, whether it’s a gas, liquid, or solid, depends on the properties of that medium. Sound is faster in liquids and solids than in gases because the molecules in liquids and solids are much closer together, allowing for more rapid transmission of mechanical waves. The speed of sound in a medium is determined by its density and elastic properties, specifically its bulk modulus. The bulk modulus measures a substance’s resistance to compression or expansion when subjected to pressure. Liquids and solids generally have higher bulk moduli than gases, making them less compressible. In a liquid, the molecules are closer together than in a gas, and they can still transmit mechanical waves. This closeness of molecules and the resistance to compression contribute to a higher sound speed in liquids than in gases.
Producing Jet Thrust
A monopropellant thruster is a basic rocket, often used for satellite propulsion and control, with a certain propellant mass stored in a pressurized storage tank. The propellant is released over time by opening the valve to flow over a catalyst bed, causing thermal decomposition and the liberation of energy from the propellant, as shown in the figure below. The flow then expands through a nozzle to produce an exit velocity and, hence, a thrust from the time rate of change of momentum of the expanding gases. The idea is that the thrust boosts the satellite to change its orbital velocity and altitude.
This is a time-varying flow problem because a propellant mass is continuously discharged from a tank, so the mass of the propellant inside the control volume is decreasing. The most general form of the continuity equation is
(38)
and in this case, it can be reduced to
(39)
where is the mass of the propellant and is the propellant flow rate. Notice that the minus sign occurs because the total propellant mass is decreasing.
Let the initial mass of the spacecraft be . Conservation of mass gives the current mass of the spacecraft, , at time later as the propellant is discharged as
(40)
where = constant. The thrust force from the propulsion system is given by Newton’s second law (force equals time rate of change of momentum), i.e.,
(41)
and so the acceleration of the spacecraft, which is not constant because of its continuously decreasing mass, is given by
(42)
Therefore,
(43)
Integrating the equation gives
(44)
which has a special name called the rocket equation.
Check Your Understanding #3 – Discharge of a monopropellant thruster
An orbiting satellite is propelled using a monopropellant thruster. The satellite has an initial mass of 5,000 kg, including its propellant mass. To give a slight corrective boost to its orbit, it opens the control valve and ejects the propellant at a constant rate of 0.1 kg/s at an effective exit velocity of 2,000 m/s. Assume that the satellite operates in a vacuum and no gravity forces are acting on it. Determine the thrust, the change in velocity of the satellite, and the new mass after 90 seconds from the start of the propulsive burn.
Show solution/hide solution.
The thrust produced can be determined from the time rate of change of momentum of the exhaust gases, i.e.,
After 90 seconds of thrusting, i.e., = 90, then using the rocket equation gives
The final mass after the burn will be
Sedimentation of Particles in a Fluid
The sedimentation of solid particles suspended in a fluid is a classic problem in fluid dynamics and is relevant to many natural[2] and industrial processes.[3] It involves understanding how suspended solid particles progressively settle in time onto a solid surface under gravity through a fluid, such as water or air, and the forces acting on them during this process. This situation is also a time-dependent flow problem, as illustrated in the figure below.
As the particles settle over time, they form a sediment layer on the surface. The concentration of particles in the sediment layer eventually becomes much higher than in the suspension above. A sharp interface or “sedimentation front” often develops as particles settle between the fluid above and the particle-rich suspension below. This front then moves downward over time as particles settle onto the bed. For tiny particles, Brownian motion (i.e., the random thermal motion of particles) can oppose or enhance sedimentation, depending on the particle size and the fluid viscosity. This behavior leads to diffusion effects, spreading the particle distribution over time.
Stokes’ drag law applies to small, independent, spherical particles settling slowly in a viscous fluid with negligible inertial effects, i.e., they move at low Reynolds numbers. According to Stokes’ law, the drag force, , on a particle of radius, , moving with velocity, , in a fluid with viscosity, , is given by
(45)
The settling of the particle involves a balance between the gravitational force pulling the particle down and the drag force resisting its motion. Buoyancy forces, , albeit small, also matter, which act upward, opposing the particle’s weight, . The difference between the gravitational force from the particle’s weight and the buoyancy force is the effective net force driving the sedimentation, i.e., the force equilibrium equation is
(46)
giving
(47)
where and are the densities of the particle and fluid, respectively. On rearrangement this gives the terminal velocity, , or the settling speed, i.e.,
(48)
In dilute suspensions, the particles will sediment independently, and the settling time from particles at a given height will be
(49)
Deviations from Stokes’ law occur for larger or more irregular particles, and corrections or empirical relationships are used to estimate terminal velocity. At higher concentrations, particle interactions also slow the overall settling rate. This phenomenon is known as “hindered settling.” The settling velocity, , in a suspension of many particles is lower than the terminal velocity of a single particle. An empirical relationship often used is
(50)
where is the volume fraction of particles in the suspension. The value of is an empirically determined exponent, typically around 4.65 for small particles in water. Using the hindered settling velocity, the time to settle from an initial height is
(51)
For fine particles, e.g., dust and aerosols, the exponent is typically lower than in water, in the range of 2 to 3, depending on the specific characteristics of the particles and the level of interactions. For larger particles or irregular shapes, the value of is closer to 1 to 2 because the effects of hindered settling are less pronounced in a less viscous medium.
In natural systems, particles may settle in a fluid with a non-uniform flow field, such as in a turbulent environment or where fluid recirculation occurs. This behavior introduces additional forces on the particles and can lead to more complex sedimentation patterns, e.g., the formation of humps or dunes. Furthermore, in many practical scenarios, the particles in suspension are not uniform in size, i.e., they are polydisperse. Larger particles will settle faster than smaller ones, leading to size segregation over time as they form the sediment bed.
Primer on Unsteady Airfoil Behavior
Unsteady airfoil theory deals with the aerodynamic behavior of airfoils and wings when subjected to time-varying boundary conditions, such as oscillations in angle of attack, various vertical and horizontal gusts, or other transient phenomena. Understanding unsteady aerodynamics is crucial for applications involving aircraft maneuvering, rotorcraft aerodynamics, and wind turbine blades.
In an unsteady flow, the wake behind the airfoil is time-varying in structure with the shedding of vorticity and circulation, thereby influencing the aerodynamic forces experienced by the airfoil. These effects change the magnitude of the lift and the phase of the lift response relative to the time-varying boundary conditions. At high angles of attack, the flow may separate from the airfoil surface, causing a rapid loss of lift and increased drag, a phenomenon called dynamic stall.
There are many ways of calculating these unsteady effects. One approach that is very instructive for those first learning the subject of unsteady airfoil behavior is the classic unsteady aerodynamics theories, including Theordorsen’s theory, which is discussed next. Unsteady panel methods discretize the airfoil surface into panels on the surface and extend into the wake. CFD methods, which are the most computationally intensive, solve the Navier-Stokes equations for unsteady flow, which can capture detailed flow phenomena such as vortex shedding and dynamic stall.
Quasi-Steady Aerodynamics
In quasi-steady flows, the behavior can be evaluated by applying the principles of steady flow under the instantaneous boundary conditions of flow tangency to the airfoil surface. For example, for a thin airfoil, the changes in the angle of attack of a harmonically oscillating airfoil can be expressed as
(52)
where is the mean angle of attack and is the angle of attack amplitude of the oscillation. If the flow were quasi-steady, in the sense that flow adjustments were to take place instantaneously, the lift coefficient would be given by
(53)
recognizing that the pitch rate, i.e., , must affect the lift coefficient because it changes the angle of attack. This equation then leads to
(54)
where is a frequency-dependent phase angle given by
(55)
Therefore, even based on quasi-steady arguments, the lift response will no longer be in phase with the angle of attack (it will lead by the angle in this case) nor have the same amplitude (it is larger than by a factor ). The upshot is that the aerodynamic forces will have different values depending on whether the angle of attack is increasing or decreasing with respect to time.
Indeed, measurements suggest that the unsteady lift on an airfoil or wing is attenuated (diminished) by unsteady effects not amplified as the quasi-steady theory would suggest, and there is a lag (not a lead) in the aerodynamic response with respect to the angle of attack. These effects are formally embodied within the theory of unsteady aerodynamic behavior, for which extensive literature exists. One of these theories is called Theordorsen’s theory.
Theodorsen’s Theory
Theodorsen’s theory[4] provides a linearized, small-disturbance solution for oscillating airfoils in an incompressible flow. The theory is based on the assumption of small oscillation amplitudes. It gives an exact analytic theory (within the assumptions and limitations) to describe the unsteady air loads on a thin airfoil. Theodorsen’s problem was obtaining the solution for the loading (bound circulation) on the airfoil surface under harmonic forcing conditions, the flow model being shown in the figure below.
The governing integral equation is
(56)
where is the downwash on the airfoil surface. This equation must be solved subject to invoking the Kutta condition at the trailing edge, i.e., . There is also a connection to be drawn between the change in circulation about the airfoil and the circulation shed into the wake. Conservation of circulation requires that
(57)
Assuming that the vorticity in the shed wake is convected at the free stream velocity, , this gives , and so
(58)
where the total circulation about the airfoil, is given by
(59)
This wake vorticity changes the downwash velocity over the airfoil, and the loads on the airfoil are also affected. Suppose the net circulation about the airfoil is changing with respect to time. In that case, circulation will be continuously shed into the wake and will continuously affect the aerodynamic loads on the airfoil. In the limit, as the changing boundary conditions become constant, the shed wake vorticity cast off the trailing edge of the airfoil becomes zero, and the remaining circulation in the wake convects downstream to infinity.
The unsteady problem posed above is certainly not trivial to solve. Still, for simple harmonic motion, the solution is given by Theodorsen in a form that represents a transfer function between the forcing (angle of attack) and the aerodynamic response (i.e., the chordwise loading distribution, lift, and pitching moment). The Theodorsen function, denoted by , is a complex-valued function that captures the instantaneous effect of the circulation shed into the wake on the unsteady aerodynamic response and is given by
(60)
where and are Hankel functions of the second kind, and is the reduced frequency. The Hankel function is defined as , with and being Bessel functions of the first and second kind, respectively. The function is an aerodynamic transfer function that relates the input (forcing) to the response.
Recognizing that each Bessel function has an argument , then the real or in-phase () part and imaginary or out-of-phase part () of the Theodorsen function can be written as
(61)
and
(62)
Despite its apparent mathematical complexity, the Theodorsen function is readily calculated numerically using MATLAB and is plotted in the figure below.
Another way of interpreting the Theodorsen function, which is perhaps more tangible in the first instance, is in terms of the amplitude, , and phase angle, , of the response, as also shown in the figure above, which are given by
(63)
MATLAB code to calculate the Theodorsen function
Show the code/hide the code
function theodorsen_function_plot_and_save
% Define the range of reduced frequencies
k = linspace(0.01, 10, 500); % Example range from 0.01 to 10 with 500 points
% Preallocate arrays for real and imaginary parts
Ck_real = zeros(size(k));
Ck_imag = zeros(size(k));
% Calculate the Theodorsen function for each k
for i = 1:length(k)
[Ck_real(i), Ck_imag(i)] = theodorsen_function(k(i));
end
% Plot the real part
figure;
subplot(2, 1, 1);
plot(k, Ck_real, ‘b-‘, ‘LineWidth’, 1.5);
xlabel(‘Reduced Frequency k’);
ylabel(‘Real Part of C(k)’);
title(‘Real Part of Theodorsen Function’);
grid on;
% Plot the imaginary part
subplot(2, 1, 2);
plot(k, Ck_imag, ‘r-‘, ‘LineWidth’, 1.5);
xlabel(‘Reduced Frequency k’);
ylabel(‘Imaginary Part of C(k)’);
title(‘Imaginary Part of Theodorsen Function’);
grid on;
% Save the data to a file
data = [k.’, Ck_real.’, Ck_imag.’];
save(‘theodorsen_function_data.txt’, ‘data’, ‘-ascii’, ‘-double’);
end
function [Ck_real, Ck_imag] = theodorsen_function(k)
% This function calculates the Theodorsen function C(k)
% and returns its real and imaginary parts.
% k is the reduced frequency
% Bessel functions of the first kind
J0 = besselj(0, k);
J1 = besselj(1, k);
% Bessel functions of the second kind
Y0 = bessely(0, k);
Y1 = bessely(1, k);
% Hankel functions of the second kind
H1 = J1 – 1i * Y1;
H0 = J0 – 1i * Y0;
% Theodorsen function C(k)
Ck = H1 ./ (H1 + 1i * H0);
% Separate into real and imaginary parts
Ck_real = real(Ck);
Ck_imag = imag(Ck);
end
Theodoresen’s theory is a frequency domain theory; even though the results for the aerodynamic response may be expressed as a function of time, the unsteady forcing must be oscillatory, which is expressed in terms of the reduced frequency. Recall that the reduced frequency is a dimensionless parameter representing the unsteadiness of the flow, defined by
(64)
Physically, the reduced frequency can be interpreted as the ratio of the oscillation frequency to the convective frequency of the flow. Airfoils respond differently to different values of . Low values imply a quasi-steady behavior, while high values indicate highly unsteady conditions. Quasi-steady means that the flow responds almost immediately to any changes in the angle of attack; therefore, the aerodynamic response will be in phase with the airfoil motion.
The effect of the Theodorsen function, which emulates a physical wake lag response, can be appreciated if a pure oscillatory variation in the angle of attack is considered, as before, i.e., ,[5] where is the mean angle of attack and is the angle of attack amplitude of the oscillation. If the flow were quasi-steady, in the sense that flow adjustments were to take place instantaneously, the lift coefficient would be given by
(65)
where for a thin-airfoil in incompressible flow. In the unsteady case, this circulatory part of the airfoil lift coefficient becomes
(66)
or in terms of the real and imaginary parts, then
(67)
The amplitude and phase lag of the lift response is then
(68)
Representative results from the Theodorsen theory are shown in the figure below in terms of lift coefficient versus angle of attack, . Notice that the Theodorsen function introduces an amplitude reduction and phase lag effect on the circulatory part of the lift response compared to the result obtained under quasi-steady conditions.
For , the steady-state lift behavior is obtained, i.e., is linearly proportional to . However, as is increased, the lift plots develop into hysteresis loops, and these loops rotate such that the amplitude of the lift response (slope of the lift curve) decreases with increasing values of reduced frequency. Notice also that these loops are circumvented in a counterclockwise direction such that the lift is lower than the steady value when is increasing with time and higher than the steady value when is decreasing with time, i.e., there is a phase lag accounted for by . Notice that for infinite reduced frequency, the circulatory part of the lift amplitude is half that at , and there is no phase lag angle.
These concepts are readily extended to more complicated boundary conditions, including the effects of pitch rate and combinations of frequencies. The results are obtained by linear superposition.
Apparent Mass Terms
It should also be appreciated that Theodorsen’s theory applies only to the “circulatory” component of the lift. Additional non-circulatory forces (and moments) are produced by inertial effects or “apparent mass” in an unsteady flow. This must be included to round out the unsteady aerodynamic theory if it is to be used in any predictive sense.
The apparent mass terms account for the pressure forces required to accelerate the fluid near the airfoil. For example, for a thin airfoil of chord moving normal to its surface at velocity , the noncirculatory fluid force, , acting on the surface is
(69)
The term is known as the apparent mass and, in this case, is given by
(70)
Therefore, the non-circulatory lift for a planning motion where the airfoil moves normal to its surface with velocity is
(71)
or in terms of lift coefficient, then
(72)
Unsteady (Dynamic) Stall
Under unsteady conditions, airfoils and wings behave differently than in steady flow. In particular, their stall characteristics are very different, and the behavior is called a dynamic stall. This phenomenon has received much research interest because of the unsteady aerodynamic effects produced, which include the shedding of a leading-edge vortex, as shown in the simulation below.
The vortex shedding increases lift and drag, giving powerful nose-down (negative) pitching moments because of the aft-moving center of pressure induced by the shedding vortex, as shown in the figure below, based on measurements and flow visualization. Dynamic stall is known for exhibiting higher levels of maximum lift, drag, and pitching moments and hysteresis effects that can lead to flutter phenomena. Dynamic stall is particularly relevant in designing and analyzing helicopter rotors, wind turbine blades, and specific aircraft configurations.
Engineers and researchers study the phenomenon of dynamic stall to understand the effects better and develop strategies to mitigate any adverse impact on the performance of aerospace systems. In particular, helicopter rotor blades often experience dynamic stall during higher-speed forward flight or during maneuvers. Indeed, the onset of dynamic stall limits the helicopter’s operational flight envelope. Wind turbine blades may encounter dynamic stall during sudden changes in wind conditions, leading to variations in power output and high blade loads that may produce concerning vibration levels. Aircraft wings can also undergo dynamic stall during aggressive maneuvers, which is a behavior that must be considered when establishing the maneuvering flight envelope for military combat aircraft.
Summary & Closure
Unsteady flows, characterized by fluctuations in fluid properties over time, are a fundamental aspect of fluid dynamics that are especially relevant in many engineering and scientific disciplines. Unlike steady flows, unsteady flows are more representative of real-world scenarios where conditions are rarely constant. The study of unsteady flows is essential for designing and optimizing advanced aerodynamic systems, propulsion technologies, and other applications where time-dependent behavior plays a critical role. The practical implications of unsteady flows extend beyond basic fluid dynamics, impacting areas such as turbulence modeling, wave propagation, and the transient behavior of fluid systems. Addressing these complexities requires sophisticated mathematical models and computational methods capable of capturing the transient phenomena. An important step in analyzing unsteady flows is their classification, often achieved using dimensionless parameters like the reduced frequency, which helps to quantify the degree of unsteadiness relative to the system’s inherent time scales. Throughout this chapter, various classic examples illustrate the principles of solving unsteady flow problems. This understanding is crucial for developing accurate predictions and effective solutions in applications where unsteady flows are a key consideration.
5-Question Self-Assessment Quickquiz
For Further Thought or Discussion
- Can you provide real-world situations where unsteady flow is crucial in fluid systems or engineering applications?
- Why is reduced frequency an essential parameter for analyzing fluid flow oscillatory or vibratory motion?
- In what scenarios might the Strouhal number be a relevant parameter to consider?
- How does turbulence contribute to the complexity of fluid dynamics, and why is it often associated with unsteadiness?
- Discuss the practical implications of unsteady flows in engineering design. How might engineers account for unsteady conditions in the design of systems like pipelines, aircraft, or water turbines?
- How might unsteady flows impact the efficiency and performance of propulsion systems, such as those in aircraft or marine vehicles?
Other Useful Online Resources
To learn more about unsteady flows, take a look at some of these online resources:
- A good video explaining the differences between steady and unsteady flows.
- An explanation of unsteady flows, including the effects of thermodynamics.
- An explanation of the conservation of mass in unsteady flow.
- A video presentation on the application of the continuity equation to unsteady flows.
- Joukowsky, N., “Über den hydraulischen Stoss in Wasserleitungsröhren.” (“On the hydraulic hammer in water supply pipes.”) Mémoires de l'Académie Impériale des Sciences de St.-Pétersbourg, Series 8, 9(5), 1-71, 1900 (in German). ↵
- For example, the sedimentation of Sahara dust is a fascinating phenomenon with significant environmental, meteorological, and climatological implications. Sahara dust, composed of fine mineral particles, is lifted into the atmosphere by strong winds and transported across vast distances, often reaching as far as the Americas and Europe. This dust plays a crucial role in various ecological and atmospheric processes. ↵
- In drug manufacturing, sedimentation separates active ingredients from suspension. In the food Industry, sedimentation is used in processes like clarifying juices, where solid particles must be removed from the liquid. ↵
- Theodorsen, T., "General Theory for Aerodynamic Instability and the Mechanism of Flutter," NACA TR 496, 1935. ↵
- Notice that this case does not represent an oscillating airfoil because the pitch-rate terms are not included. ↵