Session 5: Review of Classical Astrodynamics
|
|
- Beatrice Fitzgerald
- 5 years ago
- Views:
Transcription
1 Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve as constraints, is the change in velocity, Δv. The urose of the roulsion system is to rovide this Δv, which deends on the articular trajectory followed by the sacecraft and is determined by the orbital mechanics at lay. Because of this, the use of concets in astrodynamics becomes very relevant to sace roulsion. In this Lecture, we resent a short review of the fundamental laws in astrodynamics (Further reading: R. Battin, An Introduction to the Mathematics and Methods of Astrodynamics, AIAA Education Series, 1999). We start by analyzing the two-body roblem with central forces. Such analysis rovides the fundamentals of most orbital mechanics roblems. Consider the situation deicted in the figure below: x z m 1 r 1 r = r r 1 r y m We have two oint masses interacting through central (gravitational) forces. Under no external forces, the equations of motion for these articles are, d r 1 r r 1 m 1 = Gm 1 m dt r r 1 3 (1) d r r r 1 m = Gm 1 m dt r r 1 3 Defining the relative osition vector r = r r 1, we subtract Eqs. (1) and obtain a single equation of motion, d r r d v r + µ = 0 or + µ = 0 () dt r 3 dt r 3 where µ = G(m 1 + m ) and v = d r/dt. In most cases of interest in astrodynamics, one of the masses dominates over the other. For the articular situation in the figure, we could assume the motion of m 1 is negligible for m 1» m, therefore µ = Gm 1. This is a erfectly accetable aroximation for orbital motion about lanetary bodies. In the case of earth, µ E = 398,601 km 3 /s. This innocently looking equation () ket occuied several of the most rodigious minds in mathematics, like Gauss and Euler, ever since Newton roclaimed his gravitational law and in arallel with Leibniz, develoed infinitesimal calculus. In general, there are two ways 1
2 to tackle this equation: i) through numerical integration, by taking advantage of today s comutational ower (a luxury, that neither Gauss of Euler ever hoed to ta into) and ii) analytically, by uncovering a number of integrals of motion related to the structure of this non-linear differential equation. The analytical method, as exected, is much more enlightening and will be exlored in some detail in this lecture. We begin by taking the cross roduct of vector r with Eq. (), ( ) d v r d v d( r v) d r d( r v) r + µ = r = v = = 0 (3) dt r 3 dt dt dt dt where r r = v v = 0. From here, we see that the resulting vector r v must be constant. We define this constant as, h = r v (4) which reresents a massless angular momentum, which is conserved when no external forces are alied to the masses. The angular momentum is our first integral of motion. We now take Eq. () and calculate its cross roduct with h, d ( v ) r h + µ h = 0 (5) dt r 3 We use the vector identity r h = r ( r v) = r( r v) r v and notice that r v = rv r, where v r = dr/dt is the radial comonent of the velocity vector. The second term in Eq. (5) can be written as (excluding µ), Using this result, we can write Eq. (5) as, r ( ) h v r v d r = r = (6) r 3 r r dt r d dt This equation can be directly integrated, ( d r v ) ( ) h µ = 0 (7) dt r r v h µ = µ e (8) r where the constant (vector) of integration was selected as µ e for convenience, since later on will be identified as the eccentricity of a geometrical orbit (times µ). The eccentricity is our second integral of motion. We calculate the magnitude of the eccentricity vector, e e = ( v ) h r µ r (9)
3 Exanding the square on the right hand side, v ) ( ) h r v h = µ r µ ( ( ) v h r + 1 µr (10) Here we note that ( v h) = v h since v h. For the last term in Eq. (10), we use the vector identity v h = v ( r v) = rv rv r v, and therefore, ( v h) r = r (v v r) = r v θ = h The magnitude of the eccentricity vector becomes, e v h h 1 µ µ ( ) = + 1 E v = e (11) µ T 1 µr r h which is known as the vis-viva integral and reresents the total energy ET (er unit mass), our third and last integral of motion. In this way, without exlicitly solving the differential equation Eq. (), it is ossible to write exressions for the two-body roblem in terms of integrals of motion. To see this, use Eq. (8) to calculate the dot roduct e r, ( v h) r h e r = r = (1) µ µ r Define f as the angle between the eccentricity and radial vectors (known as the true anomaly) and solve for r, h /µ r = or r = (13) 1 + e cos f 1 + e cos f Eq. (13), known as the equation of orbit, is the solution in sace (not in time) of the differential Eq. (). It describes the relative lanar motion of the two oint masses considered in the roblem. A closer look at Eqs. (11) and (13) reveals three motion regimes that deend on the magnitude of the eccentricity vector. 0 e < 1 E T < 0 bounded ellitic e = 1 E T = 0 unbounded (zero velocity at ) arabolic e > 1 E T > 0 unbounded hyerbolic Unbounded trajectories are interesting, esecially in sace roulsion where in some instances the objective is to achieve an escae trajectory for exloration missions. However, for the moment, we focus on bounded (ellitical) orbits. Eventually, we will look at the result of erturbing these orbits through small forces rovided by electric roulsion. The lanar motion in the bounded case (0 e < 1) defines an ellitical trajectory (circular, if e = 0), as shown in the figure below. From Eq. (13), we see that when f = π/ the radius 3
4 ! r f! e a becomes h /µ. This quantity is tyically reresented by and is known as the arameter of the orbit. The eriaxis and aoaxis can be found directly from Eq. (13) with f = 0 and f = π, resectively, r = and r a = (14) 1 + e 1 e Since (twice) the semi-major axis of the orbit is a = r a + r, we find, Consequently, the orbital energy in Eq. (11) can be written as, = a(1 e ) (15) µ E T = (16) a Given the eccentricity and semi-major axis, we use Eq. (13) to find the orbital radius as a function of the true anomaly. Then, we use Eq. (11) to find the orbital velocity at that articular location. What remains is orbital timing. In other words, given an initial time and osition, find where the orbiting object will be some time afterwards. Use the conservation of angular momentum, h = r (df/dt), and substitute in the equation of orbit, Aroriately, this equation is known as Keler s equation. df µ = dt (17) (1 + e cos f) 3 The eriod of the orbit could be easily found by integrating Keler s law (orbit swees equal areas in equal u times) since the area of the ellise is πab and the semi-minor axis is b = a 1 e = h a/µ, the familiar result is, a T = π 3 (18) µ We mentioned earlier that a number of brilliant mathematicians worked intensively in solving the two-body roblem. In fact, most of that work was devoted to solve Keler s equation (17). Once the orbital timing is found, the solution of the two-body roblem is comlete. There are a number of different ways in which this equation could be written, each one with its own rocedure to find an aroximate solution. In general, Eq. (17) could be solved 4
5 numerically, but of course, little is gained that way since, after all, the original equation of motion Eq. () could also be solved numerically. At the end, aroximate solutions using iterative methods, are extremely accurate. It is of course, more interesting for our articular alication to solve the equations of motion Eq. (1) when there are forces acting on the masses (from thrusters, for instance). In articular, if we assume a force F is alied to m, the resulting equation of motion will be, d v r + µ = a c ( r, t) (19) dt r 3 where a c = F /m is the acceleration due to external forces (other than gravity) acting on the moving mass. This new equation brings no additional difficulty if the roblem is solved numerically. However, our analytical aroach is criled since it is no longer ossible to find integrals of motion (neither the angular momentum, eccentricity or energy are longer constant). It is said that the resulting trajectories are non-kelerian. Desite this, if the acceleration is a small quantity (which is tyical for the tye of maneuvers using electric roulsion), the constants of motion will drift slowly. This allows us to follow these changes and rovide accurate redictions of the actual trajectory. Not surrisingly, this lies in the realm of the method of variation of arameters. For now, let us finalize this review by writing the vectorial equation of motion, Eq. (19), as two differential equations in olar coordinates. This will be useful later on to find analytical solutions to some simle roblems where small accelerations are imarted and motion is confined to a constant orbital lane: ( ) d r dθ µ d θ dr dθ a θ r + = a r and + = (0) dt dt r dt r dt dt r 5
6 MIT OenCourseWare htt://ocw.mit.edu 16.5 Sace Proulsion Sring 015 For information about citing these materials or our Terms of Use, visit: htt://ocw.mit.edu/terms.
-aε Lecture 4. Subjects: Hyperbolic orbits. Interplanetary transfer. (1) Hyperbolic orbits
6.50 Lecture 4 ubjects: Hyerbolic orbits. Interlanetary transfer. () Hyerbolic orbits The trajectory is still described by r, but now we have ε>, so that the +! cos" radius tends to infinity at the asymtotic
More informationCentral Force Motion Challenge Problems
Central Force Motion Challenge Problems Problem 1: Ellitic Orbit A satellite of mass m s is in an ellitical orbit around a lanet of mass m which is located at one focus of the ellise. The satellite has
More informationA Closed-Form Solution to the Minimum V 2
Celestial Mechanics and Dynamical Astronomy manuscrit No. (will be inserted by the editor) Martín Avendaño Daniele Mortari A Closed-Form Solution to the Minimum V tot Lambert s Problem Received: Month
More informationSession 12 : Monopropellant Thrusters
Session 12 : Monoroellant Thrusters Electrothermal augmentation of chemical rockets was the first form of electric roulsion alied in sace vehicles. In its original imlementation, resistojets were used
More informationAN ALTERNATIVE DESCRIPTION OF SWING-BY TRAJECTORIES IN TWO AND THREE DIMENSIONS
ADVANCES IN SPACE DYNAMICS 4: CELESTIAL MECHANICS AND ASTRONAUTICS, H. K. Kuga, Editor, 1-14 (004). Instituto Nacional de Pesquisas Esaciais INPE, São José dos Camos, SP, Brazil. ISBN 85-17-0001-9 AN ALTERNATIVE
More informationINDIRECT PLANETARY CAPTURE VIA PERIODIC ORBITS ABOUT LIBRATION POINTS
INDIRECT PLANETARY CAPTURE VIA PERIODIC ORBITS ABOUT LIBRATION POINTS Li Xiangyu 1,2, Qiao Dong 1,2, Cui Pingyuan 1,2 (1. Institute of Dee Sace Exloration Technology, Beijing Institute of Technology, Beijing,
More informationWhen solving problems involving changing momentum in a system, we shall employ our general problem solving strategy involving four basic steps:
10.9 Worked Examles 10.9.1 Problem Solving Strategies When solving roblems involving changing momentum in a system, we shall emloy our general roblem solving strategy involving four basic stes: 1. Understand
More informationChapter 1 Fundamentals
Chater Fundamentals. Overview of Thermodynamics Industrial Revolution brought in large scale automation of many tedious tasks which were earlier being erformed through manual or animal labour. Inventors
More informationintegral invariant relations is not limited to one or two such
The Astronomical Journal, 126:3138 3142, 2003 December # 2003. The American Astronomical Society. All rights reserved. Printed in U.S.A. EFFICIENT ORBIT INTEGRATION BY SCALING AND ROTATION FOR CONSISTENCY
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationTrajectory Optimization of Multi-Asteroids Exploration with Low Thrust
Trans. Jaan Soc. Aero. Sace Sci. Vol. 52, No. 175,. 47 54, 2009 Technical Note Trajectory Otimization of Multi-Asteroids Exloration with Low Thrust By Kaijian ZHU, Fanghua JIANG, Junfeng LI and Hexi BAOYIN
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationCharacterizing planetary orbits and the trajectories of light in the Schwarzschild metric
St. John Fisher College Fisher Digital Publications Physics Faculty Publications Physics 4-9-200 Characterizing lanetary orbits and the trajectories of light in the Schwarzschild metric Foek T. Hioe Saint
More informationA PIEZOELECTRIC BERNOULLI-EULER BEAM THEORY CONSIDERING MODERATELY CONDUCTIVE AND INDUCTIVE ELECTRODES
Proceedings of the 6th International Conference on Mechanics and Materials in Design, Editors: J.F. Silva Gomes & S.A. Meguid, P.Delgada/Aores, 26-30 July 2015 PAPER REF: 5513 A PIEZOELECTRIC BERNOULLI-EULER
More informationFind the equation of a plane perpendicular to the line x = 2t + 1, y = 3t + 4, z = t 1 and passing through the point (2, 1, 3).
CME 100 Midterm Solutions - Fall 004 1 CME 100 - Midterm Solutions - Fall 004 Problem 1 Find the equation of a lane erendicular to the line x = t + 1, y = 3t + 4, z = t 1 and assing through the oint (,
More informationUse of Orbiting Reflectors to Decrease the Technological Challenges of Surviving the Lunar Night
IAC-11-A5.1.11 Use of Orbiting Reflectors to Decrease the Technological Challenges of Surviving the Lunar Night Russell Bewick*, Joan-Pau Sanchez and Colin R. McInnes Advanced Sace Concets Laboratory,
More informationVI. Electrokinetics. Lecture 31: Electrokinetic Energy Conversion
VI. Electrokinetics Lecture 31: Electrokinetic Energy Conversion MIT Student 1 Princiles 1.1 General Theory We have the following equation for the linear electrokinetic resonse of a nanochannel: ( ) (
More informationSession 6: Analytical Approximations for Low Thrust Maneuvers
Session 6: Analytical Approximations for Low Thrust Maneuvers As mentioned in the previous lecture, solving non-keplerian problems in general requires the use of perturbation methods and many are only
More information16.2. Infinite Series. Introduction. Prerequisites. Learning Outcomes
Infinite Series 6.2 Introduction We extend the concet of a finite series, met in Section 6., to the situation in which the number of terms increase without bound. We define what is meant by an infinite
More information9 The Theory of Special Relativity
9 The Theory of Secial Relativity Assign: Read Chater 4 of Carrol and Ostlie (2006) Newtonian hysics is a quantitative descrition of Nature excet under three circumstances: 1. In the realm of the very
More informationCHAPTER 25. Answer to Checkpoint Questions
CHAPTER 5 ELECTRIC POTENTIAL 68 CHAPTER 5 Answer to Checkoint Questions. (a) negative; (b) increase. (a) ositive; (b) higher 3. (a) rightward; (b),, 3, 5: ositive; 4: negative; (c) 3, then,, and 5 tie,
More informationQuadratic Residues, Quadratic Reciprocity. 2 4 So we may as well start with x 2 a mod p. p 1 1 mod p a 2 ±1 mod p
Lecture 9 Quadratic Residues, Quadratic Recirocity Quadratic Congruence - Consider congruence ax + bx + c 0 mod, with a 0 mod. This can be reduced to x + ax + b 0, if we assume that is odd ( is trivial
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationLATELY, numerous missions involving satellites flying in
JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 3, No., January February Secial Inclinations Allowing Minimal Drift Orbits for Formation Flying Satellites M. Sabatini University of Rome La Saienza, Rome,
More informationCHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules
CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is
More informationOn Doob s Maximal Inequality for Brownian Motion
Stochastic Process. Al. Vol. 69, No., 997, (-5) Research Reort No. 337, 995, Det. Theoret. Statist. Aarhus On Doob s Maximal Inequality for Brownian Motion S. E. GRAVERSEN and G. PESKIR If B = (B t ) t
More informationCalculation of gravity due to a vertical cylinder using a spherical harmonic series and numerical integration
CSIRO PUBISHING Exloration Geohysics htt://dx.doi.org/.7/eg43 Calculation of gravity due to a vertical cylinder using a sherical harmonic series and numerical integration Sung-Ho Na,3 Hyoungrea Rim,3,4
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationSPACE situational awareness (SSA) encompasses intelligence,
JOURNAL OF GUIDANCE,CONTROL, AND DYNAMICS Vol. 34, No. 6, November December 2011 Gaussian Sum Filters for Sace Surveillance: Theory and Simulations Joshua T. Horwood, Nathan D. Aragon, and Aubrey B. Poore
More informationIcarus 218 (2012) Contents lists available at SciVerse ScienceDirect. Icarus. journal homepage:
Icarus 8 () 4 4 Contents lists available at SciVerse ScienceDirect Icarus journal homeage: www.elsevier.com/locate/icarus The fate of sub-micron circumlanetary dust grains I: Aligned diolar magnetic fields
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Jul 1998
arxiv:cond-mat/98773v1 [cond-mat.stat-mech] 5 Jul 1998 Floy modes and the free energy: Rigidity and connectivity ercolation on Bethe Lattices P.M. Duxbury, D.J. Jacobs, M.F. Thore Deartment of Physics
More information0.1 Practical Guide - Surface Integrals. C (0,0,c) A (0,b,0) A (a,0,0)
. Practical Guide - urface Integrals urface integral,means to integrate over a surface. We begin with the stud of surfaces. The easiest wa is to give as man familiar eamles as ossible ) a lane surface
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationAvailability and Maintainability. Piero Baraldi
Availability and Maintainability 1 Introduction: reliability and availability System tyes Non maintained systems: they cannot be reaired after a failure (a telecommunication satellite, a F1 engine, a vessel
More informationReal beads on virtual strings: Charged particles on magnetic field lines
Real beads on virtual strings: Charged articles on magnetic field lines Boris Breizman and Vladimir Khudik a) Institute for Fusion Studies, The University of Texas at Austin, Texas 78712 (Received 2 March
More information3.4 Design Methods for Fractional Delay Allpass Filters
Chater 3. Fractional Delay Filters 15 3.4 Design Methods for Fractional Delay Allass Filters Above we have studied the design of FIR filters for fractional delay aroximation. ow we show how recursive or
More information16. CHARACTERISTICS OF SHOCK-WAVE UNDER LORENTZ FORCE AND ENERGY EXCHANGE
16. CHARACTERISTICS OF SHOCK-WAVE UNDER LORENTZ FORCE AND ENERGY EXCHANGE H. Yamasaki, M. Abe and Y. Okuno Graduate School at Nagatsuta, Tokyo Institute of Technology 459, Nagatsuta, Midori-ku, Yokohama,
More information8.7 Associated and Non-associated Flow Rules
8.7 Associated and Non-associated Flow Rules Recall the Levy-Mises flow rule, Eqn. 8.4., d ds (8.7.) The lastic multilier can be determined from the hardening rule. Given the hardening rule one can more
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More information16.2. Infinite Series. Introduction. Prerequisites. Learning Outcomes
Infinite Series 6. Introduction We extend the concet of a finite series, met in section, to the situation in which the number of terms increase without bound. We define what is meant by an infinite series
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationTemperature, current and doping dependence of non-ideality factor for pnp and npn punch-through structures
Indian Journal of Pure & Alied Physics Vol. 44, December 2006,. 953-958 Temerature, current and doing deendence of non-ideality factor for n and nn unch-through structures Khurshed Ahmad Shah & S S Islam
More informationLower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
More informationPHYS 301 HOMEWORK #9-- SOLUTIONS
PHYS 0 HOMEWORK #9-- SOLUTIONS. We are asked to use Dirichlet' s theorem to determine the value of f (x) as defined below at x = 0, ± /, ± f(x) = 0, - < x
More informationWhere: Where: f Wave s frequency (Hz) c Speed of light ( ms -1 ) Wavelength (m)
in a direction to both of the fields as shown in Figure 1. In wave model, the electromagnetic radiation is commonly associated with wavelength and frequency, exressed mathematically as: c f...(1) f Wave
More informationThermal Propellant Gauging System for BSS 601
5th AIAA International Communications Satellite Systems Conference (organized by APSCC) AIAA 007-3149 Thermal Proellant Gauging System for BSS 601 T. Narita. 1 JSAT Cor, 9-1 Miho-Cho, Midori-ku, Yokohama
More informationAnnouncements. p FINAL EXAM
Announcements FINAL EXAM n PHYS 131-001: Wednesday, 12/12 @ 8-9:50 am n PHYS 131-002: Monday, 12/10 @ 10-11:50 am NO New Homework! CQ6: 10 m/s 2 13.14: a) 3.0 x 10 24 kg b) 0.89 m/s 2 13.20: 1.5 x 10 4
More informationModeling and Estimation of Full-Chip Leakage Current Considering Within-Die Correlation
6.3 Modeling and Estimation of Full-Chi Leaage Current Considering Within-Die Correlation Khaled R. eloue, Navid Azizi, Farid N. Najm Deartment of ECE, University of Toronto,Toronto, Ontario, Canada {haled,nazizi,najm}@eecg.utoronto.ca
More informationRANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES
RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES AARON ZWIEBACH Abstract. In this aer we will analyze research that has been recently done in the field of discrete
More informationChurilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics
Churilova Maria Saint-Petersburg State Polytechnical University Deartment of Alied Mathematics Technology of EHIS (staming) alied to roduction of automotive arts The roblem described in this reort originated
More informationShadow Computing: An Energy-Aware Fault Tolerant Computing Model
Shadow Comuting: An Energy-Aware Fault Tolerant Comuting Model Bryan Mills, Taieb Znati, Rami Melhem Deartment of Comuter Science University of Pittsburgh (bmills, znati, melhem)@cs.itt.edu Index Terms
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More information14 th AAS/AIAA Space Flight Mechanics Conference
Paer AAS 04-157 ASPHERCAL FORMATONS NEAR THE LBRATON PONTS N THE SUN-EARTH/MOON EPHEMERS SYSTEM B. G. Marchand and K. C. Howell School of Aeronautics and Astronautics, Purdue University West Lafayette,
More informationdn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential
Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system
More informationChapter 7 Rational and Irrational Numbers
Chater 7 Rational and Irrational Numbers In this chater we first review the real line model for numbers, as discussed in Chater 2 of seventh grade, by recalling how the integers and then the rational numbers
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationAnalysis of cold rolling a more accurate method
Analysis of cold rolling a more accurate method 1.1 Rolling of stri more accurate slab analysis The revious lecture considered an aroximate analysis of the stri rolling. However, the deformation zone in
More informationUsing the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process
Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process P. Mantalos a1, K. Mattheou b, A. Karagrigoriou b a.deartment of Statistics University of Lund
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationMetrics Performance Evaluation: Application to Face Recognition
Metrics Performance Evaluation: Alication to Face Recognition Naser Zaeri, Abeer AlSadeq, and Abdallah Cherri Electrical Engineering Det., Kuwait University, P.O. Box 5969, Safat 6, Kuwait {zaery, abeer,
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationarxiv:cond-mat/ v2 25 Sep 2002
Energy fluctuations at the multicritical oint in two-dimensional sin glasses arxiv:cond-mat/0207694 v2 25 Se 2002 1. Introduction Hidetoshi Nishimori, Cyril Falvo and Yukiyasu Ozeki Deartment of Physics,
More informationCalculation of MTTF values with Markov Models for Safety Instrumented Systems
7th WEA International Conference on APPLIE COMPUTE CIENCE, Venice, Italy, November -3, 7 3 Calculation of MTTF values with Markov Models for afety Instrumented ystems BÖCÖK J., UGLJEA E., MACHMU. University
More informationNUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS
NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372
More information2016-r1 Physics 220: Worksheet 02 Name
06-r Physics 0: Worksheet 0 Name Concets: Electric Field, lines of force, charge density, diole moment, electric diole () An equilateral triangle with each side of length 0.0 m has identical charges of
More informationAnalysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia
roceedings World Geothermal Congress 00 Bali, Indonesia, 5-9 Aril 00 Analysis of ressure Transient Resonse for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia Jorge A.
More informationDimensional perturbation theory for Regge poles
Dimensional erturbation theory for Regge oles Timothy C. Germann Deartment of Chemistry, University of California, Berkeley, California 94720 Sabre Kais Deartment of Chemistry, Purdue University, West
More informationMODEL-BASED MULTIPLE FAULT DETECTION AND ISOLATION FOR NONLINEAR SYSTEMS
MODEL-BASED MULIPLE FAUL DEECION AND ISOLAION FOR NONLINEAR SYSEMS Ivan Castillo, and homas F. Edgar he University of exas at Austin Austin, X 78712 David Hill Chemstations Houston, X 77009 Abstract A
More informationPulse Propagation in Optical Fibers using the Moment Method
Pulse Proagation in Otical Fibers using the Moment Method Bruno Miguel Viçoso Gonçalves das Mercês, Instituto Suerior Técnico Abstract The scoe of this aer is to use the semianalytic technique of the Moment
More informationSolutions of the Duffing and Painlevé-Gambier Equations by Generalized Sundman Transformation
Solutions of the Duffing and Painlevé-Gambier Equations by Generalized Sundman Transformation D.K.K. Adjaï a, L. H. Koudahoun a, J. Akande a, Y.J.F. Komahou b and M. D. Monsia a 1 a Deartment of Physics,
More informationWolfgang POESSNECKER and Ulrich GROSS*
Proceedings of the Asian Thermohysical Proerties onference -4 August, 007, Fukuoka, Jaan Paer No. 0 A QUASI-STEADY YLINDER METHOD FOR THE SIMULTANEOUS DETERMINATION OF HEAT APAITY, THERMAL ONDUTIVITY AND
More informationClassical gas (molecules) Phonon gas Number fixed Population depends on frequency of mode and temperature: 1. For each particle. For an N-particle gas
Lecture 14: Thermal conductivity Review: honons as articles In chater 5, we have been considering quantized waves in solids to be articles and this becomes very imortant when we discuss thermal conductivity.
More informationNode-voltage method using virtual current sources technique for special cases
Node-oltage method using irtual current sources technique for secial cases George E. Chatzarakis and Marina D. Tortoreli Electrical and Electronics Engineering Deartments, School of Pedagogical and Technological
More informationMaximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming
Maximum Entroy and the Stress Distribution in Soft Disk Packings Above Jamming Yegang Wu and S. Teitel Deartment of Physics and Astronomy, University of ochester, ochester, New York 467, USA (Dated: August
More informationOn Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm
On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm Gabriel Noriega, José Restreo, Víctor Guzmán, Maribel Giménez and José Aller Universidad Simón Bolívar Valle de Sartenejas,
More informationSpacecraft Power System Controller Based on Neural Network
Proceedings of the 4 th International Middle East Power Systems Conference (MEPCON ), Cairo University, Egyt, December 9-2, 2, Paer ID 242. Sacecraft Power System Controller Based on Neural Network Hanaa
More informationOn Wald-Type Optimal Stopping for Brownian Motion
J Al Probab Vol 34, No 1, 1997, (66-73) Prerint Ser No 1, 1994, Math Inst Aarhus On Wald-Tye Otimal Stoing for Brownian Motion S RAVRSN and PSKIR The solution is resented to all otimal stoing roblems of
More informationOrdered and disordered dynamics in random networks
EUROPHYSICS LETTERS 5 March 998 Eurohys. Lett., 4 (6),. 599-604 (998) Ordered and disordered dynamics in random networks L. Glass ( )andc. Hill 2,3 Deartment of Physiology, McGill University 3655 Drummond
More informationChapter 6: Sound Wave Equation
Lecture notes on OPAC0- ntroduction to Acoustics Dr. Eser OLĞAR, 08 Chater 6: Sound Wave Equation. Sound Waves in a medium the wave equation Just like the eriodic motion of the simle harmonic oscillator,
More informationGeo-E2010 Advanced Soil Mechanics L Wojciech Sołowski. 07 March 2017
Geo-E2010 Advanced Soil Mechanics L Wojciech Sołowski 07 March 2017 Soil modeling: critical state soil mechanics and Modified Cam Clay model Outline 1. Refresh of the theory of lasticity 2. Critical state
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationSimple geometric interpretation of signal evolution in phase-sensitive fibre optic parametric amplifier
Simle geometric interretation of signal evolution in hase-sensitive fibre otic arametric amlifier A.A. REDYUK,,,* A.E. BEDNYAKOVA,, S.B. MEDVEDEV, M.P. FEDORUK,, AND S.K. TURITSYN,3 Novosibirsk State University,
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More informationHEAT, WORK, AND THE FIRST LAW OF THERMODYNAMICS
HET, ORK, ND THE FIRST L OF THERMODYNMIS 8 EXERISES Section 8. The First Law of Thermodynamics 5. INTERPRET e identify the system as the water in the insulated container. The roblem involves calculating
More informationAPPROXIMATIONS DES FORMULES
MÉMOIRE SUR LES APPROXIMATIONS DES FORMULES QUI SONT FONCTIONS DE TRÈS GRANDS NOMBRES (SUITE) P.S. Lalace Mémoires de l Académie royale des Sciences de Paris, year 1783; 1786. Oeuvres comlètes X,. 95 338
More informationModelling a Partly Filled Road Tanker during an Emergency Braking
Proceedings of the World Congress on Engineering and Comuter Science 217 Vol II, October 25-27, 217, San Francisco, USA Modelling a Partly Filled Road Tanker during an Emergency Braking Frank Otremba,
More informationWeather and Climate Laboratory Spring 2009
MIT OenCourseWare htt://ocw.mit.edu 12.307 Weather and Climate Laboratory Sring 2009 For information about citing these materials or our Terms of Use, visit: htt://ocw.mit.edu/terms. Thermal wind John
More informationIntroduction to Landau s Fermi Liquid Theory
Introduction to Landau s Fermi Liquid Theory Erkki Thuneberg Deartment of hysical sciences University of Oulu 29 1. Introduction The rincial roblem of hysics is to determine how bodies behave when they
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationCMSC 425: Lecture 7 Geometric Programming: Sample Solutions
CMSC 425: Lecture 7 Geometric Programming: Samle Solutions Samles: In the last few lectures, we have been discussing affine and Euclidean geometr, coordinate frames and affine transformations, and rotations.
More informationLong-term dynamics of high area-to-mass ratio objects in high-earth orbit
Available online at www.sciencedirect.com ScienceDirect Advances in Sace Research 52 (2013) 1545 1560 www.elsevier.com/locate/asr Long-term dynamics of high area-to-mass ratio objects in high-earth orbit
More informationRe-entry Protocols for Seismically Active Mines Using Statistical Analysis of Aftershock Sequences
Re-entry Protocols for Seismically Active Mines Using Statistical Analysis of Aftershock Sequences J.A. Vallejos & S.M. McKinnon Queen s University, Kingston, ON, Canada ABSTRACT: Re-entry rotocols are
More information8 STOCHASTIC PROCESSES
8 STOCHASTIC PROCESSES The word stochastic is derived from the Greek στoχαστικoς, meaning to aim at a target. Stochastic rocesses involve state which changes in a random way. A Markov rocess is a articular
More informationBayesian Spatially Varying Coefficient Models in the Presence of Collinearity
Bayesian Satially Varying Coefficient Models in the Presence of Collinearity David C. Wheeler 1, Catherine A. Calder 1 he Ohio State University 1 Abstract he belief that relationshis between exlanatory
More informationCFL Conditions for Runge-Kutta Discontinuous Galerkin Methods on Triangular Grids
CFL Conditions for Runge-Kutta Discontinuous Galerkin Methods on Triangular Grids T. Toulorge a,, W. Desmet a a K.U. Leuven, Det. of Mechanical Engineering, Celestijnenlaan 3, B-31 Heverlee, Belgium Abstract
More informationCharacteristics of Beam-Based Flexure Modules
Shorya Awtar e-mail: shorya@mit.edu Alexander H. Slocum e-mail: slocum@mit.edu Precision Engineering Research Grou, Massachusetts Institute of Technology, Cambridge, MA 039 Edi Sevincer Omega Advanced
More informationIntroduction to Arithmetic Geometry Fall 2013 Lecture #10 10/8/2013
18.782 Introduction to Arithmetic Geometry Fall 2013 Lecture #10 10/8/2013 In this lecture we lay the groundwork needed to rove the Hasse-Minkowski theorem for Q, which states that a quadratic form over
More information