Introduction Optimization via Adaptive Model Reduction Large-Scale, Constrained Optimization Conclusion Accelerating PDE-Constrained Optimization Problems using Adaptive Reduced-Order Models Matthew J. Zahr Advisor: Charbel Farhat Computational and Mathematical Engineering Stanford University Department of Aerospace and Mechanical Engineering University of Southern California Los Angeles, CA February 26, 2016 Zahr PDE-Constrained Optimization with Adaptive ROMs
• ‒ ‒ ‒ • ‒ • Introduction Optimization via Adaptive Model Reduction Large-Scale, Constrained Optimization Conclusion Multiphysics Optimization Key Player in Next-Gen Problems Current interest in computational physics reaches far beyond analysis of a single configuration of a physical system into design (shape and topology 1 ), control , and uncertainty quantification Micro-Aerial Vehicle EM Launcher Engine System 1 Emergence of additive manufacturing technologies has made topology optimization increasingly relevant, particularly in DOE. Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Optimization via Adaptive Model Reduction Large-Scale, Constrained Optimization Conclusion Topology Optimization and Additive Manufacturing 2 Emergence of AM has made TO an increasingly relevant topic AM+TO lead to highly e ffi cient designs that could not be realized previously Challenges: smooth topologies require very fine meshes and modeling of complex manufacturing process 2 MIT Technology Review , Top 10 Technological Breakthrough 2013 Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion PDE-Constrained Optimization I Goal: Rapidly solve PDE-constrained optimization problem of the form minimize J ( u , µ ) u ∈ R n u , µ ∈ R n µ subject to r ( u , µ ) = 0 where r : R n u ⇥ R n µ ! R n u is the discretized partial di ff erential equation J : R n u ⇥ R n µ ! R is the objective function u 2 R n u is the PDE state vector µ 2 R n µ is the vector of parameters red indicates a large-scale quantity, O ( mesh ) Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Nested Approach to PDE-Constrained Optimization Virtually all expense emanates from primal/dual PDE solvers Optimizer Primal PDE Dual PDE Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Nested Approach to PDE-Constrained Optimization Virtually all expense emanates from primal/dual PDE solvers Optimizer µ Primal PDE Dual PDE Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Nested Approach to PDE-Constrained Optimization Virtually all expense emanates from primal/dual PDE solvers Optimizer J ( u , µ ) Primal PDE Dual PDE Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Nested Approach to PDE-Constrained Optimization Virtually all expense emanates from primal/dual PDE solvers Optimizer µ J ( u , µ ) u Primal PDE Dual PDE Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Nested Approach to PDE-Constrained Optimization Virtually all expense emanates from primal/dual PDE solvers Optimizer d J d µ ( u , µ ) J ( u , µ ) Primal PDE Dual PDE Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Projection-Based Model Reduction to Reduce PDE Size Model Order Reduction (MOR) assumption: state vector lies in low-dimensional subspace ∂ u ∂ u r u ⇡ Φ u u r ∂ µ ⇡ Φ u ∂ µ where 2 R n u ⇥ k u is the reduced basis ⇥ φ 1 φ k u ⇤ Φ u = · · · u u u r 2 R k u are the reduced coordinates of u n u � k u Substitute assumption into High-Dimensional Model (HDM), r ( u , µ ) = 0, and project onto test subspace Ψ u 2 R n u × k u T r ( Φ u u r , µ ) = 0 Ψ u Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Connection to Finite Element Method: Hierarchical Subspaces S S - infinite-dimensional trial space Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Connection to Finite Element Method: Hierarchical Subspaces S h S S - infinite-dimensional trial space S h - (large) finite-dimensional trial space Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Connection to Finite Element Method: Hierarchical Subspaces S k h S h S S - infinite-dimensional trial space S h - (large) finite-dimensional trial space S k h - (small) finite-dimensional trial space S k h ⇢ S h ⇢ S Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Few Global, Data-Driven Basis Functions v. Many Local Ones Instead of using traditional local shape functions (e.g., FEM), use global shape functions Instead of a-priori, analytical shape functions, leverage data-rich computing environment by using data-driven modes Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Definition of Φ u : Data-Driven Reduction State-Sensitivity Proper Orthogonal Decomposition (POD) Collect state and sensitivity snapshots by sampling HDM ⇥ u ( µ 1 ) u ( µ n ) ⇤ X = u ( µ 2 ) · · · h i ∂ u ∂ u ∂ u Y = ∂ µ ( µ 1 ) ∂ µ ( µ 2 ) · · · ∂ µ ( µ n ) Use Proper Orthogonal Decomposition to generate reduced basis for each individually Φ X = POD( X ) Φ Y = POD( Y ) Concatenate and orthogonalize to get reduced-order basis ⇣h i⌘ u ( µ ∗ ) ∂ u ∂ µ ( µ ∗ ) Φ u = QR Φ X Φ Y Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion Definition of Ψ u : Minimum-Residual ROM Least-Squares Petrov-Galerkin (LSPG) 3 projection Ψ u = ∂ r ∂ u Φ u Minimum-Residual Property A ROM possesses the minimum-residual property if Ψ u r ( Φ u u r , µ ) = 0 is equivalent to the optimality condition of ( Θ � 0) minimize || r ( Φ u u r , µ ) || Θ u r ∈ R k u Implications Recover exact solution when basis not truncated (consistent 3 ) Monotonic improvement of solution as basis size increases Ensures sensitivity information in Φ u cannot degrade state approximation 4 LSPG possesses minimum-residual property 3 [Bui-Thanh et al., 2008] 4 [Fahl, 2001] Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion O ffl ine-Online Approach to Optimization . . . Schematic µ -space Breakdown of Computational E ff ort Zahr PDE-Constrained Optimization with Adaptive ROMs
Introduction Model Order Reduction Optimization via Adaptive Model Reduction Non-Quadratic Trust-Region Solver Large-Scale, Constrained Optimization Shape Optimization: Airfoil Design Conclusion O ffl ine-Online Approach to Optimization HDM HDM . . . HDM HDM Schematic µ -space HDM HDM · · · HDM HDM Breakdown of Computational E ff ort Zahr PDE-Constrained Optimization with Adaptive ROMs
Recommend
More recommend
