Spectral adjoint-based assimilation of sparse data in unsteady simulations of turbulent flows

TL;DR

This paper introduces a Fourier-based data assimilation framework for URANS simulations of turbulent flows, improving accuracy and efficiency by incorporating a divergence-free forcing term and using adjoint methods in Fourier space.

Contribution

The novel approach combines Fourier transform techniques with adjoint-based data assimilation in URANS, enabling efficient correction of turbulence model errors without specialized solvers.

Findings

Accurately reconstructs mean flow around a cylinder at Re=3900.

Improves vortex shedding frequency prediction with zeroth mode data.

Captures low-frequency flow dynamics using first mode data.

Abstract

The URANS equations provide a computationally efficient tool to simulate unsteady turbulent flows for a wide range of applications. To account for the errors introduced by the turbulence closure model, recent works have adopted data assimilation (DA) to enhance their predictive capabilities. Recognizing the challenges posed by the computational cost of 4DVar DA for unsteady flows, we propose a 3DVar DA framework that incorporates a time-discrete Fourier transform of the URANS equations, facilitating the use of the stationary discrete adjoint method in Fourier space. Central to our methodology is the introduction of a corrective, divergence-free, and unsteady forcing term, derived from a Fourier series expansion, into the URANS equations. This term aims at mitigating discrepancies in the modeled divergence of Reynolds stresses, allowing for the tuning of stationary parameters across different Fourier modes. Our implementation is built upon an extended version of the coupled URANS solver in OpenFOAM, enhanced to compute adjoint variables and gradients. This design choice ensures straightforward applicability to various flow setups and solvers, eliminating the need for specialized harmonic solvers. A gradient-based optimizer is employed to minimize discrepancies between simulated results and sparse velocity reference data. The effectiveness of our approach is demonstrated through its application to flow around a two-dimensional circular cylinder at a Reynolds number of 3900. Our results highlight the method's ability to reconstruct mean flow accurately and improve the vortex shedding frequency through the assimilation of zeroth mode data. Additionally, the assimilation of first mode data further enhances the simulation's capability to capture low-frequency dynamics of the flow, and finally, it runs efficiently by leveraging a coarse mesh.