Expanded Equations for Torque and Force on a Cylindrical Permanent Magnet Core in a Large- Gap Magnetic Suspension System

Size: px
Start display at page:

Download "Expanded Equations for Torque and Force on a Cylindrical Permanent Magnet Core in a Large- Gap Magnetic Suspension System"

Transcription

1 NASA Technical Paper 3638 Expanded Equations for Torque and Force on a Cylindrical Permanent Magnet Core in a Large- Gap Magnetic Suspension System Nelson J. Groom Langley Research Center Hampton, Virginia National Aeronautics and Space Administration Langley Research Center Hampton, Virginia February 1997

2 Aailable electronically at the following URL address: Printed copies aailable from the following: NASA Center for AeroSpace Information National Technical Information Serice (NTIS) 800 Elkridge Landing Road 5285 Port Royal Road Linthicum Heights, MD Springfield, VA (301) (703)

3 Abstract The expanded equations for torque and force on a cylindrical permanent magnet core in a large-gap magnetic suspension system are presented. The core is assumed to be uniformly magnetized, and equations are deeloped for two orientations of the magnetization ector. One orientation is parallel to the axis of symmetry, and the other is perpendicular to this axis. Fields and gradients produced by suspension system electromagnets are assumed to be calculated at a point in inertial space which coincides with the origin of the core axis system in its initial alignment. Fields at a gien point in the core are defined by expanding the fields produced at the origin as a Taylor series. The assumption is made that the fields can be adequately defined by expansion up to second-order terms. Examination of the expanded equations for the case where the magnetization ector is perpendicular to the axis of symmetry reeals that some of the second-order gradient terms proide a method of generating torque about the axis of magnetization and therefore proide the ability to produce sixdegree-of-freedom control. Introduction This paper deelops the expanded equations for torque and force on a cylindrical permanent magnet core in a large-gap magnetic suspension system. The core is assumed to be uniformly magnetized, and equations are deeloped for two orientations of the magnetization ector. One orientation is parallel to the axis of symmetry, and the other is perpendicular to this axis. Fields and gradients produced by suspension system electromagnets are assumed to be calculated at a point in inertial space which coincides with the origin of the core axis system in its initial alignment with a reference inertial axis system. Fields at a gien point in the core are defined by expanding the fields produced at the origin as a Taylor series. The assumption is made that the fields can be adequately described by expansion up to second-order terms. The expansion of the fields and gradients is presented in appendix A. The equations for torques and forces on a magnetic core that are produced by a large-gap magnetic suspension system hae been presented and discussed in a number of papers. For example, see references 1 through 6. The torques on the core are usually approximated as a function of the external or applied fields at the centroid of the core, and the forces on the core are usually approximated as a function of the gradients of the applied fields at the centroid. It is generally assumed that terms that are a function of second-order or higher gradients of the applied fields at the centroid can be neglected. In practical applications that inole large-gap magnetic suspension systems, these assumptions hae proen to be alid. For an axisymmetric core, such as a cylinder, it can be shown that if the direction of magnetization is along the axis of symmetry, then the torque about that axis, produced by the applied fields and gradients of the applied fields, is always zero (ref. 3). Various methods of oercoming this constraint, which include shaping the core and using nonuniform three-dimensional magnetization, are discussed in references 4 and 7. Howeer, examination of the expanded equations for a cylindrical core reeals that for the case of uniform magnetization perpendicular to the axis of symmetry, some of the secondorder gradient terms proide a method of generating torque about the axis of magnetization and therefore the ability to produce six-degree-of-freedom control (ref. 8). For completeness, all gradient terms for expansion of the fields up to second order are presented. Finally, instead of deeloping the torque and force equations from a set of goerning equations that are a function of core olume, core magnetization ector, and suspension system fields and gradients, appendix B presents a deelopment that begins at a more fundamental leel in an attempt to proide better insight into the origin of these equations than is commonly aailable in the literature. Symbols A area, m 2 a B B radius of core, m magnetic flux density ector, T expanded magnetic flux density ector, T [ B] matrix of field gradients, T/m [ B ] matrix of expanded field gradients, T/m F δf I l total force ector on core, N force ector on incremental olume of core, N coil current ector, A length of core, m m magnetic moment ector, A-m 2

4 δm M Q m r T δt δt T m magnetic moment ector of dipole with incremental olume, A-m 2 magnetization ector, A/m pole strength, A-m position ector, m total torque ector on core, N-m torque ector on incremental olume of core, N-m torque on incremental olume of core about core origin, N-m inertial coordinate to suspended-element coordinate ector-transformation matrix U potential energy core olume, m 3 W δ work incremental olume, m 3 x, y, z coordinates in orthogonal axis system, m θ Euler orientation for 3, 2, 1 rotation sequence, rad gradient operator Subscripts: x, y, z components along x-, y-, z-axes, respectiely ij (ij)k Matrix notations: partial deriatie of i component in j-direction partial deriatie of ij partial deriatie in k-direction [ ] T transpose of matrix row ector A bar oer a symbol indicates that it is referenced to suspended-element coordinates. Magnetic Torques and Forces The torques and forces on a cylindrical permanent magnet core are deeloped in this section by integrating the equations for torques and forces on an incremental olume of the core with magnetic moment Mδ oer the core olume. These equations are deeloped in appendix B. Figure 1 shows the cylindrical core and the core coordinate system. The core coordinate system consists of a set of orthogonal x, y, z body-fixed axes that are initially aligned with a set of orthogonal x-, y-, z-axes fixed in inertial space. In order to define the fields and gradients z, z l/2 l/2 y, y x, x a Figure 1. Core coordinate system. 2

5 at any point in the core, the fields and gradients at the origin of the core axis system are expanded as a Taylor series. It is assumed that the fields can be adequately described by expansion up to second-order terms. The expanded fields and gradients in both inertial and core coordinates are presented in appendix A. For simplicity in deeloping the equations in this section, relatie motion between the core and the reference inertial coordinate system is assumed to be zero. This assumption remoes the requirement to transform between the inertial and core coordinate systems and eliminates a significant number of components which are small relatie to the fundamental terms in the equations when small-angle assumptions are used. In particular, the transformation of second-order gradient terms from inertial to core coordinates is ery complicated, as illustrated by equation (A14). The torque on an incremental olume of the core, about the core origin, can be written as δt δt + ( r δf) (1) where δt and δf are the torque and force on the incremental olume due to the field at that location, and r x y z is the position ector of the incremental olume (fig. 2). The total torque on the core can be written as T δt [ δt + ( r δf )] (2) where the integration is oer the core. Substituting equations (B20) and (B21) results in T {( M B ) + [ r ( M ) B ]} d The total force on the core can be written as F ( M )B d The term ( M )B can be written as (ref. 5) ( M )B [ B ]M (3) (4) (5) z δ r y x Figure 2. Incremental core olume. 3

6 where T y M x B z d + M x ( B xxz B xzx ) d (12) and B B xx B xy B xz B yx B yy B yz B zx B zy B zz (6) T z M x B y d + M x ( B xyx B xxy ) d Ealuating B xzy d first (13) M M x M y M z In equation (6) the notation f ij f i j is used. Magnetization Along Axis of Symmetry (7) For orientation of the magnetization ector along the axis of symmetry (x-axis) of the permanent magnet core, the only nonzero term in M is M x. Expanding equation (5) results in Substituting equation (8) into the second part of equation (3) and expanding results in Expanding the first part of equation (3) results in B xx [ B ]M M x B xy B xz ( B xyz + B xzy) r ( M ) B M x ( B xxz B xzx) 0 M B M x B z ( B xxy + B xyx) B y (8) (9) (10) B xzy d B xz y d + B ( xx)z xy d + B ( xy)z y 2 d + B ( xz)z zy d (14) where B xz has been expanded by using a Taylor series, as detailed in appendix A, and the notation f ( ij)k ( f i / j )/ k has been used. All the integrals inoling first-order terms in equation (14) are zero. Ealuating y 2 d yields l/ y 2 d y 2 ( z) a y 2 dy dz a y 1 ( z) l/2 (15) where y 1 ( z) y 2 ( z) a 2 z 2, a is the radius of the permanent magnet core, and l is the length (fig. 1). Since the area of the face of the permanent magnet core is A πa 2 and the olume is Al, equation (15) reduces to Substituting in equation (14) results in ( a 4 4)πl ( a 4 4)πl ( a 2 4) B xzy d ( a 2 4)B ( xy)z Ealuating B xyz d in a similar manner results in dx (16) (17) The components of equation (3) become 4 T x 0 + M x ( B xzy B xyz) d (11) B xyz d ( a 2 4)B ( xy)z (18) which is equal to equation (17). Therefore the torques about the x-axis due to second-order gradients cancel

7 out and T x 0 as expected. Going next to equation (12), Finally, noting that B (ij)k B (ik)j B (jk)i... and collecting terms, the components of torque become B z d B z + ( l 2 24)B ( zx)x + ( a 2 8)B ( zy)y + ( a 2 8)B ( zz)z (19) T y T x 0 M x B z M x [( a 2 4) ( l 2 8 ) ]B ( xx)z M x ( a 2 8) ( B ( yy)z + B ) ( zz)z (25) (26) and ( B xxz B xzx ) d [ ( a 2 4) ( l 2 12) ]B ( xx)z Substituting into equation (12) results in (20) T z M x B y M x l [( 8) ( a 4 ) ]B( xx)y M x ( a 2 8) ( B ( yy)y + B ) ( yz)z (27) The force on the core, equation (4), can be ealuated in a similar manner. The components of equation (4) are T y Mx B z + ( l 2 24)B ( zx)x + ( a 2 8)B ( zy)y F x M x B xx d (28) + ( a 2 8)B ( zz)z + M x ( a 2 4) ( l 2 12) B ( xx)z Continuing to equation (13), (21) F y M x B xy d F z M x B xz d (29) (30) B y d B y + ( l 2 24)B ( yx)x + ( a 2 8)B ( yy)y and + ( a 2 8)B ( yz)z ( B xyx B xxy ) d [ ( l 2 12) ( a 2 4) ]B ( xx)y Substituting into equation (13) results in (22) (23) Expanding the integral of equation (28) results in B xx d ( B xx + B ( xx)y y + B ( xx)z z + B ( xx)x x d Since the integrals containing first-order terms in and z are zero, equation (31) reduces to F x M x B xx d M x B xx (31) x, y, (32) Ealuating equations (29) and (30) results in T z M x B y + ( l 2 24)B ( yx)x + ( a 2 8)B ( yy)y F y M x B xy (33) + ( a 2 8)B ( yz)z + M x [( l 2 12) and ( a 2 4) ]B ( xx)y (24) F z M x B xz (34) 5

8 Magnetization Perpendicular to Axis of Symmetry For orientation of the magnetization ector perpendicular to the axis of symmetry (x-axis), the only nonzero term in M is M z. Equation (5) becomes B xz [ B ]M M z B yz B zz and the second part of equation (3) becomes ( B yzz + B zzy) r ( M ) B M z ( B xzz B zzx) ( B xzy + B yzx) The first part of equation (3) becomes M B M z B y B x 0 Substituting (36) and (37) into (3) results in T x M z B y d + M z ( B zzy B yzz ) d T y M z B x d + M z ( B xzz B zzx ) d T z 0 + Mz ( B yzx B xzy ) d (35) (36) (37) (38) (39) (40) Ealuating the integrals as before and collecting terms results in T x M z B y M z ( l 2 24)B ( xx)y M z ( a 2 8) ( B ( yy)y + B ( yz)z ) T y M z B x + M z [( 3a 2 8) ( l 2 12) ]B ( xz)z + M z ( l 2 24)B ( xx)x + M z ( a 2 8)B ( xy)y T z M z [( l 2 12) ( a 2 4) ]B ( xy)z (41) (42) (43) The components of the force (eq. (4)) using equation (35) becomes Discussion of Results (44) (45) (46) Examination of equations (25) to (27), (32) to (34), (41) to (43), and (44) to (46) reeals that, for expansion of the applied fields up to second-order terms, no coupling exists between force and torque components. As stated earlier, it is generally assumed that the higher order torque terms, which are functions of second order gradients, can be neglected. Howeer, equation (43) indicates that for magnetization perpendicular to the axis of symmetry, torque about the axis of magnetization can be generated by controlling a higher order term directly, thus allowing the core to be controlled in six degrees of freedom. For a cylindrical permanent magnet core magnetized along the axis of symmetry, from equation (25), only fie-degree-of-freedom control is possible. Concluding Remarks This paper has deeloped the expanded equations for torque and force on a cylindrical permanent magnet core in a large-gap magnetic suspension system. The core was assumed to be uniformly magnetized, and equations were deeloped for two orientations of the magnetization ector. One orientation was parallel to the axis of symmetry of the core and the other was perpendicular to this axis. It is generally assumed that terms that are a function of second-order or higher gradients of the applied fields can be neglected. In practical applications inoling large-gap magnetic suspension systems, these assumptions hae proen to be alid. Howeer, in the case where the magnetization ector is perpendicular to the axis of symmetry of the core, the expanded equations indicate that torque about the magnetization ector can be produced by controlling a second-order gradient directly. This case allows the core to be controlled in six degrees of freedom whereas a cylindrical permanent magnet core magnetized along its axis of symmetry can be controlled only in fie degrees of freedom. NASA Langley Research Center Hampton, VA October 24, 1996 F x M z B xz d F y M z B yz d M z B xz M z B yz F z M z B zz d M z B zz 6

9 Appendix A Expansion of Fields and Gradients About the Nominal Operating Point of a Cylindrical Permanent Magnet Core In appendix A the fields and gradients produced by the suspension system electromagnets are expanded by using a Taylor series about the initial suspension point of the permanent magnet core. The assumption is made that the fields can be adequately described by expansion up to second-order terms. Figure 1 shows the cylindrical core and core coordinate system. The core coordinate system consists of a set of orthogonal x, y, z body-fixed axes that define the motion of the core with respect to an orthogonal x, y, z system fixed in inertial space. The core coordinate system is initially aligned with the x, y, z system. The transformation from inertial coordinates to core coordinates is gien by x y z [ T m ] (A1) where [ T m ] is the orthogonal transformation matrix for a 3, 2, 1 (z, y, x) Euler rotation sequence and is defined as x y z x where r y and is the gradient operator. Using z compact notation, each element of B in equation (A3) can be written as where and 2 B i r 2 B i B i B i r ( 1 2)r T 2 B i r r r 2 B i r ( B i x) x ( B i y) x ( B i z) x B i x B i y B i z ( B i x) y ( B i y) y ( B i z) y ( B i x) z ( B i y) z ( B i z) z (A4) (A5) (A6) Using the notation f ij f i / j and f ( ij)k ( f i / j)/ k, equations (A5) and (A6) can be written as [ T m ] cθ z cθ y sθ z cθ y sθ y ( cθ z sθ y sθ x sθ z cθ x ) ( sθ z sθ y sθ x + cθ z cθ x ) cθ y sθ x ( cθ z sθ y cθ x + sθ z sθ x ) ( sθ z sθ y cθ x cθ z sθ x ) cθ y cθ x and B i B r ix B iy B iz (A7) (A2) where sin has been shortened to s, cos has been shortened to c, and θ z, θ y, and θ x are angles of rotation about the z-, y-, and x-axes, respectiely. The field B and gradients of B produced by the suspension system electromagnets, which are fixed in the inertial frame, are calculated at the origin of the x, y, z system. Expanding B about the origin of the x, y, z system as a Taylor series, up to second order, results in The first-order gradients of B can be written as where 2 B i r 2 B ( ix)x B ( ix)y B ( ix)z B ( iy)x B ( iy)y B ( iy)z B ( iz)x B ( iz)y B ( iz)z ( B i / j) B ij B ij r r (A8) (A9) B B + ( r )B + ( 1 2) ( r ) 2 B (A3) ( B i / j) B r ( ij)x B ( ij)y B ( ij)z (A10) 7

10 The expanded fields can be expressed in core coordinates as Transforming back into core coordinates, B B + ( r )B + ( 1 2) ( r ) 2 B (A11) B x B x [ T r m ] T r ( 1 /2 )r T [ T m ] 2 B x r 2 [ T m ] T r x where r y is the displacement in core coordinates, z B [ T m ]B, and [ T m ]. Since, from equation (A1), r [ T m ] T r (A12) each element of B can be expanded in inertial coordinates by substituting equation (A12) into equation (A4), B i B i B i [ T r m ] T r ( 1 2)r T [ T m ] 2 B i r 2 [ T m ] T r (A13) B [ T m ] B y B y [ T r m ] T r ( 1 /2 )r T [ T m ] 2 B y r 2 [ T m ] T r B z B z [ T r m ] T r ( 1 /2 )r T [ T m ] 2 B z r 2 [ T m ] T r (A14) The expansion of equation (A14) can be simplified by using small-angle assumptions (ref. 6). Under smallangle assumptions, cosθ 1, sinθ θ, and products of angles are neglected. The transformation matrix [ T m ] then becomes [ T m ] 1 θ z θ y θ z 1 θ x θ y θ x 1 (A15) 8

11 Appendix B Torques and Forces on a Magnetic Dipole With Incremental Volume The torques and forces on a magnetic dipole in a steady magnetic field are identical to those on an infinitesimal current loop with the same magnetic moment (ref. 9). Therefore, the equations for torque and force on an infinitesimal current loop will be deeloped first by using the fundamental relationship for the force on a current-carrying-conductor element in a uniform, steady magnetic field. For a discussion of magnetic dipoles and infinitesimal current loops, see references 9 and 10. Infinitesimal Current Loop Consider a plane loop of conductor with steady current I located in the external, uniform, steady magnetic field B (fig. B1). In this region B B 0. The force on an element dl of the conductor is gien by the fundamental relationship (obtained from the Lorentz force law) df I dl B (B1) where df is a ector indicating magnitude and direction of force on the conductor element; I is the scalar magnitude of the current in the conductor element; dl is a ector whose magnitude equals the length of the conductor element and whose direction is in the positie direction of the current; and B is a ector indicating magnitude and direction of the flux density of the external field component. The torque on the loop can be written as [ T I r ( dl B ) ] (B2) where r is the position ector of dl and the integration is around the loop. By using the identity r ( dl B ) dl( r B) B( r dl) equation (B2) can be written as T I ( r B)dl B r dl ] [ (B3) (B4) Using Stokes s theorem and a related result (ref. 9, p. 289), the line integrals in equation (B4) can be transformed into surface integrals resulting in T I [ da ( r B) ] B ( r ) da s s (B5) where da is a ector whose magnitude is a differential area and whose direction is normal to the plane of the current loop in the sense of the right-hand rule relatie to the direction of current flow, is the gradient operator, and the integrals are oer the surface that is defined by the conductor loop. Since r is zero and ( r B ) B for constant B, equation (B5) simplifies to Taking the integral results in (B6) (B7) An infinitesimal current loop can be defined by letting A go toward zero and I go toward infinity, keeping the product IA finite. For an infinitesimal current loop, the requirement that B be uniform no longer exists. The product IA is called the magnetic moment of the loop and is designated by the letter m. Therefore equation (B7) becomes (B8) The torque T acts on the infinitesimal current loop in a direction to align the magnetic moment m with the external field B. If m and B are misaligned by the angle θ, the magnitude of the torque is (B9) To increase θ by the amount d θ, work dw must be done against the torque T resulting in an increase in potential energy du: (B10) The potential energy of an infinitesimal current loop in an external magnetic field can then be obtained by integrating equation (B10): (B11) where the constant of integration is chosen to be zero when m is perpendicular to B. The force on the infinitesimal current loop can be obtained from equation (B11). If an external force F displaces the infinitesimal current loop by the infinitesimal distance dr, then the work done dw will be equal to a decrease in potential energy, du: Therefore T I ( da B) s T IA B T m B T mbsinθ du dw T dθ mbsinθdθ U mbcosθ m B dw F dr du U dr F U ( m B) (B12) (B13) 9

12 The right-hand side of equation (B13) can be expanded as ( m B) m ( B) + ( m )B + B ( m) + ( B )m (B14) Since ( B), ( B), and ( m) are zero, equation (B13) can be written in the form F ( m )B (B15) This form is generally used in the deelopment of the equations for large-gap magnetic suspension systems (refs. 1 through 6). Magnetic Dipole With Incremental Volume The magnetic moment of a permanent magnet dipole with north and south poles separated by length l and with pole strength Q m is defined as m Q m l (B16) The magnetic moment m is a ector pointing from the south pole to the north pole. In the case of an actual magnet, Q m and l may be indefinite but m can be determined and is sufficient to specify the fields of the magnet at a large distance from it. At large distances, a magnetic dipole with magnetic moment Q m l can be treated the same as an infinitesimal current loop with magnetic moment IA and is identical in effect if Q m l IA. Therefore, in a steady magnetic field B the equations for torques and forces on a magnetic dipole with magnetic moment m are the same as equations (B8) and (B15). In theory, it can be assumed that a permanent magnet of a gien olume consists of a large number of uniformly distributed permanent magnet dipoles with incremental olumes δ which are oriented in the same direction. The magnetic moment δm of a gien dipole with incremental olume δ can be coneniently described by a quantity called the magnetization M, which is defined as the magnetic moment per unit olume. That is, (B17) The total magnetic moment m for a gien permanent magnet can then be written as (B18) where the integration is oer the olume of the permanent magnet. Magnetization is also a ector and has the same direction as m. If the permanent magnet is uniformly magnetized, that is, M is constant oer the olume of the permanent magnet, then (B19) For a discussion of magnetic dipoles and magnetization, see reference 10. The torques and forces on an incremental olume of permanent magnet material, in terms of the magnetization M, can then be written as and M m m δm δ from equations (B8), (B15), and (B17). M d M δt ( M B)δ δf ( M )Bδ (B20) (B21) 10

13 z I dl r B y x Figure B1. Plane loop of conductor with steady current I in uniform steady magnetic field B. 11

14 References 1. Basmajian, V. V.; Copeland, A. B.; and Stephens, T.: Studies Related to the Design of a Magnetic Suspension and Balance System. NASA CR-66233, Stephens, Timothy: Design, Construction, and Ealuation of a Magnetic Suspension and Balance System for Wind Tunnels. NASA CR-66903, Coert, Eugene E.; Finston, Morton; Vlajinac, Milan; and Stephens, Timothy: Magnetic Balance and Suspension Systems for Use in Wind Tunnels. Progress in Aerospace Sciences, Volume 14, D. Küchemann, P. Carriére, B. Etkin, W. Fiszdon, N. Rott, J. Smolderen, I. Tani, and W. Wuest, eds., Pergamon Press, 1973, pp Britcher, C. P.: Some Aspects of Wind Tunnel Magnetic Suspension Systems With Special Applications at Large Physical Scales. NASA CR , Groom, Nelson J.: Analytical Model of a Fie Degree of Freedom Magnetic Suspension and Positioning System. NASA TM , Groom, Nelson J.; and Britcher, Colin P.: Open-Loop Characteristics of Magnetic Suspension Systems Using Electromagnets Mounted in a Planar Array. NASA TP-3229, Goodyer, M. J.: Roll Control Techniques on Magnetic Suspension Systems. Aeronaut. Q., ol. xiii, Feb. 1967, pt. 1, pp Groom, Nelson J.: Description of the Large Gap Magnetic Suspension System (LGMSS) Ground-Based Experiment. Technology 2000, NASA CP-3109, ol. 2, 1991, pp Clemmow, P. C.: An Introduction to Electromagnetic Theory. Cambridge Uni. Press, Kraus, John D.: Electromagnetics. McGraw-Hill Book Co., Inc.,

15

16 Form Approed OMB No Public reporting burden for this collection of information is estimated to aerage 1 hour per response, including the time for reiewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reiewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Serices, Directorate for Information Operations and Reports, 1215 Jefferson Dais Highway, Suite 1204, Arlington, VA , and to the Office of Management and Budget, Paperwork Reduction Project ( ), Washington, DC AGENCY USE ONLY (Leae blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 6. AUTHOR(S) REPORT DOCUMENTATION PAGE February 1997 Technical Paper Expanded Equations for Torque and Force on a Cylindrical Permanent Magnet Core in a Large-Gap Magnetic Suspension System Nelson J. Groom WU PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) NASA Langley Research Center Hampton, VA PERFORMING ORGANIZATION REPORT NUMBER L SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) National Aeronautics and Space Administration Washington, DC SPONSORING/MONITORING AGENCY REPORT NUMBER NASA TP SUPPLEMENTARY NOTES 12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE Unclassified Unlimited Subject Category 31 Aailability: NASA CASI (301) ABSTRACT (Maximum 200 words) The expanded equations for torque and force on a cylindrical permanent magnet core in a large-gap magnetic suspension system are presented. The core is assumed to be uniformly magnetized, and equations are deeloped for two orientations of the magnetization ector. One orientation is parallel to the axis of symmetry, and the other is perpendicular to this axis. Fields and gradients produced by suspension system electromagnets are assumed to be calculated at a point in inertial space which coincides with the origin of the core axis system in its initial alignment. Fields at a gien point in the core are defined by expanding the fields produced at the origin as a Taylor series. The assumption is made that the fields can be adequately defined by expansion up to second-order terms. Examination of the expanded equations for the case where the magnetization ector is perpendicular to the axis of symmetry reeals that some of the second-order gradient terms proide a method of generating torque about the axis of magnetization and therefore proide the ability to produce six-degree-of-freedom control. 14. SUBJECT TERMS Magnetic suspension; Large-gap magnetic suspension; Magnetic leitation; Magnetic suspension model 17. SECURITY CLASSIFICATION OF REPORT 18. SECURITY CLASSIFICATION OF THIS PAGE 19. SECURITY CLASSIFICATION OF ABSTRACT Unclassified Unclassified Unclassified 15. NUMBER OF PAGES 16. PRICE CODE 20. LIMITATION OF ABSTRACT NSN Standard Form 298 (Re. 2-89) Prescribed by ANSI Std. Z A03

Simplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System

Simplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System NASA Technical Memorandum 112868 Simplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System Nelson J. Groom Langley Research Center, Hampton, Virginia June 1997 National

More information

IMPLEMENTATION OF A DECOUPLED CONTROLLER ELECTROMAGNETS MOUNTED IN A PLANAR ARRAY. NASA Langley Research Center. Hampton, VA.

IMPLEMENTATION OF A DECOUPLED CONTROLLER ELECTROMAGNETS MOUNTED IN A PLANAR ARRAY. NASA Langley Research Center. Hampton, VA. IMPLEMENTATION OF A DECOUPLED CONTROLLER FOR A MAGNETIC SUSPENSION SYSTEM USING ELECTROMAGNETS MOUNTED IN A PLANAR ARRAY D. E. Cox N. J. Groom NASA Langley Research Center Hampton, VA Abstract An implementation

More information

Volume 6 Water Surface Profiles

Volume 6 Water Surface Profiles A United States Contribution to the International Hydrological Decade HEC-IHD-0600 Hydrologic Engineering Methods For Water Resources Development Volume 6 Water Surface Profiles July 1975 Approved for

More information

Demonstration of the Coupled Evolution Rules 163 APPENDIX F: DEMONSTRATION OF THE COUPLED EVOLUTION RULES

Demonstration of the Coupled Evolution Rules 163 APPENDIX F: DEMONSTRATION OF THE COUPLED EVOLUTION RULES Demonstration of the Coupled Evolution Rules 163 APPENDIX F: DEMONSTRATION OF THE COUPLED EVOLUTION RULES Before going into the demonstration we need to point out two limitations: a. It assumes I=1/2 for

More information

( ) R kj. = y k y j. y A ( ) z A. y a. z a. Derivatives of the second order electrostatic tensor with respect to the translation of ( ) δ yβ.

( ) R kj. = y k y j. y A ( ) z A. y a. z a. Derivatives of the second order electrostatic tensor with respect to the translation of ( ) δ yβ. Supporting information Derivatives of R with respect to the translation of fragment along the y and z axis: y = y k y j (S1) z ( = z z k j) (S2) Derivatives of S with respect to the translation of fragment

More information

Electromagnetic Scattering From Arbitrarily Shaped Aperture Backed by Rectangular Cavity Recessed in Infinite Ground Plane

Electromagnetic Scattering From Arbitrarily Shaped Aperture Backed by Rectangular Cavity Recessed in Infinite Ground Plane NASA Technical Paper 3640 Electromagnetic Scattering From Arbitrarily Shaped Aperture Backed by Rectangular Cavity Recessed in Infinite Ground Plane C. R. Cockrell and Fred B. Beck Langley Research Center

More information

Closed-Form Evaluation of Mutual Coupling in a Planar Array of Circular Apertures

Closed-Form Evaluation of Mutual Coupling in a Planar Array of Circular Apertures NASA Technical Paper 3552 Closed-Form Evaluation of Mutual Coupling in a Planar Array of Circular Apertures M. C. Bailey Langley Research Center Hampton, Virginia National Aeronautics and Space Administration

More information

EJL R Sβ. sum. General objective: Reduce the complexity of the analysis by exploiting symmetry. Garth J. Simpson

EJL R Sβ. sum. General objective: Reduce the complexity of the analysis by exploiting symmetry. Garth J. Simpson e sum = EJL R Sβ General objective: Reduce the complexity of the analysis by exploiting symmetry. Specific Objectives: 1. The molecular symmetry matrix S. How to populate it.. Relationships between the

More information

2x (x 2 + y 2 + 1) 2 2y. (x 2 + y 2 + 1) 4. 4xy. (1, 1)(x 1) + (1, 1)(y + 1) (1, 1)(x 1)(y + 1) 81 x y y + 7.

2x (x 2 + y 2 + 1) 2 2y. (x 2 + y 2 + 1) 4. 4xy. (1, 1)(x 1) + (1, 1)(y + 1) (1, 1)(x 1)(y + 1) 81 x y y + 7. Homework 8 Solutions, November 007. (1 We calculate some derivatives: f x = f y = x (x + y + 1 y (x + y + 1 x = (x + y + 1 4x (x + y + 1 4 y = (x + y + 1 4y (x + y + 1 4 x y = 4xy (x + y + 1 4 Substituting

More information

Application of Model Based Parameter Estimation for RCS Frequency Response Calculations Using Method of Moments

Application of Model Based Parameter Estimation for RCS Frequency Response Calculations Using Method of Moments NASA/CR-1998-206951 Application of Model Based Parameter Estimation for RCS Frequency Response Calculations Using Method of Moments C. J. Reddy Hampton University, Hampton, Virginia March 1998 The NASA

More information

Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.

Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J. Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics

More information

The Bernoulli equation in a moving reference frame

The Bernoulli equation in a moving reference frame IOP PUBLISHING Eur. J. Phys. 3 (011) 517 50 EUROPEAN JOURNAL OF PHYSICS doi:10.1088/0143-0807/3//0 The Bernoulli equation in a moving reference frame Carl E Mungan Physics Department, US Naval Academy

More information

Characterization of Multilayer Piezoelectric Actuators for Use in Active Isolation Mounts

Characterization of Multilayer Piezoelectric Actuators for Use in Active Isolation Mounts NASA Technical Memorandum 4742 Characterization of Multilayer Piezoelectric Actuators for Use in Active Isolation Mounts Stephanie A. Wise Langley Research Center Hampton, Virginia Matthew W. Hooker Science

More information

Dynamics of Transport and Variability in the Denmark Strait Overflow

Dynamics of Transport and Variability in the Denmark Strait Overflow Approved for public release; distribution is unlimited. Dynamics of Transport and Variability in the Denmark Strait Overflow by James B. Girton Technical Report APL-UW TR 0103 August 2001 Applied Physics

More information

Report Documentation Page

Report Documentation Page Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

Multiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015

Multiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015 Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction

More information

A Simple Method for Nucleon-Nucleon Cross Sections in a Nucleus

A Simple Method for Nucleon-Nucleon Cross Sections in a Nucleus NASA/TP-1999-209125 A Simple Method for Nucleon-Nucleon Cross Sections in a Nucleus R. K. Tripathi Hampton University, Hampton, Virginia Francis A. Cucinotta Lyndon B. Johnson Space Center, Houston, Texas

More information

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

Multiple Integrals and Vector Calculus: Synopsis

Multiple Integrals and Vector Calculus: Synopsis Multiple Integrals and Vector Calculus: Synopsis Hilary Term 28: 14 lectures. Steve Rawlings. 1. Vectors - recap of basic principles. Things which are (and are not) vectors. Differentiation and integration

More information

Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications

Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications Joseph R. Gabriel Naval Undersea Warfare Center Newport, Rl 02841 Steven M. Kay University of Rhode

More information

Magnetostatics and the vector potential

Magnetostatics and the vector potential Magnetostatics and the vector potential December 8, 2015 1 The divergence of the magnetic field Starting with the general form of the Biot-Savart law, B (x 0 ) we take the divergence of both sides with

More information

REPORT DOCUMENTATION PAGE

REPORT DOCUMENTATION PAGE REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 Public reportinq burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

Lecture 5 Charge Density & Differential Charge. Sections: 2.3, 2.4, 2.5 Homework: See homework file

Lecture 5 Charge Density & Differential Charge. Sections: 2.3, 2.4, 2.5 Homework: See homework file Lecture 5 Charge Density & Differential Charge Sections: 2.3, 2.4, 2.5 Homework: See homework file Point Charge as an Approximation charge occupies a finite olume and may hae arying density a charged body

More information

0 a 3 a 2 a 3 0 a 1 a 2 a 1 0

0 a 3 a 2 a 3 0 a 1 a 2 a 1 0 Chapter Flow kinematics Vector and tensor formulae This introductory section presents a brief account of different definitions of ector and tensor analysis that will be used in the following chapters.

More information

(Refer Slide Time: 1:58 min)

(Refer Slide Time: 1:58 min) Applied Mechanics Prof. R. K. Mittal Department of Applied Mechanics Indian Institution of Technology, Delhi Lecture No. # 13 Moments and Products of Inertia (Contd.) Today s lecture is lecture thirteen

More information

Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method

Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method Manohar D. Deshpande, NASA Langley Research Center, Hmapton, VA., USA Abstract

More information

Microwave Circuit Simulator for MATLAB

Microwave Circuit Simulator for MATLAB Microwave Circuit Simulator for MATLAB Romeo D. del Rosario and Daniel C. Judy ARL-TN-184 March 2002 Approved for public release; distribution unlimited. The findings in this report are not to be construed

More information

Benefits of Force Limiting Vibration Testing

Benefits of Force Limiting Vibration Testing NASA/TM 1999-209382 Benefits of Force Limiting Vibration Testing Mark E. McNelis Glenn Research Center, Cleveland, Ohio Terry D. Scharton Jet Propulsion Laboratory, Pasadena, California August 1999 The

More information

Predicting Deformation and Strain Behavior in Circuit Board Bend Testing

Predicting Deformation and Strain Behavior in Circuit Board Bend Testing Predicting Deformation and Strain Behavior in Circuit Board Bend Testing by Ronen Aniti ARL-CR-0693 May 2012 prepared by Science and Engineering Apprenticeship Program George Washington University 2121

More information

DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION

DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION CETN-I-6 3185 DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION PURPOSE : To present a parameterized model of a directional spectrum of the sea surface using an energy spectrum and a value for the

More information

Spacecraft Dynamics and Control

Spacecraft Dynamics and Control Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 16: Euler s Equations Attitude Dynamics In this Lecture we will cover: The Problem of Attitude Stabilization Actuators Newton

More information

Classical Mechanics. Luis Anchordoqui

Classical Mechanics. Luis Anchordoqui 1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body

More information

Non-Linear Finite Element Modeling of THUNDER Piezoelectric Actuators

Non-Linear Finite Element Modeling of THUNDER Piezoelectric Actuators NASA/TM-1999-209322 ARL-TR-1977 Non-Linear Finite Element Modeling of THUNDER Piezoelectric Actuators Barmac K. Taleghani U. S. Army Research Laboratory Vehicle Technology Directorate Langley Research

More information

Grant Number: N IP To compare obtained theoretical results with NPAL experimental data.

Grant Number: N IP To compare obtained theoretical results with NPAL experimental data. Coherence of Low-Frequency Sound Signals Propagating through a Fluctuating Ocean: Analysis and Theoretical Interpretation of 2004 NPAL Experimental Data Alexander G. Voronovich NOAA/ESRL, PSD4, 325 Broadway,

More information

A GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS

A GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS A GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS Michael Thiel*, Truong Vu Bang Giang and Achim Dreher Institute of

More information

High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems

High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems Aerodynamic Issues of Unmanned Air Vehicles Fluid-Structure Interaction High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems Raymond E. Gordnier Computational Sciences

More information

CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017

CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017 CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and

More information

Section 6: PRISMATIC BEAMS. Beam Theory

Section 6: PRISMATIC BEAMS. Beam Theory Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam

More information

Eulerian Mapping Closure Approach for Probability Density Function of Concentration in Shear Flows

Eulerian Mapping Closure Approach for Probability Density Function of Concentration in Shear Flows NASA/CR-2002-211627 ICASE Report No. 2002-9 Eulerian Mapping Closure Approach for Probability Density Function of Concentration in Shear Flows Guowei He ICASE, Hampton, Virginia May 2002 The NASA STI Program

More information

HIGH PERFORMANCE CONTROLLERS BASED ON REAL PARAMETERS TO ACCOUNT FOR PARAMETER VARIATIONS DUE TO IRON SATURATION

HIGH PERFORMANCE CONTROLLERS BASED ON REAL PARAMETERS TO ACCOUNT FOR PARAMETER VARIATIONS DUE TO IRON SATURATION HIGH PERFORMANCE CONTROLLERS BASED ON REAL PARAMETERS TO ACCOUNT FOR PARAMETER VARIATIONS DUE TO IRON SATURATION Jorge G. Cintron-Rivera, Shanelle N. Foster, Wesley G. Zanardelli and Elias G. Strangas

More information

Manipulator Dynamics 2. Instructor: Jacob Rosen Advanced Robotic - MAE 263D - Department of Mechanical & Aerospace Engineering - UCLA

Manipulator Dynamics 2. Instructor: Jacob Rosen Advanced Robotic - MAE 263D - Department of Mechanical & Aerospace Engineering - UCLA Manipulator Dynamics 2 Forward Dynamics Problem Given: Joint torques and links geometry, mass, inertia, friction Compute: Angular acceleration of the links (solve differential equations) Solution Dynamic

More information

Chapter 27: Magnetic Field and Magnetic Forces

Chapter 27: Magnetic Field and Magnetic Forces Chapter 27: Magnetic Field and Magnetic Forces Iron ore found near Magnesia Compass needles align N-S: magnetic Poles North (South) Poles attracted to geographic North (South) Like Poles repel, Opposites

More information

Phys 7221 Homework # 8

Phys 7221 Homework # 8 Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with

More information

Awell-known lecture demonstration1

Awell-known lecture demonstration1 Acceleration of a Pulled Spool Carl E. Mungan, Physics Department, U.S. Naval Academy, Annapolis, MD 40-506; mungan@usna.edu Awell-known lecture demonstration consists of pulling a spool by the free end

More information

Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects

Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects Chris Harrison Centre for Marine

More information

Ocean Acoustics Turbulence Study

Ocean Acoustics Turbulence Study Ocean Acoustics Turbulence Study PI John Oeschger Coastal Systems Station 6703 West Highway 98 Panama City, FL 32407 phone: (850) 230-7054 fax: (850) 234-4886 email: OeschgerJW@ncsc.navy.mil CO-PI Louis

More information

Lecture 8. Stress Strain in Multi-dimension

Lecture 8. Stress Strain in Multi-dimension Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]

More information

REPORT DOCUMENTATION PAGE

REPORT DOCUMENTATION PAGE REPORT DOCUMENTATION PAGE Form Approved OMB NO. 0704-0188 Public Reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

A Method for Characterizing Thermoset Polyimides

A Method for Characterizing Thermoset Polyimides NASA Technical Memorandum 4707 A Method for Characterizing Thermoset Polyimides C. Ranganathaiah, D. R. Sprinkle, and R. H. Pater Langley Research Center Hampton, Virginia A. Eftekhari Hampton University

More information

Optimal Pitch Thrust-Vector Angle and Benefits for all Flight Regimes

Optimal Pitch Thrust-Vector Angle and Benefits for all Flight Regimes NASA/TM-2000-209021 Optimal Pitch Thrust-Vector Angle and Benefits for all Flight Regimes Glenn B. Gilyard Dryden Flight Research Center Edwards, California Alexander Bolonkin Senior Research Associate

More information

Stress, Strain, Mohr s Circle

Stress, Strain, Mohr s Circle Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected

More information

UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS

UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS SUBAREA I. COMPETENCY 1.0 UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS MECHANICS Skill 1.1 Calculating displacement, aerage elocity, instantaneous elocity, and acceleration in a gien frame of reference

More information

The NASA STI Program Office... in Profile

The NASA STI Program Office... in Profile NASA/CR-1998-206950 Application of Model Based Parameter Estimation for Fast Frequency Response Calculations of Input Characteristics of Cavity-Backed Aperture Antennas Using Hybrid FEM/MoM Technique C.

More information

REPORT DOCUMENTATION PAGE

REPORT DOCUMENTATION PAGE REPORT DOCUMENTATION PAGE Form Approved OMB NO. 0704-0188 Public Reporting Burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures

Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures AFOSR grant #FA9550-10-1-0061 Program manager: Dr. Les Lee PI: C.T. Sun Purdue University West Lafayette, Indiana

More information

CRS Report for Congress

CRS Report for Congress CRS Report for Congress Received through the CRS Web Order Code RS21396 Updated May 26, 2006 Summary Iraq: Map Sources Hannah Fischer Information Research Specialist Knowledge Services Group This report

More information

arxiv:physics/ v1 [physics.class-ph] 10 Nov 2006

arxiv:physics/ v1 [physics.class-ph] 10 Nov 2006 Derivation of magnetic Coulomb s law for thin, semi-infinite solenoids Masao Kitano Department of Electronic Science and Engineering, Kyoto University, Katsura, Kyoto 615-8510, Japan and CREST, Japan Science

More information

FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES

FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES 1 FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES Robert L. Street Environmental Fluid Mechanics Laboratory Department of Civil Engineering Stanford University Stanford, CA 94305-4020 650-723-4969;

More information

VLBA IMAGING OF SOURCES AT 24 AND 43 GHZ

VLBA IMAGING OF SOURCES AT 24 AND 43 GHZ VLBA IMAGING OF SOURCES AT 24 AND 43 GHZ D.A. BOBOLTZ 1, A.L. FEY 1, P. CHARLOT 2,3 & THE K-Q VLBI SURVEY COLLABORATION 1 U.S. Naval Observatory 3450 Massachusetts Ave., NW, Washington, DC, 20392-5420,

More information

Diagonal Representation of Certain Matrices

Diagonal Representation of Certain Matrices Diagonal Representation of Certain Matrices Mark Tygert Research Report YALEU/DCS/RR-33 December 2, 2004 Abstract An explicit expression is provided for the characteristic polynomial of a matrix M of the

More information

Real-Time Environmental Information Network and Analysis System (REINAS)

Real-Time Environmental Information Network and Analysis System (REINAS) Calhoun: The NPS Institutional Archive Faculty and Researcher Publications Faculty and Researcher Publications 1998-09 Real-Time Environmental Information Network and Analysis System (REINAS) Nuss, Wendell

More information

Parameterizing the Effects of Upper-Ocean Large Eddies on Air-Sea Interaction

Parameterizing the Effects of Upper-Ocean Large Eddies on Air-Sea Interaction Parameterizing the Effects of Upper-Ocean Large Eddies on Air-Sea Interaction Ming Li Horn Point Laboratory University of Maryland Center for Environmental Science 2020 Horn Point Road, Cambridge, MD 21613

More information

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models PI Isaac Ginis Graduate School of Oceanography

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

Analysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E.

Analysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E. Concept V- Hull Floor Configuration in Two Dimensions Dr. Bijan Khatib-Shahidi & Rob E. Smith 10 November 2010 : Dist A. Approved for public release Report Documentation Page Form Approved OMB No. 0704-0188

More information

Thin Film Heat Flux Sensor of Improved Design

Thin Film Heat Flux Sensor of Improved Design NASA/TM -11566 Thin Film Heat Flux Sensor of Improved Design Gus Fralick and John Wrbanek Glenn Research Center, Cleveland, Ohio Charles Blaha Akima Corporation, Cleveland, Ohio September The NASA STI

More information

Estimating Geometric Aspects of Relative Satellite Motion Using Angles-Only Measurements

Estimating Geometric Aspects of Relative Satellite Motion Using Angles-Only Measurements Estimating Geometric Aspects of Relative Satellite Motion Using Angles-Only Measurements Jason Schmidt 1 Dept of Mechanical & Aerospace Engineering, Utah State University, Logan, UT, 84322 Thomas Alan

More information

Rotational & Rigid-Body Mechanics. Lectures 3+4

Rotational & Rigid-Body Mechanics. Lectures 3+4 Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions

More information

A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any

A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»

More information

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models PI Isaac Ginis Graduate School of Oceanography

More information

Review Sheet for the Final

Review Sheet for the Final Review Sheet for the Final Math 6-4 4 These problems are provided to help you study. The presence of a problem on this handout does not imply that there will be a similar problem on the test. And the absence

More information

Multipole moments. Dipole moment. The second moment µ is more commonly called the dipole moment, of the charge. distribution and is a vector

Multipole moments. Dipole moment. The second moment µ is more commonly called the dipole moment, of the charge. distribution and is a vector Dipole moment Multipole moments The second moment µ is more commonly called the dipole moment, of the charge distribution and is a vector µ = µ x ˆx + µ y ŷ + µ z ẑ where the α component is given by µ

More information

REPORT DOCUMENTATION PAGE

REPORT DOCUMENTATION PAGE REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

Scattering of Internal Gravity Waves at Finite Topography

Scattering of Internal Gravity Waves at Finite Topography Scattering of Internal Gravity Waves at Finite Topography Peter Muller University of Hawaii Department of Oceanography 1000 Pope Road, MSB 429 Honolulu, HI 96822 phone: (808)956-8081 fax: (808)956-9164

More information

An Observational and Modeling Study of Air-Sea Fluxes at Very High Wind Speeds

An Observational and Modeling Study of Air-Sea Fluxes at Very High Wind Speeds An Observational and Modeling Study of Air-Sea Fluxes at Very High Wind Speeds Kerry Emanuel Room 54-1620, MIT 77 Massachusetts Avenue Cambridge, MA 02139 phone: (617) 253-2462 fax: (425) 740-9133 email:

More information

21.60 Worksheet 8 - preparation problems - question 1:

21.60 Worksheet 8 - preparation problems - question 1: Dynamics 190 1.60 Worksheet 8 - preparation problems - question 1: A particle of mass m moes under the influence of a conseratie central force F (r) =g(r)r where r = xˆx + yŷ + zẑ and r = x + y + z. A.

More information

11.1 Three-Dimensional Coordinate System

11.1 Three-Dimensional Coordinate System 11.1 Three-Dimensional Coordinate System In three dimensions, a point has three coordinates: (x,y,z). The normal orientation of the x, y, and z-axes is shown below. The three axes divide the region into

More information

Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO

Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO William L. Dimpfl Space Science Applications Laboratory The Aerospace Corporation

More information

Testing for Random Limit Load Versus Static Limit Load

Testing for Random Limit Load Versus Static Limit Load NASA Technical Memorandum 108542 Testing for Random Limit Load Versus Static Limit Load H.M. Lee September 1997 NASA Technical Memorandum 108542 Testing for Random Limit Load Versus Static Limit Load H.M.

More information

REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT)

REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT) AL/EQ-TP-1993-0307 REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT) John T. Mourand, John C. Crittenden, David W. Hand, David L. Perram, Sawang Notthakun Department of Chemical Engineering

More information

7a3 2. (c) πa 3 (d) πa 3 (e) πa3

7a3 2. (c) πa 3 (d) πa 3 (e) πa3 1.(6pts) Find the integral x, y, z d S where H is the part of the upper hemisphere of H x 2 + y 2 + z 2 = a 2 above the plane z = a and the normal points up. ( 2 π ) Useful Facts: cos = 1 and ds = ±a sin

More information

ICASE ON THE REMOVAL OF BOUNDARY ERRORS CAUSED BY RUNGE-KUTTA INTEGRATION OF NON-LINEAR PARTIAL DIFFERENTIAL EQUATIONS

ICASE ON THE REMOVAL OF BOUNDARY ERRORS CAUSED BY RUNGE-KUTTA INTEGRATION OF NON-LINEAR PARTIAL DIFFERENTIAL EQUATIONS DTIC pi FCTF NASA Contractor Report 194989 %-- Q^Q g7 1994 1 ICASE Report No. 94-79 ICASE ON THE REMOVAL OF BOUNDARY ERRORS CAUSED BY RUNGE-KUTTA INTEGRATION OF NON-LINEAR PARTIAL DIFFERENTIAL EQUATIONS

More information

Guide to the Development of a Deterioration Rate Curve Using Condition State Inspection Data

Guide to the Development of a Deterioration Rate Curve Using Condition State Inspection Data Guide to the Development of a Deterioration Rate Curve Using Condition State Inspection Data by Guillermo A. Riveros and Elias Arredondo PURPOSE: The deterioration of elements of steel hydraulic structures

More information

P. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE

P. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE A wavelet-based generalization of the multifractal formalism from scalar to vector valued d- dimensional random fields : from theoretical concepts to experimental applications P. Kestener and A. Arneodo

More information

Levinson s Theorem and the Nonlocal Saito Potential

Levinson s Theorem and the Nonlocal Saito Potential Journal of Mathematics Research February, 21 Levinson s Theorem and the Nonlocal Saito Potential S.B. Qadri Department of Physics, Ohio State University Columbus, Ohio 4321, USA Naval Research Laboratory,

More information

^^lo^j - L^ REPORT DOCUMENTATION PAGE. U. S. Army Research Office P.O. Box Research Triangle Park, NC

^^lo^j - L^ REPORT DOCUMENTATION PAGE. U. S. Army Research Office P.O. Box Research Triangle Park, NC REPORT DOCUMENTATION PAGE Fomi Approved 0MB NO. 0704-0188 Public Reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

SW06 Shallow Water Acoustics Experiment Data Analysis

SW06 Shallow Water Acoustics Experiment Data Analysis DISTRIBUTION STATEMENT A: Approved for public release; distribution is unlimited. SW06 Shallow Water Acoustics Experiment Data Analysis James F. Lynch MS #12, Woods Hole Oceanographic Institution Woods

More information

Testing Turbulence Closure Models Against Oceanic Turbulence Measurements

Testing Turbulence Closure Models Against Oceanic Turbulence Measurements Testing Turbulence Closure Models Against Oceanic Turbulence Measurements J. H. Trowbridge Woods Hole Oceanographic Institution Woods Hole, MA 02543 phone: 508-289-2296 fax: 508-457-2194 e-mail: jtrowbridge@whoi.edu

More information

Monte Carlo Simulation of Endlinking Oligomers

Monte Carlo Simulation of Endlinking Oligomers NASA/TM-1998-207649 Monte Carlo Simulation of Endlinking Oligomers Jeffrey A. Hinkley Langley Research Center, Hampton, Virginia Jennifer A. Young University of Virginia, Charlottesville, Virginia April

More information

Babylonian resistor networks

Babylonian resistor networks IOP PUBLISHING Eur. J. Phys. 33 (2012) 531 537 EUROPEAN JOURNAL OF PHYSICS doi:10.1088/0143-0807/33/3/531 Babylonian resistor networks Carl E Mungan 1 and Trevor C Lipscombe 2 1 Physics Department, US

More information

Attribution Concepts for Sub-meter Resolution Ground Physics Models

Attribution Concepts for Sub-meter Resolution Ground Physics Models Attribution Concepts for Sub-meter Resolution Ground Physics Models 76 th MORS Symposium US Coast Guard Academy Approved for public release distribution. 2 Report Documentation Page Form Approved OMB No.

More information

Metrology Experiment for Engineering Students: Platinum Resistance Temperature Detector

Metrology Experiment for Engineering Students: Platinum Resistance Temperature Detector Session 1359 Metrology Experiment for Engineering Students: Platinum Resistance Temperature Detector Svetlana Avramov-Zamurovic, Carl Wick, Robert DeMoyer United States Naval Academy Abstract This paper

More information

A report (dated September 20, 2011) on. scientific research carried out under Grant: FA

A report (dated September 20, 2011) on. scientific research carried out under Grant: FA A report (dated September 20, 2011) on scientific research carried out under Grant: FA2386-10-1-4150 First-principles determination of thermal properties in nano-structured hexagonal solids with doping

More information

LISA: THE LASER INTERFEROMETER SPACE ANTENNA

LISA: THE LASER INTERFEROMETER SPACE ANTENNA LISA: THE LASER INTERFEROMETER SPACE ANTENNA Masimo Tinto Jet Propulsion Laboratory, California Institute of Technology Abstract The Laser Interferometer Space Antenna (LISA) is a deep-space mission, jointly

More information

Mine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST)

Mine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST) Mine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST) Marcelo H. Garcia Department of Civil and Environmental Engineering University of Illinois at Urbana-Champaign 205 North Mathews

More information

Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)

Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS) Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS) D.N. Kurk, T.E. Beegle, S.C. Spence and R.J. Swatski Report Documentation

More information

INFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS

INFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS INFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS AMOS 2005 TECHNICAL CONFERENCE WORKSHOP 5 September, 2005 Maui, Hawaii M. Venner AFRL, Edwards AFB, CA M. Braunstein, L. Bernstein Spectral Sciences,

More information

Understanding Near-Surface and In-Cloud Turbulent Fluxes in the Coastal Stratocumulus-Topped Boundary Layers

Understanding Near-Surface and In-Cloud Turbulent Fluxes in the Coastal Stratocumulus-Topped Boundary Layers Understanding Near-Surface and In-Cloud Turbulent Fluxes in the Coastal Stratocumulus-Topped Boundary Layers Qing Wang Meteorology Department, Naval Postgraduate School Monterey, CA 93943 Phone: (831)

More information

Mixture Distributions for Modeling Lead Time Demand in Coordinated Supply Chains. Barry Cobb. Alan Johnson

Mixture Distributions for Modeling Lead Time Demand in Coordinated Supply Chains. Barry Cobb. Alan Johnson Mixture Distributions for Modeling Lead Time Demand in Coordinated Supply Chains Barry Cobb Virginia Military Institute Alan Johnson Air Force Institute of Technology AFCEA Acquisition Research Symposium

More information

Report Documentation Page

Report Documentation Page Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information