No Arabic abstract
Common modal decomposition techniques for flowfield analysis, data-driven modeling and flow control, such as proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) are usually performed in an Eulerian (fixed) frame of reference with snapshots from measurements or evolution equations. The Eulerian description poses some difficulties, however, when the domain or the mesh deforms with time as, for example, in fluid-structure interactions. For such cases, we first formulate a Lagrangian modal analysis (LMA) ansatz by a posteriori transforming the Eulerian flow fields into Lagrangian flow maps through an orientation and measure-preserving domain diffeomorphism. The development is then verified for Lagrangian variants of POD and DMD using direct numerical simulations (DNS) of two canonical flow configurations at Mach 0.5, the lid-driven cavity and flow past a cylinder, representing internal and external flows, respectively, at pre- and post-bifurcation Reynolds numbers. The LMA is demonstrated for several situations encompassing unsteady flow without and with boundary and mesh deformation as well as non-uniform base flows that are steady in Eulerian but not in Lagrangian frames. We show that LMA application to steady nonuniform base flow yields insights into flow stability and post-bifurcation dynamics. LMA naturally leads to Lagrangian coherent flow structures and connections with finite-time Lyapunov exponents (FTLE). We examine the mathematical link between FTLE and LMA by considering a double-gyre flow pattern. Dynamically important flow features in the Lagrangian sense are recovered by performing LMA with forward and backward (adjoint) time procedures.
The hydrodynamics of a liquid-vapour interface in contact with an heterogeneous surface is largely impacted by the presence of defects at the smaller scales. Such defects introduce morphological disturbances on the contact line and ultimately determine the force exerted on the wedge of liquid in contact with the surface. From the mathematical point of view, defects introduce perturbation modes, whose space-time evolution is governed by the interfacial hydrodynamic equations of the contact line. In this paper we derive the response function of the contact line to such generic perturbations. The contact line response may be used to design simplified 1+1 dimensional models accounting for the complexity of interfacial flows coupled to nanoscale defects, yet offering a more tractable mathematical framework to include thermal fluctuations and explore thermally activated contact line motion through a disordered energy landscape.
We apply supervised machine learning techniques to a number of regression problems in fluid dynamics. Four machine learning architectures are examined in terms of their characteristics, accuracy, computational cost, and robustness for canonical flow problems. We consider the estimation of force coefficients and wakes from a limited number of sensors on the surface for flows over a cylinder and NACA0012 airfoil with a Gurney flap. The influence of the temporal density of the training data is also examined. Furthermore, we consider the use of convolutional neural network in the context of super-resolution analysis of two-dimensional cylinder wake, two-dimensional decaying isotropic turbulence, and three-dimensional turbulent channel flow. In the concluding remarks, we summarize on findings from a range of regression type problems considered herein.
Modal stability analysis provides information about the long-time growth or decay of small-amplitude perturbations around a steady-state solution of a dynamical system. In fluid flows, exponentially growing perturbations can initiate departure from laminar flow and trigger transition to turbulence. Although flow of a Newtonian fluid through a pipe is linearly stable for very large values of the Reynolds number ($Re sim 10^7$), a transition to turbulence often occurs for $Re$ as low as $1500$. When a dilute polymer solution is used in the place of a Newtonian fluid, the transitional value of the Reynolds number decreases even further. Using the spectral collocation method and Oldroyd-B constitutive equation, Garg et al. (Phys. Rev. Lett. 121:024502, 2018) claimed that such a transition in viscoelastic fluids is related to linear instability. Since differential matrices in the collocation method become ill-conditioned when a large number of basis functions is used, we revisit this problem using the well-conditioned spectral integration method. We show modal stability of viscoelastic pipe flow for a broad range of fluid elasticities and polymer concentrations, including cases considered by Garg et al. Similarly, we find that plane Poiseuille flow is linearly stable for cases where Garg et al. report instability. In both channel and pipe flows, we establish the existence of spurious modes that diverge slowly with finer discretization and demonstrate that these can be mistaken for grid-independent modes if the discretization is not fine enough.
The characterization of intermittency in turbulence has its roots in the K62 theory, and if no proper definition is to be found in the literature, statistical properties of intermittency were studied and models were developed in attempt to reproduce it. The first contribution of this work is to propose a requirement list to be satisfied by models designed within the Lagrangian framework. Multifractal stochastic processes are a natural choice to retrieve multifractal properties of the dissipation. Among them, following the proposition of cite{Mandelbrot1968}, we investigate the Gaussian Multiplicative Chaos formalism, which requires the construction of a log-correlated stochastic process $X_t$. The fractional Gaussian noise of Hurst parameter $H = 0$ is of great interest because it leads to a log-correlation for the logarithm of the process.Inspired by the approximation of fractional Brownian motion by an infinite weighted sum of correlated Ornstein-Uhlenbeck processes, our second contribution is to propose a new stochastic model: $X_t = int_0^infty Y_t^x k(x) d x$, where $Y_t^x$ is an Ornstein-Uhlenbeck process with speed of mean reversion $x$ and $k$ is a kernel. A regularization of $k(x)$ is required to ensure stationarity, finite variance and logarithmic auto-correlation. A variety of regularizations are conceivable, and we show that they lead to the aforementioned multifractal models.To simulate the process, we eventually design a new approach relying on a limited number of modes for approximating the integral through a quadrature $X_t^N = sum_{i=1}^N omega_i Y_t^{x_i}$, using a conventional quadrature method. This method can retrieve the expected behavior with only one mode per decade, making this strategy versatile and computationally attractive for simulating such processes, while remaining within the proposed framework for a proper description of intermittency.
In recent years, there have been a surge in applications of neural networks (NNs) in physical sciences. Although various algorithmic advances have been proposed, there are, thus far, limited number of studies that assess the interpretability of neural networks. This has contributed to the hasty characterization of most NN methods as black boxes and hindering wider acceptance of more powerful machine learning algorithms for physics. In an effort to address such issues in fluid flow modeling, we use a probabilistic neural network (PNN) that provide confidence intervals for its predictions in a computationally effective manner. The model is first assessed considering the estimation of proper orthogonal decomposition (POD) coefficients from local sensor measurements of solution of the shallow water equation. We find that the present model outperforms a well-known linear method with regard to estimation. This model is then applied to the estimation of the temporal evolution of POD coefficients with considering the wake of a NACA0012 airfoil with a Gurney flap and the NOAA sea surface temperature. The present model can accurately estimate the POD coefficients over time in addition to providing confidence intervals thereby quantifying the uncertainty in the output given a particular training data set.