04 - Tensor Calculus - Tensor Analysis

Total Page:16

File Type:pdf, Size:1020Kb

04 - Tensor Calculus - Tensor Analysis 04 - tensor calculus - tensor analysis 04 - tensor calculus 1 tensor analysis - frechet derivative • consider smooth differentiable scalar field with scalar argument vector argument tensor argument • frechet derivative (tensor notation) scalar argument vector argument tensor argument tensor calculus 2 tensor analysis - gateaux derivative • consider smooth differentiable scalar field with scalar argument vector argument tensor argument • gateaux derivative,i.e.,frechet wrt direction (tensor notation) scalar argument vector argument tensor argument tensor calculus 3 tensor analysis - gradient • consider scalar- and vector field in domain • gradient of scalar- and vector field renders vector- and 2nd order tensor field tensor calculus 4 tensor analysis - divergence • consider vector- and 2nd order tensor field in domain • divergence of vector- and 2nd order tensor field renders scalar- and vector field tensor calculus 5 tensor analysis - laplace operator • consider scalar- and vector field in domain • laplace operator acting on scalar- and vector field renders scalar- and vector field tensor calculus 6 tensor analysis - transformation formulae • consider scalar,vector and 2nd order tensor field on • useful transformation formulae (tensor notation) tensor calculus 7 tensor analysis - transformation formulae • consider scalar,vector and 2nd order tensor field on • useful transformation formulae (index notation) tensor calculus 8 tensor analysis - integral theorems • consider scalar,vector and 2nd order tensor field on • integral theorems (tensor notation) green gauss gauss tensor calculus 9 tensor analysis - integral theorems • consider scalar,vector and 2nd order tensor field on • integral theorems (tensor notation) green gauss gauss tensor calculus 10 voigt / matrix vector notation • strain tensors as vectors in voigt notation • stress tensors as vectors in voigt notation • why are strain & stress different? check energy expression! tensor calculus 11 voigt / matrix vector notation • fourth order material operators as matrix in voigt notation • why are strain & stress different? check these expressions! tensor calculus 12.
Recommended publications
  • Arxiv:1507.07356V2 [Math.AP]
    TEN EQUIVALENT DEFINITIONS OF THE FRACTIONAL LAPLACE OPERATOR MATEUSZ KWAŚNICKI Abstract. This article discusses several definitions of the fractional Laplace operator ( ∆)α/2 (α (0, 2)) in Rd (d 1), also known as the Riesz fractional derivative − ∈ ≥ operator, as an operator on Lebesgue spaces L p (p [1, )), on the space C0 of ∈ ∞ continuous functions vanishing at infinity and on the space Cbu of bounded uniformly continuous functions. Among these definitions are ones involving singular integrals, semigroups of operators, Bochner’s subordination and harmonic extensions. We collect and extend known results in order to prove that all these definitions agree: on each of the function spaces considered, the corresponding operators have common domain and they coincide on that common domain. 1. Introduction We consider the fractional Laplace operator L = ( ∆)α/2 in Rd, with α (0, 2) and d 1, 2, ... Numerous definitions of L can be found− in− literature: as a Fourier∈ multiplier with∈{ symbol} ξ α, as a fractional power in the sense of Bochner or Balakrishnan, as the inverse of the−| Riesz| potential operator, as a singular integral operator, as an operator associated to an appropriate Dirichlet form, as an infinitesimal generator of an appropriate semigroup of contractions, or as the Dirichlet-to-Neumann operator for an appropriate harmonic extension problem. Equivalence of these definitions for sufficiently smooth functions is well-known and easy. The purpose of this article is to prove that, whenever meaningful, all these definitions are equivalent in the Lebesgue space L p for p [1, ), ∈ ∞ in the space C0 of continuous functions vanishing at infinity, and in the space Cbu of bounded uniformly continuous functions.
    [Show full text]
  • 1.2 Topological Tensor Calculus
    PH211 Physical Mathematics Fall 2019 1.2 Topological tensor calculus 1.2.1 Tensor fields Finite displacements in Euclidean space can be represented by arrows and have a natural vector space structure, but finite displacements in more general curved spaces, such as on the surface of a sphere, do not. However, an infinitesimal neighborhood of a point in a smooth curved space1 looks like an infinitesimal neighborhood of Euclidean space, and infinitesimal displacements dx~ retain the vector space structure of displacements in Euclidean space. An infinitesimal neighborhood of a point can be infinitely rescaled to generate a finite vector space, called the tangent space, at the point. A vector lives in the tangent space of a point. Note that vectors do not stretch from one point to vector tangent space at p p space Figure 1.2.1: A vector in the tangent space of a point. another, and vectors at different points live in different tangent spaces and so cannot be added. For example, rescaling the infinitesimal displacement dx~ by dividing it by the in- finitesimal scalar dt gives the velocity dx~ ~v = (1.2.1) dt which is a vector. Similarly, we can picture the covector rφ as the infinitesimal contours of φ in a neighborhood of a point, infinitely rescaled to generate a finite covector in the point's cotangent space. More generally, infinitely rescaling the neighborhood of a point generates the tensor space and its algebra at the point. The tensor space contains the tangent and cotangent spaces as a vector subspaces. A tensor field is something that takes tensor values at every point in a space.
    [Show full text]
  • Curl, Divergence and Laplacian
    Curl, Divergence and Laplacian What to know: 1. The definition of curl and it two properties, that is, theorem 1, and be able to predict qualitatively how the curl of a vector field behaves from a picture. 2. The definition of divergence and it two properties, that is, if div F~ 6= 0 then F~ can't be written as the curl of another field, and be able to tell a vector field of clearly nonzero,positive or negative divergence from the picture. 3. Know the definition of the Laplace operator 4. Know what kind of objects those operator take as input and what they give as output. The curl operator Let's look at two plots of vector fields: Figure 1: The vector field Figure 2: The vector field h−y; x; 0i: h1; 1; 0i We can observe that the second one looks like it is rotating around the z axis. We'd like to be able to predict this kind of behavior without having to look at a picture. We also promised to find a criterion that checks whether a vector field is conservative in R3. Both of those goals are accomplished using a tool called the curl operator, even though neither of those two properties is exactly obvious from the definition we'll give. Definition 1. Let F~ = hP; Q; Ri be a vector field in R3, where P , Q and R are continuously differentiable. We define the curl operator: @R @Q @P @R @Q @P curl F~ = − ~i + − ~j + − ~k: (1) @y @z @z @x @x @y Remarks: 1.
    [Show full text]
  • A Some Basic Rules of Tensor Calculus
    A Some Basic Rules of Tensor Calculus The tensor calculus is a powerful tool for the description of the fundamentals in con- tinuum mechanics and the derivation of the governing equations for applied prob- lems. In general, there are two possibilities for the representation of the tensors and the tensorial equations: – the direct (symbolic) notation and – the index (component) notation The direct notation operates with scalars, vectors and tensors as physical objects defined in the three dimensional space. A vector (first rank tensor) a is considered as a directed line segment rather than a triple of numbers (coordinates). A second rank tensor A is any finite sum of ordered vector pairs A = a b + ... +c d. The scalars, vectors and tensors are handled as invariant (independent⊗ from the choice⊗ of the coordinate system) objects. This is the reason for the use of the direct notation in the modern literature of mechanics and rheology, e.g. [29, 32, 49, 123, 131, 199, 246, 313, 334] among others. The index notation deals with components or coordinates of vectors and tensors. For a selected basis, e.g. gi, i = 1, 2, 3 one can write a = aig , A = aibj + ... + cidj g g i i ⊗ j Here the Einstein’s summation convention is used: in one expression the twice re- peated indices are summed up from 1 to 3, e.g. 3 3 k k ik ik a gk ∑ a gk, A bk ∑ A bk ≡ k=1 ≡ k=1 In the above examples k is a so-called dummy index. Within the index notation the basic operations with tensors are defined with respect to their coordinates, e.
    [Show full text]
  • Tensor Calculus and Differential Geometry
    Course Notes Tensor Calculus and Differential Geometry 2WAH0 Luc Florack March 10, 2021 Cover illustration: papyrus fragment from Euclid’s Elements of Geometry, Book II [8]. Contents Preface iii Notation 1 1 Prerequisites from Linear Algebra 3 2 Tensor Calculus 7 2.1 Vector Spaces and Bases . .7 2.2 Dual Vector Spaces and Dual Bases . .8 2.3 The Kronecker Tensor . 10 2.4 Inner Products . 11 2.5 Reciprocal Bases . 14 2.6 Bases, Dual Bases, Reciprocal Bases: Mutual Relations . 16 2.7 Examples of Vectors and Covectors . 17 2.8 Tensors . 18 2.8.1 Tensors in all Generality . 18 2.8.2 Tensors Subject to Symmetries . 22 2.8.3 Symmetry and Antisymmetry Preserving Product Operators . 24 2.8.4 Vector Spaces with an Oriented Volume . 31 2.8.5 Tensors on an Inner Product Space . 34 2.8.6 Tensor Transformations . 36 2.8.6.1 “Absolute Tensors” . 37 CONTENTS i 2.8.6.2 “Relative Tensors” . 38 2.8.6.3 “Pseudo Tensors” . 41 2.8.7 Contractions . 43 2.9 The Hodge Star Operator . 43 3 Differential Geometry 47 3.1 Euclidean Space: Cartesian and Curvilinear Coordinates . 47 3.2 Differentiable Manifolds . 48 3.3 Tangent Vectors . 49 3.4 Tangent and Cotangent Bundle . 50 3.5 Exterior Derivative . 51 3.6 Affine Connection . 52 3.7 Lie Derivative . 55 3.8 Torsion . 55 3.9 Levi-Civita Connection . 56 3.10 Geodesics . 57 3.11 Curvature . 58 3.12 Push-Forward and Pull-Back . 59 3.13 Examples . 60 3.13.1 Polar Coordinates in the Euclidean Plane .
    [Show full text]
  • Fast Higher-Order Functions for Tensor Calculus with Tensors and Subtensors
    Fast Higher-Order Functions for Tensor Calculus with Tensors and Subtensors Cem Bassoy and Volker Schatz Fraunhofer IOSB, Ettlingen 76275, Germany, [email protected] Abstract. Tensors analysis has become a popular tool for solving prob- lems in computational neuroscience, pattern recognition and signal pro- cessing. Similar to the two-dimensional case, algorithms for multidimen- sional data consist of basic operations accessing only a subset of tensor data. With multiple offsets and step sizes, basic operations for subtensors require sophisticated implementations even for entrywise operations. In this work, we discuss the design and implementation of optimized higher-order functions that operate entrywise on tensors and subtensors with any non-hierarchical storage format and arbitrary number of dimen- sions. We propose recursive multi-index algorithms with reduced index computations and additional optimization techniques such as function inlining with partial template specialization. We show that single-index implementations of higher-order functions with subtensors introduce a runtime penalty of an order of magnitude than the recursive and iterative multi-index versions. Including data- and thread-level parallelization, our optimized implementations reach 68% of the maximum throughput of an Intel Core i9-7900X. In comparison with other libraries, the aver- age speedup of our optimized implementations is up to 5x for map-like and more than 9x for reduce-like operations. For symmetric tensors we measured an average speedup of up to 4x. 1 Introduction Many problems in computational neuroscience, pattern recognition, signal pro- cessing and data mining generate massive amounts of multidimensional data with high dimensionality [13,12,9]. Tensors provide a natural representation for massive multidimensional data [10,7].
    [Show full text]
  • Part IA — Vector Calculus
    Part IA | Vector Calculus Based on lectures by B. Allanach Notes taken by Dexter Chua Lent 2015 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after lectures. They are nowhere near accurate representations of what was actually lectured, and in particular, all errors are almost surely mine. 3 Curves in R 3 Parameterised curves and arc length, tangents and normals to curves in R , the radius of curvature. [1] 2 3 Integration in R and R Line integrals. Surface and volume integrals: definitions, examples using Cartesian, cylindrical and spherical coordinates; change of variables. [4] Vector operators Directional derivatives. The gradient of a real-valued function: definition; interpretation as normal to level surfaces; examples including the use of cylindrical, spherical *and general orthogonal curvilinear* coordinates. Divergence, curl and r2 in Cartesian coordinates, examples; formulae for these oper- ators (statement only) in cylindrical, spherical *and general orthogonal curvilinear* coordinates. Solenoidal fields, irrotational fields and conservative fields; scalar potentials. Vector derivative identities. [5] Integration theorems Divergence theorem, Green's theorem, Stokes's theorem, Green's second theorem: statements; informal proofs; examples; application to fluid dynamics, and to electro- magnetism including statement of Maxwell's equations. [5] Laplace's equation 2 3 Laplace's equation in R and R : uniqueness theorem and maximum principle. Solution of Poisson's equation by Gauss's method (for spherical and cylindrical symmetry) and as an integral. [4] 3 Cartesian tensors in R Tensor transformation laws, addition, multiplication, contraction, with emphasis on tensors of second rank. Isotropic second and third rank tensors.
    [Show full text]
  • 4.1 Discrete Differential Geometry
    Spring 2015 CSCI 599: Digital Geometry Processing 4.1 Discrete Differential Geometry Hao Li http://cs599.hao-li.com 1 Outline • Discrete Differential Operators • Discrete Curvatures • Mesh Quality Measures 2 Differential Operators on Polygons Differential Properties! • Surface is sufficiently differentiable • Curvatures → 2nd derivatives 3 Differential Operators on Polygons Differential Properties! • Surface is sufficiently differentiable • Curvatures → 2nd derivatives Polygonal Meshes! • Piecewise linear approximations of smooth surface • Focus on Discrete Laplace Beltrami Operator • Discrete differential properties defined over (x) N 4 Local Averaging Local Neighborhood ( x ) of a point !x N • often coincides with mesh vertex vi • n-ring neighborhood ( v ) or local geodesic ball Nn i 5 Local Averaging Local Neighborhood ( x ) of a point !x N • often coincides with mesh vertex vi • n-ring neighborhood ( v ) or local geodesic ball Nn i Neighborhood size! • Large: smoothing is introduced, stable to noise • Small: fine scale variation, sensitive to noise 6 Local Averaging: 1-Ring (x) N x Barycentric cell Voronoi cell Mixed Voronoi cell (barycenters/edgemidpoints) (circumcenters) (circumcenters/midpoint) tight error bound better approximation 7 Barycentric Cells Connect edge midpoints and triangle barycenters! • Simple to compute • Area is 1/3 o triangle areas • Slightly wrong for obtuse triangles 8 Mixed Cells Connect edge midpoints and! • Circumcenters for non-obtuse triangles • Midpoint of opposite edge for obtuse triangles • Better approximation,
    [Show full text]
  • A Doubly Nonlocal Laplace Operator and Its Connection to the Classical Laplacian Petronela Radu University of Nebraska - Lincoln, [email protected]
    University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Publications, Department of Mathematics Mathematics, Department of 2019 A Doubly Nonlocal Laplace Operator and Its Connection to the Classical Laplacian Petronela Radu University of Nebraska - Lincoln, [email protected] Kelseys Wells University of Nebraska - Lincoln Follow this and additional works at: https://digitalcommons.unl.edu/mathfacpub Part of the Applied Mathematics Commons, and the Mathematics Commons Radu, Petronela and Wells, Kelseys, "A Doubly Nonlocal Laplace Operator and Its Connection to the Classical Laplacian" (2019). Faculty Publications, Department of Mathematics. 215. https://digitalcommons.unl.edu/mathfacpub/215 This Article is brought to you for free and open access by the Mathematics, Department of at DigitalCommons@University of Nebraska - Lincoln. It has been accepted for inclusion in Faculty Publications, Department of Mathematics by an authorized administrator of DigitalCommons@University of Nebraska - Lincoln. A DOUBLY NONLOCAL LAPLACE OPERATOR AND ITS CONNECTION TO THE CLASSICAL LAPLACIAN PETRONELA RADU, KELSEY WELLS Abstract. In this paper, motivated by the state-based peridynamic frame- work, we introduce a new nonlocal Laplacian that exhibits double nonlocality through the use of iterated integral operators. The operator introduces addi- tional degrees of flexibility that can allow for better representation of physical phenomena at different scales and in materials with different properties. We study mathematical properties of this state-based Laplacian, including connec- tions with other nonlocal and local counterparts. Finally, we obtain explicit rates of convergence for this doubly nonlocal operator to the classical Laplacian as the radii for the horizons of interaction kernels shrink to zero. 1. Introduction Physical phenomena that are beset by discontinuities in the solution, or in the domain, have been challenging to study in the context of classical partial differen- tial equations (PDEs).
    [Show full text]
  • The Electrostatic Potential Is the Coulomb Integral Transform of the Electric Charge Density Revista Mexicana De Física, Vol
    Revista Mexicana de Física ISSN: 0035-001X [email protected] Sociedad Mexicana de Física A.C. México Medina, L.; Ley Koo, E. Mathematics motivated by physics: the electrostatic potential is the Coulomb integral transform of the electric charge density Revista Mexicana de Física, vol. 54, núm. 2, diciembre, 2008, pp. 153-159 Sociedad Mexicana de Física A.C. Distrito Federal, México Available in: http://www.redalyc.org/articulo.oa?id=57028302013 How to cite Complete issue Scientific Information System More information about this article Network of Scientific Journals from Latin America, the Caribbean, Spain and Portugal Journal's homepage in redalyc.org Non-profit academic project, developed under the open access initiative ENSENANZA˜ REVISTA MEXICANA DE FISICA´ E 54 (2) 153–159 DICIEMBRE 2008 Mathematics motivated by physics: the electrostatic potential is the Coulomb integral transform of the electric charge density L. Medinaa and E. Ley Koo a;b aFacultad de Ciencias, Universidad Nacional Autonoma´ de Mexico,´ Ciudad Universitaria, Mexico´ D.F., 04510, Mexico,´ e-mail: [email protected], bInstituto de F´ısica, Universidad Nacional Autonoma´ de Mexico,´ Apartado Postal 20-364, 01000 Mexico D.F., Mexico, e-mail: eleykoo@fisica.unam.mx Recibido el 19 de octubre de 2007; aceptado el 11 de marzo de 2007 This article illustrates a practical way to connect and coordinate the teaching and learning of physics and mathematics. The starting point is the electrostatic potential, which is obtained in any introductory course of electromagnetism from the Coulomb potential and the superposition principle for any charge distribution. The necessity to develop solutions to the Laplace and Poisson differential equations is also recognized, identifying the Coulomb potential as the generating function of harmonic functions.
    [Show full text]
  • Nonlocal Exterior Calculus on Riemannian Manifolds
    The Pennsylvania State University The Graduate School Department of Mathematics NONLOCAL EXTERIOR CALCULUS ON RIEMANNIAN MANIFOLDS A Dissertation in Mathematics by Thinh Duc Le c 2013 Thinh Duc Le Submitted in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy August 2013 ii The dissertation of Thinh Duc Le was reviewed and approved* by the following: Qiang Du Verne M. Willaman Professor of Mathematics Dissertation Adviser Chair of Committee Long-Qing Chen Distinguished Professor of Material Sciences and Engineering Ping Xu Distinguished Professor of Mathematics Mathieu Stienon Associate Professor of Mathematics Svetlana Katok Director of Graduate Studies, Department of Mathematics *Signatures are on file in the Graduate School. iii Abstract Exterior calculus and differential forms are basic mathematical concepts that have been around for centuries. Variations of these concepts have also been made over the years such as the discrete exterior calculus and the finite element exterior calculus. In this work, motivated by the recent studies of nonlocal vector calculus we develop a nonlocal exterior calculus framework on Riemannian manifolds which mimics many properties of the standard (local/smooth) exterior calculus. However the key difference is that nonlocal \interactions" (functions, operators, fields,...) are not required to be smooth. Also any point/particle can interact directly with any other point/particle in the studied domain (at least in principle). Just as in the standard context, we introduce all necessary elements of exterior calculus such as forms, vector fields, exterior derivatives, etc. We point out the relationships between these elements with the known ones in (local) exterior calcu- lus, discrete exterior calculus, etc.
    [Show full text]
  • An Introduction to Space–Time Exterior Calculus
    mathematics Article An Introduction to Space–Time Exterior Calculus Ivano Colombaro 1,* , Josep Font-Segura 1 and Alfonso Martinez 1 1 Department of Information and Communication Technologies, Universitat Pompeu Fabra, 08018 Barcelona, Spain; [email protected] (J.F.-S.); [email protected] (A.M.) * Correspondence: [email protected]; Tel.: +34-93-542-1496 Received: 21 May 2019; Accepted: 18 June 2019; Published: 21 June 2019 Abstract: The basic concepts of exterior calculus for space–time multivectors are presented: Interior and exterior products, interior and exterior derivatives, oriented integrals over hypersurfaces, circulation and flux of multivector fields. Two Stokes theorems relating the exterior and interior derivatives with circulation and flux, respectively, are derived. As an application, it is shown how the exterior-calculus space–time formulation of the electromagnetic Maxwell equations and Lorentz force recovers the standard vector-calculus formulations, in both differential and integral forms. Keywords: exterior calculus, exterior algebra, electromagnetism, Maxwell equations, differential forms, tensor calculus 1. Introduction Vector calculus has, since its introduction by J. W. Gibbs [1] and Heaviside, been the tool of choice to represent many physical phenomena. In mechanics, hydrodynamics and electromagnetism, quantities such as forces, velocities and currents are modeled as vector fields in space, while flux, circulation, divergence or curl describe operations on the vector fields themselves. With relativity theory, it was observed that space and time are not independent but just coordinates in space–time [2] (pp. 111–120). Tensors like the Faraday tensor in electromagnetism were quickly adopted as a natural representation of fields in space–time [3] (pp.
    [Show full text]