Abstract:Reynolds-averaged Navier--Stokes (RANS) closure must be sensitive to the flow physics, including nonlocality and anisotropy of the effective eddy viscosity. Recent approaches used forced direct numerical simulations to probe these effects, including the macroscopic forcing method (MFM) of Mani and Park ($\textit{Phys. Rev. Fluids}$ $\textbf{6}$, 054607 (2021)) and the Green's function approach of Hamba ($\textit{Phys. Fluids}$ $\textbf{17}$, 115102 (2005)). The resulting nonlocal and anisotropic eddy viscosities are exact and relate Reynolds stresses to mean velocity gradients at all locations. They can be used to inform RANS models of the sensitivity to the mean velocity gradient and the suitability of local and isotropic approximations. However, these brute-force approaches are expensive. They force the mean velocity gradient at each point in the averaged space and measure the Reynolds stress response, requiring a separate simulation for each mean velocity gradient location. Thus, computing the eddy viscosity requires as many simulations as degrees of freedom in the averaged space, which can be cost-prohibitive for problems with many degrees of freedom. In this work, we develop an adjoint-based MFM to obtain the eddy viscosity at a given Reynolds stress location using a single simulation. This approach recovers the Reynolds stress dependence at a location of interest, such as a separation point or near a wall, on the mean velocity gradient at all locations. We demonstrate using adjoint MFM to compute the eddy viscosity for a specified wall-normal location in an incompressible turbulent channel flow using one simulation. In contrast, a brute-force approach for the same problem requires $N=144$ simulations (the number of grid points in the non-averaged coordinate direction). We show that a local approximation for the eddy viscosity would have been inappropriate.

Efficient global resolvent analysis via the one-way Navier-Stokes equations. Part 2. Optimal response

Efficient global resolvent analysis via the one-way Navier-Stokes equations. Part 1. Forced response

Data-driven resolvent analysis

An invitation to resolvent analysis

Scalable resolvent analysis for three-dimensional flows

Efficient harmonic resolvent analysis via time-stepping

Resolvent analysis of a finite wing in transonic flow

Low-dimensional modelling of high-Reynolds-number shear flows incorporating constraints from the Navier–Stokes equation

Recursive one-way Navier-Stokes equations with PSE-like cost

Wavelet-based resolvent analysis of non-stationary flows

An Adjoint Method for a High-Order Discretization of Deforming Domain Conservation Laws for Optimization of Flow Problems

Targeted computation of nonlocal closure operators via an adjoint-based macroscopic forcing method

High-Resolution Viscous Terms Discretization and ILW Solid Wall Boundary Treatment for the Navier–Stokes Equations

Sparse space-time resolvent analysis for statistically-stationary and time-varying flows

On the reconstruction of unknown driving forces from low-mode observations in the 2D Navier-Stokes Equations

A sparsity-promoting resolvent analysis for the identification of spatiotemporally-localized amplification mechanisms

Identification of cross-frequency interactions in compressible cavity flow using harmonic resolvent analysis

Well-conditioned ultraspherical and spectral integration methods for resolvent analysis of channel flows of Newtonian and viscoelastic fluids

An Efficient Numerical Method for Coupling the RANS Equations with Spalart-Allmaras Turbulence Model Equation