Showing papers in "arXiv: Fluid Dynamics in 2018"

Centrifugal Microfluidics for Biomedical Applications

Abstract: Yoon-Kyoung Cho is a full professor in Biomedical Engineering at UNIST. In this contribution she describes a fully integrated and automated lab-on-a-disc using centrifugal microfluidics to provide a "sample-in and answer-out" type of biochemical analysis solution with simple, size-reduced, and cost-efficient instrumentation. She also gives various examples of the fully integrated "lab-on-a-disc" developed for broad applications ranging from medical diagnostics to food, environment, and energy applications.

Turbulence Modeling in the Age of Data

Abstract: Data from experiments and direct simulations of turbulence have historically been used to calibrate simple engineering models such as those based on the Reynolds-averaged Navier--Stokes (RANS) equations. In the past few years, with the availability of large and diverse datasets, researchers have begun to explore methods to systematically inform turbulence models with data, with the goal of quantifying and reducing model uncertainties. This review surveys recent developments in bounding uncertainties in RANS models via physical constraints, in adopting statistical inference to characterize model coefficients and estimate discrepancy, and in using machine learning to improve turbulence models. Key principles, achievements and challenges are discussed. A central perspective advocated in this review is that by exploiting foundational knowledge in turbulence modeling and physical constraints, data-driven approaches can yield useful predictive models.

Rortex A New Vortex Vector Definition and Vorticity Tensor and Vector Decompositions

Abstract: A vortex is intuitively recognized as the rotational/swirling motion of the fluids. However, an unambiguous and universally-accepted definition for vortex is yet to be achieved in the field of fluid mechanics, which is probably one of the major obstacles causing considerable confusions and misunderstandings in turbulence research. In our previous work, a new vector quantity which is called vortex vector was proposed to accurately describe the local fluid rotation and clearly display vortical structures. In this paper, the definition of the vortex vector, named Rortex here, is revisited from the mathematical perspective. The existence of the rotational axis is proved through real Schur decomposition. Based on real Schur decomposition, a fast algorithm for calculating Rortex is also presented. In addition, new vorticity tensor and vector decompositions are introduced: the vorticity tensor is decomposed to a rigidly rotational part and an anti-symmetric deformation part, and the vorticity vector is decomposed to a rigidly rotational vector and a non-rotational vector. Several cases, including 2D Couette flow, 2D rigid rotational flow and 3D boundary layer transition on a flat plate, are studied to demonstrate the justification of the definition of Rortex. It can be observed that Rortex identifies both the precise swirling strength and the rotational axis, and thus it can reasonably represent the local fluid rotation and provide a new powerful tool for vortex dynamics and turbulence research.

Super-resolution reconstruction of turbulent flows with machine learning

Abstract: We use machine learning to perform super-resolution analysis of grossly under-resolved turbulent flow field data to reconstruct the high-resolution flow field. Two machine-learning models are developed; namely the convolutional neural network (CNN) and the hybrid Downsampled Skip-Connection Multi-Scale (DSC/MS) models. These machine-learning models are applied to two-dimensional cylinder wake as a preliminary test and show remarkable ability to reconstruct laminar flow from low-resolution flow field data. We further assess the performance of these models for two-dimensional homogeneous turbulence. The CNN and DSC/MS models are found to reconstruct turbulent flows from extremely coarse flow field images with remarkable accuracy. For the turbulent flow problem, the machine-leaning based super-resolution analysis can greatly enhance the spatial resolution with as little as 50 training snapshot data, holding great potential to reveal subgrid-scale physics of complex turbulent flows. With the growing availability of flow field data from high-fidelity simulations and experiments, the present approach motivates the development of effective super-resolution models for a variety of fluid flows.

Spectral analysis of the budget equation in turbulent channel flows at high Re

Abstract: The transport equations for velocity variances are investigated using data from DNS of incompressible channel flows at $Re_\tau$ up to 5200. Each term in the transport equation has been spectrally decomposed to expose the contribution of turbulence at different length scales to the processes governing the flow of energy in the wall-normal direction, in scale and among components. The outer-layer turbulence is dominated by very large-scale streamwise elongated modes. Away from the wall, production occurs primarily in these large-scale streamwise-elongated modes in the streamwise velocity, but dissipation occurs nearly isotropically in both velocity components and scale. For this to happen, the energy is transferred from the streamwise elongated modes to modes with a range of orientations through non-linear interactions, and then transferred to other velocity components. This allows energy to be transferred more-or-less isotropically from these large scales to the small scales at which dissipation occurs. The VLSMs also transfer energy to the wall-region resulting in a modulation of the autonomous near-wall dynamics. The near-wall energy flows are consistent with the well-known autonomous near-wall dynamics. Through the overlap region between outer and inner layer turbulence, there is a self-similar structure to the energy flows. The VLSM production occurs at spanwise scales that grow with $y$. There is transport of energy away from the wall over a range of scales that grows with $y$. And, there is transfer of energy to small dissipative scales which grow like $y^{1/4}$. Finally, the small-scale near-wall processes characterised by wavelengths less than 1000 wall units are largely Reynolds number independent, while the larger-scale outer layer process are strongly Reynolds number dependent. The interaction between them appears to be relatively simple.

Artificial Neural Networks trained through Deep Reinforcement Learning discover control strategies for active flow control

Abstract: We present the first application of an Artificial Neural Network trained through a Deep Reinforcement Learning agent to perform active flow control. It is shown that, in a 2D simulation of the Karman vortex street at moderate Reynolds number (Re = 100), our Artificial Neural Network is able to learn an active control strategy from experimenting with the mass flow rates of two jets on the sides of a cylinder. By interacting with the unsteady wake, the Artificial Neural Network successfully stabilizes the vortex alley and reduces drag by about 8%. This is performed while using small mass flow rates for the actuation, on the order of 0.5% of the mass flow rate intersecting the cylinder cross section once a new pseudo-periodic shedding regime is found. This opens the way to a new class of methods for performing active flow control.

Deep Learning of Vortex Induced Vibrations

Abstract: Vortex induced vibrations of bluff bodies occur when the vortex shedding frequency is close to the natural frequency of the structure Of interest is the prediction of the lift and drag forces on the structure given some limited and scattered information on the velocity field This is an inverse problem that is not straightforward to solve using standard computational fluid dynamics (CFD) methods, especially since no information is provided for the pressure An even greater challenge is to infer the lift and drag forces given some dye or smoke visualizations of the flow field Here we employ deep neural networks that are extended to encode the incompressible Navier-Stokes equations coupled with the structure's dynamic motion equation In the first case, given scattered data in space-time on the velocity field and the structure's motion, we use four coupled deep neural networks to infer very accurately the structural parameters, the entire time-dependent pressure field (with no prior training data), and reconstruct the velocity vector field and the structure's dynamic motion In the second case, given scattered data in space-time on a concentration field only, we use five coupled deep neural networks to infer very accurately the vector velocity field and all other quantities of interest as before This new paradigm of inference in fluid mechanics for coupled multi-physics problems enables velocity and pressure quantification from flow snapshots in small subdomains and can be exploited for flow control applications and also for system identification

Turbulent superstructures in Rayleigh-B\'enard convection

Abstract: Turbulent Rayleigh-Benard convection displays a large-scale order in the form of rolls and cells on lengths larger than the layer height once the fluctuations of temperature and velocity are removed. These turbulent superstructures are reminiscent of the patterns close to the onset of convection. They are analyzed by numerical simulations of turbulent convection in fluids at different Prandtl number ranging from 0.005 to 70 and for Rayleigh numbers up to $10^7$. For each case, we identify characteristic scales and times that separate the fast, small-scale turbulent fluctuations from the gradually changing large-scale superstructures. The characteristic scales of the large-scale patterns, which change with Prandtl and Rayleigh number, are also found to be correlated with the boundary layer dynamics, and in particular the clustering of thermal plumes at the top and bottom plates. Our analysis suggests a scale separation and thus the existence of a simplified description of the turbulent superstructures in geo- and astrophysical settings.

Sub-grid modelling for two-dimensional turbulence using neural networks

Abstract: In this investigation, a data-driven turbulence closure framework is introduced and deployed for the sub-grid modelling of Kraichnan turbulence. The novelty of the proposed method lies in the fact that snapshots from high-fidelity numerical data are used to inform artificial neural networks for predicting the turbulence source term through localized grid-resolved information. In particular, our proposed methodology successfully establishes a map between inputs given by stencils of the vorticity and the streamfunction along with information from two well-known eddy-viscosity kernels. Through this we predict the sub-grid vorticity forcing in a temporally and spatially dynamic fashion. Our study is both a-priori and a-posteriori in nature. In the former, we present an extensive hyper-parameter optimization analysis in addition to learning quantification through probability density function based validation of sub-grid predictions. In the latter, we analyse the performance of our framework for flow evolution in a classical decaying two-dimensional turbulence test case in the presence of errors related to temporal and spatial discretization. Statistical assessments in the form of angle-averaged kinetic energy spectra demonstrate the promise of the proposed methodology for sub-grid quantity inference. In addition, it is also observed that some measure of a-posteriori error must be considered during optimal model selection for greater accuracy. The results in this article thus represent a promising development in the formalization of a framework for generation of heuristic-free turbulence closures from data.

Shapes and paths of an air bubble rising in quiescent liquids

Abstract: Shapes and paths of an air bubble rising inside a liquid are investigated experimentally. About three hundred experiments are conducted in order to generate a phase plot in the Galilei and Eotvos numbers plane, which separates distinct regimes in terms of bubble behaviour. A wide range of the Galilei and Eotvos numbers are obtained by using aqueous glycerol solutions of different concentrations as the surrounding fluid, and by varying the bubble size. The dynamics is investigated in terms of shapes, topological changes and trajectories of the bubbles. Direct numerical simulations are conducted to study the bubble dynamics, which show excellent agreement with the experiments. To the best of our knowledge, this is the first time an experimentally obtained phase plot showing the distinct behaviour of an air bubble rising in a quiescent medium is reported for such a large range of Galilei and Eotvos numbers.

Quantification of Model Uncertainty in RANS Simulations: A Review

Abstract: In computational fluid dynamics simulations of industrial flows, models based on the Reynolds-averaged Navier--Stokes (RANS) equations are expected to play an important role in decades to come. However, model uncertainties are still a major obstacle for the predictive capability of RANS simulations. This review examines both the parametric and structural uncertainties in turbulence models. We review recent literature on data-free (uncertainty propagation) and data-driven (statistical inference) approaches for quantifying and reducing model uncertainties in RANS simulations. Moreover, the fundamentals of uncertainty propagation and Bayesian inference are introduced in the context of RANS model uncertainty quantification. Finally, the literature on uncertainties in scale-resolving simulations is briefly reviewed with particular emphasis on large eddy simulations.

Synthetic turbulent inflow generator using machine learning

Abstract: We propose a methodology for generating time-dependent turbulent inflow data with the aid of machine learning (ML), which has a possibility to replace conventional driver simulations or synthetic turbulent inflow generators. As for the ML model, we use an auto-encoder type convolutional neural network (CNN) with a multi-layer perceptron (MLP). For the test case, we study a fully-developed turbulent channel flow at the friction Reynolds number of ${\rm Re}_{\tau} = 180$ for easiness of assessment. The ML models are trained using a time series of instantaneous velocity fields in a single cross-section obtained by direct numerical simulation (DNS) so as to output the cross-sectional velocity field at a specified future time instant. From the a priori test in which the output from the trained ML model are recycled to the input, the spatio-temporal evolution of cross-sectional structure is found to be reasonably well reproduced by the proposed method. The turbulence statistics obtained in the a priori test are also, in general, in reasonable agreement with the DNS data, although some deviation in the flow rate was found. It is also found that the present machine-learned inflow generator is free from the spurious periodicity, unlike the conventional driver DNS in a periodic domain. As an a posteriori test, we perform DNS of inflow-outflow turbulent channel flow with the trained ML model used as a machine-learned turbulent inflow generator (MLTG) at the inlet. It is shown that the present MLTG can maintain the turbulent channel flow for a long time period sufficient to accumulate turbulent statistics, with much lower computational cost than the corresponding driver simulation. It is also demonstrated that we can obtain accurate turbulent statistics by properly correcting the deviation in the flow rate.

Performance of interFoam on the simulation of progressive waves

Abstract: The performance of interFoam (a widely-used solver within the popular open source CFD package OpenFOAM) in simulating the propagation of a nonlinear (stream function solution) regular wave is investigated in this work, with the aim of systematically documenting its accuracy. It is demonstrated that over time there is a tendency for surface elevations to increase, wiggles to appear in the free surface, and crest velocities to become (severely) overestimated. It is shown that increasing the temporal and spatial resolution can mitigate these undesirable effects, but that a relatively small Courant number is required. It is further demonstrated that the choice of discretization schemes and solver settings (often treated as a "black box" by users) can have a major impact on the results. This impact is documented, and it is shown that obtaining a "diffusive balance" is crucial to accurately propagate a surface wave over long distances without requiring exceedingly high temporal and spatial resolutions. Finally, the new code isoAdvector is compared to interFoam, which is demonstrated to produce comparably accurate results, while maintaining a sharper surface. It is hoped that the systematic documentation of the performance of the interFoam solver will enable its more accurate and optimal use, as well as increase awareness of potential shortcomings, by CFD researchers interested in the general CFD simulation of free surface waves.

Turbulent thermal superstructures in Rayleigh-B\'enard convection

Abstract: We report the observation of superstructures, i.e.\ very large-scale and long living coherent structures in highly turbulent Rayleigh-Benard convection up to Rayleigh $Ra=10^9$. We perform direct numerical simulations in horizontally periodic domains with aspect ratios up to $\Gamma=128$. In the considered $Ra$ number regime the thermal superstructures have a horizontal extend of six to seven times the height of the domain and their size is independent of $Ra$. Many laboratory experiments and numerical simulations have focused on small aspect ratio cells in order to achieve the highest possible $Ra$. However, here we show that for very high $Ra$ integral quantities such as the Nusselt number and volume averaged Reynolds number only converge to the large aspect ratio limit around $\Gamma \approx 4$, while horizontally averaged statistics such as standard deviation and kurtosis converge around $\Gamma \approx 8$, and the integral scale converges around $\Gamma \approx 32$, and the peak position of the temperature variance and turbulent kinetic energy spectra only around $\Gamma \approx 64$.

Robust flow field reconstruction from limited measurements via sparse representation

Abstract: In many applications it is important to estimate a fluid flow field from limited and possibly corrupt measurements. Current methods in flow estimation often use least squares regression to reconstruct the flow field, finding the minimum-energy solution that is consistent with the measured data. However, this approach may be prone to overfitting and sensitive to noise. To address these challenges we instead seek a sparse representation of the data in a library of examples. Sparse representation has been widely used for image recognition and reconstruction, and it is well-suited to structured data with limited, corrupt measurements. We explore sparse representation for flow reconstruction on a variety of fluid data sets with a wide range of complexity, including vortex shedding past a cylinder at low Reynolds number, a mixing layer, and two geophysical flows. In addition, we compare several measurement strategies and consider various types of noise and corruption over a range of intensities. We find that sparse representation has considerably improved estimation accuracy and robustness to noise and corruption compared with least squares methods. We also introduce a sparse estimation procedure on local spatial patches for complex multiscale flows that preclude a global sparse representation. Based on these results, sparse representation is a promising framework for extracting useful information from complex flow fields with realistic measurements.

A geometric state function for two-fluid flow in porous media

Abstract: Models that describe two-fluid flow in porous media suffer from a widely-recognized problem that the constitutive relationships used to predict capillary pressure as a function of the fluid saturation are non-unique, thus requiring a hysteretic description. As an alternative to the traditional perspec- tive, we consider a geometrical description of the capillary pressure, which relates the average mean curvature, the fluid saturation, the interfacial area between fluids, and the Euler characteristic. The state equation is formulated using notions from algebraic topology and cast in terms of measures of the macroscale state. Synchrotron-based X-ray micro-computed tomography ({\mu}CT) and high- resolution pore-scale simulation is applied to examine the uniqueness of the proposed relationship for six different porous media. We show that the geometric state function is able to characterize the microscopic fluid configurations that result from a wide range of simulated flow conditions in an averaged sense. The geometric state function can serve as a closure relationship within macroscale models to effectively remove hysteretic behavior attributed to the arrangement of fluids within a porous medium. This provides a critical missing component needed to enable a new generation of higher fidelity models to describe two-fluid flow in porous media.

Data-driven prediction of unsteady flow fields over a circular cylinder using deep learning

Abstract: Unsteady flow fields over a circular cylinder are trained and predicted using four different deep learning networks: convolutional neural networks with and without consideration of conservation laws, generative adversarial networks with and without consideration of conservation laws. Flow fields at future occasions are predicted based on information of flow fields at previous occasions. Deep learning networks are trained first using flow fields at Reynolds numbers of 100, 200, 300, and 400, while flow fields at Reynolds numbers of 500 and 3000 are predicted using the trained deep learning networks. Physical loss functions are proposed to explicitly impose information of conservation of mass and momentum to deep learning networks. An adversarial training is applied to extract features of flow dynamics in an unsupervised manner. Effects of the proposed physical loss functions, adversarial training, and network sizes on the prediction accuracy are analyzed. Predicted flow fields using deep learning networks are in favorable agreement with flow fields computed by numerical simulations.

Sub-grid scale model classification and blending through deep learning

Abstract: In this article we detail the use of machine learning for spatiotemporally dynamic turbulence model classification and hybridization for the large eddy simulations (LES) of turbulence. Our predictive framework is devised around the determination of local conditional probabilities for turbulence models that have varying underlying hypotheses. As a first deployment of this learning, we classify a point on our computational grid as that which requires the functional hypothesis, the structural hypothesis or no modeling at all. This ensures that the appropriate model is specified from \emph{a priori} knowledge and an efficient balance of model characteristics is obtained in a particular flow computation. In addition, we also utilize the conditional probability predictions of the same machine learning to blend turbulence models for another hybrid closure. Our test-case for the demonstration of this concept is given by Kraichnan turbulence which exhibits a strong interplay of enstrophy and energy cascades in the wave number domain. Our results indicate that the proposed methods lead to robust and stable closure and may potentially be used to combine the strengths of various models for complex flow phenomena prediction.

Microfluidic pumping using artificial magnetic cilia

Abstract: One of the vital functions of naturally occurring cilia is fluid transport. Biological cilia use spatially asymmetric strokes to generate a net fluid flow that can be utilized for feeding, swimming, and other functions. Biomimetic synthetic cilia with similar asymmetric beating can be useful for fluid manipulations in lab-on-chip devices. In this paper, we demonstrate the microfluidic pumping by magnetically actuated synthetic cilia arranged in multi-row arrays. We use a microchannel loop to visualize flow created by the ciliary array and to examine pumping for a range of cilia and microchannel parameters. We show that magnetic cilia can achieve flow rates of up to 11 {\mu}l/min with the pressure drop of~ 1 Pa. Such magnetic ciliary array can be useful in microfluidic applications requiring rapid and controlled fluid transport.

RANS Equations with Reynolds Stress Closure Can Be Ill-Conditioned

Abstract: Reynolds-averaged Navier-Stokes (RANS) simulations with turbulence models continue to play important roles in industrial flow simulations as high-fidelity simulations are prohibitively expensive for such flows. Commonly used linear eddy viscosity models are intrinsically unable to handle flows with non-equilibrium turbulence (e.g., flows with massive separation). Reynolds stress models, on the other hand, are plagued by their lack of robustness and stability. Recent studies found that even substituting Reynolds stresses from DNS databases (with errors below 0.5%) into RANS equations leads to grossly inaccurate velocities. Such an observation is not only disturbing for the recently emerging data-driven Reynolds stress models but also relevant for traditional, equation-based models. This observation cannot be explained by the global matrix condition number of the discretized RANS equations. In this work, we propose a metric based on local condition numbers for a priori evaluation of the stability of Reynolds stress models. Numerical tests on turbulent channel flows at various Reynolds numbers suggest that the proposed metric can adequately explain observations in previous studies, i.e., decreased model stability with increasing Reynolds number, and better stability of the implicit treatment of Reynolds stress compared to the explicit treatment.

A two-phase mixing layer between parallel gas and liquid streams: multiphase turbulence statistics and influence of interfacial instability

Abstract: The two-phase mixing layer formed between parallel gas and liquid streams is an important fundamental problem in turbulent multiphase flows. The problem is relevant to many industrial applications and natural phenomena, such as air-blast atomizers in fuel injection systems and breaking waves in the ocean. The velocity difference between the gas and liquid streams triggers an interfacial instability which can be convective or absolute depending on the stream properties and injection parameters. In the present study, a direct numerical simulation of a two-phase gas-liquid mixing layer that lie in the absolute instability regime is conducted. A dominant frequency is observed in the simulation and the numerical result agrees well with the prediction from viscous stability theory. As the interfacial wave plays a critical role in turbulence transition and development, the temporal evolution of turbulent fluctuations (such as the enstrophy) also exhibits a similar frequency. In order to investigate the statistical response of the multiphase turbulence flow, the simulation has been run for a long physical time so that time-averaging can be performed to yield the statistically converged results for Reynolds stresses and the turbulent kinetic energy (TKE) budget. An extensive mesh refinement study using from 8 million to about 4 billions cells has been carried out. The turbulent dissipation is shown to be highly demanding on mesh resolution compared to other terms in TKE budget. The results obtained with the finest mesh are shown to be not far from converged results of turbulent dissipation which allow us to obtain estimations of the Kolmogorov and Hinze scales. The computed Hinze scale is significantly larger than the size of droplets observed and does not seem to be a relevant length scale to describe the smallest size of droplets formed in atomization.

Review on the conversion of thermoacoustic power into electricity

Abstract: Thermoacoustic engines convert heat energy into high amplitude acoustic waves and subsequently into electric power. This article provides a review of the four main methods to convert the (thermo)acoustic power into electricity. First, loudspeakers and linear alternators are discussed in a section on electromagnetic devices. This is followed by sections on piezoelectric transducers, magnetohydrodynamic generators, and bidirectional turbines. Each segment provides a literature review of the given technology for the field of thermoacoustics, focusing on possible configurations, operating characteristics, output performance, and analytical and numerical methods to study the devices. This information is used as an input to discuss the performance and feasibility of each method, and to identify challenges that should be overcome for a more successful implementation in thermoacoustic engines. The work is concluded by a comparison of the four technologies, concentrating on the possible areas of application, the conversion efficiency, maximum electrical power output and more generally the suggested focus for future work in the field.

Dynamic slip wall model for large-eddy simulation

Abstract: Wall modelling in large-eddy simulation (LES) is necessary to overcome the prohibitive near-wall resolution requirements in high-Reynolds-number turbulent flows. Most existing wall models rely on assumptions about the state of the boundary layer and require a priori prescription of tunable coefficients. They also impose the predicted wall stress by replacing the no-slip boundary condition at the wall with a Neumann boundary condition in the wall-parallel directions while maintaining the no-transpiration condition in the wall-normal direction. In the present study, we first motivate and analyse the Robin (slip) boundary condition with transpiration (nonzero wall-normal velocity) in the context of wall-modelled LES. The effect of the slip boundary condition on the one-point statistics of the flow is investigated in LES of turbulent channel flow and flat-plate turbulent boundary layer. It is shown that the slip condition provides a framework to compensate for the deficit or excess of mean momentum at the wall. Moreover, the resulting nonzero stress at the wall alleviates the well-known problem of the wall-stress under-estimation by current subgrid-scale (SGS) models. Secondly, we discuss the requirements for the slip condition to be used in conjunction with wall models and derive the equation that connects the slip boundary condition with the stress at the wall. Finally, a dynamic procedure for the slip coefficients is formulated, providing a dynamic slip wall model free of a priori specified coefficients. The performance of the proposed dynamic wall model is tested in a series of LES of turbulent channel flow, non-equilibrium three-dimensional channel flow, and flat-plate turbulent boundary layer. The results show that the dynamic wall model is able to accurately predict one-point turbulence statistics for various flow configurations, Reynolds numbers, and grid resolutions.

Data-driven decomposition of the streamwise turbulence kinetic energy in boundary layers. Part 1. Energy spectra

Abstract: In wall-bounded turbulence, a multitude of coexisting turbulence structures form the streamwise velocity energy spectrum from the viscosity- to the inertia-dominated range of scales. Definite scaling-trends for streamwise spectra have remained empirically elusive, although a prominent school of thought stems from the works of Perry and Abell (J. Fluid Mech., vol. 79, 1977, pp. 785-799) and Perry et al. (J. Fluid Mech., vol. 165, 1986, pp. 163-199), which were greatly inspired by the attached-eddy hypothesis of Townsend (The Structure of Turbulent Shear Flow, Cambridge University Press, 1976). In this paper, we re-examine the turbulence kinetic energy of the streamwise velocity component in the context of the spectral decompositions of Perry and coworkers. Two universal spectral filters are derived via spectral coherence analysis of two-point velocity signals, spanning a Reynolds number range $Re_\tau \sim \mathcal{O}(10^3)$ to $\mathcal{O}(10^6)$ and form the basis for our decomposition of the logarithmic-region turbulence into stochastically wall-detached and wall-attached portions of energy. The latter is composed of scales larger than a streamwise/wall-normal ratio of $\lambda_x/z \approx 14$. If the decomposition is accepted, a $k_x^{-1}$ scaling region can only appear for $Re_\tau \geq 80\,000$, at a wall-normal position of $z^+ = 100$. Following Perry and co-workers, it is hypothesized that spectral contributions from turbulence structures other than attached eddies obscure a $k_x^{-1}$ scaling. When accepting the idea of different spectral contributions it is furthermore shown that a broad outer-spectral peak is present even at low $Re_\tau$.

A data-driven Koopman model predictive control framework for nonlinear flows

Abstract: The Koopman operator theory is an increasingly popular formalism of dynamical systems theory which enables analysis and prediction of the nonlinear dynamics from measurement data. Building on the recent development of the Koopman model predictive control framework (Korda and Mezic 2016), we propose a methodology for closed-loop feedback control of nonlinear flows in a fully data-driven and model-free manner. In the first step, we compute a Koopman-linear representation of the control system using a variation of the extended dynamic mode decomposition algorithm and then we apply model predictive control to the constructed linear model. Our methodology handles both full-state and sparse measurement; in the latter case, it incorporates the delay-embedding of the available data into the identification and control processes. We illustrate the application of this methodology on the periodic Burgers' equation and the boundary control of a cavity flow governed by the two-dimensional incompressible Navier-Stokes equations. In both examples the proposed methodology is successful in accomplishing the control tasks with sub-millisecond computation time required for evaluation of the control input in closed-loop, thereby allowing for a real-time deployment.

RANS Equations with Explicit Data-Driven Reynolds Stress Closure Can Be Ill-Conditioned

Abstract: Reynolds-averaged Navier--Stokes (RANS) simulations with turbulence closure models continue to play important roles in industrial flow simulations. However, the commonly used linear eddy viscosity models are intrinsically unable to handle flows with non-equilibrium turbulence. Reynolds stress models, on the other hand, are plagued by their lack of robustness. Recent studies in plane channel flows found that even substituting Reynolds stresses with errors below 0.5% from direct numerical simulation (DNS) databases into RANS equations leads to velocities with large errors (up to 35%). While such an observation may have only marginal relevance to traditional Reynolds stress models, it is disturbing for the recently emerging data-driven models that treat the Reynolds stress as an explicit source term in the RANS equations, as it suggests that the RANS equations with such models can be ill-conditioned. So far, a rigorous analysis of the condition of such models is still lacking. As such, in this work we propose a metric based on local condition number function for a priori evaluation of the conditioning of the RANS equations. We further show that the ill-conditioning cannot be explained by the global matrix condition number of the discretized RANS equations. Comprehensive numerical tests are performed on turbulent channel flows at various Reynolds numbers and additionally on two complex flows, i.e., flow over periodic hills and flow in a square duct. Results suggest that the proposed metric can adequately explain observations in previous studies, i.e., deteriorated model conditioning with increasing Reynolds number and better conditioning of the implicit treatment of Reynolds stress compared to the explicit treatment. This metric can play critical roles in the future development of data-driven turbulence models by enforcing the conditioning as a requirement on these models.

Periodic motion representing isotropic turbulence

Abstract: Temporally periodic solutions are extracted numerically from forced box turbulence with high symmetry. Since they are unstable to small perturbations, they are not found by forward integration but can be captured by Newton-Raphson iterations. Several periodic flows of various periods are identified for the micro-scale Reynolds number $R_{\lambda}$ between $50$ and $67$. The statistical properties of these periodic flows are compared with those of turbulent flow. It is found that the one with the longest period, which is two to three times the large-eddy-turnover time of turbulence, exhibits the same behaviour quantitatively as turbulent flow. In particular, we compare the energy spectrum, the Reynolds number dependence of the energy-dissipation rate, the pattern of the energy-cascade process, and the magnitude of the largest Lyapunov exponent. This periodic motion consists of high-activity and low-activity periods, which turbulence approaches, more often around its low-activity part, at the rate of once over a few eddy-turnover times. With reference to this periodic motion the Kaplan-York dimension and the Kolmogorov-Sinai entropy of the turbulence with high symmetry are estimated at $R_{\lambda}=67$ to be $19.7$ and $0.992$ respectively. The significance of such periodic solutions, embedded in turbulence, for turbulence analysis is discussed.

Modelling uncertainty using circulation-preserving stochastic transport noise in a 2-layer quasi-geostrophic model

Abstract: The stochastic variational approach for geophysical fluid dynamics was introduced by Holm (Proc Roy Soc A, 2015) as a framework for deriving stochastic parameterisations for unresolved scales The key feature of transport noise is that it respects the Kelvin circulation theorem This paper applies the variational stochastic parameterisation in a two-layer quasi-geostrophic model for a $\beta$-plane channel flow configuration The parameterisation is tested by comparing it with a deterministic high resolution eddy-resolving solution that has reached statistical equilibrium We describe a stochastic time-stepping scheme for the two-layer model and discuss its consistency in time Then we describe a procedure for estimating the stochastic forcing to approximate unresolved components using data from the high resolution deterministic simulation We compare an ensemble of stochastic solutions at lower resolution with the numerical solution of the deterministic model These computations quantify the uncertainty of the coarse grid computation relative to the fine grid computation The results show that the proposed parameterisation is efficient and effective for both homogeneous and heterogeneous flows, and they lay a solid foundation for data assimilation

Resolvent-analysis-based design of airfoil separation control

Abstract: We combine three-dimensional (3D) large-eddy simulations (LES) and resolvent analysis to design active separation control techniques on a NACA 0012 airfoil. Spanwise-periodic flows over the airfoil at a chord-based Reynolds number of $23,000$ and a free-stream Mach number of $0.3$ are considered at two post-stall angles of attack of $6^\circ$ and $9^\circ$. Near the leading edge, localized unsteady thermal actuation is introduced in an open-loop manner with two tunable parameters of actuation frequency and spanwise wavelength. For the most successful control case that achieves full reattachment, we observe a reduction in drag by up to $49\%$ and increase in lift by up to $54\%$. To provide physics-based guidance for the effective choice of these control input parameters, we conduct global resolvent analysis on the baseline turbulent mean flows to identify the actuation frequency and wavenumber that provide high energy amplification. The present analysis also considers the use of a temporal filter to limit the time horizon for assessing the energy amplification to extend resolvent analysis to unstable base flows. We incorporate the amplification and response mode from resolvent analysis to provide a metric that quantifies momentum mixing associated with the modal structure. By comparing this metric from resolvent analysis and the LES results of controlled flows, we demonstrate that resolvent analysis can predict the effective range of actuation frequency as well as the global response to the actuation input. Supported by the agreements between the results from resolvent analysis and LES, we believe that this study provides insights for the use of resolvent analysis in guiding future active flow control.

Tumbling, buckling, snaking: Morphological transitions of elastic filaments in shear flow

Abstract: The morphological dynamics, instabilities and transitions of elastic filaments in viscous flows underlie a wealth of biophysical processes from flagellar propulsion to intracellular streaming, and are also key to deciphering the rheological behavior of many complex fluids and soft materials. Here, we combine experiments and computational modeling to elucidate the dynamical regimes and morphological transitions of elastic Brownian filaments in a simple shear flow. Actin filaments are employed as an experimental model system and their conformations are investigated through fluorescence microscopy in microfluidic channels. Simulations matching the experimental conditions are also performed using inextensible Euler-Bernoulli beam theory and non-local slender-body hydrodynamics in the presence of thermal fluctuations, and agree quantitatively with observations. We demonstrate that filament dynamics in this system is primarily governed by a dimensionless elasto-viscous number comparing viscous drag forces to elastic bending forces, with thermal fluctuations only playing a secondary role. While short and rigid filaments perform quasi-periodic tumbling motions, a buckling instability arises above a critical flow strength. A second transition to strongly-deformed shapes occurs at a yet larger value of the elasto-viscous number and is characterized by the appearance of localized high-curvature bends that propagate along the filaments in apparent "snaking" motions. A theoretical model for the so far unexplored onset of snaking accurately predicts the transition and explains the observed dynamics. For the first time, we present a complete characterization of filament morphologies and transitions as a function of elasto-viscous number and scaled persistence length and demonstrate excellent agreement between theory, experiments and simulations.