This paper suggests reduced-order modeling using the Galerkin proper orthogonal decomposition (POD) to find approximate solutions for parabolic partial differential equations. We first transform a parabolic partial differential equation to the frequency-dependent elliptic equations using the Fourier integral transform in time. Such a frequency-domain method enables efficiently implementing a parallel computation to approximate the solutions because the frequency-variable elliptic equations have independent frequencies. Then, we introduce reduced-order modeling to determine approximate solutions of the frequency-variable elliptic equations quickly. A set of snapshots consists of the finite element solutions of the frequency-variable elliptic equations with some selected frequencies. The solutions are approximated using the general basis of the high-dimensional finite element space in a Hilbert space. reduced-order modeling employs the Galerkin POD for the snapshot subspace spanned by a set of snapshots. An orthonormal basis for the snapshot space can be easily computed using the spectral decomposition of the correlation matrix of the snapshots. Additionally, using an appropriate low-order basis of the snapshot space allows approximating the solutions of the frequency-variable elliptic equations quickly, where the approximate solutions are used for the inverse Fourier transforms to determine the approximated solutions in the time variable. Several numerical tests based on the finite element method are presented to asses the efficient performances of the suggested approaches.
Citation: Jeong-Kweon Seo, Byeong-Chun Shin. Reduced-order modeling using the frequency-domain method for parabolic partial differential equations[J]. AIMS Mathematics, 2023, 8(7): 15255-15268. doi: 10.3934/math.2023779
This paper suggests reduced-order modeling using the Galerkin proper orthogonal decomposition (POD) to find approximate solutions for parabolic partial differential equations. We first transform a parabolic partial differential equation to the frequency-dependent elliptic equations using the Fourier integral transform in time. Such a frequency-domain method enables efficiently implementing a parallel computation to approximate the solutions because the frequency-variable elliptic equations have independent frequencies. Then, we introduce reduced-order modeling to determine approximate solutions of the frequency-variable elliptic equations quickly. A set of snapshots consists of the finite element solutions of the frequency-variable elliptic equations with some selected frequencies. The solutions are approximated using the general basis of the high-dimensional finite element space in a Hilbert space. reduced-order modeling employs the Galerkin POD for the snapshot subspace spanned by a set of snapshots. An orthonormal basis for the snapshot space can be easily computed using the spectral decomposition of the correlation matrix of the snapshots. Additionally, using an appropriate low-order basis of the snapshot space allows approximating the solutions of the frequency-variable elliptic equations quickly, where the approximate solutions are used for the inverse Fourier transforms to determine the approximated solutions in the time variable. Several numerical tests based on the finite element method are presented to asses the efficient performances of the suggested approaches.
[1] | J. A. Atwell, B. B. King, Reduced-order controllers for spatially distributed systems via proper orthogonal decomposition, SIAM J. Sci. Comput., 26 (2004), 128–151. https://doi.org/10.1137/S1064827599360091 doi: 10.1137/S1064827599360091 |
[2] | J. Burkardt, M. Gunzburger, H. C. Lee, Centroidal Voronoi tessellation-based reduced-order modeling of complex systems, SIAM J. Sci. Comput., 28 (2006), 459–484. https://doi.org/10.1137/5106482750342221x doi: 10.1137/5106482750342221x |
[3] | J. Burkardt, M. Gunzburger, H. C. Lee, POD and CVT-based reduced-order modeling of Navier-Stokes flows, Comput. Meth. Appl. Mech. Eng., 196 (2006), 337–355. https://doi.org/10.1016/j.cma.2006.04.004 doi: 10.1016/j.cma.2006.04.004 |
[4] | J. Douglas, Jr., J. E. Santos, D. Sheen, Approximation of scalar waves in the space-frequency domain, Math. Models Meth. Appl. Sci., 4 (1994), 509–531. https://doi.org/10.1142/S0218202594000297 doi: 10.1142/S0218202594000297 |
[5] | J. Douglas, Jr., J. E. Santos, D. Sheen, L. S. Bennethum, Frequency domain treatment of one-dimensional scalar waves, Math. Model Meth. Appl. Sci., 3 (1993), 171–194. https://doi.org/10.1142/S0218202593000102 doi: 10.1142/S0218202593000102 |
[6] | X. Feng, D. Sheen, An elliptic regularity estimate for a problem arising from the frequency domain treatment of waves, Trans. Ams. Math. Soc., 346 (1994), 475–487. |
[7] | K. Kunisch, S. Volkwein, Control of the Burgers equation by a reduced-order approach using proper orthogonal decomposition, J. Optim. Theory Appl., 102 (1999), 345–371. https://doi.org/10.1023/A:1021732508059 doi: 10.1023/A:1021732508059 |
[8] | K. Kunisch, S. Volkwein, Galerkin proper orthogonal decomposition methods for parabolic equations, Numer. Math., 90 (2001), 117–148. https://doi.org/10.1007/s002110100282 doi: 10.1007/s002110100282 |
[9] | K. Kunisch, S. Volkwein, Galerkin proper orthogonal decomposition methods for a general equation in fluid dynamics, SIAM J. Numer. Anal., 40 (2002), 492–515. https://doi.org/10.1137/S0036142900382612 doi: 10.1137/S0036142900382612 |
[10] | C. O. Lee, J. Lee, D. Sheen, Y. Yeom, A frequency-domain parallel method for the numerical approximation of parabolic problems, Comput. Meth. Appl. Mech. Eng., 169 (1999), 19–29. https://doi.org/10.1016/S0045-7825(98)00168-6 doi: 10.1016/S0045-7825(98)00168-6 |
[11] | C. O. Lee, J. Lee, D. Sheen, Frequency domain formulation of linearized Navier-Stokes equations, Comput. Meth. Appl. Mech. Eng., 187 (2000), 351–362. https://doi.org/10.1016/S0045-7825(99)00132-2 doi: 10.1016/S0045-7825(99)00132-2 |
[12] | J. Lee, D. Sheen, An accurate numerical inversion of Laplace transforms based on the location of their poles, Comput. Math. Appl., 48 (2004), 1415-1423. https://doi.org/10.1016/j.camwa.2004.08.003 doi: 10.1016/j.camwa.2004.08.003 |
[13] | J. Lee, D. Sheen, A parallel method for backward parabolic problems based on the Laplace transformation, SIAM J. Numer. Anal., 44 (2006), 1466–1486. https://doi.org/10.1137/050624649 doi: 10.1137/050624649 |
[14] | H. C. Lee, S. W. Lee, G. R. Piao, Reduced-order modeling of Burgers equations based on centroidal Voronoi tessellation, Int. J. Numer. Anal. Model., 4 (2007), 559–583. |
[15] | R. Li, Q. Wu, S. Zhu, Proper orthogonal decomposition with SUPG-stabilized isogeometric analysis for reduced order modeling of unsteady convection-dominated convection-diffusion-reaction problems, J. Comput. Phys., 387 (2019), 280–302. https://doi.org/10.1016/j.jcp.2019.02.051 doi: 10.1016/j.jcp.2019.02.051 |
[16] | Y. C. Liang, H. P. Lee, S. P. Lim, W. Z. Lin, K. H. Lee, C. G. Wu, Proper orthogonal decomposition and its application–part 1: theory, J. Sound Vib., 252 (2002), 527–544. https://doi.org/10.1006/jsvi.2001.4041 doi: 10.1006/jsvi.2001.4041 |
[17] | J. K. Seo, B. C. Shin, Numerical solutions of Burgers equation by reduced-order modeling based on pseudo-spectral collocation method, J. Korean Soc. Ind. Appl. Math., 19 (2015), 123–135. https://doi.org/10.12941/jksiam.2015.19.123 doi: 10.12941/jksiam.2015.19.123 |
[18] | D. Sheen, I. H. Sloan, V. Thomée, A parallel method for time-discretization of parabolic problems based on contour integral representation and quadrature, Math. Comp., 69 (2000), 177-195. |
[19] | D. Sheen, I. H. Sloan, V. Thomée, A parallel method for time-discretization of parabolic equations based on Laplace transformation and quadrature, IMA J. Numer. Anal., 23 (2003), 269–299. https://doi.org/10.1093/imanum/23.2.269 doi: 10.1093/imanum/23.2.269 |
[20] | S. Zhu, L. Dede, A. Quarteroni, Isogeometric analysis and proper orthogonal decomposition for parabolic problems, Numer. Math., 135 (2017), 333–370. https://doi.org/10.1007/s00211-016-0802-5 doi: 10.1007/s00211-016-0802-5 |