HONOM EUROPEAN WORKSHOP on High Order Nonlinear Numerical Methods for Evolutionary PDEs, Bordeaux, March 18-22, 2013
|
|
- Winfred Perkins
- 5 years ago
- Views:
Transcription
1 HONOM 2013 EUROPEAN WORKSHOP on High Order Nonlinear Numerical Methods for Evolutionary PDEs, Bordeaux, March 18-22, 2013 A priori-based mesh adaptation for third-order accurate Euler simulation Alexandre Carabias, Estelle Mbinky ( ) INRIA - Tropics project Sophia-Antipolis, France ( )INRIA - Gamma Project Rocquencourt, France Alexandre.Carabias@inria.fr
2 Motivation of the study Propagation of waves in large spaces. Mesh adaptation should follow the wave propagation. Adjoint-based adaptation criterion will concentrate on zone of interest for a chosen functional output. An accurate approximation scheme needs be chosen. 2 Carabias-Mbinky A priori mesh adaption for third-order...
3 Scope of the talk 1. Interpolation error analysis 2. Reconstruction error analysis 3. CENO2 Scheme 4. Approximation error analysis 5. Optimal metric 6. Resolution of optimum 7. Numerical experiments 8. Concluding remarks 3 Carabias-Mbinky A priori mesh adaption for third-order...
4 1. Interpolation error analysis (1) Π 2 : 2D Lagrange quadratic interpolation: (u Π 2 u)(δx) P u (δx). δx = (δx, δy) defines mesh size. The generic error model is a homogeneous polynomial of degree k = 3 in 2 variables : P u (x) = i=3 i=0 a(u) i are third derivatives of u. ( k i ) a(u) i x i y k i 4 Carabias-Mbinky A priori mesh adaption for third-order...
5 1. Interpolation error error analysis (2): Modeling of the quadratic order interpolation error Local optimization problem Main Idea: (P u ) 2 3 being homogeneous with a quadratic form with respect to mesh size, let us approach locally the variations of P u by a quadratic definite positive form taken at power 3 2, i.e, for (x, y) R2, P u (x, y) ( ) 3 t (x, y) Su loc 2 (x, y) Find the optimal local metric Su loc (optimal local directions and sizes) to construct a majorant of the local error model P u. The optimal local metric Su loc is the one whose unit ball is the maximum area ellipse included in the isoline 1 of P u. 5 Carabias-Mbinky A priori mesh adaption for third-order...
6 1. Interpolation error error analysis (3): Example P u = 6 x x 2 y 6.98 xy y 3 Optimal local metric in P u in ( t (x, y) S loc u (x, y)) 3 2 Superposition of the iso-lines of P u and ( t (x, y) S loc u (x, y)) Carabias-Mbinky A priori mesh adaption for third-order...
7 1. Interpolation error error analysis (4) Construction of optimal local metrics Binary decomposition by Sylvester s theorem [Brachat-Mourrain, 1999] P u (x, y) = Real case S loc u = t Q u λ λ (rank) j=1! Q u ««α1 β Q u = 1 ; α 2 β t α1 α Q u = 2 2 β 1 β 2 λ j (α j x + β j y) 3 ; β i = 1 In the complex case, λ 2 = λ 1 and α 2 = ᾱ 1. S loc u Complex case S loc u = tc Q u λ λ 2 2 3! «; t ᾱ1 c ᾱ Q u = 2 β 1 β 2 is real symmetric in both cases Q u 7 Carabias-Mbinky A priori mesh adaption for third-order...
8 1. Interpolation error error analysis (6) Optimal estimate Find the maximal ellipse inside the P u (x) = 1 level set. S u = c opt S loc u = t R u ( µ1 0 0 µ 2 ) R u c opt = 1 real case ; c opt = complex case. ( ) 3 P u (x, y) (x, y) S 2 u (x, y) 8 Carabias-Mbinky A priori mesh adaption for third-order...
9 1. Interpolation error error analysis (7) Mesh parametrization by a Riemannian metric For any x in Ω, M(x) is a symmetric matrix. M Mesh(M) Such that any edge ab of Mesh(M) has a unit M-length: 1 0 t ab M(a + tab) ab dt = 1. M(x) = d M (x)r(x)λ(x)r t (x) with det Λ = 1. The rotation R = (e ξ, e η ) parametrises the two orthogonal stretching directions of the metric. Λ is a diagonal matrix with eigenvalues λ 1 = m η /m ξ and λ 2 = m ξ /m η where m ξ and m η represent the two directional local mesh sizes in the characteristic/stretching directions of M. λ 1 /λ 2 is the stretching ratio. d M = m 1 ξ m 1 η is the node density. 9 Carabias-Mbinky A priori mesh adaption for third-order...
10 1. Interpolation error error analysis (8) Find the optimal Riemannian metric M defined on the computational domain for smallest weighted-l 1 error ) 3 min g(x) (trace(m 12 Su M 12 2 ) dx M Ω with d(m)dx = N M opt = N ( Ω ( µ1 ν 1 = µ 2 Ω ) ( 1 κ dx κ t ν1 0 R u 0 ν1 1 ) 1 2 κ = (2g µ 1 µ 2 ) 3 5. ) R u 10 Carabias-Mbinky A priori mesh adaption for third-order...
11 2. Reconstruction error error analysis Vertex, dual cell, 2-exact quadratic reconstruction Given ū i on on each cell i of centroid G i, find the c i,α, α 2 s.t. R2 0 ū i (x) = ū i + c i,α [(X G i ) α (X G i ) α dx] Cell i α 2 R2 0ū = R2 0 ū i,i dx = u i (R2 0ūi,j u j ) 2 = Min. Cell i j N(i) 11 Carabias-Mbinky A priori mesh adaption for third-order...
12 2. Reconstruction error analysis (2) Same treatment as interpolation R 2 : 2D quadratic ENO reconstruction: i=3 u(x) R2 0 u(x) ( k i )a i x i y k i i=0 ( t (x, y) S u (x, y) ) 3 2 a(u) i are third derivatives of u. 12 Carabias-Mbinky A priori mesh adaption for third-order...
13 3. CENO2 Scheme (1) Variational statement of a model problem B(u, v) = v F(u) dω Ω Γ vf Γ (u) dγ, B is linear with respect to v. Find u V such that B(u, v) = 0 v V. CENO discrete statement (after C. Groth), vertex version V 0 = {v 0, V 0 Celli = const i vertex} Find u 0 V 0 such that B(R 0 2 u 0, v 0 ) = 0 v 0 V 0 13 Carabias-Mbinky A priori mesh adaption for third-order...
14 3. CENO2 Scheme (2) 2-exact flux integration The integral on a cell interface C ij = C i C j is split into the integrals on the two segments of C ij. On each segment C (1) ij and C (2) ij a numerical integration with two Gauss points (two Riemann solvers) is applied. 14 Carabias-Mbinky A priori mesh adaption for third-order...
15 4. A priori error analysis (1) Cf. A posteriori analysis for a large set of p-order reconstruction-based Godunov methods Barth-Larson Scalar output j(u) = (g, u). Projection π 0 : V V 0, v π 0 v C i, dual cell, π 0 v Ci = vdx/meas(c i ). C i Output error δj = (g, R 0 2 π 0 u R 0 2 u 0 ). The adjoint state u 0 V 0 is the solution of: B u (R0 2 u 0 )(R 0 2 v 0, u 0) = (g, R 0 2 v 0 ), v 0 V Carabias-Mbinky A priori mesh adaption for third-order...
16 4. A priori error analysis (2) Simpler case: B(u, v) = F (v) where B is bilinear, F is linear, for example: B(u, v) = Ω vdiv(vu)dω + Γ uvv ndγ and F (v) = Γ u BvV ndγ. B(v, u0 ) = (g, v) (discr.adj. eq.) B(R2 0u 0, v) = F (v) (discr.state eq.) B(u, v) = F (v) (cont.state eq.) (g, R2 0π 0u R2 0u 0) = B(R2 0π 0u R2 0u 0, u0 ) (discr.adj. eq.) B(R2 0π 0u, u0 ) B(R0 2 u 0, u0 ) B(R2 0π 0u, u0 ) F (u 0 ) (discr.state eq.) B(R2 0π 0u, u0 ) B(u, u 0 ) (cont.state eq.) B(R2 0π 0u u, u0 ) The error is directly expressed in terms of the reconstruction error for exact solution. 16 Carabias-Mbinky A priori mesh adaption for third-order...
17 4. A priori error analysis (3) Unsteady Euler: W = (ρ, ρu, ρv, ρe) For the case of Euler eqs, we get after some calculations: B(R2 0 π 0 W W, W0 ) 2 K q (W, W ) R2 0 π 0 u q u q dω Ω q with (u q ) q=1,8 = (W, W t ), and in which the K q (W, W ) are built from space derivatives of W and Euler fluxes derivatives with respect to dependant variable W. 17 Carabias-Mbinky A priori mesh adaption for third-order...
18 5. Optimal metric (1) Using the previous reconstruction error estimate: R 0 2 π 0u q u q ( trace(m 1/2 S q M 1/2 ) ) 3 2. Optimization problem After the a priori analysis, we have to minimise the following error: = E = d 3 2 M q ( ) 3 K q (W, W ) trace(m 1/2 S q M 1/2 2 ) = ( trace(m 1/2 SM 1/2 )) 3 2 (trace[(r M Λ M R T M) 1 2 S (RM Λ M R T M) 1 2 ] ) 3 2 dxdy dxdy dxdy with constraint d M dxdy = N. 18 Carabias-Mbinky A priori mesh adaption for third-order...
19 5. Optimal metric (2) d 3 2 M Minimize E(M) = ) (trace[(r M Λ M R T M) 1 2 S (RM Λ M R T M) ] dxdy with constraint d M dxdy = N. Optimal solution: x :S = R S ΛS R T S R M opt = R S, Λ Mopt = det(s) 3 2 (d Mopt ) 5 2 = const. d Mopt = N ( det(s) 3 5 dxdy) 1 det(s) Λ S /det(s) M opt = d Mopt R Mopt Λ Mopt R T M opt 19 Carabias-Mbinky A priori mesh adaption for third-order...
20 5. Optimal metric (3): isotropic case Minimize E(M) = d 3 2 M (trace[(λ M ) 1 2 S (λm ) 1 2 ] ) 3 2 with constraint dxdy d M dxdy = N. Optimal solution: λ Mopt = trace[s] 1 /det(s) det(s) 3 2 (d Mopt ) 5 2 = const. d Mopt = N ( det(s) 3 5 dxdy) 1 det(s) 3 5 M opt = d Mopt λ Mopt 20 Carabias-Mbinky A priori mesh adaption for third-order...
21 6. Resolution of optimum, unsteady case The continuous model giving the adapted mesh involves a state system, an adjoint system and the optimality relation giving M opt. We discretise it. The time discretisation of the metric is made of coarser time intervals. We solve it. by the Global Unsteady Fixed Point algorithm of Belme et al. JCP Carabias-Mbinky A priori mesh adaption for third-order...
22 6. Resolution of optimum, unsteady case 22 Carabias-Mbinky A priori mesh adaption for third-order...
23 7. Numerical experiments We present a first series of experiments where the propagation of an acoustical perturbation is followed by the mesh adaptator in order to minimise the mesh effort. 1 Noise wall problem 2 Diffraction problem 23 Carabias-Mbinky A priori mesh adaption for third-order...
24 7. Numerical experiments: Noise wall problem (1) Propagation of an acoustic wave from a source while observing the impact on the detector. Acoustic source : r = Ae Bln(2)[x2 +y 2 ] cos(2πf ). Functional: j(w ) = T 0 M 1 2 (p p air ) 2 dmdt. 24 Carabias-Mbinky A priori mesh adaption for third-order...
25 7. Numerical experiments: Noise wall problem (2) Density field evolving in time on uniform mesh (line 1), adapted one (middle line), with corresponding adapted meshes (last line): 25 Carabias-Mbinky A priori mesh adaption for third-order...
26 7. Numerical experiments: Noise wall problem (3) Scalar output comparaison : MUSCL (-) vs CENO ( ) schemes. 5e e-08./functional_evolution_CENO.dat u 1:2./functional_evolution_MUSCL.dat u 1:2 4e e-08 Scalar output 3e e-08 2e e-08 1e-08 5e Time 26 Carabias-Mbinky A priori mesh adaption for third-order...
27 7. Numerical experiments: Diffraction problem (1) The diffraction of an acoustic wave travelling from a source location to a detector situated in a small region under the step. Acoustic source : r = Ae Bln(2)[x2 +y 2 ] cos(2πf ). Functional: j(w ) = T 0 M 1 2 (p p air ) 2 dmdt. 27 Carabias-Mbinky A priori mesh adaption for third-order...
28 7. Numerical experiments: Diffraction problem (2) Density field evolving in time on uniform mesh (line 1), adapted one (middle line), with corresponding adapted meshes (last line): 28 Carabias-Mbinky A priori mesh adaption for third-order...
29 7. Numerical experiments: Diffraction problem (3) Our goal-oriented method ensures the accuracy of the functional output j(w ) = T 1 0 M 2 (p p air ) 2 dmdt. It is thus interesting to analyse the integrand k(t) = 1 M 2 (p p air ) 2 dm of j(w ) on the micro M for different sizes of uniform vs. adapted meshes. 29 Carabias-Mbinky A priori mesh adaption for third-order...
30 7. Numerical experiments: Diffraction problem (4) Functional time integrand calculation on different sizes of non-adapted meshes (28K, 40K, 68K) vs. adapted meshes (mean sizes: 2620, 4892, 6130). 2e e-10./Uniform40K.dat u 1:2./Uniform28K.dat u 1:2./Uniform68K.dat u 1:2 3.5e-10 3e-10./Adapt2620.dat u 1:2./Adapt4892.dat u 1:2./Adapt6130.dat u 1:2 1.6e e e-10 Functional Integrand 1.2e-10 1e-10 8e-11 Functional Integrand 2e e-10 1e-10 6e-11 4e-11 5e-11 2e Time Time The convergence order is found to be 0.6 for uniform meshes and 1.98 for the adapted one. 30 Carabias-Mbinky A priori mesh adaption for third-order...
31 7. Numerical experiments: 3rd order error vs 2nd order error The previous computations where performed with 2nd-order error, made of 2nd-order derivatives, used with the 3rd-order CENO2 scheme. In a preliminary use of 3rd-order estimate, relying on 3rd-order derivatives, we compare two calculations: - one with 3rd-order scheme and 2nd-order error, - the second with 3rd-order scheme and 3rd-order error. The test case is the translation by Euler of a density Gaussian from left to right in a rectangle, with a functional observation with a Gaussian weight at rectangle center. 31 Carabias-Mbinky A priori mesh adaption for third-order...
32 7. Numerical experiments: Use of 2nd order error The integrand of the functional is a function of time which culminate when the center of density Gaussian coincides with the center of Gaussian which is the functional weight (center of rectangle). Left, 350 equivalent nodes, right 1500 equivalent nodes. 32 Carabias-Mbinky A priori mesh adaption for third-order...
33 7. Numerical experiments: Use of 3rd order error Here now are the results with the new 3rd-order estimate. Left, 350 equivalent nodes, right 1500 equivalent nodes. 33 Carabias-Mbinky A priori mesh adaption for third-order...
34 8. Concluding remarks (1) A third-order (spatially) accurate goal-oriented mesh adaptation method has been built on the basis of: The extension of Hessian analysis to higher order interpolation of Mbinky et al.. A novel a priori analysis. The Unsteady Global Fixed-Point mesh adaptation algorithm of Belme et al.. Numerical experiments are yet only 2D and at the very beginning. 34 Carabias-Mbinky A priori mesh adaption for third-order...
35 8. Concluding remarks (2) Next studies will address: Steady test cases. Flows with singularities. 3D extension. 35 Carabias-Mbinky A priori mesh adaption for third-order...
36 36 Carabias-Mbinky A priori mesh adaption for third-order...
Third 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 informationDiscontinuous Fluctuation Distribution
Discontinuous Fluctuation Distribution Matthew Hubbard School of Computing, University of Leeds, Leeds, LS2 9JT, UK. Abstract This paper describes a new numerical scheme for the approximation of steady
More informationDDFV Schemes for the Euler Equations
DDFV Schemes for the Euler Equations Christophe Berthon, Yves Coudière, Vivien Desveaux Journées du GDR Calcul, 5 juillet 2011 Introduction Hyperbolic system of conservation laws in 2D t W + x f(w)+ y
More informationA few multilinear algebraic definitions I Inner product A, B = i,j,k a ijk b ijk Frobenius norm Contracted tensor products A F = A, A C = A, B 2,3 = A
Krylov-type methods and perturbation analysis Berkant Savas Department of Mathematics Linköping University Workshop on Tensor Approximation in High Dimension Hausdorff Institute for Mathematics, Universität
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 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 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 informationAirfoil Design Optimization Using Reduced Order Models Based on Proper Orthogonal Decomposition
Airfoil Design Optimization Using Reduced Order Models Based on Proper Orthogonal Decomposition.5.5.5.5.5.5.5..5.95.9.85.8.75.7 Patrick A. LeGresley and Juan J. Alonso Dept. of Aeronautics & Astronautics
More informationHierarchical Reconstruction for Spectral Volume Method on Unstructured Grids
Hierarchical Reconstruction for Spectral Volume Method on Unstructured Grids Zhiliang Xu, Yingjie Liu and Chi-Wang Shu April 4, 2009 Abstract The hierarchical reconstruction (HR) [, 24] is applied to a
More informationOptimality Bounds for a Variational Relaxation of the Image Partitioning Problem
Optimality Bounds for a Variational Relaxation of the Image Partitioning Problem Jan Lellmann, Frank Lenzen, Christoph Schnörr Image and Pattern Analysis Group Universität Heidelberg EMMCVPR 2011 St. Petersburg,
More informationSOLVING 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 informationAdaptive Mesh Refinement Based on a Posteriori Error Estimation. Martin Juhas
Adaptive Mesh Refinement Based on a Posteriori Error Estimation by Martin Juhas A thesis submitted in conformity with the requirements for the degree of Master of Applied Science Graduate Department of
More informationAn L L p mesh-adaptive method for computing unsteady bi-fluid flows
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2010; 84:1376 1406 Published online 15 June 2010 in Wiley Online Library (wileyonlinelibrary.com)..2954 An L L p mesh-adaptive
More informationOn the order of accuracy and numerical performance of two classes of finite volume WENO schemes
On the order of accuracy and numerical performance of two classes of finite volume WENO schemes Rui Zhang, Mengping Zhang and Chi-Wang Shu November 29, 29 Abstract In this paper we consider two commonly
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 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 information3-D Wind Field Simulation over Complex Terrain
3-D Wind Field Simulation over Complex Terrain University Institute for Intelligent Systems and Numerical Applications in Engineering Congreso de la RSME 2015 Soluciones Matemáticas e Innovación en la
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 informationA third order conservative Lagrangian type scheme on curvilinear meshes for the compressible Euler equations
A third order conservative Lagrangian type scheme on curvilinear meshes for the compressible Euler equations Juan Cheng and Chi-Wang Shu 2 Abstract Based on the high order essentially non-oscillatory (ENO)
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 informationParameterization II Some slides from the Mesh Parameterization Course from Siggraph Asia
Parameterization II Some slides from the Mesh Parameterization Course from Siggraph Asia 2008 1 Non-Convex Non Convex Boundary Convex boundary creates significant distortion Free boundary is better 2 Fixed
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationweighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces.
weighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces. joint work with (S. Osher, R. Fedkiw and M. Kang) Desired properties for surface reconstruction:
More informationLecture 1: Finite Volume WENO Schemes Chi-Wang Shu
Lecture 1: Finite Volume WENO Schemes Chi-Wang Shu Division of Applied Mathematics Brown University Outline of the First Lecture General description of finite volume schemes for conservation laws The WENO
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
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 informationOn the Construction, Comparison, and Local Characteristic Decomposition for High-Order Central WENO Schemes
Journal of Computational Physics 8, 87 09 (00) doi:0.006/jcph.00.79 On the Construction, Comparison, and Local Characteristic Decomposition for High-Order Central WENO Schemes Jianxian Qiu, and Chi-Wang
More informationUnstructured 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 informationGeneralized trace ratio optimization and applications
Generalized trace ratio optimization and applications Mohammed Bellalij, Saïd Hanafi, Rita Macedo and Raca Todosijevic University of Valenciennes, France PGMO Days, 2-4 October 2013 ENSTA ParisTech PGMO
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 informationIsogeometric analysis with Axel
Isogeometric analysis with Axel Gang Xu, Régis Duvigneau,Bernard Mourrain INRIA Sophia-Antipolis gxu@sophia.inria.fr SAGA Workshop, March 18th, 2010 Outline 1 Isogeometric toolbox in EXCITING 2 Isogeometric
More informationDiffusion Wavelets for Natural Image Analysis
Diffusion Wavelets for Natural Image Analysis Tyrus Berry December 16, 2011 Contents 1 Project Description 2 2 Introduction to Diffusion Wavelets 2 2.1 Diffusion Multiresolution............................
More informationSemi-Conservative Schemes for Conservation Laws
Semi-Conservative Schemes for Conservation Laws Rosa Maria Pidatella 1 Gabriella Puppo, 2 Giovanni Russo, 1 Pietro Santagati 1 1 Università degli Studi di Catania, Catania, Italy 2 Università dell Insubria,
More informationGEOMETRICAL CONSTRAINTS IN THE LEVEL SET METHOD FOR SHAPE AND TOPOLOGY OPTIMIZATION
1 GEOMETRICAL CONSTRAINTS IN THE LEVEL SET METHOD FOR SHAPE AND TOPOLOGY OPTIMIZATION Grégoire ALLAIRE CMAP, Ecole Polytechnique Results obtained in collaboration with F. Jouve (LJLL, Paris 7), G. Michailidis
More informationCS 565 Computer Vision. Nazar Khan PUCIT Lectures 15 and 16: Optic Flow
CS 565 Computer Vision Nazar Khan PUCIT Lectures 15 and 16: Optic Flow Introduction Basic Problem given: image sequence f(x, y, z), where (x, y) specifies the location and z denotes time wanted: displacement
More informationTVD and ENO Schemes for Multidimensional Steady and Unsteady Flows A Comparative Analysis
TVD and ENO Schemes for Multidimensional Steady and Unsteady Flows A Comparative Analysis Philippe Angot IRPHE équipe MNM, UMR CNRS 38, Universit Aix-Marseille I. & II., 2, Av. Général Leclerc, 30 03 Marseille,
More informationImplicit Matrix Representations of Rational Bézier Curves and Surfaces
Implicit Matrix Representations of Rational Bézier Curves and Surfaces Laurent Busé INRIA Sophia Antipolis, France Laurent.Buse@inria.fr GD/SPM Conference, Denver, USA November 11, 2013 Overall motivation
More informationInterpolating and approximating scattered 3D-data with hierarchical tensor product B-splines
Interpolating and approximating scattered 3D-data with hierarchical tensor product B-splines Günther Greiner Kai Hormann Abstract In this note we describe surface reconstruction algorithms based on optimization
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 informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-4: Constrained optimization Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428 June
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 informationWhen Sparsity Meets Low-Rankness: Transform Learning With Non-Local Low-Rank Constraint For Image Restoration
When Sparsity Meets Low-Rankness: Transform Learning With Non-Local Low-Rank Constraint For Image Restoration Bihan Wen, Yanjun Li and Yoram Bresler Department of Electrical and Computer Engineering Coordinated
More informationNON-OSCILLATORY HIERARCHICAL RECONSTRUCTION FOR CENTRAL AND FINITE VOLUME SCHEMES
NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION FOR CENTRAL AND FINITE VOLUME SCHEMES YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. This is the continuation of the paper central discontinuous
More informationClustering: Classic Methods and Modern Views
Clustering: Classic Methods and Modern Views Marina Meilă University of Washington mmp@stat.washington.edu June 22, 2015 Lorentz Center Workshop on Clusters, Games and Axioms Outline Paradigms for clustering
More informationShape Modeling and Geometry Processing
252-0538-00L, Spring 2018 Shape Modeling and Geometry Processing Discrete Differential Geometry Differential Geometry Motivation Formalize geometric properties of shapes Roi Poranne # 2 Differential Geometry
More informationRobust and Efficient Adaptive Moving Mesh Solution of the 2-D Euler equations
Contemporary Mathematics Robust and Efficient Adaptive Moving Mesh Solution of the -D Euler equations P. A. Zegeling, W. D. de Boer, and H. Z. Tang Abstract. In this paper we describe an adaptive moving
More information04 - Normal Estimation, Curves
04 - Normal Estimation, Curves Acknowledgements: Olga Sorkine-Hornung Normal Estimation Implicit Surface Reconstruction Implicit function from point clouds Need consistently oriented normals < 0 0 > 0
More informationFinite Elements in Analysis and Design
Finite Elements in Analysis and Design 46 (2010) 181 -- 202 Contents lists available at ScienceDirect Finite Elements in Analysis and Design journal homepage: www. elsevier. com/ locate/ finel Size gradation
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationNon-Oscillatory Third Order Fluctuation Splitting Schemes for Steady Scalar Conservation Laws
Non-Oscillatory Third Order Fluctuation Splitting Schemes for Steady Scalar Conservation Laws Matthew Hubbard School of Computing, University of Leeds, Leeds, LS 9JT, UK. Abstract This paper addresses
More informationJournal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp
Journal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp. 59 74 SOLID MECHANICS WAVE PROPAGATION DUE TO AN EMBEDDED SEISMIC SOURCE IN A GRADED HALF-PLANE WITH RELIEF PECULIARITIES.
More informationNormalized cuts and image segmentation
Normalized cuts and image segmentation Department of EE University of Washington Yeping Su Xiaodan Song Normalized Cuts and Image Segmentation, IEEE Trans. PAMI, August 2000 5/20/2003 1 Outline 1. Image
More informationA cell-centered Arbitrary Lagrangian Eulerian method for two-dimensional multimaterial problems
A cell-centered Arbitrary Lagrangian Eulerian method for two-dimensional multimaterial problems Application to Inertial Confinement Fusion modeling in the direct drive context Pierre-Henri Maire, maire@celia.u-bordeaux1.fr
More informationAnalysis of the Adjoint Euler Equations as used for Gradient-based Aerodynamic Shape Optimization
Analysis of the Adjoint Euler Equations as used for Gradient-based Aerodynamic Shape Optimization Final Presentation Dylan Jude Graduate Research Assistant University of Maryland AMSC 663/664 May 4, 2017
More informationNon-Oscillatory Hierarchical Reconstruction for Central and Finite Volume Schemes
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol., No., pp. 933-963 Commun. Comput. Phys. October 7 Non-Oscillatory Hierarchical Reconstruction for Central and Finite Volume Schemes Yingjie Liu,, Chi-Wang Shu,
More informationManagement and control of the anisotropic unstructured meshes Applications in aerodynamics Institute of Aerodynamics and FLow Technology (DLR)
] Management and control of the anisotropic unstructured meshes in aerodynamics Institute of Aerodynamics and FLow Technology (DLR) Youssef Mesri INRIA, project-team Smash, 2004 Route des Lucioles, BP
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 informationInverse and Implicit functions
CHAPTER 3 Inverse and Implicit functions. Inverse Functions and Coordinate Changes Let U R d be a domain. Theorem. (Inverse function theorem). If ϕ : U R d is differentiable at a and Dϕ a is invertible,
More informationModel Reduction for Variable-Fidelity Optimization Frameworks
Model Reduction for Variable-Fidelity Optimization Frameworks Karen Willcox Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics Massachusetts Institute of Technology Workshop
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 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 informationStereo Matching for Calibrated Cameras without Correspondence
Stereo Matching for Calibrated Cameras without Correspondence U. Helmke, K. Hüper, and L. Vences Department of Mathematics, University of Würzburg, 97074 Würzburg, Germany helmke@mathematik.uni-wuerzburg.de
More informationDebojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
More informationHigh-Order CENO Reconstruction Scheme For Three-Dimensional Unstructured Mesh
High-Order CENO Reconstruction Scheme For Three-Dimensional Unstructured Mesh by Al-Amin Aziz A thesis submitted in conformity with the requirements for the degree of Masters of Applied Science Graduate
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 informationAnisotropic 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 informationFast sweeping methods for eikonal equations on triangular meshes
Fast sweeping methods for eikonal equations on triangular meshes Jianliang Qian, Yong-Tao Zhang 2, and Hong-Kai Zhao 3 Abstract The original fast sweeping method, which is an efficient iterative method
More informationA MULTI-DOMAIN ALE ALGORITHM FOR SIMULATING FLOWS INSIDE FREE-PISTON DRIVEN HYPERSONIC TEST FACILITIES
A MULTI-DOMAIN ALE ALGORITHM FOR SIMULATING FLOWS INSIDE FREE-PISTON DRIVEN HYPERSONIC TEST FACILITIES Khalil Bensassi, and Herman Deconinck Von Karman Institute for Fluid Dynamics Aeronautics & Aerospace
More informationHigh-order, conservative, finite difference schemes for computational MHD
High-order, conservative, finite difference schemes for computational MHD A. Mignone 1, P. Tzeferacos 1 and G. Bodo 2 [1] Dipartimento di Fisica Generale, Turin University, ITALY [2] INAF Astronomic Observatory
More informationSubdivision-stabilised immersed b-spline finite elements for moving boundary flows
Subdivision-stabilised immersed b-spline finite elements for moving boundary flows T. Rüberg, F. Cirak Department of Engineering, University of Cambridge, Trumpington Street, Cambridge CB2 1PZ, U.K. Abstract
More informationSPH: Why and what for?
SPH: Why and what for? 4 th SPHERIC training day David Le Touzé, Fluid Mechanics Laboratory, Ecole Centrale de Nantes / CNRS SPH What for and why? How it works? Why not for everything? Duality of SPH SPH
More information55:148 Digital Image Processing Chapter 11 3D Vision, Geometry
55:148 Digital Image Processing Chapter 11 3D Vision, Geometry Topics: Basics of projective geometry Points and hyperplanes in projective space Homography Estimating homography from point correspondence
More informationCALIBRATION BETWEEN DEPTH AND COLOR SENSORS FOR COMMODITY DEPTH CAMERAS. Cha Zhang and Zhengyou Zhang
CALIBRATION BETWEEN DEPTH AND COLOR SENSORS FOR COMMODITY DEPTH CAMERAS Cha Zhang and Zhengyou Zhang Communication and Collaboration Systems Group, Microsoft Research {chazhang, zhang}@microsoft.com ABSTRACT
More informationThe Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
More informationStability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow
Stability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow F. Haider 1, B. Courbet 1, J.P. Croisille 2 1 Département de Simulation Numérique des Ecoulements
More informationCOMPARISON OF FINITE VOLUME HIGH-ORDER SCHEMES FOR THE TWO-DIMENSIONAL EULER EQUATIONS
ECCOMAS Congress 206 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 0 June 206
More informationReconstruction of dielectrics from experimental data via a hybrid globally convergent/adaptive inverse algorithm
Reconstruction of dielectrics from experimental data via a hybrid globally convergent/adaptive inverse algorithm Larisa Beilina Michael V. Klibanov August 17, 21 Abstract The validity of a synthesis of
More informationsmooth coefficients H. Köstler, U. Rüde
A robust multigrid solver for the optical flow problem with non- smooth coefficients H. Köstler, U. Rüde Overview Optical Flow Problem Data term and various regularizers A Robust Multigrid Solver Galerkin
More informationOn the Order of Accuracy and Numerical Performance of Two Classes of Finite Volume WENO Schemes
Commun. Comput. Phys. doi: 8/cicp.99.84s Vol. 9, No. 3, pp. 87-87 March On the Order of Accuracy and Numerical Performance of Two Classes of Finite Volume WENO Schemes Rui Zhang, Mengping Zhang and Chi-Wang
More informationEnclosures of Roundoff Errors using SDP
Enclosures of Roundoff Errors using SDP Victor Magron, CNRS Jointly Certified Upper Bounds with G. Constantinides and A. Donaldson Metalibm workshop: Elementary functions, digital filters and beyond 12-13
More informationA Piecewise Linear Finite Element Discretization of the Diffusion Equation for Arbitrary Polyhedral Grids
UCRL-PROC-213665 A Piecewise Linear Finite Element Discretization of the Diffusion Equation for Arbitrary Polyhedral Grids T. S. Bailey, M. L. Adams, B. Yang, M. R. Zika July 18, 2005 Mathematics and Computation,
More informationDense Image-based Motion Estimation Algorithms & Optical Flow
Dense mage-based Motion Estimation Algorithms & Optical Flow Video A video is a sequence of frames captured at different times The video data is a function of v time (t) v space (x,y) ntroduction to motion
More informationThree-Dimensional Sensors Lecture 6: Point-Cloud Registration
Three-Dimensional Sensors Lecture 6: Point-Cloud Registration Radu Horaud INRIA Grenoble Rhone-Alpes, France Radu.Horaud@inria.fr http://perception.inrialpes.fr/ Point-Cloud Registration Methods Fuse data
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 informationA High-Order Accurate Unstructured GMRES Solver for Poisson s Equation
A High-Order Accurate Unstructured GMRES Solver for Poisson s Equation Amir Nejat * and Carl Ollivier-Gooch Department of Mechanical Engineering, The University of British Columbia, BC V6T 1Z4, Canada
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 informationBackward facing step Homework. Department of Fluid Mechanics. For Personal Use. Budapest University of Technology and Economics. Budapest, 2010 autumn
Backward facing step Homework Department of Fluid Mechanics Budapest University of Technology and Economics Budapest, 2010 autumn Updated: October 26, 2010 CONTENTS i Contents 1 Introduction 1 2 The problem
More informationAlgorithms for arithmetic Kleinian groups
supervised by John Voight September 6, 2011 Université Bordeaux 1 Definition A quaternion algebra over a field F is a central simple algebra of dimension 4. Explicitely, if char F 2 it admits a presentation
More informationfunctional prediction
An anisotropic hp-adaptation framework for functional prediction Marco Ceze 1 and Krzysztof J. Fidkowski 2 Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 48109 This paper presents
More informationFault Tolerant Domain Decomposition for Parabolic Problems
Fault Tolerant Domain Decomposition for Parabolic Problems Marc Garbey and Hatem Ltaief Department of Computer Science, University of Houston, Houston, TX 77204 USA garbey@cs.uh.edu, ltaief@cs.uh.edu 1
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 informationA Fast Phase Space Method for Computing Creeping Rays KTH. High Frequency Wave Propagation CSCAMM, September 19-22, 2005
A Fast Phase Space Method for Computing Creeping Rays Olof Runborg Mohammad Motamed KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Geometrical optics for high-frequency waves Consider
More informationLecture 8 Object Descriptors
Lecture 8 Object Descriptors Azadeh Fakhrzadeh Centre for Image Analysis Swedish University of Agricultural Sciences Uppsala University 2 Reading instructions Chapter 11.1 11.4 in G-W Azadeh Fakhrzadeh
More informationMotion Tracking and Event Understanding in Video Sequences
Motion Tracking and Event Understanding in Video Sequences Isaac Cohen Elaine Kang, Jinman Kang Institute for Robotics and Intelligent Systems University of Southern California Los Angeles, CA Objectives!
More informationRe-rendering from a Dense/Sparse Set of Images
Re-rendering from a Dense/Sparse Set of Images Ko Nishino Institute of Industrial Science The Univ. of Tokyo (Japan Science and Technology) kon@cvl.iis.u-tokyo.ac.jp Virtual/Augmented/Mixed Reality Three
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 informationKernels and representation
Kernels and representation Corso di AA, anno 2017/18, Padova Fabio Aiolli 20 Dicembre 2017 Fabio Aiolli Kernels and representation 20 Dicembre 2017 1 / 19 (Hierarchical) Representation Learning Hierarchical
More informationNumerical Methods for PDEs. SSC Workgroup Meetings Juan J. Alonso October 8, SSC Working Group Meetings, JJA 1
Numerical Methods for PDEs SSC Workgroup Meetings Juan J. Alonso October 8, 2001 SSC Working Group Meetings, JJA 1 Overview These notes are meant to be an overview of the various memory access patterns
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 informationMulti-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