Comparison of the Reduced-Basis and POD a-posteriori Error Estimators for an Elliptic Linear-Quadratic Optimal Control
|
|
- Martin Burns
- 5 years ago
- Views:
Transcription
1 SpezialForschungsBereich F 3 Karl Franzens Universität Graz Technische Universität Graz Medizinische Universität Graz Comparison of the Reduced-Basis and POD a-posteriori Error Estimators for an Elliptic Linear-Quadratic Optimal Control T. Tonn K. Urban S. Volkwein SFB-Report No. 3 April A 8 GRAZ, HEINRICHSTRASSE 36, AUSTRIA Supported by the Austrian Science Fund (FWF)
2 SFB sponsors: Austrian Science Fund (FWF) University of Graz Graz University of Technology Medical University of Graz Government of Styria City of Graz
3 Mathematical and Computer Modelling of Dynamical Systems Vol., No., March, 3 Article Comparison of the Reduced-Basis and POD a-posteriori Error Estimators for an Elliptic Linear-Quadratic Optimal Control Problem Timo Tonn a, K. Urban a and S. Volkwein b a Institut für Numerische Mathematik, Universität Ulm, Helmholtzstrasse 8, D-8969 Ulm, Germany; b Fachbereich Mathematik und Statistik, Universität Konstanz, Universitätsstraße, D Konstanz, Germany (February ) In this article, a linear-quadratic optimal control problem governed by the Helmholtz equation is considered. For the computation of suboptimal solutions, two different model reduction techniques are compared: the reduced-basis method (RBM) and proper orthogonal decomposition (POD). By an a-posteriori error estimator for the optimal control problem the accuracy of the suboptimal solutions is ensured. The efficiency of both model reduction approaches is illustrated by a numerical example. Keywords: optimal control, a-posteriori error estimate, proper orthogonal decomposition, reduced-basis method AMS Subject Classification: 49K, 65K, 9C, 9C46. Introduction Optimal control problems for partial differential equations are often hard to tackle numerically because their discretization leads to very large scale optimization problems. Therefore, different techniques of model reduction have been developed to approximate these problems by smaller ones that are tractable with less effort. Among them, the reduced-basis (RB) and the proper orthogonal decomposition (POD) method are widely used, also in the context of nonlinear problems. Some reduced order methods like balanced truncation offer a reliable a-priori error analysis for optimal control applications, [, ]. However, for POD and RB it is not a-priori clear how far the optimal control of the reduced-order problem is from the exact one, unless its snapshots are generating a sufficiently rich space, where sufficiently rich implies that the space contains all possible snapshots. However, we are able to compensate for the lack of a-priori analysis for the POD and RB method by utilizing an a-posteriori analysis. This approach is based on a fairly standard perturbation argument to deduce how far the suboptimal control, computed on the basis of the reduced-order method, is away from the (unknown) exact one. Corresponding author. Stefan.Volkwein@uni-konstanz.de ISSN: print/issn online c Taylor & Francis DOI:.8/387395YYxxxxxxxx
4 T. Tonn, K. Urban, and S. Volkwein Domain Ω and edgelabels Impedance Melamin mm y axis Re(Z) Im(Z) x axis frequency range in Hz Figure. Two-dimensional cross section of an idealized interior Ω of the vehicle, where the boundary part Γ R consists of parts 4 and 5 ; impedance values Z = Z R + jz I for melamine mm width in the frequency range from to Hz.. The optimal control problem and optimality conditions In this section, we introduce the linear-quadratic elliptic optimal control problem, recall the associated first-order necessary optimality conditions and define a reduced problem only for the control variable. The optimal control problem Suppose that the interior of a car is simplified by the two-dimensional domain Ω plotted in Figure. The boundary Γ = Ω is split into two measurable disjunct parts Γ R and Γ N. Denoting the frequency by f, let Z f C (see Figure ) be a given complex impedance. Then, the associated sound pressure p : Ω C is governed by the Helmholtz equation p(x) k f p(x) = u b(x) for all x = (x, x ) Ω, together with the boundary conditions j p(x) ϱ ω f n = p(x) for all x Γ R, Z f j ϱ ω f p(x) n = for all x Γ N. (c) In (), ϱ =.9985 [ kg ] m is an ambient density, ωf = πf is the circle frequency 3 and k f = ω f /c is the wave number, where the constant c = [ ] m s denotes the speed of sound. The right-hand side is a simplified model for a source at the point x q = (.,.8) (e.g., a loudspeaker located at x q ) with the intensity u, u C, and shape function b(x) = exp( ( x x q ) for x = (x, x ) Ω. For the normal impedance boundary condition, let j be the imaginary unit and n denote the derivative in the outward normal direction. All other parts on the boundary are assumed to be perfectly rigid, see (c). Throughout the paper, we suppose the following assumption.
5 Mathematical and Computer Modelling of Dynamical Systems 3.4 Domain Ω, measurement points and source. x axis x axis Figure. Measurement points ( ) and location of the source ( ) for our numerical examples. Assumption : For all given frequencies f F := [ Hz, Hz], all given u C and the impedance function Z f, plotted in Figure, the problem () admits a unique solution p = p f (u). Remark : Due to the Fredholm theory, [3], we can ensure existence of a solution provided kf is not an eigenvalue of considered on Ω with Neumann and Robin boundary conditions on Γ N and Γ R, respectively. Let p i f C, i =,..., n m, n m N, be measured sound pressures at different given observation points x i Ω Γ N, i n m ; see Figure for given frequency f F. We introduce the quadratic cost functional J f (p, u) := α n m i= p(xi ) p i f + σ u u f, where u f C is a given reference intensity, z is the absolute modulus of a complex number z, α is non-negative and σ is positive. Then, given f, we consider the optimal control problem min J f (x f ) subject to (s.t.) x f = (p f, u f ) solves (). (P f ) Thus, (P f ) is a linear-quadratic optimal control problem for any frequency f F. First-order optimality conditions The first-order necessary optimality conditions to (P f ) consist of the state equation, the adjoint system for the Lagrange multiplier λ f n m λ f (x) k f λ f (x) = α (p i f p f (x i))δ xi (x) for all x Ω, i= j λ f (x) = λ f (x) for all x Γ R, ϱ ω f n Z f j λ f (x) = for all x Γ N ϱ ω f n
6 4 T. Tonn, K. Urban, and S. Volkwein and the equation for the optimal control u f σ ( u f ) u f b(x)λ (x) dx = in C, Ω where δ xi denotes the Dirac delta distribution satisfying δ xi, ϕ = ϕ(x i ) for ϕ C(Ω Γ N ) and i =,..., n m. As usual z = x jy is the complex conjugate of z = x + jy C. For more details we refer the reader to [4, 5]. The reduced problem Motivated by Assumption we introduce the reduced cost functional Ĵ : C [, ) by Ĵ f (u) = J f (p f (u), u) for u C where p f (u) again denotes the unique solution to () for a given control input u at a given frequency f. Then, (P f ) is equivalent to the reduced optimal control problem min Ĵf (u f ) s.t. u f C. ( ˆP f ) If x f = (p f, u f ) is an optimal solution to (P f ), then u f solves ( ˆP f ). On the other hand, if u f is an optimal solution to ( ˆP f ), then the pair (p f (u f ), u f ) solves (P f ). The ultimate goal is to determine the optimal control for many values f F. In order to avoid a naive and possibly costly evaluation of the elliptic linear-quadratic optimal control problem for various values of f we introduce a model reduction. This leads to a reduced-order model for ( ˆP f ). The accuracy of the reduced-order model is controlled by an a-posteriori error analysis. This is the focus of the next section. 3. A-posteriori error estimate for the optimal control problem In this section we recall some results of the a-posteriori error estimator. For more details, we refer the reader to [6, 7]. Suppose that u C is an arbitrary control input. Then, the difference u f u can be estimated without requiring knowledge of u f. Theorem 3. : Let u f be an optimal solution to ( ˆP f ). Let p f and λ f be the associated state variable and Lagrange multiplier, respectively. Suppose that u C
7 Mathematical and Computer Modelling of Dynamical Systems 5 is chosen arbitrarily, p f = p f (u) solves () and λ f = λ f (p f ) solves λ f (x) + k f λ f (x) = α n m i= (p i f p f (x i ))δ xi (x), for all x Ω, j λ f (x) = λ f (x), for all x Γ R, ϱ ω f n Z f j λ f (x) =, for all x Γ N. ϱ ω f n () Then, it follows that where ζ f = ζ f (u, λ f ) = σ ( u u ) f Ω b(x)λ f (x) dx. u f u ζf, (3) σ Proof : The claim follows from Theorem 3., Proposition 3., and Remark 3.3 in [6]. We will call the right-hand side of (3) an a-posteriori error estimate, since, in the next two sections, we shall apply it to suboptimal controls u = û f that have been computed from a reduced-order model (utilizing the reduced-basis respectively the POD method). After having computed û f, we determine the associated state ˆp f = p f (û f ) and adjoint state ˆλ f = λ f (ˆp f ). Then, we can compute ˆζ f = ζ f (û f, ˆλ f ) and its norm, s.t. (3) gives an upper bound for the distance of û f to u f. In this way, the error u f û f caused by the reduced-order model can be estimated a-posteriori. If the error is too large, we have to include more basis functions into the reduced model. 4. Reduced-order methods In this section, we introduce the reduced-order modeling for () utilizing POD and RBM. Let us emphasize that the computation of the reduced-order model is only based on the state equation (), so that no information from the optimal solution of (P f ) is needed. 4.. Proper orthogonal decomposition In this subsection, we briefly review the POD method for our problem and derive the reduced-order model. Let X be either the Lebesgue space L (Ω) supplied with the inner product ψ, ϕ L(Ω) = Ω ψϕ dx for ψ, ϕ L (Ω) or the Sobolev space H (Ω) endowed with the common inner product ψ, ϕ H (Ω) = Ω ( ψ ϕ + ψϕ) dx for ψ, ϕ H (Ω).
8 6 T. Tonn, K. Urban, and S. Volkwein Normalized eigenvalues real part imaginary part Normalized eigenvalues real part imaginary part Figure 3. Decay of the first 6 POD eigenvalues for X = L (Ω) and X = H (Ω). Note that () is a linear elliptic equation and u is a complex number. Thus, by linear superposition for any frequency f, the solution to () is given by a linear combination of p and p, where p solves () for u = and p is the solution to () for u = j. Let Ξ f := {f,..., f n } F be a given snapshot grid. We denote by p j = R(p (f j )) and p j+n = R(p (f j )) the real part of the solution p to () for u = and u = j, respectively, at frequency instances f j, j =,..., n. Let V R = span {p,..., p n } X with d R := dim V R n. Then, for given l n we consider the minimization problem n l min p j p j, ψ i ψ,...,ψ l X X ψ i X j= i= s.t. ψ i, ψ j X = δ ij for i, j l. Remark : In some applications, the mean value of the snapshots is included in the POD modeling (see, e.g., in [8]). In our numerical experiments it turns out that this approach does not give rise to better results. The solution to (4) is given by the solution of the eigenvalue problem Rψ i = n j= p j, ψ i X p j = λ i ψ i for i =,..., l, where R : X V R X is a linear, bounded, compact, self-adjoint and nonnegative operator; see [8]. Thus, there exists an orthonormal set {ψ i } dr i= of eigenfunctions and corresponding non-negative eigenvalues {λ i } dr i= satisfying Rψ i = λ i ψ i, λ λ λ d R. In [9] the dependence of {ψ i } l i= and {λ i} l i= for l n on the chosen snapshot grid {f j } n j= is investigated. Analogously, a POD basis {φ i } l i= for the snapshot space V I = span {I(p, ),..., I(p,n ), I(p, ),..., I(p,n )} X can be introduced. Since the eigenvalues for the real and for the imaginary parts decay in a similar rate (see (4) For z = x + jy C we use the notation R(z) := x and I(z) := y.
9 Mathematical and Computer Modelling of Dynamical Systems 7 Figure 3), we choose the same number of ansatz functions for the real and for the imaginary parts. In order to highlight similarities and differences of POD and RB method we stress, that the modes are linear combinations of all snapshots p j, j =,..., n. 4.. The reduced-basis method Before going into detail on the application of the reduced basis method, we reconsider the snapshot generation of the POD from above. Recall that Ξ f denotes the frequency snapshot grid which we now extend by the control, i.e. we introduce f = (f, u) F C and Ξ f := { f = (f, u) : f Ξ f, u {, j}}. Then, for POD the snapshots are the solutions p f = p f (u) of () for all f Ξ f. In contrast to this, for the RBM, the Helmholtz equation () has to be solved for only some particular values in Ξ f, where the selection procedure is stated below. As a start, let S l f := { f,..., f l } Ξ f denote the set of identified values, VR l := span {ψ i := R(p f ), i l} X and i V l I := span {φ i := I(p f i ), i l} X. Note, that for algebraic stability of the reduced order model, the bases {ψ i } l i= and {φ i} l i= are orthonormalized w.r.t. X by a Gram-Schmidt procedure. Moreover, the modes are a linear combination of just a few snapshots. Next, in order to identify frequencies in S l f, which lead to good resulting bases {ψ i } l i= and {φ i} l i=, we utilize the standard greedy procedure (c.p. []) using an a-posteriori error estimator l f to be explained later, which reads: Algorithm : : Set l max N, ε > and choose f Ξ f arbitrarily. : for l = to l max do 3: Compute ε l := max f Ξ f l f / ˆp l f X. 4: if ε l < ε or l = l max then 5: break. 6: else 7: Set f l+ := arg max f Ξ l f f / ˆpl f X. 8: Compute p f. l+ 9: Set V l+ R = VR l {R(p f )} and V l+ l+ I = VI l {I(p f )}. l+ : end if : end for Note, that ˆp l f in Algorithm denotes the reduced order solution of () using VR l jv I l as trial- and testspaces. It is well known, that l f := β f R l f X. (5) is a rigorous a-posteriori error estimator for p f ˆp l f X, where X := R( ) X + I( ) X, β(f) is the inf-sup constant of () and R l f is the residual of () for the reduced order solution ˆp l f. Obviously, for the application in Algorithm, l f should be rapidly evaluable.
10 8 T. Tonn, K. Urban, and S. Volkwein Figure 4. Inf-sup constant its approximation ˆβ f as well as f i identified by Algorithm. While this holds for the dual norm of the residual, which can be evaluated in O(l ) (c.p. []), the evaluation of is more involved than solving the actual problem (). A possible way out is the application of the so called Successive Constraint Method (SCM, c.p. []) to obtain a lower bound to that is rapidly evaluable, but at the cost of an expensive (offline) precomputation. On the other hand, in this work, we are interested in û l f, the suboptimal solution to ( ˆP f ), for which we already have an a-posteriori error estimator (c.p. Theorem 3.), s.t. l f does not necessarily need to be rigorous, rather than a good indicator. Hence, the easiest way to avoid the cost of an expensive (offline) precomputation is to use ˆβ f =. Another possibility is to compute for some frequencies (say n), only, and use the part of the SCM that determines ˆ, an upper bound to, in O(n). Figure 4 visualizes as well as ˆ obtained from computing for equidistant f i, i n =. Moreover, it shows the by Algorithm identified values f i = (fi, u i ), where f i is marked on the bottom (top) axis if u i = (u i = j). We will investigate the differences resulting from using ˆ = and ˆβ f = as approximations to in Section 5. The CPU time needed for computing ˆβ f is reported in Table. Remark : Note, that for the POD method the bases {ψ i } l i= and {φ i} l i= are computed, s.t. the mean error to the snapshots is minimized (c.p. (4)), whereas for the RBM {ψ i } l i= and {φ i} l i= are computed, s.t. the maximum error to the snapshots is minimized (c.p. Algorithm ). The a-posteriori error estimator for (ˆP f ) The a-posteriori error estimator for ( ˆP f ) is formulated in Algorithm. Algorithm : : Choose l =, l max = 6, and ε. Set flag =. : for f = to do 3: 4: repeat 5: Set l = max(, l ). 6: Calculate the suboptimal solution û l f to ( ˆP f ).
11 Mathematical and Computer Modelling of Dynamical Systems Figure 5. Run : Number l of POD ansatz functions over the frequency band F for X = L (Ω) and for X = H (Ω) with the tolerance ε = 4 for the a-posteriori error estimator. 7: Compute the solutions ˆp l f to () for u = ûl f and ˆλ l f to () for p f = ˆp l f. 8: Set ˆζ f l = σ(ûl f u (f)) Ω b(x)ˆλ l f (x) dx. 9: if ˆζ f l < ε then : Set flag =. : Return flag, l, ū l f and STOP. : else 3: Set l = l +. 4: end if 5: until l > l max 6: if flag = then 7: Increase l max and restart the algorithm. 8: break. 9: end if : end for 5. Numerical experiments In this section, we present numerical experiments. One example is constructed in such a way that the optimal control is known. Thus, we can verify quality of the estimate (3). In the second example, the (exact) optimal control is unknown. All coding is done in Matlab using routines from the PDE Toolbox with finite elements (FE). We apply a standard piecewise linear FE discretization with m = 4957 degrees of freedom. Run : Let the impedance Z f be given for the material melamine with mm depth; see Figure. We choose for the cost functional α =., σ = and ( ( π(f ) ) ( π(f ) ) ) u f = cos + j sin. The measurement points x i, i =,..., n m with n m = are depicted in Figure. The values p i f, i =,...,, are given by pi f = p f (u f )(x i), i.e., by the solution to () evaluated at x i for i =,..., and f F. Then, the optimal solution to ( ˆP f ) is u f = u (f) for f F. First we apply Algorithm for the tolerance ε = 4. In Figure 5, we present the number of l of utilized POD basis functions for the choices X = L (Ω) and X = H (Ω). It turns out that the required tolerance is achieved
12 T. Tonn, K. Urban, and S. Volkwein Figure 6. Run : Number l of RB ansatz functions over the frequency band F for X = L (Ω) and for X = H (Ω) with the tolerance ε = 4 for the a-posteriori error estimator. for all f F with l basis functions. We observe that the number of POD basis functions changes over the frequency band F. Moreover, the different choices for X also influences the value of l. The results for the reduced-basis method are given in Figure 6. Compared to the POD approximations, the RB methods needs more basis functions at the beginning, but (significantly) less for f (3, 3]. The needed CPU times are presented in Table. Note, that for RBM Table contains two columns, where the first one corresponds to ˆ and the second one to ˆ for approximating. Table. Run : CPU times in seconds. Snapshot generation with FE.9 POD computation.4 POD optimization for X = L (Ω), ε = POD optimization for X = H (Ω), ε = POD optimization for X = L (Ω), ε = POD optimization for X = H (Ω), ε = Approximation of inf-sup constant RBM computation RBM optimization for X = L (Ω), ε = RBM optimization for X = H (Ω), ε = RBM optimization for X = L (Ω), ε = RBM optimization for X = H (Ω), ε = Next, we apply Algorithm again with the smaller tolerance ε = 6. In Figure 7, the number l of utilized POD basis functions is plotted. The results for the reduced-basis method are given in Figure 8. Again, the required tolerance is achieved for all variants of the POD and RB method. As expected, more basis functions have to be included in our reduced-order modeling. Recall, that u f denotes the optimal solution to ( ˆP f ) and û l f C denotes the suboptimal solution obtained from using l basis functions. The decay of max f Ξf u f ûl f is visualized in Figure 9 for X = L (Ω) using POD and RBM. Note, that we do not present the effectivity, i.e., the ratio between the true error u f ûl f and its estimator ˆζ f l /σ, as it turns out to be bounded by. for all f Ξ f. Run : In the second example, we choose the same parameter as in the first run, but now u f =. Thus, the optimal solution to ( ˆP f ) is not known. We apply Algorithm with the tolerance ε = 6. In Figure, the number l of utilized POD basis functions is plotted. The results for the reduced-basis method are given in Figure.
13 Mathematical and Computer Modelling of Dynamical Systems Figure 7. Run : Number l of POD ansatz functions over the frequency band F for X = L (Ω) and for X = H (Ω) with the tolerance ε = 6 for the a-posteriori error estimator Figure 8. Run : Number l of RB ansatz functions over the frequency band F for X = L (Ω) and for X = H (Ω) with the tolerance ε = 6 for the a-posteriori error estimator. POD RB, RB, l Figure 9. Run : Decay of max f Ξf u f ûl f for POD and RBM in X = L (Ω). 6. Conclusion We have compared RBM and POD for a linear-quadratic optimal control problem which is constrained by the Helmholtz equation. Thus, the following conclusions hold for this particular problem only. Moreover, it should be noted that the control u C (and thus the parameter for RBM and POD, respectively) is only formally two-dimensional (real and imaginary part) since we have seen that linear superposition allows to reduce the parameter dimension. We have demonstrated that both methods are quite efficient, the differences are
14 T. Tonn, K. Urban, and S. Volkwein Figure. Run : Number l of POD ansatz functions over the frequency band F for X = L (Ω) and for X = H (Ω) with the tolerance ε = 6 for the a-posteriori error estimator Figure. Run : Number l of RB ansatz functions over the frequency band F for X = L (Ω) and for X = H (Ω) with the tolerance ε = 4 for the a-posteriori error estimator. more or less marginal. As we can see in Table, RBM is a little more efficient than POD as long as we can avoid the computation of an approximation for the inf sup constant. Also the average number of RBM modes is smaller than than the one for POD as we can deduce by the comparison of Figures 5,6 and Figures 7,8 for Run, Figures, for Run as well as the CPU timings for Algorithm in Table. On the other hand, Figure 9 shows that the true error in L (Ω) over the frequency sample Ξ f is smaller for POD. This is somehow astonishing since the POD is known to be the best approximation of the state whereas here we consider the control. It also seems that POD earlier reaches the asymptotic regime. Of course, these conclusions have to verified for different examples, in particular also for higher parameter dimensions.
15 References REFERENCES 3 [] A.C. Antoulas Approximation of Large-Scale Dynamical Systems (Advances in Design and Control) (Advances in Design and Control), Society for Industrial and Applied Mathematics, Philadelphia, PA, USA, 5. [] P. Benner and E.S. Quintana-ortí, Model reduction based on spectral projection methods, in Dimension Reduction of Large-Scale Systems Springer-Verlag, 5, pp [3] S. Reed and B. Simon Methods of Modern Mathematical Physics, I. Functional Analysis, second Academic Press, San Diego, 98. [4] S. Volkwein, Admittance identification from point-wise sound pressure measurements using reducedorder modelling, (submitted 9),. [5] S. Volkwein and A. Hepberger, Impedance identification by POD model reduction techniques, at Automatisierungstechnik 8 (8), pp [6] F. Tröltzsch and S. Volkwein, POD a-posteriori error estimates for linear-quadratic optimal control problems, Computational Optimization and Applications 44 (9), pp [7] M. Kahlbacher and S. Volkwein, POD a-posteriori error based inexact SQP method for bilinear elliptic optimal control problems, (submitted ),. [8] P. Holmes, J. Lumley, and G. Berkooz Turbulence, Coherent Structures, Dynamical Systems and Symmetry, Cambridge University Press, Cambridge, 996. [9] K. Kunisch and S. Volkwein, Galerkin Proper Orthogonal Decomposition Methods for a General Equation in Fluid Dynamics, SIAM J. Numer. Anal. (), pp [] N.C. Nguyen, Reduced-Basis Approximations and A Posteriori Error Bounds for Nonaffine and Nonlinear Partial Differential Equations: Application to Inverse Analysis, Singapore-MIT Alliance, 5. [] D. Huynh, G. Rozza, S. Sen, and A. Patera, A successive constraint linear optimization method for lower bounds of parametric coercivity and inf sup stability constants, Comptes Rendus Mathematique 345 (7), pp
Reduced-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 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 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 informationPOD for Parametric PDEs and for Optimality Systems
POD for Parametric PDEs and for Optimality Systems M. Kahlbacher, K. Kunisch, H. Müller and S. Volkwein Institute for Mathematics and Scientific Computing University of Graz, Austria DMV-Jahrestagung 26,
More informationEileen Kammann Fredi Tröltzsch Stefan Volkwein
Universität Konstanz A method of a-posteriori error estimation with application to proper orthogonal decomposition Eileen Kammann Fredi Tröltzsch Stefan Volkwein Konstanzer Schriften in Mathematik Nr.
More informationModel order reduction & PDE constrained optimization
1 Model order reduction & PDE constrained optimization Fachbereich Mathematik Optimierung und Approximation, Universität Hamburg Leuven, January 8, 01 Names Afanasiev, Attwell, Burns, Cordier,Fahl, Farrell,
More informationFINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS
Numerical Functional Analysis and Optimization, 28(7 8):957 973, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0163-0563 print/1532-2467 online DOI: 10.1080/01630560701493305 FINITE ELEMENT APPROXIMATION
More informationAn optimal control problem for a parabolic PDE with control constraints
An optimal control problem for a parabolic PDE with control constraints PhD Summer School on Reduced Basis Methods, Ulm Martin Gubisch University of Konstanz October 7 Martin Gubisch (University of Konstanz)
More informationProper Orthogonal Decomposition. POD for PDE Constrained Optimization. Stefan Volkwein
Proper Orthogonal Decomposition for PDE Constrained Optimization Institute of Mathematics and Statistics, University of Constance Joined work with F. Diwoky, M. Hinze, D. Hömberg, M. Kahlbacher, E. Kammann,
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 informationA Posteriori Error Bounds for Meshless Methods
A Posteriori Error Bounds for Meshless Methods Abstract R. Schaback, Göttingen 1 We show how to provide safe a posteriori error bounds for numerical solutions of well-posed operator equations using kernel
More informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More informationINTERPRETATION OF PROPER ORTHOGONAL DECOMPOSITION AS SINGULAR VALUE DECOMPOSITION AND HJB-BASED FEEDBACK DESIGN PAPER MSC-506
INTERPRETATION OF PROPER ORTHOGONAL DECOMPOSITION AS SINGULAR VALUE DECOMPOSITION AND HJB-BASED FEEDBACK DESIGN PAPER MSC-506 SIXTEENTH INTERNATIONAL SYMPOSIUM ON MATHEMATICAL THEORY OF NETWORKS AND SYSTEMS
More informationProper Orthogonal Decomposition for Optimal Control Problems with Mixed Control-State Constraints
Proper Orthogonal Decomposition for Optimal Control Problems with Mixed Control-State Constraints Technische Universität Berlin Martin Gubisch, Stefan Volkwein University of Konstanz March, 3 Martin Gubisch,
More informationContributors and Resources
for Introduction to Numerical Methods for d-bar Problems Jennifer Mueller, presented by Andreas Stahel Colorado State University Bern University of Applied Sciences, Switzerland Lexington, Kentucky, May
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationNumerical approximation for optimal control problems via MPC and HJB. Giulia Fabrini
Numerical approximation for optimal control problems via MPC and HJB Giulia Fabrini Konstanz Women In Mathematics 15 May, 2018 G. Fabrini (University of Konstanz) Numerical approximation for OCP 1 / 33
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationConvergence of a finite element approximation to a state constrained elliptic control problem
Als Manuskript gedruckt Technische Universität Dresden Herausgeber: Der Rektor Convergence of a finite element approximation to a state constrained elliptic control problem Klaus Deckelnick & Michael Hinze
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 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 informationAn Introduction to Variational Inequalities
An Introduction to Variational Inequalities Stefan Rosenberger Supervisor: Prof. DI. Dr. techn. Karl Kunisch Institut für Mathematik und wissenschaftliches Rechnen Universität Graz January 26, 2012 Stefan
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 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 informationPOD A-POSTERIORI ERROR BASED INEXACT SQP METHOD FOR BILINEAR ELLIPTIC OPTIMAL CONTROL PROBLEMS
ESAIM: M2AN 46 22) 49 5 DOI:.5/m2an/26 ESAIM: Mathematical Modelling and Numerical Analysis www.esaim-m2an.org POD A-POSTERIORI ERROR BASED INEXACT SQP METHOD FOR BILINEAR ELLIPTIC OPTIMAL CONTROL PROBLEMS
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationHamburger Beiträge zur Angewandten Mathematik
Hamurger Beiträge zur Angewandten Mathematik POD asis updates for nonlinear PDE control Carmen Gräßle, Martin Guisch, Simone Metzdorf, Sarina Rogg, Stefan Volkwein Nr. 16-18 August 16 C. Gräßle, M. Guisch,
More informationRESEARCH ARTICLE. A strategy of finding an initial active set for inequality constrained quadratic programming problems
Optimization Methods and Software Vol. 00, No. 00, July 200, 8 RESEARCH ARTICLE A strategy of finding an initial active set for inequality constrained quadratic programming problems Jungho Lee Computer
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 informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationReceding horizon controller for the baroreceptor loop in a model for the cardiovascular system
SpezialForschungsBereich F 32 Karl Franzens Universität Graz Technische Universität Graz Medizinische Universität Graz Receding horizon controller for the baroreceptor loop in a model for the cardiovascular
More informationProper Orthogonal Decomposition in PDE-Constrained Optimization
Proper Orthogonal Decomposition in PDE-Constrained Optimization K. Kunisch Department of Mathematics and Computational Science University of Graz, Austria jointly with S. Volkwein Dynamic Programming Principle
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 MIXED VARIATIONAL FRAMEWORK FOR THE RADIATIVE TRANSFER EQUATION
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz A MIXED VARIATIONAL FRAMEWORK FOR THE RADIATIVE TRANSFER EQUATION H. Egger M. Schlottbom
More informationA Posteriori Error Estimation for Reduced Order Solutions of Parametrized Parabolic Optimal Control Problems
A Posteriori Error Estimation for Reduced Order Solutions of Parametrized Parabolic Optimal Control Problems Mark Kärcher 1 and Martin A. Grepl Bericht r. 361 April 013 Keywords: Optimal control, reduced
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 informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationOPTIMALITY CONDITIONS AND POD A-POSTERIORI ERROR ESTIMATES FOR A SEMILINEAR PARABOLIC OPTIMAL CONTROL
OPTIMALITY CONDITIONS AND POD A-POSTERIORI ERROR ESTIMATES FOR A SEMILINEAR PARABOLIC OPTIMAL CONTROL O. LASS, S. TRENZ, AND S. VOLKWEIN Abstract. In the present paper the authors consider an optimal control
More informationRobust error estimates for regularization and discretization of bang-bang control problems
Robust error estimates for regularization and discretization of bang-bang control problems Daniel Wachsmuth September 2, 205 Abstract We investigate the simultaneous regularization and discretization of
More informationRiesz bases of Floquet modes in semi-infinite periodic waveguides and implications
Riesz bases of Floquet modes in semi-infinite periodic waveguides and implications Thorsten Hohage joint work with Sofiane Soussi Institut für Numerische und Angewandte Mathematik Georg-August-Universität
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
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 informationMEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS
MEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS KARL KUNISCH, KONSTANTIN PIEPER, AND BORIS VEXLER Abstract. A directional sparsity framework allowing for measure valued controls
More informationPriority Program 1253
Deutsche Forschungsgemeinschaft Priority Program 1253 Optimization with Partial Differential Equations Klaus Deckelnick and Michael Hinze A note on the approximation of elliptic control problems with bang-bang
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 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 informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationSpectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem
Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem L u x f x BC u x g x with the weak problem find u V such that B u,v
More informationNONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian UNDER NONHOMOGENEOUS NEUMANN BOUNDARY CONDITION
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 210, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu NONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian
More informationarxiv: v1 [math.na] 16 Feb 2016
NONLINEAR MODEL ORDER REDUCTION VIA DYNAMIC MODE DECOMPOSITION ALESSANDRO ALLA, J. NATHAN KUTZ arxiv:62.58v [math.na] 6 Feb 26 Abstract. We propose a new technique for obtaining reduced order models for
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationAn Algorithm for Determining the Active Slip Systems in Crystal Plasticity Based on the Principle of Maximum Dissipation
SpezialForschungsBereich F32 Karl Franzens Universität Graz Technische Universität Graz Medizinische Universität Graz An Algorithm for Determining the Active Slip Systems in Crystal Plasticity Based on
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationK. Krumbiegel I. Neitzel A. Rösch
SUFFICIENT OPTIMALITY CONDITIONS FOR THE MOREAU-YOSIDA-TYPE REGULARIZATION CONCEPT APPLIED TO SEMILINEAR ELLIPTIC OPTIMAL CONTROL PROBLEMS WITH POINTWISE STATE CONSTRAINTS K. Krumbiegel I. Neitzel A. Rösch
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
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 informationPerturbation Theory for Self-Adjoint Operators in Krein spaces
Perturbation Theory for Self-Adjoint Operators in Krein spaces Carsten Trunk Institut für Mathematik, Technische Universität Ilmenau, Postfach 10 05 65, 98684 Ilmenau, Germany E-mail: carsten.trunk@tu-ilmenau.de
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 Computational Approach to Study a Logistic Equation
Communications in MathematicalAnalysis Volume 1, Number 2, pp. 75 84, 2006 ISSN 0973-3841 2006 Research India Publications A Computational Approach to Study a Logistic Equation G. A. Afrouzi and S. Khademloo
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 informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationREGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS
REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS fredi tröltzsch 1 Abstract. A class of quadratic optimization problems in Hilbert spaces is considered,
More informationNONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS
Fixed Point Theory, Volume 9, No. 1, 28, 3-16 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS GIOVANNI ANELLO Department of Mathematics University
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationAnalysis in weighted spaces : preliminary version
Analysis in weighted spaces : preliminary version Frank Pacard To cite this version: Frank Pacard. Analysis in weighted spaces : preliminary version. 3rd cycle. Téhéran (Iran, 2006, pp.75.
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 informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationChapitre 1. Cours Reduced Basis Output Bounds Methodes du 1er Semestre
Chapitre 1 Cours Methodes du 1er Semestre 2007-2008 Pr. Christophe Prud'homme christophe.prudhomme@ujf-grenoble.fr Université Joseph Fourier Grenoble 1 1.1 1 2 FEM Approximation Approximation Estimation
More informationA Posteriori Estimates for Cost Functionals of Optimal Control Problems
A Posteriori Estimates for Cost Functionals of Optimal Control Problems Alexandra Gaevskaya, Ronald H.W. Hoppe,2 and Sergey Repin 3 Institute of Mathematics, Universität Augsburg, D-8659 Augsburg, Germany
More informationEmpirical 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 informationProper Orthogonal Decomposition Surrogate Models for Nonlinear Dynamical Systems: Error Estimates and Suboptimal Control
Proper Orthogonal Decomposition Surrogate Models for Nonlinear Dynamical Systems: Error Estimates and Suboptimal Control Michael Hinze 1 and Stefan Volkwein 1 Institut für Numerische Mathematik, TU Dresden,
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationSemismooth Newton Methods for an Optimal Boundary Control Problem of Wave Equations
Semismooth Newton Methods for an Optimal Boundary Control Problem of Wave Equations Axel Kröner 1 Karl Kunisch 2 and Boris Vexler 3 1 Lehrstuhl für Mathematische Optimierung Technische Universität München
More informationNumerical Solution of Heat Equation by Spectral Method
Applied Mathematical Sciences, Vol 8, 2014, no 8, 397-404 HIKARI Ltd, wwwm-hikaricom http://dxdoiorg/1012988/ams201439502 Numerical Solution of Heat Equation by Spectral Method Narayan Thapa Department
More informationKey words. preconditioned conjugate gradient method, saddle point problems, optimal control of PDEs, control and state constraints, multigrid method
PRECONDITIONED CONJUGATE GRADIENT METHOD FOR OPTIMAL CONTROL PROBLEMS WITH CONTROL AND STATE CONSTRAINTS ROLAND HERZOG AND EKKEHARD SACHS Abstract. Optimality systems and their linearizations arising in
More informationSELF-ADJOINTNESS OF SCHRÖDINGER-TYPE OPERATORS WITH SINGULAR POTENTIALS ON MANIFOLDS OF BOUNDED GEOMETRY
Electronic Journal of Differential Equations, Vol. 2003(2003), No.??, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) SELF-ADJOINTNESS
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 155 A posteriori error estimates for stationary slow flows of power-law fluids Michael Bildhauer,
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationCOMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY
COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY BIN HAN AND HUI JI Abstract. In this paper, we provide a family of compactly supported orthonormal complex wavelets with dilation
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 informationGeneralized finite element method using proper orthogonal. decomposition
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2008; 00:1 31 [Version: 2000/01/19 v2.0] Generalized finite element method using proper orthogonal decomposition W.
More informationBALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS
BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007
More informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationLinear Algebra and Dirac Notation, Pt. 2
Linear Algebra and Dirac Notation, Pt. 2 PHYS 500 - Southern Illinois University February 1, 2017 PHYS 500 - Southern Illinois University Linear Algebra and Dirac Notation, Pt. 2 February 1, 2017 1 / 14
More informationA Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems
A Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems D. Meidner and B. Vexler Abstract In this article we discuss a priori error estimates for Galerkin
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 information08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms
(February 24, 2017) 08a. Operators on Hilbert spaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2016-17/08a-ops
More informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More informationA DUALITY ALGORITHM FOR THE OBSTACLE PROBLEM
Ann. Acad. Rom. Sci. Ser. Math. Appl. ISSN 2066-6594 Vol. 5, No. -2 / 203 A DUALITY ALGORITHM FOR THE OBSTACLE PROBLEM Diana Merluşcă Abstract We consider the obstacle problem in Sobolev spaces, of order
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 informationLecture Notes 1: Vector spaces
Optimization-based data analysis Fall 2017 Lecture Notes 1: Vector spaces In this chapter we review certain basic concepts of linear algebra, highlighting their application to signal processing. 1 Vector
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationMartin Gubisch Stefan Volkwein
Universität Konstanz POD a-posteriori error analysis for optimal control problems with mixed control-state constraints Martin Gubisch Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 34, Juni 3 ISSN
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationMultipoint secant and interpolation methods with nonmonotone line search for solving systems of nonlinear equations
Multipoint secant and interpolation methods with nonmonotone line search for solving systems of nonlinear equations Oleg Burdakov a,, Ahmad Kamandi b a Department of Mathematics, Linköping University,
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More information