The proper definition of the added mass for the water entry problem

Size: px
Start display at page:

Download "The proper definition of the added mass for the water entry problem"

Transcription

1 The proper efinition of the ae mass for the water entry problem Leonaro Casetta Celso P. Pesce LIE&MO lui-structure Interaction an Offshore Mechanics Laboratory Mechanical Engineering Department Escola Politécnica University of São Paulo , São Paulo, Brazil Introuction The water entry problem is consiere from the point of view of analytic mechanics. It is usual practice to treat potential hyroynamic problems involving motion of soli boies within the frame of system ynamics. This is one whenever a finite number of generalize coorinates can be use as a proper representation for the motion of the whole flui. We refer to this kin of approach as hyro-mechanical. This is mae possible through the use of the well-known concept of ae mass. In the vertical water entry problem the generalize coorinate is the penetration epth of the impacting boy. rom another point of view, a usual way to calculate the impact force is the integration of the pressure-fiel over the wette surface of the impacting boy. The jets are then exclue from the portion of the flui that wets the boy. z or simplicity we take the case of vertical entry only. Uner the hyro-mechanical approach, one may consier only the of the flui to express the impact force as U ), (1) where M is the ae mass associate to the kinetic energy of the of the flui, jets being exclue. However, if the ae mass is interprete as a measure of the kinetic energy of the whole flui, another correct an alternative form for the impact force is given by 1 M U ) U, ()

2 where M refers to the whole flui omain. Shoul equivocal arguments be use one woul obtain the misleaing expression, 1 M U ) U. (3) The reason for the iscrepancy between expressions (1) an (3) has alreay been aresse an actually emerges from a consistent energy balance; see, e.g., the iscussion in Molin et al. (1996) or in Pesce (003, 006). In fact, via an asymptotic analysis near the contact line, it has been shown, for some particular cases, that a consierable amount 1 of kinetic energy is raine from the of the flui through the jets; Molin et al. (1996), Scolan an Korobkin (003), Casetta an Pesce (005). In other wors, the jets must be consistently consiere if energy arguments are applie in the impact force calculation. In the present work the role of the jets in the impact force calculation is iscusse. We show how one can obtain the correct impact force expression in either way, by consiering the jets, or by not consiering them. In this sense, proper efinitions of the ae mass for the problem are presente an then use in the corresponing proper forms of the Lagrange equation. Analysis consiering the of the flui omain A usual an practical approach in marine hyroynamics is to treat the water entry problem by taking the omain of analysis as the of the flui only. The whole flui omain is therefore ivie into two parts: the an the jets. As the boy penetrates the, jets are expelle from the neighborhoo of the contact line ue to a very rapi expansion of the boy wette surface. The 1 kinetic energy of the of the flui, T = M U, varies explicitly with the penetration epth of the boy, an so oes the ae mass associate to it. In this case, where an out-flux of kinetic energy oes exist from the omain uner analysis (the of flui), there is an effective loss of ae mass through the jets an one must use the extene Lagrange equation vali for systems with mass explicitly epenent on position - see Pesce (003) - that reas M z 1 U = U z. (4) The left-han sie of equation (4) must be interprete as the total force that acts on the of the flui. The first term is the impact force we want to calculate. The secon term, 1 ( ) M / z U, is the reactive force acting on the ue to the out-flux of kinetic energy. In fact, Casetta an Pesce (005) show, for the case of a generic an arbitrary shape of the contact line, that the rate of 1 / z U, being the reactive force kinetic energy raine by the jets is always equal to ( ) 3 term 1 ( / z U obtain the impact force in the form of equation (1), i.e., M 1 M ) also generally vali. Therefore, as T / z ( M /z) U =, one easily U ). (1a) 1 Actually this amount is exactly half the whole kinetic energy of the flui if the impact velocity is enforce to be constant; see Molin et al. (1996), for the case of cyliners; Casetta an Pesce (005), for the general 3D case. Such an extene form of the Lagrange equation shoul also contain a reactive force term ue to a real (physical) mass out-flux. In the water entry problem, the real mass out-flux through the jets is of secon orer, as shown by Cointe an Arman (1978) an Molin et al. (1996), for the particular an important case of a circular cyliner.

3 Analysis consiering the whole flui omain The water entry problem can also be treate by taking the whole flui omain. The ae mass 1 M, is now a measure of the kinetic energy of the whole flui omain, i.e., T = M U. In this case, the jets are inclue in the analysis omain an, obviously, there is no out-flux of kinetic energy neither an out-flux of ae mass. In other wors, there is no loss of energy from the system an this is the key point. The extene Lagrange equation, for systems with mass explicitly epenent on position is no longer applicable. One must then apply the usual form 3 of the Lagrange equation, as in Lamb (193), art. 136 an 137, i.e., = t U z. (5) The impact force expression is then obtaine from equation (5) in the form of equation (), 1 M U ) U. (a) The equivalence between impact force expressions By equating the alternative equations (1) an (), we promptly obtain M. (6) 1 ( M U ) U = ( M U ) Shoul a misleaing assumption be taken, by enforcing result woul be foun, i.e., M = M, an obviously erroneous 1 M U = 0, (7) such that 1 ( M / t) U woul be also null. This woul imply the amount of the kinetic energy in the jets to be null, Casetta an Pesce (005), an obviously false assertive. Energy balance rom equations (1) an () one can easily verify that an that T 1 U U = ( M ) 0, (8) T U = 0. (9) 3 Recall that the usual form of Lagrange equation is invariant with respect to systems with mass varying as a function of time - as is the case if the whole flui is taken as the omain; see, e.g., Pesce (003), for a etaile iscussion on this subject.

4 Both equivalent equations (8) an (9) represent the correct balance of energy between the flui an the impacting boy. The thir term in equation (8) is the time rate of the kinetic energy that fills in the jet. It is also interesting to mention that the general equation (8) can be alternatively obtaine via the classical velocity potential approach, by consiering the water entry problem as a nonlinear bounary-value problem; see Casetta an Pesce (005). Conclusion It was shown how to treat the water entry problem in two equivalent an consistent ways. If one consiers only the of the flui as the omain of analysis, such that there is an out-flux of kinetic energy to the jets, an if an analytic mechanics point of view is followe, the extene form of the Lagrange equation, for systems with mass explicitly epenent on position, must be applie to calculate the impact force. Alternatively, if the whole flui omain is taken, so that no longer exists an out-flux of kinetic energy, the usual form of the Lagrange equation must be applie. The consistency between both approaches lies on a proper omain efinition of the ae mass. Acknowlegments We acknowlege APESP, the State of São Paulo Research ounation, PhD scholarship, process n o 04/ an CNPq, The National Research Council, Research Grant n o 30450/00-5. The authors are especially grateful to Professor Alexaner Korobkin, for very interesting iscussions uring his visit to Brazil. References Casetta, L., Pesce, C.P., 005, A Noticeable Question of the Water Entry Problem: the Split of Kinetic Energy uring the Initial Stage, COBEM005, 18 th International Congress of Mechanical Engineering, Brazilian Society of Engineering an Mechanical Sciences, Ouro Preto, Brazil. Cointe, R. an Arman, 1987, Hyroynamic Impact Analysis of a Cyliner, Journal of Offshore Mechanics an Arctic Engineering, Vol. 109, pp Cointe, R., ontaine, E., Molin, B. an Scolan, Y.M., 004, On Energy Arguments Applie to the Hyroynamic Impact orce, Journal of Engineering Mathematics, Vol.48, pp Lamb, H., 193, Hyroynamics, Dover Publications, 6 th e., 738 pp. Molin, B., Cointe, R. an ontaine, E., 1996, On Energy Arguments Applie to the Slamming orce, 11 st.international Workshop on Water Waves an loating Boies, Hamburg. Pesce, C.P., 003, The Application of Lagrange Equations to Mechanical Systems with Mass Explicitly Depenent on Position, Journal of Applie Mechanics, Vol. 70, pp Pesce, C.P., 006, A Note on the Classical ree Surface Hyroynamic Impact Problem, Worl Scientific, Hyroynamics Symposium in honor to T.-Y. Wu, to appear. Scolan, Y.M. an Korobkin, A.A., 003, Energy Distribution from Vertical Impact of a Three- Dimensional Soli Boy onto the lat ree Surface of an Ieal Liqui, Journal of luis an Structures, Vol.17, pp Wu, G. X., 1998, Hyroynamic orce on a Rigi Boy uring Impact with Liqui, Journal of luis an Structures, Vol. 1, pp

5 Casetta, L., Pesce, C.P. The proper efinition of the ae mass for the water entry problem Discusser - R.C.T. Rainey: I isagree with your equation (1). If the flui has infinite epth, the flui momentum is not efine, see beginning of chapter 6 in Lamb (193). If the flui has finite epth, there is a reaction force on the seabe, which oes not ten to zero as the epth increases. Reply: irst, we agree that in the case of finite epth there woul be a reaction force on the seabe. Moreover the ae mass woul be also a function of the proximity to the seabe. Even in that finite epth case, however, the flui omain might be unboune. We note that we are treating the infinite epth case only. The momentum (an the kinetic energy) of the liqui may be efine whenever the velocity fiel is integrable (square integrable, i.e, integrable in the energy norm) over the whole liqui omain. This is the case in the present paper.

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const. G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

A note on the Mooney-Rivlin material model

A note on the Mooney-Rivlin material model A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel

More information

A. Incorrect! The letter t does not appear in the expression of the given integral

A. Incorrect! The letter t does not appear in the expression of the given integral AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

Conservation laws a simple application to the telegraph equation

Conservation laws a simple application to the telegraph equation J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Chapter Primer on Differentiation

Chapter Primer on Differentiation Capter 0.01 Primer on Differentiation After reaing tis capter, you soul be able to: 1. unerstan te basics of ifferentiation,. relate te slopes of te secant line an tangent line to te erivative of a function,.

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

PDE Notes, Lecture #11

PDE Notes, Lecture #11 PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =

More information

The continuity equation

The continuity equation Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

Numerical Integrator. Graphics

Numerical Integrator. Graphics 1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,

More information

Centrum voor Wiskunde en Informatica

Centrum voor Wiskunde en Informatica Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.

More information

Abstract A nonlinear partial differential equation of the following form is considered:

Abstract A nonlinear partial differential equation of the following form is considered: M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent

More information

On the number of isolated eigenvalues of a pair of particles in a quantum wire

On the number of isolated eigenvalues of a pair of particles in a quantum wire On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

Completely passive natural convection

Completely passive natural convection Early View publication on wileyonlinelibrary.com (issue an page numbers not yet assigne; citable using Digital Object Ientifier DOI) ZAMM Z. Angew. Math. Mech., 1 6 (2011) / DOI 10.1002/zamm.201000030

More information

7.1 Support Vector Machine

7.1 Support Vector Machine 67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

A Second Time Dimension, Hidden in Plain Sight

A Second Time Dimension, Hidden in Plain Sight A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will

More information

NON-SMOOTH DYNAMICS USING DIFFERENTIAL-ALGEBRAIC EQUATIONS PERSPECTIVE: MODELING AND NUMERICAL SOLUTIONS. A Thesis PRIYANKA GOTIKA

NON-SMOOTH DYNAMICS USING DIFFERENTIAL-ALGEBRAIC EQUATIONS PERSPECTIVE: MODELING AND NUMERICAL SOLUTIONS. A Thesis PRIYANKA GOTIKA NON-SMOOTH DYNAMICS USING DIFFERENTIAL-ALGEBRAIC EQUATIONS PERSPECTIVE: MODELING AND NUMERICAL SOLUTIONS A Thesis by PRIYANKA GOTIKA Submitte to the Office of Grauate Stuies of Texas A&M University in

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation

More information

The Hamiltonian structure of a 2-D rigid cylinder interacting dynamically with N point vortices

The Hamiltonian structure of a 2-D rigid cylinder interacting dynamically with N point vortices The Hamiltonian structure of a -D rigi cyliner interacting ynamically with N point vortices Banavara Shashikanth, Jerrol Marsen, Joel Burick, Scott Kelly Control an Dynamical Systems, 07-8, Mechanical

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

arxiv: v2 [math.ds] 26 Apr 2017

arxiv: v2 [math.ds] 26 Apr 2017 GAUSSIAN CURVATURE AND GYROSCOPES GRAHAM COX AND MARK LEVI arxiv:1607.03217v2 [math.ds] 26 Apr 2017 Abstract. We relate Gaussian curvature to the gyroscopic force, thus giving a mechanical interpretation

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

A look at Einstein s clocks synchronization

A look at Einstein s clocks synchronization A look at Einstein s clocks synchronization ilton Penha Departamento e Física, Universiae Feeral e Minas Gerais, Brasil. nilton.penha@gmail.com Bernhar Rothenstein Politehnica University of Timisoara,

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Basic Thermoelasticity

Basic Thermoelasticity Basic hermoelasticity Biswajit Banerjee November 15, 2006 Contents 1 Governing Equations 1 1.1 Balance Laws.............................................. 2 1.2 he Clausius-Duhem Inequality....................................

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

Notes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk

Notes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk Notes on Lie Groups, Lie algebras, an the Exponentiation Map Mitchell Faulk 1. Preliminaries. In these notes, we concern ourselves with special objects calle matrix Lie groups an their corresponing Lie

More information

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (

More information

This section outlines the methodology used to calculate the wave load and wave wind load values.

This section outlines the methodology used to calculate the wave load and wave wind load values. COMPUTERS AND STRUCTURES, INC., JUNE 2014 AUTOMATIC WAVE LOADS TECHNICAL NOTE CALCULATION O WAVE LOAD VALUES This section outlines the methoology use to calculate the wave loa an wave win loa values. Overview

More information

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (II)

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (II) Laplace s Equation in Cylinrical Coorinates an Bessel s Equation (II Qualitative properties of Bessel functions of first an secon kin In the last lecture we foun the expression for the general solution

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz A note on asymptotic formulae for one-imensional network flow problems Carlos F. Daganzo an Karen R. Smilowitz (to appear in Annals of Operations Research) Abstract This note evelops asymptotic formulae

More information

ME338A CONTINUUM MECHANICS

ME338A CONTINUUM MECHANICS global vs local balance equations ME338A CONTINUUM MECHANICS lecture notes 11 tuesay, may 06, 2008 The balance equations of continuum mechanics serve as a basic set of equations require to solve an initial

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Fluid Mechanics EBS 189a. Winter quarter, 4 units, CRN Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 Veihmeyer Hall.

Fluid Mechanics EBS 189a. Winter quarter, 4 units, CRN Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 Veihmeyer Hall. Flui Mechanics EBS 189a. Winter quarter, 4 units, CRN 52984. Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 eihmeyer Hall. Course Description: xioms of flui mechanics, flui statics,

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences.

. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences. S 63 Lecture 8 2/2/26 Lecturer Lillian Lee Scribes Peter Babinski, Davi Lin Basic Language Moeling Approach I. Special ase of LM-base Approach a. Recap of Formulas an Terms b. Fixing θ? c. About that Multinomial

More information

CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS

CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN

More information

Systems & Control Letters

Systems & Control Letters Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information

Chapter 2 Governing Equations

Chapter 2 Governing Equations Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Short Intro to Coordinate Transformation

Short Intro to Coordinate Transformation Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent

More information

Markov Chains in Continuous Time

Markov Chains in Continuous Time Chapter 23 Markov Chains in Continuous Time Previously we looke at Markov chains, where the transitions betweenstatesoccurreatspecifietime- steps. That it, we mae time (a continuous variable) avance in

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

DIFFERENTIAL GEOMETRY OF CURVES AND SURFACES 7. Geodesics and the Theorem of Gauss-Bonnet

DIFFERENTIAL GEOMETRY OF CURVES AND SURFACES 7. Geodesics and the Theorem of Gauss-Bonnet A P Q O B DIFFERENTIAL GEOMETRY OF CURVES AND SURFACES 7. Geoesics an the Theorem of Gauss-Bonnet 7.. Geoesics on a Surface. The goal of this section is to give an answer to the following question. Question.

More information

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E. AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 14 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering an Applie Science Gonzaga University SDC

More information

2-7. Fitting a Model to Data I. A Model of Direct Variation. Lesson. Mental Math

2-7. Fitting a Model to Data I. A Model of Direct Variation. Lesson. Mental Math Lesson 2-7 Fitting a Moel to Data I BIG IDEA If you etermine from a particular set of ata that y varies irectly or inversely as, you can graph the ata to see what relationship is reasonable. Using that

More information

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything

More information

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Bair, faculty.uml.eu/cbair University of Massachusetts Lowell 1. Pre-Einstein Relativity - Einstein i not invent the concept of relativity,

More information

Control Volume Derivations for Thermodynamics

Control Volume Derivations for Thermodynamics Control olume Derivations for Thermoynamics J. M. Powers University of Notre Dame AME 327 Fall 2003 This ocument will give a summary of the necessary mathematical operations necessary to cast the conservation

More information

1. Filling an initially porous tube under a constant head imposed at x =0

1. Filling an initially porous tube under a constant head imposed at x =0 Notes on Moving Bounary problems, Voller U o M, volle00@umn.eu. Filling an initially porous tube uner a constant hea impose at x =0 Governing equation is base on calculating the water volume lux by the

More information

2.20 Marine Hydrodynamics Lecture 3

2.20 Marine Hydrodynamics Lecture 3 2.20 Marine Hyroynamics, Fall 2018 Lecture 3 Copyright c 2018 MIT - Department of Mechanical Engineering, All rights reserve. 1.7 Stress Tensor 2.20 Marine Hyroynamics Lecture 3 1.7.1 Stress Tensor τ ij

More information

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Convective heat transfer

Convective heat transfer CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,

More information

Zachary Scherr Math 503 HW 3 Due Friday, Feb 12

Zachary Scherr Math 503 HW 3 Due Friday, Feb 12 Zachary Scherr Math 503 HW 3 Due Friay, Feb 1 1 Reaing 1. Rea sections 7.5, 7.6, 8.1 of Dummit an Foote Problems 1. DF 7.5. Solution: This problem is trivial knowing how to work with universal properties.

More information

Momentum and Energy. Chapter Conservation Principles

Momentum and Energy. Chapter Conservation Principles Chapter 2 Momentum an Energy In this chapter we present some funamental results of continuum mechanics. The formulation is base on the principles of conservation of mass, momentum, angular momentum, an

More information

Statics. There are four fundamental quantities which occur in mechanics:

Statics. There are four fundamental quantities which occur in mechanics: Statics Mechanics isabranchofphysicsinwhichwestuythestate of rest or motion of boies subject to the action of forces. It can be ivie into two logical parts: statics, where we investigate the equilibrium

More information

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK AIAA Guiance, Navigation, an Control Conference an Exhibit 5-8 August, Monterey, California AIAA -9 VIRTUAL STRUCTURE BASED SPACECRAT ORMATION CONTROL WITH ORMATION EEDBACK Wei Ren Ranal W. Bear Department

More information

Gyroscopic matrices of the right beams and the discs

Gyroscopic matrices of the right beams and the discs Titre : Matrice gyroscopique es poutres roites et es i[...] Date : 15/07/2014 Page : 1/16 Gyroscopic matrices of the right beams an the iscs Summary: This ocument presents the formulation of the matrices

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The effect of nonvertical shear on turbulence in a stably stratified medium The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:

More information

DECOMPOSITION OF POLYNOMIALS AND APPROXIMATE ROOTS

DECOMPOSITION OF POLYNOMIALS AND APPROXIMATE ROOTS DECOMPOSITION OF POLYNOMIALS AND APPROXIMATE ROOTS ARNAUD BODIN Abstract. We state a kin of Eucliian ivision theorem: given a polynomial P (x) an a ivisor of the egree of P, there exist polynomials h(x),

More information

ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS

ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS ADRIAN CONSTANTIN AND BORIS KOLEV Abstract. Accoring to the principle of least action, the spatially perioic motions of one-imensional

More information