AMS527: Numerical Analysis II
|
|
- Magnus Johns
- 5 years ago
- Views:
Transcription
1 AMS527: Numerical Analysis II A Brief Overview of Finite Element Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 1 / 25
2 Overview Basic concepts Mathematical formulation and discretization Computer implementation of FEM (next lecture) Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 2 / 25
3 Basic Concepts Finite-element method is a spatial discretization method Commonly used in computational structural mechanics for solving linear static problem Sometimes also used in computational fluid dynamics, especially in commercial codes Source: Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 3 / 25
4 Typical Steps of Finite Element Method 1 Write the equation in its weak form, integrated with test functions v(x, y) 2 Subdivide the region into triangles or quadrilaterals 3 Choose N simple trial functions T j (x, y) and look for U = U 1 T U N T N. The 1D hat functions T (x) can change to 2D pyramid functions T (x, y) 4 Produce N equations KU = F from test functions V 1,..., V N (often V j = T j ) 5 Assemble the stiffness matrix K and the load vector F. Solve KU = F. Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 4 / 25
5 Example: Poisson s Equation Strong form 2 u x 2 2 u = f (x, y) in the open set S y 2 u = 0 on the boundary C Weak form multiply strong form by test function v(x, y) and integrate over S ( ) 2 u S x 2 2 u y 2 v dx dy = fv dx dy S Integrate by left side by part, and requiring v to vanish on the boundary: ( u x S v x + u y ) v dx dy = fv dx dy y S In weak form, residual u xx + u yy + f is required to be zero only in a weighted sense, where weight is provided by v Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 5 / 25
6 Galerkin s Method Basic steps of Galerkin s method Choose a finite set of shape functions (trial functions) T i (x, y), i = 1, 2,..., n Admit approximations to u of the form U(x, y) = n j=1 U jt j (x, y) Determine n unknown U i from weak form, using n different test functions To maintain symmetry, require test functions to be the same as trial functions, i.e., v = T i (x, y) Then weak form becomes n T j U j T n i S x x + j=1 j=1 = ft i dx dy S U j T j y T i y dx dy Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 6 / 25
7 Galerkin s Method Cont d Rewrite in matrix form into linear system Ku = f where K is stiffness matrix and f is load vector On left-hand side, we obtain ( Ti K ij = x S T j x + T i y ) T j dx dy y Stiffness matrix K is symmetric in that K ij = K ji. In addition, K is positive definite (later) On right-hand side, we have f i = S ft i dx dy Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 7 / 25
8 Rayleigh-Ritz Method For Poisson s equation, one can obtain same approximation with the minimum principle to minimize [ ( ) 1 u 2 P(u) = + 1 ( ) u 2 fu] dx dy 2 x 2 y S over all functions that satisfy essential boundary condition u = 0 Approximate u in the form U(x, y) = n j=1 U jt j (x, y), then P(u) = 1 2 ut Ku u T f and we obtain the same linear system Ku = f The derivation based on minimum principle is called Rayleigh-Ritz Method From energy P(u), it is obvious that K is symmetric and positive definite (SPD) Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 8 / 25
9 Finite Element Shape Functions Desired properties of shape (trial) functions Function T i must be capable of approximating true u(x, y) and must vanish at boundary Entries K ij and f i must be convenient to compute (ease of differentiation and integration) K should be sparse and well-conditioned Typical shape functions are piecewise polynomials over triangles or quadrilaterals (known as elements) For triangles, piecewise linear (P 1 ), quadratic (P 2 ), and cubic (P 3 ) For quadrilaterals, piecewise bilinear (Q 1 ) and biquadratic (Q 2 ) Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 9 / 25
10 Examples of Element Types Source: Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 10 / 25
11 More on Shape Functions Finite element shape functions are Lagrangian interpolants Each shape function is associated with a vertex It has value 1 at the vertex and vanishes at all other vertices Question: What are the meanings of the coefficients of the shape functions? There are two mathematically equivalent views of shape functions View of each shape function over whole spatial domain View of shape functions within each element Mathematical formulation is based on former, but computational implementation is typically based on latter for convenience Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 11 / 25
12 Examples of Linear Shape Functions Linear (P1) Bilinear (Q1) Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 12 / 25
13 Quadratic Shape Functions for Triangles Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 13 / 25
14 Biquadratic Shape Functions for Quadrilaterals Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 14 / 25
15 Stiffness Matrix To implement finite-element method for Poisson s equation, key is to obtain linear system where K ij = f i = S S ( Ti x ft i dx dy Ku = f T j x + T i y ) T j dx dy y On triangles, each T j is pyramid over jth vertex v j, and ( Ti T j K ij = {e v i e v j e} e x x + T ) i T j dx dy y y f i = ft i dx dy {e v i e v j e} e Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 15 / 25
16 Element Stiffness Matrices Within each triangle, let U = a + bx + cy, compute ( NI N J (k e ) IJ = x x + N ) I N J dx dy, I, J = 1, 2, 3 y y e where N I is shape function w.r.t. I th vertex. Matrix k e = [(k e ) IJ ], I, J = 1, 2, 3 is element stiffness matrix For Laplacian equation, it can be shown that c 2 + c 3 c 3 c 2 k e = c 3 c 1 + c 3 c 1 1, where c i = 2 tan θ c 2 c 1 c 1 + c i 2 Question: Is k e symmetric and positive definite? Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 16 / 25
17 Mass Matrix and Element Mass Matrices Besides stiffness matrix, mass matrix [M ij ] is very common in finite element methods M ij = T i T j dx dy S = T i T j dx dy {e v i e v j e} Its computation involves element mass matrix m e with (m e ) IJ = N I N J dx dy, I, J = 1, 2, 3 e e Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 17 / 25
18 General Procedure for Elemental Computation In general, N I / x and N I / y is obtained from chain rule ] ] ] [ x ξ x η y ξ y η } {{ } J T [ NI x N I y [ x2 x where J = 1 x 3 x 1 y 2 y 1 y 3 y 1 coordinates = [ NI ξ N I η ]. ξ and η are called natural Requires computing Jacobian matrix and derivatives of shape functions Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 18 / 25
19 Example Quadrature Rules Over Triangle Requires quadrature rules over elements Source: Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 19 / 25
20 Example Quadrature Rules Over Triangle Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 20 / 25
21 Assembling Element Matrices For each triangle, k e is 3 3 matrix, each of whose rows (columns) correspond to a vertex Each vertex of triangle has global vertex ID. Example: triangle 1: nodes 3, 1, 2 triangle 2: nodes 3, 4, 2... triangle m: nodes 5, 4, n Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 21 / 25
22 Assembling Element Matrices Map k e for element e to an n n matrix K e, based on local to global mapping of vertex IDs, e.g., (k e ) 22 (k e ) 23 (k e ) 21 (k e ) 32 (k e ) 33 (k e ) 31 (k e ) 12 (k e ) 13 (k e ) 11 K 1 = Then, K = K 1 + K K m Without boundary condition, K is singular! After applying boundary condition, K is SPD In iterative linear solvers involving only matrix-vector multiplications, Kx = K 1 x + K 2 x + + K m x, so K need not be assembled explicitly, leading to matrix-free linear solvers Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 22 / 25
23 Role of Assembler Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 23 / 25
24 Summary of Key Aspects for Implementing FEM Input of finite element methods Vertex coordinates, used for computing Jacobian matrix Element connectivity, used for local to global mapping in assembling matrix Boundary conditions Elemental computation Shape functions, first and second derivatives Numerical quadrature rules and Gauss points Assemble element matrices into stiffness matrix Solution of sparse linear system Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 24 / 25
25 Accuracy and Convergence of Finite Elements When using degree p piecewise polynomials, basic theory of FEM can be stated as: The finite element method converges if p 1, and its error is O(h 2p ), where h is largest edge length. The answers are exact for solutions of degree p (checked by patch test) In addition, finite-element method gives optimal solution U that is closest to exact solution u in that it minimizes the energy E = (U u) T K(U u) within all feature solutions in function space of trial functions Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 25 / 25
f xx + f yy = F (x, y)
Application of the 2D finite element method to Laplace (Poisson) equation; f xx + f yy = F (x, y) M. R. Hadizadeh Computer Club, Department of Physics and Astronomy, Ohio University 4 Nov. 2013 Domain
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 20: Sparse Linear Systems; Direct Methods vs. Iterative Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 26
More informationWhen implementing FEM for solving two-dimensional partial differential equations, integrals of the form
Quadrature Formulas in Two Dimensions Math 57 - Finite Element Method Section, Spring Shaozhong Deng, PhD (shaodeng@unccedu Dept of Mathematics and Statistics, UNC at Charlotte When implementing FEM for
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationAdaptive numerical methods
METRO MEtallurgical TRaining On-line Adaptive numerical methods Arkadiusz Nagórka CzUT Education and Culture Introduction Common steps of finite element computations consists of preprocessing - definition
More informationFEM Convergence Requirements
19 FEM Convergence Requirements IFEM Ch 19 Slide 1 Convergence Requirements for Finite Element Discretization Convergence: discrete (FEM) solution approaches the analytical (math model) solution in some
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 36
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 36 In last class, we have derived element equations for two d elasticity problems
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 5: Sparse Linear Systems and Factorization Methods Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis I 1 / 18 Sparse
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 3, 2017, Lesson 1
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 3, 2017, Lesson 1 1 Politecnico di Milano, February 3, 2017, Lesson 1 2 Outline
More informationIMA Preprint Series #
WEB-SPLINE FINITE ELEMENTS By Yanping Cao, Olga Terlyga, Jon Van Laarhoven Jianbao Wu, Guangri Xue, Ping Zhang and Thomas Grandine (mentor) IMA Preprint Series # 2133-2 ( September 2006 ) INSTITUTE FOR
More informationIndex. C m (Ω), 141 L 2 (Ω) space, 143 p-th order, 17
Bibliography [1] J. Adams, P. Swarztrauber, and R. Sweet. Fishpack: Efficient Fortran subprograms for the solution of separable elliptic partial differential equations. http://www.netlib.org/fishpack/.
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 24
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 24 So in today s class, we will look at quadrilateral elements; and we will
More informationChapter 13. Boundary Value Problems for Partial Differential Equations* Linz 2002/ page
Chapter 13 Boundary Value Problems for Partial Differential Equations* E lliptic equations constitute the third category of partial differential equations. As a prototype, we take the Poisson equation
More informationIsogeometric Collocation Method
Chair for Computational Analysis of Technical Systems Faculty of Mechanical Engineering, RWTH Aachen University Isogeometric Collocation Method Seminararbeit By Marko Blatzheim Supervisors: Dr. Stefanie
More informationToday. Motivation. Motivation. Image gradient. Image gradient. Computational Photography
Computational Photography Matthias Zwicker University of Bern Fall 009 Today Gradient domain image manipulation Introduction Gradient cut & paste Tone mapping Color-to-gray conversion Motivation Cut &
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationComputational Fluid Dynamics - Incompressible Flows
Computational Fluid Dynamics - Incompressible Flows March 25, 2008 Incompressible Flows Basis Functions Discrete Equations CFD - Incompressible Flows CFD is a Huge field Numerical Techniques for solving
More informationNumerical Methods for PDEs : Video 11: 1D FiniteFebruary Difference 15, Mappings Theory / 15
22.520 Numerical Methods for PDEs : Video 11: 1D Finite Difference Mappings Theory and Matlab February 15, 2015 22.520 Numerical Methods for PDEs : Video 11: 1D FiniteFebruary Difference 15, Mappings 2015
More informationSet No. 1 IV B.Tech. I Semester Regular Examinations, November 2010 FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 Answer any FIVE Questions All Questions carry equal marks
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationComputational Physics PHYS 420
Computational Physics PHYS 420 Dr Richard H. Cyburt Assistant Professor of Physics My office: 402c in the Science Building My phone: (304) 384-6006 My email: rcyburt@concord.edu My webpage: www.concord.edu/rcyburt
More information1 Exercise: Heat equation in 2-D with FE
1 Exercise: Heat equation in 2-D with FE Reading Hughes (2000, sec. 2.3-2.6 Dabrowski et al. (2008, sec. 1-3, 4.1.1, 4.1.3, 4.2.1 This FE exercise and most of the following ones are based on the MILAMIN
More informationNumerical integration of polynomials and discontinuous functions on irregular convex polygons and polyhedrons
Comput Mech (2011) 47:535 554 DOI 10.1007/s00466-010-0562-5 ORIGINAL PAPER Numerical integration of polynomials and discontinuous functions on irregular convex polygons and polyhedrons S. E. Mousavi N.
More informationFoundations of Analytical and Numerical Field Computation
Foundations of Analytical and Numerical Field Computation Stephan Russenschuck, CERN-AT-MEL Stephan Russenschuck CERN, TE-MCS, 1211 Geneva, Switzerland 1 Permanent Magnet Circuits 2 Rogowski profiles Pole
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 informationAdaptivity For Meshfree Point Collocation Methods In Linear Elastic Solid Mechanics
University of South Carolina Scholar Commons Theses and Dissertations 2015 Adaptivity For Meshfree Point Collocation Methods In Linear Elastic Solid Mechanics Joshua Wayne Derrick University of South Carolina
More informationBarycentric Finite Element Methods
University of California, Davis Barycentric Finite Element Methods N. Sukumar UC Davis Workshop on Generalized Barycentric Coordinates, Columbia University July 26, 2012 Collaborators and Acknowledgements
More informationBarycentric Finite Element Methods
University of California, Davis Barycentric Finite Element Methods N. Sukumar University of California at Davis SIAM Conference on Geometric Design and Computing November 8, 2007 Collaborators and Acknowledgements
More informationDocumentation for Numerical Derivative on Discontinuous Galerkin Space
Documentation for Numerical Derivative on Discontinuous Galerkin Space Stefan Schnake 204 Introduction This documentation gives a guide to the syntax and usage of the functions in this package as simply
More informationNumerical integration of polynomials and discontinuous functions on irregular convex polygons and polyhedrons
Noname manuscript No. (will be inserted by the editor) Numerical integration of polynomials and discontinuous functions on irregular convex polygons and polyhedrons S. E. Mousavi N. Sukumar Received: date
More information3D NURBS-ENHANCED FINITE ELEMENT METHOD
7th Workshop on Numerical Methods in Applied Science and Engineering (NMASE 8) Vall de Núria, 9 a 11 de enero de 28 c LaCàN, www.lacan-upc.es 3D NURBS-ENHANCED FINITE ELEMENT METHOD R. Sevilla, S. Fernández-Méndez
More informationMicroprocessor Thermal Analysis using the Finite Element Method
Microprocessor Thermal Analysis using the Finite Element Method Bhavya Daya Massachusetts Institute of Technology Abstract The microelectronics industry is pursuing many options to sustain the performance
More informationA Hybrid Geometric+Algebraic Multigrid Method with Semi-Iterative Smoothers
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 013; 00:1 18 Published online in Wiley InterScience www.interscience.wiley.com). A Hybrid Geometric+Algebraic Multigrid Method with
More informationCommunications in Applied Mathematics and Computational Science
Communications in Applied Mathematics and Computational Science Volume 1 No. 1 2006 A COMPARISON OF THE EXTENDED FINITE ELEMENT METHOD WITH THE IMMERSED INTERFACE METHOD FOR ELLIPTIC EQUATIONS WITH DISCONTINUOUS
More informationThe numerical simulation of complex PDE problems. A numerical simulation project The finite element method for solving a boundary-value problem in R 2
Universidad de Chile The numerical simulation of complex PDE problems Facultad de Ciencias Físicas y Matemáticas P. Frey, M. De Buhan Year 2008 MA691 & CC60X A numerical simulation project The finite element
More informationElement Quality Metrics for Higher-Order Bernstein Bézier Elements
Element Quality Metrics for Higher-Order Bernstein Bézier Elements Luke Engvall and John A. Evans Abstract In this note, we review the interpolation theory for curvilinear finite elements originally derived
More informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More informationA new 8-node quadrilateral spline finite element
Journal of Computational and Applied Mathematics 195 (2006) 54 65 www.elsevier.com/locate/cam A new 8-node quadrilateral spline finite element Chong-Jun Li, Ren-Hong Wang Institute of Mathematical Sciences,
More informationModule 1: Introduction to Finite Element Analysis. Lecture 4: Steps in Finite Element Analysis
25 Module 1: Introduction to Finite Element Analysis Lecture 4: Steps in Finite Element Analysis 1.4.1 Loading Conditions There are multiple loading conditions which may be applied to a system. The load
More informationAdaptive boundary element methods in industrial applications
Adaptive boundary element methods in industrial applications Günther Of, Olaf Steinbach, Wolfgang Wendland Adaptive Fast Boundary Element Methods in Industrial Applications, Hirschegg, September 30th,
More informationWhat is Multigrid? They have been extended to solve a wide variety of other problems, linear and nonlinear.
AMSC 600/CMSC 760 Fall 2007 Solution of Sparse Linear Systems Multigrid, Part 1 Dianne P. O Leary c 2006, 2007 What is Multigrid? Originally, multigrid algorithms were proposed as an iterative method to
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 1: Course Overview; Matrix Multiplication Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis I 1 / 21 Outline 1 Course
More informationFinite Element Methods
Chapter 5 Finite Element Methods 5.1 Finite Element Spaces Remark 5.1 Mesh cells, faces, edges, vertices. A mesh cell is a compact polyhedron in R d, d {2,3}, whose interior is not empty. The boundary
More informationNumerical Analysis I - Final Exam Matrikelnummer:
Dr. Behrens Center for Mathematical Sciences Technische Universität München Winter Term 2005/2006 Name: Numerical Analysis I - Final Exam Matrikelnummer: I agree to the publication of the results of this
More informationParallel Implementations of Gaussian Elimination
s of Western Michigan University vasilije.perovic@wmich.edu January 27, 2012 CS 6260: in Parallel Linear systems of equations General form of a linear system of equations is given by a 11 x 1 + + a 1n
More informationNumerical Methods in Physics Lecture 2 Interpolation
Numerical Methods in Physics Pat Scott Department of Physics, Imperial College November 8, 2016 Slides available from http://astro.ic.ac.uk/pscott/ course-webpage-numerical-methods-201617 Outline The problem
More informationFinite element solution of multi-scale transport problems using the least squares based bubble function enrichment
Finite element solution of multi-scale transport problems using the least squares based bubble function enrichment A. Yazdani a, V. Nassehi b1 a Cranfield University, School of Applied Sciences, Cranfield,
More informationIntroduction to Finite Element Analysis using ANSYS
Introduction to Finite Element Analysis using ANSYS Sasi Kumar Tippabhotla PhD Candidate Xtreme Photovoltaics (XPV) Lab EPD, SUTD Disclaimer: The material and simulations (using Ansys student version)
More informationAn introduction to mesh generation Part IV : elliptic meshing
Elliptic An introduction to mesh generation Part IV : elliptic meshing Department of Civil Engineering, Université catholique de Louvain, Belgium Elliptic Curvilinear Meshes Basic concept A curvilinear
More informationGalerkin Projections Between Finite Element Spaces
Galerkin Projections Between Finite Element Spaces Ross A. Thompson Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements
More informationlecture 10: B-Splines
9 lecture : -Splines -Splines: a basis for splines Throughout our discussion of standard polynomial interpolation, we viewed P n as a linear space of dimension n +, and then expressed the unique interpolating
More informationTHE MORTAR FINITE ELEMENT METHOD IN 2D: IMPLEMENTATION IN MATLAB
THE MORTAR FINITE ELEMENT METHOD IN D: IMPLEMENTATION IN MATLAB J. Daněk, H. Kutáková Department of Mathematics, University of West Bohemia, Pilsen MECAS ESI s.r.o., Pilsen Abstract The paper is focused
More informationMeshless Modeling, Animating, and Simulating Point-Based Geometry
Meshless Modeling, Animating, and Simulating Point-Based Geometry Xiaohu Guo SUNY @ Stony Brook Email: xguo@cs.sunysb.edu http://www.cs.sunysb.edu/~xguo Graphics Primitives - Points The emergence of points
More informationA comparison of accuracy and computational efficiency between the Finite Element Method and the Isogeometric Analysis for twodimensional
A comparison of accuracy and computational efficiency between the Finite Element Method and the Isogeometric Analysis for twodimensional Poisson Problems Per Ståle Larsen Master of Science in Physics and
More informationFinite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More informationIN FINITE element simulations usually two computation
Proceedings of the International Multiconference on Computer Science and Information Technology pp. 337 342 Higher order FEM numerical integration on GPUs with OpenCL 1 Przemysław Płaszewski, 12 Krzysztof
More informationAdaptive Isogeometric Analysis by Local h-refinement with T-splines
Adaptive Isogeometric Analysis by Local h-refinement with T-splines Michael Dörfel 1, Bert Jüttler 2, Bernd Simeon 1 1 TU Munich, Germany 2 JKU Linz, Austria SIMAI, Minisymposium M13 Outline Preliminaries:
More informationFebruary 23 Math 2335 sec 51 Spring 2016
February 23 Math 2335 sec 51 Spring 2016 Section 4.1: Polynomial Interpolation Interpolation is the process of finding a curve or evaluating a function whose curve passes through a known set of points.
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationTopology Optimization of Two Linear Elastic Bodies in Unilateral Contact
2 nd International Conference on Engineering Optimization September 6-9, 2010, Lisbon, Portugal Topology Optimization of Two Linear Elastic Bodies in Unilateral Contact Niclas Strömberg Department of Mechanical
More informationf ( 1,η ) 1,2 2,2 2,1
Chapter 5 Mesh Generation and Assembly 5.1 Introduction There are several reasons for the popularity of nite element methods. Large code segments can be implemented for a wide class of problems. The software
More informationSpace Filling Curves and Hierarchical Basis. Klaus Speer
Space Filling Curves and Hierarchical Basis Klaus Speer Abstract Real world phenomena can be best described using differential equations. After linearisation we have to deal with huge linear systems of
More informationFinite element algorithm with adaptive quadtree-octree mesh refinement
ANZIAM J. 46 (E) ppc15 C28, 2005 C15 Finite element algorithm with adaptive quadtree-octree mesh refinement G. P. Nikishkov (Received 18 October 2004; revised 24 January 2005) Abstract Certain difficulties
More informationGeometric Modeling Assignment 3: Discrete Differential Quantities
Geometric Modeling Assignment : Discrete Differential Quantities Acknowledgements: Julian Panetta, Olga Diamanti Assignment (Optional) Topic: Discrete Differential Quantities with libigl Vertex Normals,
More informationARE BILINEAR QUADRILATERALS BETTER THAN LINEAR TRIANGLES?
ARE BILINEAR QUADRILATERALS BETTER THAN LINEAR TRIANGLES? E. F. D AZEVEDO y Abstract. This paper compares the theoretical effectiveness of bilinear approximation over quadrilaterals with linear approximation
More informationMethod of Finite Elements I
Institute of Structural Engineering Page 1 Held by Prof. Dr. E. Chatzi, Dr. P. Steffen Assistants: Adrian Egger (HIL E 13.3), Harry Mylonas (HIL H33.1), Konstantinos Tatsis (HIL H33.1) Lectures homepage:
More informationMotivation Patch tests Numerical examples Conclusions
Weakening the tight coupling between geometry and simulation in isogeometric analysis: from sub- and super- geometric analysis to Geometry Independent Field approximation (GIFT) Elena Atroshchenko, Gang
More informationAS 5850 Finite Element Analysis
AS 5850 Finite Element Analysis INTRODUCTION Prof. IIT Madras Goal of engineering computations Perform analysis and design of physical systems and processes subjected to imposed conditions (or loads) and
More informationIterative Algorithms I: Elementary Iterative Methods and the Conjugate Gradient Algorithms
Iterative Algorithms I: Elementary Iterative Methods and the Conjugate Gradient Algorithms By:- Nitin Kamra Indian Institute of Technology, Delhi Advisor:- Prof. Ulrich Reude 1. Introduction to Linear
More informationA meshfree weak-strong form method
A meshfree weak-strong form method G. R. & Y. T. GU' 'centre for Advanced Computations in Engineering Science (ACES) Dept. of Mechanical Engineering, National University of Singapore 2~~~ Fellow, Singapore-MIT
More information3D Modeling Parametric Curves & Surfaces. Shandong University Spring 2013
3D Modeling Parametric Curves & Surfaces Shandong University Spring 2013 3D Object Representations Raw data Point cloud Range image Polygon soup Surfaces Mesh Subdivision Parametric Implicit Solids Voxels
More informationAn Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (2000/2001)
An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (000/001) Summary The objectives of this project were as follows: 1) Investigate iterative
More informationImage Registration Lecture 4: First Examples
Image Registration Lecture 4: First Examples Prof. Charlene Tsai Outline Example Intensity-based registration SSD error function Image mapping Function minimization: Gradient descent Derivative calculation
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationFinite Element Modeling Techniques (2) دانشگاه صنعتي اصفهان- دانشكده مكانيك
Finite Element Modeling Techniques (2) 1 Where Finer Meshes Should be Used GEOMETRY MODELLING 2 GEOMETRY MODELLING Reduction of a complex geometry to a manageable one. 3D? 2D? 1D? Combination? Bulky solids
More information(x,y ) 3 (x,y ) (x,y ) φ 1
Chapter 4 Finite Element Approximation 4. Introduction Our goal in this chapter is the development of piecewise-polynomial approximations U of a two- or three-dimensional function u. For this purpose,
More informationInvestigation of a Robust Method for Connecting Dissimilar 3D Finite Element Models. David M. Trujillo 1. December 2005
Investigation of a Robust Method for Connecting Dissimilar 3D Finite Element Models by David M. Trujillo 1 December 2005 1 Consultant, TRUCOMP, Fountain Valley, California trucomp@earthlink.net Abstract
More informationTransfinite Interpolation Based Analysis
Transfinite Interpolation Based Analysis Nathan Collier 1 V.M. Calo 1 Javier Principe 2 1 King Abdullah University of Science and Technology 2 International Center for Numerical Methods in Engineering
More informationLinear Finite Element Methods
Chapter 3 Linear Finite Element Methods The finite element methods provide spaces V n of functions that are piecewise smooth and simple, and locally supported basis function of these spaces to achieve
More informationShape Functions, Derivatives, and Integration
Shape Functions, Derivatives, and Integration CHAPTER 6 6. Introduction In the previous chapter we found that the quasi-harmonic equation created a weak form that contained derivatives of the dependent
More informationCubic B-spline Solution of Two-point Boundary Value Problem Using HSKSOR Iteration
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 11 (2017), pp. 7921-7934 Research India Publications http://www.ripublication.com Cubic B-spline Solution of Two-point Boundary
More informationCITY AND GUILDS 9210 UNIT 135 MECHANICS OF SOLIDS Level 6 TUTORIAL 15 - FINITE ELEMENT ANALYSIS - PART 1
Outcome 1 The learner can: CITY AND GUILDS 9210 UNIT 135 MECHANICS OF SOLIDS Level 6 TUTORIAL 15 - FINITE ELEMENT ANALYSIS - PART 1 Calculate stresses, strain and deflections in a range of components under
More informationFinite Element Implementation
Chapter 8 Finite Element Implementation 8.1 Elements Elements andconditions are the main extension points of Kratos. New formulations can be introduced into Kratos by implementing a new Element and its
More informationPolynomial and Rational Functions. Copyright Cengage Learning. All rights reserved.
2 Polynomial and Rational Functions Copyright Cengage Learning. All rights reserved. 2.1 Quadratic Functions Copyright Cengage Learning. All rights reserved. What You Should Learn Analyze graphs of quadratic
More informationInformation Coding / Computer Graphics, ISY, LiTH. Splines
28(69) Splines Originally a drafting tool to create a smooth curve In computer graphics: a curve built from sections, each described by a 2nd or 3rd degree polynomial. Very common in non-real-time graphics,
More informationConcept of Curve Fitting Difference with Interpolation
Curve Fitting Content Concept of Curve Fitting Difference with Interpolation Estimation of Linear Parameters by Least Squares Curve Fitting by Polynomial Least Squares Estimation of Non-linear Parameters
More informationComputer Graphics. Lecture 14 Bump-mapping, Global Illumination (1)
Computer Graphics Lecture 14 Bump-mapping, Global Illumination (1) Today - Bump mapping - Displacement mapping - Global Illumination Radiosity Bump Mapping - A method to increase the realism of 3D objects
More informationPuffin User Manual. March 1, Johan Hoffman and Anders Logg.
Puffin User Manual March 1, 2006 Johan Hoffman and Anders Logg www.fenics.org Visit http://www.fenics.org/ for the latest version of this manual. Send comments and suggestions to puffin-dev@fenics.org.
More informationA New Approach for Solving Partial Differential Equations Based on B-splines and Transfinite Interpolation
A New Approach for Solving Partial Differential Equations Based on B-splines and Transfinite Interpolation Yuanjie Liu Key Laboratory of Mathematics Mechanization AMSS Chinese Academy of Sciences, 100190,
More informationEXTENSION. a 1 b 1 c 1 d 1. Rows l a 2 b 2 c 2 d 2. a 3 x b 3 y c 3 z d 3. This system can be written in an abbreviated form as
EXTENSION Using Matrix Row Operations to Solve Systems The elimination method used to solve systems introduced in the previous section can be streamlined into a systematic method by using matrices (singular:
More information3D Modeling Parametric Curves & Surfaces
3D Modeling Parametric Curves & Surfaces Shandong University Spring 2012 3D Object Representations Raw data Point cloud Range image Polygon soup Solids Voxels BSP tree CSG Sweep Surfaces Mesh Subdivision
More informationMaple as an Instructional Tool
Maple as an Instructional Tool Terence A. Weigel 1 Abstract Computer algebra software such as Maple is an important component of the engineer s toolkit, much as are Matlib, MathCAD and Excel. However,
More informationIll-Posed Problems with A Priori Information
INVERSE AND ILL-POSED PROBLEMS SERIES Ill-Posed Problems with A Priori Information V.V.Vasin andalageev HIV SPIII Utrecht, The Netherlands, 1995 CONTENTS Introduction 1 CHAPTER 1. UNSTABLE PROBLEMS 1 Base
More informationComputer Graphics Curves and Surfaces. Matthias Teschner
Computer Graphics Curves and Surfaces Matthias Teschner Outline Introduction Polynomial curves Bézier curves Matrix notation Curve subdivision Differential curve properties Piecewise polynomial curves
More informationSpectral(Finite) Volume Method for Conservation Laws on Unstructured Grids
Journal of Computational Physics 179, 665 697 (2002) doi:10.1006/jcph.2002.7082 Spectral(Finite) Volume Method for Conservation Laws on Unstructured Grids II. Extension to Two-Dimensional Scalar Equation
More informationFor each question, indicate whether the statement is true or false by circling T or F, respectively.
True/False For each question, indicate whether the statement is true or false by circling T or F, respectively. 1. (T/F) Rasterization occurs before vertex transformation in the graphics pipeline. 2. (T/F)
More informationReconstructing High-Order Surfaces for Meshing
Engineering with Computers manuscript No. (will be inserted by the editor) Reconstructing High-Order Surfaces for Meshing Xiangmin Jiao, Duo Wang Department of Applied Mathematics, Stony Brook University,
More informationLinear Interpolating Splines
Jim Lambers MAT 772 Fall Semester 2010-11 Lecture 17 Notes Tese notes correspond to Sections 112, 11, and 114 in te text Linear Interpolating Splines We ave seen tat ig-degree polynomial interpolation
More informationThis is the main idea of the evolution Galerkin (EG) methods, which evolve the initial data using the bicharacteristic cone and then project them onto
Finite volume evolution Galerkin methods for multidimensional hyperbolic problems M. Lukacova { Medvid'ova 1 3, K. W. Morton 2,G.Warnecke 1 1 Institut fur Analysis und Numerik, Otto-von-Guericke-Universitat
More information