Algebraic Flux Correction Schemes for High-Order B-Spline Based Finite Element Approximations
|
|
- Amber Kristina Chase
- 6 years ago
- Views:
Transcription
1 Algebraic Flux Correction Schemes for High-Order B-Spline Based Finite Element Approximations Numerical Analysis group Matthias Möller 7th DIAM onderwijs- en onderzoeksdag, 5th November 2014
2 Outline 1 Motivation Finite elements in a nutshell 2 Introduction to B-splines Knot insertion (h-refinement) Order elevation (p-refinement) Geometric mapping Definition of ansatz spaces Properties of B-splines 3 Constrained data projection Numerical examples 4 Constrained transport Numerical examples M. Möller (NA group) AFC for B-Spline based FEM 2 / 34
3 Finite elements in a nutshell Strong problem: find u C k (Ω) such that Lu = f in Ω + bc s Weighted residual formulation: find u V such that w[lu f] dx = 0 w W Ω Boundary conditions: V (trial) and W (test spaces) contain essential bc s natural bc s are incorporated via integration by parts M. Möller (NA group) AFC for B-Spline based FEM 3 / 34
4 Finite elements in a nutshell, cont d Galerkin finite elements: choose finite-dimensional spaces V h := {ϕ j } V and W h := {φ i } W and find u h = j u jϕ j V h such that φ i [Lu h f] dx = 0 i = 1,..., dim(w h ) Ω h neglecting complications due to bc s this yields [ ] φ i Lϕ j dx u j = φ i f dx i = 1,..., dim(w h ) Ω h Ω h j M. Möller (NA group) AFC for B-Spline based FEM 4 / 34
5 Problem Poor approximation of discontinuities/steep gradients if standard Galerkin methods are used without proper stabilization L 2 -projection convective transport M. Möller (NA group) AFC for B-Spline based FEM 5 / 34
6 Problem Poor approximation of discontinuities/steep gradients if standard Galerkin methods are used with proper stabilization L 2 -projection convective transport M. Möller (NA group) AFC for B-Spline based FEM 5 / 34
7 Algebraic Flux Correction Methodology based on algebraic design criteria to derive robust and accurate high-resolution finite element schemes for Constrained data projection [5] Convection-dominated transport processes [3, 5, 6, 7, 10, 11] Anisotropic diffusion processes [3, 9] Processes with maximum-packing limit [4]... Many success stories published... ;-) M. Möller (NA group) AFC for B-Spline based FEM 6 / 34
8 Algebraic Flux Correction Methodology based on algebraic design criteria to derive robust and accurate high-resolution finite element schemes for Constrained data projection [5] Convection-dominated transport processes [3, 5, 6, 7, 10, 11] Anisotropic diffusion processes [3, 9] Processes with maximum-packing limit [4]... Many success stories published... ;-) mostly for P 1 and Q 1 finite elements ;-( M. Möller (NA group) AFC for B-Spline based FEM 6 / 34
9 Algebraic Flux Correction, cont d Extension of AFC to P1 / Q 1 elements [12]: P 2 elements [2]: CR-FE satisfy the necessary prerequisites [...] fail completely [...] yielding overdiffusive approximate solutions. RT-FE provides an accurate resolution [...] if the integral mean value based variant is adopted. In summary, algebraic flux correction for quadratic finite elements seems to be feasible but gives rise to many challenging open problems. Objective: to extend AFC to high-order B-spline basis functions M. Möller (NA group) AFC for B-Spline based FEM 7 / 34
10 B-splines in a nutshell Define knot vector Ξ = (ξ 1, ξ 2,..., ξ n+p+1 ) as a sequence of non-decreasing coordinates in the parameter space Ω 0 = [0, 1]: ξ i R is the i th knot with index i = 1, 2,..., n + p + 1 p is the polynomial order of the B-splines n is the number of B-spline functions Cox-de Boor recursion formula for N i,p : Ω 0 R p = 0 : N i,0 = { 1 if ξi ξ < ξ i+1 0 otherwise define 0 0 = 0 p > 0 : N i,p = ξ ξ i N i,p 1 + ξ i+p+1 ξ N i+1,p 1 ξ i+p ξ i ξ i+p+1 ξ i+1 M. Möller (NA group) AFC for B-Spline based FEM 8 / 34
11 Knot insertion (h-refinement) Ξ = (0, 0, 1, 1), p = 1 Ξ = (0, 0, 0.5, 1, 1), p = one element N 0,1 N 1, N 0,1 N 1,1 N 2,1 two elements, C 0 -continuity In general we have C p m i -continuity across element boundaries, where m i is the multiplicity of the value of ξ i in the knot vector. M. Möller (NA group) AFC for B-Spline based FEM 9 / 34
12 Knot insertion, cont d Ξ = (0, 0, 1, 1), p = 1 Ξ = (0, 0, 0.5, 0.5, 1, 1), p = one element N 0,1 N 1, N 0,1 N 1,1 N 2,1 N 3,1 two elements, discontinuity In general we have C p m i -continuity across element boundaries, where m i is the multiplicity of the value of ξ i in the knot vector. M. Möller (NA group) AFC for B-Spline based FEM 10 / 34
13 Order elevation 1 Ξ = (0, 0, 1, 1), p = 1 1 Ξ = (0, 0, 0, 1, 1, 1), p = one element N 0,1 N 1, one element N 0,2 N 1,2 N 2,2 In contrast to standard Lagrange finite element basis functions, the B-spline functions never become negative over their support. M. Möller (NA group) AFC for B-Spline based FEM 11 / 34
14 Nonuniform continuity at element boundaries th -order B-spline functions Ξ = (0, 0, 0, 0, 0, 1/3, 2/3, 2/3, 2/3, 1, 1, 1, 1, 1) }{{}}{{}}{{}}{{} discontinuity C 3 C 1 discontinuity N 0,4 N 1,4 N 2,4 N 3,4 N 4,4 N 5,4 N 6,4 N 7,4 N 8,4 M. Möller (NA group) AFC for B-Spline based FEM 12 / 34
15 Nonuniform continuity at element boundaries First derivatives of 4 th -order B-spline functions Ξ = (0, 0, 0, 0, 0, 1/3, 2/3, 2/3, 2/3, 1, 1, 1, 1, 1) }{{}}{{}}{{}}{{} discontinuity C 3 C 1 discontinuity N 0,4 N 1,4 N 2,4 N 3,4 N 4,4 N 5,4 N 6,4 N 7,4 N 8,4 M. Möller (NA group) AFC for B-Spline based FEM 13 / 34
16 Geometric mapping Define mapping between parameter space Ω 0 = [0, 1] and the computational domain Ω using control points p i R d, d = 1, 2, 3 G : Ω 0 Ω, G = n N i,p p i i=1 Examples: G : [0, 1] [a, b] R G : [0, 1] curve in R 2 or R 3 (linear mapping) (work in progress) M. Möller (NA group) AFC for B-Spline based FEM 14 / 34
17 Ansatz spaces Construct ansatz space from B-spline basis functions V h (Ω 0, p, Ξ, G) = span{ϕ i (x) = N i,p G 1 (x)} Approximate the solution the standard way u(x) u h (x) = n ϕ i (x)u i, i=1 x Ω Approximate fluxes by Fletcher s group formulation, e.g., v(x)u(x) (vu) h (x) = n ϕ i (x)(v i u i ), i=1 x Ω M. Möller (NA group) AFC for B-Spline based FEM 15 / 34
18 Properties of B-splines Derivative of p th order B-spline is a B-spline of order p 1 N i,p = p p N i,p 1 N i+1,p 1 ξ i+p ξ i ξ i+p+1 ξ i+1 B-splines form a partition of unity. That is, for all ξ [a, b] n N i,p = 1 i=1 n N i,p = 0 i=1 AFC-1: Edge-wise flux decomposition b n c ij = ϕ i ϕ j dx, c ij = 0 a j=1 n j=1 c ij u j = j i c ij (u j u i ) M. Möller (NA group) AFC for B-Spline based FEM 16 / 34
19 Properties of B-splines, cont d B-splines of order p have compact support supp N i,p = [ξ i, ξ i+p+1 ), i = 1,..., n B-splines are strictly positive over the interior of their support N i,p > 0 for ξ (ξ i, ξ i+p+1 ), i = 1,..., n AFC-2: positive consistent and lumped mass matrices b n m ij = ϕ i ϕ j dx > 0 m i = m ij > 0 a j=1 M. Möller (NA group) AFC for B-Spline based FEM 17 / 34
20 Constrained data projection Find u L 2 ([a, b]) s.t. b a w(u f)dx = 0 w L 2 ([a, b]) Consistent L 2 -projection n b m ij u H j = ϕ i f dx j=1 a Lumped L 2 -projection m i u L i = b a ϕ i f dx Constrained L 2 -projection [5] u i = u L i + 1 α ij fij H, m i 0 α ij = α ji 1, fji H = fij H j i M. Möller (NA group) AFC for B-Spline based FEM 18 / 34
21 Symmetric flux limiting algorithm [5] 1 Prelimited raw antidiffusive fluxes { fij H m ij (u H i u H j = ), if (uh i u H j )(ul i u L j ) > 0 0, otherwise (!) 2 Bounds and antidiffusive increments Q ± i = max min j i (0, u L i u L j ), P ± i = max min j i (0, f ij ) 3 Nodal and edge-wise limiting coefficients R ± i = m iq ± i P ± i, α ij = { min(r + i, R j ), if f ij > 0 min(r + j, R i ), if f ij < 0 M. Möller (NA group) AFC for B-Spline based FEM 19 / 34
22 Test case: Semi-ellipse of McDonald p = 1, n = 32 1 L 1 -/L 2 -errors u H u 1 = u L u 1 = u u 1 = u H u 2 = u L u 2 = u u 2 = M. Möller (NA group) AFC for B-Spline based FEM 20 / 34
23 Test case: semi-ellipse of McDonald p = 2, n = L 1 -/L 2 -errors u H u 1 = u L u 1 = u u 1 = u H u 2 = u L u 2 = u u 2 = M. Möller (NA group) AFC for B-Spline based FEM 21 / 34
24 Test case: semi-ellipse of McDonald p = 3, n = 32 1 L 1 -/L 2 -errors u H u 1 = u L u 1 = u u 1 = u H u 2 = u L u 2 = u u 2 = M. Möller (NA group) AFC for B-Spline based FEM 22 / 34
25 Test case: semi-ellipse of McDonald p = 4, n = 32 1 L 1 -/L 2 -errors u H u 1 = u L u 1 = u u 1 = u H u 2 = u L u 2 = u u 2 = M. Möller (NA group) AFC for B-Spline based FEM 23 / 34
26 Thought experiment: What are ideal knots? M. Möller (NA group) AFC for B-Spline based FEM 24 / 34
27 Thought experiment: What are ideal knots? L 1 -/L 2 -errors u H u 1 = u u 1 = u H u 2 = u u 2 = Deduce ideal knots from nodal correction factors or residual Ξ = (0, 0, 0, 0, 0.1, 0.2, 0.2, 0.2, 0.2,0.3, 0.4, 0.4, 0.4, 0.4,0.5, 0.6, 0.6, 0.6, 0.6,0.7, 0.7, 0.7,0.8, 0.8, 0.8, 0.8,0.9, 1, 1, 1, 1) Use smoothness indicator [8] to avoid peak clipping. M. Möller (NA group) AFC for B-Spline based FEM 24 / 34
28 Constrained transport b Find u HD 1 ([a, b]) s.t. w(vu) x + dw x u x dx = 0 j=1 j=1 a w H 1 0 ([a, b]) Galerkin scheme Discrete upwind scheme n (k ij + s ij )u H j = 0 n (k ij + d ij + s ij )u L j = 0 b k ij = v j ϕ i ϕ jdx, a s ij = d b a ϕ iϕ jdx, d ij = max(k ij, 0, k ji ) High-resolution TVD-type scheme [11] n ij + d ij + s ij )u j=1(k j + α ij f ij = 0, j i f ji = f ij M. Möller (NA group) AFC for B-Spline based FEM 25 / 34
29 Test case: steady convection-diffusion v d = p = 1, n = p = 2, n = p = 3, n = u H u 1 = u L u 1 = u u 1 = u H u 1 = u L u 1 = u u 1 = u H u 1 = u L u 1 = u u 1 = M. Möller (NA group) AFC for B-Spline based FEM 26 / 34
30 Summary Algebraic flux correction concept has been generalized to higher-order approximations based on B-spline bases Original lowest-order approximation is naturally included Nodal correction factors/residual provide information to locally reduce inter-element continuity by increasing knot multiplicity Peak clipping at smooth extrema is prevented by locally deactivating the flux limiter using the smoothness indicator [8] M. Möller (NA group) AFC for B-Spline based FEM 27 / 34
31 Current and future research Analysis for general geometric mappings G : Ω 0 Ω Extension to multi-dimensions by tensor-product construction G(ξ, η, ζ) = n m l N i,p N j,q (η)n k,r (ζ)p ijk i=1 j=1 k=1 Embed local B-spline based AFC-scheme into global outer Galerkin method with unstructured quad/hexa macro mesh Exploit potential of fully structured data per macro element A. Jaeschke (COSSE-MSc), S-R. Janssen (DD: AM-LR) M. Möller (NA group) AFC for B-Spline based FEM 28 / 34
32 Outlook: Isogeometric Analysis [1] Poor approximation of curved boundaries (with low-order FEs) Ω Ω h IgA approach adopts the same (hierarchical) B-spline, NURBS, etc. basis functions for the approximate solution u h u and for exactly representing the geometry Ω h = Ω Our activities in IgA: bi-weekly PhD-seminar (participants from LR) PhD-candidate in CSC-15 and/or EU-project (?) M. Möller (NA group) AFC for B-Spline based FEM 29 / 34
33 References I T.J.R. Hughes, J.A. Cottrell, and Y. Bazilevs. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Computer Methods in Applied Mechanics and Engineering, 194: , D. Kuzmin. On the design of algebraic flux correction schemes for quadratic finite elements. Journal of Computational and Applied Mathematics, 218(1):79 87, Special Issue: Finite Element Methods in Engineering and Science (FEMTEC 2006) Special Issue: Finite Element Methods in Engineering and Science (FEMTEC 2006). M. Möller (NA group) AFC for B-Spline based FEM 30 / 34
34 References II D. Kuzmin. Linearity-preserving flux correction and convergence acceleration for constrained Galerkin schemes. Journal of Computational and Applied Mathematics, 236(9): , D. Kuzmin and Y. Gorb. A flux-corrected transport algorithm for handling the close-packing limit in dense suspensions. Journal of Computational and Applied Mathematics, 236(18): , {FEMTEC} 2011: 3rd International Conference on Computational Methods in Engineering and Science, May 913, M. Möller (NA group) AFC for B-Spline based FEM 31 / 34
35 References III D. Kuzmin, M. Möller, J.N. Shadid, and M. Shashkov. Failsafe flux limiting and constrained data projections for equations of gas dynamics. Journal of Computational Physics, 229(23): , D. Kuzmin, M. Möller, and S. Turek. Multidimensional FEM-FCT schemes for arbitrary time stepping. International Journal for Numerical Methods in Fluids, 42(3): , D. Kuzmin, M. Möller, and S. Turek. High-resolution FEMFCT schemes for multidimensional conservation laws. Computer Methods in Applied Mechanics and Engineering, 193(4547): , M. Möller (NA group) AFC for B-Spline based FEM 32 / 34
36 References IV D. Kuzmin and F. Schieweck. A parameter-free smoothness indicator for high-resolution finite element schemes. Central European Journal of Mathematics, 11(8): , D. Kuzmin, M.J. Shashkov, and D. Svyatskiy. A constrained finite element method satisfying the discrete maximum principle for anisotropic diffusion problems. Journal of Computational Physics, 228(9): , D. Kuzmin and S. Turek. Flux correction tools for finite elements. Journal of Computational Physics, 175(2): , M. Möller (NA group) AFC for B-Spline based FEM 33 / 34
37 References V D. Kuzmin and S. Turek. High-resolution FEM-TVD schemes based on a fully multidimensional flux limiter. Journal of Computational Physics, 198(1): , M. Möller. Algebraic flux correction for nonconforming finite element discretizations of scalar transport problems. Computing, 95(5): , M. Möller (NA group) AFC for B-Spline based FEM 34 / 34
Recent Developments in Isogeometric Analysis with Solid Elements in LS-DYNA
Recent Developments in Isogeometric Analysis with Solid Elements in LS-DYNA Liping Li David Benson Attila Nagy Livermore Software Technology Corporation, Livermore, CA, USA Mattia Montanari Nik Petrinic
More informationEfficient Multigrid based solvers for Isogeometric Analysis
Efficient Multigrid based solvers for Isogeometric Analysis R. Tielen, M. Möller and C. Vuik Delft Institute of Applied Mathematics (DIAM) Numerical Analysis / 22 Isogeometric Analysis (IgA) Extension
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 informationFrequency Domain Analysis for Isogeometric Element in LS-DYNA
Frequency Domain Analysis for Isogeometric Element in LS-DYNA Liping Li 1, Yun Huang 1, Zhe Cui 1, Stefan Hartmann 2, David J. Benson 3 1 Livermore Software Technology Corporation, Livermore, CA, USA 2
More informationIsogeometric Analysis for Compressible Flows with Application in Turbomachinery
Isogeometric Analysis for Compressible Flows with Application in Turbomachinery by Andrzej Jaeschke to obtain the degree of Master of Science at the Delft University of Technology, to be defended publicly
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 informationThe NURBS and GeoPDEs packages Octave software for research on IGA
The NURBS and GeoPDEs packages Octave software for research on IGA Rafael Vázquez IMATI Enrico Magenes, Pavia Consiglio Nazionale delle Ricerche R. Vázquez (IMATI-CNR Italy) The NURBS and GeoPDEs packages
More informationIsogeometric Analysis of Fluid-Structure Interaction
Isogeometric Analysis of Fluid-Structure Interaction Y. Bazilevs, V.M. Calo, T.J.R. Hughes Institute for Computational Engineering and Sciences, The University of Texas at Austin, USA e-mail: {bazily,victor,hughes}@ices.utexas.edu
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 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 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 informationMAFELAP Mini-Symposium: Advanced FEM methodologies and isogeometric analysis. Organisers: Matthias Möller and Bernd Simeon
MAFELAP 2016 Conference on the Mathematics of Finite Elements and Applications 14 17 June 2016 Mini-Symposium: Advanced FEM methodologies and isogeometric analysis Organisers: Matthias Möller and Bernd
More informationParallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes
Parallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes Stefan Vater 1 Kaveh Rahnema 2 Jörn Behrens 1 Michael Bader 2 1 Universität Hamburg 2014 PDES Workshop 2 TU München Partial
More informationAMS527: 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 informationNodal 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 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 informationMODELLING STIFFENED LIGHTWEIGHT STRUCTURES WITH ISOGEOMETRIC ANALYSIS VIA MORTAR METHODS
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationA Study of Isogeometric Analysis for Scalar Convection-Diffusion Equations
A Study of Isogeometric Analysis for Scalar Convection-Diffusion Equations Volker John a,b,, Liesel Schumacher c, a Weierstrass Institute for Applied Analysis and Stochastics, Leibniz Institute in Forschungsverbund
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 informationALE METHODS FOR DETERMINING STATIONARY SOLUTIONS OF METAL FORMING PROCESSES
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2000 Barcelona, 11-14 September 2000 c ECCOMAS ALE METHODS FOR DETERMINING STATIONARY SOLUTIONS OF METAL FORMING PROCESSES
More informationOn the high order FV schemes for compressible flows
Applied and Computational Mechanics 1 (2007) 453-460 On the high order FV schemes for compressible flows J. Fürst a, a Faculty of Mechanical Engineering, CTU in Prague, Karlovo nám. 13, 121 35 Praha, Czech
More informationICES 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 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 informationPolynomial Splines over locally refined box-partitions
Polynomial Splines over locally refined box-partitions Tor Dokken*, Tom Lyche +, Kjell Fredrik Pettersen* *SINTEF and + University of Oslo The work is partly funded by the European Union through The TERRIFIC
More informationAN ISOGEOMETRIC REISSNER-MINDLIN SHELL WITH LAGRANGE BASIS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationNUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.
Convergent Science White Paper COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. This document contains information that is proprietary to Convergent Science. Public dissemination of this document
More informationAN A PRIORI MODEL REDUCTION FOR ISOGEOMETRIC BOUNDARY ELEMENT METHOD
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
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 informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationNonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws Ivan Christov Bojan Popov Department of Mathematics, Texas A&M University, College Station, Texas
More informationTopology optimization of heat conduction problems
Topology optimization of heat conduction problems Workshop on industrial design optimization for fluid flow 23 September 2010 Misha Marie Gregersen Anton Evgrafov Mads Peter Sørensen Technical University
More informationA Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods
A Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods Scott A. Sarra, Derek Sturgill Marshall University, Department of Mathematics, One John Marshall Drive, Huntington
More informationA new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws
HYP 2012, Padova A new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws Argiris I. Delis & Ioannis K. Nikolos (TUC)
More informationNonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws Ivan Christov 1,* Bojan Popov 1 Peter Popov 2 1 Department of Mathematics, 2 Institute for Scientific
More informationAPPROXIMATING PDE s IN L 1
APPROXIMATING PDE s IN L 1 Veselin Dobrev Jean-Luc Guermond Bojan Popov Department of Mathematics Texas A&M University NONLINEAR APPROXIMATION TECHNIQUES USING L 1 Texas A&M May 16-18, 2008 Outline 1 Outline
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationICES REPORT February The Institute for Computational Engineering and Sciences The University of Texas at Austin Austin, Texas 78712
ICES REPORT 3-03 February 03 Isogeometric Collocation: Cost Comparison with Galerkin Methods and Extension to Adaptive Hierarchical NURBS Discretizations by Dominik Schillinger, John A. Evans, Alessandro
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationThird Order WENO Scheme on Three Dimensional Tetrahedral Meshes
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 2-4, pp. 86-848 Commun. Comput. Phys. February 29 Third Order WENO Scheme on Three Dimensional Tetrahedral Meshes Yong-Tao Zhang 1, and Chi-Wang Shu
More informationf 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 informationRobust Poisson Surface Reconstruction
Robust Poisson Surface Reconstruction V. Estellers, M. Scott, K. Tew, and S. Soatto Univeristy of California, Los Angeles Brigham Young University June 2, 2015 1/19 Goals: Surface reconstruction from noisy
More informationComputergrafik. 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 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 informationCollocation and optimization initialization
Boundary Elements and Other Mesh Reduction Methods XXXVII 55 Collocation and optimization initialization E. J. Kansa 1 & L. Ling 2 1 Convergent Solutions, USA 2 Hong Kong Baptist University, Hong Kong
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 informationUNSTRUCTURED GRIDS ON NURBS SURFACES. The surface grid can be generated either in a parameter. surfaces. Generating grids in a parameter space is
UNSTRUCTURED GRIDS ON NURBS SURFACES Jamshid Samareh-Abolhassani 1 Abstract A simple and ecient computational method is presented for unstructured surface grid generation. This method is built upon an
More informationHigh Order Weighted Essentially Non-Oscillatory Schemes for Convection. Dominated Problems. Chi-Wang Shu 1
High Order Weighted Essentially Non-Oscillatory Schemes for Convection Dominated Problems Chi-Wang Shu Division of Applied Mathematics, Brown University, Providence, Rhode Island 09 ABSTRACT High order
More informationChapter 6. Petrov-Galerkin Formulations for Advection Diffusion Equation
Chapter 6 Petrov-Galerkin Formulations for Advection Diffusion Equation In this chapter we ll demonstrate the difficulties that arise when GFEM is used for advection (convection) dominated problems. Several
More informationComputergrafik. Matthias Zwicker Universität Bern Herbst 2016
Computergrafik Matthias Zwicker Universität Bern Herbst 2016 Today Curves NURBS Surfaces Parametric surfaces Bilinear patch Bicubic Bézier patch Advanced surface modeling 2 Piecewise Bézier curves Each
More informationOn the thickness of discontinuities computed by THINC and RK schemes
The 9th Computational Fluid Dynamics Symposium B7- On the thickness of discontinuities computed by THINC and RK schemes Taku Nonomura, ISAS, JAXA, Sagamihara, Kanagawa, Japan, E-mail:nonomura@flab.isas.jaxa.jp
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 informationInteractive modelling and simulation using blending techniques. Aleksander Pedersen
Interactive modelling and simulation using blending techniques Aleksander Pedersen NUC R&D Group in mathematical and geometrical modeling, numerical simulations, programming and visualization, Narvik University
More informationarxiv: v2 [cs.na] 27 Sep 2013
Isogeometric analysis: an overview and computer implementation aspects Vinh Phu Nguyen a,1,, Stéphane P.A. Bordas a,2, Timon Rabczuk b,3 a School of Engineering, Institute of Mechanics and Advanced Materials,
More informationDomain Decomposition and hp-adaptive Finite Elements
Domain Decomposition and hp-adaptive Finite Elements Randolph E. Bank 1 and Hieu Nguyen 1 Department of Mathematics, University of California, San Diego, La Jolla, CA 9093-011, USA, rbank@ucsd.edu. Department
More informationAn 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 informationSimulation in Computer Graphics. Particles. Matthias Teschner. Computer Science Department University of Freiburg
Simulation in Computer Graphics Particles Matthias Teschner Computer Science Department University of Freiburg Outline introduction particle motion finite differences system of first order ODEs second
More informationCOMPUTER AIDED GEOMETRIC DESIGN. Thomas W. Sederberg
COMPUTER AIDED GEOMETRIC DESIGN Thomas W. Sederberg January 31, 2011 ii T. W. Sederberg iii Preface This semester is the 24 th time I have taught a course at Brigham Young University titled, Computer Aided
More informationCurrent Status of Isogeometric Analysis in LS-DYNA
Current Status of Isogeometric Analysis in LS-DYNA Stefan Hartmann Developer Forum, September 24 th, 2013, Filderstadt, Germany Development in cooperation with: D.J. Benson: Professor of Applied Mechanics,
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 informationIsogeometric Simulation and Shape Optimization with Applications to Electrical Machines
Isogeometric Simulation and Shape Optimization with Applications to Electrical Machines Peter Gangl, Ulrich Langer, Angelos Mantzaflaris, Rainer Schneckenleitner G+S Report No. 75 September 2018 Isogeometric
More informationRay-casting Algebraic Surfaces using the Frustum Form. Eurographics 2008 Crete, Thursday April 17.
Ray-casting Algebraic Surfaces using the Frustum Form Martin Reimers Johan Seland Eurographics 2008 Crete, Thursday April 17. Algebraic Surfaces Zero set of polynomial f : R 3 R f (x, y, z) = f ijk x i
More informationImportant Properties of B-spline Basis Functions
Important Properties of B-spline Basis Functions P2.1 N i,p (u) = 0 if u is outside the interval [u i, u i+p+1 ) (local support property). For example, note that N 1,3 is a combination of N 1,0, N 2,0,
More informationDEVELOPMENTS OF PARALLEL CURVED MESHING FOR HIGH-ORDER FINITE ELEMENT SIMULATIONS
DEVELOPMENTS OF PARALLEL CURVED MESHING FOR HIGH-ORDER FINITE ELEMENT SIMULATIONS By Qiukai Lu A Thesis Submitted to the Graduate Faculty of Rensselaer Polytechnic Institute in Partial Fulfillment of the
More informationIsogeometric Analysis of Solids & Structures Small and Finite Deformation Applications
Trends & Challenges in Computational Mechanics Padua, Italy - September 12-14, 2011 Isogeometric Analysis of Solids & Structures Small and Finite Deformation Applications Robert L. Taylor Department of
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 informationA Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere
A Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere Ram Nair Computational and Information Systems Laboratory (CISL) National Center for Atmospheric Research
More informationA NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION OF FLOW DOMAINS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 2018, Glasgow, UK A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION
More informationNIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011
NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011 First-Order Hyperbolic System Method If you have a CFD book for hyperbolic problems, you have a CFD book for all problems.
More informationICES REPORT January D. Schillinger, L. Dede, M.A. Scott, J.A. Evans, M.J. Borden, E. Rank, T.J.R. Hughes
ICES REPORT 12-05 January 2012 An Isogeometric Design-through-analysis Methodology based on Adaptive Hierarchical Refinement of NURBS, Immersed Boundary Methods, and T-spline CAD Surfaces by D. Schillinger,
More informationLiterature Report. Daniël Pols. 23 May 2018
Literature Report Daniël Pols 23 May 2018 Applications Two-phase flow model The evolution of the momentum field in a two phase flow problem is given by the Navier-Stokes equations: u t + u u = 1 ρ p +
More informationTopology optimization in B-spline space
Topology optimization in B-spline space Xiaoping Qian Mechanical, Materials, and Aerospace Engineering Department Illinois Institute of Technology Chicago, IL 60062, USA Email: qian@iit.edu Highlights
More informationCURVILINEAR MESH GENERATION IN 3D
CURVILINEAR MESH GENERATION IN 3D Saikat Dey, Robert M. O'Bara 2 and Mark S. Shephard 2 SFA Inc. / Naval Research Laboratory, Largo, MD., U.S.A., dey@cosmic.nrl.navy.mil 2 Scientific Computation Research
More informationA story about Non Uniform Rational B-Splines. E. Shcherbakov
A story about Non Uniform Rational B-Splines E. Shcherbakov Speakers 09-06: B-spline curves (W. Dijkstra) 16-06: NURBS (E. Shcherbakov) 30-06: B-spline surfaces (M. Patricio) Seminar 16-06-2004 2 Outline
More informationSPLINE-BASED MESHING TECHNIQUES FOR INDUSTRIAL APPLICATIONS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK SPLINE-BASED MESHING TECHNIQUES FOR INDUSTRIAL APPLICATIONS
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 informationIsoGeometric Analysis: Be zier techniques in Numerical Simulations
IsoGeometric Analysis: Be zier techniques in Numerical Simulations Ahmed Ratnani IPP, Garching, Germany July 30, 2015 1/ 1 1/1 Outline Motivations Computer Aided Design (CAD) Zoology Splines/NURBS Finite
More informationThe Level Set Method applied to Structural Topology Optimization
The Level Set Method applied to Structural Topology Optimization Dr Peter Dunning 22-Jan-2013 Structural Optimization Sizing Optimization Shape Optimization Increasing: No. design variables Opportunity
More informationRational Bezier Curves
Rational Bezier Curves Use of homogeneous coordinates Rational spline curve: define a curve in one higher dimension space, project it down on the homogenizing variable Mathematical formulation: n P(u)
More informationDaubechies Wavelets and Interpolating Scaling Functions and Application on PDEs
Daubechies Wavelets and Interpolating Scaling Functions and Application on PDEs R. Schneider F. Krüger TUB - Technical University of Berlin November 22, 2007 R. Schneider F. Krüger Daubechies Wavelets
More informationA NUMERICAL SOLUTION OF ONE DIMENSIONAL HEAT EQUATION USING CUBIC B-SPLINE BASIS FUNCTIONS
A NUMERICAL SOLUTION OF ONE DIMENSIONAL HEAT EQUATION USING CUBIC B-SPLINE BASIS FUNCTIONS 1 SHARANJEET DHAWAN, 2 SHEO KUMAR 1 Department of mathematics 2 Dr. B.R.Ambedkar National Institute of Technology,
More informationThe multi-level hp-method for three-dimensional problems: dynamically changing high-order mesh refinement with arbitrary hanging nodes
The multi-level hp-method for three-dimensional problems: dynamically changing high-order mesh refinement with arbitrary hanging nodes Nils Zander 1, Tino Bog 1, Mohamed Elhaddad 1, Felix Frischmann 1,
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 informationAPPM/MATH Problem Set 4 Solutions
APPM/MATH 465 Problem Set 4 Solutions This assignment is due by 4pm on Wednesday, October 16th. You may either turn it in to me in class on Monday or in the box outside my office door (ECOT 35). Minimal
More informationFluent User Services Center
Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume
More informationAlmost Curvature Continuous Fitting of B-Spline Surfaces
Journal for Geometry and Graphics Volume 2 (1998), No. 1, 33 43 Almost Curvature Continuous Fitting of B-Spline Surfaces Márta Szilvási-Nagy Department of Geometry, Mathematical Institute, Technical University
More informationNumerical simulation techniques for the efficient and accurate treatment of local fluidic transport processes together with chemical reactions
Numerical simulation techniques for the efficient and accurate treatment of local fluidic transport processes together with chemical reactions Stefan Turek,, Chair of Applied Mathematics and Numerics,
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
More informationTo 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 informationAsynchronous OpenCL/MPI numerical simulations of conservation laws
Asynchronous OpenCL/MPI numerical simulations of conservation laws Philippe HELLUY 1,3, Thomas STRUB 2. 1 IRMA, Université de Strasbourg, 2 AxesSim, 3 Inria Tonus, France IWOCL 2015, Stanford Conservation
More informationReal time Ray-Casting of Algebraic Surfaces
Real time Ray-Casting of Algebraic Surfaces Martin Reimers Johan Seland Center of Mathematics for Applications University of Oslo Workshop on Computational Method for Algebraic Spline Surfaces Thursday
More informationFEM techniques for interfacial flows
FEM techniques for interfacial flows How to avoid the explicit reconstruction of interfaces Stefan Turek, Shu-Ren Hysing (ture@featflow.de) Institute for Applied Mathematics University of Dortmund Int.
More informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationAdaptive-Mesh-Refinement Pattern
Adaptive-Mesh-Refinement Pattern I. Problem Data-parallelism is exposed on a geometric mesh structure (either irregular or regular), where each point iteratively communicates with nearby neighboring points
More informationMetafor 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 informationModelling of Isogeometric Analysis for Plane Stress Problem Using Matlab GUI Norliyana Farzana Binti Zulkefli, Ahmad Razin Bin Zainal Abidin
Modelling of Isogeometric Analysis for Plane Stress Problem Using Matlab GUI Norliyana Farzana Binti Zulkefli, Ahmad Razin Bin Zainal Abidin Faculty of Civil Engineering, Universiti Teknologi Malaysia,
More informationJournal of Engineering Research and Studies E-ISSN
Journal of Engineering Research and Studies E-ISS 0976-79 Research Article SPECTRAL SOLUTIO OF STEADY STATE CODUCTIO I ARBITRARY QUADRILATERAL DOMAIS Alavani Chitra R 1*, Joshi Pallavi A 1, S Pavitran
More informationAn Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems
An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems Long Chen University of California, Irvine chenlong@math.uci.edu Joint work with: Huayi Wei (Xiangtan University),
More informationIsogeometric Analysis (IGA) Part I (MS62) and Part II (MS78)
Isogeometric Analysis (IGA) Part I (MS62) and Part II (MS78) July 11, 4:00-6:00 PM Challenges in Isogeometric Analysis (IGA) Tor Dokken, Oslo, Norway Locally Refined B-splines, Tom Lyche, University of
More informationLab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.
More information