2.1. Linearisation of the incompressible Navier–Stokes equations

The incompressible NS equations are given by

(2.1a)$$\begin{gather} \frac{\partial \boldsymbol{u}}{\partial t}+\boldsymbol{u}\boldsymbol{\cdot}\boldsymbol{\nabla}\boldsymbol{u} ={-}\boldsymbol{\nabla} p+\frac{1}{Re_\tau}\boldsymbol{\nabla}\boldsymbol{\cdot}(\boldsymbol{\nabla} \boldsymbol{u}+\boldsymbol{\nabla}\boldsymbol{u}^{\rm T}), \end{gather}$$

(2.1b)$$\begin{gather}\boldsymbol{\nabla}\boldsymbol{\cdot}\boldsymbol{u} = 0, \end{gather}$$

where $Re_\tau = {u_{\tau } h}/{\nu }$ is the friction Reynolds number, $u_{\tau }$ is the friction velocity, $h$ is the half-channel width, $\nu$ is the kinetic viscosity and the superscript ${}^{\rm T}$ denotes transpose. The linearised NS equations with respect to the fluctuation velocity $\boldsymbol {u}^{\prime }$ and pressure $p^{\prime }$ are expressed by

(2.2a)$$\begin{gather} \frac{\partial \boldsymbol{u}^\prime}{\partial t}+\boldsymbol{U} \boldsymbol{\cdot}\boldsymbol{\nabla}\boldsymbol{u}^\prime +\boldsymbol{u}^\prime\boldsymbol{\cdot}\boldsymbol{\nabla} \boldsymbol{U}+\boldsymbol{\nabla} p^\prime -\frac{1}{Re_\tau}\boldsymbol{\nabla}\boldsymbol{\cdot} (\boldsymbol{\nabla}\boldsymbol{u}^\prime+{\boldsymbol{\nabla}\boldsymbol{u}^\prime}^{\rm T}) = \boldsymbol{f}_0^{\prime}, \end{gather}$$

(2.2b)$$\begin{gather}\boldsymbol{\nabla}\boldsymbol{\cdot}\boldsymbol{u}^\prime = 0, \end{gather}$$

where ${}^\prime$ denotes the fluctuation variable, $U$ is the mean velocity and the nonlinear forcing term is defined as

(2.3)\begin{equation} \boldsymbol{f}_0^{\prime} ={-} \left( \boldsymbol{u}^\prime\boldsymbol{\cdot}\boldsymbol{\nabla}\boldsymbol{u}^\prime - \overline{\overline{\boldsymbol{u}^\prime\boldsymbol{\cdot}\boldsymbol{\nabla}\boldsymbol{u}^\prime}}\right) \end{equation}

with the double overlines denoting the averaged value along the temporal and uniform spatial directions.

On the other hand, the nonlinear term on the right-hand side of (2.2) can be further decomposed into two parts, i.e. the eddy-viscosity terms that are linear functions of $\boldsymbol {u}^{\prime }$ and the forcing $\boldsymbol {f}^{\prime }$ (Illingworth et al. Reference Illingworth, Monty and Marusic2018; Morra et al. Reference Morra, Semeraro, Henningson and Cossu2019; Towne et al. Reference Towne, Lozano-Durán and Yang2020). The linearised NS equations considering the eddy-viscosity model are expressed by

(2.4a)$$\begin{gather} \frac{\partial \boldsymbol{u}^\prime}{\partial t}+\boldsymbol{U}\boldsymbol{\cdot}\boldsymbol{\nabla}\boldsymbol{u}^\prime +\boldsymbol{u}^\prime\boldsymbol{\cdot}\boldsymbol{\nabla} \boldsymbol{U}+\boldsymbol{\nabla} p^\prime -\frac{1}{Re_\tau}\boldsymbol{\nabla}\boldsymbol{\cdot} \left[\frac{\nu_t+\nu}{\nu}(\boldsymbol{\nabla}\boldsymbol{u}^\prime+{\boldsymbol{\nabla}\boldsymbol{u}^\prime}^{\rm T})\right] = \boldsymbol{f}^{\prime}, \end{gather}$$

(2.4b)$$\begin{gather}\boldsymbol{\nabla}\boldsymbol{\cdot}\boldsymbol{u}^\prime = 0, \end{gather}$$

where the forcing term is defined as

(2.5)\begin{equation} \boldsymbol{f}^{\prime}= \boldsymbol{f}_0^{\prime} -\frac{1}{Re_\tau}\boldsymbol{\nabla}\boldsymbol{\cdot} \left[\frac{\nu_t}{\nu}(\boldsymbol{\nabla}\boldsymbol{u}^\prime+{\boldsymbol{\nabla}\boldsymbol{u}^\prime}^{\rm T})\right]. \end{equation}

The eddy viscosity $\nu _t$ can be calculated from the semi-empirical expression of Cess (Reference Cess1958) as reported by Reynolds & Hussain (Reference Reynolds and Hussain1972)

(2.6) \begin{align} \nu_t=\frac{\nu}{2} \left\{ 1+\frac{\kappa^2 Re_\tau ^2}{9} (2\tilde{y}-\tilde{y}^2)^2 (3-4\tilde{y}+2\tilde{y}^2)^2 \left[1-{\rm exp}\left(\frac{-Re_\tau (1-\left|1-\tilde{y}\right|)}{A}\right)\right]^2 \right\}^{1/2} -\frac{\nu}{2}, \end{align}

where $\tilde{y}=y/h$, the constants $\kappa =0.426$ and $A=25.4$. For the flow state estimation problems that are focused on in this study, the inclusion of the eddy-viscosity terms provides improved results (Illingworth et al. Reference Illingworth, Monty and Marusic2018; Amaral et al. Reference Amaral, Cavalieri, Martini, Jordan and Towne2021). Thus, we will also adopt the eddy-viscosity-modelled linearised formulations (2.4) in this study.

Taking the Fourier transformation to (2.4) in the homogeneous spatial directions, the linearised NS equations in each spatial scale $\boldsymbol {k}_{s}$ are obtained. For instance, in the fully developed turbulent channel flow, the Fourier transformation is taken in the streamwise $(x)$ and spanwise $(z)$ directions, where $\boldsymbol {k}_{s} = (k_x , k_z)$. The linearised NS equations are written in a discretised state-space form with $N_y$ points in the wall-normal direction as (Luhar, Sharma & McKeon Reference Luhar, Sharma and McKeon2014a,Reference Luhar, Sharma and McKeonb; Amaral et al. Reference Amaral, Cavalieri, Martini, Jordan and Towne2021)

(2.7a)$$\begin{gather} \boldsymbol{\mathsf{M}} \frac{\partial \hat{\boldsymbol{q}}_{\boldsymbol{k}_{s}}(t)}{\partial t}= \boldsymbol{\mathsf{A}}_{\boldsymbol{k}_{s}}\hat{\boldsymbol{q}}_{\boldsymbol{k}_{s}}(t) + \boldsymbol{\mathsf{B}} \hat{\boldsymbol{f}}_{\boldsymbol{k}_{s}}(t), \end{gather}$$

(2.7b)$$\begin{gather}\hat{\boldsymbol{m}}_{\boldsymbol{k}_{s}}(t) = \boldsymbol{\mathsf{C}} \hat{\boldsymbol{q}}_{\boldsymbol{k}_{s}}(t) + \hat{\boldsymbol{n}}_{\boldsymbol{k}_{s}}(t), \end{gather}$$

where $\hat {}$ denotes the variable after Fourier transformation, the state $\hat {\boldsymbol {q}}_{\boldsymbol {k}_{s}}(t) = [\hat {u} , \hat {v} , \hat {w} , \hat {p}]^{\rm T}$ and the operators are defined as

(2.8a–c)\begin{equation} \boldsymbol{\mathsf{M}} = \left[\begin{matrix} \boldsymbol{\mathsf{I}}_{3N_y} & 0 \\ 0 & 0 \end{matrix} \right], {\quad}\boldsymbol{\mathsf{A}}_{\boldsymbol{k}_{s}} = \left[\begin{matrix} \mathcal{L}_{\boldsymbol{k}_{s}} & -\boldsymbol{\nabla} ^{\rm T} \\ \boldsymbol{\nabla} & 0 \end{matrix} \right], {\quad} \boldsymbol{\mathsf{B}} = \left[ \begin{matrix} \boldsymbol{\mathsf{I}}_{3N_y} \\ 0 \end{matrix} \right], \end{equation}

where $\boldsymbol{\mathsf{C}}$ is the spatial observation tensor, $\hat {\boldsymbol {n}}_{\boldsymbol {k}_{s}}(t)$ is the measurement noise, $\boldsymbol {\nabla }$ is the Hamiltonian operator and $\mathcal {L}_{\boldsymbol {k}_{s}}$ is the spatial linearised NS operator. The expressions of $\mathcal {L}_{\boldsymbol {k}_{s}}$, $\boldsymbol {\nabla }$ and $\boldsymbol{\mathsf{C}}$ are provided in Appendix A. Further taking Fourier transformations along the statistically homogeneous temporal direction, the linearised NS equations at each spatio-temporal scale $\boldsymbol {k} = (\boldsymbol {k}_{s},\omega )$ are obtained, as expressed by

(2.9)\begin{equation} \hat{\boldsymbol{q}}_{\boldsymbol{k}} = \boldsymbol{\mathsf{H}}_{\boldsymbol{k}} \boldsymbol{\cdot} \hat{\boldsymbol{f}}_{\boldsymbol{k}}, \end{equation}

where $\boldsymbol{\mathsf{H}}_{\boldsymbol {k}}$ is the resolvent operator that is expressed as

(2.10)\begin{equation} \boldsymbol{\mathsf{H}}_{\boldsymbol{k}} = \left(-{\rm i}\omega \boldsymbol{\mathsf{M}} - \boldsymbol{\mathsf{A}}_{\boldsymbol{k}_{s}} \right)^{{-}1}, \end{equation}

with ${\rm i}=\sqrt {-1}$. From (2.9), $\hat {\boldsymbol {q}}_{\boldsymbol {k}}$ is regarded as the response of the input forcing $\hat {\boldsymbol {f}}_{\boldsymbol {k}}$ through the resolvent operator $\boldsymbol{\mathsf{H}}_{\boldsymbol {k}}$. Taking singular value decomposition of the transfer function $\boldsymbol{\mathsf{H}}_{\boldsymbol {k}}$ yields the resolvent input and output modes ordered by the gains, where the low-order modes can be utilised to investigate the predominant coherent structures assuming the low-rank behaviour, see, e.g. Moarref et al. (Reference Moarref, Sharma, Tropp and McKeon2013) and Sharma & McKeon (Reference Sharma and McKeon2013). Meanwhile, only the linearised relationship (2.9) built by the resolvent analysis will be utilised in this study.

Based on the resolvent formulation (2.9), the RBE that estimates the flow state from limited measurements was proposed by Towne et al. (Reference Towne, Lozano-Durán and Yang2020), and is further improved (Martini et al. Reference Martini, Cavalieri, Jordan, Towne and Lesshafft2020) by considering the forcing statistics. In the following, the RBE approach is further generalised to cases where the measurements are arbitrarily sampled.

2.2. Generalised RBE in cases with arbitrary temporal measurements

In this section, the case where the information of the turbulence field is partially available from measurements that are noise contaminated and arbitrarily sampled in time is considered. Before the estimation, Fourier transformations are conducted in statistically uniform spatial directions of the flow field, such as the streamwise and spanwise directions in fully developed turbulent channel flows, leading to the decomposed measurements at separated spatial scales $\boldsymbol {k}_{s}$. Meanwhile, since the flow signals can be unresolved in time, the Fourier transformation will not be conducted along the temporal direction. In this section, the subscript $_{\boldsymbol {k}_{s}}$ denoting the spatial scale of a flow quantity is omitted for brevity.

Consider the MIMO estimation case where the measured data are sampled at $N_{\check {T}}$ time instances with inadequate resolution. Meanwhile, to fully resolve the temporal information, $N_T$ snapshots are actually needed $(N_T > N_{\check {T}})$. The amounts of response and observation of different types or at different locations at a given spatial scale $\boldsymbol {k}_{s}$ are denoted as $N_{q}$ and $N_{m}$, respectively. Without loss of generality, a sequence of measured data $m(j)$, where $j \in [1,N_m N_{\check {T}}]$, is taken into consideration. The linear estimation of the flow state $q(i)$, where $i \in [1,N_q N_{\check {T}}]$ is based on all the available measurements $m(j)$, can be expressed by the convolution operation as

(2.11) \begin{equation} q_{e}(i) = \sum_{j=1}^{N_m N_{\check{T}}} \mathscr{h}(i,j) m(j), \end{equation}

where the subscript ${}_{e}$ denotes the estimated quantity and $\mathscr {h}(i,j)$ is the convolution kernel. Equation (2.11) can also be expressed in matrix form by

(2.12)\begin{equation} \boldsymbol{q}_{e} = \boldsymbol{\mathsf{h}} \boldsymbol{\cdot} \boldsymbol{m}, \end{equation}

where $\boldsymbol{\mathsf{h}}$ is a two-dimensional matrix consisting of $\mathscr {h}(i,j)$ at the $(i,j)$ location.

Denoting the estimation error $\varepsilon (i) = q(i) - q_{{e}}(i)$, the cost function $\mathcal {J}(i)$ is defined as the MSE that is expressed by

(2.13)\begin{equation} \mathcal{J}(i)=\left\langle \left| \varepsilon(i) \right|^2\right\rangle = \left\langle \left[ q(i) - \sum_{j=1}^{N_m N_{\check{T}}} \mathscr{h}(i,j) m(j) \right] {\cdot} \overline{\left[ q(i) - \sum_{k=1}^{N_m N_{\check{T}}} \mathscr{h}(i,k) m(k) \right]} \right\rangle, \end{equation}

where the overline denotes the complex conjugate and $\langle \rangle$ denotes ensemble average. To minimise the MSE of the estimated signal with optimised transfer kernel $\mathscr {h}(i,j) = \mathscr {h}^{R}(i,j) + {\rm i} \mathscr {h}^{I}(i,j)$, the partial derivatives of the cost function $\mathcal {J}(i)$ to any $\mathscr {h}^{R}(i,j)$ or $\mathscr {h}^{I}(i,j)$ with $j \in [1,N_m N_{\check {T}}]$ should be zero, i.e.

(2.14)\begin{equation} \left.\begin{gathered} \frac{\partial \mathcal{J} (i) }{\partial \mathscr{h}^{R(I)}(i,j)} = \left\langle -2 \left[ q(i)-\sum_{j=1}^{N_m N_{\check{T}}} \mathscr{h}(i,k) m(k) \right] \overline{ m(j) } \right\rangle_{R(I)} = 0,\\ \Rightarrow \quad \left\langle q(i) \overline{m(j)} \right\rangle - \sum_{j=1}^{N_m N_{\check{T}}} \mathscr{h}(i,k) \left\langle m(k) \overline{m(j)} \right\rangle_{R(I)} = 0. \end{gathered}\right\} \end{equation}

The optimal transfer kernel can be obtained by solving the system of linear equations (2.14). The solution can be expressed as

(2.15)\begin{equation} \boldsymbol{\mathsf{h}} = \left\langle \boldsymbol{q} \boldsymbol{m}^{{\ast}} \right\rangle \boldsymbol{\cdot} \left\langle \boldsymbol{m} \boldsymbol{m}^{{\ast}} \right\rangle^{{-}1}, \end{equation}

where the superscript ${}^{\ast }$ denotes the Hermitian transpose.

Equation (2.15) provides the expression of the optimal linear transfer kernel in terms of the covariance matrices of the flow state and measured signal. However, since such spatio-temporal covariance matrices are not directly available, the expressions of the optimal transfer kernel as a function of the forcing CSD tensors, which can be provided by the forcing models, should be further derived.

Considering the resolvent formulation, the expressions of $\boldsymbol {q}$ and $\boldsymbol {m}$ in the physical space as functions of the forcing $\hat {\boldsymbol {f}}_{j}$ at $\omega _j$, where $j \in [1,N_T]$, in the Fourier space are expressed as

(2.16a)$$\begin{gather} \boldsymbol{q} = \frac{1}{N_T}\sum_{j=1}^{N_T} \check{\zeta}_{\boldsymbol{q},j} \boldsymbol{\cdot} \hat{\boldsymbol{q}}_{j} = \frac{1}{N_T}\sum_{j=1}^{N_T} \check{\zeta}_{\boldsymbol{q},j} \boldsymbol{\cdot} \left( \boldsymbol{\mathsf{H}}_j \boldsymbol{\cdot} \hat{\boldsymbol{f}}_{j}\right) , \end{gather}$$

(2.16b)$$\begin{gather}\boldsymbol{m} = \frac{1}{N_T}\sum_{j=1}^{N_T} \check{\zeta}_{\boldsymbol{m},j} \boldsymbol{\cdot} \hat{\boldsymbol{m}}_{j} = \frac{1}{N_T} \sum_{j=1}^{N_T} \check{\zeta}_{\boldsymbol{m},j} \boldsymbol{\cdot} \left( \boldsymbol{\mathsf{C}} \boldsymbol{\cdot} \boldsymbol{\mathsf{H}}_j \boldsymbol{\cdot} \hat{\boldsymbol{f}}_{j} + \hat{\boldsymbol{n}}_{j} \right) , \end{gather}$$

where $\hat {\boldsymbol {q}}_{j}$, $\hat {\boldsymbol {m}}_{j}$ and $\hat {\boldsymbol {n}}_{j}$ are the Fourier transforms of response, measurement and sensor noise at $\omega _j$. Here, $\check {\zeta }_{\boldsymbol {q}(\boldsymbol {m}),j} = \textrm {diag} [ \overbrace {\check {\xi }_{j}, \ldots, \check {\xi }_{j} }^{N_{q(m)}} ]$, $\check {\xi }_{j} = \textrm {e}^{\textrm {i} \omega _j \check {\boldsymbol {t}}}$, $\check {\boldsymbol {t}} \in \mathbb {R}^{N_{\check {T}}}$ is the vector consisting of the sampled time instances.

Algorithm 1 GRBE

Using (2.16a,b), the spatio-temporal covariance matrices can be predicted with the forcing models by

(2.17a)$$\begin{gather} \left\langle \boldsymbol{q} \boldsymbol{m}^{{\ast}} \right\rangle_{p} = \sum_j \check{\zeta}_{\boldsymbol{q},j} \boldsymbol{\mathsf{S}}_{\boldsymbol{q}\boldsymbol{m},j} \check{\zeta}_{\boldsymbol{m},j}^{{\ast}} = \sum_j \check{\zeta}_{\boldsymbol{q},j} \boldsymbol{\mathsf{H}}_{j} \boldsymbol{\mathsf{S}}_{\boldsymbol{f}\boldsymbol{f},j} \boldsymbol{\mathsf{H}}_{\boldsymbol{m},j}^{{\ast}} \check{\zeta}_{\boldsymbol{m},j}^{{\ast}}, \end{gather}$$

(2.17b)$$\begin{gather}\left\langle \boldsymbol{m} \boldsymbol{m}^{{\ast}} \right\rangle_{p} = \sum_j \check{\zeta}_{\boldsymbol{m},j} \boldsymbol{\mathsf{S}}_{\boldsymbol{m}\boldsymbol{m},j} \check{\zeta}_{\boldsymbol{m},j}^{{\ast}} = \sum_j \check{\zeta}_{\boldsymbol{m},j} \left( \boldsymbol{\mathsf{H}}_{\boldsymbol{m},j} \boldsymbol{\mathsf{S}}_{\boldsymbol{f}\boldsymbol{f},j} \boldsymbol{\mathsf{H}}_{\boldsymbol{m},j}^{{\ast}} + \boldsymbol{\mathsf{S}}_{\boldsymbol{n}\boldsymbol{n},j} \right) \check{\zeta}_{\boldsymbol{m},j}^{{\ast}}, \end{gather}$$

where $\boldsymbol{\mathsf{S}}_{\boldsymbol {q}\boldsymbol {m},j}=\langle \hat {\boldsymbol {q}}_{j} \hat {\boldsymbol {m}}_{j}^{\ast }\rangle$, $\boldsymbol{\mathsf{S}}_{\boldsymbol {m}\boldsymbol {m},j}=\langle \hat {\boldsymbol {m}}_{j} \hat {\boldsymbol {m}}_{j}^{\ast }\rangle$, $\boldsymbol{\mathsf{S}}_{\boldsymbol {f}\boldsymbol {f},j}=\langle \hat {\boldsymbol {f}}_{j} \hat {\boldsymbol {f}}_{j}^{\ast }\rangle$ and $\boldsymbol{\mathsf{S}}_{\boldsymbol {n}\boldsymbol {n},j}=\langle \hat {\boldsymbol {n}}_{j} \hat {\boldsymbol {n}}_{j}^{\ast }\rangle$. Substituting (2.17a,b) into (2.15), the optimal transfer kernel that considers the forcing and noise statistics is expressed by

(2.18a)\begin{align} \boldsymbol{\mathsf{h}}&=\left\langle \boldsymbol{q} \boldsymbol{m}^{{\ast}} \right\rangle_{p} \boldsymbol{\cdot} \left\langle \boldsymbol{m} \boldsymbol{m}^{{\ast}} \right\rangle_{p}^{{-}1} \end{align}

(2.18b)\begin{align} &=\left( \sum_j \check{\zeta}_{\boldsymbol{q},j} \boldsymbol{\mathsf{S}}_{\boldsymbol{q}\boldsymbol{m},j} \check{\zeta}_{\boldsymbol{m},j}^{{\ast}} \right) \boldsymbol{\cdot} \left( \sum_j \check{\zeta}_{\boldsymbol{m},j} \boldsymbol{\mathsf{S}}_{\boldsymbol{m}\boldsymbol{m},j} \check{\zeta}_{\boldsymbol{m},j}^{{\ast}}\right) ^{{-}1} \end{align}

(2.18c)\begin{align} &=\left[ \sum_j \check{\zeta}_{\boldsymbol{q},j} \boldsymbol{\mathsf{H}}_{j} \boldsymbol{\mathsf{S}}_{\boldsymbol{f}\boldsymbol{f},j} \boldsymbol{\mathsf{H}}_{\boldsymbol{m},j}^{{\ast}} \check{\zeta}_{\boldsymbol{m},j}^{{\ast}} \right] \boldsymbol{\cdot} \left[ \sum_j \check{\zeta}_{\boldsymbol{m},j} \left( \boldsymbol{\mathsf{H}}_{\boldsymbol{m},j} \boldsymbol{\mathsf{S}}_{\boldsymbol{f}\boldsymbol{f},j} \boldsymbol{\mathsf{H}}_{\boldsymbol{m},j}^{{\ast}} + \boldsymbol{\mathsf{S}}_{\boldsymbol{n}\boldsymbol{n},j} \right) \check{\zeta}_{\boldsymbol{m},j}^{{\ast}}\right] ^{{-}1}. \end{align}

From the above derivations, the transfer kernel $\boldsymbol{\mathsf{h}}$ is the optimal linear estimator of the response $\boldsymbol {q}$ when the real statistics $\boldsymbol{\mathsf{S}}_{\boldsymbol {ff},j}$ and $\boldsymbol{\mathsf{S}}_{\boldsymbol {nn},j}$ are incorporated, which can be obtained from another existing database. More importantly, the transfer kernel defined in (2.18) also constructs an open framework that enables the incorporation of the spatio-temporal forcing models that predict the forcing CSD tensors all over the frequency domain at each spatial scale $\boldsymbol {k}_{s}$ when the real statistics of forcing are unknown. The current models that provide spatio-temporal prediction of the forcing include the baseline model (B-model) (Hwang & Cossu Reference Hwang and Cossu2010; Madhusudanan et al. Reference Madhusudanan, Illingworth and Marusic2019), W-model and $\lambda$-model (Gupta et al. Reference Gupta, Madhusudanan, Wan, Illingworth and Juniper2021), which assume that the forcing is white in time and uncorrelated in space with certain forms of spatial energy distributions along the wall-normal height. These forcing models are reviewed in Appendix B. The new approach generalises the RBE to the cases with arbitrary sampled measurements in time, and thus it is denoted as the GRBE in this study. The procedure of the GRBE is provided in Algorithm 1.

2.3. Expected estimation error from the GRBE

In § 2.2, the theoretically optimal transfer kernel of arbitrarily sampled measurements is derived in (2.18). On the other hand, the expected statistics of the estimation error can also be theoretically derived. Recall (2.13) that describes the linear relationship between the observations and estimations, the matrix form holds

(2.19)\begin{align} \left\langle \boldsymbol{\varepsilon} \boldsymbol{\varepsilon}^{{\ast}} \right\rangle & = \left\langle \left( \boldsymbol{q} - \boldsymbol{\mathsf{h}} \boldsymbol{\cdot} \boldsymbol{m} \right) \boldsymbol{\cdot} \left( \boldsymbol{q} - \boldsymbol{\mathsf{h}} \boldsymbol{\cdot} \boldsymbol{m} \right)^{{\ast}} \right\rangle \nonumber\\ & = \left\langle \boldsymbol{q} \boldsymbol{q}^{{\ast}} \right\rangle + \boldsymbol{\mathsf{h}} \boldsymbol{\cdot} \left\langle \boldsymbol{m} \boldsymbol{m}^{{\ast}} \right\rangle \boldsymbol{\cdot} \boldsymbol{\mathsf{h}}^{{\ast}} - \left\langle \boldsymbol{q} \boldsymbol{m}^{{\ast}} \right\rangle \boldsymbol{\cdot} \boldsymbol{\mathsf{h}}^{{\ast}} - \boldsymbol{\mathsf{h}} \boldsymbol{\cdot} \left\langle \boldsymbol{m} \boldsymbol{q}^{{\ast}} \right\rangle . \end{align}

When the actual forcing and noise statistics are incorporated into the transfer kernel $\boldsymbol{\mathsf{h}}$, the minimum expected estimation error can be obtained, as calculated by

(2.20)\begin{align} \left\langle \boldsymbol{\varepsilon} \boldsymbol{\varepsilon}^{{\ast}} \right\rangle_{OPT} & = \left\langle \boldsymbol{q} \boldsymbol{q}^{{\ast}} \right\rangle + \left\langle {\boldsymbol{q}} {\boldsymbol{m}}^{ {\ast}} \right\rangle \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle^{{-}1} \boldsymbol{\cdot} \left\langle \boldsymbol{m} \boldsymbol{m}^{{\ast}} \right\rangle \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle^{{-}1} \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{q}}^{ {\ast}} \right\rangle\nonumber\\ &\quad- \left\langle \boldsymbol{q} \boldsymbol{m}^{{\ast}} \right\rangle \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle^{{-}1} \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{q}}^{ {\ast}} \right\rangle - \left\langle {\boldsymbol{q}} {\boldsymbol{m}}^{ {\ast}} \right\rangle \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle^{{-}1} \boldsymbol{\cdot} \left\langle \boldsymbol{m} \boldsymbol{q}^{{\ast}} \right\rangle \nonumber\\ & = \left\langle \boldsymbol{q} \boldsymbol{q}^{{\ast}} \right\rangle - \left\langle \boldsymbol{q} \boldsymbol{m}^{{\ast}} \right\rangle \left\langle \boldsymbol{m} \boldsymbol{m}^{{\ast}} \right\rangle^{{-}1} \left\langle \boldsymbol{m} \boldsymbol{q}^{{\ast}} \right\rangle, \end{align}

where the subscript ‘$_{OPT}$’ denotes the variable from the optimal linear estimation. Considering that $\langle \boldsymbol {q} \boldsymbol {m}^{\ast } \rangle \langle \boldsymbol {m} \boldsymbol {m}^{\ast } \rangle ^{-1} \langle \boldsymbol {m} \boldsymbol {q}^{\ast } \rangle = \langle \boldsymbol {q}_{e} \boldsymbol {q}_{e}^{\ast } \rangle$ according to (2.12) and (2.18a), (2.20) indicates the orthogonality principle for the optimal linear estimation, i.e. $\langle \boldsymbol {q}_{e} \boldsymbol {\cdot } \boldsymbol {\varepsilon }^{\ast } \rangle = \boldsymbol{\mathsf{O}}$. Considering (2.16), the covariance matrix of the estimated error is further expressed by

(2.21)\begin{align} \left\langle \boldsymbol{\varepsilon} \boldsymbol{\varepsilon}^{{\ast}} \right\rangle_{OPT} &= \left[ \sum_j \check{\zeta}_{\boldsymbol{q},j} \boldsymbol{\mathsf{H}}_{j} \boldsymbol{\mathsf{S}}_{\boldsymbol{f}\boldsymbol{f},\boldsymbol{j}} \boldsymbol{\mathsf{H}}_{j}^{{\ast}} \check{\zeta}_{\boldsymbol{q},j}^{{\ast}} \right] - \left[ \sum_j \check{\zeta}_{\boldsymbol{q},j} \boldsymbol{\mathsf{H}}_{j} \boldsymbol{\mathsf{S}}_{\boldsymbol{f}\boldsymbol{f},\boldsymbol{j}} \boldsymbol{\mathsf{H}}_{\boldsymbol{m},j}^{{\ast}} \check{\zeta}_{\boldsymbol{m},j}^{{\ast}} \right]\nonumber\\ &\quad \boldsymbol{\cdot} \left[\sum_j \check{\zeta}_{\boldsymbol{m},j} \left( \boldsymbol{\mathsf{H}}_{\boldsymbol{m},j} \boldsymbol{\mathsf{S}}_{\boldsymbol{f}\boldsymbol{f},\boldsymbol{j}} \boldsymbol{\mathsf{H}}_{\boldsymbol{m},j}^{{\ast}} + \boldsymbol{\mathsf{S}}_{\boldsymbol{n}\boldsymbol{n},\boldsymbol{j}} \right) \check{\zeta}_{\boldsymbol{m},j}^{{\ast}} \right]^{{-}1} \boldsymbol{\cdot} \left[ \sum_j \check{\zeta}_{\boldsymbol{q},j} \boldsymbol{\mathsf{H}}_{j} \boldsymbol{\mathsf{S}}_{\boldsymbol{f}\boldsymbol{f},\boldsymbol{j}} \boldsymbol{\mathsf{H}}_{\boldsymbol{m},j}^{{\ast}} \check{\zeta}_{\boldsymbol{m},j}^{{\ast}} \right]. \end{align}

On the other hand, when the forcing models rather than the real flow statistics are incorporated in the GRBE estimator, the orthogonal principle as in (2.20) for the optimal linear estimation no longer holds. Instead, the expected error energy from the modelled GRBE estimator should be directly derived from (2.19) by considering (2.18a), i.e.

(2.22)\begin{align} \left\langle {\boldsymbol{\varepsilon}} {\boldsymbol{\varepsilon}}^{{\ast}} \right\rangle_{MOD} & = \left\langle \boldsymbol{q} \boldsymbol{q}^{{\ast}} \right\rangle + \boldsymbol{\mathsf{h}}_{\boldsymbol{q},{MOD}} \boldsymbol{\cdot} \left\langle \boldsymbol{m} \boldsymbol{m}^{{\ast}} \right\rangle \boldsymbol{\cdot} \boldsymbol{\mathsf{h}}_{\boldsymbol{q},{MOD}}^{{\ast}} - \left\langle \boldsymbol{q} \boldsymbol{m}^{{\ast}} \right\rangle \boldsymbol{\cdot} \boldsymbol{\mathsf{h}}_{\boldsymbol{q},{MOD}}^{{\ast}} - \boldsymbol{\mathsf{h}}_{\boldsymbol{q},{ MOD}} \boldsymbol{\cdot} \left\langle \boldsymbol{m} \boldsymbol{q}^{{\ast}} \right\rangle\nonumber\\ & = \left\langle \boldsymbol{q} \boldsymbol{q}^{{\ast}} \right\rangle + \left\langle {\boldsymbol{q}} {\boldsymbol{m}}^{ {\ast}} \right\rangle_{p} \boldsymbol{\cdot}\left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle_{p}^{{-}1} \boldsymbol{\cdot} \left\langle \boldsymbol{m} \boldsymbol{m}^{{\ast}} \right\rangle \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle_{p}^{{-}1} \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{q}}^{ {\ast}} \right\rangle_{p}\nonumber\\ &\quad - \left\langle \boldsymbol{q} \boldsymbol{m}^{{\ast}} \right\rangle \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle_{p}^{{-}1} \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{q}}^{ {\ast}} \right\rangle_{p} - \left\langle {\boldsymbol{q}} {\boldsymbol{m}}^{ {\ast}} \right\rangle_{p} \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle_{p}^{{-}1} \boldsymbol{\cdot} \left\langle \boldsymbol{m} \boldsymbol{q}^{{\ast}} \right\rangle, \end{align}

where $\langle {\boldsymbol {q}} {\boldsymbol {m}}^{ \ast } \rangle _{p}$ and $\langle {\boldsymbol {m}} {\boldsymbol {m}}^{ \ast } \rangle _{p}$ are defined in (2.17) and the subscript ‘$_{MOD}$’ denotes the variable corresponding to the linear estimation informed by the modelled forcing. Substituting the estimation error from optimal linear estimation in (2.20) to (2.22), the estimation error from the modelled forcing can be expressed as

(2.23)\begin{align} \left\langle {\boldsymbol{\varepsilon}} {\boldsymbol{\varepsilon}}^{{\ast}} \right\rangle_{MOD} & = \left\langle {\boldsymbol{\varepsilon}} {\boldsymbol{\varepsilon}}^{{\ast}} \right\rangle_{OPT} + \left[ \left\langle {\boldsymbol{q}} {\boldsymbol{m}}^{ {\ast}} \right\rangle \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle^{{-}1} - \left\langle {\boldsymbol{q}} {\boldsymbol{m}}^{ {\ast}} \right\rangle_{p} \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle_{p}^{{-}1} \right] \nonumber\\ &\quad \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle \boldsymbol{\cdot} \left[ \left\langle {\boldsymbol{q}} {\boldsymbol{m}}^{ {\ast}} \right\rangle \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle^{{-}1} - \left\langle {\boldsymbol{q}} {\boldsymbol{m}}^{ {\ast}} \right\rangle_{p} \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle_{p}^{{-}1} \right]^{{\ast}}\nonumber\\ & = \left\langle {\boldsymbol{\varepsilon}} {\boldsymbol{\varepsilon}}^{{\ast}} \right\rangle_{OPT} + \left( \boldsymbol{\mathsf{h}}_{\boldsymbol{q},{OPT}} - \boldsymbol{\mathsf{h}}_{\boldsymbol{q},{MOD}} \right) \boldsymbol{\cdot} \left\langle {\boldsymbol{m}} {\boldsymbol{m}}^{{\ast}} \right\rangle \boldsymbol{\cdot} \left( \boldsymbol{\mathsf{h}}_{\boldsymbol{q},{OPT}} - \boldsymbol{\mathsf{h}}_{\boldsymbol{q},{MOD}} \right)^{{\ast}}, \end{align}

where $\boldsymbol{\mathsf{h}}_{\boldsymbol {q},{OPT}}=\langle {\boldsymbol {q}} {\boldsymbol {m}}^{ \ast } \rangle \boldsymbol {\cdot } \langle {\boldsymbol {m}} {\boldsymbol {m}}^{\ast } \rangle ^{-1}$. On the right-hand side of (2.23), the extra term $[ ( \boldsymbol{\mathsf{h}}_{\boldsymbol {q},{OPT}} - \boldsymbol{\mathsf{h}}_{\boldsymbol {q},{MOD}} ) \boldsymbol {\cdot } \langle {\boldsymbol {m}} {\boldsymbol {m}}^{\ast } \rangle \boldsymbol {\cdot } ( \boldsymbol{\mathsf{h}}_{\boldsymbol {q},{OPT}} - \boldsymbol{\mathsf{h}}_{\boldsymbol {q},{MOD}} )^{\ast } ]$ is positive semi-definite. Thus, the estimation error energy from the modelled forcing will be minimised when $\boldsymbol{\mathsf{h}}_{\boldsymbol {q},{MOD}} = \boldsymbol{\mathsf{h}}_{\boldsymbol {q},{OPT}}$.

3.1. The complex Ginzburg–Landau equation

The complex Ginzburg–Landau equation is a one-dimensional model that mimics a real turbulence system (Bagheri, Brandt & Henningson Reference Bagheri, Brandt and Henningson2009; Towne, Schmidt & Colonius Reference Towne, Schmidt and Colonius2018; Towne et al. Reference Towne, Lozano-Durán and Yang2020) and facilitates the analysis of the properties of the linear operator. The equation can be written in the same form as (2.7a), with the linear operators $\boldsymbol{\mathsf{A}}$, $\boldsymbol{\mathsf{B}}$ and $\boldsymbol{\mathsf{M}}$ expressed by

(3.1)\begin{equation} \left.\begin{gathered} \boldsymbol{\mathsf{A}} ={-} \nu \frac{\partial}{\partial x} + \gamma \frac{{\partial}^2}{\partial x^2} + \mu (x), \\ \boldsymbol{\mathsf{M}} = \boldsymbol{\mathsf{B}} = \boldsymbol{\mathsf{I}}. \end{gathered}\right\} \end{equation}

The parameters in (3.1) are set as $\mu (x) = (\mu _0 - c_u^2) + \mu _2 x^2 / 2$, $\nu = 2+0.2\textrm {i}$, $\gamma = 1-\textrm {i}$, $\nu _0 = 0.3$, $c_u = 0.2$ and $\mu _2 = -0.01$, which are identical to those in Towne et al. (Reference Towne, Schmidt and Colonius2018). The spatial domain is discretised with second-order upwind scheme at evenly distributed grid nodes with $N_x = 221$ in a computational domain $x \in [ -85,85 ]$, where the grid size $\Delta x = 0.77$. Following Towne et al. (Reference Towne, Schmidt and Colonius2018), the input forcing that excites the system is realised by setting the value at each discrete time step (with a constant step size of 0.05) and grid node to be a random complex number with uniformly distributed phase from 0 to $2 {\rm \pi}$ and normally distributed amplitude with unit variance. The forcing is then low-pass filtered using a tenth-order finite-impulse-response filter with a cutoff frequency equal to $60\,\%$ of the Nyquist frequency $\omega _{s}$ equal to $20 {\rm \pi}$. The spatial scope of forcing is limited in the interior of the computational domain with an exponential envelope of $\textrm {exp} [ -( x/L )^p ]$, where $L = 60$ and $p = 10$. The equations are integrated in time with the Crank–Nicolson scheme (Crank & Nicolson Reference Crank and Nicolson1947). After the computational result achieves a statistical stationary state, the flow data at the measurement location are sampled for estimation.

3.1.1. Non-causal estimation results

In this section, three cases in terms of different sampling time intervals are considered in this section, namely the temporally resolved case (TR), the randomly sampled case (RS) and the temporally unresolved case (TU), with $\Delta t = 0.5$, $\Delta t = \textrm {rand}(0.5,8.0)$ that is uniformly distributed between 0.5 and 8.0, and $\Delta t = 8.0$, respectively. The same sampling period of $T = 200$ is applied for each case. To construct the estimator of the GRBE in (2.18) and RBE in (2.24), the incorporated CSD tensor of the forcing is set to be $\boldsymbol{\mathsf{S}}_{\boldsymbol {ff}} = \boldsymbol{\mathsf{I}}$ at each frequency.

In the TR case with $\Delta t = 0.5$, the temporal information of the measured signals is resolved in time. In this case, the flow information at the energy-containing temporal scales can be well recovered by the Fourier transform. According to the derivations in § 2.4, the RBE results should be the same as the GRBE in this case. Figure 1(a,d) depicts the estimation of the flow data at $x = 10$ using the measurements at $x = 0$. It is observed that the prediction results of the RBE are the same as those of the GRBE, which match pretty well with the reference numerical results. To quantitatively analyse the dependence of estimation accuracy on the distance between the measured and predicted signals, the relative MSE at the $i$th location is defined as

(3.2)\begin{equation} {\rm Err} = \frac{\left\langle\sum_{j=1}^{N_{\check{T}}} \left| \varepsilon(i,j) \right| ^2 \right\rangle}{\left\langle \sum_{j=1}^{N_{\check{T}}} \left| u(i,j) \right| ^2 \right\rangle}, \end{equation}

where $\varepsilon (i,j) = u(i,j)-u_{e}(i,j)$ is the error at the $i$th location and $j$th time instance. Figure 2(a) shows the relative MSEs of the estimations at different spatial locations in the TR case. The relative MSEs from the RBE and GRBE overlap with each other at all of the considered locations.

Figure 1. Estimation of the flow states at $x = 10$ using the measurements at $x = 0$ with the GRBE (a–c) and RBE (d–f) in cases TR (a,d), RS (b,e) and TU (c,f). The solid and dashed curves denote the reference numerical results and the estimated ones, respectively.

Figure 2. Relative MSEs of the estimations with measurements at $x=0$ in cases TR ($a$), RS ($b$) and TU ($c$). The translucent grey curves denote the spatial distribution of the relative fluctuation energy from the reference numerical simulation.

In the RS case, the sampling time steps are set to be uniformly distributed in the range of $0.5 \leqslant t \leqslant 8.0$. From figure 1(b,e), the GRBE results match well with the reference numerical results, while those from the RBE deviate obviously from the reference data. In figure 2(b), the relative MSEs of the RBE increase rapidly when moving away from the measurement location, while those of the GRBE are much lower at corresponding locations. Similar phenomena are observed in the TU case, where the estimation accuracy of the GRBE is also much lower than the RBE.

3.1.2. Causal estimation results

In Martini et al. (Reference Martini, Cavalieri, Jordan, Towne and Lesshafft2020), the casual RBE estimator is proposed by truncating the resolvent-based kernel (2.26), which can be explicitly expressed in (2.30) with the weighting function $w(\tau )=0$ when $\tau > t$, where $t$ is the estimation time step. On the other hand, the causal GRBE estimator can be realised by setting $\check {\boldsymbol {t}}$ containing only the time instances with $\tau = t$ and $\tau \leqslant t$ when constructing the Fourier bases $\check {\zeta }_{\boldsymbol {q},j}$ and $\check {\zeta }_{\boldsymbol {m},j}$ in (2.18), respectively.

Figures 3 and 4 show the estimated instantaneous flow states and relative MSEs of the one-dimensional Ginzburg–Landau equation with measurements located at $x=0$. The sampling time interval is set as $\Delta t = 0.5$, which is the same as that in the TR case. Since the $H_2$-optimal estimation provides optimal causal results with white-in-time forcing (Doyle et al. Reference Doyle, Glover, Khargonekar and Francis1989; Zhou & Doyle Reference Zhou and Doyle1998), which is met in this case, its results are used for the validation of the causal GRBE herein. From figures 3 and 4, the causal GRBE results overlap with the $H_2$-optimal estimation results at all of the tested locations. On the other hand, the estimation errors from the causal RBE results are obviously larger than those from the causal GRBE results at $x < 0$. Even in the downstream region with $x \geqslant 0$, the estimation errors of the causal RBE estimator are still subtly larger than those of the causal GRBE estimator, as in figure 4(b). Such observations are consistent with those in Martini et al. (Reference Martini, Cavalieri, Jordan, Towne and Lesshafft2020), where the causal RBE results are found to underestimate the flow state compared with the causal Kalman filter results.

Figure 3. Causal estimation of the flow states at $x = -1.5$ (a,c,e) and $x=10$ (b,d,f) using the measurements at $x = 0$ with the $H_2$-optimal estimation (a,b), GRBE (c,d) and RBE (e,f). The solid and dashed curves denote the reference numerical results and the estimated ones, respectively.

Figure 4. Relative MSEs of the causal estimations with measurements at $x=0$ ($a$), with zoomed results near $x=0$ ($b$). The translucent grey curves denote the spatial distribution of the relative fluctuation energy from the reference numerical simulation.

3.2. Estimation of the turbulent channel flows

In this section, validations of the GRBE in estimating the turbulent channel flows with $Re_{\tau } = 186$, 547 and 934 will be conducted, where the RBE results are also presented.

3.2.1. Descriptions of the DNS database

The code that is used to generate the widely used open-source DNS database (Del Alamo & Jiménez Reference Del Alamo and Jiménez2003; Hoyas & Jiménez Reference Hoyas and Jiménez2008) for incompressible turbulent channel flows with $Re_{\tau}=186$, 547 and 934 is applied to compute the DNS results in this study. Details of the DNS set-ups are listed in table 1. To provide temporally resolved results, the sampling time intervals $\delta t^+=\delta t {\cdot } u_{\tau }^2 / \nu$ are set as 2.13, 4.81 and 5.09 for cases with $Re_{\tau } = 186$, 547 and 934, respectively. The normalised total simulation time $(u_\tau {\cdot } T)/h$ is larger than 5.0 in each case to obtain statistically converged results. To process the DNS datasets, the computational domain is divided into blocks with a spatial size of ${L_x}/h = 4{\rm \pi}$, ${L_z}/h = {\rm \pi}$ and ${L_y}/h = 2$ in each case. Considering that the fully developed channel flow is statistically symmetric about the centreline $y=h$, the DNS data that are mirrored about the centreline are treated as another set of blocks besides the original one. The time periods of the blocks are set as $240\delta t$, $160\delta t$ and $160\delta t$ for cases with $Re_{\tau } = 186$, 547 and 934, respectively, with $75\,\%$ overlap in the temporal direction. The rectangular window function is utilised for spectral analysis. With the above set-ups for data processing, the numbers of blocks are 408, 352 and 252 for cases with $Re_{\tau } = 186$, 547 and 934, respectively. The DNS datasets utilised in this study have been validated in Ying et al. (Reference Ying, Liang, Li and Fu2023) by comparisons with the time-averaged and r.m.s. profiles of the velocities from the open-source DNS database (Del Alamo & Jiménez Reference Del Alamo and Jiménez2003; Hoyas & Jiménez Reference Hoyas and Jiménez2008). In the following, the DNS results with $Re_{\tau } = 934$ will be mainly discussed, supplemented by the investigations on the estimation results with $Re_{\tau } = 186$ and 547.

Table 1. Parameters of the incompressible channel DNS set-ups.

3.2.2. Comparisons of the estimation results from the GRBE and RBE

Like the case settings in the validations with the Ginzburg–Landau equation, three cases are set with different temporal resolutions, namely the TR, RS and TU cases. In the TR case, the sampled data are temporally resolved with $\Delta t^+ = \delta t^+$, which corresponds to the case where the measured data are finely sampled in time and well restored. In the RS case, the sampling time intervals are uniformly distributed between $\Delta t^+ = \delta t^+$ and $32 \delta t^+$, as shown in figure 5. In this case, the temporal information is partially utilised, which corresponds to the case where the initially time-resolved measured data are missing arbitrarily. As for the TU case, the sampling interval is set as $\Delta t^+ = 32 \delta t^+$. Besides the sampling temporal intervals, the types of observed quantities may also influence the estimation accuracy. Meanwhile, in this section, we only show the estimation results from the observations of both the wall shear stress and pressure at the wall for brevity. On the other hand, detailed discussions on the impacts of choices of measurement types and the inclusion of forcing models on the estimation accuracy of the GRBE will be presented in the following § 4. The estimators of both RBE and GRBE are informed by the real flow statistics. With the above settings, the GRBE and RBE estimators are informed by abundant information, including all the considered wall quantities, i.e. wall shear stress and pressure, and exact flow statistics. Such choice enables us to adequately compare the performance of the two methods with given sampling time intervals of the measurement while excluding the influence from other factors.

Figure 5. Sampling time intervals in the RS case of turbulent channel flow with $Re_{\tau }=934$.

Intuitive insights into the estimation results can be found in the instantaneous flow field at $y^+ = 100$ with $Re_{\tau } = 934$ shown in figure 6. It is observed that the estimation results from the GRBE and RBE are the same in the TR case. The estimated results from both approaches recover well the representative flow structures in the DNS results. With the increase of the sampling intervals, the estimation accuracy of the RBE rapidly deteriorates, where some non-physical small-scale spanwise streaks are observed, as in figure 6( f,g) in the RS and TU cases. In the meantime, although the flow structures estimated by the GRBE also become blurred a little in cases RS and TU, which is attributed to the decrease of measurement information, they still provide credible estimations of the large-scale structures.

Figure 6. Instantaneous flow state at $y^+ = 100$ from the DNS results ($a$), the GRBE results (b–d) and the RBE results (e–g). The sampling time intervals are set as temporally resolved (in the TR case) (b,e), randomly sampled (in the RS case) (c,f) and temporally unresolved (in the TU case) (d,g). The values shown in the figures are normalised by the maximum velocity fluctuation value in the DNS result.

To quantitatively investigate the estimation accuracy of the GRBE and RBE results, the relative MSE for the estimated turbulent channel flow is defined as

(3.3)\begin{equation} {\rm Err}(y)=\frac{\left\langle \sum_{k_x}\sum_{k_z} \sum_t \hat{\varepsilon}(k_x,k_z,y,t) \bar{\hat{\varepsilon}}(k_x,k_z,y,t)\right\rangle }{ \left\langle \sum_{k_x}\sum_{k_z} \sum_t \hat{u}(k_x,k_z,y,t)\bar{\hat{u}}(k_x,k_z,y,t)\right\rangle }. \end{equation}

The profiles of the relative MSEs from the GRBE and RBE with different $Re_{\tau }$ and sampling time intervals are shown in figure 7. In the TR case where $\Delta t = \delta t$, the curves denoting the estimation error from the GRBE and RBE overlap with each other, which again demonstrates that the GRBE is equivalent to RBE when the sampling time is resolved in time, as pointed out in § 2.4. When the sampling intervals are randomly distributed in the RS case, the estimation errors of the GRBE along the height slightly increase compared with those in the TR case due to the reduction of measurement information. Meanwhile, the estimation errors of the RBE in the RS case are much larger in comparison with those in the TR case. When the sampling time interval increases to $32 \delta t^+$ in the TU case, the estimation errors of the RBE results become even larger, while those from the GRBE maintain a relatively low magnitude. For instance, compared with the TR cases, the relative MSE from the RBE at $y^+ = 100$ with $Re_{\tau }=934$ increases by $43.7\,\%$ in the TU case, while that from the GRBE increases by only $10.1\,\%$.

Figure 7. Comparison of the relative MSEs from the RBE and GRBE with different measurement sampling intervals with $Re_{\tau } = 186$ ($a$), 547 ($b$) and 934 ($c$).

From the above comparisons between the estimation results from the RBE and GRBE, the GRBE is capable of providing a satisfying estimation for the turbulent channel flows with different sampling time intervals when the measurements are wall shear stress and pressure with real flow statistics informing the estimator.

4.1. Impacts of the abundance of measurement information on the estimation accuracy

The impact of different observation types, including wall shear stress and pressure, is preliminarily studied in Amaral et al. (Reference Amaral, Cavalieri, Martini, Jordan and Towne2021), leading to the conclusion that the estimation accuracy is improved with the increasing amount of the measured quantity types. Meanwhile, the problem becomes more complex when the abundance degree of temporal information is introduced. The variations of estimation accuracy with the measurement abundance at different flow scales also await further examination.

To investigate the impact of the abundance of temporal measurement information, three cases with different sampling temporal intervals are set, namely TU for the temporally unresolved case, RS for the randomly sampled case and TR for the temporally resolved case. Note that the sampling intervals in the TU case in this section are further enlarged to be $240\delta t$, $160\delta t$ and $160\delta t$ for turbulent channel flows with $Re_{\tau } = 186$, 547 and 934, respectively, which are all the same as the lengths of the DNS data in temporal direction for each block with corresponding $Re_{\tau }$, as presented in § 3.2.1. This means that the flow information at only one snapshot is used for each TU case. On the other hand, to investigate the impacts of the sensor types on the estimation accuracy, different combinations of wall quantities, including the streamwise shear stress $\tau _u$, spanwise shear stress $\tau _w$ and pressure $p$, are also considered in this study. All of the to-be-tested combinations of the sampling intervals and sensor types in this study are summarised in table 2. In this section, the estimations with wall measurements informed by the real flow statistics are conducted.

Table 2. Nomenclature of the cases testing the impacts of the sampling time intervals and sensor types on estimating the turbulent channel flow in this study.

Instantaneous snapshots of the DNS result and the estimations at $y^+ = 100$ with $Re_{\tau }=934$ are depicted in figure 8. Different types of sensors with various extents of temporal information are considered. In case TR$_{u,w,p}$ where the measured flow information is the most abundant, the GRBE estimator provides the most accurate results, where the flow structures recover well the representative ones in the DNS results. With the decrease of either the sensor types or temporal information, the estimated small flow structures become blurred. In case TU$_u$ with the least measured information, only the large-scale structures with long streaky structures could be observed. Interestingly, it is found that the flow structures are still well estimated in case TU$_{u,w,p}$ where the flow information at only one snapshot is used for the estimation. Compared with case TR$_{u}$ with measurements of $N_x \times N_z \times N_T$, the results from case TU$_{u,w,p}$ with $3 \times N_x \times N_z$ provide comparable flow information with a much lower requirement of measured data.

Figure 8. Estimated instantaneous flow state at $y^+ = 100$ with measurements at the wall from the GRBE. The DNS result ($a$) is presented as a reference. The measured quantities are set as $(\tau _u,\tau _w,p)$ (b–d), $(\tau _u,\tau _w)$ (e–g) and $\tau _u$ (h–j). The sampling time intervals are set as temporally resolved (b,e,h), randomly sampled (c,f,i) and temporally unresolved (d,g,j). The values shown in the figures are normalised by the maximum velocity fluctuation value in the DNS result.

Next, we further investigate the impacts of the measurement abundance through the profiles of the relative MSEs. Figure 9 shows the profiles of the relative MSEs along the height. For all the tested $Re_{\tau }$, it is observed that the estimation errors generally increase with the enlarged distance from the measurement location, i.e. the wall. At each height, the results in TR cases are the most accurate among those with the same sensor types. Meanwhile, case TU$_{u,w,p}$ that only needs the measured data at one instant has an even smaller error compared with case TR$_{u}$ that needs temporally resolved measurements when $y^+ < 100$. The Reynolds number effects on the GRBE accuracy can also be observed from figure 9. At a given height $y^+$, the estimation errors decrease with the enlargement of $Re_{\tau }$ for each case. For instance, at $y^+ = 100$, the estimation errors in case TR$_{u,w,p}$ are equal to 0.88, 0.60 and 0.48 when $Re_{\tau } = 186$, 547 and 934, respectively. In the meantime, an approximate plateau of the relative MSEs in the buffer zone is observed when $Re_{\tau } = 547$ and 934, which is more obvious with the higher $Re_{\tau }$. When $Re_{\tau } = 934$, the estimation errors even slightly decrease with the increase of height for cases with measurements of $\tau _u$.

Figure 9. Profiles of the relative MSEs of the streamwise velocity fluctuations from the GRBE in turbulent channel flows with $Re_{\tau }=186$ ($a$), 547 ($b$) and 934 ($c$).

To further excavate the underlying mechanism of the differences in the estimation accuracy with different $Re_{\tau }$, the relative estimation error energy as a function of $y$ and the flow scale $\lambda _x$ is defined by

(4.1)\begin{equation} {\rm Err}_x(\lambda_x , y) =\frac{\left\langle \sum_{k_z} \sum_t \hat{\varepsilon}(k_x,k_z,y,t) \bar{\hat{\varepsilon}}(k_x,k_z,y,t)\right\rangle }{ \left\langle \sum_{k_z} \sum_t \hat{u}(k_x,k_z,y,t)\bar{\hat{u}}(k_x,k_z,y,t)\right\rangle }, \end{equation}

where $\lambda _x = 2 {\rm \pi}/ k_x$. As shown in figure 10, with the increasing flow scale, the estimation errors at a given height decrease monotonically. In the meantime, the plateaus in the estimation errors that are observed in the overall error energy profiles in figure 9(b,c) are much more obvious for the VLSMs ($\lambda _x \geqslant 3h$) (Adrian, Meinhart & Tomkins Reference Adrian, Meinhart and Tomkins2000; Lee & Sung Reference Lee and Sung2013; Lee et al. Reference Lee, Lee, Choi and Sung2014) in the buffer zone. Especially, for the VLSMs with $\lambda _x = 4 {\rm \pi}$, the values of the estimation errors are around 0.18 at the plateaus at $y^+ = 40 \sim 100$ in case TR$_{u,w,p}$ with $Re_{\tau } = 934$. Meanwhile, for the motions with small scales of $\lambda _x \leqslant 3h$, the plateaus are not found in the buffer zone, which indicates that the nonlinear effects are stronger for such small-scale motions.

Figure 10. Relative MSEs of the streamwise velocity fluctuations from the GRBE in turbulent channel flows as functions of the flow scale $\lambda _x^+$ and wall-normal distance $y^+$ in cases ${\rm TR}_{u,v,w}$ (a,c,e) and ${\rm TU}_u$ (b,d,f) with $Re_{\tau} = 934$ (a,b), 547 (c,d), and 186 (e,f). The blue curves denote the profiles of the relative MSEs for flow motions with given $\lambda^{+}_{x}$, while the red ones denote the contours of ${\rm Err}_x = 0.2$, 0.4, 0.6 and 0.8, respectively.

In the following, deeper investigations on the GRBE estimator are conducted by studying the two-dimensional energy spectra of the estimation and the corresponding energy distributions of the errors as functions of $\lambda _x$ and $\lambda _z$. Similar to (4.1), the two-dimensional relative energy distribution of the estimation error at a given height $y_o$ is defined by

(4.2)\begin{equation} {\rm Err}_{xz}(\lambda_x , \lambda_z) =\frac{\left\langle \sum_t \hat{\varepsilon}(k_x,k_z,y_o,t) \bar{\hat{\varepsilon}}(k_x,k_z,y_o,t)\right\rangle }{ \left\langle \sum_t \hat{u}(k_x,k_z,y_o,t)\bar{\hat{u}}(k_x,k_z,y_o,t)\right\rangle }, \end{equation}

where $\lambda _x = 2 {\rm \pi}/ k_x$. The premultiplied spectra of estimated velocity fluctuation energy and the corresponding relative MSEs of $u$ at $y^+ = 40$ and 100 are shown in figures 11 and 12, respectively. To quantify the spectral region with relatively high estimation accuracy, the contours of $\textrm {Err}_{xz}(\lambda _x , \lambda _z) = 0.25$ and 0.75 in the corresponding figures are highlighted with white solid and dashed curves, respectively. Meanwhile, the scope of the self-similar attached eddies, i.e. $\lambda _x \geqslant 2\lambda _z$, $\lambda _x \geqslant 10y$, and $\lambda _z \geqslant 2y$, as suggested by Lozano-Durán, Flores & Jiménez (Reference Lozano-Durán, Flores and Jiménez2012), Hwang (Reference Hwang2015) and Cheng et al. (Reference Cheng, Li, Lozano-Durán and Liu2019), are also highlighted with black dot-dashed lines.

Figure 11. Premultiplied energy spectra ($1\textrm {st}$ and $3\textrm {rd}$ panels) and the relative error energies ($2\textrm {nd}$ and $4\textrm {th}$ panels) of $u$ at $y^+ = 40$ using the measurements at the wall. The sensor types are $\tau _u$, $\tau _w$, and $p$ (a–d), $\tau _u$ (e–h), $\tau _w$ (i–l) and $p$ (m–p). The black dot-dashed lines in each figure denote the range of the attached eddies. In the $1\textrm {st}$ and $3\textrm {rd}$ panels, the solid and dashed curves denote the contours of the energy spectra from DNS and estimations, respectively. In the $2\textrm {nd}$ and $4\textrm {th}$ panels, the contours of $\textrm {Err}_{xz} = 0.25$ and 0.75 are highlighted with white solid and dashed curves, respectively.

Figure 12. Same as figure 11, but at $y^+ = 100$.

From figures 11 and 12, estimated premultiplied energy spectra in the dashed curves match relatively well with the DNS results within the spectral region corresponding to the attached eddies when the sensor types are $\lbrace \tau _u,\tau _w,p \rbrace$ or $\lbrace \tau _u \rbrace$. Correspondingly, the relative energies of the estimation errors are generally smaller for the flow portions with larger wavelengths. On the other hand, the spectral regions where $\textrm {Err}_{xz} \leqslant 0.75$ at $y^+ = 40$ are much larger than those at $y^+ = 100$, which means that the estimation errors of the small-scale structures are significantly enlarged. Meanwhile, the extents of the regions where $\textrm {Err}_{xz} \leqslant 0.25$ for the large-scale flows do not shrink so obviously when the estimation layers lifts from $y^+ = 40$ to $y^+ = 100$. This phenomenon is consistent with the plateaus for the VLSMs existing at the buffer zones observed in figure 10.

The sensor types are found to have significant impacts on the energy distributions of the estimation errors in figures 11 and 12. Comparing the results in cases TR$_{u,w,p}$ and TR$_{u}$, obvious increases of the estimation errors are observed in the spectral region of $10y < \lambda _x < 2\lambda _z$. At $(\lambda _x / h , \lambda _z / h) = (0.36{\rm \pi}, {\rm \pi})$ on $y^+ = 40$, the value of the error energy is 0.52 and 0.94 in cases TR$_{u,w,p}$ and TR$_{u}$, respectively, with an increasing rate of $81.5\,\%$. On the other hand, at $(\lambda _x / h,\lambda _z / h) = (4{\rm \pi}, 0.5{\rm \pi} )$ inside the scope of the attached eddies, the error energy increases by $20.0\,\%$ in case TR$_u$ in comparison with that in case TR$_{u,w,p}$, which is not so dramatic as that at $(\lambda _x / h,\lambda _z / h) = (0.36{\rm \pi}, {\rm \pi})$. Besides the estimation results with measurements of $\lbrace \tau _u , \tau _w , p \rbrace$ and $\lbrace \tau _u \rbrace$, those with $\lbrace \tau _w \rbrace$ and $\lbrace p \rbrace$ are also shown in figures 11 and 12 to study the characteristics of the GRBE in cross-quantity estimation. In case TR$_{w}$, the region where $\textrm {Err}_{xz}(\lambda _x,\lambda _z) \leqslant 0.75$ that is surrounded by the white dashed lines is comparable to that in case TR$_{u}$, which roughly matches the range of the attached eddies. However, for large-scale flow motions, the estimation errors in case TR$_{w}$ are generally larger than 0.25, which is larger than those in TR$_{u}$. Thus, although the GRBE is capable of providing fair estimations of the small-scale structures of $u$ with measurements of $\tau _w$, the large-scale structures cannot be well estimated from cross-quantity estimations. In case TR$_{p}$ with wall pressure measurements, the estimation errors are obviously lower than those in the TR$_{u}$ case for the small-scale flows, which is beyond the scope of the attached eddies. For $(\lambda _x / h,\lambda _z / h) = (0.36{\rm \pi}, {\rm \pi})$ at $y^+ = 40$, the error energy in case TR$_{p}$ is $0.75$, which is $20\,\%$ lower than that in case TR$_{u}$. Meanwhile, for the estimations of large-scale structures, the GRBE in the TR$_{p}$ case is not as good as those in other cases.

Besides the sensor types, the abundance of the measurement temporal information is indeed another factor that influences the performance of the GRBE. Compared with cases TR$_{u,w,p}$ and TU$_{u,w,p}$, the reduction of the measured information in time reduces the estimation accuracy for all the considered spatial scales. Unlike the sensor types, which have significant influences on the patterns of the energy distributions of errors, the variations of the temporal information of measurements do not obviously change the error energy distributions. For instance, at $(\lambda _x,\lambda _z) = (0.36{\rm \pi}, {\rm \pi})$ and $(4{\rm \pi}, 0.5{\rm \pi} )$ on $y^+ = 40$ that are outside and inside the scopes of attached eddies, the value of the relative error energy increases by $46.2\,\%$ and $53.9\,\%$ in cases TU$_{u,w,p}$ compared with those in the TR$_{u,w}$ case, respectively, which are relatively close to each other.

To find out the performance of the GRBE with the variations of sampling intervals, we gradually increase the sampling time interval. The variations of the relative MSEs with the enlargement of the sampling intervals are shown in figure 13(a), while the increasing rate of the relative MSEs compared with those in the TR cases are shown in figure 13(b). For each combination of the sensor types, the estimation error does not obviously increase when $\Delta t / \delta t \leqslant 5$. The relative deviations compared with the TR cases are no greater than $5\,\%$ until $\Delta t / \delta t$ increases to $16$ ($\Delta t^+ = 81.41$). When $16 \leqslant \Delta t / \delta t \leqslant 80$, a relatively obvious increase in the estimation error with the enlarging sampling interval is observed. When the sampling intervals increase over $80 \delta t$, the estimation errors become almost unchanged with a further increase of $\Delta t$. For the estimations with measurements including $\tau _u$, the largest relative MSEs in the case with $\Delta t = 160 \delta t$ are around $15\,\%$ larger than the temporally resolved cases. Meanwhile, the errors of the cross-quantity estimation with measurements of $\lbrace \tau _w \rbrace$ and $\lbrace p\rbrace$ increase by around $8\,\%$ in TU cases compared with the TR ones. The relatively small increase in the cross-quantity estimation errors should be attributed to the fact that their magnitudes are already large in the TR cases, as shown in figure 13(a). In the meantime, it is noticed that the estimation error in case TU$_{u,w,p}$ is almost the same as that in the TR$_{u}$ case with a discrepancy of less than $1\,\%$.

Figure 13. Variations of the estimation error energies with measurement time intervals. ($a$) Values of the relative MSEs. ($b$) Variations of the relative MSEs compared with those in the TR cases with the corresponding measured quantities.

To further investigate the impact of sampling temporal intervals on the estimation accuracy, the increasing rates of the error energies compared with those in the TR cases as functions of measurement time intervals $\Delta t$ and the streamwise wavelength $\lambda _x$ are depicted in figure 14. To quantify the impact of the increasing sampling temporal intervals on the estimation accuracy of flow motions with different scales, the threshold sampling interval $\Delta t_{T}$ is defined herein as the maximum sampling interval with which the increasing rate of $\textrm {Err}_x$ compared with the TR cases is smaller than $5\,\%$. For the attached eddies with $\lambda _x \geqslant 10y$, the threshold temporal interval $\Delta t_{T}$ is generally larger than $20 \delta t$, where $\Delta t_{T}^+ \geqslant 101.76$. On the other hand, when $3y \leqslant \lambda _x < 10y$, $\Delta t_{T}$ decreases with decreasing wavelength, which indicates that finer temporal resolution is required to provide effective information of the smaller flow scales. When $\lambda _x = 3y$ with the sensor of $\lbrace \tau _u , \tau _w , p \rbrace$, $\Delta t_{T}$ decreases to the minimum value of approximately $10 \delta t$. When the flow scales further decrease to $\lambda _x < 3y$, $\Delta t_{T}$ dramatically increases instead, which is attributed to the fact that the estimation errors with $\lambda _x < 3y$ are already near 1.0 in the TR cases, as shown in figure 10.

Figure 14. The increasing rate of the estimation error energies compared with those in the TR case as functions of sampling time intervals $\Delta t$ and streamwise wavelengths $\lambda _x$ at $y^+ = 100$ with wall measurements of $\lbrace \tau _u , \tau _w , p \rbrace$ ($a$), $\lbrace \tau _u \rbrace$ ($b$), $\lbrace \tau _w \rbrace$ ($c$) and $\lbrace p \rbrace$ ($d$). The black dashed and solid lines denote the contours of $5\,\%$ and $15\,\%$, respectively.

From the above discussions, it can be concluded that the temporal information has limited impacts on the estimation accuracy, with a maximum reduction of around $13.0\,\%$ of the estimation error in the TR case compared with the TU cases at $y^+ = 100$ with measurements of wall shear stress and pressure. When the sampling time intervals $\Delta t = 16 \delta t$ ($\Delta t^+ = 81.41$), the increasing rates of the estimation errors are around $5\,\%$ compared with the TR results. Such findings indicate that the measurements with sparse samplings can be effectively used to provide reasonable estimations of the flow state with the GRBE. On the other hand, the inclusion of more types of observations can be efficient in improving the estimation accuracy. In the TU cases without temporal measurement information, the estimation error reduces by $12.3\,\%$ when the sensors for $\tau _w$ and $p$ are added besides that for $\tau _u$. Consequently, the relative MSEs at $y^+ = 100$ in the TU$_{u,w,p}$ case with the most sensor types but least temporal information are almost the same as those in the TR$_u$ case with the least sensor type but most temporal information.

Besides the above conclusions, we would also like to discuss the tolerated estimation error from the GRBE when applied in real engineering scenarios. Generally speaking, the feasibility of the GRBE depends on specific requirements on the estimation accuracy for a given range of flow scales. For instance, in flow control problems, the combination of sensors and actuators that are both limited to one wall-normal location provides a comparable control performance to the case where the limited-in-space actuators are informed by the entire flow information (Oehler & Illingworth Reference Oehler and Illingworth2021). This is due to the fact that the linear actuators are also limited by the flow coherence like the linear estimator, which means that the energy-containing LSMs could be better controlled. Since the GRBE provides reliable estimations of the large-scale structures in the inner layer, it can already well excavate the potentials of the linear actuator in controlling the large-scale turbulent flows. On the other hand, for the estimation problem of small-scale structures outside the scale range of the attached eddies, which is unresolved yet for all the current linear and nonlinear estimators (e.g. Guastoni et al. Reference Guastoni, Güemes, Ianiro, Discetti, Schlatter, Azizpour and Vinuesa2021) to the best of the authors’ knowledge, the GRBE is also not able to provide reliable results. Meanwhile, since GRBE has been validated to be the optimal linear estimator, the comprehensive discussions in this section are meaningful to evaluate the reliability and quantify the uncertainties of the linear estimation results in engineering applications.

4.2. Performance of the forcing models in flow state estimation

In the above, the performance of the GRBE estimators informed by the real flow statistics under different cases is discussed. Thus, no matter how the performance varies in different cases, they are already the optimal linear estimators with given measurements. In practice, the flow statistics are not always available when estimating the flow state, where the forcing models that predict the spatio-temporal statistical properties of turbulent flows are needed to inform the GRBE estimator. Thus, it is important to clarify their performance when informing the linear estimator of GRBE. In the following, three existing forcing models that provide spatio-temporal descriptions of the forcing statistics, including the B-model (Hwang & Cossu Reference Hwang and Cossu2010; Madhusudanan et al. Reference Madhusudanan, Illingworth and Marusic2019), W-model and $\lambda$-model (Gupta et al. Reference Gupta, Madhusudanan, Wan, Illingworth and Juniper2021), will be applied to inform the GRBE estimator. Mathematical descriptions of these three models are provided in Appendix B. Estimation results with the wall measurements will be used for the discussions in this section.

Figure 15 shows the estimation errors of different forcing models with different sensor types along the wall-normal direction. Moreover, to provide intuitive results, figures 16, 17 and 18 show the instantaneous flow states estimated by GRBE estimators that are informed by the B-model, W-model and $\lambda$-model, respectively. In the TU case, the W-model-informed GRBE provides the most accurate results, which are closer to the optimal linear estimated results than those from the other models. In case TU$_u$, the estimation error in the W-model-informed result at $y^+ = 100$ is $13.3\,\%$ larger than that from the optimal linear estimations, as in figure 15(b). From figure 17(d,g,j), it is also observed that the estimated large-scale structures from the W-model match well with the DNS results in the TU cases. Such good agreement between the W-model-informed result and the optimal linear result in case TU$_u$ indicates that the W-model can effectively predict the spatial coherence of the streamwise velocity between the wall and the higher region. Compared with the estimation errors from the W-model and $\lambda$-model, those from the B-model are obviously larger in the near-wall region of $y^+ \leqslant 10$. Meanwhile, with the increase of estimation height, the errors from the $\lambda$-model rapidly increase and exceed those from the B-model where $y^+ \geqslant 15$. On the other hand, the estimation error from each forcing model does not obviously decrease with the increase of sensor types. When converting the sensor types from $\lbrace \tau _u \rbrace$ to $\lbrace \tau _u, \tau _w, p \rbrace$, the estimation errors at $y^+ = 100$ decrease by $0.7\,\%$, $3.9\,\%$ and $8.0\,\%$ for the B-model, W-model and $\lambda$-model, respectively, which are lower than the value of $13.3\,\%$ by the optimal linear estimator. The relatively low efficiency of the forcing-model-informed GRBE estimator to improve the estimation accuracy with an increasing number of measured variables is attributed to the fact that all the tested forcing models cannot effectively predict the covariance between different quantities.

Figure 15. Variations of the estimation errors from the B-model (a,d,g), W-model (b,e,h) and $\lambda$-model (c,f,i) with the wall-normal height in the TU (a–c), RS (d–f) and TR (g–i) cases. The scattered symbols and translucent curves denote the estimation results from the GRBE estimator incorporating the forcing models and real flow statistics, respectively.

Figure 16. Estimations of the flow state with the B-model-informed estimator. The depicted figures are reference DNS result ($a$); estimations with measurements of wall shear stress and pressure (b–d), wall shear stress (e–g) and streamwise wall shear stress (h–j). The measurements are temporally resolved (b,e,h), randomly sampled (c,f,i) and temporally unresolved (d,g,j). The values shown in the figures are normalised by the maximum velocity fluctuation value in the DNS result.

Figure 17. Same as figure 16, but with the W-model-informed estimator.

Figure 18. Same as figure 16, but with the $\lambda$-model-informed estimator.

With the increase of the temporal information of measurements, it is interesting to find that the estimation errors from the W-model- and $\lambda$-model-informed estimators increase accordingly, the tendencies of which are opposite to those of the optimal linear estimator. In the TR$_{ u, w, p}$ case, where the estimation error is minimised in the optimal linear estimator, it is instead maximised for the W-model results, which are $48.9\,\%$ larger than that in case TU$_{u, w, p}$. A similar non-physical phenomenon is observed in the $\lambda$-model results. From figures 17(a) and 18(a), the fluctuation energy of small-scale structures is incorrectly amplified by the W- and $\lambda$-models in the TR cases compared with the estimated results in the TU cases. Such non-physical increase of estimation error with the increase of temporal measurement information and sensor types indicates that the W-model and $\lambda$-model cannot effectively predict the temporal properties, although the W-model works well in predicting the spatial coherence of the streamwise velocity. On the other hand, the estimation errors of the B-model results decrease monotonically with the increase of the measurement abundance. Compared with the relative MSE in case TU$_{u,w,p}$ at $y^+ = 100$, that in case TR$_{u,w,p}$ decreases by $13.2\,\%$ from the B-model results.

In summary, the tested forcing models perform differently in cases with different measurement abundance. In case TU$_u$ with the least measurement information, the W-model performs the best. When the sensor types increase from $\lbrace \tau _u \rbrace$ to $\lbrace \tau _u, \tau _w, p \rbrace$, the estimation accuracy of the tested forcing models is not obviously improved. Moreover, with the increase of temporal information, the estimation errors from the W-model and $\lambda$-model non-physically increase. Such tendencies are opposite to those in the optimal linear estimation results, indicating that the temporal properties predicted by the W-model and $\lambda$-model cannot effectively reflect these of the real turbulence. A forcing model that properly models the forcing in the temporal domain will effectively improve the linear estimation accuracy of the GRBE without the requirement of the real flow statistics, which is to be explored in the future.