REDUCED BASIS METHOD
|
|
- Arthur Grant
- 5 years ago
- Views:
Transcription
1 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 France and Division of Applied Maths, Brown University CEMRACS SUMMER SCHOOL Luminy, July 2016
2 This talk deals with linear constructive approximation methods specifically tailored to approximate the functions in F (or in a set close to F see PBDW).
3 Small Kolomogorov manifold This talk deals with linear constructive approximation methods specifically tailored to approximate the functions in F (or in a set close to F see PBDW).
4 Small Kolomogorov manifold Let F be a compact subset of a functional Banach space X This talk deals with linear constructive approximation methods specifically tailored to approximate the functions in F (or in a set close to F see PBDW).
5 Small Kolomogorov manifold Let F be a compact subset of a functional Banach space X An example can be a set of parametric functions F = {u(., µ),µ2 D} This talk deals with linear constructive approximation methods specifically tailored to approximate the functions in F (or in a set close to F see PBDW).
6 Small Kolomogorov manifold Let F be a compact subset of a functional Banach space X An example can be a set of parametric functions F = {u(., µ),µ2 D} This talk deals with linear constructive approximation methods specifically tailored to approximate the functions in F (or in a set close to F see PBDW). For this, we consider approximations on appropriate n-dimensional spaces X n of relatively small dimension under the hypothesis that the Kolmogorov n-width of F is small.
7
8 The Kolmogorov n-width of F in X, is defined as d n (F, X )= inf X n X dim(xn)applen sup u2f inf ku vk X, ( v2x n And the assumption is that d n (F, X ) is small.
9 The Kolmogorov n-width of F in X, is defined as d n (F, X )= inf X n X dim(xn)applen sup u2f inf ku vk X, ( v2x n And the assumption is that d n (F, X ) is small. In order to satisfy the constraint on the computing time, we consider continuous linear approximations J n : X! X n on appropriate n-dimensional spaces X n of relatively small dimension. A first bound on the approximation error is then supku u2f J n [u]k X apple (1 + n )sup u2f inf ku vk X, (1) v2x n where n =sup '2X kj n [']k X k'k X (2) is the norm of the operator J n : X! X n and is known as the Lebesgue constant.
10 The Kolmogorov n-width of F in X, is defined as d n (F, X )= inf X n X dim(xn)applen sup u2f inf ku vk X, ( v2x n And the assumption is that d n (F, X ) is small. What is the proper definition of X n In order to satisfy the constraint on the computing time, we consider continuous linear approximations J n : X! X n on appropriate n-dimensional spaces X n of relatively small dimension. A first bound on the approximation error is then supku u2f J n [u]k X apple (1 + n )sup u2f inf ku vk X, (1) v2x n where n =sup '2X kj n [']k X k'k X (2) is the norm of the operator J n : X! X n and is known as the Lebesgue constant.
11 EIM/GEIM In 2004 with M. Barrault, N. C. Nguyen, and A. T. Patera, we proposed n generic approach : the Empirical Interpolation Method: Application to cient Reduced-Basis Discretization Of Partial Di erential Equations, that as proven successful This approach allows to determine an empirical optimal set of interpolation points and/or set of interpolating functions. In 2013, with Olga Mula, we have generalized it (GEIM) to include more general output from the functions we want to interpolate : not only pointwize values but also some moments.
12 The first generating function is ' 1 = arg max '2F k'( )k L1 ( ), the associated interpolation point satisfies x 1 = arg max x2 ' 1 (x), we then set q 1 = ' 1 ( )/' 1 (x 1 ) and B11 1 = 1. Then the construction proceeds by induction : assume the nested sets of interpolation points M 1 = {x 1,...,x M 1 },M apple M max, and the associated nested sets of basis functions {q 1,...,q M 1 } are given 1. We first solve the interpolation problem : Find such that I M 1 ['( )] = 1) where Mmax M is some given upper bound fixed a priori MX 1 j=1 M 1,j [']q j, I M 1 ['( )](x i )='(x i ), i =1,...,M 1, that allows to define the M 1,j ['], 1 apple j apple M 1, as it can be proven indeed that the (M 1) (M 1) matrix of running entry q j (x i ) is invertible, actually it is lower triangular with unity diagonal.
13 We then set (note indeed that h M i (x j )= ij ). 8' 2 F, " M 1 (') =k' I M 1 [']k L1 ( ), and define and x M = arg max x2 ' M = arg max '2F " M 1('), ' M (x) J M 1 [' M ](x), we finally set r M (x) =' M (x) J M 1 [' M (x)], q M = r M /r M (x M ) and B M ij = q j (x i ), 1 apple i, j apple M The Lagrangian functions that can be used to build the interpolation operator I M in X M = span {' i, 1 apple i apple M} = span {q i, 1 apple i apple M} over the set of points M = {x i, 1 apple i apple M} verify for any given M, I M [u( )] = MX u(x i )h M i ( ) i=1 where h M i ( )= MX q j ( )[B M ] j=1 1 ji
14 For GEIM, we assume now that we do not have access to the values of ' 2 F at points in easily, but, on the contrary, that we have a dictionary of linear forms 2 assumed to be continuous in some sense, e.g. in L 2 ( ) with norm 1 the application of which over each ' 2 F is easy. Our extension consists in defining ' 1, ' 2,..., ' M and a family of associated linear forms 1, 2,..., M such that the following generalized interpolation process (our GEIM) is well defined : J M ['] = MX j=1 j ' j, such that 8i =1,...,M, i(j M [']) = i (') (1) Note that the GEIM reduces to the EIM when the dictionary is composed of dirac masses, defined in the dual space of C 0 ( ).
15
16 Note that obviously everything has to be implemented on a computer nd thus discretized!!
17 Note that obviously everything has to be implemented on a computer nd thus discretized!! Note also that for some reasons you may want to use your proper set f basis/interpolating function in your preferredspace X N that may come from ntuition, previous knowledge, or POD/SVD.
18 Note that obviously everything has to be implemented on a computer nd thus discretized!! Note also that for some reasons you may want to use your proper set f basis/interpolating function in your preferredspace X N that may come from ntuition, previous knowledge, or POD/SVD.
19 Note that obviously everything has to be implemented on a computer nd thus discretized!! Note also that for some reasons you may want to use your proper set f basis/interpolating function in your preferredspace X N that may come from ntuition, previous knowledge, or POD/SVD. Journal: Communications on Pure and Applied Analysis - A general multipurpose interpolation procedure: the magic points Yvon Maday - Ngoc Cuong Nguyen - Anthony T. Patera - George S. H. Pau Pages: , Volume 8, Issue 1, January 2009 doi: /cpaa
20
21 The important thing is to measure to which extent the EIM/GEIM misses the optimality of the best approximation in the best optimal discrete space suggested by the definition of the Kolmogorov n-width.
22 Formula suggests that n plays an important role in the result and it is therefore important to discuss its behavior as n increases. First of all, n depends both on the choices of the interpolating functions and interpolation points. We have proven (YM-Mula-Patera-Yano) that n =1/ n,where n = inf '2X n sup 2Span{ 0,..., n 1 } h', i X,X 0 k'k X k k X 0.
23 Formula suggests that n plays an important role in the result and it is therefore important to discuss its behavior as n increases. First of all, n depends both on the choices of the interpolating functions and interpolation points. We have proven (YM-Mula-Patera-Yano) that n =1/ n,where n = inf '2X n sup 2Span{ 0,..., n 1 } h', i X,X 0 k'k X k k X 0. GEIM interpreted as an oblic projection
24
25
26
27 We are now working on the justification of this nice behavior.
28 The next point is to justify the choice of the interpolating functions. It can be proven that Lemma For any n 1, the nth interpolating function ' n verifies k' n P n (' n )k X 1+ n max f2f kf P n(f)k X. Which allows to use the frame weak greedy of the papers by P. Binev, A. Cohen, W. Dahmen, R.A. DeVore, G. Petrova, and P. Wojtaszczyk, and R. A. DeVore, G. Petrova, and P. Wojtaszczyk, to analyse the convergence properties of our algorithm Work done with O. Mula and G. Turinici (SINUM (2016)
29 In a nutshell, in the case where we have a Hilbert framework, our result states that Theorem If ( n ) 1 n=1 is a monotonically increasing sequence then i) if d n apple C 0 n for any n 1, then n apple C 0 nn,with n := n, if n 2. ii) if d n apple C 0 e c 1n for n 1 and C 0 1, then n apple C 0 ne c 2n,with n := p 2 n, if n 2. Work done with O. Mula and G. Turinici (SINUM (2016)
30 Another application : fluid flow
31 Another application : fluid flow
32
33
34 Application of EIM for Reduced Basis Methods for non affine and nonlinear PDE: see Efficient reduced-basis treatment of nonaffine and nonlinear partial differential equations MA Grepl, Y Maday, NC Nguyen, AT Patera ESAIM (2007) or the two recent books Certified Reduced Basis Methods for Parametrized Partial Differential Equations Authors: Hesthaven, Jan S, Rozza, Gianluigi, Stamm, Benjamin Reduced Basis Methods for Partial Differential Equations An Introduction Authors: Quarteroni, Alfio, Manzoni, Andrea, Negri, Federico
35 More about (G)EIM
36 In a recent paper, with J.P. Argaud, B. Bouriquet, H. Gong and O. Mula, we have introduced the GEIM to monitor nuclear reactor. The challenge here is that we do not have access inside the core of the reactor and thus have only the possibility to place the captors inside the surrounding region.
37 The two group di usion equation in matrix notation reads A(µ)' = 1 k eff F (µ)' Where µ is the parameters set, e.g. D,, f. A and F are 2 2 matrix and ' is a 2-element column vector: r D A(µ) = 1 r +( 1 a + 1!2 s ) 0 1!2 s r D 2 r + 2 a 1 1 f 1 2 f F (µ) = 2 1 f 2 2 f '1 ' = Where D i, i =1, 2 is called the di usion coe cient of each group; i a, i = 1, 2 is the absorption cross section of each group; ' i,i=1, 2 is the neutron flux of each group; 1!2 s is called the removal cross section from group 1 to group 2; 1 f, i =1, 2 is the fission source term of each group; i, i =1, 2 is called the fission spectrum of each group; finally k eff is the e ective multiplication factor, also the eigenvalue of equation. ' 2
38 In 1D, this looks like
39 In order you are convinced that the Kolmogorov dimension is small
40 First classical EIM
41
42 Then with a restriction on the position of the points
43
44 Similar results in two dimensions The first twenty interpolation points distribution (constrained in fuel region)
45 Similar results in two dimensions
46 Similar results in two dimensions
47 What about noisy data We clearly see the effect of the Lebesgue constant
48 What about noisy data : how to minimize the effect of the Lebesgue constant
49 Application of GEIM in data assimilation and monitoring
50 Application of GEIM in data assimilation and monitoring
51 Application of GEIM in data assimilation and monitoring
52 Incorporating the model error : Parametrized-Background Data-Weak (PBDW) formulation with A.T. Patera, J. D. Penn and M. Yano The PBDW formulation integrates a parametrized mathematical model and M experimental observations associated with the configuration C to estimate the true field u true [C] as well as any desired output l out (u true [C]) 2 C for given output functional l out. We first introduce a sequence of background spaces that reflect our (prior) best knowledge, Z 1 Z Nmax U; here the second ellipsis indicates that we may consider the sequence of length N max as resulting from a truncation of an infinite sequence. Our goal is to choose the background spaces such that In words, we choose the background spaces such that the most dominant physics that we anticipate to encounter for various system configurations is well represented for a relatively small N.
53 Incorporating the model error : Parametrized-Background Data-Weak (PBDW) formulation with A.T. Patera, J. D. Penn and M. Yano Concerning the form of the noises (e m ) m, we make the following three assumptions:
54 We first associate with each observation functional `om 2 U 0 an observable function,
55
56 Analysis
57 Analysis
58 Experimental settings
59
60
61 Data-Driven Empirical Enrichment of the Background Space
62
63
64 Conclusion
65 Conclusion
66 Conclusion EIM and GEIM have been used now extensively in order to implement easily the reduced basis method
67 Conclusion EIM and GEIM have been used now extensively in order to implement easily the reduced basis method From its initial magic point statement to empirical interpolation we have been now able (thanks to the work initiated by devore on greedy algorithms) to state it on more firm ground
68 Conclusion EIM and GEIM have been used now extensively in order to implement easily the reduced basis method From its initial magic point statement to empirical interpolation we have been now able (thanks to the work initiated by devore on greedy algorithms) to state it on more firm ground There has been some nice generalizations like the matrix form in the group of A. Quartering and G. Rozza or operator interpolation M. Drohmann, B. Haasdonk, M. Ohlberger
69 Conclusion EIM and GEIM have been used now extensively in order to implement easily the reduced basis method From its initial magic point statement to empirical interpolation we have been now able (thanks to the work initiated by devore on greedy algorithms) to state it on more firm ground There has been some nice generalizations like the matrix form in the group of A. Quartering and G. Rozza or operator interpolation M. Drohmann, B. Haasdonk, M. Ohlberger There is still some missing point on the Lebesgue constant behavior.. on which we are passive but it works still empirically (see however the interpretation of the greedy algorithm in this respect)
70 Conclusion EIM and GEIM have been used now extensively in order to implement easily the reduced basis method From its initial magic point statement to empirical interpolation we have been now able (thanks to the work initiated by devore on greedy algorithms) to state it on more firm ground There has been some nice generalizations like the matrix form in the group of A. Quartering and G. Rozza or operator interpolation M. Drohmann, B. Haasdonk, M. Ohlberger There is still some missing point on the Lebesgue constant behavior.. on which we are passive but it works still empirically (see however the interpretation of the greedy algorithm in this respect) The pure interpolation process can be enriched by more data as in the PBDW method with A.T. Patera, J. Penn and M. Yano
71 Conclusion EIM and GEIM have been used now extensively in order to implement easily the reduced basis method From its initial magic point statement to empirical interpolation we have been now able (thanks to the work initiated by devore on greedy algorithms) to state it on more firm ground There has been some nice generalizations like the matrix form in the group of A. Quartering and G. Rozza or operator interpolation M. Drohmann, B. Haasdonk, M. Ohlberger There is still some missing point on the Lebesgue constant behavior.. on which we are passive but it works still empirically (see however the interpretation of the greedy algorithm in this respect) The pure interpolation process can be enriched by more data as in the PBDW method with A.T. Patera, J. Penn and M. Yano
72 The PBDW formulation is endowed with the following characteristics: Weak formulation. Actionable a priori theory. The weak formulation facilitates the construction of a priori error estimates Background space that best reflects our (prior) best knowledge of the phenomenon under consideration Design of quasi-optimal set of observations from a library of experimentally realizable observations in order to maximize the stability of the data assimilation. Correction of unmodeled physics with uncertainty. Online computational cost is O(M). We may realize real-time state estimation Simple non-intrusive implementation and generality. The mathematical model appears only in the offline stage. Recently a work of P. Binev, A. Cohen, W. Dahmen, R. DeVore, G. Petrova, P. Wojtaszczyk has analyzed this approach in terms of optimal recovery and our algorithm is optimal in this frame. They have extended it to the Multi-Space Case.
73 Reduced Basis method for a bifurcation study of a thermal convection problem Henar Herrero, Yvon Maday, Francisco Pla
74 RB-4-RB : Formulation of the problem A sketch of the domain and the physical situation is shown here The domain is a rectangle of depth d and width L. The domain contains a fluid that is heated from below, so that on the bottom plate a temperature T0 is imposed and on the upper plate the temperature is where T 1 = T 0 T = T 0 d is the vertical temperature gradient.
75 Formulation of the problem The equations governing the system are the incompressible Navier-Stokes equations with the Boussinesq approximation coupled with a heat equation
76 Bifurcation diagram
77 Stationary equations : Legendre collocation method
78 Greedy procedure In this calculations we have considered expansions of order n=35 in the x-direction and m=13 in the z-direction.
79 Greedy procedure
80 Greedy procedure
81 Greedy procedure
82 Galerkin procedure for each branch
83 Galerkin procedure for each branch
84 Galerkin procedure for each branch
85 A new concept of reduced basis approximation for convection dominated problems Nicolas Cagniart, Yvon Maday, Benjamin Stamm
86 And if the Kolmogorov n-width is not small? Try to better look at the set of solutions
87 A case where the dimension is not small 1D viscous Burger equation Snapshots of the solution to the unsteady viscous burger equation with u 0 = sin(x), =4, =0.04
88 A case where the dimension is not small 1D viscous Burger equation Snapshots of the solution to the unsteady viscous burger equation with u 0 = + sin(x), =4, =0.04
89 Nevertheless all solutions look alike what we have done above is to center the solutions these are the u(x n,t n )
90 and the dimension diminishes largely!! Eigenvalues of the POD decomposition of the original set of snapshots (in red) and of the centered set of snapshots (in green)
91 and the dimension diminishes largely!! Reconstruction of a snapshot (blue) using 3 POD modes. Left figure is in the centered case. Right figure is the uncentered case
92 The discrete scheme we want to mimic u n+1 (,µ)=u n (,µ) dt u n (,µ)u n x(,µ)+dt u n xx(,µ) leads to At each time step, we are looking for a set of reduced coordinates { i n} and a translation parameter n such that our true solution u N (,µ,t n )iswell apprxoximated by : u N (,µ,t n )= X i n i (µ) i ( n ) (1) We replace the search of the real and s as in the usual galerkin method. That is, we are trying to minimize the residual. min n+1 2 min n+1 2R X i n+1 i i( n+1 ) u n dt u n u n x + dt u n xx 2 (2) Complexity.. how to compute and solve??
93 We have assumed periodic boundary conditions (for the sake of simplicity) thus we need to compute online the following terms : 8 R 8, 8i, j, >< i ( ) j ( ) R 8, 8i, j, p >: i ( ) j ( )( p ) x ( ) R 8, 8i, j, ( i) x ( )( j ) x ( ) for unknown value (to be optimized) of For a su ciently small time step, we expect to be of order t c where c is some local characteristic velocity. We have chosen the following method: precompute the scalar products for a predefined small set of discrete values of in [ t c max, t c max ], where c max is the maximum expected shock speed during the simulation. using the regularity of the scalar product, we use spline interpolation to get and approximate value for any in [ t c max, t c max ]
94 The RB discrete scheme thus iterates between the evaluation of the proper best translation definition of the coefficient n+1 i n+1 and the proper X i n+1 i i( n+1 ) Mean reconstruction error w.r.t number of POD basis used
95 extensions.. what needs to be done - non periodic (superposition of basic space and convection space) - higher dimensions (POD representation of the translations ) see 1) - better fitting (add the derivatives of the POD functions) - better fitting (replace least square with L 1, see 2) approximation) 1) Iollo, A., Lombardi, D.: Advection modes by optimal mass transfer. Physical Review E 89(2), (2014) 2) Technical paper: Roxana Crisovan, Rémi Abgrall, David Amsallem Robust Model Reduction by L1-norm Minimization and Approximation via Dictionaries: Application to Linear and Nonlinear Hyperbolic Problems
Empirical Interpolation Methods
Empirical Interpolation Methods Yvon Maday Laboratoire Jacques-Louis Lions - UPMC, Paris, France IUF and Division of Applied Maths Brown University, Providence USA Doctoral Workshop on Model Reduction
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 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 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 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 informationReduced Modeling in Data Assimilation
Reduced Modeling in Data Assimilation Peter Binev Department of Mathematics and Interdisciplinary Mathematics Institute University of South Carolina Challenges in high dimensional analysis and computation
More informationThe Generalized Empirical Interpolation Method: Stability theory on Hilbert spaces with an application to the Stokes equation
The Generalized Empirical Interpolation Method: Stability theory on Hilbert spaces with an application to the Stokes equation The MIT Faculty has made this article openly available. Please share how this
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 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 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 informationrecent developments of approximation theory and greedy algorithms
recent developments of approximation theory and greedy algorithms Peter Binev Department of Mathematics and Interdisciplinary Mathematics Institute University of South Carolina Reduced Order Modeling in
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 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 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 informationDEIM-BASED PGD FOR PARAMETRIC NONLINEAR MODEL ORDER REDUCTION
VI International Conference on Adaptive Modeling and Simulation ADMOS 213 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) DEIM-BASED PGD FOR PARAMETRIC NONLINEAR MODEL ORDER REDUCTION JOSE
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 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 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 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 informationSurrogate gravitational waveform models
Department of Physics, University of Maryland Theoretical astrophysics seminar Cornell University, November 30, 2013 Previous works - Field, Galley, Herrmann, Hesthaven, Ochsner, Tiglio (Phys. Rev. Lett.
More informationPBDW: a non-intrusive Reduced Basis Data Assimilation Method and its application to outdoor Air Quality Models
PBDW: a non-intrusive Reduced Basis Data Assimilation Method and its application to outdoor Air Quality Models Janelle K Hammond, Rachida Chakir, Frédéric Bourquin, Yvon Maday To cite this version: Janelle
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 informationA Reduced-Basis Approach to the Reduction of Computations in Multiscale Models & Uncertainty Quantification
A Reduced-Basis Approach to the Reduction of Computations in Multiscale Models & Uncertainty Quantification Sébastien Boyaval1,2,3 1 Univ. Paris Est, Laboratoire d hydraulique Saint-Venant (EDF R&D Ecole
More informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Methods for Nonlinear Systems Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu 1 / 65 Outline 1 Nested Approximations 2 Trajectory PieceWise Linear (TPWL)
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 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 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 informationMasayuki Yano. Appointments and research experience Assistant Professor, Institute for Aerospace Studies, University of Toronto 10/2015 present
Education Masayuki Yano University of Toronto Institute for Aerospace Studies 4925 Dufferin St, Toronto, M3H 5T6, Canada E-mail: myano@utias.utoronto.ca Web: http://arrow.utias.utoronto.ca/~myano/ Massachusetts
More informationarxiv: v1 [math.na] 8 Mar 2016
arxiv:1603.02652v1 [math.na] 8 Mar 2016 ROBUST MODEL REDUCTION OF HYPERBOLIC PROBLEMS BY L 1 -NORM MINIMIZATION AND DICTIONARY APPROXIMATION R. ABGRALL (1), D. AMSALLEM (2) AND R. CRISOVAN (1) (1): INSTITUT
More informationLocal null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain
Local null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain Nicolás Carreño Université Pierre et Marie Curie-Paris 6 UMR 7598 Laboratoire
More informationA general, multipurpose interpolation procedure: the magic points
A general, multipurpose interpolation procedure: the magic points Yvon aday, Ngoc Cuong Nguyen, Anthony T. Patera, George S.H. Pau To cite this version: Yvon aday, Ngoc Cuong Nguyen, Anthony T. Patera,
More informationConvergence Rates for Greedy Algorithms in Reduced Basis Methods
Convergence Rates for Greedy Algorithms in Reduced Basis Methods Peter Binev, Albert Cohen, Wolfgang Dahmen, Ronald DeVore, Guergana Petrova, and Przemyslaw Wojtaszczyk October 4, 2010 Abstract The reduced
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 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 informationPOD-DEIM APPROACH ON DIMENSION REDUCTION OF A MULTI-SPECIES HOST-PARASITOID SYSTEM
POD-DEIM APPROACH ON DIMENSION REDUCTION OF A MULTI-SPECIES HOST-PARASITOID SYSTEM G. Dimitriu Razvan Stefanescu Ionel M. Navon Abstract In this study, we implement the DEIM algorithm (Discrete Empirical
More informationA-Posteriori Error Estimation for Parameterized Kernel-Based Systems
A-Posteriori Error Estimation for Parameterized Kernel-Based Systems Daniel Wirtz Bernard Haasdonk Institute for Applied Analysis and Numerical Simulation University of Stuttgart 70569 Stuttgart, Pfaffenwaldring
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
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 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 informationProper Orthogonal Decomposition (POD)
Intro Results Problem etras Proper Orthogonal Decomposition () Advisor: Dr. Sorensen CAAM699 Department of Computational and Applied Mathematics Rice University September 5, 28 Outline Intro Results Problem
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationarxiv: v2 [math.na] 15 May 2014
A nonintrusive Reduced Basis Method applied to aeroacoustic simulations arxiv:1401.3684v2 [math.na] 15 May 2014 March 26, 2018 Fabien Casenave 1, Alexandre Ern 1, and Tony Lelièvre 1,2 1 Université Paris-Est,
More informationPOD/DEIM 4DVAR Data Assimilation of the Shallow Water Equation Model
nonlinear 4DVAR 4DVAR Data Assimilation of the Shallow Water Equation Model R. Ştefănescu and Ionel M. Department of Scientific Computing Florida State University Tallahassee, Florida May 23, 2013 (Florida
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 informationPOD-based Reduced Order Models for Parameterized PDEs
POD-based Reduced Order Models for Parameterized PDEs Florian De Vuyst florian.de-vuyst@ecp.fr Ecole Centrale Paris Laboratoire Mathématiques Appliquées aux Systèmes July 2008 the 1st Outline 1. Setting
More informationA reduced basis approach for modelling the movement of nuclear reactor control rods
Alberto Sartori Nuclear Engineering Division Department of Energy Politecnico di Milano Milano, 20156, Italy Email: alberto.sartori@polimi.it Antonio Cammi Nuclear Engineering Division Department of Energy
More informationThe Certified Reduced-Basis Method for Darcy Flows in Porous Media
The Certified Reduced-Basis Method for Darcy Flows in Porous Media S. Boyaval (1), G. Enchéry (2), R. Sanchez (2), Q. H. Tran (2) (1) Laboratoire Saint-Venant (ENPC - EDF R&D - CEREMA) & Matherials (INRIA),
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationStability of the Parareal Algorithm
Stability of the Parareal Algorithm Gunnar Andreas Staff and Einar M. Rønquist Norwegian University of Science and Technology Department of Mathematical Sciences Summary. We discuss the stability of the
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 reduced basis method and ROM-based optimization for batch chromatography
MAX PLANCK INSTITUTE MedRed, 2013 Magdeburg, Dec 11-13 A reduced basis method and ROM-based optimization for batch chromatography Yongjin Zhang joint work with Peter Benner, Lihong Feng and Suzhou Li Max
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 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 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 informationInverse problems in tumor growth modelling
MC2 Damiano Lombardi Angelo Iollo Thierry Colin Olivier Saut Marseille, 18-09-09 Outline 1) Introduction: phenomenological model for tumor growth models 2) Models based on mixture theory: a simplified
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 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 informationREDUCED ORDER METHODS
REDUCED ORDER METHODS Elisa SCHENONE Legato Team Computational Mechanics Elisa SCHENONE RUES Seminar April, 6th 5 team at and of Numerical simulation of biological flows Respiration modelling Blood flow
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 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 informationApplication of POD-DEIM Approach on Dimension Reduction of a Diffusive Predator-Prey System with Allee effect
Application of POD-DEIM Approach on Dimension Reduction of a Diffusive Predator-Prey System with Allee effect Gabriel Dimitriu 1, Ionel M. Navon 2 and Răzvan Ştefănescu 2 1 The Grigore T. Popa University
More informationRecent advances on nonlinear Reduced Order Modelling for stability and bifurcations problems in incompressible fluid dynamics
Recent advances on nonlinear Reduced Order Modelling for stability and bifurcations problems in incompressible fluid dynamics Giuseppe Pitton, Gianluigi Rozza SISSA, International School for Advanced Studies,
More informationOptimal rotary control of the cylinder wake using POD reduced order model p. 1/30
Optimal rotary control of the cylinder wake using POD reduced order model Michel Bergmann, Laurent Cordier & Jean-Pierre Brancher Michel.Bergmann@ensem.inpl-nancy.fr Laboratoire d Énergétique et de Mécanique
More informationOn a Data Assimilation Method coupling Kalman Filtering, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model
On a Data Assimilation Method coupling, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model 2016 SIAM Conference on Uncertainty Quantification Basile Marchand 1, Ludovic
More informationSplitting Techniques in the Face of Huge Problem Sizes: Block-Coordinate and Block-Iterative Approaches
Splitting Techniques in the Face of Huge Problem Sizes: Block-Coordinate and Block-Iterative Approaches Patrick L. Combettes joint work with J.-C. Pesquet) Laboratoire Jacques-Louis Lions Faculté de Mathématiques
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More informationMIT 9.520/6.860, Fall 2017 Statistical Learning Theory and Applications. Class 19: Data Representation by Design
MIT 9.520/6.860, Fall 2017 Statistical Learning Theory and Applications Class 19: Data Representation by Design What is data representation? Let X be a data-space X M (M) F (M) X A data representation
More informationTowards Reduced Order Modeling (ROM) for Gust Simulations
Towards Reduced Order Modeling (ROM) for Gust Simulations S. Görtz, M. Ripepi DLR, Institute of Aerodynamics and Flow Technology, Braunschweig, Germany Deutscher Luft und Raumfahrtkongress 2017 5. 7. September
More informationApplication of POD-DEIM Approach for Dimension Reduction of a Diffusive Predator-Prey System with Allee Effect
Application of POD-DEIM Approach for Dimension Reduction of a Diffusive Predator-Prey System with Allee Effect Gabriel Dimitriu 1, Ionel M. Navon 2 and Răzvan Ştefănescu 2 1 The Grigore T. Popa University
More informationOn Riesz-Fischer sequences and lower frame bounds
On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition
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 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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationMath 108b: Notes on the Spectral Theorem
Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator
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 informationReal Analysis Qualifying Exam May 14th 2016
Real Analysis Qualifying Exam May 4th 26 Solve 8 out of 2 problems. () Prove the Banach contraction principle: Let T be a mapping from a complete metric space X into itself such that d(tx,ty) apple qd(x,
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
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 informationProper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems. Stefan Volkwein
Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems Institute for Mathematics and Scientific Computing, Austria DISC Summerschool 5 Outline of the talk POD and singular value decomposition
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationA NEW ITERATIVE METHOD FOR THE SPLIT COMMON FIXED POINT PROBLEM IN HILBERT SPACES. Fenghui Wang
A NEW ITERATIVE METHOD FOR THE SPLIT COMMON FIXED POINT PROBLEM IN HILBERT SPACES Fenghui Wang Department of Mathematics, Luoyang Normal University, Luoyang 470, P.R. China E-mail: wfenghui@63.com ABSTRACT.
More informationarxiv: v1 [math.na] 10 Oct 2017
submitted to Journal of Computational Physics, October 2017 A Reduced Basis Technique for Long-Time Unsteady Turbulent Flows Lambert Fick 1, Yvon Maday 2,3, Anthony T Patera 4, Tommaso Taddei 2 arxiv:1710.03569v1
More informationJOINT SOURCE AND SENSOR PLACEMENT FOR SOUND FIELD CONTROL BASED ON EMPIRICAL INTERPOLATION METHOD Hongo, Bunkyo-ku, Tokyo , Japan
JOINT SOURCE AND SENSOR PLACEMENT FOR SOUND FIELD CONTROL BASED ON EMPIRICAL INTERPOLATION METHOD Shoichi Koyama 1,2, Gilles Chardon 3, and Laurent Daudet 2 1 The University of Tokyo, Graduate School of
More informationGreedy bases are best for m-term approximation
Greedy bases are best for m-term approximation Witold Bednorz Institute of Mathematics, Warsaw University Abstract We study the approximation of a subset K in a Banach space X by choosing first basis B
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 informationSubmitted to SISC; revised August A SPACE-TIME PETROV-GALERKIN CERTIFIED REDUCED BASIS METHOD: APPLICATION TO THE BOUSSINESQ EQUATIONS
Submitted to SISC; revised August 2013. A SPACE-TIME PETROV-GALERKI CERTIFIED REDUCED BASIS METHOD: APPLICATIO TO THE BOUSSIESQ EQUATIOS MASAYUKI YAO Abstract. We present a space-time certified reduced
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 informationReduced basis method applied to a convective stability problem
Herrero et al. Journal of Mathematics in Industry (208) 8: https://doi.org/0.86/s3362-08-0043-6 R E S E A R C H Open Access Reduced basis method applied to a convective stability problem Henar Herrero,
More informationLow Rank Approximation Lecture 3. Daniel Kressner Chair for Numerical Algorithms and HPC Institute of Mathematics, EPFL
Low Rank Approximation Lecture 3 Daniel Kressner Chair for Numerical Algorithms and HPC Institute of Mathematics, EPFL daniel.kressner@epfl.ch 1 Sampling based approximation Aim: Obtain rank-r approximation
More informationStochastic methods for solving partial differential equations in high dimension
Stochastic methods for solving partial differential equations in high dimension Marie Billaud-Friess Joint work with : A. Macherey, A. Nouy & C. Prieur marie.billaud-friess@ec-nantes.fr Centrale Nantes,
More informationImprovement of Reduced Order Modeling based on Proper Orthogonal Decomposition
ICCFD5, Seoul, Korea, July 7-11, 28 p. 1 Improvement of Reduced Order Modeling based on Proper Orthogonal Decomposition Michel Bergmann, Charles-Henri Bruneau & Angelo Iollo Michel.Bergmann@inria.fr http://www.math.u-bordeaux.fr/
More informationI = i 0,
Special Types of Matrices Certain matrices, such as the identity matrix 0 0 0 0 0 0 I = 0 0 0, 0 0 0 have a special shape, which endows the matrix with helpful properties The identity matrix is an example
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 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 informationA Reduced Basis Method for Evolution Schemes with Parameter-Dependent Explicit Operators
B. Haasdonk a M. Ohlberger b G. Rozza c A Reduced Basis Method for Evolution Schemes with Parameter-Dependent Explicit Operators Stuttgart, March 29 a Instite of Applied Analysis and Numerical Simulation,
More informationLeast-Squares Finite Element Methods
Pavel В. Bochev Max D. Gunzburger Least-Squares Finite Element Methods Spri ringer Contents Part I Survey of Variational Principles and Associated Finite Element Methods 1 Classical Variational Methods
More informationMultivariate Predictions of Local Parametric Reduced-Order Models Characteristics and Their Applications
Multivariate Predictions of Local Parametric Reduced-Order Models Characteristics and Their Applications Azam Moosavi, a Razvan Stefanescu b1 and Adrian Sandu a a Computational Science Laboratory, Department
More informationAn Empirical Chaos Expansion Method for Uncertainty Quantification
An Empirical Chaos Expansion Method for Uncertainty Quantification Melvin Leok and Gautam Wilkins Abstract. Uncertainty quantification seeks to provide a quantitative means to understand complex systems
More informationSeptember 22, :50 WSPC/INSTRUCTION FILE KNP Boussinesq revised extended
Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company REDUCED BASIS APPROXIMATIO AD A POSTERIORI ERROR ESTIMATIO FOR THE PARAMETRIZED USTEADY BOUSSIESQ EQUATIOS DAVID
More informationA perturbation-based reduced basis method for parametric eigenvalue problems
A perturbation-based reduced basis method for parametric eigenvalue problems Eric Cancès 1, Virginie Ehrlacher 1, David Gontier 2, Damiano Lombardi 3 and Antoine Levitt 1 1. CERMICS, Ecole des Ponts ParisTech
More information