No Arabic abstract
Wakes of aircraft and automobiles with relatively flat slanted aftbodies are often characterized by a streamwise-oriented vortex pair, whose strength affects drag and other crucial performance parameters. We examine the stability characteristics of the vortex pair emerging over an abstraction comprised of a streamwise-aligned cylinder terminated with an upswept plane. The Reynolds number is fixed at 5000 and the upsweep angle is increased from 20deg to 32deg. At 20deg, the LES yields a steady streamwise-oriented vortex pair, and the global modes are also stable. At 32deg, the LES displays unsteady flow behavior. Linear analysis of the mean flow reveals different unstable modes. The lowest oscillation frequency is an antisymmetric mode, which is attached to the entire slanted base. At the highest frequency, the mode is symmetric and has the same rotational orientation as the mean vortex pair. Its support is prominent in the rear part of the slanted base and spreads relatively rapidly downstream with prominent helical structures. A receptivity analysis of low- and high-frequency modes suggests the latter holds promise to affect the vortical flow, providing a potential starting point for a control strategy to modify the vortex pair.
We report the results of a complete modal and nonmodal linear stability analysis of the electrohydrodynamic flow (EHD) for the problem of electroconvection in the strong injection region. Convective cells are formed by Coulomb force in an insulating liquid residing between two plane electrodes subject to unipolar injection. Besides pure electroconvection, we also consider the case where a cross-flow is present, generated by a streamwise pressure gradient, in the form of a laminar Poiseuille flow. The effect of charge diffusion, often neglected in previous linear stability analyses, is included in the present study and a transient growth analysis, rarely considered in EHD, is carried out. In the case without cross-flow, a non-zero charge diffusion leads to a lower linear stability threshold and thus to a more unstable low. The transient growth, though enhanced by increasing charge diffusion, remains small and hence cannot fully account for the discrepancy of the linear stability threshold between theoretical and experimental results. When a cross-flow is present, increasing the strength of the electric field in the high-$Re$ Poiseuille flow yields a more unstable flow in both modal and nonmodal stability analyses. Even though the energy analysis and the input-output analysis both indicate that the energy growth directly related to the electric field is small, the electric effect enhances the lift-up mechanism. The symmetry of channel flow with respect to the centerline is broken due to the additional electric field acting in the wall-normal direction. As a result, the centers of the streamwise rolls are shifted towards the injector electrode, and the optimal spanwise wavenumber achieving maximum transient energy growth increases with the strength of the electric field.
We perform a three-dimensional, short-wavelength stability analysis on the numerically simulated two-dimensional flow past a circular cylinder for Reynolds numbers in the range $50le Rele300$; here, $Re = U_{infty}D/ u$ with $U_infty$, $D$ and $ u$ being the free-stream velocity, the diameter of the cylinder and the kinematic viscosity of the fluid, respectively. For a given $Re$, inviscid local stability equations from the geometric optics approach are solved on three distinct closed fluid particle trajectories (denoted as orbits 1, 2 & 3) for purely transverse perturbations. The inviscid instability on orbits 1 & 2, which are symmetric counterparts of one another, is shown to undergo bifurcations at $Reapprox50$ and $Reapprox250$. Upon incorporating finite-wavenumber, finite-Reynolds number effects to compute corrected local instability growth rates, the inviscid instability on orbits 1 & 2 is shown to be suppressed for $Relesssim262$. Orbits 1 & 2 are thus shown to exhibit a synchronous instability for $Regtrsim262$, which is remarkably close to the critical Reynolds number for the mode-B secondary instability. Further evidence for the connection between the local instability on orbits 1 & 2, and the mode-B secondary instability, is provided via a comparison of the growth rate variation with span-wise wavenumber between the local and global stability approaches. In summary, our results strongly suggest that the three-dimensional short-wavelength instability on orbits 1 & 2 is a possible mechanism for the emergence of the mode B secondary instability.
A rigorous derivation and validation for linear fluid-structure-interaction (FSI) equations for a rigid-body-motion problem is performed in an Eulerian framework. We show that the added-stiffness terms arising in the formulation of Fanion et al. (2000) vanish at the FSI interface in a first-order approximation. Several numerical tests with rigid-body motion are performed to show the validity of the derived formulation by comparing the time evolution between the linear and non-linear equations when the base flow is perturbed by identical small-amplitude perturbations. In all cases both the growth rate and angular frequency of the instability matches within $0.1%$ accuracy. The derived formulation is used to investigate the phenomenon of symmetry breaking for a rotating cylinder with an attached splitter-plate. The results show that the onset of symmetry breaking can be explained by the existence of a zero-frequency linearly unstable mode of the coupled fluid-structure-interaction system. Finally, the structural sensitivity of the least stable eigenvalue is studied for an oscillating cylinder, which is found to change significantly when the fluid and structural frequencies are close to resonance.
In this article we consider the linear stability of the two-dimensional flow induced by the linear stretching of a surface in the streamwise direction. The basic flow is a rare example of an exact analytical solution of the Navier-Stokes equations. Using results from a large Reynolds number asymptotic study and a highly accurate spectral numerical method we show that this flow is linearly unstable to disturbances in the form of Tollmien-Schlichting waves. Previous studies have shown this flow is linearly stable. However, our results show that this is only true for G{o}rtler-type disturbances.
Modeling realistic fluid and plasma flows is computationally intensive, motivating the use of reduced-order models for a variety of scientific and engineering tasks. However, it is challenging to characterize, much less guarantee, the global stability (i.e., long-time boundedness) of these models. The seminal work of Schlegel and Noack (JFM, 2015) provided a theorem outlining necessary and sufficient conditions to ensure global stability in systems with energy-preserving, quadratic nonlinearities, with the goal of evaluating the stability of projection-based models. In this work, we incorporate this theorem into modern data-driven models obtained via machine learning. First, we propose that this theorem should be a standard diagnostic for the stability of projection-based and data-driven models, examining the conditions under which it holds. Second, we illustrate how to modify the objective function in machine learning algorithms to promote globally stable models, with implications for the modeling of fluid and plasma flows. Specifically, we introduce a modified trapping SINDy algorithm based on the sparse identification of nonlinear dynamics (SINDy) method. This method enables the identification of models that, by construction, only produce bounded trajectories. The effectiveness and accuracy of this approach are demonstrated on a broad set of examples of varying model complexity and physical origin, including the vortex shedding in the wake of a circular cylinder.