Yuhui Yin , Chenhui Kou , Shengkun Jia , Lu Lu , Xigang Yuan , Yiqing Luo
{"title":"PCDMD: Physics-constrained dynamic mode decomposition for accurate and robust forecasting of dynamical systems with imperfect data and physics","authors":"Yuhui Yin , Chenhui Kou , Shengkun Jia , Lu Lu , Xigang Yuan , Yiqing Luo","doi":"10.1016/j.cpc.2024.109303","DOIUrl":null,"url":null,"abstract":"<div><p>The dynamic mode decomposition (DMD) method has attracted widespread attention as a representative modal-decomposition method and can build a predictive model. However, DMD may give predicted results that deviate from physical reality in some scenarios, such as dealing with translation problems or noisy data. Here, we propose a physics-constrained DMD (PCDMD) method to address this issue. The proposed PCDMD method first employs a data-driven model using DMD, then calculates the residual of the physical equations, and finally corrects the predicted results using Kalman filter and gain coefficients. In this way, the PCDMD method can integrate the physics-informed equations with the data-driven model generated by DMD. Numerical experiments are conducted using PCDMD, including the Allen–Cahn, advection-diffusion, Burgers' equations and lid-driven cavity flow. The results demonstrate that the proposed PCDMD method can reduce the reconstruction and prediction errors by 1%-10% by incorporating physical constraints. Regarding noisy datasets and imperfect physical constraints, PCDMD can still ensure that the predicted results satisfy the physical constraints, thereby reducing errors.</p></div><div><h3>Program summary</h3><p><em>Program title:</em> PCDMD</p><p><em>Dataset link:</em> <span>https://github.com/YinYuhuiTJU/PCDMD</span><svg><path></path></svg></p><p><em>Licensing provisions:</em> MIT</p><p><em>Programming language:</em> Python</p></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"304 ","pages":"Article 109303"},"PeriodicalIF":7.2000,"publicationDate":"2024-07-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer Physics Communications","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0010465524002261","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
The dynamic mode decomposition (DMD) method has attracted widespread attention as a representative modal-decomposition method and can build a predictive model. However, DMD may give predicted results that deviate from physical reality in some scenarios, such as dealing with translation problems or noisy data. Here, we propose a physics-constrained DMD (PCDMD) method to address this issue. The proposed PCDMD method first employs a data-driven model using DMD, then calculates the residual of the physical equations, and finally corrects the predicted results using Kalman filter and gain coefficients. In this way, the PCDMD method can integrate the physics-informed equations with the data-driven model generated by DMD. Numerical experiments are conducted using PCDMD, including the Allen–Cahn, advection-diffusion, Burgers' equations and lid-driven cavity flow. The results demonstrate that the proposed PCDMD method can reduce the reconstruction and prediction errors by 1%-10% by incorporating physical constraints. Regarding noisy datasets and imperfect physical constraints, PCDMD can still ensure that the predicted results satisfy the physical constraints, thereby reducing errors.
期刊介绍:
The focus of CPC is on contemporary computational methods and techniques and their implementation, the effectiveness of which will normally be evidenced by the author(s) within the context of a substantive problem in physics. Within this setting CPC publishes two types of paper.
Computer Programs in Physics (CPiP)
These papers describe significant computer programs to be archived in the CPC Program Library which is held in the Mendeley Data repository. The submitted software must be covered by an approved open source licence. Papers and associated computer programs that address a problem of contemporary interest in physics that cannot be solved by current software are particularly encouraged.
Computational Physics Papers (CP)
These are research papers in, but are not limited to, the following themes across computational physics and related disciplines.
mathematical and numerical methods and algorithms;
computational models including those associated with the design, control and analysis of experiments; and
algebraic computation.
Each will normally include software implementation and performance details. The software implementation should, ideally, be available via GitHub, Zenodo or an institutional repository.In addition, research papers on the impact of advanced computer architecture and special purpose computers on computing in the physical sciences and software topics related to, and of importance in, the physical sciences may be considered.