<p>In this paper, we propose the <i>Frequency Reduced-Basis</i> method: a reduced-basis method to solve time-dependent partial differential equations based on the Laplace transform. Unlike traditional approaches, we begin by applying the Laplace transform to the evolution problem, yielding a time-independent boundary value problem that depends on the complex Laplace parameter. First, in an offline stage, we appropriately sample the Laplace parameter and solve the collection of problems using the finite element method. Next, we apply a Proper Orthogonal Decomposition (POD) to this collection of solutions in order to obtain a reduced basis that is of dimension much smaller than that of the original solution space. This reduced basis, in turn, is then used to solve the evolution problem using any suitable time-stepping method. A key insight to justify the formulation of the method relies on Hardy spaces of analytic functions. By applying the Paley-Wiener theorem we can then define an isometry between the solution of the time-dependent problem and its Laplace transform. As a consequence of this result, one may conclude that computing a POD with samples taken in the Laplace domain produces an exponentially accurate reduced basis for the time-dependent problem. Numerical experiments characterising the performance of the method, in terms of accuracy and speed-up, are included for a variety of relevant time-evolution equations.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

The Frequency Reduced-Basis method: reduced order models for time-dependent problems using the Laplace transform

  • Ricardo Reyes

摘要

In this paper, we propose the Frequency Reduced-Basis method: a reduced-basis method to solve time-dependent partial differential equations based on the Laplace transform. Unlike traditional approaches, we begin by applying the Laplace transform to the evolution problem, yielding a time-independent boundary value problem that depends on the complex Laplace parameter. First, in an offline stage, we appropriately sample the Laplace parameter and solve the collection of problems using the finite element method. Next, we apply a Proper Orthogonal Decomposition (POD) to this collection of solutions in order to obtain a reduced basis that is of dimension much smaller than that of the original solution space. This reduced basis, in turn, is then used to solve the evolution problem using any suitable time-stepping method. A key insight to justify the formulation of the method relies on Hardy spaces of analytic functions. By applying the Paley-Wiener theorem we can then define an isometry between the solution of the time-dependent problem and its Laplace transform. As a consequence of this result, one may conclude that computing a POD with samples taken in the Laplace domain produces an exponentially accurate reduced basis for the time-dependent problem. Numerical experiments characterising the performance of the method, in terms of accuracy and speed-up, are included for a variety of relevant time-evolution equations.