Reduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling
|
|
- Nickolas Cunningham
- 5 years ago
- Views:
Transcription
1 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 de Lausanne In collaboration with Alfio Quarteroni (EPFL & Politecnico di Milano) Supported by the ERC-Mathcard Project (ERC-2008-AdG ), the Swiss National Science Foundation (Project ), and the Emil Aaltonen Foundation Model Reduction for Complex Dynamical Systems, TU Berlin, December 2-4, 2010
2 1 Introduction 2 Navier-Stokes equations 3 Reduced basis approximation 4 Application in cardiovascular modelling
3 Challenges in modelling the human cardiovascular system Human cardiovascular system is a complex flow network of different spatial and temporal scales. When investigating fluid flow processes the flow geometries are changing over time. The geometric variation causes a strong nonlinearity in the equations. Medical professionals are interested in accurate simulation of spatial quantities, such as wall shear stresses at the location of a possible pathology. Computational costs can become unacceptably high, especially if the objective is to model the entire network, and strategies to reduce numerical efforts and model order are being developed.
4 Parametric incompressible Navier-Stokes equations (steady case) We consider the following model problem: For a given parameter vector µ D R P, find U(µ) X s.t. a(u(µ),v ; µ) = f (V ; µ) V X, µ D where U := (u,p) and V := (v,q) consist of the velocity field and the pressure, the product space X = V Q [H 1 (Ω)] 2 L 2 (Ω), and the problem consists of a linear part a 0 and a nonlinear (quadratic in U) part a 1 : a(u,v ; µ) := a 0 (U,V ; µ) + a 1 (U,U,V ; µ) U,V X, µ D
5 Parametric incompressible Navier-Stokes equations (steady case) We consider the following model problem: For a given parameter vector µ D R P, find U(µ) X s.t. a(u(µ),v ; µ) = f (V ; µ) V X, µ D where U := (u,p) and V := (v,q) consist of the velocity field and the pressure, the product space X = V Q [H 1 (Ω)] 2 L 2 (Ω), and the problem consists of a linear part a 0 and a nonlinear (quadratic in U) part a 1 : a(u,v ; µ) := a 0 (U,V ; µ) + a 1 (U,U,V ; µ) U,V X, µ D For example, if the parameter is simply µ = ν (fluid viscosity), we have a 0 (U,V ; µ) = [µ u : v p div(v) q div(u)] dω Ω a 1 (U,W,V ) = v (u )w dω Ω + appropriate boundary conditions.
6 Parametric incompressible Navier-Stokes equations (steady case) We consider the following model problem: For a given parameter vector µ D R P, find U(µ) X s.t. a(u(µ),v ; µ) = f (V ; µ) V X, µ D where U := (u,p) and V := (v,q) consist of the velocity field and the pressure, the product space X = V Q [H 1 (Ω)] 2 L 2 (Ω), and the problem consists of a linear part a 0 and a nonlinear (quadratic in U) part a 1 : a(u,v ; µ) := a 0 (U,V ; µ) + a 1 (U,U,V ; µ) U,V X, µ D For example, if the parameter is simply µ = ν (fluid viscosity), we have a 0 (U,V ; µ) = [µ u : v p div(v) q div(u)] dω Ω a 1 (U,W,V ) = v (u )w dω Ω + appropriate boundary conditions. Typically we are interested in linear functionals of the field solutions (outputs) s(µ) := l(u(µ)), i.e. need to find a reduced model s(µ) that has is within certified tolerance of the actual outputs: s(µ) s(µ) < TOL for all µ D.
7 Finite element approximation to the Navier-Stokes solution Starting from initial guess U 0, solve at each step k of a FP iteration for U k s.t. until convergence. a 0 (U k,v ; µ) + a 1 (U k 1,U k,v ) = f (V ) V V Q Stable discretization with P 2 /P 1 FE spaces for velocity and pressure V h := {v C(Ω,R d ) : v K [P 2 (K)] 2, K T h } V Q h := {q C(Ω,R) : q K P 1 (K), K T h } Q. Galerkin projection in FE space: solve at each step k for U k h s.t. until convergence. a 0 (U k h,v h; µ) + a 1 (U k 1 h,u k h,v h) = f (V h ) V h V h Q h Similar approach for the Newton s method...
8 Reduced basis approximation of the finite element solution 1 Assumption: parametric manifold of FE solutions M h X h is 1) low dimensional and 2) depends smoothly on µ (valid for small Reynolds number) 2 Choose a representative set of parameter values µ 1,..., µ N 3 Snapshot solutions u h (µ 1 ),...,u h (µ N ) span a subspace V N h for the velocity and p h (µ 1 ),...,p h (µ N ) span a subspace Q N h for the pressure 4 Galerkin reduced basis: given µ D, find U N h (µ) X N h s.t. a 0 (U k,n h,vh N ; µ) + a 1(U k 1,N h,uh k,v h N ) = f (V h N ) V h N X h N 5 Adaptive sampling procedure (greedy algorithm) for the choice of µ 1,..., µ N M = {U(µ) X ; µ D} M h = {U h (µ) X h ; µ D}
9 Reduced basis approximation of the finite element solution 1 Assumption: parametric manifold of FE solutions M h X h is 1) low dimensional and 2) depends smoothly on µ (valid for small Reynolds number) 2 Choose a representative set of parameter values µ 1,..., µ N 3 Snapshot solutions u h (µ 1 ),...,u h (µ N ) span a subspace V N h for the velocity and p h (µ 1 ),...,p h (µ N ) span a subspace Q N h for the pressure 4 Galerkin reduced basis: given µ D, find U N h (µ) X N h s.t. a 0 (U k,n h,vh N ; µ) + a 1(U k 1,N h,uh k,v h N ) = f (V h N ) V h N X h N 5 Adaptive sampling procedure (greedy algorithm) for the choice of µ 1,..., µ N Reliability / accuracy? 1 is based on the quality of the sampling 2 relies on computable and rigorous a posteriori error estimator N (µ): M = {U(µ) X ; µ D} M h = {U h (µ) X h ; µ D} X N h = span{u h(µ i ), i = 1,...,N} U h (µ) U N h (µ) X N (µ), s(µ) s N (µ) s N (µ) = l X h N(µ)
10 A posteriori error estimation of the reduced basis approximation (Veroy-Patera 2005) If τ N (µ) < 1 and β h (µ) > 0 there exists a unique solution U h (µ) s.t. Here: U h (µ) Uh N (µ) X N (µ) =: β [ h(µ) 1 ] 1 τ ρ(µ) N (µ) β h (µ) is the Babuska inf-sup constant that needs to be estimated da(u h (µ); µ)(w,v ) inf sup = β h (µ) > β 0 > 0 W X h V X h W V for the Fréchet derivative of a(u,w,v ) w.r.t first argument at U h ρ(µ) is a Sobolev embedding constant that needs to be estimated τ N (µ) := 2ρ(µ)ε N (µ) β h (µ) 2, where ε N (µ) := f ( ; µ) a(uh N, ; µ) X h is the RB residual
11 A posteriori error estimation of the reduced basis approximation (Veroy-Patera 2005) If τ N (µ) < 1 and β h (µ) > 0 there exists a unique solution U h (µ) s.t. Here: U h (µ) Uh N (µ) X N (µ) =: β [ h(µ) 1 ] 1 τ ρ(µ) N (µ) β h (µ) is the Babuska inf-sup constant that needs to be estimated da(u h (µ); µ)(w,v ) inf sup = β h (µ) > β 0 > 0 W X h V X h W V for the Fréchet derivative of a(u,w,v ) w.r.t first argument at U h ρ(µ) is a Sobolev embedding constant that needs to be estimated τ N (µ) := 2ρ(µ)ε N (µ) β h (µ) 2, where ε N (µ) := f ( ; µ) a(uh N, ; µ) X h is the RB residual Note: for large viscosity we obtain the Stokes equations and the estimator simplifies to U h (µ) Uh N (µ) N(µ) = ε } N(µ) key ingredients β h (µ)
12 Generalization of reduced basis method to parametric geometries Parametrized formulation on fixed reference domain (Rozza 2009) Evaluate output of interest s(µ) = l(u h (µ)) s.t. U h = (u h (µ),p(µ)) V h ( Ω) Q h ( Ω) solves a(u h (µ),v h ; µ) = f (V h ; µ) V h X ( Ω) v a((v,p),(w,q); µ) = ν ij (x, µ) w dω pχ ij (x, µ) w j dω qχ ij (x, µ) v j dω, Ω x i x j Ω x i Ω x i The parametrized (original) domain Ω(µ) is the image of a reference domain Ω through a parametric mapping T ( ; µ) : Ω Ω(µ) One possible parametrization using free-form deformations (L.-Rozza 2009) Transformation tensors (J T = J T (x, µ) = Jacobian of T (x, µ)) ν(x, µ) = J 1 T νo J T T J T and χ(x, µ) = J 1 T χo J T Problem reduced to a parametric PDEs system on Ω (reference domain)
13 Reduced basis offline/online computational framework Offline stage involves precomputation of structures required for the certified error estimate and choice of the reduced basis functions. Online stage has complexity only depending on N and allows evaluation of output s(µ) for any µ D with a certified error bounds.
14 Shape Optimization of Aorto-Coronaric Bypass Grafts Shape optimization of cardiovascular geometries helps to avoid post-surgical complications Local fluid patterns (vorticity) and wall shear stress are strictly related to the development of cardiovascular diseases
15 Shape Optimization of Aorto-Coronaric Bypass Grafts Shape optimization of cardiovascular geometries helps to avoid post-surgical complications Local fluid patterns (vorticity) and wall shear stress are strictly related to the development of cardiovascular diseases Shape optimization problem [Agoshkov-Quarteroni-Rozza 2006] minj(ω; v) s.t. ν v + p = f in Ω v = 0 in Ω v = v g on Γ D := Ω \ Γ out, p n + ν v n = 0 on Γ out J o (Ω;v) = v 2 dω, Ω df J o (Ω;v) = ν v Ω n tdγ o Pictures taken from: Lei et al., J Vasc. Surg. 25(1997), ; Loth et al., Annu. Rev. Fluid Mech. 40(2008),
16 Shape optimization of aorto-coronaric bypass grafts A possible free-form deformation approach [Manzoni-Quarteroni-Rozza 2010] Several analyses show a deep impact of the graft-artery diameter ratio Φ and anastomotic angle α on shear stress and vorticity distributions Oscillatory shear stress with different graft-artery diameter ratios Φ and anastomotic angles α. Picture taken from F.L. Xiong, C.K. Chong, Med. Eng. & Phys. 30(2008),
17 Shape optimization of aorto-coronaric bypass grafts A possible free-form deformation approach [Manzoni-Quarteroni-Rozza 2010] Several analyses show a deep impact of the graft-artery diameter ratio Φ and anastomotic angle α on shear stress and vorticity distributions Oscillatory shear stress with different graft-artery diameter ratios Φ and anastomotic angles α. Picture taken from F.L. Xiong, C.K. Chong, Med. Eng. & Phys. 30(2008), In order to get a low-dimensional FFD parametrization we need to maximize the influence of the control points by placing them close to the sensitive regions 8 parameters (7 vertical and 1 horizontal displacements) to control the anastomotic angle, the graft-artery diameter ratio, the upper side, the lower wall
18 Shape optimization of aorto-coronaric bypass grafts! N (µ) RB approximation space construction RB space dimension N Number of FE dof N v + N p Lattice FFD control points P i,j 5 6 Number of design variables P 8 Number of RB functions N 22 Error tolerance RB greedy εtol RB Affine operator components Q 222 Error estimation (energy norm) for RB space construction (greedy procedure) and selected snapshots µ µ 1 µ 2 µ 3 µ 4 µ 6 µ 7 µ Reduction in linear system dimension 500:1 Computational speedup (single flow simulation) 107 Reduction in parametric complexity w.r.t. explicit nodal deformation 102:1
19 Shape optimization of aorto-coronaric bypass grafts Vorticity Minimization (downfield region) Automatic iterative minimization procedure (sequential quadratic programming) Vorticity evaluation by using the reduced basis velocity at each step [Manzoni-Quarteroni-Rozza 2010] Optimized bypass anastomosis and Stokes flow (velocity magnitude and pressure) Optimal (black) and unperturbed (grey) configurations of FFD parametrization Vorticity magnitude for the unperturbed (left) and optimal (right) configuration
20 Reduction of time-dependent Navier-Stokes (ongoing work) One-dimensional prototype problem (viscous Burgers equation): find u(µ) L 2 (0,T ;V ) C 0 ([0,T ];L 2 (Ω)) s.t. d dt Ω (u(µ),v) + µ u x (µ)v x dx 1 2 u 2 (µ)v x dx = f (v) Ω where Ω = (0,1) and V H 1 (Ω). v V Time discretization with implicit Euler = time-discrete equations Spatial discretization with FEM + reduced basis reduction as before Stability constant ρ N := inf v V da(u N,k h )(v,v) not necessarily positive! v X A posteriori estimator (Nguyen-Rozza-Patera 2009) for k = 1,..., t T ( ) uh k (µ) un,k h (µ) k t µ N (µ) := k m=1 εn 2 (tm ; µ) m 1 j=1 (1 + tρ N(t j ; µ)) k m=1 (1 + tρ N(t m ; µ)) BUT the error bound grows exponentially in time
21 Conclusions Reduced basis methods a reliable MOR method for parametric PDEs Parameters can also describe the (variable) flow geometry Certified error bounds for spatial outputs of reduced field variables Extensions to noncoercive (Stokes) and nonlinear (Navier-Stokes) cases Future work Time-dependent Navier-Stokes, improved error estimates Reduction of coupled multiphysics problems Parameter identification and inverse problems
22 References V.I. Agoshkov, A. Quarteroni and G. Rozza. Shape design in aorto-coronaric bypass using perturbation theory. SIAM J. Num. Anal. 44(1), , D.B.P. Huynh, D.J. Knezevic, Y. Chen, J.S. Hesthaven, and A.T. Patera. A natural-norm successive constraint method for inf-sup lower bounds. Comput. Methods Appl. Mech. Engrg. 199: , T. L. and G. Rozza. Parametric free-form shape design with PDE models and reduced basis method. Comput. Methods Appl. Mech. Engrg. 199(23-24): , A. Manzoni, A. Quarteroni, G. Rozza. Shape optimization for viscous flows by reduced basis methods and free-form deformation, J. Comp. Phys., submitted, A. Quarteroni and G. Rozza. Numerical solution of parametrized Navier-Stokes equations by reduced basis methods. Numer. Methods Partial Differential Equations, 23(4): , N.-C. Nguyen, G. Rozza and A.T. Patera. Reduced basis approximation and a posteriori error estimation for the time-dependent viscous Burgers equation. Calcolo 46(3): , G. Rozza, D.B.P. Huynh, A.T. Patera. Reduced basis approximation and a posteriori error estimation for affinely parametrized elliptic coercive PDEs. Arch. Comput. Methods Engrg.,15: , G. Rozza and K. Veroy. On the stability of the reduced basis method for Stokes equations in parametrized domains. Comput. Methods Appl. Mech. Engrg., 196(7): , K. Veroy and A.T. Patera. Certified real-time solution of the parametrized steady incompressible Navier-Stokes equations; rigorous reduced-basis a posteriori error bounds. Int. J. Numer. Methods Fluids 47(89), , 2005.
Some 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 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 informationModel 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 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 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 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 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 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 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 informationReduced basis method for multi-parameter dependent steady Navier Stokes equations: applications to natural convection in a cavity
Reduced basis method for multi-parameter dependent steady Navier Stokes equations: applications to natural convection in a cavity Simone Deparis, Gianluigi Rozza Ecole Polytechnique Fédérale de Lausanne,
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 informationDownloaded 05/29/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 34, No. 2, pp. A1187 A1213 c 2012 Society for Industrial and Applied Mathematics A REDUCED BASIS MODEL WITH PARAMETRIC COUPLING FOR FLUID-STRUCTURE INTERACTION PROBLEMS TONI LASSILA,
More informationModel 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 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 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 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. 22.2013 July 2013 Model order reduction in fluid
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 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 informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
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 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 informationA model order reduction framework for parametrized nonlinear PDE constrained optimization
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 11.2015 Mai 2015 A model order reduction framework
More informationHIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 1, March 2008 HIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IOAN TELEAGA AND JENS
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 informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
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 informationOne-Shot Domain Decomposition Methods for Shape Optimization Problems 1 Introduction Page 565
1 One-Shot Domain Decomposition Methods for Shape 2 Optimization Problems 3 Rongliang Chen 1 and Xiao-Chuan Cai 2 4 1 College of Mathematics and Econometrics, Hunan University, Changsha, Hunan 410082,
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
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 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 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 informationREDUCED BASIS METHOD
REDUCED BASIS METHOD approximation of PDE s, interpolation, a posteriori estimates Yvon Maday Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie, Paris, France, Institut Universitaire de
More informationNumerical modelling of shear-thinning non-newtonian flows in compliant vessels
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2007; 00:1 [Version: 2002/09/18 v1.01] Numerical modelling of shear-thinning non-newtonian flows in compliant vessels M.
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 informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
More informationModel reduction of parametrized aerodynamic flows: discontinuous Galerkin RB empirical quadrature procedure. Masayuki Yano
Model reduction of parametrized aerodynamic flows: discontinuous Galerkin RB empirical quadrature procedure Masayuki Yano University of Toronto Institute for Aerospace Studies Acknowledgments: Anthony
More informationA Petrov-Galerkin reduced basis approximation of the Stokes equation in parametrized geometries
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 35.2014 August 2014 A Petrov-Galerkin reduced basis
More informationReceived *****; accepted after revision Presented by. arxiv: v1 [math.ap] 23 Apr 2015
, Received *****; accepted after revision +++++ Presented by arxiv:1504.06131v1 [math.ap] 23 Apr 2015 1 Simultaneous Empirical Interpolation and Reduced Basis method for non-linear problems Cécile Daversin,
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 informationCoupling Stokes and Darcy equations: modeling and numerical methods
Coupling Stokes and Darcy equations: modeling and numerical methods Marco Discacciati NM2PorousMedia 24 Dubrovnik, September 29, 24 Acknowledgment: European Union Seventh Framework Programme (FP7/27-23),
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 informationneed to improve the plant availability enhancing the performance for Nuclear Power Plants (NPPs) + a more and more ability to follow grid demands
need to improve the plant availability enhancing the performance for Nuclear Power Plants (NPPs) + a more and more ability to follow grid demands Classic control techniques Example: Proportional-Integral-
More informationNumerical Solution of Parametrized Navier Stokes Equations by Reduced Basis Methods
Numerical Solution of Parametrized Navier Stokes Equations by Reduced Basis Methods Alfio Quarteroni, 1,2 Gianluigi Rozza 2, * 1 Modellistica e Calcolo Scientifico (MOX), Dipartimento di Matematica F.Brioschi,
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationFourier Series Development for Solving a Steady Fluid Structure Interaction Problem
Int. J. Contemp. Math. Sciences, Vol. 8, 203, no. 2, 75-84 HIKARI Ltd, www.m-hikari.com Fourier Series Development for Solving a Steady Fluid Structure Interaction Problem Ibrahima Mbaye University of
More informationReduced Basis Methods for Nonlinear Parametrized Partial Differential Equations
Reduced Basis Methods for Nonlinear Parametrized Partial Differential Equations M. Grepl Aachen Institute for Advanced Study in Computational Engineering Science AICES Progress since 2009 Overview of 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 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 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 informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationToni Lassila 1, Andrea Manzoni 3, Alfio Quarteroni 1,2 and Gianluigi Rozza 3
ESAIM: M2AN 47 (2013) 1107 1131 DOI: 10.1051/m2an/2012059 ESAIM: Mathematical Modelling and Numerical Analysis www.esaim-m2an.org BOUNDARY CONTROL AND SHAPE OPTIMIZATION FOR THE ROBUST DESIGN OF BYPASS
More informationGreedy control. Martin Lazar University of Dubrovnik. Opatija, th Najman Conference. Joint work with E: Zuazua, UAM, Madrid
Greedy control Martin Lazar University of Dubrovnik Opatija, 2015 4th Najman Conference Joint work with E: Zuazua, UAM, Madrid Outline Parametric dependent systems Reduced basis methods Greedy control
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA REDUCED BASIS METHOD FOR THE HAMILTON-JACOBI-BELLMAN EQUATION WITH APPLICATION TO THE EUROPEAN UNION EMISSION TRADING SCHEME
A REDUCED BASIS METHOD FOR THE HAMILTO-JACOBI-BELLMA EQUATIO WITH APPLICATIO TO THE EUROPEA UIO EMISSIO TRADIG SCHEME SEBASTIA STECK AD KARSTE URBA Abstract. This paper draws on two sources of motivation:
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationReduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems. Laura Iapichino Stefan Trenz Stefan Volkwein
Universität Konstanz Reduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems Laura Iapichino Stefan Trenz Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 347, Dezember 2015
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
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 informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationReduced Basis Method for Explicit Finite Volume Approximations of Nonlinear Conservation Laws
Proceedings of Symposia in Applied Mathematics Reduced Basis Method for Explicit Finite Volume Approximations of Nonlinear Conservation Laws Bernard Haasdonk and Mario Ohlberger Abstract. The numerical
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 informationPALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations
2013 SIAM Conference On Computational Science and Engineering Boston, 27 th February 2013 PALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations U. Villa,
More informationGreedy Control. Enrique Zuazua 1 2
Greedy Control Enrique Zuazua 1 2 DeustoTech - Bilbao, Basque Country, Spain Universidad Autónoma de Madrid, Spain Visiting Fellow of LJLL-UPMC, Paris enrique.zuazua@deusto.es http://enzuazua.net X ENAMA,
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
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 informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationSpace time finite element methods in biomedical applications
Space time finite element methods in biomedical applications Olaf Steinbach Institut für Angewandte Mathematik, TU Graz http://www.applied.math.tugraz.at SFB Mathematical Optimization and Applications
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 informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
More informationA Multiphysics Strategy for Free Surface Flows
A Multiphysics Strategy for Free Surface Flows Edie Miglio, Simona Perotto, and Fausto Saleri MOX, Modeling and Scientific Computing, Department of Mathematics, Politecnico of Milano, via Bonardi 9, I-133
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationReduced-Basis Approximations and A Posteriori Error Bounds for Nonaffine and Nonlinear Partial Differential Equations: Application to Inverse Analysis
Reduced-Basis Approximations and A Posteriori Error Bounds for Nonaffine and Nonlinear Partial Differential Equations: Application to Inverse Analysis by Nguyen Ngoc Cuong B.Eng., HCMC University of Technology
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationTowards parametric model order reduction for nonlinear PDE systems in networks
Towards parametric model order reduction for nonlinear PDE systems in networks MoRePas II 2012 Michael Hinze Martin Kunkel Ulrich Matthes Morten Vierling Andreas Steinbrecher Tatjana Stykel Fachbereich
More informationAN ADAPTIVE REDUCED BASIS COLLOCATION METHOD BASED ON PCM ANOVA DECOMPOSITION FOR ANISOTROPIC STOCHASTIC PDES
AN ADAPTIVE REDUCED BASIS COLLOCATION METHOD BASED ON PCM ANOVA DECOMPOSITION FOR ANISOTROPIC STOCHASTIC PDES Heyrim Cho & Howard Elman 2 Department of Mathematics, University of Maryland, College Park,
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationNOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7 Number Pages 3 39 c Institute for Scientific Computing and Information NOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationGoal-oriented error estimation for the reduced basis method, with application to sensitivity analysis
Goal-oriented error estimation for the reduced basis method, with application to sensitivity analysis Alexandre Janon Maëlle Nodet Clémentine Prieur July 10, 2014 Abstract The reduced basis method is a
More informationRICE UNIVERSITY. Nonlinear Model Reduction via Discrete Empirical Interpolation by Saifon Chaturantabut
RICE UNIVERSITY Nonlinear Model Reduction via Discrete Empirical Interpolation by Saifon Chaturantabut A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy
More informationImplementation of 3D Incompressible N-S Equations. Mikhail Sekachev
Implementation of 3D Incompressible N-S Equations Mikhail Sekachev Navier-Stokes Equations The motion of a viscous incompressible fluid is governed by the Navier-Stokes equations u + t u = ( u ) 0 Quick
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 informationarxiv: v1 [math.na] 30 Jul 2018
The Effort of Increasing Reynolds Number in Projection-Based Reduced Order Methods: from Laminar to Turbulent Flows Saddam Hijazi, Shafqat Ali, Giovanni Stabile, Francesco Ballarin and Gianluigi Rozza
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationA NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES EQUATIONS
MATHEMATICS OF COMPUTATION Volume 74, Number 250, Pages 543 554 S 0025-5718(04)01711-9 Article electronically published on August 31, 2004 A NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationA generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem
A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem Dave McCormick joint work with James Robinson and José Rodrigo Mathematics and Statistics Centre for Doctoral Training University
More informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
More informationElectronic Transactions on Numerical Analysis Volume 32, 2008
Electronic Transactions on Numerical Analysis Volume 32, 2008 Contents 1 On the role of boundary conditions for CIP stabilization of higher order finite elements. Friedhelm Schieweck. We investigate the
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationThe Generalized Empirical Interpolation Method: Analysis of the convergence and application to data assimilation coupled with simulation
The Generalized Empirical Interpolation Method: Analysis of the convergence and application to data assimilation coupled with simulation Y. Maday (LJLL, I.U.F., Brown Univ.) Olga Mula (CEA and LJLL) G.
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More information