Squirming Sphere in an Ideal Fluid
|
|
- Linda Malone
- 5 years ago
- Views:
Transcription
1 Squirming Sphere in an Ideal Fluid D Rufat June 8, 9 Introduction The purpose of this wor is to derive abalytically and with an explicit formula the motion due to the variation of the surface of a squirming body in a fluid. Such a derivation has been done before, but only in two dimensions, and there is no literature on any attempts to generalize to three dimensions, which is cosiderably more involved algebraically. But utilizing software such as Mathematica, the derivation is quite straightforward. The primary motivation behind this wor is the search for innovative propulsion mechanisms. Surface variations may not lead to a very high speed, but they provide higher maneuverability and efficiency and are the primary source of motion of aquatic animals. Two types of fluid are studied - an Eulerian (potential fluid where all viscous effects are ignored and inertia is the only factor, and a Stoesian (creeping fluid where all inertial effects are ignored and viscous friction is the only factor. We consider only those cases because they linearize the equations governing the fluid motion, and permit analytic solutions. Hydromechanical Connections Consider a body immersed in a fluid. We regard its shape as a shape manifold M. And the configuration manifold Q = M G represents the complete position and orientation of the body in space. Q = M G can be thought of as a trivial bundle with a group structure G, and an action of G on Q given by Φ h (r, g = (r, hg for (r, g Q and h G. An element of the tangent space is given by q = (ṙ, ġ T Q. We are intererested in computing a connection on the configuration bundle Q that will allow us to compute the fiberwise translations resulting from cyclical changes in shape.
2 A connection is a vector valued one-form Γ : T Q g on Q. We follow (Kelly, 998 in defining two types of connections: Mechanical connection Γ mech : (q, q I (qj(q, q where I : g g is the loced inertia tensor such that I(qξ, η = ξ Q (q, η Q (q KE and I(r, g = Ad g I loc(rad g and J : T Q g is the momentum map. Stoes connection Γ Stoes : (q, q V (qk(q, q where K : T Q g is a momentum map K(q, q, ξ = F (q, q, ξ Q (q associated with a dissipative force F, and V : g g is a viscosity tensor such that V(qξ, η = F (ξ Q (q, η Q (q and V(r, g = Ad g V loc(rad g Kelly uses the following local expressions Γ mech = Ad g (g ġ + A mech ṙ Γ mech = Ad g (g ġ + A Stoes ṙ and for zero initial momentum g ġ = A mech ṙ g ġ = A Stoes ṙ
3 3 z F (θ, ϕ θ y ϕ x Figure : Coordinate system. which completely describe the evolution of a Potential and Stoesian swimmers respecitlvey. 3 Squirming Sphere Consider a three dimensional nearly spherical body whose surface is described by the following equation: r = F (θ, ϕ; s,..., s n Such a body has already been studied previously numerically (Rufat, 5, and the current method that allows us to obtain explicit formulas for the displacement as will be later seen, are valuable tools for validating numerical simulations. As an example we consider the propulsive movement of a roughly spherical device whose boundary shape is modulated by a small amount. We parametrize the boundary in terms of the shape variables s i (t, and express it as a sum of an unperturbed base shape and a perturbed part: The configuration manifold is F (θ; s = + ɛ (s cos θ + s cos 3θ Q = R SE(3
4 4 z F(θ, s s s Figure : Boundary perturbations. 3. Potential Flow Irrotational potential flow is best described in terms of the velocity potential, which satisfies with the Neumann boundary condition φ = where u is the prescribed boundary motion. φ n = u n ( For an axisymmetric problem, such as this one, the solution of the Laplace equation can be expressed as a series φ(r, θ = = ( A r + B r + P (cos θ where P (x are the Legendre polynomials satisfying the orthogonality condition: ˆ + P m (xp n (xdx = n + δ mn
5 5 We require that u as r, and therefore we must tae A = for. φ(r, θ = = a P (cos θ r + a Φ The problem of solving for φ is reduced to finding the coefficients a. First, in order to enforce the boundary conditions on the perturbed surface, we need to compute the normal n at each point. The two tangential components on the surface are = t = F (F ˆr = θ θ ˆr + F ˆθ t = F (F ˆr = ˆr + F sin θ ˆϕ ϕ ϕ and therefore the normal component can be evaluated as the cross product of the two n = t t = F sin θ ˆr F F θ sin θ ˆθ F F ϕ ˆϕ Note that we need not normalize the expression for n since it occurs on both sides of the boundary condition equation (. Expanding the spherical coordinates (r, θ, ϕ of n in terms of ɛ: n = n ( +ɛn ( +ɛ n ( (c[θ]s + c[3θ]s = + ɛ s[θ]s + 3s[3θ]s + ɛ (c[θ]s + c[3θ]s (c[θ]s + c[3θ]s (s[θ]s + 3s[3θ]s where c[θ] = cos θ and s[θ] = sin θ. Next, we expand the potential in a power series in terms of the small perturbation paramater ɛ. φ = φ ( + ɛφ ( + ɛ φ ( ɛ i φ (i +... where φ (i = a (i Φ = Plugging bac in, we expand both sides in terms of ɛ φ n = u n = i= ɛ i = ɛ i u (i i= a (i f (i
6 6 and equate terms of equal power = a (i f (i = u (i where the first few terms are = = η = 3 ( 3η 3 = ( 3η 5η 3 4 = 5 ( 3 3η + 35η = 3 ( 5η 7η η = 7 ( 5 5η + 35η 4 3η = ( 35η 35η η 5 49η 7 where η = cos θ By the linearity of Laplace s equation (using superposition one can write following Kirchoff, Then, up to second order in ɛ φ = φ z ż + φ s ṡ + +φ s ṡ φ z = P ( r + ɛ 9P s r + 9P s 7r 3 6P 3s 5r 4 48P 4s 35r ( 5 + ɛ (P 3339s 499s 45r 78P s s 735r 3 + 6P ( 3 38s 45s 95r P 4s s 3475r 5 6P ( 5 9s 9s 637r 6 64P 6s s r 7 5P 7s 43r 8 ( P φ s = ɛ 3r 4P ( 9r 3 + ɛ 4P s 3P s 5r 35r + 8P s r 3 + P 3s 5r 4 5P 4s 55r 5 64P 5s 63r 6 ( 3P φ s = ɛ r P ( 3 5r 4 + ɛ 34P s 39P s 35r 35r 3P s 35r 3 + P 3s 5r 4 + 8P 4s 95r 5 8P 5s 89r 6 856P 6s 67r 7 The total inetic energy is given by
7 7 L = = π ˆ F ˆ π ˆ dv φ dφ = π 5 ż ɛ π 5 +ɛ π 5 ˆ π ˆ ˆ sin θdθ r dr φ dη r dr φ ( 49 5 s ż + 6 ( 5ṡż 3 89ṡ + 4 5ṡ + ( s ṡ s ṡ ( 576 ż + 75 s s ż + O(ɛ 3 This Lagrangian is invariant under translation along z, so its momentum map is J = π 5 ż ɛπ 5 ( 49 5 s ż + 3 +ɛ π 5ṡ 5 ( s ṡ 68 ( s ṡ + 75 s s ż +O(ɛ 3 The loced inertia tensor is given by Therefore, I(s, s : ξ π ( 5 ξ ɛ 49 ( s + ɛ 75 s s + O(ɛ 3 Γ mech = I J = dz ɛ 3 ( 69 5 ds + ɛ 55 s ṡ s ṡ 49 5 s 3 + O(ɛ 3 5ṡ = dz ɛ 3 ( 8 5 ds + ɛ 5 s ṡ s ṡ + O(ɛ 3 Self-propulsion of the sphere requires that its motion remain horizontal with respet to this connection, therefore 3. Stoes Flow ż = ɛ 3 ( 68 + ɛ 5ṡ 45 s ṡ 8 5 s ṡ Thes solution of creeping-motion can be reduce in two dimensions to the biharmonic equation in terms of the stream function 4 ψ = which satisfies the no-slip and no-penetration boundary condition on the surface
8 8 u = (u r, u θ = ( F, Following (Leal, 7 p.46 the general solution is ψ = (A r +3 + B r + + C r + D r Q (η = where Q (x are polynomial functions defined in terms of the Legendre polynomials Q (x = and satisfiying the orthogonality condition ˆ x P (xdx ˆ + Q m (xq n (x x dx = n(n + (n + δ mn Again we want the velocity to vanish at infinity, so we will neglect the first two terms ψ = (a r + b r Q (η = u r = r ψ η u θ η = r ψ r Again we will exploit the linearity of the equation by using Kirchoff s method to express teh streamfunction as ψ = ψ z ż + ψ s ṡ + +ψ s ṡ ψ z = ( r 3r [ ( s Q + ɛ Q 3 [ ( +ɛ Q 3r ( 448s + 45s 5 ( +Q 3 6 ( 38s 35s 95r ( 48s +Q 7 43r 7 536s 43r 5 ( + Q 9s r + 9rs + Q 4 ( 4s 5r 3 s 5r 7r 4 4s 7r s + 85s ( + Q 87 45r 49 s s + 97s s 49r 4 ( 38s + 45s ( 95r 3 + Q 4 764s s 695r 4 836s s 695r ( ( 48 9s + Q 54s 5 637r 5 3 ( 9s 8s 637r 3 + Q ( 7s ( 4s s + Q 6 r 6 6s s r 4 4 7s ] 4r ] Due to the lac of time, I am unable to continue with the calculations. The Mathematica code that I use seems to have trouble computing the ψ s and ψ s components. This could be related
9 9 the 3D equivalent of the Stoe s paradox, or could simply be a bug in the code. In any case, it would be important to carry this calculation to completion at some point in the future in order to compare the displacement resulting from creeping flow to that in potential flow. 4 Discussion A calculation similar to the one done above was carried by Kelly for a two dimensional circle where the algebra is less involved. Figure 3 compares the results. In all cases the shape trajectory is a unit circle (s, s shape space. potential 3d potential d z stoes d.3.. s s Figure 3: Geometric phase for a particular gate. TODO: stoes 3d. Explicitly, we quote the results. For the the potential flow case, Kelly obtained: and for Stoe s case ż = ɛ (s ṡ s ṡ + O(ɛ 3 ż = ɛ ( 4 s ṡ + s ṡ + O(ɛ If we compare those, to our result
10 Figure 4: Ilustration of axisymmetroc squirming. ż = ɛ 3 ( 68 + ɛ 5ṡ 45 s ṡ 8 5 s ṡ ( The first thing that we notice is the presence of O(ɛ in the 3d case (our solution. However, this term has no effect on the net displacement and it only results in a more pronounced oscillation of the center of mass during each cycle. Also we notice that the body achieves higher displacementin 3d. This simplicity with which such explicit equations such as can be obtained is a testament to the power of the geometric approach to mechanics. References Bloch A, 7 Nonholonomic Mechanics and Control (Springer Holm D D, 8 GEOMETRIC MECHANICS: Dynamics and Symmetry illustrated edition edition (Imperial College Press Kelly S D, 998 The Mechanics and Control of Robotic Locomotion with Appliations to Aquatic Vehicles Ph.D. thesis California Institute of Technology Kohr M, Pop I, 4 Viscous Incompressible Flow for Low Reynolds Numbers (Advances in Boundary Elements illustrated edition edition (WIT Press (UK Leal L G, 7 Advanced Transport Phenomena: Fluid Mechanics and Convective Transport Processes st edition (Cambridge University Press Marsden J E, Ratiu T S, Introduction to Mechanics and Symmetry: A Basic Exposition of Classical Mechanical Systems nd edition (Springer Mason R, Burdic J, 999 Propulsion and control of deformable bodies in an ideal fluid in 999 IEEE International Conference on Robotics and Automation, 999. Proceedings, volume
11 Melli J B, Rowley C W, Rufat D S, 6 Motion planning for an articulated body in a perfect planar fluid SIAM Journal on Applied Dynamical Systems 5 65 Nair S, Kanso E, 7 Hydrodynamically coupled rigid bodies Journal of Fluid Mechanics Rufat D S, 5 Locomotion of a Three Dimensional Variable Shape Body in Inviscid Fluid Undergraduate independent wor Princeton University
Chapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationThe Squirmer model and the Boundary Element Method
Hauptseminar : Active Matter The Squirmer model and the Boundary Element Method Miru Lee April 6, 017 Physics epartment, University of Stuttgart 1 Introduction A micro scale swimmer exhibits directional
More informationHomework Two. Abstract: Liu. Solutions for Homework Problems Two: (50 points total). Collected by Junyu
Homework Two Abstract: Liu. Solutions for Homework Problems Two: (50 points total). Collected by Junyu Contents 1 BT Problem 13.15 (8 points) (by Nick Hunter-Jones) 1 2 BT Problem 14.2 (12 points: 3+3+3+3)
More informationLocomotion of Articulated Bodies in a Perfect Fluid
J. Nonlinear Sci. Vol. 15: pp. 255 289 (2005) DOI: 10.1007/s00332-004-0650-9 2005 Springer Science+Business Media, Inc. Locomotion of Articulated Bodies in a Perfect Fluid E. Kanso, 1 J. E. Marsden, 1
More informationPassive locomotion via normal-mode coupling in a submerged spring mass system
J. Fluid Mech. (9), vol. 64, pp. 5 5. c Cambridge University Press 9 doi:.7/s999357 5 Passive locomotion via normal-mode coupling in a submerged spring mass system EVA KANSO AND PAUL K NEWTON, Department
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationReview of Lagrangian Mechanics and Reduction
Review of Lagrangian Mechanics and Reduction Joel W. Burdick and Patricio Vela California Institute of Technology Mechanical Engineering, BioEngineering Pasadena, CA 91125, USA Verona Short Course, August
More informationPassive locomotion via normal mode. coupling in a submerged spring-mass system
Under consideration for publication in J. Fluid Mech. Passive locomotion via normal mode coupling in a submerged spring-mass system By E V A K A N S O A N D P A U L K N E W T O N,2 Department of Aerospace
More informationControlling Mechanical Systems by Active Constraints. Alberto Bressan. Department of Mathematics, Penn State University
Controlling Mechanical Systems by Active Constraints Alberto Bressan Department of Mathematics, Penn State University 1 Control of Mechanical Systems: Two approaches by applying external forces by directly
More informationi.e. the conservation of mass, the conservation of linear momentum, the conservation of energy.
04/04/2017 LECTURE 33 Geometric Interpretation of Stream Function: In the last class, you came to know about the different types of boundary conditions that needs to be applied to solve the governing equations
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationAvailable online at ScienceDirect. IFAC-PapersOnLine (2016)
Available online at www.sciencedirect.com ScienceDirect IFAC-PapersOnLine 49-18 (216) 988 993 Geometric Controllability of The Purcell s Swimmer and its Symmetrized Cousin Sudin Kadam Ravi N. Banavar Systems
More informationGlobal Formulations of Lagrangian and Hamiltonian Dynamics on Embedded Manifolds
1 Global Formulations of Lagrangian and Hamiltonian Dynamics on Embedded Manifolds By Taeyoung Lee, Melvin Leok, and N. Harris McClamroch Mechanical and Aerospace Engineering, George Washington University,
More informationM3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011
M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011 Text for the course: Professor Darryl D Holm 25 October 2011 Imperial College London d.holm@ic.ac.uk http://www.ma.ic.ac.uk/~dholm/ Geometric Mechanics
More informationGFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability
GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What
More informationOffshore Hydromechanics Module 1
Offshore Hydromechanics Module 1 Dr. ir. Pepijn de Jong 4. Potential Flows part 2 Introduction Topics of Module 1 Problems of interest Chapter 1 Hydrostatics Chapter 2 Floating stability Chapter 2 Constant
More informationSteering the Chaplygin Sleigh by a Moving Mass
Steering the Chaplygin Sleigh by a Moving Mass Jason M. Osborne Department of Mathematics North Carolina State University Raleigh, NC 27695 jmosborn@unity.ncsu.edu Dmitry V. Zenkov Department of Mathematics
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
More informationFluid dynamics is the study of flows of fluid. However, before we begin discussing the dynamics of
Hydrodynamics of Newtonian Fluids Cameron Stewart Universität Stuttgart cameron.n.stewart@gmail.com April 5, 2017 Abstract In this review we provide an overview of the hydrodynamics of viscous, Newtonian
More informationUse conserved quantities to reduce number of variables and the equation of motion (EOM)
Physics 106a, Caltech 5 October, 018 Lecture 8: Central Forces Bound States Today we discuss the Kepler problem of the orbital motion of planets and other objects in the gravitational field of the sun.
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationElectromagnetism: Worked Examples. University of Oxford Second Year, Part A2
Electromagnetism: Worked Examples University of Oxford Second Year, Part A2 Caroline Terquem Department of Physics caroline.terquem@physics.ox.ac.uk Michaelmas Term 2017 2 Contents 1 Potentials 5 1.1 Potential
More informationLocomotion of Articulated Bodies in a Perfect Fluid
Locomotion of Articulated Bodies in a Perfect Fluid Eva Kanso, 1 Jerrold E. Marsden, 1 Clarence W. Rowley 2 and Juan Melli-Huber 2 1 Control and Dynamical Systems 2 Mechanical and Aerospace Engineering
More information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More informationCOMPLETE ALL ROUGH WORKINGS IN THE ANSWER BOOK AND CROSS THROUGH ANY WORK WHICH IS NOT TO BE ASSESSED.
BSc/MSci EXAMINATION PHY-304 Time Allowed: Physical Dynamics 2 hours 30 minutes Date: 28 th May 2009 Time: 10:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from section
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationAn Introduction to Differential Flatness. Will Burgoyne May 7, 2007 ME598 Geometric Mechanics
An Introduction to Differential Flatness Will Burgoyne May 7, 2007 ME598 Geometric Mechanics Definitions Underdetermined System a system of m ODEs with n dependent variables where m
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationCDS 205 Final Project: Incorporating Nonholonomic Constraints in Basic Geometric Mechanics Concepts
CDS 205 Final Project: Incorporating Nonholonomic Constraints in Basic Geometric Mechanics Concepts Michael Wolf wolf@caltech.edu 6 June 2005 1 Introduction 1.1 Motivation While most mechanics and dynamics
More informationMath 575-Lecture 19. In this lecture, we continue to investigate the solutions of the Stokes equations.
Math 575-Lecture 9 In this lecture, we continue to investigate the solutions of the Stokes equations. Energy balance Rewrite the equation to σ = f. We begin the energy estimate by dotting u in the Stokes
More informationTheory of Vibrations in Stewart Platforms
Theory of Vibrations in Stewart Platforms J.M. Selig and X. Ding School of Computing, Info. Sys. & Maths. South Bank University London SE1 0AA, U.K. (seligjm@sbu.ac.uk) Abstract This article develops a
More informationAutonomous Underwater Vehicles: Equations of Motion
Autonomous Underwater Vehicles: Equations of Motion Monique Chyba - November 18, 2015 Departments of Mathematics, University of Hawai i at Mānoa Elective in Robotics 2015/2016 - Control of Unmanned Vehicles
More informationClassical Mechanics III (8.09) Fall 2014 Assignment 3
Classical Mechanics III (8.09) Fall 2014 Assignment 3 Massachusetts Institute of Technology Physics Department Due September 29, 2014 September 22, 2014 6:00pm Announcements This week we continue our discussion
More informationGeneral Solution of the Incompressible, Potential Flow Equations
CHAPTER 3 General Solution of the Incompressible, Potential Flow Equations Developing the basic methodology for obtaining the elementary solutions to potential flow problem. Linear nature of the potential
More informationWKB Approximation in 3D
1 WKB Approximation in 3D We see solutions ψr of the stationary Schrodinger equations for a spinless particle of energy E: 2 2m 2 ψ + V rψ = Eψ At rst, we just rewrite the Schrodinger equation in the following
More informationp + µ 2 u =0= σ (1) and
Master SdI mention M2FA, Fluid Mechanics program Hydrodynamics. P.-Y. Lagrée and S. Zaleski. Test December 4, 2009 All documentation is authorized except for reference [1]. Internet access is not allowed.
More informationElectric fields in matter
Electric fields in matter November 2, 25 Suppose we apply a constant electric field to a block of material. Then the charges that make up the matter are no longer in equilibrium: the electrons tend to
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationComplex functions in the theory of 2D flow
Complex functions in the theory of D flow Martin Scholtz Institute of Theoretical Physics Charles University in Prague scholtz@utf.mff.cuni.cz Faculty of Transportation Sciences Czech Technical University
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationDynamics and Stability application to submerged bodies, vortex streets and vortex-body systems
Dynamics and Stability application to submerged bodies, vortex streets and vortex-body systems Eva Kanso University of Southern California CDS 140B Introduction to Dynamics February 5 and 7, 2008 Fish
More informationQUALIFYING EXAMINATION, Part 1. Solutions. Problem 1: Mathematical Methods. r r r 2 r r2 = 0 r 2. d 3 r. t 0 e t dt = e t
QUALIFYING EXAMINATION, Part 1 Solutions Problem 1: Mathematical Methods (a) For r > we find 2 ( 1 r ) = 1 ( ) 1 r 2 r r2 = 1 ( 1 ) r r r 2 r r2 = r 2 However for r = we get 1 because of the factor in
More informationHamiltonian aspects of fluid dynamics
Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline
More informationFaraday Instability on Elastic Interfaces
Faraday Instability on Elastic Interfaces A Project Report Submitted to the Faculty of the WORCESTER POLYTECHNIC INSTITUTE In partial fulfillment of the requirements for the Degree of Master of Science
More informationAssignment 2. Goldstein 2.3 Prove that the shortest distance between two points in space is a straight line.
Assignment Goldstein.3 Prove that the shortest distance between two points in space is a straight line. The distance between two points is given by the integral of the infinitesimal arclength: s = = =
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More information1. Comparison of stability analysis to previous work
. Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing
More informationParticle self-diffusiophoresis near solid walls and interfaces
Particle self-diffusiophoresis near solid walls and interfaces Darren G. CROWDY Tel.: +44(0)2075948587; Fax: +44(0)2075948517; Email: d.crowdy@imperial.ac.uk Department of Mathematics, Imperial College
More informationApplications in Fluid Mechanics
CHAPTER 8 Applications in Fluid 8.1 INTRODUCTION The general topic of fluid mechanics encompasses a wide range of problems of interest in engineering applications. The most basic definition of a fluid
More informationThe Kelvin transformation in potential theory and Stokes flow
IMA Journal of Applied Mathematics (2009 74, 427 438 doi:10.1093/imamat/hxn027 Advance Access publication on September 9, 2008 The Kelvin transformation in potential theory and Stokes flow GEORGE DASSIOS
More informationContents. Dynamics and control of mechanical systems. Focus on
Dynamics and control of mechanical systems Date Day 1 (01/08) Day 2 (03/08) Day 3 (05/08) Day 4 (07/08) Day 5 (09/08) Day 6 (11/08) Content Review of the basics of mechanics. Kinematics of rigid bodies
More informationPart II. Fluid Dynamics II. Year
Part II Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 Paper 2, Section II 38C 41 An initially unperturbed two-dimensional inviscid jet in h < y < h has uniform speed U
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More informationF11AE1 1. C = ρν r r. r u z r
F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationSolution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3
: Solution Set Two Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein October 7, 2015 Contents 1 Problem #1: Projectile Motion. 2 1.1 Cartesian Coordinates....................................
More information7 Curvilinear coordinates
7 Curvilinear coordinates Read: Boas sec. 5.4, 0.8, 0.9. 7. Review of spherical and cylindrical coords. First I ll review spherical and cylindrical coordinate systems so you can have them in mind when
More informationLecture 17 - Gyroscopes
Lecture 17 - Gyroscopes A Puzzle... We have seen throughout class that the center of mass is a very powerful tool for evaluating systems. However, don t let yourself get carried away with how useful it
More informationProblem 1. Mathematics of rotations
Problem 1. Mathematics of rotations (a) Show by algebraic means (i.e. no pictures) that the relationship between ω and is: φ, ψ, θ Feel free to use computer algebra. ω X = φ sin θ sin ψ + θ cos ψ (1) ω
More informationA Lagrangian approach to the kinematic dynamo
1 A Lagrangian approach to the kinematic dynamo Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc/ 5 March 2001 with Allen
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More informationA note on interior vs. boundary-layer damping of surface waves in a circular cylinder
J. Fluid Mech. (1998), vol. 364, pp. 319 323. Printed in the United Kingdom c 1998 Cambridge University Press 319 A note on interior vs. boundary-layer damping of surface waves in a circular cylinder By
More informationDIFFERENTIAL GEOMETRY HW 12
DIFFERENTIAL GEOMETRY HW 1 CLAY SHONKWILER 3 Find the Lie algebra so(n) of the special orthogonal group SO(n), and the explicit formula for the Lie bracket there. Proof. Since SO(n) is a subgroup of GL(n),
More informationAE301 Aerodynamics I UNIT B: Theory of Aerodynamics
AE301 Aerodynamics I UNIT B: Theory of Aerodynamics ROAD MAP... B-1: Mathematics for Aerodynamics B-: Flow Field Representations B-3: Potential Flow Analysis B-4: Applications of Potential Flow Analysis
More information3.5 Vorticity Equation
.0 - Marine Hydrodynamics, Spring 005 Lecture 9.0 - Marine Hydrodynamics Lecture 9 Lecture 9 is structured as follows: In paragraph 3.5 we return to the full Navier-Stokes equations (unsteady, viscous
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
More informationGeometric Control of Particle Manipulation in a Two-Dimensional Fluid
Submitted, 2009 Conference on Decision and Control (CDC) http://www.cds.caltech.edu/~murray/papers/or+09-cdc.html Geometric Control of Particle Manipulation in a Two-Dimensional Fluid Yizhar Or, Joris
More informationGeometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics
Geometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics Harris McClamroch Aerospace Engineering, University of Michigan Joint work with Taeyoung Lee (George Washington University) Melvin
More informationLecture 1: Introduction to Linear and Non-Linear Waves
Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types
More informationHyperviscous regularization of the Navier Stokes equation and the motion of slender swimmers
Hyperviscous regularization of the Navier Stokes equation and the motion of slender swimmers Giulio G. Giusteri Dipartimento di Matematica e Fisica N. Tartaglia Università Cattolica del Sacro Cuore International
More informationPhys 7221 Homework # 8
Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with
More informationÜbungen zu RT2 SS (4) Show that (any) contraction of a (p, q) - tensor results in a (p 1, q 1) - tensor.
Übungen zu RT2 SS 2010 (1) Show that the tensor field g µν (x) = η µν is invariant under Poincaré transformations, i.e. x µ x µ = L µ νx ν + c µ, where L µ ν is a constant matrix subject to L µ ρl ν ση
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationAnalysis of a Hinge-Connected Flapping Plate with an Implemented Torsional Spring Model
Analysis of a Hinge-Connected Flapping Plate with an Implemented Torsional Spring Model Zach Gaston 1, Hui Wan 2 and Haibo Dong 3 Department of Mechanical & Materials Engineering, Wright State University,
More informationROLLING OF A SYMMETRIC SPHERE ON A HORIZONTAL PLANE WITHOUT SLIDING OR SPINNING. Hernán Cendra 1. and MARÍA ETCHECHOURY
Vol. 57 (2006) REPORTS ON MATHEMATICAL PHYSICS No. 3 ROLLING OF A SYMMETRIC SPHERE ON A HORIZONTAL PLANE WITHOUT SLIDING OR SPINNING Hernán Cendra 1 Universidad Nacional del Sur, Av. Alem 1253, 8000 Bahía
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationThe Navier-Stokes equations for an incompressible fluid in indicial notation are. i u i = 0 (9.1)
Chapter 9 Viscous flow The Navier-Stokes equations for an incompressible fluid in indicial notation are i u i = 0 9.1 ρ t + u j j u i = i p + µ 2 j u i 9.2 When these are applied to practical problems,
More informationarxiv: v1 [cs.sy] 21 Mar 2018
Modelling and controllability of the motion of a slender, flexible micro-swimmer Sudin Kadam and Ravi N. Banavar arxiv:83.8795v [cs.sy] 2 Mar 28 Abstract The mechanism of swimming at very low Reynolds
More informationThe Astrojax Pendulum and the N-Body Problem on the Sphere: A study in reduction, variational integration, and pattern evocation.
The Astrojax Pendulum and the N-Body Problem on the Sphere: A study in reduction, variational integration, and pattern evocation. Philip Du Toit 1 June 005 Abstract We study the Astrojax Pendulum and the
More information06. Lagrangian Mechanics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 06. Lagrangian Mechanics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationBefore you begin read these instructions carefully.
MATHEMATICAL TRIPOS Part IB Friday, 7 June, 2013 1:30 pm to 4:30 pm PAPER 4 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question
More informationarxiv: v2 [math-ph] 14 Apr 2008
Exact Solution for the Stokes Problem of an Infinite Cylinder in a Fluid with Harmonic Boundary Conditions at Infinity Andreas N. Vollmayr, Jan-Moritz P. Franosch, and J. Leo van Hemmen arxiv:84.23v2 math-ph]
More informationScattering. March 20, 2016
Scattering March 0, 06 The scattering of waves of any kind, by a compact object, has applications on all scales, from the scattering of light from the early universe by intervening galaxies, to the scattering
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationarxiv: v1 [physics.flu-dyn] 1 Feb 2015
Stokesian swimmers and active particles B. U. Felderhof Institut für Theorie der Statistischen Physik RWTH Aachen University Templergraben 55 52056 Aachen Germany arxiv:1502.00242v1 [physics.flu-dyn] 1
More informationNoether s Theorem. 4.1 Ignorable Coordinates
4 Noether s Theorem 4.1 Ignorable Coordinates A central recurring theme in mathematical physics is the connection between symmetries and conservation laws, in particular the connection between the symmetries
More informationExercise 9, Ex. 6.3 ( submarine )
Exercise 9, Ex. 6.3 ( submarine The flow around a submarine moving at at velocity V can be described by the flow caused by a source and a sink with strength Q at a distance a from each other. V x Submarine
More informationHigher Orders Instability of a Hollow Jet Endowed with Surface Tension
Mechanics and Mechanical Engineering Vol. 2, No. (2008) 69 78 c Technical University of Lodz Higher Orders Instability of a Hollow Jet Endowed with Surface Tension Ahmed E. Radwan Mathematics Department,
More informationQuestion 1: Spherical Pendulum
Question 1: Spherical Pendulum Consider a two-dimensional pendulum of length l with mass M at its end. It is easiest to use spherical coordinates centered at the pivot since the magnitude of the position
More informationGeometric Control of Particle Manipulation in a Two-Dimensional Fluid
To appear, 2009 Conference on Decision and Control (CDC) http://www.cds.caltech.edu/~murray/papers/or+09-cdc.html Geometric Control of Particle Manipulation in a Two-Dimensional Fluid Yizhar Or, Joris
More informationOn mechanical control systems with nonholonomic constraints and symmetries
ICRA 2002, To appear On mechanical control systems with nonholonomic constraints and symmetries Francesco Bullo Coordinated Science Laboratory University of Illinois at Urbana-Champaign 1308 W. Main St,
More informationMassachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004
Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve
More informationNonconservativity and noncommutativity in locomotion
Eur. Phys. J. Special Topics 224, 3141 3174 (215) c EDP Sciences, Springer-Verlag 215 DOI: 1.114/epjst/e215-585-y THE EUROPEAN PHYSICAL JOURNAL SPECIAL TOPICS Regular Article Nonconservativity and noncommutativity
More informationWeek 7: Integration: Special Coordinates
Week 7: Integration: Special Coordinates Introduction Many problems naturally involve symmetry. One should exploit it where possible and this often means using coordinate systems other than Cartesian coordinates.
More information