φ(r, θ, t) = a 2 U(t) cos θ. (7.1)

Similar documents
Homework Two. Abstract: Liu. Solutions for Homework Problems Two: (50 points total). Collected by Junyu

Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics

Fluid Mechanics Prof. T. I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay

Department of Mechanical Engineering

Applied Fluid Mechanics

1. Fluid Dynamics Around Airfoils

SELF-PROPELLED ANGUILLIFORM SWIMMING: SIMULTANEOUS SOLUTION OF THE TWO-DIMENSIONAL NAVIER STOKES EQUATIONS AND NEWTON S LAWS OF MOTION

F11AE1 1. C = ρν r r. r u z r

(3) BIOMECHANICS of LOCOMOTION through FLUIDS

Math 575-Lecture Failure of ideal fluid; Vanishing viscosity. 1.1 Drawbacks of ideal fluids. 1.2 vanishing viscosity

Water is sloshing back and forth between two infinite vertical walls separated by a distance L: h(x,t) Water L

- Marine Hydrodynamics. Lecture 14. F, M = [linear function of m ij ] [function of instantaneous U, U, Ω] not of motion history.

FLUID MECHANICS. Chapter 9 Flow over Immersed Bodies

All that begins... peace be upon you

SPC Aerodynamics Course Assignment Due Date Monday 28 May 2018 at 11:30

Detailed Outline, M E 521: Foundations of Fluid Mechanics I

Vortex motion. Wasilij Barsukow, July 1, 2016

Offshore Hydromechanics Module 1

MAE 101A. Homework 7 - Solutions 3/12/2018

AA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow

7 EQUATIONS OF MOTION FOR AN INVISCID FLUID

PEMP ACD2505. M.S. Ramaiah School of Advanced Studies, Bengaluru

Flight Vehicle Terminology

UNIVERSITY of LIMERICK

Friction Factors and Drag Coefficients

Fluid Dynamics Problems M.Sc Mathematics-Second Semester Dr. Dinesh Khattar-K.M.College

equation 4.1 INTRODUCTION

OFFSHORE HYDROMECHANICS OE 4620-d

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

Fluid Dynamics Exercises and questions for the course

2. Mass, Force and Acceleration

BLUFF-BODY AERODYNAMICS

Marine Hydrodynamics Prof.TrilochanSahoo Department of Ocean Engineering and Naval Architecture Indian Institute of Technology, Kharagpur

Lecture-4. Flow Past Immersed Bodies

Student name: This is a closed book examination. You are allowed 1 sheet of 8.5 x 11 paper with notes.

Given a stream function for a cylinder in a uniform flow with circulation: a) Sketch the flow pattern in terms of streamlines.

External Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.

Lab Reports Due on Monday, 11/24/2014

3.8 The First Law of Thermodynamics and the Energy Equation

AE301 Aerodynamics I UNIT B: Theory of Aerodynamics

Passive locomotion via normal-mode coupling in a submerged spring mass system

Fundamentals of Fluid Dynamics: Waves in Fluids

Rotational Kinetic Energy

5.1 Fluid momentum equation Hydrostatics Archimedes theorem The vorticity equation... 42

Some Basic Plane Potential Flows

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

R.N. Govardhan and J.H. Arakeri*

1. Introduction, tensors, kinematics

3.5 Vorticity Equation

Control Volume Analysis For Wind Turbines

(Jim You have a note for yourself here that reads Fill in full derivation, this is a sloppy treatment ).

The Divergence Theorem Stokes Theorem Applications of Vector Calculus. Calculus. Vector Calculus (III)

Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation

Lifting Airfoils in Incompressible Irrotational Flow. AA210b Lecture 3 January 13, AA210b - Fundamentals of Compressible Flow II 1

Passive locomotion via normal mode. coupling in a submerged spring-mass system

Physics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top

The Bernoulli Equation

FLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes

GFD 2013 Lecture 4: Shallow Water Theory

Vortex Induced Vibrations

AERODYNAMICS STUDY NOTES UNIT I REVIEW OF BASIC FLUID MECHANICS. Continuity, Momentum and Energy Equations. Applications of Bernouli s theorem

Pressure in stationary and moving fluid Lab- Lab On- On Chip: Lecture 2

Forces Part 1: Newton s Laws

Effect of wavelength of fish-like undulation of a hydrofoil in a free-stream flow

Lecture 9 - Rotational Dynamics

CEE 3310 Control Volume Analysis, Oct. 7, D Steady State Head Form of the Energy Equation P. P 2g + z h f + h p h s.

Vorticity Equation Marine Hydrodynamics Lecture 9. Return to viscous incompressible flow. N-S equation: v. Now: v = v + = 0 incompressible

ν δ - 1 -

Momentum (Newton s 2nd Law of Motion)

Exam Question 6/8 (HL/OL): Circular and Simple Harmonic Motion. February 1, Applied Mathematics: Lecture 7. Brendan Williamson.

6. The Momentum Equation

Tutorial 10. Boundary layer theory

Mathematical Modelling and 3D Simulation of a Virtual Robotic Fish

PC 1141 : AY 2012 /13

UNIVERSITY OF MALTA G.F. ABELA JUNIOR COLLEGE

On the generation of a reverse Von Karman street for the controlled cylinder wake in the laminar regime

Simulations of optimized anguilliform swimming

Module 2: External Flows Lecture 12: Flow Over Curved Surfaces. The Lecture Contains: Description of Flow past a Circular Cylinder

Classical Mechanics. FIG. 1. Figure for (a), (b) and (c). FIG. 2. Figure for (d) and (e).

Aerodynamics. High-Lift Devices

1. Introduction - Tutorials

Symmetries 2 - Rotations in Space

The vorticity field. A dust devil

A Biologically Inspired Computational Study of Flow Past Tandem Flapping Foils

Class 11 Physics NCERT Exemplar Solutions Motion in a Straight Line

In this section, mathematical description of the motion of fluid elements moving in a flow field is

Bluff Body, Viscous Flow Characteristics ( Immersed Bodies)

τ du In his lecture we shall look at how the forces due to momentum changes on the fluid and viscous forces compare and what changes take place.

Iran University of Science & Technology School of Mechanical Engineering Advance Fluid Mechanics

4.1 LAWS OF MECHANICS - Review


Wings and Bodies in Compressible Flows

FLOW SEPARATION. Aerodynamics Bridge-Pier Design Combustion Chambers Human Blood Flow Building Design Etc.

MAE 3130: Fluid Mechanics Lecture 7: Differential Analysis/Part 1 Spring Dr. Jason Roney Mechanical and Aerospace Engineering

Chapter 10. Rotation of a Rigid Object about a Fixed Axis

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Concept Question: Normal Force

11.1 Introduction Galilean Coordinate Transformations

MECHANICAL PROPERTIES OF FLUIDS

Performance. 5. More Aerodynamic Considerations

Transcription:

BioFluids Lectures 7-8: Slender Fish Added Mass for Lateral Motion At high Reynolds number, most of the effort required in swimming is pushing water out of the way, that is our energy goes in providing kinetic energy for a sizable volume of fluid. For a thin, streamlined fish this effect for the sideways movements during active swimming is more significant than combatting the inevitable drag due to the net forward motion. To a good approximation the flow can be treated as inviscid, and irrotational outside the wake region. We consider a fish thin in the z-direction swimming in the negative x-direction. From observation it performs small amplitude motions in the z-direction, and in so doing, propels itself forwards surprisingly fast. Now when it moves its body sideways, it sets some fluid in motion. Thus its lateral movements seem to carry with them a certain amount of added mass of water. We begin by formalising this idea. Consider a circular cylinder of radius a moving with velocity Ut) through fluid of density ρ. If we ignore separation, and treat the flow as potential, then we can write u = φ where in terms of cylindrical polar coordinates r, θ, z), φr, θ, t) = a 2 Ut) cos θ. 7.1) r The fluid velocity at the cylinder surface is U cos θ, U sin θ, ) and slips over the solid boundary which moves with U cos θ, U sin θ, ). and the associated pressure field can be found from the Bernoulli equation p p = ρ φ t 1 2 ρ u 2 = a 2 ρ du cos θ 1 dt r 2 ρu 2. 7.2) Now integrating the pressure around the cylinder, we get the force per unit length F, exerted by the cylinder on the fluid. It is in the x-direction with magnitude F = πa 2 ρ du dt. 7.3) Clearly we can regard F as accelerating a mass per unit length m = πa 2 ρ of fluid. We can also calculate the total kinetic energy of the fluid, E = 1 2 ρ φ 2 dv = 1 2 ρ φ φ n ds = 1 2 ρu 2 πa 2 = 1 2 mu 2, 7.4) r>a r=a which once again suggests an apparent added mass m of fluid moving with the cylinder. More generally, we can show that a body of any shape moving with velocity U gives rise to fluid kinetic energy E = 1 2 M iju i U j where M is the added mass tensor, 7.5) and depends only on the instantaneous body shape. To a good approximation the appropriate added mass m for a given body cross-section can be found from two-dimensional 1

theory i.e. by finding the added mass for a cylinder of the local cross-section. This approximation requires the wavelength of the body movements to be significantly greater than the cross-section depth, s. Given this, it is only really the depth s of the cross-section which is important. Potential flow calculations show that if we write m = 1 4 πs2 ρβ then β 1 for most shapes, 7.6) the relation being exact β = 1) for ellipses of any aspect ratio. Small amplitude Elongated Body Theory We now investigate the consequences of Lighthill s elongated body theory for the design of a fish which is to swim efficiently. We have argued that the motion of a slender fish could be well understood by considering the sum of the added mass effects of small amplitude lateral motions. We consider a fish in the x, y) plane, with a depth sx) in the y-direction for < x < L. We assume it undergoes displacements z = hx, t) as a result of which it swims with a constant velocity U,, ). To begin with, we work in a frame in which the fish has zero mean motion, so that the fluid at infinity moves past it with uniform speed U. Then the kinematic boundary condition equality of the normal components of fluid and solid velocity) requires = D Dh z h) = w = Dt Dt h t + U h x, 7.7) where wx, t) is the z-velocity of the fluid at the fish surface. Note that for an inviscid flow this is different from the z-velocity of the fish, which is W = h t. 7.8) Then using the local cylinder approximation, the total z-momentum of the fluid added mass is mx)wx, t). As this is time-dependent, the fish must exert a force F z in the z-direction, which is F z = Dmw). Dt 7.9) Now this force is provided by the fish, which has to work at a rate W F z in order to produce it. Thus the total rate of working of the fish over its entire length is P = = = d dt W F z dx = t + U x [ h t t + U ) mw) dx x ) dx mw w t t dx ) [ dx + U L. t ) t 1 2 mw2 2 7.1)

The first term on the RHS of 7.1) clearly has zero time-average for any quasiperiodic motion, and so will not contribute directly to any mechanism enabling steady swimming against a non-zero drag. Now in 7.6) we argued that the added mass mx) = 1 4 πρs2, where sx) is the total depth of the fish cross-section. In particular, this means that m = at the front of the fish. Many fish have a large caudal fin at their rear, and so m will be sizable there. Even if their shape tapers to a point in a nicely streamlined shape, we know some separation will occur, and a wake region will commence at some effective non-zero sx). Intriguingly, it appears as though an important factor in swimming is the the value of muww at x = L, and already we might surmise that a large caudal fin might be a good idea. The term muww can be interpreted as the rate of working required to shed lateral momentum into the wake. We now derive an alternative expression for the work-rate P, Comparison with the result 7.1) then enables us to calculate the forwards thrust of the fish, which must balance the drag in the wake, at least on average. Let us now work in a stationary frame, through which the fish swims with speed U. The sideways velocity components w and W are the same as before. Then if the fish exerts a force T in the swim direction, the rate at which this force does work is T U. This is included in the total workrate of the fish, P, along with the rate of increase of kinetic energy of the fluid and the energy flux into the wake. Consistently, we only consider the z-velocity of the fluid, so that [ P = T U + d L 1 dt 2 mw2 dx + [U 12 mw2 L. 7.11) Combining the two expressions 7.11) and 7.1), we obtain UT = U [mww 12 w) + d ) h mw dt t w dx 7.12) so that using 7.7) we obtain the time-dependent thrust T [ T = [mww 12 w) d L dt x dx. 7.13) Denoting time-averages by, we obtain the mean thrust and power from 7.13) and 7.11) T = [ [ mww 1 2 w2 ), P = muww. 7.14) Now this mean thrust must be positive, T, as it must balance the mean drag, D, which we discussed in lecture 2. It is quadratic in the swim speed U, and proportional to the cross-sectional area in the direction of swim. This cross-section is increased by lateral movements, so swimming strategies must balance the increased thrust and increased drag obtained by lateral movements. Some fish adopt a swim/glide cycle, accelerating with rapid wiggles, followed by a long flat glide. 3

Estimating the efficiency of swimming We can estimate the efficiency, η, of a swimming cycle by comparing UT and P η = Useful work rate Total work rate = UT P = 1 1 w 2 2 ww ). 7.15) For consistency, we have η 1, so that w and W are positively correlated. We thus see that the most efficient swimming occurs when w W or Uh x h t. However, that leads to small mean thrust and limits the swim speed. Let s take a simple travelling wave as an example. If h = h sin kx V t L), then wl, t) = V U)kh cos kv t and W L, t) = kh V coskv t). 7.16) Then the mean thrust T and η are given by T = 1 2 ml)k2 h 2 V U)V 1 2 V U)) = 1 4 ml)k2 h 2 V 2 U 2 ), 7.17) and η = 1 1 2 V U) 2 /2 V V U)/2 = U + V 2V. 7.18) We see that efficient swimming occurs when the swim speed equals the wave speed, U = V, but this leads to zero thrust and so cannot balance the drag. In practice, we expect the swim speed to be close to, but less than, the wave speed. Evolution towards a hydrodynamically perfect fish. It is clear that conditions at the tail, or caudal fin are vital when swimming. It is therefore no surprise that efficient swimmers have developed a large caudal fin, with resultant large ml). A further evolutionary development has been the change from anguilliform to carangiform or thunniform) swimming. The anguilliform eel-like) mode has essentially the entire body undulating, whereas in carangiform mode the travelling wave has sizable amplitude only towards the rear of the body. This has a potential disadvantage known as recoil. As the caudal fin flaps, it generates a moment about the centre of mass which could cause the fish nose to bend, increasing the drag and cancelling out any advantage. This effect can be minimised by morphological changes, in which the start position of the travelling wave corresponds to a thin cross-section, and very low mx) at that point. The high forward mass leads to very limited recoil of the front body. Strictly speaking, we should include a sideways motion and rotation in the theory to date, using a displacement h x, t) = hx, t) + xω t) + W t), 7.19) rather than the h we have used so far. Ω and W should then be determined by Newton s laws for the entire fish. However, W and Ω are not large in small amplitude motion. The 4

theory can also be extended to large amplitude displacements without too much effort, but we won t investigate that here. Observations of a variety of fish shapes and behaviour indicate two changes as the swimming mode becomes more carangiform. Simultaneously, there is a reduction in extent of the dorsal fin, which is much used in anguilliform motion, a decrease in the minimum cross-section forward of the caudal fin and thirdly a thickening of the front portion of the body. All these changes can be understood in terms of the above theory. A further evolutionary change is the development of the lunate tail. The shedding of vorticity into the wake takes place at the tips of the caudal fin; there is no particular advantage in having a continuous tail, which can increase the drag a little. Some of the fastest swimmers have a caudal fin which resembles swept-back wings. Some of the preceding theory requires mild modification for such designs. So what is the caudal fin doing in these fast swimmers such as the tuna? As it waggles it sheds vortices into the wake, from alternate sides, which at first sight look like a Karman vortex street. However, the sense of each vortex is opposite from that of the street, so that there is a rearward flux of momentum in the wake. Schooling: Schools of fish can use the wake structure to improve their swimming efficiency, in a manner reminscent of slip-streaming. While it would be a mistake for two fish to swim one behind each other, net benefit can derive from a diamond shape pattern, as in the figure. Other advantages of schooling relate to interactions with predators. Another way of understanding what is actually happening in the wake of one of these fast swimming fish is to regard the caudal fin as acting like a wing, shedding vorticity in order to generate circulation and a resultant lift force. As the tail oscillates, this force always has a forward component the fish, quite literally, flies through the water. This suggests we should begin to consider lift generation and animal flight... which we shall do next time. There be four things which I know not: The way of an eagle in the air; the way of a serpent upon a rock; the way of a fish in the midst of the sea; and the way of a man with a maid. Proverbs 3:18-19. 5