MA8502 Numerical solution of partial differential equations. The Poisson problem: Mixed Dirichlet/Neumann boundary conditions along curved boundaries
|
|
- Lindsay Montgomery
- 5 years ago
- Views:
Transcription
1 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 NTNU, N-749 Trondheim, Norwa All rights reserved
2 Model problem ˆ Γ 4 Γ F ˆΓ 4 ˆΓ Γ ˆΓ Figure : The phsical domain = (, L ) (, L ) can be viewed as coming from the reference domain = (, ) (, ) through the one-to-one mapping F :. We consider here the numerical solution of the two-dimensional Poisson equation u = f in, with the following mied boundar conditions specified along : u Γ = { Dirichlet u n = u { Homogeneous Γ n = Neumann Γ4 conditions { u Non-homogeneous n = g() Γ Neumann Here, u n = u n is the derivative of u in the outward normal direction. We first formulate the corresponding weak formulation. Let X be the function space where X = {v H () v Γ = }. }{{} Onl Dirichlet conditions The weak form is then: Find u X such that u v d d = fv d d + ( u n)v ds = 4 k= Γ k u n v ds = ( u n) v ds v X, gv ds v X.
3 where We ma also formulate this problem abstractl as: Find u X such that a(u, v) = l(v) v X, a(u, v) = l(v) = u v d d, fv d d + gv ds. Γ } {{} Non-zero Neumann conditions We now proceed with a discretization based on high order polnomials. We first define the discrete space X N = {v X v F P()}. The discrete problem (with quadrature) is then: Find u N X N such that a N (u N, v) = l N (v) v X N, where a N (, ) and l N ( ) denote evaluation of the bilinear form and linear form b Gauss- Lobatto Legendre (GLL) quadrature. The net step is to derive a sstem of algebraic equations. To this end, we have to define a basis for our discrete space, X N : v X N, v F = ˆv, ˆv(ξ, η) = m= n= v mn l m (ξ) l n (η) where l m (ξ), m =,..., N are the usual one-dimensional Lagrangian interpolants through the GLL points. We remark that dim(x N ) = N(N + ). With this nodal, tensor-product basis, we can epress the numerical solution (in the reference variables) as Note that, since u Γ =, û N (ξ, η) = u mn l m (ξ)l n (η). () m= n= As usual, we choose test functions û N (ξ, ) = = u m =, m =,..., N. ˆv(ξ, η) = l i (ξ)l j (η) i =,..., N j =,..., N. ()
4 Using () and (), we are now read to give an eplicit epression for the bilinear form a N (u N, v). First, we note that ( L û N ˆv a(u N, v) = u N v d d = L ξ ξ + L ) û N ˆv dξ dη. L η η Using () and appling GLL quadrature, we arrive at (for this particular test function) with a N (u N, v) =  pq = m= n=  im ˆBjn + L ) ˆBim  jn u mn, L L ( L ρ α D αp D αq (the one-dimensional Laplace operator on (, )), α= ˆB pq = ρ p δ pq (the one-dimensional mass matri on (, )). In a similar fashion, we can derive an eplicit epression for the linear form l N (v). First, we note that l(v) = fv d d + gv ds = L L ˆf ˆv dξ dη + L ĝ ˆv dξ. 4 ˆ Again, using () and appling GLL quadrature, we arrive at l N (v) = L L 4 We note that α= β= ρ α ρ β f αβ l i (ξ α )l j (ξ β ) + L ρ α g α l i (ξ α )l j () α= = L L ρ i ρ j f ij + L 4 ρ ig i δ jn, i =,..., N, j =,..., N. }{{} onl contribution along J = L L 4 J S = L f ij = f((ξ i ), (ξ j )) g i = g((ξ i )) : Jacobian (volume) : Jacobian (surface) 4
5 ˆ Γ 4 Γ F ˆΓ 4 ˆΓ Γ ˆΓ Model problem The same problem as model problem, ecept that is now deformed. The same boundar conditions. The same abstract formulation. The domain is approimated using an isoparametric representation. The bilinear form a(u N, v) can now be epressed as a(u N, v) = ( ˆ ˆv) T G ˆ ûn dξ dη (see earlier notes), and a N (u N, v) is again obtained using GLL quadrature when evaluating the integral. The linear form l(v) is given as l(v) = fv d d + gv ds Γ = ˆf ˆvJ dξ dη + ĝ ˆvJ S dξ, ˆ where J = ξ J S = s ξ η η ξ Jacobian (volume), Jacobian (surface). s 5
6 Setting ˆv(ξ, η) = l i (ξ)l j (η) and appling GLL quadrature gives l N (v) = = α= β= α= β= ρ α ρ β J αβ f αβ l i (ξ α )l j (ξ β ) + ρ α ρ β J αβ f αβ δ iα δ jβ + ρ α Jα S g α l i (ξ α )l j () α= ρ α Jα S g α δ iα δ jn α= = ρ i ρ j J ij f ij + ρ i Ji S g i δ jn }{{} onl, i =,..., N, j =,..., N. contribution along Computationall, we let N P N (), N P N (), J = N N ξ η N η (( J S = ds dξ = N ξ N ξ, ˆΓ) + ( N ξ ) / ˆΓ). s ˆ Note also that f ij = f((ξ i, ξ j ), (ξ i, ξ j )), g i = g(s i ) = g(s(ξ i )), where (ξ i, ξ j ) = ij (ξ i, ξ j ) = ij } computed using Gordon-Hall s si coordinates ( in, in ) 6
7 Numerical results Let us conclude this section with some actual numerical results. We consider the solution of the Laplace equation in a 9 degree bend; see Figure. Γ 4 Γ Γ Figure : The computational domain corresponds to an aismmetric bend. Two of the edges correspond to a quarter of two concentric circles with radii r = and r =. We solve the Laplace equation (i.e., f = ) using a spectral approimation based on high order polnomials. The boundar conditions along Γ i, i =,,, 4 are as described in Model problem and Model problem. In order to approimate the geometr b high order polnomials, we first distribute N + GLL points along each of the four edges. The points are distributed radiall along Γ and Γ 4, and in an angular fashion along Γ and (i.e., according to arc length). Once the boundar points are distributed, we generate all the interior points using the Gordon-Hall algorithm described earlier. This will give us the phsical coordinates corresponding to all the GLL points on the reference domain = (, ). Figure and Figure 4 show the intermediate stages in this construction, as well as final set of points { mn } N m,n=. We can now compute N corresponding to an point (ξ, η) via our high order tensor-product basis N (ξ, η) = m= n= mn l m (ξ)l n (η). 7
8 Figure : The Gordon-Hall algorithm. Left plot: preserving Γ and (the mapping F ξ ). Right plot: preserving Γ and Γ 4 (the mapping F η ) Figure 4: The Gordon-Hall algorithm. Left plot: preserving the four corners (the mapping F ξη ). Right plot: the final mapping F = F ξ + F η F ξη. 8
9 We first impose a constant flu g = along the boundar. The numerical solution is depicted in Figure 5. As epected, we observe eponential convergence for this problem with regular data and geometr (represented using an isoparametric mapping) Error N Figure 5: Numerical results when imposing a constant flu g = along. Left plot: contour lines of the numerical solution, u N, to the Laplace problem. Right plot: the error u u N in the discrete L -norm as a function of the polnomial degree, N, used in the numerical approimation. Net, we impose a trigonometric variation in the heat flu g along the boundar. The numerical solution is depicted in Figure Figure 6: Numerical results when imposing a variable flu, g, along. Contour lines of the numerical solution, u N, to the Laplace problem. 9
10 Eercise Consider the curve Γ given b the function () = cos() for [, π]. We would like to compute the length L of this curve. The length can be epressed as L = ds, where ds is an infinitesimal length element along the curve Γ. Γ. Suggest a reparametrization of the curve on the form ((ξ), (ξ)), for ξ [, ]. Here, ξ = corresponds to the left end point of the curve, i.e. the point (, ), while ξ = corresponds to the right end point of the curve, i.e., the point (π, ).. Use the reparametrization to epress L as an integral with respect to ξ instead of s.. Compute the integral numericall using GLL quadrature.
Applications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
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 informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
More informationChapter 11 Three-Dimensional Stress Analysis. Chapter 11 Three-Dimensional Stress Analysis
CIVL 7/87 Chapter - /39 Chapter Learning Objectives To introduce concepts of three-dimensional stress and strain. To develop the tetrahedral solid-element stiffness matri. To describe how bod and surface
More informationElliptic Equations. Chapter Definitions. Contents. 4.2 Properties of Laplace s and Poisson s Equations
5 4. Properties of Laplace s and Poisson s Equations Chapter 4 Elliptic Equations Contents. Neumann conditions the normal derivative, / = n u is prescribed on the boundar second BP. In this case we have
More informationEvaluation of Flux Integral Outputs for the Reduced Basis Method
Evaluation of Flux Integral Outputs for the Reduced Basis Method Jens L. Eftang and Einar M. Rønquist 20th October 2009 Abstract In this paper, we consider the evaluation of flux integral outputs from
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationPROBLEM SET 6. E [w] = 1 2 D. is the smallest for w = u, where u is the solution of the Neumann problem. u = 0 in D u = h (x) on D,
PROBLEM SET 6 UE ATE: - APR 25 Chap 7, 9. Questions are either directl from the text or a small variation of a problem in the text. Collaboration is oka, but final submission must be written individuall.
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 informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More information5.6. Differential equations
5.6. Differential equations The relationship between cause and effect in phsical phenomena can often be formulated using differential equations which describe how a phsical measure () and its derivative
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 informationLecture 3: Asymptotic Methods for the Reduced Wave Equation
Lecture 3: Asymptotic Methods for the Reduced Wave Equation Joseph B. Keller 1 The Reduced Wave Equation Let v (t,) satisfy the wave equation v 1 2 v c 2 = 0, (1) (X) t2 where c(x) is the propagation speed
More informationVibration Analysis of Isotropic and Orthotropic Plates with Mixed Boundary Conditions
Tamkang Journal of Science and Engineering, Vol. 6, No. 4, pp. 7-6 (003) 7 Vibration Analsis of Isotropic and Orthotropic Plates with Mied Boundar Conditions Ming-Hung Hsu epartment of Electronic Engineering
More information44 CHAPTER 3. CAVITY SCATTERING
44 CHAPTER 3. CAVITY SCATTERING For the TE polarization case, the incident wave has the electric field parallel to the x 3 -axis, which is also the infinite axis of the aperture. Since both the incident
More informationThe Control-Volume Finite-Difference Approximation to the Diffusion Equation
The Control-Volume Finite-Difference Approimation to the Diffusion Equation ME 448/548 Notes Gerald Recktenwald Portland State Universit Department of Mechanical Engineering gerr@mepdedu ME 448/548: D
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 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 informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
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 informationA Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations
Int. Conference on Boundary and Interior Layers BAIL 6 G. Lube, G. Rapin Eds c University of Göttingen, Germany, 6 A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations. Introduction
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationSOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4
SOLUTIONS TO SELECTED POBLEMS FOM ASSIGNMENTS 3, 4 Problem 5 from Assignment 3 Statement. Let be an n-dimensional bounded domain with smooth boundary. Show that the eigenvalues of the Laplacian on with
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationExercise solutions: concepts from chapter 5
1) Stud the oöids depicted in Figure 1a and 1b. a) Assume that the thin sections of Figure 1 lie in a principal plane of the deformation. Measure and record the lengths and orientations of the principal
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationNST1A: Mathematics II (Course A) End of Course Summary, Lent 2011
General notes Proofs NT1A: Mathematics II (Course A) End of Course ummar, Lent 011 tuart Dalziel (011) s.dalziel@damtp.cam.ac.uk This course is not about proofs, but rather about using different techniques.
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More information3.3.1 Linear functions yet again and dot product In 2D, a homogenous linear scalar function takes the general form:
3.3 Gradient Vector and Jacobian Matri 3 3.3 Gradient Vector and Jacobian Matri Overview: Differentiable functions have a local linear approimation. Near a given point, local changes are determined by
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationOn the Stability of the Flux Reconstruction Schemes on Quadrilateral Elements for the Linear Advection Equation
OI 10.1007/s10915-015-0102-8 On the Stability of the Flux Reconstruction Schemes on Quadrilateral Elements for the Linear Advection Equation Abhishe Sheshadri 1 Antony Jameson 1 Received: 25 March 2015
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
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 informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationBoundary value problems for the elliptic sine Gordon equation in a semi strip
Boundary value problems for the elliptic sine Gordon equation in a semi strip Article Accepted Version Fokas, A. S., Lenells, J. and Pelloni, B. 3 Boundary value problems for the elliptic sine Gordon equation
More informationFirst Order Equations
10 1 Linear and Semilinear Equations Chapter First Order Equations Contents 1 Linear and Semilinear Equations 9 Quasilinear Equations 19 3 Wave Equation 6 4 Sstems of Equations 31 1 Linear and Semilinear
More informationThe 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 informationINFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS
INFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS Michael Newman Department of Aerospace Engineering Texas A&M University 34 TAMU 744D H.R. Bright Building College Station,
More informationAsymptotics for evolution problems with nonlocal diffusion. Julio D. Rossi
Asymptotics for evolution problems with nonlocal diffusion Julio D. Rossi IMDEA Matemáticas, C-IX, Campus Cantoblanco, UAM, Madrid, Spain and Departamento de Matemática, FCEyN UBA (1428) Buenos Aires,
More informationCOMPUTING FRACTIONAL LAPLACIANS ON COMPLEX-GEOMETRY DOMAINS: ALGORITHMS AND SIMULATIONS
SIAM J. SCI. COMPUT. Vol. 39, No. 4, pp. A3 A344 c 7 Society for Industrial and Applied Mathematics COMPUTING FRACTIONAL LAPLACIANS ON COMPLEX-GEOMETRY DOMAINS: ALGORITHMS AND SIMULATIONS FANGYING SONG,
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationSteady and unsteady diffusion
Chapter 5 Steady and unsteady diffusion In this chapter, we solve the diffusion and forced convection equations, in which it is necessary to evaluate the temperature or concentration fields when the velocity
More information16.5. Maclaurin and Taylor Series. Introduction. Prerequisites. Learning Outcomes
Maclaurin and Talor Series 6.5 Introduction In this Section we eamine how functions ma be epressed in terms of power series. This is an etremel useful wa of epressing a function since (as we shall see)
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationNONLINEAR STRUCTURAL DYNAMICS USING FE METHODS
NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.
More informationOn explicit integration of two non-holonomic problems
On explicit integration of two non-holonomic problems Alexey V. Borisov 1 1 Institute of Computer Sciences, Izhevsk, Russia Generalized Chaplygin systems Equations of motion of the generalized Chaplygin
More informationarxiv:gr-qc/ v1 6 Sep 2006
Introduction to spectral methods Philippe Grandclément Laboratoire Univers et ses Théories, Observatoire de Paris, 5 place J. Janssen, 995 Meudon Cede, France This proceeding is intended to be a first
More informationA Method for Geometry Optimization in a Simple Model of Two-Dimensional Heat Transfer
A Method for Geometr Optimization in a Simple Model of Two-Dimensional Heat Transfer X. Peng, K. Niakhai B. Protas Jul, 3 arxiv:37.48v [math.na] 4 Jul 3 Abstract This investigation is motivated b the problem
More informationTMA4220: Programming project - part 1
TMA4220: Programming project - part 1 TMA4220 - Numerical solution of partial differential equations using the finite element method The programming project will be split into two parts. This is the first
More informationis the two-dimensional curl of the vector field F = P, Q. Suppose D is described by a x b and f(x) y g(x) (aka a type I region).
Math 55 - Vector alculus II Notes 4.4 Green s Theorem We begin with Green s Theorem: Let be a positivel oriented (parameterized counterclockwise) piecewise smooth closed simple curve in R and be the region
More informationSurface integrals, Divergence theorem of Gauss
c jbquig-ucd Februar 13, 3 Chapter 6 Surface integrals, Divergence theorem of Gauss introduction A solid D R 3 might be described b an inequalit, see chapter5. A surface R 3 is often described b a single
More information4.7. Newton s Method. Procedure for Newton s Method HISTORICAL BIOGRAPHY
4. Newton s Method 99 4. Newton s Method HISTORICAL BIOGRAPHY Niels Henrik Abel (18 189) One of the basic problems of mathematics is solving equations. Using the quadratic root formula, we know how to
More informationIn particular, if the initial condition is positive, then the entropy solution should satisfy the following positivity-preserving property
1 3 4 5 6 7 8 9 1 11 1 13 14 15 16 17 18 19 1 3 4 IMPLICIT POSITIVITY-PRESERVING HIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR CONSERVATION LAWS TONG QIN AND CHI-WANG SHU Abstract. Positivity-preserving
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 informationAircraft Structures Structural & Loading Discontinuities
Universit of Liège Aerospace & Mechanical Engineering Aircraft Structures Structural & Loading Discontinuities Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/
More informationFitting Integrands to Basic Rules. x x 2 9 dx. Solution a. Use the Arctangent Rule and let u x and a dx arctan x 3 C. 2 du u.
58 CHAPTER 8 Integration Techniques, L Hôpital s Rule, and Improper Integrals Section 8 Basic Integration Rules Review proceres for fitting an integrand to one of the basic integration rules Fitting Integrands
More informationJournal of Computational Physics
Journal of Computational Physics 23 (211) 8573 862 Contents lists available at SciVerse ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp An eigen-based high-order
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationWBHS Algebra 2 - Final Exam
Class: _ Date: _ WBHS Algebra 2 - Final Eam Multiple Choice Identify the choice that best completes the statement or answers the question. Describe the pattern in the sequence. Find the net three terms.
More informationFinite elements for plates and shells. Advanced Design for Mechanical System LEC 2008/11/04
Finite elements for plates and shells Sstem LEC 008/11/04 1 Plates and shells A Midplane The mid-surface of a plate is plane, a curved geometr make it a shell The thickness is constant. Ma be isotropic,
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
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 informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
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 informationChapter 4. Two-Dimensional Finite Element Analysis
Chapter 4. Two-Dimensional Finite Element Analysis general two-dimensional boundary-value problems 4.1 The Boundary-Value Problem 2nd-order differential equation to consider α α β φ Ω (4.1) Laplace, Poisson
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 informationA new class of pseudodifferential operators with mixed homogenities
A new class of pseudodifferential operators with mixed homogenities Po-Lam Yung University of Oxford Jan 20, 2014 Introduction Given a smooth distribution of hyperplanes on R N (or more generally on a
More informationAlf Emil Løvgren. Reduced basis modeling of hierarchical flow systems. Doctoral Theses at NTNU, 2005:212 Alf Emil Løvgren
Alf Emil Løvgren Doctoral Theses at NTNU, 2005:212 Alf Emil Løvgren Reduced basis modeling of hierarchical flow systems ISBN 82-471-7322-0 (printed ver.) ISBN 82-471-7321-2 (electronic ver.) ISSN 1503-8181
More informationAn anisotropic, superconvergent nonconforming plate finite element
Journal of Computational and Applied Mathematics 220 (2008 96 0 www.elsevier.com/locate/cam An anisotropic, superconvergent nonconforming plate finite element Shaochun Chen a,liyin a, Shipeng Mao b, a
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains
Department of Mathematics Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains Samira Nikkar and Jan Nordström LiTH-MAT-R--/--SE Department of Mathematics
More informationA Unstructured Nodal Spectral-Element Method for the Navier-Stokes Equations
Commun. Comput. Phys. doi: 10.4208/cicp.070111.140711a Vol. 12, o. 1, pp. 315-336 July 2012 A Unstructured odal Spectral-Element Method for the avier-stokes Equations Lizhen Chen 1, Jie Shen 1,2 and Chuanju
More informationFast and accurate methods for the discretization of singular integral operators given on surfaces
Fast and accurate methods for the discretization of singular integral operators given on surfaces James Bremer University of California, Davis March 15, 2018 This is joint work with Zydrunas Gimbutas (NIST
More informationTHE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER QUADRILATERAL MESHES. Zhong-Ci Shi and Xuejun Xu.
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 9, Number, January 2008 pp. 63 82 THE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER
More informationStress-strain relations
SICLLY INDRMIN SRSS SYSMS staticall determinate stress sstem simple eample of this is a bar loaded b a weight, hanging in tension. he solution for the stress is simpl W/ where is the cross sectional area.
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More informationFast tensor-product solvers. Part I: Partially deformed three-dimensional domains
Fast tensor-product solvers. Part I: Partially deformed three-dimensional domains Tormod Bjøntegaard, Yvon Maday, Einar M Rønquist To cite this version: Tormod Bjøntegaard, Yvon Maday, Einar M Rønquist.
More informationFitting Integrands to Basic Rules
6_8.qd // : PM Page 8 8 CHAPTER 8 Integration Techniques, L Hôpital s Rule, and Improper Integrals Section 8. Basic Integration Rules Review proceres for fitting an integrand to one of the basic integration
More informationPolynomial approximation and Splines
Polnomial approimation and Splines 1. Weierstrass approimation theorem The basic question we ll look at toda is how to approimate a complicated function f() with a simpler function P () f() P () for eample,
More informationTowards Higher-Order Schemes for Compressible Flow
WDS'6 Proceedings of Contributed Papers, Part I, 5 4, 6. ISBN 8-867-84- MATFYZPRESS Towards Higher-Order Schemes for Compressible Flow K. Findejs Charles Universit, Facult of Mathematics and Phsics, Prague,
More informationσ = F/A. (1.2) σ xy σ yy σ zy , (1.3) σ xz σ yz σ zz The use of the opposite convention should cause no problem because σ ij = σ ji.
Cambridge Universit Press 978-0-521-88121-0 - Metal Forming: Mechanics Metallurg, Third Edition Ecerpt 1 Stress Strain An understing of stress strain is essential for analzing metal forming operations.
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationσ = F/A. (1.2) σ xy σ yy σ zx σ xz σ yz σ, (1.3) The use of the opposite convention should cause no problem because σ ij = σ ji.
Cambridge Universit Press 978-1-107-00452-8 - Metal Forming: Mechanics Metallurg, Fourth Edition Ecerpt 1 Stress Strain An understing of stress strain is essential for the analsis of metal forming operations.
More informationDirect and Iterative Solution of the Generalized Dirichlet-Neumann Map for Elliptic PDEs on Square Domains
Direct and Iterative Solution of the Generalized Dirichlet-Neumann Map for Elliptic PDEs on Square Domains A.G.Sifalakis 1, S.R.Fulton, E.P. Papadopoulou 1 and Y.G.Saridakis 1, 1 Applied Mathematics and
More informationUseful Mathematics. 1. Multivariable Calculus. 1.1 Taylor s Theorem. Monday, 13 May 2013
Useful Mathematics Monday, 13 May 013 Physics 111 In recent years I have observed a reticence among a subpopulation of students to dive into mathematics when the occasion arises in theoretical mechanics
More informationnorges teknisk-naturvitenskapelige universitet Semi-Lagrangian exponential integrators for the incompressible Navier-Stokes equations
norges teknisk-naturvitenskapelige universitet Semi-Lagrangian eponential integrators for the incompressible Navier-Stokes equations b Elena Celledoni, Bawfeh Kingsle Kometa and Olivier Verdier preprint
More information3.6. Implicit Differentiation. Implicitly Defined Functions
3.6 Implicit Differentiation 205 3.6 Implicit Differentiation 5 2 25 2 25 2 0 5 (3, ) Slope 3 FIGURE 3.36 The circle combines the graphs of two functions. The graph of 2 is the lower semicircle an passes
More informationMultiple Integration
Contents 7 Multiple Integration 7. Introduction to Surface Integrals 7. Multiple Integrals over Non-rectangular Regions 7. Volume Integrals 4 7.4 Changing Coordinates 66 Learning outcomes In this Workbook
More informationCONSERVATION OF MASS AND BALANCE OF LINEAR MOMENTUM
CONSERVATION OF MASS AND BALANCE OF LINEAR MOMENTUM Summary of integral theorems Material time derivative Reynolds transport theorem Principle of conservation of mass Principle of balance of linear momentum
More information18.303: Introduction to Green s functions and operator inverses
8.33: Introduction to Green s functions and operator inverses S. G. Johnson October 9, 2 Abstract In analogy with the inverse A of a matri A, we try to construct an analogous inverse  of differential
More informationEnergy Preserving Numerical Integration Methods
Energy Preserving Numerical Integration Methods Department of Mathematics and Statistics, La Trobe University Supported by Australian Research Council Professorial Fellowship Collaborators: Celledoni (Trondheim)
More informationA GLOBAL COLLOCATION METHOD FOR TWO-DIMENSIONAL RECTANGULAR DOMAINS
JOURNAL OF MECHANICS OF MATERIALS AND STRUCTURES Vol., No., A GLOBAL COLLOCATION METHOD FOR TWO-DIMENSIONAL RECTANGULAR DOMAINS CHRISTOPHER G. PROVATIDIS This paper proposes the use of a global collocation
More informationA UNSTRUCTURED NODAL SPECTRAL-ELEMENT METHOD FOR THE NAVIER-STOKES EQUATIONS
A UNSTRUCTURED NODAL SPECTRAL-ELEMENT METHOD FOR THE NAVIER-STOKES EQUATIONS LIZHEN CHEN 1 JIE SHEN 1,2 CHUANJU XU 1,3 Abstract. An unstructured nodal spectral-element method for the Navier-Stokes equations
More informationChapter 1. Condition numbers and local errors in the boundary element method
Chapter 1 Condition numbers and local errors in the boundary element method W. Dijkstra, G. Kakuba and R. M. M. Mattheij Eindhoven University of Technology, Department of Mathematics and Computing Science,
More informationA Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations
A Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations Institute for Mathematics Applied to the Geosciences (IMAGe) National Center for Atmospheric Research (NCAR) Boulder CO
More information