  Understanding, predicting and controlling laminar-turbulent boundary-layer
transition is crucial for the next generation aircraft design. However, in real
flight experiments, or wind tunnel tests, often only sparse sensor measurements
can be collected at fixed locations. Thus, in developing reduced models for
predicting and controlling the flow at the sensor locations, the main challenge
is in accounting for how the surrounding field of unobserved variables
interacts with the observed variables at the fixed sensor locations. This makes
the Mori-Zwanzig (MZ) formalism a natural choice, as it results in the
Generalized Langevin Equations which provides a framework for constructing
non-Markovian reduced-order models that includes the effects the unresolved
variables have on the resolved variables. These effects are captured in the so
called memory kernel and orthogonal dynamics. In this work, we explore the
data-driven MZ formulations to two boundary layer flows obtained from DNS data;
a low speed incompressible flow; and a high speed compressible flow over a
flared cone at Mach 6. An array of "sensors" are placed near the surface of the
solid boundary, and the MZ operators are learned and the predictions are
compared to the Extended Dynamic Mode Decomposition (EDMD), both using delay
embedded coordinates. Further comparisons are made with Long Short-Term Memory
(LSTM) and a regression based projection framework using neural networks for
the MZ operators. First we compare the effects of including delay embedded
coordinates with EDMD and Mori based MZ and provide evidence that using both
memory and delay embedded coordinates minimizes generalization errors on the
relevant time scales. Next, we provide numerical evidence that the data-driven
regression based projection MZ model performs best with respect to the
prediction accuracy (minimum generalization error) on the relevant time scales.
