Showing papers in "arXiv: Fluid Dynamics in 2017"
••
TL;DR: In this article, the authors describe wall-bounded turbulence as a deterministic high-dimensional dynamical system of interacting coherent structures, defined as eddies with enough internal dynamics to behave relatively autonomously from any remaining incoherent part of the flow.
Abstract: This article discusses the description of wall-bounded turbulence as a deterministic high-dimensional dynamical system of interacting coherent structures, defined as eddies with enough internal dynamics to behave relatively autonomously from any remaining incoherent part of the flow. The guiding principle is that randomness is not a property, but a methodological choice of what to ignore in the flow, and that a complete understanding of turbulence, including the possibility of control, requires that it be kept to a minimum. After briefly reviewing the underlying low-order statistics of flows at moderate Reynolds numbers, the article examines what two-point statistics imply for the decomposition of the flow into individual eddies. Intense eddies are examined next, including their temporal evolution, and shown to satisfy many of the properties required for coherence. In particular, it is shown that coherent structures larger than the Corrsin scale are a natural consequence of the shear. In wall-bounded turbulence, they can be classified into coherent dispersive waves and transient bursts. The former are found in the viscous layer near the wall and as very-large structures spanning the boundary layer thickness. Although they are shear-driven, these waves have enough internal structure to maintain a uniform advection velocity. Conversely, bursts exist at all scales, are characteristic of the logarithmic layer, and interact almost linearly with the shear. While the waves require a wall to determine their length scale, the bursts are essentially independent from it. The article concludes with a brief review of our present theoretical understanding of turbulent structures, and with a list of open problems and future perspectives.
268 citations
••
TL;DR: In this paper, the authors proposed a general dynamic reduced-order modeling framework for typical experimental data: time-resolved sensor data and optional non-time-resolution PIV snapshots.
Abstract: We propose a general dynamic reduced-order modeling framework for typical experimental data: time-resolved sensor data and optional non-time-resolved PIV snapshots. This framework contains four steps. First, the sensor signals are lifted to a dynamic feature space. Second, we identify a sparse human-interpretable nonlinear dynamical system for the feature state based on the sparse identification of nonlinear dynamics (SINDy). Third, if PIV snapshots are available, a local linear mapping from the feature state to velocity fields is shown to be orders of magnitudes more accurate than optimal modal expansions of the same order. Fourth, a generalized feature-based modal decomposition identifies coherent structures that are most dynamically correlated with the linear and nonlinear interaction terms in the sparse model, adding interpretability. Steps 1 and 2 define a black-box model. Optional steps 3 and 4 lift the black-box dynamics to a 'gray-box' model of the coherent structures, if non-time-resolved full-state data is available. This gray-box modeling strategy is successfully applied to the transient and post-transient laminar cylinder wake, and compares favorably with a POD model. We foresee numerous applications of this highly flexible modeling strategy, including estimation, prediction and control. Moreover, the feature space may be based on intrinsic coordinates, which are unaffected by a key challenge of modal expansion: the slow change of low-dimensional coherent structures with changing geometry and varying parameters.
144 citations
••
TL;DR: In this paper, a unified framework was proposed to describe the dynamics of fast liquid jets, called micro-jets, produced within cavitation bubbles experiencing an aspherical collapse.
Abstract: Fast liquid jets, called micro-jets, are produced within cavitation bubbles experiencing an aspherical collapse. Here we review micro-jets of different origins, scales and appearances, and propose a unified framework to describe their dynamics by using an anisotropy parameter $\zeta$, representing a dimensionless measure of the liquid momentum at the collapse point (Kelvin impulse). This parameter is rigorously defined for various jet drivers, including gravity and nearby boundaries. Combining theoretical considerations with hundreds of high-speed visualisations of bubbles collapsing near a rigid surface, near a free surface or in variable gravity, we classify the jets into three distinct regimes: weak, intermediate and strong. Weak jets ($\zeta 0.1$) pierce the bubble early during the collapse. The dynamics of the jets is analysed through key observables, such as the jet impact time, jet speed, bubble displacement, bubble volume at jet impact and vapour-jet volume. We find that, upon normalising these observables to dimensionless jet parameters, they all reduce to straightforward functions of $\zeta$, which we can reproduce numerically using potential flow theory. An interesting consequence of this result is that a measurement of a single observable, such as the bubble displacement, suffices to estimate any other parameter, such as the jet speed. Remarkably, the dimensionless parameters of intermediate and weak jets only depend on $\zeta$, not on the jet driver. In the same regime, the jet parameters are found to be well approximated by power-laws of $\zeta$, which we explain through analytical arguments.
130 citations
••
TL;DR: In this article, a complete PIML framework for predictive turbulence modeling, including learning Reynolds stress discrepancy function, predicting Reynolds stresses in different flows, and propagating to mean flow fields, is introduced.
Abstract: Although an increased availability of computational resources has enabled high-fidelity simulations of turbulent flows, the RANS models are still the dominant tools for industrial applications. However, the predictive capabilities of RANS models are limited by potential inaccuracy driven by hypotheses in the Reynolds stress closure. Recently, a Physics-Informed Machine Learning (PIML) approach has been proposed to learn the functional form of Reynolds stress discrepancy in RANS simulations based on available data. It has been demonstrated that the learned discrepancy function can be used to improve Reynolds stresses in different flows where data are not available. However, owing to a number of challenges, the improvements have been demonstrated only in the Reynolds stress prediction but not in the corresponding propagated quantities of interest. In this work, we introduce the procedures toward a complete PIML framework for predictive turbulence modeling, including learning Reynolds stress discrepancy function, predicting Reynolds stresses in different flows, and propagating to mean flow fields. The process of Reynolds stress propagation and predictive accuracy of the propagated velocity field are investigated. To improve the learning-prediction performance, the input features are enriched based on an integrity basis of invariants. The fully developed turbulent flow in a square duct is used as the test case. The discrepancy model is trained on flow fields obtained from several Reynolds numbers and evaluated on a duct flow at a Reynolds number higher than any of the training cases. The predicted Reynolds stresses are propagated to velocity field through RANS equations. Numerical results show excellent predictive performances in both Reynolds stresses and their propagated velocities, demonstrating the merits of the PIML approach in predictive turbulence modeling.
102 citations
••
TL;DR: In this article, the structure of turbulence in jets in the subsonic, transonic, and supersonic regimes was examined by LES data and resolvent analysis of the mean flow.
Abstract: Informed by LES data and resolvent analysis of the mean flow, we examine the structure of turbulence in jets in the subsonic, transonic, and supersonic regimes. Spectral (frequency-space) proper orthogonal decomposition is used to extract energy spectra and decompose the flow into energy-ranked coherent structures. The educed structures are generally well predicted by the resolvent analysis. Over a range of low frequencies and the first few azimuthal mode numbers, these jets exhibit a low-rank response characterized by Kelvin-Helmholtz (KH) type wavepackets associated with the annular shear layer up to the end of the potential core and that are excited by forcing in the very-near-nozzle shear layer. These modes too the have been experimentally observed before and predicted by quasi-parallel stability theory and other approximations--they comprise a considerable portion of the total turbulent energy. At still lower frequencies, particularly for the axisymmetric mode, and again at high frequencies for all azimuthal wavenumbers, the response is not low rank, but consists of a family of similarly amplified modes. These modes, which are primarily active downstream of the potential core, are associated with the Orr mechanism. They occur also as sub-dominant modes in the range of frequencies dominated by the KH response. Our global analysis helps tie together previous observations based on local spatial stability theory, and explains why quasi-parallel predictions were successful at some frequencies and azimuthal wavenumbers, but failed at others.
100 citations
•
TL;DR: The proposed CNN-based approximation procedure has a profound impact on the parametric design of bluff bodies and the feedback control of separated flows.
Abstract: We present an efficient deep learning technique for the model reduction of the Navier-Stokes equations for unsteady flow problems. The proposed technique relies on the Convolutional Neural Network (CNN) and the stochastic gradient descent method. Of particular interest is to predict the unsteady fluid forces for different bluff body shapes at low Reynolds number. The discrete convolution process with a nonlinear rectification is employed to approximate the mapping between the bluff-body shape and the fluid forces. The deep neural network is fed by the Euclidean distance function as the input and the target data generated by the full-order Navier-Stokes computations for primitive bluff body shapes. The convolutional networks are iteratively trained using the stochastic gradient descent method with the momentum term to predict the fluid force coefficients of different geometries and the results are compared with the full-order computations. We attempt to provide a physical analogy of the stochastic gradient method with the momentum term with the simplified form of the incompressible Navier-Stokes momentum equation. We also construct a direct relationship between the CNN-based deep learning and the Mori-Zwanzig formalism for the model reduction of a fluid dynamical system. A systematic convergence and sensitivity study is performed to identify the effective dimensions of the deep-learned CNN process such as the convolution kernel size, the number of kernels and the convolution layers. Within the error threshold, the prediction based on our deep convolutional network has a speed-up nearly four orders of magnitude compared to the full-order results and consumes an insignificant fraction of computational resources. The proposed CNN-based approximation procedure has a profound impact on the parametric design of bluff bodies and the feedback control of separated flows.
82 citations
••
TL;DR: In this article, a new definition of a new vector quantity called vortex vector is proposed to describe the local fluid rotation in 2D and 3D flows, and the magnitude of the vortex vector was defined as the rotational part of vorticity in the direction of vortex vector, which is twice of the minimum angular velocity of fluid around the point among all azimuth in the plane perpendicular to vortex vector.
Abstract: Vortex is ubiquitous in nature. However, there is not a consensus on the vortex definition in fluid dynamics. Lack of mathematical definition has caused considerable confusions in visualizing and understanding the coherent vortical structures in turbulence. According to previous study, it is realized that vortex is not the vorticity tube and vorticity should be decomposed into a rotational part which is the vortex vector and a non-rotational part which is the shear. In this paper, several new concepts such as fluid rotation of local point, the direction of fluid rotation axis and the strength of fluid rotation are proposed by investigating the kinematics of fluid element in the 2D and 3D flows. A definition of a new vector quantity called vortex vector is proposed to describe the local fluid rotation. The direction of the vortex vector is defined as the direction of local fluid rotation axis. The velocity components in the plane orthogonal to the vortex vector have zero derivatives along the vortex vector direction. The magnitude of the vortex vector is defined as the rotational part of vorticity in the direction of the vortex vector, which is the twice of the minimum angular velocity of fluid around the point among all azimuth in the plane perpendicular to vortex vector. According to the definition of the vortex vector, vortex is defined as a connected flow region where the magnitude of the vortex vector at each point is larger than zero. The new definition for the vortex vector and vortex follows three principles: 1. Local in quantity, 2. Galilean invariant, 3. Unique. The definitions are carefully checked by DNS and LES examples which clearly show that the new defined vortex vector and vortex can fully represent the complex structures of vortices in turbulence.
82 citations
••
TL;DR: In this paper, a series of high-resolution three dimensional smulations of two initial conditions with eight independent codes were used to investigate the turbulent Richtmyer-Meshkov instability.
Abstract: Turbulent Richtmyer--Meshkov instability (RMI) is investigated through a series of high resolution three dimensional smulations of two initial conditions with eight independent codes. The simulations are initialised with a narrowband perturbation such that instability growth is due to non-linear coupling/backscatter from the energetic modes, thus generating the lowest expected growth rate from a pure RMI. By independently assessing the results from each algorithm, and computing ensemble averages of multiple algorithms, the results allow a quantification of key flow properties as well as the uncertainty due to differing numerical approaches. A new analytical model predicting the initial layer growth for a multimode narrowband perturbation is presented, along with two models for the linear and non-linear regime combined. Overall, the growth rate exponent is determined as $\theta=0.292 \pm 0.009$, in good agreement with prior studies; however, the exponent is decaying slowly in time. $\theta$ is shown to be relatively insensitive to the choice of mixing layer width measurement. The asymptotic integral molecular mixing measures $\Theta=0.792\pm 0.014$, $\Xi=0.800 \pm 0.014$ and $\Psi=0.782\pm 0.013$ which are lower than some experimental measurements but within the range of prior numerical studies. The flow field is shown to be persistently anisotropic for all algorithms, at the latest time having between 49\% and 66\% higher kinetic energy in the shock parallel direction compared to perpendicular and does not show any return to isotropy. The plane averaged volume fraction profiles at different time instants collapse reasonably well when scaled by the integral width, implying that the layer can be described by a single length scale and thus a single $\theta$. Quantitative data given for both ensemble averages and individual algorithms provide useful benchmark results for future research.
79 citations
••
TL;DR: In this article, the authors demonstrate how liquid-infused surfaces can reduce turbulent drag significantly in Taylor-Couette flow, achieving a reduction in the amount of turbulent drag exceeding 35%.
Abstract: Experiments are presented that demonstrate how liquid-infused surfaces can reduce turbulent drag significantly in Taylor-Couette flow. The test liquid was water, and the test surface was composed of square microscopic grooves measuring 100 $\mu$m to 800 $\mu$m, filled with alkane liquids with viscosities from 0.3 to 1.4 times that of water. We achieve drag reduction exceeding 35\%, four times higher than previously reported for liquid-infused surfaces in turbulent flow. The level of drag reduction increased with viscosity ratio, groove width, fluid area fraction, and Reynolds number. The optimum groove width was given by $w^+ \approx 35$.
66 citations
••
TL;DR: In this paper, an analytical model is developed for the prediction of noise radiated by an aerofoil with leading edge serration in a subsonic turbulent stream, which makes use of the Fourier expansion and Schwarzschild techniques in order to solve a set of coupled differential equations iteratively and express the far-field sound power spectral density in terms of the statistics of incoming turbulent upwash velocity.
Abstract: An analytical model is developed for the prediction of noise radiated by an aerofoil with leading edge serration in a subsonic turbulent stream. The model makes use of the Fourier Expansion and Schwarzschild techniques in order to solve a set of coupled differential equations iteratively and express the far-field sound power spectral density in terms of the statistics of incoming turbulent upwash velocity. The model has shown that the primary noise reduction mechanism is due to the destructive interference of the scattered pressure induced by the leading edge serrations. It has also shown that in order to achieve significant sound reduction, the serration must satisfy two geometrical criteria related to the serration sharpness and hydrodynamic properties of the turbulence. A parametric study has been carried out and it is shown that serrations can reduce the overall sound pressure level at most radiation angles, particularly at downstream angles close to the aerofoil surface. The sound directivity results have also shown that the use of leading edge serration does not particularity change the dipolar pattern of the far-field noise at low frequencies, but it changes the cardioid directivity pattern associated with radiation from straight-edge scattering at high frequencies to a tilted dipolar pattern.
65 citations
••
TL;DR: In this paper, numerical simulations of a turbulent channel flow over an hyper-elastic wall were performed using a one-continuum formulation, using a monolithic velocity field for both the fluid and solid phase, allowing the use of a fully Eulerian formulation.
Abstract: We perform numerical simulations of a turbulent channel flow over an hyper-elastic wall In the fluid region the flow is governed by the incompressible Navier-Stokes (NS) equations, while the solid is a neo-Hookean material satisfying the incompressible Mooney-Rivlin law The multiphase flow is solved with a one-continuum formulation, using a monolithic velocity field for both the fluid and solid phase, which allows the use of a fully Eulerian formulation The simulations are carried out at Reynolds bulk $Re=2800$ and examine the effect of different elasticity and viscosity of the deformable wall We show that the skin friction increases monotonically with the material elastic modulus The turbulent flow in the channel is affected by the moving wall even at low values of elasticity since non-zero fluctuations of vertical velocity at the interface influence the flow dynamics The near-wall streaks and the associated quasi-streamwise vortices are strongly reduced near a highly elastic wall while the flow becomes more correlated in the spanwise direction, similarly to what happens for flows over rough and porous walls As a consequence, the mean velocity profile in wall units is shifted downwards when shown in logarithmic scale, and the slope of the inertial range increases in comparison to that for the flow over a rigid wall We propose a correlation between the downward shift of the inertial range, its slope and the wall-normal velocity fluctuations at the wall, extending results for the flow over rough walls We finally show that the interface deformation is determined by the fluid fluctuations when the viscosity of the elastic layer is low, while when this is high the deformation is limited by the solid properties
••
TL;DR: In this paper, the Koopman mode decomposition (KMD) was used to study the dynamics of the lid-driven flow in a two-dimensional square cavity based on theorems related to the spectral theory of the koopman operator, which is a data-analysis technique which is often used to extract the spatio-temporal patterns of complex flows.
Abstract: The Koopman Mode Decomposition (KMD) is a data-analysis technique which is often used to extract the spatio-temporal patterns of complex flows. In this paper, we use KMD to study the dynamics of the lid-driven flow in a two-dimensional square cavity based on theorems related to the spectral theory of the Koopman operator. We adapt two algorithms, from the classical Fourier and power spectral analysis, to compute the discrete and continuous spectrum of the Koopman operator for the post-transient flows. Properties of the Koopman operator spectrum are linked to the sequence of flow regimes occurring between $Re=10000$ and $Re=30000$, and changing the flow nature from steady to aperiodic. The Koopman eigenfunctions for different flow regimes, including flows with mixed spectra, are constructed using the assumption of ergodicity in the state space. The associated Koopman modes show remarkable robustness even as the temporal nature of the flow is changing substantially. We observe that KMD outperforms the Proper Orthogonal Decomposition in reconstruction of the flows with strong quasi-periodic components.c features are present in the flow.
••
TL;DR: In this article, the authors present detailed observations of the shock waves emitted at the collapse of single cavitation bubbles using simultaneous time-resolved shadowgraphy and hydrophone pressure measurements.
Abstract: We present detailed observations of the shock waves emitted at the collapse of single cavitation bubbles using simultaneous time-resolved shadowgraphy and hydrophone pressure measurements. The geometry of the bubbles is systematically varied from spherical to very non-spherical by decreasing their distance to a free or rigid surface or by modulating the gravity-induced pressure gradient aboard parabolic flights. The non-spherical collapse produces multiple shocks that are clearly associated with different processes, such as the jet impact and the individual collapses of the distinct bubble segments. For bubbles collapsing near a free surface, the energy and timing of each shock are measured separately as a function of the anisotropy parameter $\zeta$, which represents the dimensionless equivalent of the Kelvin impulse. For a given source of bubble deformation (free surface, rigid surface or gravity), the normalized shock energy depends only on $\zeta$, irrespective of the bubble radius $R_{0}$ and driving pressure $\Delta p$. Based on this finding, we develop a predictive framework for the peak pressure and energy of shock waves from non-spherical bubble collapses. Combining statistical analysis of the experimental data with theoretical derivations, we find that the shock peak pressures can be estimated as jet impact-induced hammer pressures, expressed as $p_{h} = 0.45\left(\rho c^{2}\Delta p\right)^{1/2} \zeta^{-1}$ at $\zeta > 10^{-3}$. The same approach is found to explain the shock energy quenching as a function of $\zeta^{-2/3}$.
••
TL;DR: In this paper, the authors present results from large eddy simulations (LES) of extended wind-farms for several turbine configurations with a range of different spanwise and streamwise spacing combinations.
Abstract: We present results from large eddy simulations (LES) of extended wind-farms for several turbine configurations with a range of different spanwise and streamwise spacing combinations. The results show that for wind-farms arranged in a staggered configuration with spanwise spacings in the range $\approx[3.5,8]$D, where $D$ is the turbine diameter, the power output in the fully developed regime depends primarily on the geometric mean of the spanwise and streamwise turbine spacings. In contrast, for the aligned configuration the power output in the fully developed regime strongly depends on the streamwise turbine spacing and shows weak dependence on the spanwise spacing. Of interest to the rate of wake recovery, we find that the power output is well correlated with the vertical kinetic energy flux, which is a measure of how much kinetic energy is transferred into the wind-turbine region by the mean flow. A comparison between the aligned and staggered configurations reveals that the vertical kinetic energy flux is more localized along turbine rows for aligned wind-farms than for staggered ones. This additional mixing leads to a relatively fast wake recovery for aligned wind-farms.
••
TL;DR: In this article, the authors solved the Bhatnagar-gross-Krook equation in simplified porous media and identified two key factors that can explain Klinkenberg's experimental results: the tortuous flow path and the non-unitary tangential momentum accommodation coefficient for the gas-surface interaction.
Abstract: The apparent gas permeability of the porous medium is an important parameter in the prediction of unconventional gas production, which was first investigated systematically by Klinkenberg in 1941 and found to increase with the reciprocal mean gas pressure (or equivalently, the Knudsen number). Although the underlying rarefaction effects are well-known, the reason that the correction factor in Klinkenberg's famous equation decreases when the Knudsen number increases has not been fully understood. Most of the studies idealize the porous medium as a bundle of straight cylindrical tubes, however, according to the gas kinetic theory, this only results in an increase of the correction factor with the Knudsen number, which clearly contradicts Klinkenberg's experimental observations. Here, by solving the Bhatnagar-Gross-Krook equation in simplified (but not simple) porous media, we identify, for the first time, two key factors that can explain Klinkenberg's experimental results: the tortuous flow path and the non-unitary tangential momentum accommodation coefficient for the gas-surface interaction. Moreover, we find that Klinkenberg's results can only be observed when the ratio between the apparent and intrinsic permeabilities is $\lesssim30$; at large ratios (or Knudsen numbers) the correction factor increases with the Knudsen number. Our numerical results could also serve as benchmarking cases to assess the accuracy of macroscopic models and/or numerical schemes for the modeling/simulation of rarefied gas flows in complex geometries over a wide range of gas rarefaction.
••
TL;DR: In this article, a self-similar adverse pressure gradient (APG) turbulent boundary layer (TBL) at the verge of separation is simulated using direct numerical simulation (DNS) on a flat surface with farfield boundary conditions tailored to apply the desired pressure gradient.
Abstract: The statistical properties are presented for the direct numerical simulation (DNS) of a self-similar adverse pressure gradient (APG) turbulent boundary layer (TBL) at the verge of separation. The APG TBL has a momentum thickness based Reynolds number range from $Re_{\delta_2}=570$ to $13800$, with a self-similar region from $Re_{\delta_2} = 10000$ to $12300$. Within this domain the average non-dimensional pressure gradient parameter $\beta=39$, where for a unit density $\beta=\delta_1 P_e^\prime / \tau_w$, with $\delta_1$ the displacement thickness, $\tau_w$ the mean shear stress at the wall, and $P_e^\prime$ the farfield pressure gradient. This flow is compared to previous zero pressure gradient (ZPG) and mild APG TBL ($\beta=1$) results of similar Reynolds number. All flows are generated via the DNS of a TBL on a flat surface with farfield boundary conditions tailored to apply the desired pressure gradient. The conditions for self-similarity, and the appropriate length and velocity scales are derived. The mean and Reynolds stress profiles are shown to collapse when non-dimensionalised on the basis of these length and velocity scales. As the pressure gradient increases the flow has properties less like a ZPG TBL and more akin to a free shear layer.
••
TL;DR: In this paper, the authors studied laminar and turbulent dynamos in chiral magnetohydrodynamics with an extended set of equations that accounts for an additional contribution to the electric current due to the chiral magnetic effect (CME).
Abstract: Using direct numerical simulations (DNS), we study laminar and turbulent dynamos in chiral magnetohydrodynamics (MHD) with an extended set of equations that accounts for an additional contribution to the electric current due to the chiral magnetic effect (CME). This quantum phenomenon originates from an asymmetry between left- and right-handed relativistic fermions in the presence of a magnetic field and gives rise to a chiral dynamo. We show that the magnetic field evolution proceeds in three stages: (1) a small-scale chiral dynamo instability; (2) production of chiral magnetically driven turbulence and excitation of a large-scale dynamo instability due to a new chiral effect (alpha_mu effect); and (3) saturation of magnetic helicity and magnetic field growth controlled by a conservation law for the total chirality. The $\alpha_\mu$ effect becomes dominant at large fluid and magnetic Reynolds numbers and is not related to kinetic helicity. The growth rate of the large-scale magnetic field and its characteristic scale measured in the numerical simulations agree well with theoretical predictions based on mean-field theory. The previously discussed two-stage chiral magnetic scenario did not include stage (2) during which the characteristic scale of magnetic field variations can increase by many orders of magnitude. Based on the findings from numerical simulations, the relevance of the CME and the chiral effects revealed in the relativistic plasma of the early universe and of proto-neutron stars are discussed.
•
TL;DR: The novel DDF-ROM framework is built on general ideas of spatial filtering and optimization and is independent of (restrictive) phenomenological arguments and is significantly more accurate than the standard projection ROM.
Abstract: We propose a data-driven filtered reduced order model (DDF-ROM) framework for the numerical simulation of fluid flows. The novel DDF-ROM framework consists of two steps: (i) In the first step, we use explicit ROM spatial filtering of the nonlinear PDE to construct a filtered ROM. This filtered ROM is low-dimensional, but is not closed (because of the nonlinearity in the given PDE). (ii) In the second step, we use data-driven modeling to close the filtered ROM, i.e., to model the interaction between the resolved and unresolved modes. To this end, we use a quadratic ansatz to model this interaction and close the filtered ROM. To find the new coefficients in the closed filtered ROM, we solve an optimization problem that minimizes the difference between the full order model data and our ansatz. We emphasize that the new DDF-ROM is built on general ideas of spatial filtering and optimization and is independent of (restrictive) phenomenological arguments.
We investigate the DDF-ROM in the numerical simulation of a 2D channel flow past a circular cylinder at Reynolds number $Re=100$. The DDF-ROM is significantly more accurate than the standard projection ROM. Furthermore, the computational costs of the DDF-ROM and the standard projection ROM are similar, both costs being orders of magnitude lower than the computational cost of the full order model. We also compare the new DDF-ROM with modern ROM closure models in the numerical simulation of the 1D Burgers equation. The DDF-ROM is more accurate and significantly more efficient than these ROM closure models.
••
TL;DR: In this article, the authors explored the initiation of thin-film hydraulic jumps and found that the jump is formed when surface tension and viscous forces balance the momentum in the film and gravity plays no significant role.
Abstract: For more than a century, it has been believed that all hydraulic jumps are created due to gravity. However, we found that thin-film hydraulic jumps are not induced by gravity. This study explores the initiation of thin-film hydraulic jumps. For circular jumps produced by the normal impingement of a jet onto a solid surface, we found that the jump is formed when surface tension and viscous forces balance the momentum in the film and gravity plays no significant role. Experiments show no dependence on the orientation of the surface and a scaling relation balancing viscous forces and surface tension collapses the experimental data. Experiments on thin film planar jumps in a channel also show that the predominant balance is with surface tension, although for the thickness of the films we studied gravity also played a role in the jump formation. A theoretical analysis shows that the downstream transport of surface tension energy is the previously neglected, critical ingredient in these flows and that capillary waves play the role of gravity waves in a traditional jump in demarcating the transition from the supercritical to subcritical flow associated with these jumps.
••
TL;DR: In this paper, the Cauchy invariants equation and formula are shown to be duals of the Hodge dual of the Lagrangian map and the Jacobian of the ideal flow.
Abstract: Cauchy invariants are now viewed as a powerful tool for investigating the Lagrangian structure of three-dimensional (3D) ideal flow (Frisch & Zheligovsky, Commun. Math. Phys., vol. 326, 2014, pp. 499-505, Podvigina et al., J. Comput. Phys., vol. 306, 2016, pp. 320-342). Looking at such invariants with the modern tools of differential geometry and of geodesic flow on the space SDiff of volume-preserving transformations (Arnold, Ann. Inst. Fourier, vol. 16, 1966, pp. 319-361), all manners of generalisations are here derived. The Cauchy invariants equation and the Cauchy formula, relating the vorticity and the Jacobian of the Lagrangian map, are shown to be two expressions of this Lie-advection invariance, which are duals of each other (specifically, Hodge dual). Actually, this is shown to be an instance of a general result, which holds for flow both in flat (Euclidean) space and in a curved Riemannian space: any Lie-advection invariant p-form which is exact (i.e. is a differential of a (p-1)-form) has an associated Cauchy invariants equation and a Cauchy formula. This constitutes a new fundamental result in linear transport theory, providing a Lagrangian formulation of Lie advection for some classes of differential forms. The result has a broad applicability: examples include the magnetohydrodynamics (MHD) equations and various extensions thereof, discussed by Lingam et al. (Phys. Lett. A, vol. 380, 2016, pp. 2400-2406), and include also the equations of Tao (2016, arXiv:1606.08481 [math.AP]), Euler equations with modified Biot-Savart law, displaying finite-time blow-up. Our main result is also used for new derivations, and several new results, concerning local helicity-type invariants for fluids and MHD flow in flat or curved spaces of arbitrary dimension.
••
TL;DR: In this paper, the authors demonstrate that appropriate distortions to the velocity profile lead to a complete collapse of turbulence and subsequently friction losses are reduced by as much as 95% in pipe flow.
Abstract: Turbulence is the major cause of friction losses in transport processes and it is responsible for a drastic drag increase in flows over bounding surfaces. While much effort is invested into developing ways to control and reduce turbulence intensities, so far no methods exist to altogether eliminate turbulence if velocities are sufficiently large. We demonstrate for pipe flow that appropriate distortions to the velocity profile lead to a complete collapse of turbulence and subsequently friction losses are reduced by as much as 95%. Counterintuitively, the return to laminar motion is accomplished by initially increasing turbulence intensities or by transiently amplifying wall shear. The usual measures of turbulence levels, such as the Reynolds number (Re) or shear stresses, do not account for the subsequent relaminarization. Instead an amplification mechanism measuring the interaction between eddies and the mean shear is found to set a threshold below which turbulence is suppressed beyond recovery.
••
TL;DR: In this article, the influence that the kinematics of 2D airfoils have on the aerodynamic forces is investigated using Direct Numerical Simulations and a force decomposition algorithm.
Abstract: The influence that the kinematics of pitching and heaving 2D airfoils have on the aerodynamic forces is investigated using Direct Numerical Simulations and a force decomposition algorithm. Large amplitude motions are considered (of the order of one chord), with moderate Reynolds numbers and reduced frequencies of order 1, varying the mean pitch angle and the phase shift between the pitching and heaving motions. Our results show that the surface vorticity contribution (viscous effects) to the aerodynamic force is negligible compared to the contributions from the body motion (fluid inertia) and the vorticity within the flow (circulation). For the range of parameters considered here, the latter tends to be instantaneously oriented in the direction normal to the chord of the airfoil. Based on the results discussed in the paper, a reduced order model for the instantaneous aerodynamic force is proposed, taking advantage of the force decomposition and the chord-normal orientation of the contribution from vorticity within the flow to the total aerodynamic force. The predictions of the proposed model are compared to those of a similar model from the literature, showing a noticeable improvement on the prediction of the mean thrust, and a smaller improvement on the prediction of mean lift and the instantaneous force coefficients.
••
TL;DR: In this paper, the effects of finite flat porous extensions to semi-infinite impermeable flat plates in an attempt to control trailing-edge noise through bio-inspired adaptations are investigated.
Abstract: This paper investigates the effects of finite flat porous extensions to semi-infinite impermeable flat plates in an attempt to control trailing-edge noise through bio-inspired adaptations. Specifically the problem of sound generated by a gust convecting in uniform mean steady flow scattering off the trailing edge and permeable-impermeable junction is considered. This setup supposes that any realistic trailing-edge adaptation to a blade would be sufficiently small so that the turbulent boundary layer encapsulates both the porous edge and the permeable-impermeable junction, and therefore the interaction of acoustics generated at these two discontinuous boundaries is important. The acoustic problem is tackled analytically through use of the Wiener-Hopf method. A two-dimensional matrix Wiener-Hopf problem arises due to the two interaction points (the trailing edge and the permeable-impermeable junction). This paper discusses a new iterative method for solving this matrix Wiener-Hopf equation which extends to further two-dimensional problems in particular those involving analytic terms that exponentially grow in the upper or lower half planes. This method is an extension of the commonly used "pole removal" technique and avoids the needs for full matrix factorisation. Convergence of this iterative method to an exact solution is shown to be particularly fast when terms neglected in the second step are formally smaller than all other terms retained. The final acoustic solution highlights the effects of the permeable-impermeable junction on the generated noise, in particular how this junction affects the far-field noise generated by high-frequency gusts by creating an interference to typical trailing-edge scattering. This effect results in partially porous plates predicting a lower noise reduction than fully porous plates when compared to fully impermeable plates.
••
TL;DR: In this article, the authors studied the mechanisms of dispersion in the laminar flow through the pore space of a $3$-dimensional porous medium and quantified systematically in the framework of a time-domain random walk for the motion of marked elements (particles) of the transported material quantity.
Abstract: This paper studies the mechanisms of dispersion in the laminar flow through the pore space of a $3$-dimensional porous medium. We focus on pre-asymptotic transport prior to the asymptotic hydrodynamic dispersion regime, in which solute motion may be described by the average flow velocity and a hydrodynamic dispersion coefficient. High performance numerical flow and transport simulations of solute breakthrough at the outlet of a sand-like porous medium evidence marked deviations from the hydrodynamic dispersion paradigm and identifies two distinct regimes. The first regime is characterized by a broad distribution of advective residence times in single pores. The second regime is characterized by diffusive mass transfer into low-velocity region in the wake of solute grains. These mechanisms are quantified systematically in the framework of a time-domain random walk for the motion of marked elements (particles) of the transported material quantity. The model is parameterized with the characteristics of the porous medium under consideration and captures both pre-asymptotic regimes. Macroscale transport is described by an integro-differential equation for solute concentration, whose memory kernels are given in terms of the distribution of mean pore velocities and trapping times. This approach quantifies the physical non-equilibrium caused by a broad distribution of mass transfer time scales, both advective and diffusive, on the representative elementary volume (REV). Thus, while the REV indicates the scale at which medium properties like porosity can be uniquely defined, this does not imply that transport can be characterized by hydrodynamic dispersion.
•
TL;DR: It is demonstrated that in this framework, certain wave-like solutions exhibit low-rank structure and thus, can be efficiently compressed using relatively few global basis in the Lagrangian computational domain.
Abstract: Foundations of a new projection-based model reduction approach for convection dominated nonlinear fluid flows are summarized. In this method the evolution of the flow is approximated in the Lagrangian frame of reference. Global basis functions are used to approximate both the state and the position of the Lagrangian computational domain. It is demonstrated that in this framework, certain wave-like solutions exhibit low-rank structure and thus, can be efficiently compressed using relatively few global basis. The proposed approach is successfully demonstrated for the reduction of several simple but representative problems.
••
TL;DR: In this article, a model of thin layer turbulence that follows the evolution of the two-dimensional motions along the horizontal directions coupled to a single Fourier mode along the vertical direction was investigated.
Abstract: We investigate a model of thin layer turbulence that follows the evolution of the two-dimensional motions ${\bf u}_{_{2D}} (x,y)$ along the horizontal directions $(x,y)$ coupled to a single Fourier mode along the vertical direction ($z$) of the form ${\bf u}_q (x, y, z)=[v_x(x,y) \sin(qz), v_y(x,y)\sin(qz), v_z(x,y)\cos(qz)\, ]$, reducing thus the system to two coupled, two-dimensional equations. The reduced dimensionality of the model allows a thorough investigation of the transition from a forward to an inverse cascade of energy as the thickness of the layer $H=\pi/q$ is varied. Starting from a thick layer and reducing its thickness it is shown that two critical heights are met (i) one for which the forward unidirectional cascade (similar to three-dimensional turbulence) transitions to a bidirectional cascade transferring energy to both small and large scales and (ii) one for which the bidirectional cascade transitions to a unidirectional inverse cascade when the layer becomes very thin (similar to two-dimensional turbulence). The two critical heights are shown to have different properties close to criticality that we are able to analyze with numerical simulations for a wide range of Reynolds numbers and aspect ratios.
••
TL;DR: In this paper, the role of small elastic scales in the mechanism of self-sustained EIT is investigated by introducing global artificial diffusion in the hyperbolic transport equation for polymers.
Abstract: The goal of the present study is: (i) to demonstrate the two-dimensional nature of the elasto-inertial instability in elasto-inertial turbulence (EIT), (ii) to identify the role of the bi-dimensional instability in three-dimensional EIT flows and (iii) to establish the role of the small elastic scales in the mechanism of self-sustained EIT. Direct numerical simulations of FENE-P fluid flows are performed in two- and three-dimensional channels. The Reynolds number is set to $\mathrm{Re}_\tau = 85$ which is sub-critical for 2D flows but beyond transition for 3D ones. The polymer properties correspond to those of typical dilute polymer solutions and two moderate Weissenberg numbers, $\mathrm{Wi}_\tau = 40, 100$, are considered. The simulation results show that sustained turbulence can be observed in 2D sub-critical flows, confirming the existence of a bi-dimensional elasto-inertial instability. The same type of instability is also observed in 3D simulations where both Newtonian and elasto-inertial turbulent structures co-exist. Depending on the Wi number, one type of structure can dominate and drive the flow. For large Wi values, the elasto-inertial instability tends to prevail over the Newtonian turbulence. This statement is supported by (i) the absence of the typical Newtonian near-wall vortices and (ii) strong similarities between two- and three-dimensional flows when considering larger Wi numbers. The role of the small elastic scales is investigated by introducing global artificial diffusion in the hyperbolic transport equation for polymers. The study results show that the introduction of large polymer diffusion in the system strongly damps a significant part of the elastic scales that are necessary to feed turbulence, eventually leading to the flow laminarization. A sufficiently high Schmidt number is necessary to allow self-sustained turbulence to settle.
••
TL;DR: In this article, the authors used volume-of-fluid (VOF) method to perform three-dimensional numerical simulations of droplet formation of Newtonian fluids in microfluidic T-junction devices.
Abstract: We used volume-of-fluid (VOF) method to perform three-dimensional numerical simulations of droplet formation of Newtonian fluids in microfluidic T-junction devices. To evaluate the performance of the VOF method we examined the regimes of drop formation and determined droplet size as a function of system parameters. Comparison of the simulation results with four sets of experimental data from the literature showed good agreement, validating the VOF method. Motivated by the lack of adequate studies investigating the influence of viscosity ratio ({\lambda}) on the generated droplet size, we mapped the dependence of drop volume on capillary number (0.001 1. In addition, we find that at a given capillary number, the size of droplets does not vary appreciably when {\lambda} 1. We develop an analytical model for predicting droplet size that includes a viscosity-dependent breakup time for the dispersed phase. This improved model successfully predicts the effects of viscosity ratio observed in simulations. Results from this study are useful for the design of lab-on-chip technologies and manufacture of microfluidic emulsions, where there is a need to know how system parameters influence droplet size.
•
TL;DR: In this paper, the effect of airfoil thickness on the onset of dynamic stall was investigated using large eddy simulations at chord-based Reynolds number of 200,000, where the 3-D Navier Stokes solver, FDL3DI was used with a sixth-order compact finite difference scheme for spatial discretization, second-order implicit time integration, and discriminating filters to remove unresolved wavenumbers.
Abstract: Effect of airfoil thickness on onset of dynamic stall is investigated using large eddy simulations at chord-based Reynolds number of 200,000. Four symmetric NACA airfoils of thickness-to-chord ratios of 9%, 12%, 15%, and 18% are studied. The 3-D Navier Stokes solver, FDL3DI is used with a sixth-order compact finite difference scheme for spatial discretization, second-order implicit time integration, and discriminating filters to remove unresolved wavenumbers. A constant-rate pitch-up maneuver is studied with the pitching axis located at the airfoil quarter chord point. Simulations are performed in two steps. In the first step, the airfoil is kept static at a prescribed angle of attack ($=4^\circ$). In the second step, a ramp function is used to smoothly increase the pitch rate from zero to the selected value and then the pitch rate is held constant until the angle of attack goes past the lift stall point. Comparisons against XFOIL for the static simulations show good agreement in predicting the transition location. FDL3DI predicts two-stage transition for thin airfoils (9% and 12%), which is not observed in the XFOIL results. The dynamic simulations show that the onset of dynamic stall is marked by the bursting of the laminar separation bubble (LSB) in all cases. However, for the thickest airfoil tested, the reverse flow region spreads over most of the airfoil and reaches the LSB location immediately before the LSB bursts and dynamic stall begins, suggesting that stall could be triggered by the separated turbulent boundary layer. The results suggest that the boundary between different classifications of dynamic stall, particularly leading edge stall versus trailing edge stall are blurred. The dynamic stall onset mechanism changes gradually from one to the other with a gradual change in some parameters, in this case, airfoil thickness.
••
TL;DR: This study combines asymptotic predictions, laboratory and numerical results, and experimental constraints to build a heuristic framework for cross-comparison between a broad range of rotating convection studies, demonstrating that laboratory studies may soon be able to characterize geophysically-relevant flow regimes.
Abstract: Many geophysical and astrophysical phenomena are driven by turbulent fluid dynamics, containing behaviors separated by tens of orders of magnitude in scale. While direct simulations have made large strides toward understanding geophysical systems, such models still inhabit modest ranges of the governing parameters that are difficult to extrapolate to planetary settings. The canonical problem of rotating Rayleigh-Benard convection provides an alternate approach - isolating the fundamental physics in a reduced setting. Theoretical studies and asymptotically-reduced simulations in rotating convection have unveiled a variety of flow behaviors likely relevant to natural systems, but still inaccessible to direct simulation. In lieu of this, several new large-scale rotating convection devices have been designed to characterize such behaviors. It is essential to predict how this potential influx of new data will mesh with existing results. Surprisingly, a coherent framework of predictions for extreme rotating convection has not yet been elucidated. In this study, we combine asymptotic predictions, laboratory and numerical results, and experimental constraints to build a heuristic framework for cross-comparison between a broad range of rotating convection studies. We categorize the diverse field of existing predictions in the context of asymptotic flow regimes. We then consider the physical constraints that determine the points of intersection between flow behavior predictions and experimental accessibility. Applying this framework to several upcoming devices demonstrates that laboratory studies may soon be able to characterize geophysically-relevant flow regimes. These new data may transform our understanding of geophysical and astrophysical turbulence, and the conceptual framework developed herein should provide the theoretical infrastructure needed for meaningful discussion of these results.