Transfinite Interpolation Based Analysis

Size: px
Start display at page:

Download "Transfinite Interpolation Based Analysis"

Transcription

1 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 7 June 2010 Collier (KAUST) Transfinite Element Analysis 7 June / 26

2 Motivation - Collier (KAUST) Transfinite Element Analysis 7 June / 26

3 Overview 1 Transfinite Element Analysis Basis Function Construction Sample Basis Functions 2 Application 1: Terrain Modeling Automatic Refinement Insert Features 3 Application 2: Finite Elements with Simplified Meshing Current Challenges 4 Conclusions Collier (KAUST) Transfinite Element Analysis 7 June / 26

4 Interpolate Features 1 u1 d1 u d2 u2 2 1 f1(u1) F (u) f2(u2) 2 d3 3 u3 f3(u3) 3 i 1 d 1 F (u) = f 1 (u 1 ) 1 d d f 2 (u 2 ) 1 d 3 d d f 3 (u 3 ) 1 d 3 d d d 3 n f W i n f = f i (u i ) nf = f i (u i )Ŵ i (d 1, d 2, d 3 ) j W j i 1 d 2 1 d 3 Collier (KAUST) Transfinite Element Analysis 7 June / 26

5 Basis Function Construction n f F (u) = f i (u i )Ŵi(d) Each footprint can then be parameterized by its own one-dimensional basis, B i,k n b f i (u) = B i,k (u)c i,k which substituted into the surface approximation becomes n f n b F (u) = B i,k (u)c i,k Ŵ i (d) i k n f n b i k ( ) = B i,k (u)ŵ i (d) C i,k = N j C j i k }{{} j N i,k Collier (KAUST) Transfinite Element Analysis 7 June / 26

6 Sample Basis Functions Linear Footprints, Convex Polygons Collier (KAUST) Transfinite Element Analysis 7 June / 26

7 Sample Basis Functions Gradients N N x N y Node 0 Node 1 Node 2 Collier (KAUST) Transfinite Element Analysis 7 June / 26

8 Derivative Control: Ribbons Not only can we construct functions which interpolate features, but we can also control the derivative as the feature is crossed. 1 u1 d1 u d2 u2 2 1 f1(u1) F (u) f2(u2) 2 d3 3 u3 f3(u3) 3 Collier (KAUST) Transfinite Element Analysis 7 June / 26

9 Bi-cubic: 28 vs 16 Collier (KAUST) Transfinite Element Analysis 7 June / 26

10 Non-standard Basis Functions Collier (KAUST) Transfinite Element Analysis 7 June / 26

11 Overview 1 Transfinite Element Analysis Basis Function Construction Sample Basis Functions 2 Application 1: Terrain Modeling Automatic Refinement Insert Features 3 Application 2: Finite Elements with Simplified Meshing Current Challenges 4 Conclusions Collier (KAUST) Transfinite Element Analysis 7 June / 26

12 Projection onto Function Space Currently using L 2 projections. Find U A such that argmin U A Z H(U A ) 2 where Z is a function to approximate, and H is the interpolation which is controlled by parameters U A. Collier (KAUST) Transfinite Element Analysis 7 June / 26

13 Footprint Refinement Each footprint can be refined locally due to the surface interpolating each footprint. Collier (KAUST) Transfinite Element Analysis 7 June / 26

14 Footprint Refinement Tensor Product (Old way) Transfinite (New way) Collier (KAUST) Transfinite Element Analysis 7 June / 26

15 Footprint Refinement Other degrees of freedom can then be obtained via a surface L 2 projection G 1 Cubic G 2 Quintic Collier (KAUST) Transfinite Element Analysis 7 June / 26

16 Footprint Insertion Combine image processing with projective interpolation Collier (KAUST) Transfinite Element Analysis 7 June / 26

17 Footprint Insertion Combine image processing with projective interpolation Collier (KAUST) Transfinite Element Analysis 7 June / 26

18 Footprint Insertion C 0 Continuity C 1 Continuity Collier (KAUST) Transfinite Element Analysis 7 June / 26

19 Footprint Insertion Collier (KAUST) Transfinite Element Analysis 7 June / 26

20 Overview 1 Transfinite Element Analysis Basis Function Construction Sample Basis Functions 2 Application 1: Terrain Modeling Automatic Refinement Insert Features 3 Application 2: Finite Elements with Simplified Meshing Current Challenges 4 Conclusions Collier (KAUST) Transfinite Element Analysis 7 June / 26

21 Motivation - Collier (KAUST) Transfinite Element Analysis 7 June / 26

22 Laplace Equation Find θ(x) such that with boundary conditions θ = 0 x Ω θ x=0 = 0, n θ x=1 = 1, n θ y=0,1 = 0 Collier (KAUST) Transfinite Element Analysis 7 June / 26

23 Integration Check the relative change in the stiffness matrix diagonal terms as refining the integration mesh. Gaussian Quadrature 4x4 rule used with uniform refinements. Collier (KAUST) Transfinite Element Analysis 7 June / 26

24 Discontinuity B 5 (x, y) d dx B 5(x, y) d dy B 5(x, y) Collier (KAUST) Transfinite Element Analysis 7 June / 26

25 Discontinuity - Removed by Assembly B 5 (x, y) d dx B 5(x, y) d dy B 5(x, y) Pre-Assembly Post-Assembly Collier (KAUST) Transfinite Element Analysis 7 June / 26

26 Conclusion Transfinite Interpolation is a powerful, flexible way to define a basis Grants the freedom to place degrees of freedom at will without propagating them to unintentional regions of the domain Ties in directly to CAD programs whose geometries are watertight and requires no meshing in the traditional sense Current work is focused in accurate integration techniques Collier (KAUST) Transfinite Element Analysis 7 June / 26

Motivation Patch tests Numerical examples Conclusions

Motivation 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 information

Computational Fluid Dynamics - Incompressible Flows

Computational 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 information

Adaptive numerical methods

Adaptive 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 information

3D Modeling Parametric Curves & Surfaces. Shandong University Spring 2013

3D 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 information

Finite 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 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 information

AMS527: Numerical Analysis II

AMS527: Numerical Analysis II 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 Overview Basic concepts Mathematical

More information

Dual Interpolants for Finite Element Methods

Dual Interpolants for Finite Element Methods Dual Interpolants for Finite Element Methods Andrew Gillette joint work with Chandrajit Bajaj and Alexander Rand Department of Mathematics Institute of Computational Engineering and Sciences University

More information

Contents. I The Basic Framework for Stationary Problems 1

Contents. 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 information

Barycentric Finite Element Methods

Barycentric 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 information

3D Modeling Parametric Curves & Surfaces

3D 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 information

HYBRID DIRECT AND ITERATIVE SOLVER FOR H ADAPTIVE MESHES WITH POINT SINGULARITIES

HYBRID DIRECT AND ITERATIVE SOLVER FOR H ADAPTIVE MESHES WITH POINT SINGULARITIES HYBRID DIRECT AND ITERATIVE SOLVER FOR H ADAPTIVE MESHES WITH POINT SINGULARITIES Maciej Paszyński Department of Computer Science, AGH University of Science and Technology, Kraków, Poland email: paszynsk@agh.edu.pl

More information

Parameterization. Michael S. Floater. November 10, 2011

Parameterization. Michael S. Floater. November 10, 2011 Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point

More information

Finite Element Analysis Using Creo Simulate 4.0

Finite Element Analysis Using Creo Simulate 4.0 Introduction to Finite Element Analysis Using Creo Simulate 4.0 Randy H. Shih SDC PUBLICATIONS Better Textbooks. Lower Prices. www.sdcpublications.com Powered by TCPDF (www.tcpdf.org) Visit the following

More information

2D Object Definition (1/3)

2D Object Definition (1/3) 2D Object Definition (1/3) Lines and Polylines Lines drawn between ordered points to create more complex forms called polylines Same first and last point make closed polyline or polygon Can intersect itself

More information

Shape modeling Modeling technique Shape representation! 3D Graphics Modeling Techniques

Shape modeling Modeling technique Shape representation! 3D Graphics   Modeling Techniques D Graphics http://chamilo2.grenet.fr/inp/courses/ensimag4mmgd6/ Shape Modeling technique Shape representation! Part : Basic techniques. Projective rendering pipeline 2. Procedural Modeling techniques Shape

More information

CS 372: Computational Geometry Lecture 3 Line Segment Intersection

CS 372: Computational Geometry Lecture 3 Line Segment Intersection CS 372: Computational Geometry Lecture 3 Line Segment Intersection Antoine Vigneron King Abdullah University of Science and Technology September 9, 2012 Antoine Vigneron (KAUST) CS 372 Lecture 3 September

More information

Lecture 3.2 Methods for Structured Mesh Generation

Lecture 3.2 Methods for Structured Mesh Generation Lecture 3.2 Methods for Structured Mesh Generation 1 There are several methods to develop the structured meshes: Algebraic methods, Interpolation methods, and methods based on solving partial differential

More information

Set 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 information

Implicit Surfaces & Solid Representations COS 426

Implicit Surfaces & Solid Representations COS 426 Implicit Surfaces & Solid Representations COS 426 3D Object Representations Desirable properties of an object representation Easy to acquire Accurate Concise Intuitive editing Efficient editing Efficient

More information

The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a

The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a coordinate system and then the measuring of the point with

More information

2D Model For Steady State Temperature Distribution

2D Model For Steady State Temperature Distribution 2D Model For Steady State Temperature Distribution Finite Element Method Vinh University of Massachusetts Dartmouth December 14, 21 Introduction Advisor Dr. Nima Rahbar: Civil Engineering Project Objective

More information

Parameterization of triangular meshes

Parameterization of triangular meshes Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to

More information

Know it. Control points. B Spline surfaces. Implicit surfaces

Know it. Control points. B Spline surfaces. Implicit surfaces Know it 15 B Spline Cur 14 13 12 11 Parametric curves Catmull clark subdivision Parametric surfaces Interpolating curves 10 9 8 7 6 5 4 3 2 Control points B Spline surfaces Implicit surfaces Bezier surfaces

More information

Numerical integration of polynomials and discontinuous functions on irregular convex polygons and polyhedrons

Numerical 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 information

Numerical integration of polynomials and discontinuous functions on irregular convex polygons and polyhedrons

Numerical 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 information

Chapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling

Chapter 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 information

Def De orma f tion orma Disney/Pixar

Def De orma f tion orma Disney/Pixar Deformation Disney/Pixar Deformation 2 Motivation Easy modeling generate new shapes by deforming existing ones 3 Motivation Easy modeling generate new shapes by deforming existing ones 4 Motivation Character

More information

PetIGA. A Framework for High Performance Isogeometric Analysis. Santa Fe, Argentina. Knoxville, United States. Thuwal, Saudi Arabia

PetIGA. A Framework for High Performance Isogeometric Analysis. Santa Fe, Argentina. Knoxville, United States. Thuwal, Saudi Arabia PetIGA A Framework for High Performance Isogeometric Analysis Lisandro Dalcin 1,3, Nathaniel Collier 2, Adriano Côrtes 3, Victor M. Calo 3 1 Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)

More information

B-Spline Polynomials. B-Spline Polynomials. Uniform Cubic B-Spline Curves CS 460. Computer Graphics

B-Spline Polynomials. B-Spline Polynomials. Uniform Cubic B-Spline Curves CS 460. Computer Graphics CS 460 B-Spline Polynomials Computer Graphics Professor Richard Eckert March 24, 2004 B-Spline Polynomials Want local control Smoother curves B-spline curves: Segmented approximating curve 4 control points

More information

Adaptive Isogeometric Analysis by Local h-refinement with T-splines

Adaptive 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 information

Modeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures

Modeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures University of California, Davis Modeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures N. Sukumar UC Davis Rutgers University,

More information

3D Modeling techniques

3D Modeling techniques 3D Modeling techniques 0. Reconstruction From real data (not covered) 1. Procedural modeling Automatic modeling of a self-similar objects or scenes 2. Interactive modeling Provide tools to computer artists

More information

Generative Part Structural Analysis Fundamentals

Generative Part Structural Analysis Fundamentals CATIA V5 Training Foils Generative Part Structural Analysis Fundamentals Version 5 Release 19 September 2008 EDU_CAT_EN_GPF_FI_V5R19 About this course Objectives of the course Upon completion of this course

More information

Subdivision Surfaces

Subdivision Surfaces Subdivision Surfaces 1 Geometric Modeling Sometimes need more than polygon meshes Smooth surfaces Traditional geometric modeling used NURBS Non uniform rational B-Spline Demo 2 Problems with NURBS A single

More information

Conforming Vector Interpolation Functions for Polyhedral Meshes

Conforming Vector Interpolation Functions for Polyhedral Meshes Conforming Vector Interpolation Functions for Polyhedral Meshes Andrew Gillette joint work with Chandrajit Bajaj and Alexander Rand Department of Mathematics Institute of Computational Engineering and

More information

Digital Geometry Processing Parameterization I

Digital Geometry Processing Parameterization I Problem Definition Given a surface (mesh) S in R 3 and a domain find a bective F: S Typical Domains Cutting to a Disk disk = genus zero + boundary sphere = closed genus zero Creates artificial boundary

More information

PROPERTIES OF NATURAL ELEMENT COORDINATES ON ANY POLYHEDRON

PROPERTIES OF NATURAL ELEMENT COORDINATES ON ANY POLYHEDRON PROPRTIS OF NATURAL LMNT COORDINATS ON ANY POLYHDRON P. Milbradt and T. Fröbel Institute of Computer Science in Civil ngineering, Univercity of Hanover, 3067, Hanover, Germany; PH (+49) 5-76-5757; FAX

More information

3D NURBS-ENHANCED FINITE ELEMENT METHOD

3D 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 information

Shape Modeling. Differential Geometry Primer Smooth Definitions Discrete Theory in a Nutshell. CS 523: Computer Graphics, Spring 2011

Shape Modeling. Differential Geometry Primer Smooth Definitions Discrete Theory in a Nutshell. CS 523: Computer Graphics, Spring 2011 CS 523: Computer Graphics, Spring 2011 Shape Modeling Differential Geometry Primer Smooth Definitions Discrete Theory in a Nutshell 2/15/2011 1 Motivation Geometry processing: understand geometric characteristics,

More information

High Order Nédélec Elements with local complete sequence properties

High Order Nédélec Elements with local complete sequence properties High Order Nédélec Elements with local complete sequence properties Joachim Schöberl and Sabine Zaglmayr Institute for Computational Mathematics, Johannes Kepler University Linz, Austria E-mail: {js,sz}@jku.at

More information

Metafor FE Software. 2. Operator split. 4. Rezoning methods 5. Contact with friction

Metafor FE Software. 2. Operator split. 4. Rezoning methods 5. Contact with friction ALE simulations ua sus using Metafor eao 1. Introduction 2. Operator split 3. Convection schemes 4. Rezoning methods 5. Contact with friction 1 Introduction EULERIAN FORMALISM Undistorted mesh Ideal for

More information

Subdivision Surfaces

Subdivision Surfaces Subdivision Surfaces 1 Geometric Modeling Sometimes need more than polygon meshes Smooth surfaces Traditional geometric modeling used NURBS Non uniform rational B-Spline Demo 2 Problems with NURBS A single

More information

Parameterization of Meshes

Parameterization of Meshes 2-Manifold Parameterization of Meshes What makes for a smooth manifold? locally looks like Euclidian space collection of charts mutually compatible on their overlaps form an atlas Parameterizations are

More information

Numerical Methods for PDEs : Video 11: 1D FiniteFebruary Difference 15, Mappings Theory / 15

Numerical 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 information

GEOMETRIC LIBRARY. Maharavo Randrianarivony

GEOMETRIC LIBRARY. Maharavo Randrianarivony GEOMETRIC LIBRARY Maharavo Randrianarivony During the last four years, I have maintained a numerical geometric library. The constituting routines, which are summarized in the following list, are implemented

More information

Meshless Modeling, Animating, and Simulating Point-Based Geometry

Meshless 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 information

B-spline Curves. Smoother than other curve forms

B-spline Curves. Smoother than other curve forms Curves and Surfaces B-spline Curves These curves are approximating rather than interpolating curves. The curves come close to, but may not actually pass through, the control points. Usually used as multiple,

More information

Finite Element Methods

Finite 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 information

implicit surfaces, approximate implicitization, B-splines, A- patches, surface fitting

implicit surfaces, approximate implicitization, B-splines, A- patches, surface fitting 24. KONFERENCE O GEOMETRII A POČÍTAČOVÉ GRAFICE ZBYNĚK ŠÍR FITTING OF PIECEWISE POLYNOMIAL IMPLICIT SURFACES Abstrakt In our contribution we discuss the possibility of an efficient fitting of piecewise

More information

C 1 Quintic Spline Interpolation Over Tetrahedral Partitions

C 1 Quintic Spline Interpolation Over Tetrahedral Partitions C 1 Quintic Spline Interpolation Over Tetrahedral Partitions Gerard Awanou and Ming-Jun Lai Abstract. We discuss the implementation of a C 1 quintic superspline method for interpolating scattered data

More information

A popular method for moving beyond linearity. 2. Basis expansion and regularization 1. Examples of transformations. Piecewise-polynomials and splines

A popular method for moving beyond linearity. 2. Basis expansion and regularization 1. Examples of transformations. Piecewise-polynomials and splines A popular method for moving beyond linearity 2. Basis expansion and regularization 1 Idea: Augment the vector inputs x with additional variables which are transformation of x use linear models in this

More information

ME 345: Modeling & Simulation. Introduction to Finite Element Method

ME 345: Modeling & Simulation. Introduction to Finite Element Method ME 345: Modeling & Simulation Introduction to Finite Element Method Examples Aircraft 2D plate Crashworthiness 2 Human Heart Gears Structure Human Spine 3 F.T. Fisher, PhD Dissertation, 2002 Fluid Flow

More information

An introduction to mesh generation Part IV : elliptic meshing

An 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 information

Fall CSCI 420: Computer Graphics. 4.2 Splines. Hao Li.

Fall CSCI 420: Computer Graphics. 4.2 Splines. Hao Li. Fall 2014 CSCI 420: Computer Graphics 4.2 Splines Hao Li http://cs420.hao-li.com 1 Roller coaster Next programming assignment involves creating a 3D roller coaster animation We must model the 3D curve

More information

NUMERICAL COUPLING BETWEEN DEM (DISCRETE ELEMENT METHOD) AND FEA (FINITE ELEMENTS ANALYSIS).

NUMERICAL COUPLING BETWEEN DEM (DISCRETE ELEMENT METHOD) AND FEA (FINITE ELEMENTS ANALYSIS). NUMERICAL COUPLING BETWEEN DEM (DISCRETE ELEMENT METHOD) AND FEA (FINITE ELEMENTS ANALYSIS). Daniel Schiochet Nasato - ESSS Prof. Dr. José Roberto Nunhez Unicamp Dr. Nicolas Spogis - ESSS Fabiano Nunes

More information

Efficient Storage and Processing of Adaptive Triangular Grids using Sierpinski Curves

Efficient Storage and Processing of Adaptive Triangular Grids using Sierpinski Curves Efficient Storage and Processing of Adaptive Triangular Grids using Sierpinski Curves Csaba Attila Vigh, Dr. Michael Bader Department of Informatics, TU München JASS 2006, course 2: Numerical Simulation:

More information

Parallel multi-frontal solver for isogeometric finite element methods on GPU

Parallel multi-frontal solver for isogeometric finite element methods on GPU Parallel multi-frontal solver for isogeometric finite element methods on GPU Maciej Paszyński Department of Computer Science, AGH University of Science and Technology, Kraków, Poland email: paszynsk@agh.edu.pl

More information

Numerical Algorithms

Numerical Algorithms Chapter 10 Slide 464 Numerical Algorithms Slide 465 Numerical Algorithms In textbook do: Matrix multiplication Solving a system of linear equations Slide 466 Matrices A Review An n m matrix Column a 0,0

More information

Parametric Surfaces. Michael Kazhdan ( /657) HB , FvDFH 11.2

Parametric Surfaces. Michael Kazhdan ( /657) HB , FvDFH 11.2 Parametric Surfaces Michael Kazhdan (601.457/657) HB 10.6 -- 10.9, 10.1 FvDFH 11.2 Announcement OpenGL review session: When: Wednesday (10/1) @ 7:00-9:00 PM Where: Olin 05 Cubic Splines Given a collection

More information

Chapter 3 Numerical Methods

Chapter 3 Numerical Methods Chapter 3 Numerical Methods Part 1 3.1 Linearization and Optimization of Functions of Vectors 1 Problem Notation 2 Outline 3.1.1 Linearization 3.1.2 Optimization of Objective Functions 3.1.3 Constrained

More information

PROGRAMMING OF MULTIGRID METHODS

PROGRAMMING OF MULTIGRID METHODS PROGRAMMING OF MULTIGRID METHODS LONG CHEN In this note, we explain the implementation detail of multigrid methods. We will use the approach by space decomposition and subspace correction method; see Chapter:

More information

PATCH TEST OF HEXAHEDRAL ELEMENT

PATCH TEST OF HEXAHEDRAL ELEMENT Annual Report of ADVENTURE Project ADV-99- (999) PATCH TEST OF HEXAHEDRAL ELEMENT Yoshikazu ISHIHARA * and Hirohisa NOGUCHI * * Mitsubishi Research Institute, Inc. e-mail: y-ishi@mri.co.jp * Department

More information

E and. L q. AE q L AE L. q L

E and. L q. AE q L AE L. q L STRUTURL NLYSIS [SK 43] EXERISES Q. (a) Using basic concepts, members towrds local axes is, E and q L, prove that the equilibrium equation for truss f f E L E L E L q E q L With f and q are both force

More information

WESTFÄLISCHE WILHELMS-UNIVERSITÄT MÜNSTER

WESTFÄLISCHE WILHELMS-UNIVERSITÄT MÜNSTER MÜNSTER Adaptive Discretization of Liftings for Curvature Regularization SIAM Conference on Imaging Science 2016, Albuquerque, NM Ulrich Hartleif, Prof. Dr. Benedikt Wirth May 23, 2016 Outline MÜNSTER

More information

Nodal Basis Functions for Serendipity Finite Elements

Nodal Basis Functions for Serendipity Finite Elements Nodal Basis Functions for Serendipity Finite Elements Andrew Gillette Department of Mathematics University of Arizona joint work with Michael Floater (University of Oslo) Andrew Gillette - U. Arizona Nodal

More information

Von Neumann Analysis for Higher Order Methods

Von Neumann Analysis for Higher Order Methods 1. Introduction Von Neumann Analysis for Higher Order Methods Von Neumann analysis is a widely used method to study how an initial wave is propagated with certain numerical schemes for a linear wave equation

More information

Interpolation & Polynomial Approximation. Cubic Spline Interpolation II

Interpolation & Polynomial Approximation. Cubic Spline Interpolation II Interpolation & Polynomial Approximation Cubic Spline Interpolation II Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University

More information

A New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes

A New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes A New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes Sanjay Kumar Khattri Department of Mathematics, University of Bergen, Norway sanjay@mi.uib.no http://www.mi.uib.no/ sanjay Abstract. Mesh

More information

Optimization of Execution Time, Energy Consumption, and Accuracy during Finite Element Method Simulations

Optimization of Execution Time, Energy Consumption, and Accuracy during Finite Element Method Simulations Optimization of Execution Time, Energy Consumption, and Accuracy during Finite Element Method Simulations Maciej Woźniak, Marcin Łoś, Leszek Siwik, Darusz Król, Maciej Paszyński Department of Computer

More information

Computergrafik. Matthias Zwicker. Herbst 2010

Computergrafik. Matthias Zwicker. Herbst 2010 Computergrafik Matthias Zwicker Universität Bern Herbst 2010 Today Curves NURBS Surfaces Parametric surfaces Bilinear patch Bicubic Bézier patch Advanced surface modeling Piecewise Bézier curves Each segment

More information

Curves and Surfaces 1

Curves and Surfaces 1 Curves and Surfaces 1 Representation of Curves & Surfaces Polygon Meshes Parametric Cubic Curves Parametric Bi-Cubic Surfaces Quadric Surfaces Specialized Modeling Techniques 2 The Teapot 3 Representing

More information

Anisotropic quality measures and adaptation for polygonal meshes

Anisotropic quality measures and adaptation for polygonal meshes Anisotropic quality measures and adaptation for polygonal meshes Yanqiu Wang, Oklahoma State University joint work with Weizhang Huang Oct 2015, POEMs, Georgia Tech. 1 / 23 Consider the following function

More information

Splines. Parameterization of a Curve. Curve Representations. Roller coaster. What Do We Need From Curves in Computer Graphics? Modeling Complex Shapes

Splines. Parameterization of a Curve. Curve Representations. Roller coaster. What Do We Need From Curves in Computer Graphics? Modeling Complex Shapes CSCI 420 Computer Graphics Lecture 8 Splines Jernej Barbic University of Southern California Hermite Splines Bezier Splines Catmull-Rom Splines Other Cubic Splines [Angel Ch 12.4-12.12] Roller coaster

More information

Introduction to Computer Graphics

Introduction to Computer Graphics Introduction to Computer Graphics 2016 Spring National Cheng Kung University Instructors: Min-Chun Hu 胡敏君 Shih-Chin Weng 翁士欽 ( 西基電腦動畫 ) Data Representation Curves and Surfaces Limitations of Polygons Inherently

More information

When implementing FEM for solving two-dimensional partial differential equations, integrals of the form

When 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 information

This is an author-deposited version published in: Eprints ID: 4362

This is an author-deposited version published in:   Eprints ID: 4362 This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral

More information

f xx + f yy = F (x, y)

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 information

Isogeometric Analysis Application to Car Crash Simulation

Isogeometric Analysis Application to Car Crash Simulation Isogeometric Analysis Application to Car Crash Simulation S. Bouabdallah 2, C. Adam 2,3, M. Zarroug 2, H. Maitournam 3 De Vinci Engineering Lab, École Supérieure d Ingénieurs Léonard de Vinci 2 Direction

More information

COMPUTER AIDED ENGINEERING DESIGN (BFF2612)

COMPUTER AIDED ENGINEERING DESIGN (BFF2612) COMPUTER AIDED ENGINEERING DESIGN (BFF2612) BASIC MATHEMATICAL CONCEPTS IN CAED by Dr. Mohd Nizar Mhd Razali Faculty of Manufacturing Engineering mnizar@ump.edu.my COORDINATE SYSTEM y+ y+ z+ z+ x+ RIGHT

More information

The Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations

The Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving

More information

Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications

Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/

More information

Radiosity. Johns Hopkins Department of Computer Science Course : Rendering Techniques, Professor: Jonathan Cohen

Radiosity. Johns Hopkins Department of Computer Science Course : Rendering Techniques, Professor: Jonathan Cohen Radiosity Radiosity Concept Global computation of diffuse interreflections among scene objects Diffuse lighting changes fairly slowly across a surface Break surfaces up into some number of patches Assume

More information

CS-184: Computer Graphics

CS-184: Computer Graphics CS-184: Computer Graphics Lecture #12: Curves and Surfaces Prof. James O Brien University of California, Berkeley V2007-F-12-1.0 Today General curve and surface representations Splines and other polynomial

More information

CS-184: Computer Graphics. Today

CS-184: Computer Graphics. Today CS-84: Computer Graphics Lecture #5: Curves and Surfaces Prof. James O Brien University of California, Berkeley V25F-5-. Today General curve and surface representations Splines and other polynomial bases

More information

To appear in Computer Methods in Applied Mechanics and Engineering. Isogeometric shape optimization of photonic crystals via Coons patches

To appear in Computer Methods in Applied Mechanics and Engineering. Isogeometric shape optimization of photonic crystals via Coons patches To appear in Computer Methods in Applied Mechanics and Engineering Isogeometric shape optimization of photonic crystals via Coons patches Xiaoping Qian Department of Mechanical, Materials and Aerospace

More information

An Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm

An Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm An Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm P.A. Sherar, C.P. Thompson, B. Xu, B. Zhong Abstract A new method is presented to deal with shape optimization problems.

More information

1 Exercise: Heat equation in 2-D with FE

1 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 information

REPORT DOCUMENTATION PAGE

REPORT DOCUMENTATION PAGE REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

Introduction p. 1 What Is Geometric Modeling? p. 1 Computer-aided geometric design Solid modeling Algebraic geometry Computational geometry

Introduction p. 1 What Is Geometric Modeling? p. 1 Computer-aided geometric design Solid modeling Algebraic geometry Computational geometry Introduction p. 1 What Is Geometric Modeling? p. 1 Computer-aided geometric design Solid modeling Algebraic geometry Computational geometry Representation Ab initio design Rendering Solid modelers Kinematic

More information

Mesh Processing Pipeline

Mesh Processing Pipeline Mesh Smoothing 1 Mesh Processing Pipeline... Scan Reconstruct Clean Remesh 2 Mesh Quality Visual inspection of sensitive attributes Specular shading Flat Shading Gouraud Shading Phong Shading 3 Mesh Quality

More information

Documentation for Numerical Derivative on Discontinuous Galerkin Space

Documentation 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 information

SOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS

SOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS SOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS JIAN LIANG AND HONGKAI ZHAO Abstract. In this paper we present a general framework for solving partial differential equations on manifolds represented

More information

G 2 Interpolation for Polar Surfaces

G 2 Interpolation for Polar Surfaces 1 G 2 Interpolation for Polar Surfaces Jianzhong Wang 1, Fuhua Cheng 2,3 1 University of Kentucky, jwangf@uky.edu 2 University of Kentucky, cheng@cs.uky.edu 3 National Tsinhua University ABSTRACT In this

More information

Adaptive boundary element methods in industrial applications

Adaptive 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 information

Computer Graphics Curves and Surfaces. Matthias Teschner

Computer 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 information

ICES REPORT March Isogeometric Finite Element Data Structures based on Bezier Extraction of NURBS

ICES REPORT March Isogeometric Finite Element Data Structures based on Bezier Extraction of NURBS ICES REPORT 0-08 March 200 Isogeometric Finite Element Data Structures based on Bezier Extraction of NURBS by Michael J. Borden, Michael A. Scott, John A. Evans, and Thomas J.R. Hughes The Institute for

More information

Adaptive Surface Modeling Using a Quadtree of Quadratic Finite Elements

Adaptive Surface Modeling Using a Quadtree of Quadratic Finite Elements Adaptive Surface Modeling Using a Quadtree of Quadratic Finite Elements G. P. Nikishkov University of Aizu, Aizu-Wakamatsu 965-8580, Japan niki@u-aizu.ac.jp http://www.u-aizu.ac.jp/ niki Abstract. This

More information

Chapter 7: Computation of the Camera Matrix P

Chapter 7: Computation of the Camera Matrix P Chapter 7: Computation of the Camera Matrix P Arco Nederveen Eagle Vision March 18, 2008 Arco Nederveen (Eagle Vision) The Camera Matrix P March 18, 2008 1 / 25 1 Chapter 7: Computation of the camera Matrix

More information

CS 523: Computer Graphics, Spring Shape Modeling. Differential Geometry of Surfaces

CS 523: Computer Graphics, Spring Shape Modeling. Differential Geometry of Surfaces CS 523: Computer Graphics, Spring 2011 Shape Modeling Differential Geometry of Surfaces Andrew Nealen, Rutgers, 2011 2/22/2011 Differential Geometry of Surfaces Continuous and Discrete Motivation Smoothness

More information

Multi-Mesh CFD. Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech

Multi-Mesh CFD. Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech Multi-Mesh CFD Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech cjroy@vt.edu May 21, 2014 CCAS Program Review, Columbus, OH 1 Motivation Automated

More information