MAE 210B. Homework Solution #6 Winter Quarter, U 2 =r U=r 2 << 1; ) r << U : (1) The boundary conditions written in polar coordinates,

Similar documents
Equilibria of a cylindrical plasma

General Relativity Homework 5

PH126 Exam I Solutions

Quantum Mechanics I - Session 5

7.2.1 Basic relations for Torsion of Circular Members

Electric Potential and Gauss s Law, Configuration Energy Challenge Problem Solutions

15. SIMPLE MHD EQUILIBRIA

Physics Courseware Physics II Electric Field and Force

r cos, and y r sin with the origin of coordinate system located at

Conservation of Linear Momentum using RTT

Section 5: Magnetostatics

Example

Homework Set 3 Physics 319 Classical Mechanics

Solutions to Problems : Chapter 19 Problems appeared on the end of chapter 19 of the Textbook

Errors in Nobel Prize for Physics (3) Conservation of Energy Leads to Probability Conservation of Parity, Momentum and so on

rt () is constant. We know how to find the length of the radius vector by r( t) r( t) r( t)

That is, the acceleration of the electron is larger than the acceleration of the proton by the same factor the electron is lighter than the proton.

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018

GRAVITATION. Einstein Classes, Unit No. 102, 103, Vardhman Ring Road Plaza, Vikas Puri Extn., New Delhi -18 PG 1


ENGI 4430 Non-Cartesian Coordinates Page xi Fy j Fzk from Cartesian coordinates z to another orthonormal coordinate system u, v, ˆ i ˆ ˆi

Class #16 Monday, March 20, 2017

Much that has already been said about changes of variable relates to transformations between different coordinate systems.

556: MATHEMATICAL STATISTICS I

Chapter 28: Magnetic Field and Magnetic Force. Chapter 28: Magnetic Field and Magnetic Force. Chapter 28: Magnetic fields. Chapter 28: Magnetic fields

Chapters 5-8. Dynamics: Applying Newton s Laws

Physics 122, Fall December 2012

, and the curve BC is symmetrical. Find also the horizontal force in x-direction on one side of the body. h C

Moment. F r F r d. Magnitude of moment depends on magnitude of F and the length d

Do not turn over until you are told to do so by the Invigilator.

of Technology: MIT OpenCourseWare). (accessed MM DD, YYYY). License: Creative Commons Attribution- Noncommercial-Share Alike.

PHY2061 Enriched Physics 2 Lecture Notes. Gauss Law

( )( )( ) ( ) + ( ) ( ) ( )

That is, the acceleration of the electron is larger than the acceleration of the proton by the same factor the electron is lighter than the proton.

Microscopic Momentum Balances

EN40: Dynamics and Vibrations. Midterm Examination Thursday March

arxiv: v1 [physics.pop-ph] 3 Jun 2013

Electric Potential. Outline. Potential Energy per Unit Charge. Potential Difference. Potential Energy Difference. Quiz Thursday on Chapters 23, 24.

2. Radiation Field Basics I. Specific Intensity

( ) ( )( ) ˆ. Homework #8. Chapter 27 Magnetic Fields II.

Applied Aerodynamics

Electrostatics (Electric Charges and Field) #2 2010

2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum

Physics 181. Assignment 4

MATH 417 Homework 3 Instructor: D. Cabrera Due June 30. sin θ v x = v r cos θ v θ r. (b) Then use the Cauchy-Riemann equations in polar coordinates

c n ψ n (r)e ient/ h (2) where E n = 1 mc 2 α 2 Z 2 ψ(r) = c n ψ n (r) = c n = ψn(r)ψ(r)d 3 x e 2r/a0 1 πa e 3r/a0 r 2 dr c 1 2 = 2 9 /3 6 = 0.

The study of the motion of a body along a general curve. the unit vector normal to the curve. Clearly, these unit vectors change with time, u ˆ

Gauss Law. Physics 231 Lecture 2-1

Stress, Cauchy s equation and the Navier-Stokes equations

Two Dimensional Inertial Flow of a Viscous Fluid in a Corner

Kinematics of rigid bodies

Geometry of the homogeneous and isotropic spaces

I. CONSTRUCTION OF THE GREEN S FUNCTION

Welcome to Physics 272

Jerk and Hyperjerk in a Rotating Frame of Reference

Classical Mechanics Homework set 7, due Nov 8th: Solutions

VECTOR MECHANICS FOR ENGINEERS: STATICS

MODULE 5a and 5b (Stewart, Sections 12.2, 12.3) INTRO: In MATH 1114 vectors were written either as rows (a1, a2,..., an) or as columns a 1 a. ...

A Level Exam-style Practice Paper

ESTIMATION MODELS USING MATHEMATICAL CONCEPTS AND NEWTON S LAWS FOR CONIC SECTION TRAJECTORIES ON EARTH S SURFACE

A Crash Course in (2 2) Matrices

d 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c

1) Consider an object of a parabolic shape with rotational symmetry z

Physics 107 HOMEWORK ASSIGNMENT #15

THE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2

Thomas Whitham Sixth Form Mechanics in Mathematics. Rectilinear Motion Dynamics of a particle Projectiles Vectors Circular motion

EGN 3353C Fluid Mechanics

Pressure in the Average-Atom Model

On the force and torque on systems of rigid bodies: A remark on an integral formula due to Howe

To Feel a Force Chapter 7 Static equilibrium - torque and friction

Complex Heat Transfer Dimensional Analysis

CHAPTER 2 DERIVATION OF STATE EQUATIONS AND PARAMETER DETERMINATION OF AN IPM MACHINE. 2.1 Derivation of Machine Equations

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 9 Solutions [Theorems of Gauss and Stokes]

THE NAVIER-STOKES EQUATION: The Queen of Fluid Dynamics. A proof simple, but complete.

Lecture 23. Representation of the Dirac delta function in other coordinate systems

An Exact Solution of Navier Stokes Equation

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

ESCI 342 Atmospheric Dynamics I Lesson 3 Fundamental Forces II

EELE 3331 Electromagnetic I Chapter 4. Electrostatic fields. Islamic University of Gaza Electrical Engineering Department Dr.

Physics 1114: Unit 5 Hand-out Homework (Answers)

But for simplicity, we ll define significant as the time it takes a star to lose all memory of its original trajectory, i.e.,

Math 209 Assignment 9 Solutions

Radian Measure CHAPTER 5 MODELLING PERIODIC FUNCTIONS

Polar Coordinates. a) (2; 30 ) b) (5; 120 ) c) (6; 270 ) d) (9; 330 ) e) (4; 45 )

SPH4UI 28/02/2011. Total energy = K + U is constant! Electric Potential Mr. Burns. GMm

AP Physics 1 - Circular Motion and Gravitation Practice Test (Multiple Choice Section) Answer Section

Lecture 8 - Gauss s Law

Vectors, Vector Calculus, and Coordinate Systems

Physics 121 Hour Exam #5 Solution

Answers to test yourself questions

On the integration of the equations of hydrodynamics

Today s Plan. Electric Dipoles. More on Gauss Law. Comment on PDF copies of Lectures. Final iclicker roll-call

15 Solving the Laplace equation by Fourier method

Homework Set 4 Physics 319 Classical Mechanics. m m k. x, x, x, x T U x x x x l 2. x x x x. x x x x

2. Plane Elasticity Problems

A Lattice Energy Calculation for LiH

Particle Systems. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell

is the instantaneous position vector of any grid point or fluid

THE HERMITE POLYNOMIAL & QUANTIZATION OF THE HARMONIC OSCILLATOR

Diffusion and Transport. 10. Friction and the Langevin Equation. Langevin Equation. f d. f ext. f () t f () t. Then Newton s second law is ma f f f t.

Transcription:

MAE B Homewok Solution #6 Winte Quate, 7 Poblem a Expecting a elocity change of oe oe a aial istance, the conition necessay fo the ow to be ominate by iscous foces oe inetial foces is O( y ) O( ) = = << ; ) << : () b The bounay conitions witten in pola cooinates, Dene the steam function > : = ; x; = ; () > : = ; ; = : (3) ; = x ; (4) such that continuity is satise. We know that the steam function satises bihamonic equation fo two-imensional Stokes ow ( ) = : (5) The bounay conitions (taking efeence steam function = at the walls) ae > : = ; = ; > : = ; = ; = x; ; (6) = : (7) c Fom equation an B.Cs, we see that the aiables which can inuence the steam function ae = F (; ; ; ): (8)

The numbe of aiables with imensions ( = [L T x ], = [LT x ] & = [L] ) ae n = 3 an the numbe of inepenent imensions ([L] & [T ]) ae k =, whence, the numbe of goup is n x k =. Equation 5, witten in pola cooinates is = f (; ) (9) 4 4 4 3 + 4 + 4 + 4 x 3 3 3 x Substituting = f (; ) in the aboe equation, we get + 4 4 + 3 = () f i + f + f = : () Assuming the solution of the fom f exp(), we get = ~i; ~i, theefoe the solution fo f is gien by f () = A sin + B cos + C sin + D cos ; () f () = (C x B) sin + (A + D) cos x D sin + C cos ; (3) f () = ; f () = x; f () = ; f () = : (4) It's tiial to see B = an D = x(a + ) using the lowe wall bounay conitions, an the othe bounay conitions leas to Soling fo A an C, we get f () = A sin + C sin x (A + ) cos = ; (5) f () = C sin x (A + ) sin + C cos = cos : (6) A = x Afte a eaangement, we n x sin ; B = ; C = x sin cos x sin ; D = sin x sin (7) f (; ) = x sin [ sin sin( x ) x ( x ) sin ] : (8) Velocities (o not epen on ) ae f (), = xf followe fom thei enitions, x sin = x x sin fsin [sin( x ) x cos( x )] + [sin x ( x ) cos ]g ; (9) f sin sin( x ) x ( x ) sin g : ()

Pessue can iectly be obtaine fom steamfunction using the elation! = x only fo two-imensional ow) as shown below, e z (tue = p = ; ) p = z ( e z ); () + + witing the aboe equation in e an e iection gies Integating the st equation, p = ( ) = p = x ( ) = = f + f ; () (f + f ); (3) (f + f ): (4) p(; ) = x (f + f ) + G(); ) p = x (f i + f G ) + an substituting p= into the secon equation gies (5) Thus G = f i + f + f = ; ) G = constant = p : (6) p(; ) x p = x (f + f ): (7) Since f () = x3c sin x(a+3d) cos +D sin xc cos we n f +f = xc sin xd cos, theefoe p(; ) x p = (C sin + D cos ); (8) sin + sin sin( x ) p(; ) x p = x sin : (9) e Pessue on the wipe = is p(; ) x p = sin x sin : (3) To calculate the iscous stesses on the wall, we ene a nomal ecto on the wipe n = xe such that the iscous stess on that iection is y n = x e x e ealuate at =. = + ; = + : (3) 3

It is easy to see that is ientically zeo = x f f + = (3) an theefoe only the tangential iscous stess contibutes to the foce, which is = (f + f ) = The iscous stess on the wipe =, p = cos x sin y = n= cos x sin cos( x ) x sin : (33) x sin : (34) We see that both the pessue an the iscous stess acting on the wipe ay as x an hence the foce will ay as ln. Theefoe the foce acting on the wipe fom the oigin to any nite aius will be logaithmically innite. But in eal life, the wipe can nee make complete contact with the plane suface an theefoe a nite stess at the oigin, esulting in a nite foce (which will epen on the gap with an the geomety of intesection). This poblem was st sole by G.I. Taylo as a moel fo scape an he shows that the nomal an tangential stess on the scape euces to that of Lubication theoy when <<. Poblem A moe inole poblem woul be the Stokes ow aoun a ui sphee of abitay iscosity, whee both insie an outsie motion of the sphee has to be sole simultaneously. But hee the iscosity of insie ui is negligible, thus it is enough to consie only outsie ow. We stat with the enition of Stokes steam function, sin ; = x sin The fee steam bounay conition at innity an the no penetation conition at = a still hols fo this poblem as in ow oe a soli sphee, but the non-slip conition nees to be eplace with the stess fee conition in the tangential iection, = at = a i.e., =a = because at = a, 8. + = ; ) =a (35) = (36) Taylo, G.I., 96. On scaping iscous ui fom a plane suface. Miszellangen e Angewanten Mechanik (Festschift Walte Tollmien), pp.33-35. Taylo, G.I., 97. The Scientic Papes of Si Geoey Ingam Taylo: Volume 4, Mechanics of Fluis: Miscellaneous Papes (Vol. 4). Cambige niesity Pess. 4

Stokes steam function o not satisfy bihamonic equation, but satises the following equation, with bounay conitions, E 4 = " sin + = a : = ; sin # = (37) = (38)! : = sin : (39) Assuming the solution of the fom = f ( ) sin, the equation becomes x! f = ; (4) Substituting f ( ) into the aboe Eule equation gies us the oots = x; ; ; 4. The geneal solution is an the bounay conitions ae A f ( ) = = a : f = ; + B + C + D 4 (4) = (4) f! : f = : (43) pon soling, we get A =, B = xa=, C = = an D =. Theefoe Hence the elocity, f ( ) = = ( x a ); (44) ( x a) sin : (45) ( x a) cos ; (46) Equation fo pessue is E = = x x a sin : (47) p = ; (48) f x f a sin = sin (49) 5

hence (see Acheson, pp. 4) Integating of st equation, p = sin (E a ) = p = x sin (E a ) = p(; ) = x a cos + G(); ) p a G = sin + 3 cos ; (5) 3 sin : (5) an substituting into the secon equation gies G= = ) G = p. Theefoe the pessue is p(; ) = p x a cos ; p(a; ) x p = x a (5) cos : (53) To calculate iscous foce, we ene a nomal ecto on the sphee n = e, theefoe y n = ( ; ; ). But by enition = at = a, only the nomal iscous stess contibutes to ag. = = cos : (54) =a =a a Dag foce is calculate in the iection of motion so that net stess acting on the sphee in that iection is xp(a; ) + p + j =a cos an with the ieential element on the suface of the sphee = a sin, F D =Z Z = 3aZ xp(a; ) + p + =a Z a sin cos ; (55) sin cos ; (56) F D = 4 a: (57) =3 of F D comes fom iscous stess an =3 fom pessue ag. Poblem 3 The conition neee fo the ow to be ominate by iscosity is O( y ) O( ) R =R R=R << ; ) R << : (58) In this low Re limit, we can use the Stokes equation = xp + in spheical pola cooinates. Due to the symmety of the poblem, we hae = = =, thus we hae (; ) an p(; ). The momentum an momentum euces to p = ; p = ; ) p = p : (59) 6

No centifugal foce exists in Stokes ow, hence the constant pessue. The momentum equation then becomes, Bounay conitions ae + sin sin x sin = : (6) = R : = R sin ; (6)! : = : (6) Looking at the bounay conition, we assume the solution of the fom = f ( ) sin. With this equation (6) euces to f + f = ; f (R) = R; f () = : (63) The geneal solution is B f ( ) = A + It's tiial to see A = an B = R 3. Hence the elocity is = R sin R The only non-zeo stess component which can contibute to toque is = : (64) = x3 sin R : (65) 3 ; (66) =R = x3 sin : (67) Thus, the toque on the suface is obtaine fom the pouct of stess an the aial istance R sin fom the axis of otation, integate oe the suface with the ieential suface element = R sin, T =Z Z = x6 R 3 Z R sin R sin ; (68) sin 3 ; (69) T = x8 R 3 : (7) 7