DOCSLIB.ORG

4 Rotation Coordinate Systems and the Equations of Motion

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
4 Rotation Coordinate Systems and the Equations of Motion EOSC 512 2019 4 Rotation Coordinate Systems and the Equations of Motion 4.1 Rates of change revisited We have now derived the Navier-Stokes equations in an inertial (non-accelerating) frame of reference for which Newton’s third law is valid. However, in oceanography and meteorology it is more natural to put ourselves in an earth-fixed coordinate frame that is rotating with the planet (Only natural – this is the frame from which we observe motions in the ocean and the atmosphere). However, rotating reference frames are not inertial. It is therefore nexessary to examine how the equations of motion must be altered to take this into account. Consider two observers. one fixed in inertial space and one fixed in our rotating frame. Consider how the rate of change of a vector fixed in length will appear to both observers if the vector rotates with the rotation of the Earth. Where ⌦is the angular velocity of the rotating frame. The earth has a rotational period of one day, so ⌦= 5 1 7.29 10− s− . ⇥ This introduces three interesting questions: ˆ Dˆi 1. What rate of change of i ( Dt ) is seen by the inertial observer? By the rotating observer? 2. Now let ˆi be an arbitrary vector A~ = Ajˆij. How are the rates of change of this vector seen by the inertial and DA~ DA~ rotating observer related? That is, how does Dt relate to Dt ? inertial rotating ⇣ ⌘ ⇣ ⌘ ~ 3. Now choose A to be the velocity seen in the inertial frame (~u )inertial. How is the rate of change of the inertial velocity as seen in inertial space related to the acceleration as seen in the rotating coordinate frame? We will answer these questions indiviually below: 1. Consider the motion of the unit vector ˆi in a time ∆t. It will swing around the rotation axis under the unfluence of the rotation, moving at right andles to itself and to the rotation vector as shown: Page 36 EOSC 512 2019 The di↵erence between the two vectors, ∆~i can be approximated – for small ∆t – as the chord generated in that time between the perpendicular distance from the rotation axis to the tip of the unit vector at t = t0 & the perpendicular distance from the rotation axis to the tip of the unit vector at t = t0 +∆t. Top-down view: Here, the radius of the circle traced out by the tip of the unit vector is r =sinφ ˆi and the angle of swing is the product of the angular velocity of rotation and the “swinging time”, ✓ =⌦∆t. Thus, the chord length is r✓ = ˆi sin φ⌦t. ~ The direction of ∆i is perpendicular to both the unit vector and the rotation vector since the change can only alter its direction and not its length. Thus, the change in the unit vector in time ∆t is: ~ ⌦~ ˆi ∆ˆi = ˆi sin φ⌦t ⇥ (4.1) ⌦~ ˆi ⇥ Where ˆi sin φ⌦= ⌦~ ˆi by the definition of the cross product. This leaves ∆~ˆi =∆t⌦ ˆi.Dividingby∆t ⇥ ⇥ and taking the limit ∆t 0, ! Dˆi = ⌦~ ˆi (4.2) Dt ⇥ !inertial This is the rate of change of the unit vector ˆi seen by the inertial observer! Notice that if the observer is rotating with ⌦, the see no change at all, so: Dˆi = 0 (4.3) Dt !rotating 2. Let A~ be represented as three components along the coordinate axes with unit vectors ˆi1, ˆi2, and ˆi3. Assuming the coordinate axes are still rotating with angular velocity ⌦. The individual component of A~ along each coordinate axis is the shadow of the vector cast along that axis. It Page 37 EOSC 512 2019 is a scalar whose value and rate of change is seen the same by both the inertial and rotating observers. The inertial observer also sees the rate of change of the unit vectors. Thus, DA~ DA Dˆi = i ˆi + A Dt Dt i i Dt !inertial DAi (4.4) = ˆi + A ⌦ ˆi Dt i i ⇥ DA ⇣ ⌘ = i ˆi +⌦ A~ Dt i ⇥ DAi ˆ ~ Physically, Dt ii is the rate of change of A seen by the rotating obeserver. He/she is only aware of the increase or decrease of the individual scalar coordinates along the coordinate unit vectors that the observer sees as constant. This can be interpreted as a stretching of the components of A~.⌦ A~ is an additional component ⇥ of the rate of change of A~ seen only by the inertial observer. This is the swinging of A~ due to the rotation. Together, DA~ DA~ = +⌦ A~ (4.5) Dt Dt ⇥ !inertial !rotating 3. Taking the previous equation and using ~x as the position of the fluid element, D~x DAx~ = +⌦ ~x (4.6) Dt Dt ⇥ ✓ ◆inertial !rotating The derivative of position with respect to time is velocity, so: ~u = ~u +⌦ ~x (4.7) inertial rotating ⇥ Similarly, if we choose A~ to be the velocity seen in the inertial frame ~u inertial, D~u D~u inertial = inertial + ⌦~ ~u Dt Dt ⇥ inertial ✓ ◆inertial ✓ ◆rotating D = ~u rotating + ⌦~ ~x + ⌦~ ~u rotating + ⌦~ ~x Dt ⇥ rotating ⇥ ⇥ ✓ ⇣ ⌘◆ ⇣ ⌘ D~u rotating D~x = + ⌦~ + ⌦~ ~u rotating + ⌦~ ⌦~ ~x (4.8) Dt rotating ⇥ Dt rotating ⇥ ⇥ ⇥ ✓ ◆ ✓ ◆ ⇣ ⌘ D~u inertial D~u rotating = +2⌦~ ~u rotating + ⌦~ ⌦~ ~x Dt inertial Dt rotating ⇥ ⇥ ⇥ ✓ ◆ ✓ ◆ ⇣ ⌘ 2 3 1 | {z } | {z } | {z } 1 The acceleration that an observer in the rotating frame would see. ⌘ 2 The coriolis acceleration. ⌘ 3 The (familiar?) centripetal accerlation. This can be written as the gradient of a scalar, ⌦~ ⌦~ ~x = ⌘ ⇥ ⇥ 1 (⌦ ~x )2. This is the acceleration an object feels resulting from the change of direction (but⇣ not⌘ the − 2 r ⇥ change in speed!) associated with uniform circular motion. Has magnitude ⌦ 2r where r is the distance | | between the rotation axis and the tip of the vector. It is directed toward the centre of the circle (hence centripetal!). Page 38 EOSC 512 2019 4.2 The Coriolis acceleration/force We have just seen that the acceleration in the inertial frame is related to the acceleration that an observer would see in the rotating frame with two additional contributions: The centripetal accerlation (acceleration directed towards the centre of the circle traced out by the tip of the rotating vector) & the Coriolis acceleration. The so-called Coriolis acceleration acts perpendicular to the velocity as seen in the rotating frame. It represents the acceleration due to the swinging of the velocity vector by the rotation vector. A factor of two enters because the rotation vector swings the velocity ~u around, but also the velocity ⌦~ ~x which is not seen in the rotating rotating ⇥ frame. This term will increase if the position vector ~x increases, giving rise to a second factor of ⌦~ ~u in the ⇥ rotating acceleration. As we will see, the Coriolis acceleration is a dominant (lead order) term in the dynamics of large-scale, low frequency motion in the atmospheric and oceans. This includes general banded atmospheric circulations such as trade winds and westerlies and large-scale ocean gyre circulations. The Coriolis acceleration, like all accelerations, is produced by a force in the direction of the acceleration. First we assume we have uniform velocity seen in the rotating system (which you will note is counterclockwise rotation). This implies that there will be no acceleration perceived in the system. In order for the flow to continue on a straight path as in the figure above, a force must be applied from the right. If that force is removed, the fluid element can no longer continue in a straight line and will veer to the right. This gives rise to the sensation that a “force” is acting to push the fluid to the right. Of course, since there are no forces in the inertial frame, then there are no acceleration in the inertial frame. For this reason, both the centripedal force and the Coriolis force are sometimes called pseudoforces as they are experienced only in the frame of reference of the body undergoing circular motion. It should be emphasized that there is nothing fictitious about these forces or accelerations – they are real – they are just not recognized as acceleration by an observer in the rotating frame. An alternative, more heurisitic way to think of the Coriolis e↵ect is the tendency for objects in motion (viewed in a rotating frame of reference) to follow curved paths. They deflect to the right for counter-clockwise rotation. Consider an observer on Earth watching the motion of a fluid element in the atmosphere. Because the Earth is rotating, points on the Earth’s surface travel eastward with a speed given by the Earth’s Page 39 EOSC 512 2019 angular velocity. Further, because the Earth is spherical, the distance a point travels in one day at low latitudes is bigger compared to higher latitudes. As a results, the eastward velocity of a point on the Earth’s surface varies with latitude: Fast at the equator and slower as you move towards the poles. Example 4.1. If a fluid parcel at the equator has a northward inertial velocity and no other forces acting, it also has an eastward component of velocity of approximately 444 m/s. To a rotating observer at 49°, it also has an eastward component of velocity due to di↵erence in the angular velocity between the equator and 49°. To the observer at 49°, in some time ∆t, the parcel has moved north, but also eastwards by an amount (444 − Page 40 EOSC 512 2019 300) m/s∆t. That is, the observer at 49° sees the parcel travel a curved path that is deflected to the right of its initial trajectory northwards. We can heuristically think of this as arising because the parcel is moving eastward faster than the ground over which it is flying as it moves northward (in a sense, it gets ahead of the planet).
Load more

Recommended publications
  • Multidisciplinary Design Project Engineering Dictionary Version 0.0.2
    Multidisciplinary Design Project Engineering Dictionary Version 0.0.2 February 15, 2006 . DRAFT Cambridge-MIT Institute Multidisciplinary Design Project This Dictionary/Glossary of Engineering terms has been compiled to compliment the work developed as part of the Multi-disciplinary Design Project (MDP), which is a programme to develop teaching material and kits to aid the running of mechtronics projects in Universities and Schools. The project is being carried out with support from the Cambridge-MIT Institute undergraduate teaching programe. For more information about the project please visit the MDP website at http://www-mdp.eng.cam.ac.uk or contact Dr. Peter Long Prof. Alex Slocum Cambridge University Engineering Department Massachusetts Institute of Technology Trumpington Street, 77 Massachusetts Ave. Cambridge. Cambridge MA 02139-4307 CB2 1PZ. USA e-mail: [email protected] e-mail: [email protected] tel: +44 (0) 1223 332779 tel: +1 617 253 0012 For information about the CMI initiative please see Cambridge-MIT Institute website :- http://www.cambridge-mit.org CMI CMI, University of Cambridge Massachusetts Institute of Technology 10 Miller’s Yard, 77 Massachusetts Ave. Mill Lane, Cambridge MA 02139-4307 Cambridge. CB2 1RQ. USA tel: +44 (0) 1223 327207 tel. +1 617 253 7732 fax: +44 (0) 1223 765891 fax. +1 617 258 8539 . DRAFT 2 CMI-MDP Programme 1 Introduction This dictionary/glossary has not been developed as a definative work but as a useful reference book for engi- neering students to search when looking for the meaning of a word/phrase. It has been compiled from a number of existing glossaries together with a number of local additions.
    [Show full text]
  • Newtonian Mechanics Is Most Straightforward in Its Formulation and Is Based on Newton’S Second Law
    CLASSICAL MECHANICS D. A. Garanin September 30, 2015 1 Introduction Mechanics is part of physics studying motion of material bodies or conditions of their equilibrium. The latter is the subject of statics that is important in engineering. General properties of motion of bodies regardless of the source of motion (in particular, the role of constraints) belong to kinematics. Finally, motion caused by forces or interactions is the subject of dynamics, the biggest and most important part of mechanics. Concerning systems studied, mechanics can be divided into mechanics of material points, mechanics of rigid bodies, mechanics of elastic bodies, and mechanics of fluids: hydro- and aerodynamics. At the core of each of these areas of mechanics is the equation of motion, Newton's second law. Mechanics of material points is described by ordinary differential equations (ODE). One can distinguish between mechanics of one or few bodies and mechanics of many-body systems. Mechanics of rigid bodies is also described by ordinary differential equations, including positions and velocities of their centers and the angles defining their orientation. Mechanics of elastic bodies and fluids (that is, mechanics of continuum) is more compli- cated and described by partial differential equation. In many cases mechanics of continuum is coupled to thermodynamics, especially in aerodynamics. The subject of this course are systems described by ODE, including particles and rigid bodies. There are two limitations on classical mechanics. First, speeds of the objects should be much smaller than the speed of light, v c, otherwise it becomes relativistic mechanics. Second, the bodies should have a sufficiently large mass and/or kinetic energy.
    [Show full text]
  • Chapter 04 Rotational Motion
    Chapter 04 Rotational Motion P. J. Grandinetti Chem. 4300 P. J. Grandinetti Chapter 04: Rotational Motion Angular Momentum Angular momentum of particle with respect to origin, O, is given by l⃗ = ⃗r × p⃗ Rate of change of angular momentum is given z by cross product of ⃗r with applied force. p m dl⃗ dp⃗ = ⃗r × = ⃗r × F⃗ = ⃗휏 r dt dt O y Cross product is defined as applied torque, ⃗휏. x Unlike linear momentum, angular momentum depends on origin choice. P. J. Grandinetti Chapter 04: Rotational Motion Conservation of Angular Momentum Consider system of N Particles z m5 m 2 Rate of change of angular momentum is m3 ⃗ ∑N l⃗ ∑N ⃗ m1 dL d 훼 dp훼 = = ⃗r훼 × dt dt dt 훼=1 훼=1 y which becomes m4 x ⃗ ∑N dL ⃗ net = ⃗r훼 × F dt 훼 Total angular momentum is 훼=1 ∑N ∑N ⃗ ⃗ L = l훼 = ⃗r훼 × p⃗훼 훼=1 훼=1 P. J. Grandinetti Chapter 04: Rotational Motion Conservation of Angular Momentum ⃗ ∑N dL ⃗ net = ⃗r훼 × F dt 훼 훼=1 Taking an earlier expression for a system of particles from chapter 1 ∑N ⃗ net ⃗ ext ⃗ F훼 = F훼 + f훼훽 훽=1 훽≠훼 we obtain ⃗ ∑N ∑N ∑N dL ⃗ ext ⃗ = ⃗r훼 × F + ⃗r훼 × f훼훽 dt 훼 훼=1 훼=1 훽=1 훽≠훼 and then obtain 0 > ⃗ ∑N ∑N ∑N dL ⃗ ext ⃗ rd ⃗ ⃗ = ⃗r훼 × F + ⃗r훼 × f훼훽 double sum disappears from Newton’s 3 law (f = *f ) dt 훼 12 21 훼=1 훼=1 훽=1 훽≠훼 P.
    [Show full text]
  • Analytical Mechanics
    A Guided Tour of Analytical Mechanics with animations in MAPLE Rouben Rostamian Department of Mathematics and Statistics UMBC [email protected] December 2, 2018 ii Contents Preface vii 1 An introduction through examples 1 1.1 ThesimplependulumàlaNewton ...................... 1 1.2 ThesimplependulumàlaEuler ....................... 3 1.3 ThesimplependulumàlaLagrange.. .. .. ... .. .. ... .. .. .. 3 1.4 Thedoublependulum .............................. 4 Exercises .......................................... .. 6 2 Work and potential energy 9 Exercises .......................................... .. 12 3 A single particle in a conservative force field 13 3.1 The principle of conservation of energy . ..... 13 3.2 Thescalarcase ................................... 14 3.3 Stability....................................... 16 3.4 Thephaseportraitofasimplependulum . ... 16 Exercises .......................................... .. 17 4 TheKapitsa pendulum 19 4.1 Theinvertedpendulum ............................. 19 4.2 Averaging out the fast oscillations . ...... 19 4.3 Stabilityanalysis ............................... ... 22 Exercises .......................................... .. 23 5 Lagrangian mechanics 25 5.1 Newtonianmechanics .............................. 25 5.2 Holonomicconstraints............................ .. 26 5.3 Generalizedcoordinates .......................... ... 29 5.4 Virtual displacements, virtual work, and generalized force....... 30 5.5 External versus reaction forces . ..... 32 5.6 The equations of motion for a holonomic system . ...
    [Show full text]
  • Fundamental Governing Equations of Motion in Consistent Continuum Mechanics
    Fundamental governing equations of motion in consistent continuum mechanics Ali R. Hadjesfandiari, Gary F. Dargush Department of Mechanical and Aerospace Engineering University at Buffalo, The State University of New York, Buffalo, NY 14260 USA [email protected], [email protected] October 1, 2018 Abstract We investigate the consistency of the fundamental governing equations of motion in continuum mechanics. In the first step, we examine the governing equations for a system of particles, which can be considered as the discrete analog of the continuum. Based on Newton’s third law of action and reaction, there are two vectorial governing equations of motion for a system of particles, the force and moment equations. As is well known, these equations provide the governing equations of motion for infinitesimal elements of matter at each point, consisting of three force equations for translation, and three moment equations for rotation. We also examine the character of other first and second moment equations, which result in non-physical governing equations violating Newton’s third law of action and reaction. Finally, we derive the consistent governing equations of motion in continuum mechanics within the framework of couple stress theory. For completeness, the original couple stress theory and its evolution toward consistent couple stress theory are presented in true tensorial forms. Keywords: Governing equations of motion, Higher moment equations, Couple stress theory, Third order tensors, Newton’s third law of action and reaction 1 1. Introduction The governing equations of motion in continuum mechanics are based on the governing equations for systems of particles, in which the effect of internal forces are cancelled based on Newton’s third law of action and reaction.
    [Show full text]
  • Energy and Equations of Motion in a Tentative Theory of Gravity with a Privileged Reference Frame
    1 Energy and equations of motion in a tentative theory of gravity with a privileged reference frame Archives of Mechanics (Warszawa) 48, N°1, 25-52 (1996) M. ARMINJON (GRENOBLE) Abstract- Based on a tentative interpretation of gravity as a pressure force, a scalar theory of gravity was previously investigated. It assumes gravitational contraction (dilation) of space (time) standards. In the static case, the same Newton law as in special relativity was expressed in terms of these distorted local standards, and was found to imply geodesic motion. Here, the formulation of motion is reexamined in the most general situation. A consistent Newton law can still be defined, which accounts for the time variation of the space metric, but it is not compatible with geodesic motion for a time-dependent field. The energy of a test particle is defined: it is constant in the static case. Starting from ”dust‘, a balance equation is then derived for the energy of matter. If the Newton law is assumed, the field equation of the theory allows to rewrite this as a true conservation equation, including the gravitational energy. The latter contains a Newtonian term, plus the square of the relative rate of the local velocity of gravitation waves (or that of light), the velocity being expressed in terms of absolute standards. 1. Introduction AN ATTEMPT to deduce a consistent theory of gravity from the idea of a physically privileged reference frame or ”ether‘ was previously proposed [1-3]. This work is a further development which is likely to close the theory. It is well-known that the concept of ether has been abandoned at the beginning of this century.
    [Show full text]
  • Chapter 6 the Equations of Fluid Motion
    Chapter 6 The equations of fluid motion In order to proceed further with our discussion of the circulation of the at- mosphere, and later the ocean, we must develop some of the underlying theory governing the motion of a fluid on the spinning Earth. A differen- tially heated, stratified fluid on a rotating planet cannot move in arbitrary paths. Indeed, there are strong constraints on its motion imparted by the angular momentum of the spinning Earth. These constraints are profoundly important in shaping the pattern of atmosphere and ocean circulation and their ability to transport properties around the globe. The laws governing the evolution of both fluids are the same and so our theoretical discussion willnotbespecifictoeitheratmosphereorocean,butcanandwillbeapplied to both. Because the properties of rotating fluids are often counter-intuitive and sometimes difficult to grasp, alongside our theoretical development we will describe and carry out laboratory experiments with a tank of water on a rotating table (Fig.6.1). Many of the laboratory experiments we make use of are simplified versions of ‘classics’ that have become cornerstones of geo- physical fluid dynamics. They are listed in Appendix 13.4. Furthermore we have chosen relatively simple experiments that, in the main, do nor require sophisticated apparatus. We encourage you to ‘have a go’ or view the atten- dant movie loops that record the experiments carried out in preparation of our text. We now begin a more formal development of the equations that govern the evolution of a fluid. A brief summary of the associated mathematical concepts, definitions and notation we employ can be found in an Appendix 13.2.
    [Show full text]
  • Application of Newton's Second
    Chapter 8 Applications of Newton’s Second Law 8.1 Force Laws ............................................................................................................... 1 8.1.1 Hooke’s Law ..................................................................................................... 1 8.2.2 Principle of Equivalence: ................................................................................ 5 8.2.3 Gravitational Force near the Surface of the Earth ....................................... 5 8.2.4 Electric Charge and Coulomb’s Law ............................................................. 6 Example 8.1 Coulomb’s Law and the Universal Law of Gravitation .................. 7 8.3 Contact Forces ......................................................................................................... 7 8.3.1 Free-body Force Diagram ................................................................................... 9 Example 8.2 Tug-of-War .......................................................................................... 9 8.3.3 Normal Component of the Contact Force and Weight ............................... 11 8.3.4 Static and Kinetic Friction ............................................................................ 14 8.3.5 Modeling ......................................................................................................... 16 8.4 Tension in a Rope .................................................................................................. 16 8.4.1 Static Tension in a Rope ...............................................................................
    [Show full text]
  • Stress, Cauchy's Equation and the Navier-Stokes Equations
    Chapter 3 Stress, Cauchy’s equation and the Navier-Stokes equations 3.1 The concept of traction/stress • Consider the volume of fluid shown in the left half of Fig. 3.1. The volume of fluid is subjected to distributed external forces (e.g. shear stresses, pressures etc.). Let ∆F be the resultant force acting on a small surface element ∆S with outer unit normal n, then the traction vector t is defined as: ∆F t = lim (3.1) ∆S→0 ∆S ∆F n ∆F ∆ S ∆ S n Figure 3.1: Sketch illustrating traction and stress. • The right half of Fig. 3.1 illustrates the concept of an (internal) stress t which represents the traction exerted by one half of the fluid volume onto the other half across a ficticious cut (along a plane with outer unit normal n) through the volume. 3.2 The stress tensor • The stress vector t depends on the spatial position in the body and on the orientation of the plane (characterised by its outer unit normal n) along which the volume of fluid is cut: ti = τij nj , (3.2) where τij = τji is the symmetric stress tensor. • On an infinitesimal block of fluid whose faces are parallel to the axes, the component τij of the stress tensor represents the traction component in the positive i-direction on the face xj = const. whose outer normal points in the positive j-direction (see Fig. 3.2). 6 MATH35001 Viscous Fluid Flow: Stress, Cauchy’s equation and the Navier-Stokes equations 7 x3 x3 τ33 τ22 τ τ11 12 τ21 τ τ 13 23 τ τ 32τ 31 τ 31 32 τ τ τ 23 13 τ21 τ τ τ 11 12 22 33 x1 x2 x1 x2 Figure 3.2: Sketch illustrating the components of the stress tensor.
    [Show full text]
  • Lagrangian Mechanics - Wikipedia, the Free Encyclopedia Page 1 of 11
    Lagrangian mechanics - Wikipedia, the free encyclopedia Page 1 of 11 Lagrangian mechanics From Wikipedia, the free encyclopedia Lagrangian mechanics is a re-formulation of classical mechanics that combines Classical mechanics conservation of momentum with conservation of energy. It was introduced by the French mathematician Joseph-Louis Lagrange in 1788. Newton's Second Law In Lagrangian mechanics, the trajectory of a system of particles is derived by solving History of classical mechanics · the Lagrange equations in one of two forms, either the Lagrange equations of the Timeline of classical mechanics [1] first kind , which treat constraints explicitly as extra equations, often using Branches [2][3] Lagrange multipliers; or the Lagrange equations of the second kind , which Statics · Dynamics / Kinetics · Kinematics · [1] incorporate the constraints directly by judicious choice of generalized coordinates. Applied mechanics · Celestial mechanics · [4] The fundamental lemma of the calculus of variations shows that solving the Continuum mechanics · Lagrange equations is equivalent to finding the path for which the action functional is Statistical mechanics stationary, a quantity that is the integral of the Lagrangian over time. Formulations The use of generalized coordinates may considerably simplify a system's analysis. Newtonian mechanics (Vectorial For example, consider a small frictionless bead traveling in a groove. If one is tracking the bead as a particle, calculation of the motion of the bead using Newtonian mechanics) mechanics would require solving for the time-varying constraint force required to Analytical mechanics: keep the bead in the groove. For the same problem using Lagrangian mechanics, one Lagrangian mechanics looks at the path of the groove and chooses a set of independent generalized Hamiltonian mechanics coordinates that completely characterize the possible motion of the bead.
    [Show full text]
  • Model Comparison of DBD-PA-Induced Body Force in Quiescent Air and Separated Flow Over NACA0015
    Model comparison of DBD-PA-induced body force in quiescent air and separated flow over NACA0015 Di Chen, Kengo Asada, Satoshi Sekimoto, Kozo Fujii Tokyo University of Science, Katsushika, Tokyo 125-8585, Japan Hiroyuki Nishida Tokyo University of Agriculture and Technology, Koganei, Tokyo 184-8588, Japan Numerical simulations of plasma flows induced by dielectric barrier discharge plasma actuators (DBD-PA) are conducted with two different body-force models: Suzen-Huang (S- H) model and drift-diffusion (D-D) model. The induced flow generated in quiescent air over a flat plate in continuous actuation and the PA-based flow control effect with burst actuation in separated flow over NACA0015 are studied. In the comparative study, the body-force field and the induced velocity field are firstly investigated in the quiescent field to see the spatial difference and the temporal difference in a single discharge cycle. The D-D body force is computed with flush-mounted and bulge configuration of the exposed electrode, which is operated at the peak-to-peak AC voltage of 7kV and 10kV. The D-D models generate momentarily higher body force in the positive-going phase of the AC power, but activate smaller flow region than the S-H model with Dc = 0.0117, which is given by the experiment beforehand at 7kV.[21] The local induced velocity of the D-D bulge case at 7kV measured in the downstream flow has the best agreement with the experimental result.[36] The maximum wall-parallel induced velocity in the S-H case with Dc = 0.0117 is consistent with that in the experiment, however, the local induced velocity is relatively high with different flow structure.
    [Show full text]
  • Dynamic Meteorology - Introduction
    Dynamic Meteorology - Introduction Atmospheric dynamics – the study of atmospheric motions that are associated with weather and climate We will consider the atmosphere to be a continuous fluid medium, or continuum. Each “point” in the atmosphere will be made up of a large number of molecules, with certain properties. These properties are assumed to be continuous functions of position and time. The basic laws of thermodynamics and fluid mechanics can be expressed in terms of partial differential equations, with space and time as independent variables and the atmospheric properties as dependent variables. Physical Dimensions and Units Dimensional homogeneity – all terms in the equations that describe the atmosphere must have the same physical dimensions (units) The four base units we will use are: From these we will also use the following derived units: Fundamental Forces Newton’s Second Law: F=ma In atmospheric science it is typical to consider the force per unit mass acting on the atmosphere: Force = a mass In order to understand atmospheric motion (accelerations) we need to know what forces act on the atmosphere. € What are the fundamental forces of interest in atmospheric science? Body (or volume) force – a force that acts on the center of mass of a fluid parcel Surface force – a force that acts across the boundary separating a fluid parcel from its surroundings. The magnitude of surface forces are independent of the mass of the parcel. What are examples of body and surface forces? Pressure Gradient Force Example: Real-time weather map The force exerted on the left face of this air parcel due to pressure is: pA = pδyδz The force exerted on the right face of this air parcel due to pressure is: € % ∂p ( −(p + δp)δyδz = −' p + δx*δ yδz & ∂x ) The net force exerted by pressure on this air parcel is the sum of these forces and is equal to: € ∂p − δxδyδz ∂x By dividing by the mass of the air parcel (ρδxδyδz) we get the force per unit mass due to changes in pressure (i.e.
    [Show full text]