Spectral Element Methods (SEM) May 12, D Wave Equations (Strong and Weak form)
|
|
- Felicity Rich
- 5 years ago
- Views:
Transcription
1 Spectral Element Methods (SEM) May 12, D Wave Equations (Strong and Weak form) Strong form (PDEs): ρ 2 t s = T + F in (1) Multiply a global test function W(x) to both sides and integrate over to obtain the weak form: ρw s i dv = W ˆn T i dσ WT i d + WF i d (2) 2 Meshing Σ Grid into spectral elements (iteration index ispec). Each element will have GLL grid points in 3 directions (index α,β,γ). Globally number all grid points (index I) and establish the proection (α,β,γ; ispec) (I) (3) Notice this proection is one-to-one for grid points inside any element (valence = 0), and multipleto-one for grid points on the boundaries that are shared by elements (valence 1). 2.1 Interpolation of Shapes For boundary and volume elements, we not only need to access location of grid points of the element, but also the location of any arbitrary point inside the element. We can interpolate the shape of the element by given anchors x a. 1. For boundary element: x(ξ,η) = N a=1 N a (ξ,η)x a (4) and the Jacobian matrix associated with this transformation J b (ξ,η) = ( x ξ x η ) (5) 1
2 and the scalar Jacobian J b = J b. The inverse transform is given by ( ξ x η x ) = J 1 b and the normal to the boundary is described by (6) n(ξ,η) = x x ξ η x x ξ η (7) 2. For volume element: x(ξ,η,ζ) = N a=1 N a (ξ,η,ζ)x a (8) J a (ξ,η,ζ) = ( x(ξ,η,ζ) (ξ,η,ζ) ) (9) J a = J a (10) (ξ,η,ζ) x(ξ,η,ζ) = J 1 a N a (ξ,η) and N a (ξ,η,ζ) are the 2-D and 3-D shape functions. They are double or triple products of degree 1 or 2 Lagrangian polynomials. For example, degree 1 Lagrangian polynomials are N 1 (ξ) = 1(1 + ξ) and N 2 2(ξ) = 1 (1 ξ) for given anchors at 1 and Interpolation of Function Field We interpolate function field on GLL points that satisfy: (11) (1 ξ 2 ) P N (ξ) = 0 (12) and f(x(ξ,η)) Σe = αβ f(x(ξ,η,ζ)) e = αβγ f αβ l α (ξ)l β (η) (13) f αβγ l α (ξ)l β (η)l γ (ζ) 3.1 Derivatives i f(x(ξ,η,ζ)) = αβγ [ f αβγ l α (ξ)l β (η)l γ (ζ) ξ + l α (ξ) l β (η)l γ (ζ) η + l α (ξ)l β (η) l γ (ζ) ζ ] x i x i x i 2
3 Notice ξ η ζ x i, x i, x i values can be calculated from the 2D and 3D shape functions. Specifically, to compute derivatives on GLL points: [ ] [ ] [ ] i f αβγ = f βγ l (ξ α ) αβγ i ξ+ f αγ l (η β ) i η αβγ + f αβ l (ζ γ ) i ζ αβγ (15) (14) 3.2 Integration Integration quadrature f(x)dx = Σ e αβ e f(x)dx = αβγ ω α ω β f αβ J αβ b (16) ω α ω β ω γ f αβγ J αβγ a (17) 4 Global Test Functions Define global test functions W I (x), such that { W I l α (ξ)l β (η)l γ (ζ) if I e, and I e = (α,β,γ) (x) e = (18) 0 if I e If I e = (α,β,γ) has zero valence, then W I (x) is simply a 3-D local Lagrangian function extended to the whole space; otherwise it consists several pieces (valence+1) of local lagrangian function (edge ones). Therefore, for integration, one needs to loop over spectral elements, then all the GLL points, and adds contributions to the corresponding global grid point. Notice W I e (α,β,ζ ) = δ αα δ ββ δ ζζ (19) this simplifies the integration results. 3
4 5 Application to the Wave Equation 5.1 LHS ρw I s i (t) dx = ρw I e s i (t) dx (20) e e = ω α ω β ω γ ρ α β γ J α β γ s α β γ i δ αα δ ββ δ ζζ (21) e α β γ = ω α ω β ω γ ρ αβγ αβγ J αβγ s i (t) (22) e 5.2 Spatial Derivatives For example, i s (x) αβγ e = [ = i [ s α β γ α β γ l α (ξ)l β (η)l γ (ζ) s βγ l (ξ α )] i ξ(ξ α,η β,ζ γ ) ] αβγ e + [ + [ s αγ l (η β )] i η(ξ α,η β,ζ γ ) s αβ l (ζ γ )] i ζ(ξ α,η β,ζ γ ) (23) in the SEM code, hprime(α,) = l (ξ α ). The stress tensor is given by: T kl = C kli i s R kl (24) For isotropic and no pre-stress case, T kl = (κ 2 3 µ)δ klǫ ii + µ(ǫ kl + ǫ lk ) (25) 5.3 Boundary terms Σ W Iˆn T i dσ = e = e Σ e W I t i dσ e ω α ω β t αβ i J αβ (26) Where Ith global grid point is on Σ e elements, with GLL index (α,β). Normal stress t i is zero for free surface, and continuous for any internal boundaries. (For fluid-solid coupling problems, 4
5 we may solve both side independently, and the interchange of two fields will occur through this boundary condition.) For regional problems, other boundaries are absorbing, where t i = ˆn T i = ρα(ˆn ṡ )ˆn i + ρβ(δ i ˆn iˆn )ṡ (27) 5.4 Volume integral terms W I T i d = W I e T i d e e e = [ l α ξl β l γ + l α lβ ηl γ + l α l β lγ ζ]t i dx e e = {[ ] ω lα (ξ ) ξ(ξ,η β,ζ γ )T βγ i J βγ ω β ω γ e [ ] + ω lβ (η ) η(ξ α,η,ζ γ )T αγ i J αγ ω α ω γ + [ ω lγ (ζ ) ζ(ξ α,η β,ζ )T αβ i J αβ ] ω α ω β } (28) in the SEM code, sums over are named as temp[x y z][1 2 3]. 5.5 Forcing terms 1. For point force F i (x,t) = f i (t)δ(x x 0 ), W I F i d = W I (x 0 )f i (t) f i (t) if x 0 I = 0 if x 0 other grid point (29) l α (ξ(x 0 ))l β (ξ(x 0 ))l γ (ξ(x 0 ))f i (t) else x 0 e if x 0 has local parameter {ispec;ξ,η,ζ}, then compute the Lagrange polynomial at the source location. 2. For Moment tensor point force F i (x,t) = M i δ(x x 0 )g(t), W I F i d = M i W I (x 0 )g(t) (30) 5
6 Since W = W x and define G ik (ξ,η,ζ) = M i x, M i W I (x 0 ) = G ik W (ξ 0 ) = r,t,v = r,t,v l r (ξ 0 )l t (η 0 )l v (ζ 0 )G ik (ξ r,η t,ζ v ) ξk (l α (ξ r )l β (η t )l γ (ζ v )) l r (ξ 0 )l t (η 0 )l v (ζ 0 ) [G i1 (ξ r,η t,ζ v )l α(ξ r )l β (η t )l γ (ζ v )) + ] (31) 3. For body force field F i (x,t), W I F i d = e = e e W I e F i d e ω α ω β ω γ F αβγ i (t)j αβγ (32) 4. For surface force field F i (x,t) = Σ 0 τ i (x,t)δ(x x 0 )dσ 0, (monopoles on Σ 0 ) [ ] W I F i d = W I (x) τ i (x,t)δ(x x 0 )dσ 0 d (33) Σ 0 = τ i (x 0,t)W I (x 0 )dσ 0 (34) Σ 0 to compute this in practice, loop over surface elements, and then for each surface GLL point I = (α,β), Contr. to I th test function + = ω α ω β τ αβ J αβ (35) 5. For surface double couple force field F p (x,t) = Σ 0 u i (x,t)n C ipq q δ(x x 0 )dσ 0, assume that Σ 0 do not overlap with the surface Σ that encompasses, then [ ] W I F p d = W I (x) u i (x,t)n C ipq q δ dσ 0 d Σ [ 0 ] = W I (x)u i (x,t)n C ipq q δ(x x 0 )d W I (x) Σ 0 ( = q ui (t)n C ipq W I) (x 0 )dσ 0 (36) Σ 0 We define the corresponding moment density tensor m pq (t) = u i (t)n C ipq on Σ 0. In the case of a fault in an isotropic medium where u is perpendicular to n, m pq = µ(u p n q + u q n p ) is the standard moment density tensor that describes a seismic source. Now rewrite the above 6
7 expression: W I F i d = [ q m pq (t)w I + m pq(t) q W I ]dσ 0 Σ 0 (37) In practice, we store f p (t) = q m pq and m pq (t) in advance for each grid point of the surface Σ 0, then the above two terms can be calculated with ease in runtime. 5.6 Assemble Both Sides M s(t) = B + T + F(t) (38) We can update the acceleration at time t by s(t) = (M) 1 (B + T + F(t)) (39) Notice that M is diagonal (I I), B,T,F, s(t) are all vectors (I 1). Try loop over ispec and then (α,β,γ) to assemble the contribution of each grid point (in an element) into global matrices. This is true for all matrices, except the forcing term in which the contribution is added by looping over only elements that contain the forces and related GLL points. 6 Time Marching Schemes Newmark scheme Predictor: d n+1 = d n + v n t an ( t) 2 v n+1 = v n an t a n+1 = 0 (40) Corrector: a n+1 = (M) 1 (B + T + F n+1 ) v n+1 = v n an+1 t d n+1 = d n+1 (41) 7
The Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
More information8 A pseudo-spectral solution to the Stokes Problem
8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation / 59 Agenda Agenda
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationMethod in 1D. Introduction of the Basic Concepts and Comparison to Optimal FD Operators
The Spectral Element Method in 1D Introduction of the Basic Concepts and Comparison to Optimal FD Operators Seismology and Seismics Seminar 17.06.2003 Ludwigs-Maximilians-Universität München by Bernhard
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationFSI with Application in Hemodynamics Analysis and Simulation
FSI with Application in Hemodynamics Analysis and Simulation Mária Lukáčová Institute of Mathematics, University of Mainz A. Hundertmark (Uni-Mainz) Š. Nečasová (Academy of Sciences, Prague) G. Rusnáková
More information1. Basic Operations Consider two vectors a (1, 4, 6) and b (2, 0, 4), where the components have been expressed in a given orthonormal basis.
Questions on Vectors and Tensors 1. Basic Operations Consider two vectors a (1, 4, 6) and b (2, 0, 4), where the components have been expressed in a given orthonormal basis. Compute 1. a. 2. The angle
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationSPECIAL RELATIVITY AND ELECTROMAGNETISM
SPECIAL RELATIVITY AND ELECTROMAGNETISM MATH 460, SECTION 500 The following problems (composed by Professor P.B. Yasskin) will lead you through the construction of the theory of electromagnetism in special
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationImplementation of the dg method
Implementation of the dg method Matteo Cicuttin CERMICS - École Nationale des Ponts et Chaussées Marne-la-Vallée March 6, 2017 Outline Model problem, fix notation Representing polynomials Computing integrals
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationPolynomial form of the Hilbert Einstein action
1 Polynomial form of the Hilbert Einstein action M. O. Katanaev Steklov Mathematical Institute, Gubkin St. 8, Moscow, 119991, Russia arxiv:gr-qc/0507026v1 7 Jul 2005 6 July 2005 Abstract Configuration
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationShape Function Generation and Requirements
Shape Function Generation and Requirements Requirements Requirements (A) Interpolation condition. Takes a unit value at node i, and is zero at all other nodes. Requirements (B) Local support condition.
More informationMath. 460, Sec. 500 Fall, Special Relativity and Electromagnetism
Math. 460, Sec. 500 Fall, 2011 Special Relativity and Electromagnetism The following problems (composed by Professor P. B. Yasskin) will lead you through the construction of the theory of electromagnetism
More informationseismic wave propagation through spherical sections including anisotropy and visco-elastic dissipation 3D heterogeneities of all Earth model
seismic wave propagation through spherical sections including anisotropy and visco-elastic dissipation 3D heterogeneities of all Earth model properties inversion of seismic observables for 3D Earth structure
More informationWRT in 2D: Poisson Example
WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx
More informationHomework 7-8 Solutions. Problems
Homework 7-8 Solutions Problems 26 A rhombus is a parallelogram with opposite sides of equal length Let us form a rhombus using vectors v 1 and v 2 as two adjacent sides, with v 1 = v 2 The diagonals of
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationSpinor Representation of Conformal Group and Gravitational Model
Spinor Representation of Conformal Group and Gravitational Model Kohzo Nishida ) arxiv:1702.04194v1 [physics.gen-ph] 22 Jan 2017 Department of Physics, Kyoto Sangyo University, Kyoto 603-8555, Japan Abstract
More informationSimulating Thin Shells with MPM
Simulating Thin Shells with MPM Biswajit Banerjee Center for the Simulation of Accidental Fires and xplosions University of Utah March 14, 2005 Outline The Problem. Shell Theory. MPM Formulation. Results.
More informationEnergy and potential enstrophy flux constraints in quasi-geostrophic models
Energy and potential enstrophy flux constraints in quasi-geostrophic models Eleftherios Gkioulekas University of Texas Rio Grande Valley August 25, 2017 Publications K.K. Tung and W.W. Orlando (2003a),
More informationMathematics for Engineers. Numerical mathematics
Mathematics for Engineers Numerical mathematics Integers Determine the largest representable integer with the intmax command. intmax ans = int32 2147483647 2147483647+1 ans = 2.1475e+09 Remark The set
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More informationAlgorithms for polynomial multidimensional spectral factorization and sum of squares
Algorithms for polynomial multidimensional spectral factorization and sum of squares H.L. Trentelman University of Groningen, The Netherlands. Japan November 2007 Algorithms for polynomial multidimensional
More informationMA8502 Numerical solution of partial differential equations. The Poisson problem: Mixed Dirichlet/Neumann boundary conditions along curved boundaries
MA85 Numerical solution of partial differential equations The Poisson problem: Mied Dirichlet/Neumann boundar conditions along curved boundaries Fall c Einar M. Rønquist Department of Mathematical Sciences
More informationLet M be a Riemannian manifold with Levi-Civita connection D. For X, Y, W Γ(TM), we have R(X, Y )Z = D Y D X Z D X D Y Z D [X,Y ] Z,
Let M be a Riemannian manifold with Levi-Civita connection D. For X, Y, W Γ(TM), we have R(X, Y )Z = D Y D X Z D X D Y Z D [X,Y ] Z, where R is the curvature tensor of D. We have Proposition 1. For any
More informationSpectral-Element and Adjoint Methods in Seismology
Geophys. J. Int. (), Spectral-Element and Adjoint Methods in Seismology Jeroen Tromp Seismological Laboratory, California Institute of Technology, Pasadena, California 91125, USA 13 September 26 SUMMARY
More informationChapter 7. Kinematics. 7.1 Tensor fields
Chapter 7 Kinematics 7.1 Tensor fields In fluid mechanics, the fluid flow is described in terms of vector fields or tensor fields such as velocity, stress, pressure, etc. It is important, at the outset,
More informationThe effect of asymmetric large-scale dissipation on energy and potential enstrophy injection in two-layer quasi-geostrophic turbulence
The effect of asymmetric large-scale dissipation on energy and potential enstrophy injection in two-layer quasi-geostrophic turbulence Eleftherios Gkioulekas (1) and Ka-Kit Tung (2) (1) Department of Mathematics,
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationA framework for monetary-policy analysis Overview Basic concepts: Goals, targets, intermediate targets, indicators, operating targets, instruments
Eco 504, part 1, Spring 2006 504_L6_S06.tex Lars Svensson 3/6/06 A framework for monetary-policy analysis Overview Basic concepts: Goals, targets, intermediate targets, indicators, operating targets, instruments
More informationHamiltonian Formulation of Piano String Lagrangian Density with Riemann-Liouville Fractional Definition
Hamiltonian Formulation of Piano String Lagrangian Density with Riemann-Liouville Fractional Definition Emad K. Jaradat 1 1 Department of Physics, Mutah University, Al-Karak, Jordan Email: emad_jaradat75@yahoo.com;
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationWave and Elasticity Equations
1 Wave and lasticity quations Now let us consider the vibrating string problem which is modeled by the one-dimensional wave equation. Suppose that a taut string is suspended by its extremes at the points
More informationCOURSE Numerical integration of functions
COURSE 6 3. Numerical integration of functions The need: for evaluating definite integrals of functions that has no explicit antiderivatives or whose antiderivatives are not easy to obtain. Let f : [a,
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationAn Introduction to Kaluza-Klein Theory
An Introduction to Kaluza-Klein Theory A. Garrett Lisi nd March Department of Physics, University of California San Diego, La Jolla, CA 993-39 gar@lisi.com Introduction It is the aim of Kaluza-Klein theory
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 6124 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #6a due Thursday, October 25, 2012 Notes for lectures 14 and 15: Calculus on smooth
More informationMCE503: Modeling and Simulation of Mechatronic Systems. Modeling of Continuous Beams: Finite-Mode Approach
MCE503: Modeling and Simulation of Mechatronic Systems MCE503 p./23 Modeling of Continuous Beams: Finite-Mode Approach Reading: KMR Chapter 0 Cleveland State University Mechanical Engineering Hanz Richter,
More informationVectors. Three dimensions. (a) Cartesian coordinates ds is the distance from x to x + dx. ds 2 = dx 2 + dy 2 + dz 2 = g ij dx i dx j (1)
Vectors (Dated: September017 I. TENSORS Three dimensions (a Cartesian coordinates ds is the distance from x to x + dx ds dx + dy + dz g ij dx i dx j (1 Here dx 1 dx, dx dy, dx 3 dz, and tensor g ij is
More informationMidterm Exam. Math 460, Fall October 2015
Midterm Exam Math 460, Fall 2015 14 October 2015 1. This exam has 4 questions and 6 pages for a total of 110 points. Make sure you have all pages before you begin. 2. As a form of extra credit, the grade
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationCohomology and Vector Bundles
Cohomology and Vector Bundles Corrin Clarkson REU 2008 September 28, 2008 Abstract Vector bundles are a generalization of the cross product of a topological space with a vector space. Characteristic classes
More informationPartititioned Methods for Multifield Problems
Partititioned Methods for Multifield Problems Joachim Rang, 22.7.215 22.7.215 Joachim Rang Partititioned Methods for Multifield Problems Seite 1 Non-matching matches usually the fluid and the structure
More information10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections )
c Dr. Igor Zelenko, Fall 2017 1 10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections 7.2-7.4) 1. When each of the functions F 1, F 2,..., F n in right-hand side
More informationModel order reduction of mechanical systems subjected to moving loads by the approximation of the input
Model order reduction of mechanical systems subjected to moving loads by the approximation of the input Alexander Vasilyev, Tatjana Stykel Institut für Mathematik, Universität Augsburg ModRed 2013 MPI
More informationMA106 Linear Algebra lecture notes
MA106 Linear Algebra lecture notes Lecturers: Diane Maclagan and Damiano Testa 2017-18 Term 2 Contents 1 Introduction 3 2 Matrix review 3 3 Gaussian Elimination 5 3.1 Linear equations and matrices.......................
More informationLecture notes: GE 263 Computational Geophysics The spectral element method
Lecture notes: GE 263 Computational Geophysics The spectral element method Jean-Paul Ampuero Abstract The spectral element method (SEM) is a high order numerical method for solving partial differential
More informationBackground. Background. C. T. Kelley NC State University tim C. T. Kelley Background NCSU, Spring / 58
Background C. T. Kelley NC State University tim kelley@ncsu.edu C. T. Kelley Background NCSU, Spring 2012 1 / 58 Notation vectors, matrices, norms l 1 : max col sum... spectral radius scaled integral norms
More informationf dr. (6.1) f(x i, y i, z i ) r i. (6.2) N i=1
hapter 6 Integrals In this chapter we will look at integrals in more detail. We will look at integrals along a curve, and multi-dimensional integrals in 2 or more dimensions. In physics we use these integrals
More informationFinite Element simulations of a phase-field model for mode-iii fracture. by V. Kaushik
Finite Element simulations of a phase-field model for mode-iii fracture. by V. Kaushik 1 Motivation The motivation behind this model is to understand the underlying physics behind branching of cracks during
More informationMathematical Background
CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationGeneral Curvilinear Ocean Model (GCOM): Enabling Thermodynamics
General Curvilinear Ocean Model (GCOM): Enabling Thermodynamics M. Abouali, C. Torres, R. Walls, G. Larrazabal, M. Stramska, D. Decchis, and J.E. Castillo AP0901 09 General Curvilinear Ocean Model (GCOM):
More informationElectricity & Magnetism Qualifier
Electricity & Magnetism Qualifier For each problem state what system of units you are using. 1. Imagine that a spherical balloon is being filled with a charged gas in such a way that the rate of charge
More informationParticle-based Fluids
Particle-based Fluids Particle Fluids Spatial Discretization Fluid is discretized using particles 3 Particles = Molecules? Particle approaches: Molecular Dynamics: relates each particle to one molecule
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 34 Outline 1 Infinite Order Finite Differences
More informationLecture 35: The Inertia Tensor
Lecture 35: The Inertia Tensor We found last time that the kinetic energy of a rotating obect was: 1 Trot = ωω i Ii where i, ( I m δ x x x i i, k, i, k So the nine numbers represented by the I i tell us
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More information1.8. Zeroth-Order Tensors (Scalars). A scalar is a single function (i.e., one component) which is invariant under
1.8. Zeroth-Order Tensors (Scalars). A scalar is a single function (i.e., one component) which is invariant under changes of the coordinate systems. We deal with rectangular coordinate systems only. Thus
More informationFinite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.
More informationThe second-order 1D wave equation
C The second-order D wave equation C. Homogeneous wave equation with constant speed The simplest form of the second-order wave equation is given by: x 2 = Like the first-order wave equation, it responds
More informationAS Level / Year 1 Edexcel Further Maths / CP1
AS Level / Year 1 Edexcel Further Maths / CP1 March 2018 Mocks 2018 crashmaths Limited 1 (a) n k(2k 3) = 2 k 2 k=1 k=1 n n Use linearity 3 k k=1 n k=1 k(2k 3) = 2 n 6 (n +1)(2n +1) 3 n 2 (n +1) Uses standard
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationMechanics Physics 151
Mechanics Physics 151 Lecture 4 Continuous Systems and Fields (Chapter 13) What We Did Last Time Built Lagrangian formalism for continuous system Lagrangian L Lagrange s equation = L dxdydz Derived simple
More informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
More informationOn the hydrodynamic diffusion of rigid particles
On the hydrodynamic diffusion of rigid particles O. Gonzalez Introduction Basic problem. Characterize how the diffusion and sedimentation properties of particles depend on their shape. Diffusion: Sedimentation:
More informationPhysics of Continuous media
Physics of Continuous media Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 October 26, 2012 Deformations of continuous media If a body is deformed, we say that the point which originally had
More informationOn divergence-free reconstruction schemes for CED and MHD
On divergence-free reconstruction schemes for CED and MHD Praveen Chandrashekar praveen@math.tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore-560065, India
More information1 Ordinary Differential Equations
Ordinary Differential Equations.0 Mathematical Background.0. Smoothness Definition. A function f defined on [a, b] is continuous at ξ [a, b] if lim x ξ f(x) = f(ξ). Remark Note that this implies existence
More informationFDTD for 1D wave equation. Equation: 2 H Notations: o o. discretization. ( t) ( x) i i i i i
FDTD for 1D wave equation Equation: 2 H = t 2 c2 2 H x 2 Notations: o t = nδδ, x = iδx o n H nδδ, iδx = H i o n E nδδ, iδx = E i discretization H 2H + H H 2H + H n+ 1 n n 1 n n n i i i 2 i+ 1 i i 1 = c
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationCosmology. April 13, 2015
Cosmology April 3, 205 The cosmological principle Cosmology is based on the principle that on large scales, space (not spacetime) is homogeneous and isotropic that there is no preferred location or direction
More informationAxioms of Adaptivity (AoA) in Lecture 2 (sufficient for optimal convergence rates)
Axioms of Adaptivity (AoA) in Lecture 2 (sufficient for optimal convergence rates) Carsten Carstensen Humboldt-Universität zu Berlin 2018 International Graduate Summer School on Frontiers of Applied and
More informationExtra Problems and Examples
Extra Problems and Examples Steven Bellenot October 11, 2007 1 Separation of Variables Find the solution u(x, y) to the following equations by separating variables. 1. u x + u y = 0 2. u x u y = 0 answer:
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationTopology Optimization Using the SIMP Method
Fabian Wein Introduction Introduction About this document This is a fragment of a talk given interally Intended for engineers and mathematicians SIMP basics Detailed introduction (based on linear elasticity)
More informationSpectral Elements. Introduction recalling the elastic wave equation
Spectral Elements Introduction recalling the elastic wave equation The spectral-element method: General concept domain mapping from space-continuous to space-discrete time extrapolation Gauss-Lobatto-Legendre
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationMaria Cameron Theoretical foundations. Let. be a partition of the interval [a, b].
Maria Cameron 1 Interpolation by spline functions Spline functions yield smooth interpolation curves that are less likely to exhibit the large oscillations characteristic for high degree polynomials Splines
More informationA STABLE PENALTY METHOD FOR THE COMPRESSIBLE NAVIER STOKES EQUATIONS: III. MULTIDIMENSIONAL DOMAIN DECOMPOSITION SCHEMES
SIAM J. SCI. COMPUT. Vol. 20, No., pp. 62 93 c 998 Society for Industrial and Applied Mathematics A STABLE PENALTY METHOD FOR THE COMPRESSIBLE NAVIER STOKES EQUATIONS: III. MULTIDIMENSIONAL DOMAIN DECOMPOSITION
More informationPhysics 325: General Relativity Spring Final Review Problem Set
Physics 325: General Relativity Spring 2012 Final Review Problem Set Date: Friday 4 May 2012 Instructions: This is the third of three review problem sets in Physics 325. It will count for twice as much
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationCurved Spacetime I. Dr. Naylor
Curved Spacetime I Dr. Naylor Last Week Einstein's principle of equivalence We discussed how in the frame of reference of a freely falling object we can construct a locally inertial frame (LIF) Space tells
More informationComputational Methods in Optimal Control Lecture 8. hp-collocation
Computational Methods in Optimal Control Lecture 8. hp-collocation William W. Hager July 26, 2018 10,000 Yen Prize Problem (google this) Let p be a polynomial of degree at most N and let 1 < τ 1 < τ 2
More information