Ivan Zanardi, Alberto Padovan, Daniel J. Bodony, Marco Panesi
{"title":"热化学非平衡气体混合物的Petrov-Galerkin模型还原","authors":"Ivan Zanardi, Alberto Padovan, Daniel J. Bodony, Marco Panesi","doi":"10.1016/j.jcp.2025.113999","DOIUrl":null,"url":null,"abstract":"<div><div>State-specific thermochemical collisional models are crucial to accurately describe the physics of systems involving nonequilibrium plasmas, but they are also computationally expensive and impractical for large-scale, multi-dimensional simulations. Historically, computational cost has been mitigated by using empirical and physics-based arguments to reduce the complexity of the governing equations. However, the resulting models are often inaccurate and they fail to capture the important features of the original physics. Additionally, the construction of these models is often impractical, as it requires extensive user supervision and time-consuming parameter tuning. In this paper, we address these issues through an easily implementable and computationally efficient model reduction pipeline based on the Petrov-Galerkin projection of the nonlinear kinetic equations onto a low-dimensional subspace. Our approach is justified by the observation that kinetic systems in thermal nonequilibrium tend to exhibit low-rank dynamics that rapidly drive the state towards a low-dimensional subspace that can be exploited for reduced-order modeling. Furthermore, despite the nonlinear nature of the governing equations, we observe that the dynamics of these systems evolve on subspaces that can be accurately identified using the linearized equations about thermochemical equilibrium steady states, and we shall see that this allows us to significantly reduce the cost associated with the construction of the model. The approach is demonstrated on two distinct thermochemical systems: a rovibrational collisional model for the O<sub>2</sub>-O system, and a vibrational collisional model for the combined O<sub>2</sub>-O and O<sub>2</sub>-O<sub>2</sub> systems. Our method achieves high accuracy, with relative errors of less than 1% for macroscopic quantities (i.e., moments) and 10% for microscopic quantities (i.e., energy levels population), while also delivering excellent compression rates and speedups, outperforming existing state-of-the-art techniques.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"533 ","pages":"Article 113999"},"PeriodicalIF":3.8000,"publicationDate":"2025-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Petrov-Galerkin model reduction for thermochemical nonequilibrium gas mixtures\",\"authors\":\"Ivan Zanardi, Alberto Padovan, Daniel J. Bodony, Marco Panesi\",\"doi\":\"10.1016/j.jcp.2025.113999\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><div>State-specific thermochemical collisional models are crucial to accurately describe the physics of systems involving nonequilibrium plasmas, but they are also computationally expensive and impractical for large-scale, multi-dimensional simulations. Historically, computational cost has been mitigated by using empirical and physics-based arguments to reduce the complexity of the governing equations. However, the resulting models are often inaccurate and they fail to capture the important features of the original physics. Additionally, the construction of these models is often impractical, as it requires extensive user supervision and time-consuming parameter tuning. In this paper, we address these issues through an easily implementable and computationally efficient model reduction pipeline based on the Petrov-Galerkin projection of the nonlinear kinetic equations onto a low-dimensional subspace. Our approach is justified by the observation that kinetic systems in thermal nonequilibrium tend to exhibit low-rank dynamics that rapidly drive the state towards a low-dimensional subspace that can be exploited for reduced-order modeling. Furthermore, despite the nonlinear nature of the governing equations, we observe that the dynamics of these systems evolve on subspaces that can be accurately identified using the linearized equations about thermochemical equilibrium steady states, and we shall see that this allows us to significantly reduce the cost associated with the construction of the model. The approach is demonstrated on two distinct thermochemical systems: a rovibrational collisional model for the O<sub>2</sub>-O system, and a vibrational collisional model for the combined O<sub>2</sub>-O and O<sub>2</sub>-O<sub>2</sub> systems. Our method achieves high accuracy, with relative errors of less than 1% for macroscopic quantities (i.e., moments) and 10% for microscopic quantities (i.e., energy levels population), while also delivering excellent compression rates and speedups, outperforming existing state-of-the-art techniques.</div></div>\",\"PeriodicalId\":352,\"journal\":{\"name\":\"Journal of Computational Physics\",\"volume\":\"533 \",\"pages\":\"Article 113999\"},\"PeriodicalIF\":3.8000,\"publicationDate\":\"2025-04-10\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Journal of Computational Physics\",\"FirstCategoryId\":\"101\",\"ListUrlMain\":\"https://www.sciencedirect.com/science/article/pii/S0021999125002827\",\"RegionNum\":2,\"RegionCategory\":\"物理与天体物理\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q2\",\"JCRName\":\"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0021999125002827","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
Petrov-Galerkin model reduction for thermochemical nonequilibrium gas mixtures
State-specific thermochemical collisional models are crucial to accurately describe the physics of systems involving nonequilibrium plasmas, but they are also computationally expensive and impractical for large-scale, multi-dimensional simulations. Historically, computational cost has been mitigated by using empirical and physics-based arguments to reduce the complexity of the governing equations. However, the resulting models are often inaccurate and they fail to capture the important features of the original physics. Additionally, the construction of these models is often impractical, as it requires extensive user supervision and time-consuming parameter tuning. In this paper, we address these issues through an easily implementable and computationally efficient model reduction pipeline based on the Petrov-Galerkin projection of the nonlinear kinetic equations onto a low-dimensional subspace. Our approach is justified by the observation that kinetic systems in thermal nonequilibrium tend to exhibit low-rank dynamics that rapidly drive the state towards a low-dimensional subspace that can be exploited for reduced-order modeling. Furthermore, despite the nonlinear nature of the governing equations, we observe that the dynamics of these systems evolve on subspaces that can be accurately identified using the linearized equations about thermochemical equilibrium steady states, and we shall see that this allows us to significantly reduce the cost associated with the construction of the model. The approach is demonstrated on two distinct thermochemical systems: a rovibrational collisional model for the O2-O system, and a vibrational collisional model for the combined O2-O and O2-O2 systems. Our method achieves high accuracy, with relative errors of less than 1% for macroscopic quantities (i.e., moments) and 10% for microscopic quantities (i.e., energy levels population), while also delivering excellent compression rates and speedups, outperforming existing state-of-the-art techniques.
期刊介绍:
Journal of Computational Physics thoroughly treats the computational aspects of physical problems, presenting techniques for the numerical solution of mathematical equations arising in all areas of physics. The journal seeks to emphasize methods that cross disciplinary boundaries.
The Journal of Computational Physics also publishes short notes of 4 pages or less (including figures, tables, and references but excluding title pages). Letters to the Editor commenting on articles already published in this Journal will also be considered. Neither notes nor letters should have an abstract.