Model Reduction Methods
|
|
- Tiffany Crawford
- 5 years ago
- Views:
Transcription
1 Martin A. Grepl 1, Gianluigi Rozza 2 1 IGPM, RWTH Aachen University, Germany 2 MATHICSE - CMCS, Ecole Polytechnique Fédérale de Lausanne, Switzerland Summer School "Optimal Control of PDEs" Cortona (Italy), July 12-17, 21
2 rbmit: - Reduced Basis Software Library in Matlab(C) Environment developed at MIT by DBP Huynh, CN Nguyen, G Rozza and AT Patera based on extensive use of Matlab ToolBoxes like Symbolic, PDEs, Optimization The user must describe the problem. The input can be separated into three parts: The User Input geometry: Ω o(µ) is defined by providing points coordinates, straight/curvy edges describing all regions Ω k o(µ) material properties: coefficients are provided for differential operator in each region Ω k o(µ) and for boundary conditions. parameter control and settings: parameter domain D, reference parameters and other RB information (e.g. N max)
3 rbmit Users Interface Geometry Parameters PDE and BC Output
4 Problem Formulation, Offline and Online Steps The Problem Formulation Step Domain Decomposition/geometric transformations: coefficients Θ q(µ) are generated for each sub-domain (coupled with material µ-properties) FE mesh is generated, discrete FE stiffness matrices/vectors are assembled for each sub-domain to form the µ-independent components The RB Offline Step By a greedy algorithm the RB parameter sample set is obtained FE/RB matrices are saved in order to be used by the Online Step The RB Online Step Given µ D, the RB Online Evaluator returns output and error bound Online_RB (probname, µ, outputname,...): µ s N N (µ), s N (µ) The RB Visualizer renders the field variable(s) and provides error bounds Vis_RB (probname, µ): µ Ω, u N N (x; µ) for all x in Ω o(µ)
5 Problem Formulation, Offline and Online Steps The Problem Formulation Step Domain Decomposition/geometric transformations: coefficients Θ q(µ) are generated for each sub-domain (coupled with material µ-properties) FE mesh is generated, discrete FE stiffness matrices/vectors are assembled for each sub-domain to form the µ-independent components The RB Offline Step By a greedy algorithm the RB parameter sample set is obtained FE/RB matrices are saved in order to be used by the Online Step The RB Online Step Given µ D, the RB Online Evaluator returns output and error bound Online_RB (probname, µ, outputname,...): µ s N N (µ), s N (µ) The RB Visualizer renders the field variable(s) and provides error bounds Vis_RB (probname, µ): µ Ω, u N N (x; µ) for all x in Ω o(µ)
6 Problem Formulation, Offline and Online Steps The Problem Formulation Step Domain Decomposition/geometric transformations: coefficients Θ q(µ) are generated for each sub-domain (coupled with material µ-properties) FE mesh is generated, discrete FE stiffness matrices/vectors are assembled for each sub-domain to form the µ-independent components The RB Offline Step By a greedy algorithm the RB parameter sample set is obtained FE/RB matrices are saved in order to be used by the Online Step The RB Online Step Given µ D, the RB Online Evaluator returns output and error bound Online_RB (probname, µ, outputname,...): µ s N N (µ), s N (µ) The RB Visualizer renders the field variable(s) and provides error bounds Vis_RB (probname, µ): µ Ω, u N N (x; µ) for all x in Ω o(µ)
7 rbmit Example: the Thermal Fin problem Engineering aspects Heat sink designed for cooling of high-density electronic components Shaded domain due to assumed periodicity and symmetry (multi-fin sink) Flowing air is modelled though a simple convection HT coefficient: to compute temperature at the base of the spreader Physical and geometrical parametrization µ 1 = Bi = h c dper/ κ fin Biot number µ 1 [.1,.5] µ 2 = L = L/ d per nondimensional fin height µ 2 [2, 8] µ 3 = κ = κ sp/ κ fin spreader/fin conductivity ratio µ 3 [1, 1]
8 rbmit Example: the Thermal Fin problem Modeling: temperature u o(µ) over Ω o(µ) satisfies the steady heat equation Output: average temperature over the base of the spreader (component to be cooled, being the hottest location in the system)» «µ3 u x o i µ 3 x o j o(µ) = in Ω 1 o {z } κ 1 o ij» «1 u x o i 1 x o j o(µ) = in Ω 2 o (µ 2) {z }, µ2, µ , µ2 3 2, µ2 2 κ 2 o ij n o i κ 1 u o o ij (µ) = 1 on Γ x o 1 o j n o i κ 2 u o o ij (µ) + (µ x 1 )u o = on Γ R = Γ o 5 Γ 6 o j n o i κ o ij u x o j o(µ) = on Γ \ (Γ o 1 Γ R ) Z output: T o av(µ) = 2 u o(µ) Γ o 1, 3 5 (, ) 3 2, , ,
9 rbmit Example: the Thermal Fin problem Example of geometry and field visualizations provided by rbmit package Example of initial geometry, domain decompostion, FE mesh and RB solution visualization for a thermal fin problem
10 rbmit Example: the Thermal Fin problem Approximation property # of mesh nodes N 4198 # of RB functions N 1 Reduced Basis vs Finite Elements RB online.13s (N = 7) evaluation time.15s (N = 13) FEM sol. µ s N (µ) 1.96s Reduction of 4:1 in linear system dimension Online evaluation 5 6% of the FEM computational cost μ1 RB output/error bars -[s N (µ) s N (µ), s N (µ) + s N (µ)] as a function of µ 1 for µ 2 = 2, µ 3 = 1 and N = 6. RB temperature field for different choices of parameters: µ = (.5, 2, 1), µ = (.5, 2, 5), µ = (.1, 2, 1).
11 Geometrical Parametrization: the goal The parametrized (original) domain Ω o (µ) is the image of a fixed (reference) domain Ω through a map T ( ; µ) : Ω Ω o (µ) In order to recover the affine parameter dependence, the parametric map T ( ; µ) has to be an affine map. The rbmit software allows to deal with more complex configurations by means of automatically built affine mappings
12 Geometrical Parametrization: Domain Decomposition Domain decomposition: definition Original Domain Ω o (µ), u e o Xe o (Ω o(µ)) Ω o (µ) = K dom k=1 Ωk o (µ) ; Reference domain Ω, u e X e (Ω) Ω = K dom k=1 Ωk, common configuration where Ω = Ω o (µ ref ) for µ ref D. For Ω k, Ω k o (µ) we choose in R2 triangles, elliptical triangles and curvy triangles. In R 3 we choose parallelepipeds (and in theory tetrahedra). Connectivity requirement: subdomain intersections must be an entire edge, a vertex, or null.
13 Geometrical Parametrization: Domain Decomposition Domain decomposition: definition Original Domain Ω o (µ), u e o Xe o (Ω o(µ)) Ω o (µ) = K dom k=1 Ωk o (µ) ; Reference domain Ω, u e X e (Ω) Ω = K dom k=1 Ωk, common configuration where Ω = Ω o (µ ref ) for µ ref D. For Ω k, Ω k o (µ) we choose in R2 triangles, elliptical triangles and curvy triangles. In R 3 we choose parallelepipeds (and in theory tetrahedra). Connectivity requirement: subdomain intersections must be an entire edge, a vertex, or null.
14 Geometrical Parametrization: Affine Mappings Require where Ω k o (µ) = T aff,k (Ω k ; µ), 1 k K dom, T aff,k (x; µ) = C aff,k (µ) + G aff,k (µ)x, µ D is an invertible affine mapping from Ω k onto Ω k o (µ). Further require T aff,k (x; µ) = T aff,k (x; µ), x Ω k Ω k, 1 k, k K dom, µ D to ensure a continuous piecewise-affine global mapping T aff ( ; µ) from Ω onto Ω o (µ). It follows that for w o H 1 (Ω o(µ)), w o T aff = H 1 (Ω).
15 Geometrical Parametrization: Affine Mappings Require where Ω k o (µ) = T aff,k (Ω k ; µ), 1 k K dom, T aff,k (x; µ) = C aff,k (µ) + G aff,k (µ)x, µ D is an invertible affine mapping from Ω k onto Ω k o (µ). Further require T aff,k (x; µ) = T aff,k (x; µ), x Ω k Ω k, 1 k, k K dom, µ D to ensure a continuous piecewise-affine global mapping T aff ( ; µ) from Ω onto Ω o (µ). It follows that for w o H 1 (Ω o(µ)), w o T aff = H 1 (Ω).
16 Geometrical Parametrization: Affine Mappings Elliptical Triangles: definition Outwards: Inwards: { }} { O(µ) = (x cen o1, xcen o2 ( )T ) cos φ(µ) sin φ(µ) Q rot (µ) = sin φ(µ) cos φ(µ) S(µ) = diag(ρ 1 (µ), ρ 2 (µ))
17 Geometrical Parametrization: Affine Mappings Elliptical Triangles: constraints Given x 2 o (µ), x3 o (µ), find x1 o (µ), x4 o (µ) ( T aff,1&2 ) } (i) produce desired elliptical arc µ D; (ii) satisfy internal angle criterion these conditions ensure continuous invertible mappings. Explicit recipes for admissible x 1 o(µ) (Inwards case) and x 4 o(µ) (Outwards case) are readily obtained.
18 Geometrical Parametrization: Affine Mappings Elliptical Triangles: example (CinS triangulation) Ω o (µ) : µ = (µ 1, µ 2,...) D [.8, 1.2] 2...
19 Affine Mappings Elliptical Triangles: example (CinS triangulation) Ω = Ωo (µref = (1, 1)) Ωo (µ = (.8, 1.2)) 1 2
20 Geometrical Parametrization: Affine Mappings Curvy Triangles: definition Outwards: Inwards: { }} { O(µ) = (x cen o1, xcen o2 ( )T ) cos φ(µ) sin φ(µ) Q rot (µ) = sin φ(µ) cos φ(µ) S(µ) = diag(ρ 1 (µ), ρ 2 (µ))
21 Geometrical Parametrization: Affine Mappings Curvy Triangles: constraints Given x 2 o (µ), x3 o (µ), find x1 o (µ), x4 o (µ) ( T aff,1&2 ) } (i) produce desired curvy arc µ D; (ii) satisfy internal angle criterion these conditions ensure continuous invertible mappings. Quasi-explicit recipes for admissible x 1 o (µ) and x4 o (µ) can (sometimes) be obtained in the convex/concave case.
22 Geometrical Parametrization: Affine Mappings Elliptical Triangles: example (Cosine triangulation) (say) Ω o (µ) : µ = (µ 1,...) D [ 1 6, 1 2 ]...
23 Geometrical Parametrization: Affine Mappings Elliptical Triangles: example (Cosine triangulation) Ω = Ω o (µ ref = 1 3 ) Ω o(µ = 1 2 )
24 Geometrical Parametrization: Tensor Transformations Transformation: Formulation on original domain (IR 2 ) For w, v H 1 (Ω o (µ)) a o (w, v; µ) = Kdom k=1 Ω ko (µ) [ w x o1 w x o2 w u e o (µ) H1 (Ω o(µ)) ] K k oij (µ) v x o1 v x o2 v where K k o : D R3 3, SPD for 1 k K dom (note K k o affine in x o is also permissible). We consider the scalar case; the vector case (linear elasticity) admits an analogous treatment.
25 Geometrical Parametrization: Tensor Transformations Transformation: Formulation on reference domain For w, v H 1 (Ω) u e (µ) H 1 (Ω) a(w, v; µ) = Kdom k=1 Ω k [ w x 1 w x 2 w ] K k ij (µ) v x 1 v x 2 v K k (µ) = det G aff,k (µ) D(µ)K k o (µ)dt (µ), and D(µ) = (Gaff,k ) 1. 1
26 Geometrical Parametrization: Tensor Transformations Transformation: Affine form Expand a(w, v; µ) = K 1 11 }{{ (µ) w v +... } Ω 1 x 1 x }{{ 1 } Θ 1 (µ) a 1 (w,v) We can often greatly reduce the requisite Q. with as many as Q = 4K terms. Achtung! Many interesting problems are not affine (or require Q very large). For example, K k o (x; µ) for general x dependence; and nonzero Neumann conditions on curvy Ω yield non-affine a(, ; µ).
27 Geometrical Parametrization: Tensor Transformations Transformation: Affine form Expand a(w, v; µ) = K 1 11 }{{ (µ) w v +... } Ω 1 x 1 x }{{ 1 } Θ 1 (µ) a 1 (w,v) We can often greatly reduce the requisite Q. with as many as Q = 4K terms. Achtung! Many interesting problems are not affine (or require Q very large). For example, K k o (x; µ) for general x dependence; and nonzero Neumann conditions on curvy Ω yield non-affine a(, ; µ).
28 RB Simulation / Output Evaluation Geometries, Software: Geometrical Problems Flow Control and Optimal Design with RB Methods Basic parametrized geometries Automatic Geometrical Parametrization: Basic setup: AirfoilAirfoil of the NACAExample 4-digits family (symmetric c I in parallel flow Airfoil of the NACA 4-digits family (symmetric case) Ω L u µ1 x2 Thickness distribution 4-digits profile x2 = µ1 ( Ω x x1.126x1.352x x1.121x Airfoil geometry description for the rbmit software.1.2 x1 x2 «x1 x2 «= = 1 «« ±µ1 /2 «±µ1 /2 « t2.2969t.126t2.352t t6.121t8! 1 t2.2969t.126t2.352t t6.121t8!.8 1, t [,, t [.3, 1] Complexity (potential flow): # of subdomains K = 42/84 # of bilinear forms Q = 16/2.3]
29 RB Simulation / Output Evaluation Geometries, Software: Geometrical Problems Flow Control and Optimal Design with RB Methods Basic parametrized geometries Automatic Geometrical Parametrization: Basic setup: AirfoilAirfoil of the NACAExample 4-digits family (symmetric c I in parallel flow Airfoil of the NACA 4-digits family (symmetric case) Ω L u µ1 x2 Thickness distribution 4-digits profile x2 = µ1 ( Ω x x1.126x1.352x x1.121x Airfoil geometry description for the rbmit software.1.2 x1 x2 «x1 x2 «= = 1 «« ±µ1 /2 «±µ1 /2 « t2.2969t.126t2.352t t6.121t8! 1 t2.2969t.126t2.352t t6.121t8!.8 1, t [,, t [.3, 1] Complexity (potential flow): # of subdomains K = 42/84 # of bilinear forms Q = 16/2.3]
30 RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Potential Flows: Output Evaluation Parameters: thickness µ 1 [4, 24], angle of attack µ 2 [.1, π/16] and profiles distance (µ 3, µ 4) [.85, 1] [1.2, 1.6] Output: pressure coefficient c p = p p in 1 2 ρ u in 2 around the profiles
31 RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Thermal Flows: Output Evaluation Thermal Flows around an airfoil with ground effect included Physical and geometrical parametrization: thickness µ 1 [4, 24], ground distance µ 2 [1.5, 3], Peclet number µ 3 [1, 1] The thermal boundary layer on the ground and on profile becomes thinner and more separated for higher Peclet numbers
32 RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Potential/Thermal Flows: Computational Costs Potential Flow mesh nodes N 3693 subdomains 163 affine operator components Q a 12 affine rhs components Q f 3 reduced basis functions N 2 t online F E 17.15s t online RB.8s speedup 195 RB online time as funcion of N (potential flow) Thermal Flow mesh nodes N 6727 subdomains 88 affine operator components Q a 42 reduced basis functions N pr 7 reduced basis functions N du 61 t online F E 16.97s t online RB.26s speedup 65 Average outflow temperature w.r. to ground distance µ 2 (thermal flow)
33 RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Potential Flow: Shape Optimization Problem Airfoil inverse design problem Z 1 1/2 min p(s, µ) p target(s) ds«2 + λ ˆα(µ) 5 2, µ D Z Z s.t. u v dω o = fv dω o v H 1 (Ω o(µ)) Ω o(µ) u = on Γ out, Ω o(µ) u u = 1 on Γin, n Choose target airfoil (ex: NACA4412) and compute pressure distribution p target on its surface using the Bernoulli equation (p = p 1 2 u 2 ) Objective: find small perturbation of reference airfoil NACA12 s.t. pressure distribution on the airfoil surface is close to p target Add penalty term to enforce the constraint on the angle of attack (AOA = 5 ) n = elsewhere
34 RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Potential Flow: Shape Optimization Problem Geometrical parametrization: Free-Form Deformation Techniques Geometrical parameters µ 1,..., µ P are chosen as the perturbations of a (small) lattice of FFD control points
35 RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Potential Flow: Shape Optimization Problem Pressure distributions and computational cost (online solution of the parametric PDE) Finite element method Reduced basis method CPU time (seconds) Real time barrier (response time < 1 s) 1 1 (a) Inverse design (b) Target airfoil Number of mesh nodes x 1 4 (c) Computational costs Number of mesh nodes N 843 Lattice of FFD control points 6 4 Number of shape parameters 8 Number of reduced basis functions N 52 Error tolerance for RB greedy ε RB tol 1 4 Number of affine expansion terms Q a 8 Error tolerance for EIM greedy ε EIM tol Results from T. Lassila, G. Rozza, Comput. Methods Appl. Mech. Engrg. 199 (21) Reduction of 5:1 in parametric complexity compared to explicit nodal deformation Reduction of 2:1 in linear system dimension
36 Thermal Flows: Optimal Design of Airfoils RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Optimal heat exchange problem» u target 1 min µ D s.t. Z Ω o(µ) Z 2 u dγ + λ [α(µ) α ] 2, Γ out Γ out ε u v + v Z b u dω o = fv dω o u = on Γout, u = T on Γin Γfree, n u = T 1 on Γ surf, u = T 2 on airfoil Objective: find airfoil shape and vertical position s.t. average temperature over outflow equals u target and angle of attack equals α Heat exchange of an airfoil in exterior flow with b = [1; ] and ε =.2 is considered Penalty term enforces the constraint on the angle of attack (AOA = α ) Ω o(µ)
37 Thermal Flows: Optimal Design of Airfoils RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods α = 7, u target = 4.1 α = 5, u target = 4.5 Number of mesh nodes N Lattice of FFD control points 6 6 Number of shape parameters 8 Number of reduced basis functions N 36 Error tolerance for RB greedy ε RB tol 1 5 Number of affine expansion terms Q a 18 Error tolerance for EIM greedy ε EIM tol 1 4 Reduction of 1:1 in parametric complexity compared to explicit nodal deformation Reduction of 436:1 in linear system dimension Results from G. Rozza, T. Lassila, A. Manzoni, Proc. of Icosahom Conference, 29, in press
38 RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Environmental Flows: Control of Air Pollution Control of Air pollution min J(µ) = y(u(µ u)) z d µ D ZD 2 dω, o(µ g) Z s.t. (ν(µ p) y v + V(µ p) y) dω o = Ω o(µ g) y = on ΓN, n y = on ΓD Goal: to regulate the pollutant emission by industrial plants in order to keep the pollutant level below a fixed threshold over an observation area Parameters: control input µ u, which define the control function u = u(µ u); physical input µ p (e.g. viscosity, advection velocity); geometrical input µ g (domain configuration) Z Ω o(µ g) u(µ p)v dω o
39 RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Environmental Flows: Control of Air Pollution Control input µ u : variable emission rate Use of reduced basis method for the efficient solution of the state and the adjoint parametrized problems Results from A. Quarteroni, G. Rozza, A, Quaini, Advances in Numerical Mathematics, 27, p
40 RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Environmental Flows: Control of Air Pollution Physical Input µ p : wind direction Use of reduced basis method for the efficient solution of the state and the adjoint parametrized problems Results from A. Quarteroni, G. Rozza, A, Quaini, Advances in Numerical Mathematics, 27, p
41 RB Simulation / Output Evaluation Flow Control and Optimal Design with RB Methods Environmental Flows: Control of Air Pollution Geometrical Input µ g : parametrized domain Use of reduced basis method for the efficient solution of the state and the adjoint parametrized problems Results from A. Quarteroni, G. Rozza, A, Quaini, Advances in Numerical Mathematics, 27, p
Model Reduction Methods
Martin A. Grepl 1, Gianluigi Rozza 2 1 IGPM, RWTH Aachen University, Germany 2 MATHICSE - CMCS, Ecole Polytechnique Fédérale de Lausanne, Switzerland Summer School "Optimal Control of PDEs" Cortona (Italy),
More informationReduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling
Reduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling Toni Lassila, Andrea Manzoni, Gianluigi Rozza CMCS - MATHICSE - École Polytechnique Fédérale
More informationReduced basis approximation and a posteriori error estimation for Stokes flows in parametrized geometries: roles of the inf-sup stability constants
Reduced basis approximation and a posteriori error estimation for Stokes flows in parametrized geometries: roles of the inf-sup stability constants The MIT Faculty has made this article openly available.
More informationReduced basis method and domain decomposition for elliptic problems in networks and complex parametrized geometries
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 30.2015 November 2015 Reduced basis method and domain
More informationSome recent developments on ROMs in computational fluid dynamics
Some recent developments on ROMs in computational fluid dynamics F. Ballarin 1, S. Ali 1, E. Delgado 2, D. Torlo 1,3, T. Chacón 2, M. Gómez 2, G. Rozza 1 1 mathlab, Mathematics Area, SISSA, Trieste, Italy
More informationReduced order modelling in multi-parametrized system: applications to inverse problems and optimal control
Reduced order modelling in multi-parametrized system: applications to inverse problems and optimal control Gianluigi Rozza in collaboration with Federico Negri, Andrea Manzoni SISSA Mathlab EPFL - MATHICSE
More informationA Natural Norm Successive Constraint Method for Inf-Sup Lower Bounds
A Natural Norm Successive Constraint Method for Inf-Sup Lower Bounds D.B.P. Huynh a, D.J. Knezevic b, Y. Chen c, J.S. Hesthaven c, A.T. Patera,b a Mechanical Engineering Department, National University
More informationToni Lassila 1, Andrea Manzoni 1 and Gianluigi Rozza 1
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ON THE APPROIMATION OF STABILITY FACTORS FOR GENERAL PARAMETRIZED PARTIAL DIFFERENTIAL
More informationEfficient model reduction of parametrized systems by matrix discrete empirical interpolation
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 02.2015 February 2015 Efficient model reduction of
More informationUNITEXT La Matematica per il 3+2
UNITEXT La Matematica per il 3+2 Volume 92 Editor-in-chief A. Quarteroni Series editors L. Ambrosio P. Biscari C. Ciliberto M. Ledoux W.J. Runggaldier More information about this series at http://www.springer.com/series/5418
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationMATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 05.2013 February 2013 (NEW 10.2013 A weighted empirical
More informationModel reduction for multiscale problems
Model reduction for multiscale problems Mario Ohlberger Dec. 12-16, 2011 RICAM, Linz , > Outline Motivation: Multi-Scale and Multi-Physics Problems Model Reduction: The Reduced Basis Approach A new Reduced
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationThe Reduced Basis Method:
The Reduced Basis Method: Application to Problems in Fluid and Solid Mechanics ++ Aachen Institute for Advanced Study in Computational Engineering Science K. Veroy-Grepl 31 January 2014 AICES & MOR I The
More informationBounding Stability Constants for Affinely Parameter-Dependent Operators
Bounding Stability Constants for Affinely Parameter-Dependent Operators Robert O Connor a a RWTH Aachen University, Aachen, Germany Received *****; accepted after revision +++++ Presented by Abstract In
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationA Two-Step Certified Reduced Basis Method
A Two-Step Certified Reduced Basis Method The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published Publisher Eftang,
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationEmpirical Interpolation of Nonlinear Parametrized Evolution Operators
MÜNSTER Empirical Interpolation of Nonlinear Parametrized Evolution Operators Martin Drohmann, Bernard Haasdonk, Mario Ohlberger 03/12/2010 MÜNSTER 2/20 > Motivation: Reduced Basis Method RB Scenario:
More informationCertified Reduced Basis Approximation for Parametrized Partial Differential Equations and Applications
Certified Reduced Basis Approximation for Parametrized Partial Differential Equations and Applications Alfio Quarteroni,, Gianluigi Rozza and Andrea Manzoni Modelling and Scientific Computing, CMCS, Mathematics
More informationMultiobjective PDE-constrained optimization using the reduced basis method
Multiobjective PDE-constrained optimization using the reduced basis method Laura Iapichino1, Stefan Ulbrich2, Stefan Volkwein1 1 Universita t 2 Konstanz - Fachbereich Mathematik und Statistik Technischen
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationChapitre 1. Cours Reduced Basis Output Bounds Methodes du 1er Semestre
Chapitre 1 Cours Methodes du 1er Semestre 2007-2008 Pr. Christophe Prud'homme christophe.prudhomme@ujf-grenoble.fr Université Joseph Fourier Grenoble 1 1.1 1 2 FEM Approximation Approximation Estimation
More informationA SEAMLESS REDUCED BASIS ELEMENT METHOD FOR 2D MAXWELL S PROBLEM: AN INTRODUCTION
A SEAMLESS REDUCED BASIS ELEMENT METHOD FOR 2D MAXWELL S PROBLEM: AN INTRODUCTION YANLAI CHEN, JAN S. HESTHAVEN, AND YVON MADAY Abstract. We present a reduced basis element method (RBEM) for the time-harmonic
More informationReduced basis approximation and a posteriori error estimation for stress intensity factors
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2006; 00:1 6 [Version: 2002/09/18 v2.02] Reduced basis approximation and a posteriori error estimation for stress intensity
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationAn Extended Finite Element Method for a Two-Phase Stokes problem
XFEM project An Extended Finite Element Method for a Two-Phase Stokes problem P. Lederer, C. Pfeiler, C. Wintersteiger Advisor: Dr. C. Lehrenfeld August 5, 2015 Contents 1 Problem description 2 1.1 Physics.........................................
More informationReduced basis method for efficient model order reduction of optimal control problems
Reduced basis method for efficient model order reduction of optimal control problems Laura Iapichino 1 joint work with Giulia Fabrini 2 and Stefan Volkwein 2 1 Department of Mathematics and Computer Science,
More informationASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM C. Canuto 1, L. F. Pavarino 2, and A. B. Pieri 3 1 Introduction Discontinuous Galerkin (DG) methods for partial differential
More informationCertified reduced basis approximation for parametrized partial differential equations and applications
Journal of Mathematics in Industry (2011) 1:3 DOI 10.1186/2190-5983-1-3 RESEARCH OpenAccess Certified reduced basis approximation for parametrized partial differential equations and applications Alfio
More informationA CERTIFIED TRUST REGION REDUCED BASIS APPROACH TO PDE-CONSTRAINED OPTIMIZATION
A CERTIFIED TRUST REGIO REDUCED BASIS APPROACH TO PDE-COSTRAIED OPTIMIZATIO ELIZABETH QIA, MARTI GREPL, KARE VEROY, AD KARE WILLCOX Abstract. Parameter optimization problems constrained by partial differential
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationModel reduction of parametrized aerodynamic flows: stability, error control, and empirical quadrature. Masayuki Yano
Model reduction of parametrized aerodynamic flows: stability, error control, and empirical quadrature Masayuki Yano University of Toronto Institute for Aerospace Studies Acknowledgments: Anthony Patera;
More informationCertified Reduced Basis Methods for Parametrized Partial Differential Equations
Jan S Hesthaven, Gianluigi Rozza, Benjamin Stamm Certified Reduced Basis Methods for Parametrized Partial Differential Equations Monograph Springer Briefs in Mathematics June 23, 2015 Springer To our
More informationReduced Basis Method for Parametric
1/24 Reduced Basis Method for Parametric H 2 -Optimal Control Problems NUMOC2017 Andreas Schmidt, Bernard Haasdonk Institute for Applied Analysis and Numerical Simulation, University of Stuttgart andreas.schmidt@mathematik.uni-stuttgart.de
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationMixed Hybrid Finite Element Method: an introduction
Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola
More informationA reduced basis method with an exact solution certificate and spatio-parameter adaptivity: application to linear elasticity
A reduced basis method with an exact solution certificate and spatio-parameter adaptivity: application to linear elasticity Masayuki Yano Abstract We present a reduced basis method for parametrized linear
More informationNonlinear Model Reduction for Uncertainty Quantification in Large-Scale Inverse Problems
Nonlinear Model Reduction for Uncertainty Quantification in Large-Scale Inverse Problems Krzysztof Fidkowski, David Galbally*, Karen Willcox* (*MIT) Computational Aerospace Sciences Seminar Aerospace Engineering
More informationA certified trust region reduced basis approach to PDE-constrained optimization
J U E 2 0 1 6 P R E P R I T 4 5 3 A certified trust region reduced basis approach to PDE-constrained optimization Elizabeth Qian *, Martin Grepl, Karen Veroy and Karen Willcox * Institut für Geometrie
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationReduced basis approximation and a posteriori error bounds for 4D-Var data assimilation
M A R C H 2 0 1 7 P R E P R I T 4 6 3 Reduced basis approximation and a posteriori error bounds for 4D-Var data assimilation Mark Kärcher, Sébastien Boyaval, Martin A. Grepl and Karen Veroy Institut für
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationMath 302 Test 1 Review
Math Test Review. Given two points in R, x, y, z and x, y, z, show the point x + x, y + y, z + z is on the line between these two points and is the same distance from each of them. The line is rt x, y,
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationChapter 3. Formulation of FEM for Two-Dimensional Problems
Chapter Formulation of FEM for Two-Dimensional Problems.1 Two-Dimensional FEM Formulation Many details of 1D and 2D formulations are the same. To demonstrate how a 2D formulation works we ll use the following
More informationEfficient BDF time discretization of the Navier Stokes equations with VMS LES modeling in a High Performance Computing framework
Introduction VS LES/BDF Implementation Results Conclusions Swiss Numerical Analysis Day 2015, Geneva, Switzerland, 17 April 2015. Efficient BDF time discretization of the Navier Stokes equations with VS
More informationReduced Basis Approximation and A Posteriori Error. Estimation for the Time-Dependent Viscous Burgers Equation
Calcolo manuscript o. (will be inserted by the editor) Reduced Basis Approximation and A Posteriori Error Estimation for the Time-Dependent Viscous Burgers Equation goc-cuong guyen 1, Gianluigi Rozza 1,,
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationFinite Element Modeling of Electromagnetic Systems
Finite Element Modeling of Electromagnetic Systems Mathematical and numerical tools Unit of Applied and Computational Electromagnetics (ACE) Dept. of Electrical Engineering - University of Liège - Belgium
More informationBake, shake or break - and other applications for the FEM. 5: Do real-life experimentation using your FEM code
Bake, shake or break - and other applications for the FEM Programming project in TMA4220 - part 2 by Kjetil André Johannessen TMA4220 - Numerical solution of partial differential equations using the finite
More informationBuilding a finite element program in MATLAB Linear elements in 1d and 2d
Building a finite element program in MATLAB Linear elements in 1d and 2d D. Peschka TU Berlin Supplemental material for the course Numerische Mathematik 2 für Ingenieure at the Technical University Berlin,
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationExercise 1 (Formula for connection 1-forms) Using the first structure equation, show that
1 Stokes s Theorem Let D R 2 be a connected compact smooth domain, so that D is a smooth embedded circle. Given a smooth function f : D R, define fdx dy fdxdy, D where the left-hand side is the integral
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More informationBridging methods for coupling atomistic and continuum models
Bridging methods for coupling atomistic and continuum models Santiago Badia 12, Pavel Bochev 1, Max Gunzburger 3, Richard Lehoucq 1, and Michael Parks 1 1 Sandia National Laboratories, Computational Mathematics
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationEvaluation of Flux Integral Outputs for the Reduced Basis Method
Evaluation of Flux Integral Outputs for the Reduced Basis Method Jens L. Eftang and Einar M. Rønquist 20th October 2009 Abstract In this paper, we consider the evaluation of flux integral outputs from
More informationNodal O(h 4 )-superconvergence of piecewise trilinear FE approximations
Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen
More informationFinite Element Method based investigation of IPMSM losses
Finite Element Method based investigation of IPMSM losses Martin Schmidtner 1, Prof. Dr. -Ing. Carsten Markgraf 1, Prof. Dr. -Ing. Alexander Frey 1 1. Augsburg University of Applied Sciences, Augsburg,
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationSTRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD
Shape optimization 1 G. Allaire STRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD G. ALLAIRE CMAP, Ecole Polytechnique with F. JOUVE and A.-M. TOADER 1. Introduction 2. Setting of the problem 3. Shape differentiation
More informationCertified Reduced Basis Methods for Parametrized Distributed Optimal Control Problems
Certified Reduced Basis Methods for Parametrized Distributed Optimal Control Problems Mark Kärcher, Martin A. Grepl and Karen Veroy Preprint o. 402 July 2014 Key words: Optimal control, reduced basis method,
More informationOn Surface Meshes Induced by Level Set Functions
On Surface Meshes Induced by Level Set Functions Maxim A. Olshanskii, Arnold Reusken, and Xianmin Xu Bericht Nr. 347 Oktober 01 Key words: surface finite elements, level set function, surface triangulation,
More informationA Fast Fourier transform based direct solver for the Helmholtz problem. AANMPDE 2017 Palaiochora, Crete, Oct 2, 2017
A Fast Fourier transform based direct solver for the Helmholtz problem Jari Toivanen Monika Wolfmayr AANMPDE 2017 Palaiochora, Crete, Oct 2, 2017 Content 1 Model problem 2 Discretization 3 Fast solver
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationReduced-Order Greedy Controllability of Finite Dimensional Linear Systems. Giulia Fabrini Laura Iapichino Stefan Volkwein
Universität Konstanz Reduced-Order Greedy Controllability of Finite Dimensional Linear Systems Giulia Fabrini Laura Iapichino Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 364, Oktober 2017 ISSN
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationarxiv: v1 [math.na] 9 Oct 2018
PRIMAL-DUAL REDUCED BASIS METHODS FOR CONVEX MINIMIZATION VARIATIONAL PROBLEMS: ROBUST TRUE SOLUTION A POSTERIORI ERROR CERTIFICATION AND ADAPTIVE GREEDY ALGORITHMS arxiv:1810.04073v1 [math.na] 9 Oct 2018
More informationDimensionality reduction of parameter-dependent problems through proper orthogonal decomposition
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 01.2016 January 2016 (New 25.05.2016) Dimensionality
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More information3.0 FINITE ELEMENT MODEL
3.0 FINITE ELEMENT MODEL In Chapter 2, the development of the analytical model established the need to quantify the effect of the thermal exchange with the dome in terms of a single parameter, T d. In
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
More informationPeng Chen 1, Alfio Quarteroni 1,2 and Gianluigi Rozza 3
ESAI: 2AN 48 (214 943 953 DOI: 1.151/m2an/213128 ESAI: athematical odelling and Numerical Analysis www.esaim-m2an.org A WEIGHTED EPIRICAL INTERPOLATION ETHOD: APRIORI CONVERGENCE ANALYSIS AND APPLICATIONS
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationMulti-Point Constraints
Multi-Point Constraints Multi-Point Constraints Multi-Point Constraints Single point constraint examples Multi-Point constraint examples linear, homogeneous linear, non-homogeneous linear, homogeneous
More informationRobust solution of Poisson-like problems with aggregation-based AMG
Robust solution of Poisson-like problems with aggregation-based AMG Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Paris, January 26, 215 Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More informationSome improvements of Xfem for cracked domains
Some improvements of Xfem for cracked domains E. Chahine 1, P. Laborde 2, J. Pommier 1, Y. Renard 3 and M. Salaün 4 (1) INSA Toulouse, laboratoire MIP, CNRS UMR 5640, Complexe scientifique de Rangueil,
More informationOn the design of higher-order FEM satisfying the discrete maximum principle
On the design of higher-order FEM satisfying the discrete maximum principle Dmitri Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund, Germany
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationR T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T
2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce
More informationComparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems
Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Nikola Kosturski and Svetozar Margenov Institute for Parallel Processing, Bulgarian Academy of Sciences Abstract.
More informationIntrinsic finite element modeling of a linear membrane shell problem
arxiv:3.39v [math.na] 5 Mar Intrinsic finite element modeling of a linear membrane shell problem Peter Hansbo Mats G. Larson Abstract A Galerkin finite element method for the membrane elasticity problem
More informationParameterized Partial Differential Equations and the Proper Orthogonal D
Parameterized Partial Differential Equations and the Proper Orthogonal Decomposition Stanford University February 04, 2014 Outline Parameterized PDEs The steady case Dimensionality reduction Proper orthogonal
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More information