Riemannian Geometry

Total Page:16

File Type:pdf, Size:1020Kb

Riemannian Geometry Riemannian Geometry it is a draft of Lecture Notes of H.M. Khudaverdian. Manchester, 20-th May 2011 Contents 1 Riemannian manifolds 4 1.1 Manifolds. Tensors. (Recalling) . 4 1.2 Riemannian manifold|manifold equipped with Riemannian metric . 7 1.2.1 ∗ Pseudoriemannian manifold . 11 1.3 Scalar product. Length of tangent vectors and angle between vectors. Length of the curve . 11 1.3.1 Length of the curve . 12 1.4 Riemannian structure on the surfaces embedded in Euclidean space . 15 1.4.1 Internal and external coordinates of tangent vector . 15 1.4.2 Explicit formulae for induced Riemannian metric (First Quadratic form) . 17 1.4.3 Induced Riemannian metrics. Examples. 20 1.4.4 ∗Induced metric on two-sheeted hyperboloid embedded in pseudo-Euclidean space. 27 1.5 Isometries of Riemanian manifolds. 28 1.5.1 Examples of local isometries . 29 1.6 Volume element in Riemannian manifold . 31 1.6.1 Volume of parallelepiped . 31 1.6.2 Invariance of volume element under changing of coor- dinates . 33 1.6.3 Examples of calculating volume element . 34 1 2 Covariant differentiaion. Connection. Levi Civita Connec- tion on Riemannian manifold 36 2.1 Differentiation of vector field along the vector field.—Affine connection . 36 2.1.1 Definition of connection. Christoffel symbols of con- nection . 37 2.1.2 Transformation of Christoffel symbols for an arbitrary connection . 39 2.1.3 Canonical flat affine connection . 40 2.1.4 ∗ Global aspects of existence of connection . 43 2.2 Connection induced on the surfaces . 44 2.2.1 Calculation of induced connection on surfaces in E3. 45 2.3 Levi-Civita connection . 47 2.3.1 Symmetric connection . 47 2.3.2 Levi-Civita connection. Theorem and Explicit formulae 47 2.3.3 Levi-Civita connection on 2-dimensional Riemannian manifold with metric G = adu2 + bdv2. 49 2.3.4 Example of the sphere again . 49 2.4 Levi-Civita connection = induced connection on surfaces in E3 50 3 Parallel transport and geodesics 52 3.1 Parallel transport . 52 3.1.1 Definition . 52 3.1.2 ∗Parallel transport is a linear map . 53 3.1.3 Parallel transport with respect to Levi-Civita connection 54 3.2 Geodesics . 54 3.2.1 Definition. Geodesic on Riemannian manifold. 54 3.2.2 Un-parameterised geodesic . 55 3.2.3 Geodesics on surfaces in E3 . 56 3.3 Geodesics and Lagrangians of "free" particle on Riemannian manifold. 58 3.3.1 Lagrangian and Euler-Lagrange equations . 58 3.3.2 Lagrangian of "free" particle . 58 3.3.3 Equations of geodesics and Euler-Lagrange equations . 59 3.3.4 Examples of calculations of Christoffel symbols and geodesics using Lagrangians. 60 3.3.5 Variational principe and Euler-Lagrange equations . 62 3.4 Geodesics and shortest distance. 63 2 3.4.1 Again geodesics for sphere and Lobachevsky plane . 65 4 Surfaces in E3 67 4.1 Formulation of the main result. Theorem of parallel transport over closed curve and Theorema Egregium . 67 4.1.1 GaußTheorema Egregium . 69 4.2 Derivation formulae . 70 4.2.1 ∗Gauss condition (structure equations) . 72 4.3 Geometrical meaning of derivation formulae. Weingarten op- erator and second quadratic form in terms of derivation for- mulae. 73 4.3.1 Gaussian and mean curvature in terms of derivation formulae . 75 4.4 Examples of calculations of derivation formulae and curvatures for cylinder, cone and sphere . 76 4.5 ∗Proof of the Theorem of parallel transport along closed curve. 80 5 Curvtature tensor 84 5.1 Definition . 84 5.1.1 Properties of curvature tensor . 85 5.2 Riemann curvature tensor of Riemannian manifolds. 86 5.3 yCurvature of surfaces in E3.. Theorema Egregium again . 87 5.4 Relation between Gaussian curvature and Riemann curvature tensor. Straightforward proof of Theorema Egregium . 88 5.4.1 ∗Proof of the Proposition (5.25) . 90 5.5 Gauss Bonnet Theorem . 94 6 Appendices 97 6.1 ∗Integrals of motions and geodesics. 97 6.1.1 ∗Integral of motion for arbitrary Lagrangian L(x; x_) . 97 6.1.2 ∗Basic examples of Integrals of motion: Generalised momentum and Energy . 97 6.1.3 ∗Integrals of motion for geodesics . 98 6.1.4 ∗Using integral of motions to calculate geodesics . 100 6.2 Induced metric on surfaces. 101 6.2.1 Recalling Weingarten operator . 101 6.2.2 Second quadratic form . 102 6.2.3 Gaussian and mean curvatures . 103 3 6.2.4 Examples of calculation of Weingarten operator, Sec- ond quadratic forms, curvatures for cylinder, cone and sphere. 103 1 Riemannian manifolds 1.1 Manifolds. Tensors. (Recalling) I recall briefly basics of manifolds and tensor fields on manifolds. An n-dimensional manifold is a space such that in a vicinity of any point one can consider local coordinates fx1; : : : ; xng (charts). One can consider different local coordinates. If both coordinates fx1; : : : ; xng, fy1; : : : ; yng are defined in a vicinity of the given point then they are related by bijective transition functions (functions defined on domains in Rn and taking values in Rn). 8 > 10 10 1 n >x = x (x ; : : : ; x ) > 0 0 <>x2 = x2 (x1; : : : ; xn) >::: > − 0 − 0 >xn 1 = xn 1 (x1; : : : ; xn) :> 0 0 xn = xn (x1; : : : ; xn) We say that manifold is differentiable or smooth if transition functions are diffeomorphisms, i.e. they are smooth and rank of Jacobian is equal to k, i.e. 0 0 0 0 1 @x1 @x1 @x1 1 2 ::: n B @x 0 @x 0 @x 0 C B @x2 @x2 ::: @x2 C det B @x1 @x2 @xn C =6 0 (1.1) @ ::: A 0 0 0 @xn @xn @xn @x1 @x2 ::: @xn A good example of manifold is an open domain D in n-dimensional vector space Rn. Cartesian coordinates on Rn define global coordinates on D. On the other hand one can consider an arbitrary local coordinates in different domains in Rn. E.g. one can consider polar coordinates fr; 'g in a domain D = fx; y : y > 0g of R2 (or in other domain of R2) defined by standard formulae: ( x = r cos ' ; y = r sin ' 4 ! ( ) @x @x − @r @' cos ' r sin ' det @y @y = det = r (1.2) @r @' sin ' r cos ' or one can consider spherical coordinates fr; θ; 'g in a domain D = fx; y; z : x > 0; y > 0; z > 0g of R3 (or in other domain of R3) defined by standard for- mulae 8 <>x = r sin θ cos ' y = r sin θ sin ' ; :> z = r cos θ , 0 1 0 1 @x @x @x − B @r @θ @' C sin θ cos ' r cos θ cos ' r sin θ sin ' @ @y @y @y A @ A 2 det @r @θ @' = det sin θ sin ' r cos θ sin ' r sin θ cos ' = r sin θ @z @z @z cos θ −r sin θ 0 @r @θ @' (1.3) Choosing domain where polar (spherical) coordinates are well-defined we have to be award that coordinates have to be well-defined and transition functions (1.1) have to be diffeomorphisms. Examples of manifolds: Rn, Circle S1, Sphere S2, in general sphere Sn, torus S1 × S1, cylinder, cone, . We also have to recall briefly what are tensors on manifold. Tensors on Manifold For every point p on manifold M one can consider tangent vector space TpM|the space of vectors tangent to the manifold at the point M. i @ Tangent vector A = A @xi . Under changing of coordinates it transforms as follows: m0 @ @x @ 0 @ A = Ai = Ai = Am @xi @xi @xm0 @xm0 Hence i0 0 @x Ai = Ai (1.4) @xi ∗ Consider also cotangent space TpM (for every point p on manifold M)| space of linear functions on tangent vectors, i.e. space of 1-forms which sometimes are called covectors.: 5 i One-form (covector) ! = !idx transforms as follows m m @x m0 m0 ! = ! dx = ! dx = ! 0 dx : m m @xm0 m Hence @xm ! 0 = ! : (1.5) m @xm0 m Tensors: One can consider contravariant tensors of the rank p @ @ @ T = T i1i2:::ip ⊗ ⊗ · · · ⊗ @xi1 @xi2 @xik with components fT i1i2:::ik g. One can consider covariant tensors of the rank q j1 ⊗ j2 ⊗ jq S = Sj1j2:::jq dx dx : : : dx f g with components Sj1j2:::jq . One can also consider mixed tensors: i i :::i @ @ @ 1 2 p ⊗ ⊗ · · · ⊗ ⊗ j1 ⊗ j2 ⊗ jq Q = Qj j :::j dx dx : : : dx 1 2 q @xi1 @xi2 @xik ( ) p f i1i2:::ip g with components Qj j :::j . We call these tensors tensors of the type . (1 )2 q q p Tensors of the type are called contravariant tensors of the rank p. 0 ( ) 0 Tensors of the type are called covariant tensors of the rank q. q Having in mind( (1.4)) and (1.5) we come to the rule of transformation for p tensors of the type : q 0 0 0 0 0 0 i1 i2 ip j1 j2 jq i1i2:::ip @x @x @x @x @x @x i1i2:::ip Q 0 0 0 = ::: 0 0 ::: 0 Q (1.6) j j :::j i i i j j j j1j2:::jq 1 2 q @x 1 @x 2 @x p @x 1 @x 2 @x q ( ) p E.g. if S is a covariant tensor of rank 2 (tensor of the type ) then ik q @xi @xk S 0 0 = S : (1.7) i k @xi0 @xk0 ik 6 ( ) 1 If Ai is a tensor of rank (linear operator on T M) then k 1 p i0 k i0 @x @x i A 0 = A k @xi @xk0 k Remark Transformations formulae (1.4)|(1.7) define vectors, covectors and in generally any tensor fields in components.
Recommended publications
  • General Relativity Fall 2018 Lecture 10: Integration, Einstein-Hilbert Action
    General Relativity Fall 2018 Lecture 10: Integration, Einstein-Hilbert action Yacine Ali-Ha¨ımoud (Dated: October 3, 2018) σ σ HW comment: T µν antisym in µν does NOT imply that Tµν is antisym in µν. Volume element { Consider a LICS with primed coordinates. The 4-volume element is dV = d4x0 = 0 0 0 0 dx0 dx1 dx2 dx3 . If we change coordinates, we have 0 ! @xµ d4x0 = det d4x: (1) µ @x Now, the metric components change as 0 0 0 0 @xµ @xν @xµ @xν g = g 0 0 = η 0 0 ; (2) µν @xµ @xν µ ν @xµ @xν µ ν since the metric components are ηµ0ν0 in the LICS. Seing this as a matrix operation and taking the determinant, we find 0 !2 @xµ det(g ) = − det : (3) µν @xµ Hence, we find that the 4-volume element is q 4 p 4 p 4 dV = − det(gµν )d x ≡ −g d x ≡ jgj d x: (4) The integral of a scalar function f is well defined: given any coordinate system (even if only defined locally): Z Z f dV = f pjgj d4x: (5) We can only define the integral of a scalar function. The integral of a vector or tensor field is mean- ingless in curved spacetime. Think of the integral as a sum. To sum vectors, you need them to belong to the same vector space. There is no common vector space in curved spacetime. Only in flat spacetime can we define such integrals. First parallel-transport the vector field to a single point of spacetime (it doesn't matter which one).
    [Show full text]
  • Vector Calculus and Multiple Integrals Rob Fender, HT 2018
    Vector Calculus and Multiple Integrals Rob Fender, HT 2018 COURSE SYNOPSIS, RECOMMENDED BOOKS Course syllabus (on which exams are based): Double integrals and their evaluation by repeated integration in Cartesian, plane polar and other specified coordinate systems. Jacobians. Line, surface and volume integrals, evaluation by change of variables (Cartesian, plane polar, spherical polar coordinates and cylindrical coordinates only unless the transformation to be used is specified). Integrals around closed curves and exact differentials. Scalar and vector fields. The operations of grad, div and curl and understanding and use of identities involving these. The statements of the theorems of Gauss and Stokes with simple applications. Conservative fields. Recommended Books: Mathematical Methods for Physics and Engineering (Riley, Hobson and Bence) This book is lazily referred to as “Riley” throughout these notes (sorry, Drs H and B) You will all have this book, and it covers all of the maths of this course. However it is rather terse at times and you will benefit from looking at one or both of these: Introduction to Electrodynamics (Griffiths) You will buy this next year if you haven’t already, and the chapter on vector calculus is very clear Div grad curl and all that (Schey) A nice discussion of the subject, although topics are ordered differently to most courses NB: the latest version of this book uses the opposite convention to polar coordinates to this course (and indeed most of physics), but older versions can often be found in libraries 1 Week One A review of vectors, rotation of coordinate systems, vector vs scalar fields, integrals in more than one variable, first steps in vector differentiation, the Frenet-Serret coordinate system Lecture 1 Vectors A vector has direction and magnitude and is written in these notes in bold e.g.
    [Show full text]
  • Mathematical Theorems
    Appendix A Mathematical Theorems The mathematical theorems needed in order to derive the governing model equations are defined in this appendix. A.1 Transport Theorem for a Single Phase Region The transport theorem is employed deriving the conservation equations in continuum mechanics. The mathematical statement is sometimes attributed to, or named in honor of, the German Mathematician Gottfried Wilhelm Leibnitz (1646–1716) and the British fluid dynamics engineer Osborne Reynolds (1842–1912) due to their work and con- tributions related to the theorem. Hence it follows that the transport theorem, or alternate forms of the theorem, may be named the Leibnitz theorem in mathematics and Reynolds transport theorem in mechanics. In a customary interpretation the Reynolds transport theorem provides the link between the system and control volume representations, while the Leibnitz’s theorem is a three dimensional version of the integral rule for differentiation of an integral. There are several notations used for the transport theorem and there are numerous forms and corollaries. A.1.1 Leibnitz’s Rule The Leibnitz’s integral rule gives a formula for differentiation of an integral whose limits are functions of the differential variable [7, 8, 22, 23, 45, 55, 79, 94, 99]. The formula is also known as differentiation under the integral sign. H. A. Jakobsen, Chemical Reactor Modeling, DOI: 10.1007/978-3-319-05092-8, 1361 © Springer International Publishing Switzerland 2014 1362 Appendix A: Mathematical Theorems b(t) b(t) d ∂f (t, x) db da f (t, x) dx = dx + f (t, b) − f (t, a) (A.1) dt ∂t dt dt a(t) a(t) The first term on the RHS gives the change in the integral because the function itself is changing with time, the second term accounts for the gain in area as the upper limit is moved in the positive axis direction, and the third term accounts for the loss in area as the lower limit is moved.
    [Show full text]
  • Chapter 13 Curvature in Riemannian Manifolds
    Chapter 13 Curvature in Riemannian Manifolds 13.1 The Curvature Tensor If (M, , )isaRiemannianmanifoldand is a connection on M (that is, a connection on TM−), we− saw in Section 11.2 (Proposition 11.8)∇ that the curvature induced by is given by ∇ R(X, Y )= , ∇X ◦∇Y −∇Y ◦∇X −∇[X,Y ] for all X, Y X(M), with R(X, Y ) Γ( om(TM,TM)) = Hom (Γ(TM), Γ(TM)). ∈ ∈ H ∼ C∞(M) Since sections of the tangent bundle are vector fields (Γ(TM)=X(M)), R defines a map R: X(M) X(M) X(M) X(M), × × −→ and, as we observed just after stating Proposition 11.8, R(X, Y )Z is C∞(M)-linear in X, Y, Z and skew-symmetric in X and Y .ItfollowsthatR defines a (1, 3)-tensor, also denoted R, with R : T M T M T M T M. p p × p × p −→ p Experience shows that it is useful to consider the (0, 4)-tensor, also denoted R,givenby R (x, y, z, w)= R (x, y)z,w p p p as well as the expression R(x, y, y, x), which, for an orthonormal pair, of vectors (x, y), is known as the sectional curvature, K(x, y). This last expression brings up a dilemma regarding the choice for the sign of R. With our present choice, the sectional curvature, K(x, y), is given by K(x, y)=R(x, y, y, x)but many authors define K as K(x, y)=R(x, y, x, y). Since R(x, y)isskew-symmetricinx, y, the latter choice corresponds to using R(x, y)insteadofR(x, y), that is, to define R(X, Y ) by − R(X, Y )= + .
    [Show full text]
  • Applications of the Hodge Decomposition to Biological Structure and Function Modeling
    Applications of the Hodge Decomposition to Biological Structure and Function Modeling Andrew Gillette joint work with Chandrajit Bajaj and John Luecke Department of Mathematics, Institute of Computational Engineering and Sciences University of Texas at Austin, Austin, Texas 78712, USA http://www.math.utexas.edu/users/agillette university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 1/32 University Introduction Molecular dynamics are governed by electrostatic forces of attraction and repulsion. These forces are described as the solutions of a PDE over the molecular surfaces. Molecular surfaces may have complicated topological features affecting the solution. The Hodge Decomposition relates topological properties of the surface to solution spaces of PDEs over the surface. ∼ (space of forms) = (solutions to ∆u = f 6≡ 0) ⊕ (non-trivial deRham classes)university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 2/32 University Outline 1 The Hodge Decomposition for smooth differential forms 2 The Hodge Decomposition for discrete differential forms 3 Applications of the Hodge Decomposition to biological modeling university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 3/32 University Outline 1 The Hodge Decomposition for smooth differential forms 2 The Hodge Decomposition for discrete differential forms 3 Applications of the Hodge Decomposition to biological modeling university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 4/32 University Differential Forms Let Ω denote a smooth n-manifold and Tx (Ω) the tangent space of Ω at x.
    [Show full text]
  • Hodge Theory
    HODGE THEORY PETER S. PARK Abstract. This exposition of Hodge theory is a slightly retooled version of the author's Harvard minor thesis, advised by Professor Joe Harris. Contents 1. Introduction 1 2. Hodge Theory of Compact Oriented Riemannian Manifolds 2 2.1. Hodge star operator 2 2.2. The main theorem 3 2.3. Sobolev spaces 5 2.4. Elliptic theory 11 2.5. Proof of the main theorem 14 3. Hodge Theory of Compact K¨ahlerManifolds 17 3.1. Differential operators on complex manifolds 17 3.2. Differential operators on K¨ahlermanifolds 20 3.3. Bott{Chern cohomology and the @@-Lemma 25 3.4. Lefschetz decomposition and the Hodge index theorem 26 Acknowledgments 30 References 30 1. Introduction Our objective in this exposition is to state and prove the main theorems of Hodge theory. In Section 2, we first describe a key motivation behind the Hodge theory for compact, closed, oriented Riemannian manifolds: the observation that the differential forms that satisfy certain par- tial differential equations depending on the choice of Riemannian metric (forms in the kernel of the associated Laplacian operator, or harmonic forms) turn out to be precisely the norm-minimizing representatives of the de Rham cohomology classes. This naturally leads to the statement our first main theorem, the Hodge decomposition|for a given compact, closed, oriented Riemannian manifold|of the space of smooth k-forms into the image of the Laplacian and its kernel, the sub- space of harmonic forms. We then develop the analytic machinery|specifically, Sobolev spaces and the theory of elliptic differential operators|that we use to prove the aforementioned decom- position, which immediately yields as a corollary the phenomenon of Poincar´eduality.
    [Show full text]
  • Curvature of Riemannian Manifolds
    Curvature of Riemannian Manifolds Seminar Riemannian Geometry Summer Term 2015 Prof. Dr. Anna Wienhard and Dr. Gye-Seon Lee Soeren Nolting July 16, 2015 1 Motivation Figure 1: A vector parallel transported along a closed curve on a curved manifold.[1] The aim of this talk is to define the curvature of Riemannian Manifolds and meeting some important simplifications as the sectional, Ricci and scalar curvature. We have already noticed, that a vector transported parallel along a closed curve on a Riemannian Manifold M may change its orientation. Thus, we can determine whether a Riemannian Manifold is curved or not by transporting a vector around a loop and measuring the difference of the orientation at start and the endo of the transport. As an example take Figure 1, which depicts a parallel transport of a vector on a two-sphere. Note that in a non-curved space the orientation of the vector would be preserved along the transport. 1 2 Curvature In the following we will use the Einstein sum convention and make use of the notation: X(M) space of smooth vector fields on M D(M) space of smooth functions on M 2.1 Defining Curvature and finding important properties This rather geometrical approach motivates the following definition: Definition 2.1 (Curvature). The curvature of a Riemannian Manifold is a correspondence that to each pair of vector fields X; Y 2 X (M) associates the map R(X; Y ): X(M) ! X(M) defined by R(X; Y )Z = rX rY Z − rY rX Z + r[X;Y ]Z (1) r is the Riemannian connection of M.
    [Show full text]
  • Vector Calculus Operators
    Vector Calculus Operators July 6, 2015 1 Div, Grad and Curl The divergence of a vector quantity in Cartesian coordinates is defined by: −! @ @ @ r · ≡ x + y + z (1) @x @y @z −! where = ( x; y; z). The divergence of a vector quantity is a scalar quantity. If the divergence is greater than zero, it implies a net flux out of a volume element. If the divergence is less than zero, it implies a net flux into a volume element. If the divergence is exactly equal to zero, then there is no net flux into or out of the volume element, i.e., any field lines that enter the volume element also exit the volume element. The divergence of the electric field and the divergence of the magnetic field are two of Maxwell's equations. The gradient of a scalar quantity in Cartesian coordinates is defined by: @ @ @ r (x; y; z) ≡ ^{ + |^ + k^ (2) @x @y @z where ^{ is the unit vector in the x-direction,| ^ is the unit vector in the y-direction, and k^ is the unit vector in the z-direction. Evaluating the gradient indicates the direction in which the function is changing the most rapidly (i.e., slope). If the gradient of a function is zero, the function is not changing. The curl of a vector function in Cartesian coordinates is given by: 1 −! @ @ @ @ @ @ r × ≡ ( z − y )^{ + ( x − z )^| + ( y − x )k^ (3) @y @z @z @x @x @y where the resulting function is another vector function that indicates the circulation of the original function, with unit vectors pointing orthogonal to the plane of the circulation components.
    [Show full text]
  • A New Definition of the Representative Volument Element in Numerical
    UNIVERSITE´ DE MONTREAL´ A NEW DEFINITION OF THE REPRESENTATIVE VOLUMENT ELEMENT IN NUMERICAL HOMOGENIZATION PROBLEMS AND ITS APPLICATION TO THE PERFORMANCE EVALUATION OF ANALYTICAL HOMOGENIZATION MODELS HADI MOUSSADDY DEPARTEMENT´ DE GENIE´ MECANIQUE´ ECOLE´ POLYTECHNIQUE DE MONTREAL´ THESE` PRESENT´ EE´ EN VUE DE L'OBTENTION DU DIPLOME^ DE PHILOSOPHIÆ DOCTOR (GENIE´ MECANIQUE)´ AVRIL 2013 © Hadi Moussaddy, 2013. UNIVERSITE´ DE MONTREAL´ ECOLE´ POLYTECHNIQUE DE MONTREAL´ Cette th`eseintitul´ee : A NEW DEFINITION OF THE REPRESENTATIVE VOLUMENT ELEMENT IN NUMERICAL HOMOGENIZATION PROBLEMS AND ITS APPLICATION TO THE PERFORMANCE EVALUATION OF ANALYTICAL HOMOGENIZATION MODELS pr´esent´ee par : MOUSSADDY Hadi en vue de l'obtention du dipl^ome de : Philosophiæ Doctor a ´et´ed^ument accept´eepar le jury d'examen constitu´ede : M. GOSSELIN Fr´ed´eric, Ph.D., pr´esident M. THERRIAULT Daniel, Ph.D., membre et directeur de recherche M. LE´VESQUE Martin, Ph.D., membre et codirecteur de recherche M. LESSARD Larry, Ph.D., membre M. BOHM¨ Helmut J., Techn. Doct., membre iii To my parents Safwan and Aisha To my wife Nour To my whole family... iv ACKNOWLEDGEMENTS I would like to thank Pr. Daniel Therriault and Pr. Martin L´evesque for their guidance throughout this project. Your insightful criticism guided me to complete a rigorous work, while gaining scientific maturity. Also, I want to thank Pr. B¨ohm, Pr. Lessard and Pr. Henry for being my committee members and Pr. Gosselin for being the Jury president. I owe special gratitude to the members of Laboratory for Multiscale Mechanics (LM2) research group. I would like to paricularly acknowledge the contributions of Maryam Pah- lavanPour by conducting parts of the analytical modeling included in this thesis, and for numerous fruitful dicussions all along this project.
    [Show full text]
  • Electric Flux Density, Gauss's Law, and Divergence
    www.getmyuni.com CHAPTER3 ELECTRIC FLUX DENSITY, GAUSS'S LAW, AND DIVERGENCE After drawinga few of the fields described in the previous chapter and becoming familiar with the concept of the streamlines which show the direction of the force on a test charge at every point, it is difficult to avoid giving these lines a physical significance and thinking of them as flux lines. No physical particle is projected radially outward from the point charge, and there are no steel tentacles reaching out to attract or repel an unwary test charge, but as soon as the streamlines are drawn on paper there seems to be a picture showing``something''is present. It is very helpful to invent an electric flux which streams away symmetri- cally from a point charge and is coincident with the streamlines and to visualize this flux wherever an electric field is present. This chapter introduces and uses the concept of electric flux and electric flux density to solve again several of the problems presented in the last chapter. The work here turns out to be much easier, and this is due to the extremely symmetrical problems which we are solving. www.getmyuni.com 3.1 ELECTRIC FLUX DENSITY About 1837 the Director of the Royal Society in London, Michael Faraday, became very interested in static electric fields and the effect of various insulating materials on these fields. This problem had been botheringhim duringthe past ten years when he was experimentingin his now famous work on induced elec- tromotive force, which we shall discuss in Chap. 10.
    [Show full text]
  • Semi-Riemannian Manifold Optimization
    Semi-Riemannian Manifold Optimization Tingran Gao William H. Kruskal Instructor Committee on Computational and Applied Mathematics (CCAM) Department of Statistics The University of Chicago 2018 China-Korea International Conference on Matrix Theory with Applications Shanghai University, Shanghai, China December 19, 2018 Outline Motivation I Riemannian Manifold Optimization I Barrier Functions and Interior Point Methods Semi-Riemannian Manifold Optimization I Optimality Conditions I Descent Direction I Semi-Riemannian First Order Methods I Metric Independence of Second Order Methods Optimization on Semi-Riemannian Submanifolds Joint work with Lek-Heng Lim (The University of Chicago) and Ke Ye (Chinese Academy of Sciences) Riemannian Manifold Optimization I Optimization on manifolds: min f (x) x2M where (M; g) is a (typically nonlinear, non-convex) Riemannian manifold, f : M! R is a smooth function on M I Difficult constrained optimization; handled as unconstrained optimization from the perspective of manifold optimization I Key ingredients: tangent spaces, geodesics, exponential map, parallel transport, gradient rf , Hessian r2f Riemannian Manifold Optimization I Riemannian Steepest Descent x = Exp (−trf (x )) ; t ≥ 0 k+1 xk k I Riemannian Newton's Method 2 r f (xk ) ηk = −∇f (xk ) x = Exp (−tη ) ; t ≥ 0 k+1 xk k I Riemannian trust region I Riemannian conjugate gradient I ...... Gabay (1982), Smith (1994), Edelman et al. (1998), Absil et al. (2008), Adler et al. (2002), Ring and Wirth (2012), Huang et al. (2015), Sato (2016), etc. Riemannian Manifold Optimization I Applications: Optimization problems on matrix manifolds I Stiefel manifolds n Vk (R ) = O (n) =O (k) I Grassmannian manifolds Gr (k; n) = O (n) = (O (k) × O (n − k)) I Flag manifolds: for n1 + ··· + nd = n, Flag = O (n) = (O (n ) × · · · × O (n )) n1;··· ;nd 1 d I Shape Spaces k×n R n f0g =O (n) I .....
    [Show full text]
  • Deep Riemannian Manifold Learning
    Deep Riemannian Manifold Learning Aaron Lou1∗ Maximilian Nickel2 Brandon Amos2 1Cornell University 2 Facebook AI Abstract We present a new class of learnable Riemannian manifolds with a metric param- eterized by a deep neural network. The core manifold operations–specifically the Riemannian exponential and logarithmic maps–are solved using approximate numerical techniques. Input and parameter gradients are computed with an ad- joint sensitivity analysis. This enables us to fit geodesics and distances with gradient-based optimization of both on-manifold values and the manifold itself. We demonstrate our method’s capability to model smooth, flexible metric structures in graph embedding tasks. 1 Introduction Geometric domain knowledge is important for machine learning systems to interact, process, and produce data from inherently geometric spaces and phenomena [BBL+17, MBM+17]. Geometric knowledge often empowers the model with explicit operations and components that reflect the geometric properties of the domain. One commonly used construct is a Riemannian manifold, the generalization of standard Euclidean space. These structures enable models to represent non-trivial geometric properties of data [NK17, DFDC+18, GDM+18, NK18, GBH18]. However, despite this capability, the current set of usable Riemannian manifolds is surprisingly small, as the core manifold operations must be provided in a closed form to enable gradient-based optimization. Worse still, these manifolds are often set a priori, as the techniques for interpolating between manifolds are extremely restricted [GSGR19, SGB20]. In this work, we study how to integrate a general class of Riemannian manifolds into learning systems. Specifically, we parameterize a class of Riemannian manifolds with a deep neural network and develop techniques to optimize through the manifold operations.
    [Show full text]