skip to main content

This content will become publicly available on July 10, 2023

Title: Long-wavelength equations of motion for thin double vorticity layers
We consider the time evolution in two spatial dimensions of a double vorticity layer consisting of two contiguous, infinite material fluid strips, each with uniform but generally differing vorticity, embedded in an otherwise infinite, irrotational, inviscid incompressible fluid. The potential application is to the wake dynamics formed by two boundary layers separating from a splitter plate. A thin-layer approximation is constructed where each layer thickness, measured normal to the common centre curve, is small in comparison with the local radius of curvature of the centre curve. The three-curve equations of contour dynamics that fully describe the double-layer dynamics are expanded in the small thickness parameter. At leading order, closed nonlinear initial-value evolution equations are obtained that describe the motion of the centre curve together with the time and spatial variation of each layer thickness. In the special case where the layer vorticities are equal, these equations reduce to the single-layer equation of Moore ( Stud. Appl. Math. , vol. 58, 1978, pp. 119–140). Analysis of the linear stability of the first-order equations to small-amplitude perturbations shows Kelvin–Helmholtz instability when the far-field fluid velocities on either side of the double layer are unequal. Equal velocities define a circulation-free double vorticity layer, more » for which solution of the initial-value problem using the Laplace transform reveals a double pole in transform space leading to linear algebraic growth in general, but there is a class of interesting initial conditions with no linear growth. This is shown to agree with the long-wavelength limit of the full linearized, three-curve stability equations. « less
Authors:
; ; ;
Award ID(s):
1706934
Publication Date:
NSF-PAR ID:
10326441
Journal Name:
Journal of Fluid Mechanics
Volume:
942
ISSN:
0022-1120
Sponsoring Org:
National Science Foundation
More Like this
  1. We study the evolution of a 2D vortex layer at high Reynolds number. Vortex layer flows are characterized by intense vorticity concentrated around a curve. In addition to their intrinsic interest, vortex layers are relevant configurations because they are regularizations of vortex sheets. In this paper, we consider vortex layers whose thickness is proportional to the square-root of the viscosity. We investigate the typical roll-up process, showing that crucial phases in the initial flow evolution are the formation of stagnation points and recirculation regions. Stretching and folding characterizes the following stage of the dynamics, and we relate these events to the growth of the palinstrophy. The formation of an inner vorticity core, with vorticity intensity growing to infinity for larger Reynolds number, is the final phase of the dynamics. We display the inner core's self-similar structure, with the scale factor depending on the Reynolds number. We reveal the presence of complex singularities in the solutions of Navier–Stokes equations; these singularities approach the real axis with increasing Reynolds number. The comparison between these singularities and the Birkhoff–Rott singularity seems to suggest that vortex layers, in the limit $Re\rightarrow \infty$ , behave differently from vortex sheets.
  2. Abstract In this article, the recently discovered phenomenon of delayed Hopf bifurcations (DHB) in reaction–diffusion partial differential equations (PDEs) is analysed in the cubic Complex Ginzburg–Landau equation, as an equation in its own right, with a slowly varying parameter. We begin by using the classical asymptotic methods of stationary phase and steepest descents on the linearized PDE to show that solutions, which have approached the attracting quasi-steady state (QSS) before the Hopf bifurcation remain near that state for long times after the instantaneous Hopf bifurcation and the QSS has become repelling. In the complex time plane, the phase function of the linearized PDE has a saddle point, and the Stokes and anti-Stokes lines are central to the asymptotics. The non-linear terms are treated by applying an iterative method to the mild form of the PDE given by perturbations about the linear particular solution. This tracks the closeness of solutions near the attracting and repelling QSS in the full, non-linear PDE. Next, we show that beyond a key Stokes line through the saddle there is a curve in the space-time plane along which the particular solution of the linear PDE ceases to be exponentially small, causing the solution of the non-linearmore »PDE to diverge from the repelling QSS and exhibit large-amplitude oscillations. This curve is called the space–time buffer curve. The homogeneous solution also stops being exponentially small in a spatially dependent manner, as determined also by the initial data and time. Hence, a competition arises between these two solutions, as to which one ceases to be exponentially small first, and this competition governs spatial dependence of the DHB. We find four different cases of DHB, depending on the outcomes of the competition, and we quantify to leading order how these depend on the main system parameters, including the Hopf frequency, initial time, initial data, source terms, and diffusivity. Examples are presented for each case, with source terms that are a uni-modal function, a smooth step function, a spatially periodic function and an algebraically growing function. Also, rich spatio-temporal dynamics are observed in the post-DHB oscillations. Finally, it is shown that large-amplitude source terms can be designed so that solutions spend substantially longer times near the repelling QSS, and hence, region-specific control over the delayed onset of oscillations can be achieved.« less
  3. Abstract

    We formulate the two-dimensional gravity-capillary water wave equations in a spatially quasi-periodic setting and present a numerical study of solutions of the initial value problem. We propose a Fourier pseudo-spectral discretization of the equations of motion in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on a torus. We adopt a conformal mapping formulation and employ a quasi-periodic version of the Hilbert transform to determine the normal velocity of the free surface. Two methods of time-stepping the initial value problem are proposed, an explicit Runge–Kutta (ERK) method and an exponential time-differencing (ETD) scheme. The ETD approach makes use of the small-scale decomposition to eliminate stiffness due to surface tension. We perform a convergence study to compare the accuracy and efficiency of the methods on a traveling wave test problem. We also present an example of a periodic wave profile containing vertical tangent lines that is set in motion with a quasi-periodic velocity potential. As time evolves, each wave peak evolves differently, and only some of them overturn. Beyond water waves, we argue that spatial quasi-periodicity is a natural setting to study the dynamics of linear and nonlinear waves, offering a third option to the usual modeling assumptionmore »that solutions either evolve on a periodic domain or decay at infinity.

    « less
  4. Solovjovs, Sergejs (Ed.)
    In the present paper, we summarize the results of the research devoted to the problem of stability of the fluid flow moving in a channel with flexible walls and interacting with the walls. The walls of the vessel are subject to traveling waves. Experimental data show that the energy of the flowing fluid can be transferred and consumed by the structure (the walls), inducing “traveling wave flutter.” The problem of stability of fluid-structure interaction splits into two parts: (a) stability of fluid flow in the channel with harmonically moving walls and (b) stability of solid structure participating in the energy exchange with the flow. Stability of fluid flow, the main focus of the research, is obtained by solving the initial boundary value problem for the stream function. The main findings of the paper are the following: (i) rigorous formulation of the initial boundary problem for the stream function, ψ x , y , t , induced by the fluid-structure interaction model, which takes into account the axisymmetric pattern of the flow and “no-slip” condition near the channel walls; (ii) application of a double integral transformation (the Fourier transformation and Laplace transformation) to both the equation and boundary and initial conditions,more »which reduces the original partial differential equation to a parameter-dependent ordinary differential equation; (iii) derivation of the explicit formula for the Fourier transform of the stream function, ψ ˜ k , y , t ; (iv) evaluation of the inverse Fourier transform of ψ ˜ k , y , t and proving that reconstruction of ψ x , y , t can be obtained through a limiting process in the complex k -plane, which allows us to use the Residue theorem and represent the solution in the form of an infinite series of residues. The result of this research is an analytical solution describing blood flowing through a channel with flexible walls that are being perturbed in the form of a traveling wave.« less
  5. The effect of slip surfaces on the laminar–turbulent separatrix of plane Poiseuille flow is studied by direct numerical simulation. In laminar flows, the inclusion of the slip surfaces results in a drag reduction of over 10 %, which is in good agreement with previous studies and the theory of laminar slip flows. Turbulence lifetimes, the likelihood that turbulence is sustained, is investigated for transitional flows with various slip lengths. We show that slip surfaces decrease the likelihood of sustained turbulence compared to the no-slip case, and the likelihood is further decreased as slip length is increased. A more deterministic analysis of the effects of slip surfaces on a transition to turbulence is performed by using nonlinear travelling-wave solutions to the Navier–Stokes equations, also known as exact coherent solutions. Two solution families, dubbed P3 and P4, are used since their lower-branch solutions are embedded on the boundary of the basin of attraction of laminar and turbulent flows (Park & Graham, J. Fluid Mech. , vol. 782, 2015, pp. 430–454). Additionally, they exhibit distinct flow structures – the P3 and P4 are denoted as core mode and critical-layer mode, respectively. Distinct effects of slip surfaces on the solutions are observed by the skin-friction evolution, linearmore »growth rate and phase-space projection of transitional trajectories. The slip surface appears to modify the transition dynamics very little for the core mode, but quite considerably for the critical-layer mode. Most importantly, the slip surface promotes different transition dynamics – an early and bypass-like transition for the core mode and a delayed and H- or K-type-like transition for the critical-layer mode. We explain these distinct transition dynamics based on spatio-temporal and quadrant analyses. It is found that slip surfaces promote the prevalence of strong wall-toward motions (sweep-like events) near vortex cores close to the channel centre, inducing an early transition, while long sustained ejection events are present in the region of the $\unicode[STIX]{x1D6EC}$ -shaped vortex cores close to the critical layer, resulting in a delayed transition. This should motivate flow control strategies to fully exploit these distinct transition dynamics for transition to turbulence.« less