用神经网络改进数据同化,让冲击波模拟更稳定。
Neural ensemble Kalman filter: Data assimilation for compressible flows with shocks
- 将神经网络嵌入集合卡尔曼滤波,映射流场到参数空间。
- 在激波附近避免了传统方法的虚假振荡和非物理解。
- 适合高精度模拟含冲击波的可压缩流体问题的研究者。
含冲击波的可压缩流体数据同化(DA)极具挑战性,因传统方法常在不确定冲击波附近产生虚假振荡和非物理解。本文聚焦集合卡尔曼滤波(EnKF),发现其性能不佳源于冲击波附近可能出现双峰预报分布,违背了EnKF假设的近高斯分布。为此提出新型神经EnKF:将冲击波流场的预报集合映射至深度神经网络(NN)的参数空间(权重与偏置),并在该空间中进行数据同化。非线性映射能编码流场中的陡峭与平滑特征。通过物理信息迁移学习强制神经网络参数在预报集合中平滑变化,可有效避免虚假振荡与非物理解。在无粘伯格方程、Sod激波管及二维爆炸波等系列数值实验中验证了该方法的有效性。
原文摘要 · Abstract (English)
Data assimilation (DA) for compressible flows with shocks is challenging because many classical DA methods generate spurious oscillations and nonphysical features near uncertain shocks. We focus here on the ensemble Kalman filter (EnKF). We show that the poor performance of the EnKF may be attributed to the bimodal forecast distribution that can arise in the vicinity of an uncertain shock location; this violates the assumptions underpinning the EnKF, which assume a forecast which is close to Gaussian. To address this issue we introduce the new neural EnKF. The basic idea is to systematically embed neural function approximations within ensemble DA by mapping the forecast ensemble of shocked flows to the parameter space (weights and biases) of a deep neural network (NN) and to subsequently perform DA in that space. The nonlinear mapping encodes sharp and smooth flow features in an ensemble of NN parameters. Neural EnKF updates are therefore well-behaved only if the NN parameters vary smoothly within the neural representation of the forecast ensemble. We show that such a smooth variation of network parameters can be enforced via physics-informed transfer learning, and demonstrate that in so-doing the neural EnKF avoids the spurious oscillations and nonphysical features that plague the EnKF. The applicability of the neural EnKF is demonstrated through a series of systematic numerical experiments with the inviscid Burgers' equation, the Sod shock tube, and a two-dimensional blast wave.
Thank you to arXiv for use of its open access interoperability. PaperDance 不是 arXiv 官方产品;中文卡片由大模型生成,请以原文为准。