{"title":"A time-relaxation reduced order model for the turbulent channel flow","authors":"Ping-Hsuan Tsai , Paul Fischer , Traian Iliescu","doi":"10.1016/j.jcp.2024.113563","DOIUrl":null,"url":null,"abstract":"<div><div>Regularized reduced order models (Reg-ROMs) are stabilization strategies that leverage spatial filtering to alleviate the spurious numerical oscillations generally displayed by the classical Galerkin ROM (G-ROM) in under-resolved numerical simulations of turbulent flows. In this paper, we propose a new Reg-ROM, the time-relaxation ROM (TR-ROM), which filters the marginally resolved scales. We compare the new TR-ROM with the two other Reg-ROMs in current use, i.e., the Leray ROM (L-ROM) and the evolve-filter-relax ROM (EFR-ROM) and one eddy viscosity model, the mixing-length model, in the numerical simulation of the turbulent channel flow at <span><math><mi>R</mi><msub><mrow><mi>e</mi></mrow><mrow><mi>τ</mi></mrow></msub><mo>=</mo><mn>180</mn></math></span> and <span><math><mi>R</mi><msub><mrow><mi>e</mi></mrow><mrow><mi>τ</mi></mrow></msub><mo>=</mo><mn>395</mn></math></span> in both the reproduction and the predictive regimes. For each Reg-ROM, we investigate two different filters: (i) the differential filter (DF), and (ii) the higher-order algebraic filter (HOAF). In our numerical investigation, we monitor the Reg-ROM performance with respect to the ROM dimension, <em>N</em>, and the filter order. We also perform sensitivity studies of the three Reg-ROMs with respect to the time interval, relaxation parameter, and filter radius. The numerical results yield the following conclusions: (i) All three Reg-ROMs are significantly more accurate than the G-ROM. (ii) All three Reg-ROMs are more accurate than the ROM projection in terms of Reynolds stresses. (iii) With the optimal parameter values, the new TR-ROM yields more accurate results than the L-ROM and the EFR-ROM in all tests. (iv) The new TR-ROM is more accurate than the mixing-length ROM. (v) For most <em>N</em> values, DF yields the most accurate results for all three Reg-ROMs. (vi) The optimal parameters trained in the reproduction regime are also optimal for the predictive regime for most <em>N</em> values, demonstrating the Reg-ROM predictive capabilities. (vii) All three Reg-ROMs are sensitive to the filter order and the filter radius, and the EFR-ROM and the TR-ROM are sensitive to the relaxation parameter. (viii) The optimal range for the filter radius and the effect of relaxation parameter are similar for the two <span><math><mi>R</mi><msub><mrow><mi>e</mi></mrow><mrow><mi>τ</mi></mrow></msub></math></span> values.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"521 ","pages":"Article 113563"},"PeriodicalIF":3.8000,"publicationDate":"2024-11-06","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/S0021999124008118","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
Regularized reduced order models (Reg-ROMs) are stabilization strategies that leverage spatial filtering to alleviate the spurious numerical oscillations generally displayed by the classical Galerkin ROM (G-ROM) in under-resolved numerical simulations of turbulent flows. In this paper, we propose a new Reg-ROM, the time-relaxation ROM (TR-ROM), which filters the marginally resolved scales. We compare the new TR-ROM with the two other Reg-ROMs in current use, i.e., the Leray ROM (L-ROM) and the evolve-filter-relax ROM (EFR-ROM) and one eddy viscosity model, the mixing-length model, in the numerical simulation of the turbulent channel flow at and in both the reproduction and the predictive regimes. For each Reg-ROM, we investigate two different filters: (i) the differential filter (DF), and (ii) the higher-order algebraic filter (HOAF). In our numerical investigation, we monitor the Reg-ROM performance with respect to the ROM dimension, N, and the filter order. We also perform sensitivity studies of the three Reg-ROMs with respect to the time interval, relaxation parameter, and filter radius. The numerical results yield the following conclusions: (i) All three Reg-ROMs are significantly more accurate than the G-ROM. (ii) All three Reg-ROMs are more accurate than the ROM projection in terms of Reynolds stresses. (iii) With the optimal parameter values, the new TR-ROM yields more accurate results than the L-ROM and the EFR-ROM in all tests. (iv) The new TR-ROM is more accurate than the mixing-length ROM. (v) For most N values, DF yields the most accurate results for all three Reg-ROMs. (vi) The optimal parameters trained in the reproduction regime are also optimal for the predictive regime for most N values, demonstrating the Reg-ROM predictive capabilities. (vii) All three Reg-ROMs are sensitive to the filter order and the filter radius, and the EFR-ROM and the TR-ROM are sensitive to the relaxation parameter. (viii) The optimal range for the filter radius and the effect of relaxation parameter are similar for the two values.
期刊介绍:
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.