Yifei Guan, Pedram Hassanzadeh, Tapio Schneider, Oliver Dunbar, Daniel Zhengyu Huang, Jinlong Wu, Ignacio Lopez-Gomez
{"title":"Online learning of eddy-viscosity and backscattering closures for geophysical turbulence using ensemble Kalman inversion","authors":"Yifei Guan, Pedram Hassanzadeh, Tapio Schneider, Oliver Dunbar, Daniel Zhengyu Huang, Jinlong Wu, Ignacio Lopez-Gomez","doi":"arxiv-2409.04985","DOIUrl":null,"url":null,"abstract":"Different approaches to using data-driven methods for subgrid-scale closure\nmodeling have emerged recently. Most of these approaches are data-hungry, and\nlack interpretability and out-of-distribution generalizability. Here, we use\n{online} learning to address parametric uncertainty of well-known physics-based\nlarge-eddy simulation (LES) closures: the Smagorinsky (Smag) and Leith\neddy-viscosity models (1 free parameter) and the Jansen-Held (JH)\nbackscattering model (2 free parameters). For 8 cases of 2D geophysical\nturbulence, optimal parameters are estimated, using ensemble Kalman inversion\n(EKI), such that for each case, the LES' energy spectrum matches that of direct\nnumerical simulation (DNS). Only a small training dataset is needed to\ncalculate the DNS spectra (i.e., the approach is {data-efficient}). We find the\noptimized parameter(s) of each closure to be constant across broad flow regimes\nthat differ in dominant length scales, eddy/jet structures, and dynamics,\nsuggesting that these closures are {generalizable}. In a-posteriori tests based\non the enstrophy spectra and probability density functions (PDFs) of vorticity,\nLES with optimized closures outperform the baselines (LES with standard Smag,\ndynamic Smag or Leith), particularly at the tails of the PDFs (extreme events).\nIn a-priori tests, the optimized JH significantly outperforms the baselines and\noptimized Smag and Leith in terms of interscale enstrophy and energy transfers\n(still, optimized Smag noticeably outperforms standard Smag). The results show\nthe promise of combining advances in physics-based modeling (e.g., JH) and\ndata-driven modeling (e.g., {online} learning with EKI) to develop\ndata-efficient frameworks for accurate, interpretable, and generalizable\nclosures.","PeriodicalId":501270,"journal":{"name":"arXiv - PHYS - Geophysics","volume":"2 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-09-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - PHYS - Geophysics","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.04985","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
Different approaches to using data-driven methods for subgrid-scale closure
modeling have emerged recently. Most of these approaches are data-hungry, and
lack interpretability and out-of-distribution generalizability. Here, we use
{online} learning to address parametric uncertainty of well-known physics-based
large-eddy simulation (LES) closures: the Smagorinsky (Smag) and Leith
eddy-viscosity models (1 free parameter) and the Jansen-Held (JH)
backscattering model (2 free parameters). For 8 cases of 2D geophysical
turbulence, optimal parameters are estimated, using ensemble Kalman inversion
(EKI), such that for each case, the LES' energy spectrum matches that of direct
numerical simulation (DNS). Only a small training dataset is needed to
calculate the DNS spectra (i.e., the approach is {data-efficient}). We find the
optimized parameter(s) of each closure to be constant across broad flow regimes
that differ in dominant length scales, eddy/jet structures, and dynamics,
suggesting that these closures are {generalizable}. In a-posteriori tests based
on the enstrophy spectra and probability density functions (PDFs) of vorticity,
LES with optimized closures outperform the baselines (LES with standard Smag,
dynamic Smag or Leith), particularly at the tails of the PDFs (extreme events).
In a-priori tests, the optimized JH significantly outperforms the baselines and
optimized Smag and Leith in terms of interscale enstrophy and energy transfers
(still, optimized Smag noticeably outperforms standard Smag). The results show
the promise of combining advances in physics-based modeling (e.g., JH) and
data-driven modeling (e.g., {online} learning with EKI) to develop
data-efficient frameworks for accurate, interpretable, and generalizable
closures.