Semicompressible Ocean Thermodynamics and Boussinesq Energy Conservation
|
|
- Drusilla Woods
- 5 years ago
- Views:
Transcription
1 Article Semicompressible Ocean Thermoynamics an Boussinesq Energy Conservation William K. Dewar 1, *, Joseph Schoonover 1, Trevor McDougall 2 an Rupert Klein 3 1 Department of Earth, Ocean, an Atmospheric Science, Floria State University, Tallahassee, FL 32309, USA; schoonover.numerics@gmail.com 2 School of Mathematics an Statistics, UNSW, Syney, NSW 2052, Australia; trevor.mcougall@unsw.eu.au 3 Institut für Mathematik, Freie Universität, Berlin 14195, Germany; rupert.klein@math.fu-berlin.e * Corresponence: wewar@fsu.eu; Tel.: Acaemic Eitor: Pavel S. Berloff Receive: 29 Januray 2016 ; Accepte: 24 March 2016 ; Publishe: 8 April 2016 Abstract: Equations more accurate than the Boussinesq set that still filter out soun were recently introuce. While these equations were shown to have a consistent potential energy, their thermoynamical behavior an associate implications were not fully analyze. These shortcomings are remeie in the present note that argues both sets are fully consistent from a thermoynamic perspective. It is further argue that both sets remain computationally competitive with the Boussinesq set. Keywors: thermoynamics; compressibility; equations of motion 1. Introuction Dewar et al. [1] recently evelope two sets of ocean ynamical equations that inclue some egree of compressibility an exclue soun waves. They were argue to be simultaneously more accurate than the Boussinesq equations an computationally competitive with them. These were referre to as the Type I an Type II Semicompressible equations. In subsequent communications with Prof. R. Klein, he correctly pointe out that the thermoynamic consistency of these equations ha not yet been emonstrate. Thermoynamic consistency in this context means that total energy is conserve an that an entropy variable exists for the system that obeys the secon thermoynamic law. This remark was mae in comparison to pseuoincompressible equations erive in [2] for a general equation of state in which a much more thorough thermoynamic analysis ha been unertaken. The purpose of this short communication is to analyze similar characteristics of both types of the semicompressible equations. 2. The Semicompressible Equations As iscusse in [1], the Type I equations are ρ ( t u + u u + 2Ω u) = p [ρ ρ r (z) + pρ r c 2 s ρ ]gk + ρ ν visc u t ρ + (ρ u) = 0 ρ = ρ(s A, Θ, P ) (1a) t Θ = Θo t S A = S o A where ensity, ρ, is evaluate using Conservative Temperature, Θ, Absolute Salinity, S A an the hyrostatic pressure, P, etermine from the reference ensity profile, ρ r (at most a function of z) accoring to (1b) Fluis 2016, 1, 9; oi: /fluis
2 Fluis 2016, 1, 9 2 of 8 z P = ρ r g (2) The quantity ν visc is viscosity an ν = 1/ρ. De Szoeke an Samelson [3] have shown that the hyrostatic Boussinesq equations are equivalent to the compressible hyrostatic Navier-Stokes equations written in pressure coorinates. Equation (1a) improves upon the classic hyrostatic balance in the hyrostatic limit through the presence of the soun spee. The ensity approximation it represents is a more accurate approximation of true ensity than ρ, ensity evaluate at the hyrostatic pressure. Equation (1) is also more flexible in that they can be use in non-hyrostatic settings. The notation X o in Equation (1) enotes the non-avective contributions to X. Absolute Salinity, being a purely conserve quantity (ignoring chemistry, see [4]), the form is clear S o A = ν F S (3) where F S is the iffusive flux of Absolute Salinity. We purposely leave the form for Θ o unspecifie, as part of the present exercise is to etermine it. Other notation is stanar. The equations in Equation (1) satisfy the energy equation t (ρ (K + h + I )) + (u(ρ (K + h + I )) + pu) = pr + ρ ν visc K ρ ɛ (4) with potential energy playe by the quantity h + I P ρ ρ P r = P = ρρ r b I gρ r P (5) The quantity ɛ = ν visc u 2 is the viscous issipation of kinetic energy, b I is the Type I buoyancy, the quantity R R = Θ ν t Θ ν S A t S A (6) enotes iffusive terms an Equation (1b) emonstrates the presence of compressibility in these equations. The Type II equations are where ρ II ( t u + u u + 2Ω u) = p (ρ ρ II )gk p c 2 sii t ρ II + (ρ II u) = (ρ II u) = 0 ρ = ρ(s A, Θ, P ) t Θ = Θo t S A = S o A ρ II (z) = ρ II (0)e z o g c 2 z sii gk + ρ II ν visc u (7a) (7b) (8) is either a chosen reference ensity profile yieling the soun spee profile c 2 s II or a profile constructe from a chosen soun spee profile, c 2 s II. In either case, ensity, ρ II, is a function of epth only an the hyrostatic pressure is compute from it z P = ρ II g (9)
3 Fluis 2016, 1, 9 3 of 8 Compressibility appears in these equations in Equation (7b) where the time inepenence of ρ II has been recognize. The Type II equations satisfy the energy constraint t (ρ II(K + h + II )) + (uρ II(K + h + II ) + pu) = ρii ν visc K ρ II ɛ (10) where the role of potential energy is now playe by h + II = P ρ ρ P II ρ 2 P = II b II gρ II P (11) The quantity b II is the Type II buoyancy. The avantages an isavantages of these equations are iscusse in [1]. Extracting Equations (1) an (7) from the full Navier Stokes equations is akin to substituting the Boussinesq ensity for the full ensity in most places, an retaining a reference ensity profile where neee for potential energy consistency. 3. Semicompressible Thermoynamics We now consier how to augment (1) an (7) with a consistent thermoynamic behavior. The First Law of Thermoynamics is t h ν t P = T t η + µ t S A = 1 ρ F Q + ɛ (12) where T is temperature, h(η, S A, P) is specific enthalpy, η is specific entropy, specific volume is enote by ν = 1/ρ, P is pressure an µ is the relative chemical potential of salt in seawater. The quantity ɛ is the heat of viscous energy issipation, an F Q is a generalize heat flux. Enthalpy is efine by h = e + Pν (13) where e is specific internal energy. Enthalpy from Equation (12) is seen to epen naturally on the variables η, S A, an P, an in this form plays the role of a thermoynamic potential from which all thermoynamic variables can be obtaine by ifferentiation. For example, erivatives of enthalpy with respect to its natural variables are T, µ, an ν, respectively. The erivation of Equation (12) is iscusse extensively in the TEOS-10 manual ([4]; see also [5,6]). Equation (12) is exact. The semicompressible equations consier moifications to the thermoynamic variables cause mainly by the ifference between static an full pressure, an we analyze Equation (12) from that perspective. In what follows, we will aopt entropy, salinity, an pressure as state variables, an will often cast pressure epenency in terms of static pressure. For example, h = h (S A, η, P ) + h P (S A, η, P )p + O(p 2 ) (14) where P is static pressure an p the ynamic pressure, given by p = P P. The superscript ( ) will enote quantities evaluate at the static pressure. Specific volume, ensity, etc. will be written in a manner similar to Equation (14). Note, accoring to Equation (12), Equation (14) is equivalent to h = h (S A, η, P ) + ν p (15)
4 Fluis 2016, 1, 9 4 of Type I Thermoynamics We first evelop what is a consistent statement of Type I thermoynamics. consiering the form of Equation (12) to O(p 2 ); i.e., We begin by where t (h + ν p) = (ν + p P ν) t P + ν t p + T+ t η + µ+ t S A + O(p 2 ) (16) T + = T + p P T (17) an similarly for µ +. Evaluating the erivative of the static pressure, Equation (16) becomes t (h + ν p) = (ν + p P ν)( wgρ r) + ν t p + T+ t η + µ+ t S A (18) While Equation (18) is literally correct to O(p 2 ), it will turn out that a less accurate approximation, to O(p), is require on the right han sie of Equation (12) T + t η + µ+ t S A = ν F Q + ɛ (19) The reason for this will become clear later, here we stress only that it assures thermoynamic consistency of the Type I equations Type I Mechanical Energy Equation The mechanical energy equation erive from Equation (1a) is ρ t K = u p wg(ρ ρ r + ρ r p ρ c 2 ) + ρ ν visc K ρ ɛ (20) s where K = (u u)/2 is the kinetic energy ensity. Introucing the gravitational potential Φ = gz an using Equation (18) ρ ( t K + Φ + h + pν ) = u p + t p + ρ ν visc K ρ ɛ + ρ (T + t η + µ+ t S A) (21) so where Using Equation (12) P h = ν (22) P P h 1 P = h o + νp = h o + P + ρ r b P (23) ρ r g b = g (ρ ρ r ) ρ (24) The quantity h o is recognize as the potential enthalpy efine by [6], an the last integral in Equation (23) is the so-calle ynamic enthalpy of [7]. The quantity is a constant reference surface pressure. The mile integral in Equation (23) reuces upon inspection to Φ. Thus, Equation (21) becomes ρ t (K + h o + h + I ) = pu + ρ ν visc K ρ ɛ + ρ (T + t η + µ+ t S A) (25)
5 Fluis 2016, 1, 9 5 of 8 Lastly, we exploit Equation (19) to obtain ρ t (K + h o + h + I ) = pu + ρ ν visc K F Q (26) which is in conservative form. This implies the Type I equations have a well-forme energy principle, provie they are augmente by Equations (18) an (19) as their thermoynamic equations. Note that the entropy in Equation (19) can be rewritten as ρ t η = FQ µ + F S T + + (F Q µ + F S ) ( 1 T + ) F S ( µ+ T + ) + ρ ɛ T + (27) which is entropy non-conservation in its familiar form. Hence, the Type I equations also have a consistent entropy variable. At this point, a practical ifficulty is that entropy is not a typical ocean moelling variable. Instea, McDougall s Conservative Temperature Θ = h o c o P (28) McDougall ([6]) is preferable an so it is necessary to work out a Conservative Temperature equation to complement (27). From Equation (18), one fins t h ν t P = 1 ρ F Q p P T t η p P µ t S A + ɛ (29) By efinition h o = h(η, S A, ) (30) so t h o = θ t η + µ o t S A (31) where θ is potential temperature an µ o = µ(η, S A, ). Potential enthalpy is thus relate to entropy via t η = t h o µ o θ t S A (32) From Equation (23) t h o = c o p t Θ = 1 ρ F Q (T + θ) t η (µ+ µ o ) t S A + ɛ (33) so using Equation (32) t Θ = ν F Q + ɛ c o P (1 + (T+ θ) ) θ (( T+ θ )µ o µ + ) t S A c o P (1 + (T+ θ) θ ) (34) The unerline terms in Equation (34) are quite small, accoring to [8]. Neglecting them leas to the statement that Conservative Temperature is very accurately portraye as a conserve quantity; however, all of those terms must be inclue in orer to ensure full energy conservation. It is seen, therefore, that the Type I equations yiel energy conservation, but the cost is the insertion of a large quantity (c o PΘ) into the energy equation. A similar result was labelle in [7] as a crippling isavantage of the Boussinesq-Bernoulli equation. It might thus be more useful to work with a partial equation when consiering energy; i.e., one that oesn t explicitly inclue Conservative Temperature, but in any case, the full potential enthalpy Equation (34) is require. We have also foun the Type I Semicompressible equations have a proper entropy variable an conserve total energy, an so are
6 Fluis 2016, 1, 9 6 of 8 fully thermoynamically consistent. These equations are very similar to the pseuo-incompressible equation set for general equations of state in [2] Type II Equations: Potential, Kinetic, an Thermal Energies The Type II equations are quite similar to the Type I equations except for the appearance of the reference ensity profile, rather than the Boussinesq ensity, in the mass conservation equation an in front of the momentum acceleration. Most of the analysis carries over in a manner similar to the Type I analysis, with the one important istinction that the velocity ivergence iagnose from mass conservation has no iabatic contributions u II = gw c 2 sii Therefore, the mechanical energy equation for the system analogous to Equation (21) is wg ρ II K = pu + p t c 2 wg(ρ ρ II + sii p c 2 sii (35) ) + ρ II ν visc K ρ II ɛ (36) with the unerline terms cancelling. Taking the material erivative of Equation (11) an substituting in Equation (36) yiels P [ ] ρ II t K = pu + ρ bii II P S A gρ r t S A (37) P [ ] bii + ρ II P η gρ r t η ρ II t h+ II + ρii ν visc K ρ II ɛ To put Equation (37) in conservative form requires the potential enthalpy equation t h = t h o + ν P t P + which is equivalent to Equation (33). This leas to ν P P S A t S ν A + η P t η (38) ρ II t (K + h+ II + co P Θ) = pu F Q + ρ II ν visc K (39) The Type II Conservative Temperature equation becomes t Θ = ν II F Q + ɛ (( T+ θ )µ o µ + ) t S A c o P T+ θ (40) The most significant ifference relative to the Type I set is in Equation (35), where the iabatic contributions to expansion are absent. Like for Type I, it is still necessary to retain all terms for Conservative Temperature, rather than the approximate form recommene in [6]. We conclue that the Type II equations also have a consistent mechanical an thermal energy structure. Our Type II equations are very similar to the thermoynamically consistent anelastic set erive in [9] A Comparison to the Boussinesq Set To ensure full energetic consistency for the semicompressible equations, it has been necessary to carry along some aitional non-conservative effects in the calculations. We now analyze the full energetics behavior of the Boussinesq set to gauge any aitional computational buren of the semicompressible set relative to them. If we aopt Conservative Temperature, salinity, an pressure as
7 Fluis 2016, 1, 9 7 of 8 the thermoynamic coorinates an follow the proceure outline in [7], we arrive at the Boussinesq version of the Bernoulli equation t (K + p + h B) = t p + P [ ] b P P [ S A gρ o t S A + Θ b gρ o ] P t Θ ɛ (41) (see Equation (48) in [7]) where b = g(ρ ρ o )/ρ is the Boussinesq buoyancy an ρ o is the Boussinesq reference ensity. Dynamic enthalpy as efine in [7] is h B = 1 P bp (42) gρ o Returning to a form of the potential enthalpy consistent with the Boussinesq approximation h (S A, Θ, P ) = h o + 1 P (1 + b )P (43) ρ o g where h o is potential enthalphy, an taking its material erivative yiels t h ν t P = c o P [ P t Θ + b S A gρ o P [ b + Θ gρ o ] P t S A (44) ] P t Θ = 1 ρ o F Q + ɛ Using Equation (44) to eliminate the pressure integrals in Equation (41) turns it into the conservative Boussinesq energy equation t (K + p + h B + c o P Θ) = t p + ρii ν visc K 1 ρ o F Q (45) The Boussinesq Conservative Temperature equation is given by t Θ = ν o F Q + ɛ (( T θ )µ o µ + ) t S A c o P T θ (46) The point of this section is that to obtain a conservative energy equation out of the Boussinesq set, it is necessary to carry along the same aitional iffusive terms as for the semi-compressible equations. A similar result has recently been erive in [10,11]. This is not one in any Boussinesq ocean circulation moel of which we are aware. 4. Summary The Type I an II semicompressible equations iscusse in [1], which are analogues of the pseuo-incompressible an anelastic equations known in meteorology, have unergone an examination of their thermoynamic structure. To summarize, the Types I an II equations are fully consistent from an energetic an thermoynamic perspective. To obtain this, it is necessary to inclue a more complete equation for Conservative Temperature as a part of the equation sets than is normally one. Interestingly, it turns out the Boussinesq equations must also carry analogous quantities to be in purely conservative energetic form. We thus see the semicompressible equations remain computationally competitive with the Boussinesq equations while amitting some egree of compressibility. Other results here inclue the associate efinitions of the ynamic enthalpy variables belonging to the sets. From a simple computer throughput perspective, this analysis recommens Type II relative to Type I because they are slightly less computationally intensive. However, the Type II equations
8 Fluis 2016, 1, 9 8 of 8 represent the oceanic ensity structure by a single vertical profile. Any such profile cannot perform well globally, so the cost of the simplicity is a less accurate ensity representation. This promises to be the most troubling when computing sea level changes via the continuity equation. The aitional accuracy provie by the Boussinesq ensity in the Type I set might well compensate for the (slight) aitional complexity. Any further exploration with these equation sets shoul bear in min the above points. Acknowlegments: This research was make possible by grants from the National Science Founation Grants OCE an OCE an the Gulf of Mexico Research Initiative through consortium funing to the University of Miami for the Consortium for Avance Research on Transport of Hyrocarbon in the Environment (CARTHE). William Young is thanke for many interesting conversations that ae significantly to the manuscript. Author Contributions: William K. Dewar an Rupert Klein conceive the research. Trevor McDougall an Joseph Schoonover contribute substantially to the thermoynamical analysis. William K. Dewar wrote the paper. Conflicts of Interest: The authors eclare no conflict of interest. References 1. Dewar, W.; Schoonover, J.; McDougall, T.; Young, W. Semicompressible ocean ynamics. J. Phys. Oceanogr. 2015, 45, Klein, R.; Pauluis, O. Thermoynamic consistency of a psueoincompressible approximation for general equations of state. J. Atmos. Sci. 2012, 69, De Szoeke, R.; Samelson, R. The uality between the Boussinesq an non-boussinesq hyrostatic equations of motion. J. Phys. Oceanogr. 2002, 32, Intergovernmental Oceanographic Commission (IOC); Scientific Committee on Oceanic Research (SCOR); International Association for the Physical Sciences of the Oceans (IAPSO). The International Thermoynamic Equation of Seawater 2010: Calculation an Use of Thermoynamic Properties; 2010 Manuals an Guies No. 56, UNESCO (English); Intergovernmental Oceanographic Commission: West Perth, Australia, 2010; p Lanau, L.M.; Lifshitz, E.M. Flui Mechanics, 1st e.; Pergamon Press: Oxfor, UK, McDougall, T. Potential enthalpy: A conservative oceanic variable for evaluating heat content an heat fluxes. J. Phys. Oceanogr. 2003, 33, Young, W. Dynamic enthalpy, conservative temperature an the seawater Boussinesq approximation. J. Phys. Oceanogr. 2010, 40, Graham, F.S.; McDougall, T.J. Quantifying the Nonconservative Prouction of Conservative Temperature, Potential Temperature, an Entropy. J. Phys. Oceanogr. 2013, 43, Pauluis, O. Thermoynamic consistency of the anelastic approximatio for a moist atmosphere. J. Atmos. Sci. 2008, 65, Tailleux, R. On the energetics of stratifie turbulent mixing, irreversible thermoynamics, Boussinesq moels an the ocean heat engine controversy. J. Flui Mech. 2009, 638, Tailleux, R. Irreversible compressible work an available potential energy issipation in turbulent stratifie fluis. Phys. Scr. 2013, oi: / /2013/t155/ c 2016 by the authors; licensee MDPI, Basel, Switzerlan. This article is an open access article istribute uner the terms an conitions of the Creative Commons by Attribution (CC-BY) license (
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 informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationConvective 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 informationOptimized 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 informationCompletely 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 informationThe 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 informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationSources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences
Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences
More informationThe 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 informationConservation 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 informationA Modification of the Jarque-Bera Test. for Normality
Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam
More informationSchrö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 informationfv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n
Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie
More informationAPPROXIMATE 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'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationPhysics 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 informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationCentrum 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 informationExercise 4 - Hydraulic Systems
Exercise 4 - Hyraulic Systems 4.1 Hyraulic Systems Hyraulic systems are, in general, escribe by the Navier-Stokes equations as you might have learne in flui ynamics courses. In orer to simplify the moeling
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationSeparation 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 informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationHow the potentials in different gauges yield the same retarded electric and magnetic fields
How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department
More informationLecture 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 informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationinflow 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 informationAsymptotics of a Small Liquid Drop on a Cone and Plate Rheometer
Asymptotics of a Small Liqui Drop on a Cone an Plate Rheometer Vincent Cregan, Stephen B.G. O Brien, an Sean McKee Abstract A cone an a plate rheometer is a laboratory apparatus use to measure the viscosity
More informationMath 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 informationMomentum 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 informationApplication 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 informationSection 2.7 Derivatives of powers of functions
Section 2.7 Derivatives of powers of functions (3/19/08) Overview: In this section we iscuss the Chain Rule formula for the erivatives of composite functions that are forme by taking powers of other functions.
More informationTotal Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*
51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*
More informationA 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 informationCONSERVATION 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 informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationθ 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 informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationA 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 informationChapter 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 informationMagnetic helicity evolution in a periodic domain with imposed field
PHYSICAL REVIEW E 69, 056407 (2004) Magnetic helicity evolution in a perioic omain with impose fiel Axel Branenburg* Norita, Blegamsvej 17, DK-2100 Copenhagen Ø, Denmark William H. Matthaeus University
More informationEquilibrium in Queues Under Unknown Service Times and Service Value
University of Pennsylvania ScholarlyCommons Finance Papers Wharton Faculty Research 1-2014 Equilibrium in Queues Uner Unknown Service Times an Service Value Laurens Debo Senthil K. Veeraraghavan University
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationGlobal Solutions to the Coupled Chemotaxis-Fluid Equations
Global Solutions to the Couple Chemotaxis-Flui Equations Renjun Duan Johann Raon Institute for Computational an Applie Mathematics Austrian Acaemy of Sciences Altenbergerstrasse 69, A-44 Linz, Austria
More informationNotes on the function gsw_internal_energy(sa,ct,p)
Notes on gsw_internal_energy Notes on the function gsw_internal_energy(sa,ct,p) The function gsw_internal_energy(sa,ct,p) evaluates the specific internal energy of seawater as a function of Absolute Salinity,
More informationDiagonalization 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 informationTime-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 informationAssignment 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 informationarxiv:nlin/ v1 [nlin.cd] 21 Mar 2002
Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical
More informationImpurities in inelastic Maxwell models
Impurities in inelastic Maxwell moels Vicente Garzó Departamento e Física, Universia e Extremaura, E-671-Baajoz, Spain Abstract. Transport properties of impurities immerse in a granular gas unergoing homogenous
More informationAgmon 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 informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationThe maximum sustainable yield of Allee dynamic system
Ecological Moelling 154 (2002) 1 7 www.elsevier.com/locate/ecolmoel The maximum sustainable yiel of Allee ynamic system Zhen-Shan Lin a, *, Bai-Lian Li b a Department of Geography, Nanjing Normal Uni ersity,
More informationLectures - 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 informationNonlinear Lagrangian equations for turbulent motion and buoyancy in inhomogeneous flows
Nonlinear Lagrangian equations for turbulent motion an buoyancy in inhomogeneous flows Stefan Heinz a) Delft University of Technology, Faculty of Applie Physics, Section Heat Transfer, Lorentzweg 1, 68
More informationG j dq i + G j. q i. = a jt. and
Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationComputing 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 informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More information05 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 informationON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS
ON ISENTROPIC APPROXIMATIONS FOR COMPRESSILE EULER EQUATIONS JUNXIONG JIA AND RONGHUA PAN Abstract. In this paper, we first generalize the classical results on Cauchy problem for positive symmetric quasilinear
More informationThermodynamic basis for a variational model for crystal growth
PHYSICL REIEW E OLUME 60, NUMBER 1 JULY 1999 Thermoynamic basis for a variational moel for crystal growth Bayar K. Johnson* an Robert F. Sekerka Carnegie Mellon University, Pittsburgh, Pennsylvania 1513
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationq = F If we integrate this equation over all the mass in a star, we have q dm = F (M) F (0)
Astronomy 112: The Physics of Stars Class 4 Notes: Energy an Chemical Balance in Stars In the last class we introuce the iea of hyrostatic balance in stars, an showe that we coul use this concept to erive
More informationThe influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section
IOP Conference Series: Materials Science an Engineering PAPER OPEN ACCESS The influence of the equivalent hyraulic iameter on the pressure rop preiction of annular test section To cite this article: A
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationLINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form
LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1 We consier ifferential equations of the form y + a()y = b(), (1) y( 0 ) = y 0, where a() an b() are functions. Observe that this class of equations inclues equations
More informationarxiv: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 informationThe Generalized Incompressible Navier-Stokes Equations in Besov Spaces
Dynamics of PDE, Vol1, No4, 381-400, 2004 The Generalize Incompressible Navier-Stokes Equations in Besov Spaces Jiahong Wu Communicate by Charles Li, receive July 21, 2004 Abstract This paper is concerne
More informationCOUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND
More information(2012) , ISBN
Cregan, V. an O'Brien, Stephen B.G. an McKee, Sean (2012) Asymptotics of a small liqui rop on a cone an plate rheometer. In: Progress in Inustrial Mathematics at ECMI 2010. Mathematics in Inustry: Progress
More informationA 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 informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationLagrangian and Hamiltonian Dynamics
Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion
More informationNovel entropy stable schemes for 1D and 2D fluid equations
Novel entropy stable schemes for D an D flui equations Eitan Tamor an Weigang Zhong Department of Mathematics, Center for Scientific Computation an Mathematical Moeling (CSCAMM, Institute of Physical Science
More informationOn the dynamics of the electron: Introduction, 1, 9
On the ynamics of the electron: Introuction,, 9 Henri oincaré Introuction. It seems at first that the aberration of light an relate optical an electrical phenomena will provie us with a means of etermining
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS
ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationLecture 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 informationFormulation of statistical mechanics for chaotic systems
PRAMANA c Inian Acaemy of Sciences Vol. 72, No. 2 journal of February 29 physics pp. 315 323 Formulation of statistical mechanics for chaotic systems VISHNU M BANNUR 1, an RAMESH BABU THAYYULLATHIL 2 1
More informationChapter 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 informationGeneralization 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 informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More informationConstruction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems
Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu
More informationChapter 11: Feedback and PID Control Theory
Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by
More information1.2 - Stress Tensor Marine Hydrodynamics Lecture 3
13.021 Marine Hyroynamics, Fall 2004 Lecture 3 Copyright c 2004 MIT - Department of Ocean Engineering, All rights reserve. 1.2 - Stress Tensor 13.021 Marine Hyroynamics Lecture 3 Stress Tensor τ ij:. The
More informationTopological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem
Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075
More informationSYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d
SYNTHESIS ON THE ASSESSMENT OF HIPS ONTRATION OEFFIIENT Marius MILEA Technical University Gheorghe Asachi of Iași, Machine Manufacturing an Inustrial Management Faculty, B-ul D Mangeron, nr 59, A, 700050,
More informationAbstract 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 informationFluid 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 informationState observers and recursive filters in classical feedback control theory
State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent
More informationarxiv: v1 [physics.class-ph] 20 Dec 2017
arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.
More informationContinuum Mechanics Lecture 4 Fluid dynamics
Continuum Mechanics Lecture 4 Flui ynamics Prof. http://www.itp.uzh.ch/~teyssier Outline - Flui kinematics - Mass an momentum conservation laws - The energy equation - Real fluis - Ieal fluis - Incompressible
More informationFundamental Laws of Motion for Particles, Material Volumes, and Control Volumes
Funamental Laws of Motion for Particles, Material Volumes, an Control Volumes Ain A. Sonin Department of Mechanical Engineering Massachusetts Institute of Technology Cambrige, MA 02139, USA August 2001
More informationPerfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs
Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Ashish Goel Michael Kapralov Sanjeev Khanna Abstract We consier the well-stuie problem of fining a perfect matching in -regular bipartite
More information