Efficient three-dimensional variational data assimilation of multi-plane PIV data

Abstract

We perform three-dimensional variational data assimilation (3DVar) using a discrete adjoint approach to optimise the time-averaged momentum equations. The experimental data consist of sparse stereoscopic particle image velocimetry (PIV) measurements collected along $12$ cross-stream planes in the wake of a vehicle-like bluff body at a Reynolds number $Re_L = 5.64 \times 10^5$ based on the streamwise body length. Adjoint localisation is proposed and implemented to reduce the memory footprint of the discrete adjoint method for spatially-varying control variables in 3DVar by confining the control variable space to a user-defined subdomain. Restricting the control variable to $12$ % of the full control space yields a maximum reduction in peak memory of $64$ %, while producing assimilated fields of comparable fidelity with respect to mean velocity and the optimised momentum forcing field. The localised adjoint case improves upon the baseline Spalart--Allmaras turbulence model and recovers the correct asymmetric topology of the complex three-dimensional (3D) recirculation bubble. The assimilated Reynolds shear stress agrees well with the experiment, and the assimilated mean pressure is shown to be physically consistent when correlated with the in-plane vorticity fields. A data efficiency study is also performed, in which the number of planes provided for assimilation is progressively reduced, demonstrating that the data coverage must extend at least to the end of the primary recirculation bubble to adequately constrain the near-wake dynamics. The efficiency that adjoint localisation affords is crucial for assimilating sparse, experimental data for 3D separated flows on fine meshes that can tackle industrial problems of interest.