CEE 262A H YDRODYNAMICS

CEE 262A HYDRODYNAMICS Lecture 15 Unsteady solutions to the Navier-Stokes equation

Effects of unsteadiness Up to now we have assumed “ Steady State ” Poiseulle-Couette Flow H i.e. It takes a time >> for flow to reach a steady The “generic” estimate of the timescale for viscous effects to act over a lengthscale L is or we can say in time t the distance over which viscous forces are important is

3 1 2 Stokes' first problem “Infinite” Moving Plate (i) Flow is invariant in x1 direction (ii) No x2 component

(iii) From (i) x1-Mom: x3-Mom: BC’s IC Need to solve

Can we convert PDE into an ODE by defining a new variable which combines x and t ? This is a self-similar solution Let and (a) Convert and to system U0

- (9) BC’s - (10) Substitution into the governing PDE Want ODE in f and Need to get rid of t or - (8)

Let Choose

From (10) From (9) • By Definition (see, e.g. Abramowitz & Stegun)

is scaled by , the viscous length at time t • By Definition: Profile has same shape but height over which velocity varies grows like t1/2 As we argued before

E.g. u1/U0=0.05 when This formula involving the error function looks like this:

What it looks like in terms of x3 (n=1 m2/s) at different times (every 0.1 s) Note that if we had scaled the coordinate x3 for each profile by (2nt)1/2, we get the single “self-similar” profile seen on the previous slide

In terms of the vorticity: What we have calculated is the diffusion of a spike (sheet) of vorticity initially located at x3=0 into the fluid above, If we calculate the circulation around a contour that looks like this (but infinite in the x3 direction): C x3 A U0 Thus, the total vorticity is constant l0

For x3=0 the flux = 0 except for t = 0 All the vorticity enters the fluid at t = 0 due to the impulsive onset of plate motion What about the vertical flux of vorticity from the wall into the fluid? Vertical Flux Summary: Setting the plate in motion imparts a fixed amount of vorticity to the fluid – This vorticity diffuses into the fluid, spreading vertically such that the layer containing a bulk of the vorticity is thick at time t

Unsteady Couette flow The analysis of Stokes’ first problem showed a flow that was self-similar such that when viewed with a vertical scale that grew with time, the flow always looked the same. What happens if we confine the flow in a finite depth region? Qualitative: If t << D2/n - erf solution near x3=D, zero elsewhere If t >> D2/n - steady Couette flow as before

Exact method: Use of separation of variables/Fourier series This part is the steady solution, it satisfies the top and bottom BCs, but not the initial condition (at least by itself) • This part is the deviation from the steady solution. It has two properties: • With steady solution it satisfies the initial conditions. • It satisfies the BCs that are left after we take away the steady solution, i.e Note: Both parts satisfy the p.d.e

Thus, we must solve subject to: Approaches: Laplace transforms, numerical methods, or Fourier series

We use a Fourier sine series with the amplitude of each Fourier component assumed to be a function of time: • Why this form? (which works because of linearity) • The sine terms each satisfy the top and bottom boundary conditions on u’ • Any function of x3 can be described as a sum of these sines. This property is known as “completeness”. • The sines are “orthogonal” (we’ll see why this is useful). This means that:

If we substitute our expansions into the governing PDE we find that How do we find Un(0)? At t = 0, we have the condition

Thus we need to find the Un(0) to make the sum = the negative of the steady state solution. This is done by multiplying both sides of this expression by sin(mpx3/D) and integrating from 0 to D, i.e. using the orthogonality condition:

Thus, when all is said and done we find that t t t =0 to 0.3

Again we assume so that Stokes' second problem A flow closely related to the flow of Stokes' first problem is Stokes' second problem! Now we consider the flow induced above a flat plate at x3=0 that oscillates according to

First, what happens? If we move the plate to the right, we impart negative vorticity to the fluid – This diffuses upwards Now, when the plate reverses direction. A new front begins to diffuse having opposite sign vorticity. Thus as the plate oscillates for some time, we expect that the vorticity from successive fronts eventually cancels out one-another far enough away from the boundary. Near the boundary, we expect to see only the current forcing’s ‘front’. old front new front

To confirm this picture, let’s look at the forced or periodic soln. : We choose dimensionless variables using To create a parameter – free eqn. we set Expressing the fact that we expect the dominant length scale will be the distance vorticity can diffuse in one period Thus with Real part of complex #

We now set (and drop *) u is complex The P.D.E. gives We can eliminate one of these values of k by insisting that Eliminate since unbounded.

Here

These two components look like this (pic in page V39)

CEE 262A H YDRODYNAMICS