Critical Layer Analysis of Stuart Vortices in a Plane Jet

Asymptotic techniques are used to model quasi-steady-state vortices in the plane (Bickley) inviscid jet. A nonlinear critical layer analysis is used to find a family of steady-state 
finite amplitude two-dimensional vortices which are based on the Stuart vortex.


Introduction
A very basic flow occurring in both nature and engineering applications is the wake behind a bluff body, and it has been claimed [1] that the plane (Bickley) jet, which has a sech 2  velocity profile, provides a very good approximation to such a wake.For this velocity profile, there exist both symmetric (sinuous) and antisymmetric (varicose) instability modes [2,3] with critical layers centered on each of the two inflection points of the flow, as well as other modes with critical layers centered on the nose [4][5][6].For twodimensional perturbations, Kelly [7] pointed out that since the wavenumber of the neutral varicose mode was onehalf that of the neutral sinuous mode, it was possible to have a resonant interaction between the two modes.Some experiments have hinted that an interaction may indeed take place but have not attempted to quantify it.The authors of [8] conducted careful experiments on small deficit (turbulent) wakes and found that the development of some aspects of the flow was dependent on initial conditions, which they attributed to interactions between the varicose and sinuous modes, and other experiments, such as [9,10], have also suggested that these interactions may take place, as have preliminary numerical simulations [11].Kelly [7] explored the idea of this interaction theoretically, using an unsteady viscous critical layer and assuming a modal interaction of the Stuart-Watson type with a cubic nonlinearity, but found that there was no interaction of this type because the coefficients of the coupling terms were identically zero.However, [12] employed instead an unsteady nonlinear critical layer, in which the nonlinear effects were confined to a thin layer in which the nonlinear and unsteady terms were of the same order of magnitude, and found that using this approach there was indeed a resonant interaction between the modes.The authors of [12] used an  1/2 critical layer, meaning that the width of the critical layer was O( 1/2 ), where  denoted the order of magnitude of the instability wave amplitude, and such a critical layer is standard for both steady [13,14] and unsteady [15,16] two-dimensional critical layer analyses and is also the scaling employed in the present study.For threedimensional disturbances, Mallier [17,18] found that there was also an interaction between the varicose and sinuous modes.
As with any shear layer, when an instability mode is superimposed on a plane jet, the shear layer rolls up and a quasi-equilibrium state emerges in the form of large vortices.What is distinct about the plane jet is that two rows of vortices develop rather than just one, with a row centered on each of the critical layers.These large scale vortices are intriguing for a number of reasons, not least because it appears that an ordered state emerges from disorder, and the aim of the present study is to use asymptotic techniques to model them.Because the two instability modes can interact, our analysis must include both modes.

Mathematical Problems in Engineering
For a mixing layer with a single row of vortices, a single critical layer, and a single instability mode, two alternative models emerged in the 1960s, Stuart vortices [19] and Benney-Bergeron-Davis (BBD) vortices [13,14,[20][21][22][23][24].Stuart found an exact nonlinear steady-state solution of the inviscid Euler equations, for which the stream function can be cast in the form  = ln[cosh  −  cos ], where  is a small parameter.This solution represents an infinite row of corotating vortices, and when  = ±1, the vortices become point vortices, but when  = 0, we recover unidirectional shear flow.In the form above, this stream function obeys ∇ 2  = (1 −  2 ) −2 , which is a form of Liouville's equation.Since the vorticity is given by  = −∇ 2 , the vorticity distribution for the Stuart vortex is The stability of Stuart vortices has been studied both numerically and analytically [25][26][27][28].By contrast, the BBD vortices, which marked the genesis of nonlinear critical layer theory, were found using matched asymptotic expansions, posing an outer solution away from the critical layer and matching it to an inner solution close to the critical layer.Recently, Caillol and coworkers have used an analysis similar to the BBD theory to model steady multipolar planar vortices [29] and also the intensification or weakening of a rapidly rotating vortex of a hurricane or tornado type by considering the critical-layer-like interaction between a helical neutral vortex Rossby wave and an azimuthally mean vortex [29][30][31].For time-independent inviscid two-dimensional flow, it can be shown that the vorticity must be a function of the stream function, so that with Liouville's equation being a well-known specific example and the sinh-Poisson equation [32] another.It has been shown [33] that there is no single-valued solution to (2) for the plane jet with critical layers centered on the inflection points, so we must seek a series solution.While the vorticity for the Stuart vortex was given by (1), the BBD theory instead invoked the Prandtl-Batchelor theorem [34], which says that, for steady flow at high Reynolds numbers, the vorticity distribution inside closed streamlines must be uniform.Because of this, BBD argued that the vorticity should be uniform inside the vortex cores, so that the function F was taken to be constant there.
One early problem with the BBD theory was that the vorticity was discontinuous across the edges of the cat's eyes, which BBD claimed would be smoothed out by a viscous layer at the edge of the cat's eyes.However, the analysis at the corners of the cat's eye led to a Wiener-Hopf problem, and Brown and Stewartson [35] later showed that this Wiener-Hopf problem did not have a solution with the correct behavior, meaning that the discontinuity could not be removed.Haberman [23,24] showed that the problems with the BBD solutions could be removed if certain mean flow distortions were added to the problem, and we include those same mean flow corrections in our analysis here as series about the critical layers.In Haberman's solutions, the velocity and vorticity were continuous, but gradients of vorticity were not, and the same will hold true here, although it is possible to smooth out the discontinuity in the gradient of vorticity by utilizing a thin viscous layer at the edge of the vortex cores.
A second problem with the BBD solutions was also pointed out by Stewartson [16,36].In order for vorticity to homogenize inside closed streamlines, as it is assumed to do in the BBD solution where the Prandtl-Batchelor theorem was invoked, some viscosity must be present.Stewartson argued that the cat's eyes would decay before the vorticity had become fully homogenized, and indeed it is known that such a situation occurs for the Blasius boundary layer [37,38].However, since structures resembling the BBD solutions were reported in simulations such as those of [39], it was widely assumed that the BBD solutions must be correct.Klaassen and Peltier [40] later pointed out that the cat's eyes structures observed in simulations also resembled Stuart vortices, and this can be seen upon close examination of the original paper of [39].Although the vorticity contours in [39] resemble both the BBD solutions and Stuart vortices, the authors describe the "topography" of the vorticity field as consisting of a high rounded hill (the core) along with high ridges (the braids), with the core and braids separated by low passes.It is this hill-like structure of the vorticity inside the core, which the authors of [39] attribute to the coalescence of the vorticity into relatively compact cores, that suggests that the flow inside the core is the Stuart vortex rather than the BBD solutions, for which the vorticity field would be a plateau rather than a round hill.Similar results can be found in the simulations of [41] that calculated the quasi-equilibrium states in order to study their stability.In their Figure 3, they present surface plots of the vorticity field.As in [39], the vorticity field consists of a smooth round hill (the core) along with high ridges (the braids), with the two separated by a trough, which they attribute to the entrainment of irrotational fluid into the core.
In later simulations, Sommeria et al. [42] confirmed that the vorticity field in the quasi-equilibrium states resembled the Stuart vortex rather than the homogenized cores of the BBD solutions and suggested that the Stuart vortex was a preferred state on entropy grounds.In their Figures 2-4, they show the rollup of vorticity.Because their contours are labelled, it is possible to discern from their simulations that after the climax, the vorticity decays to a smooth state resembling a round hill.In their Figures 8-10, the authors of [42] plot a scatter-plot of the vorticity against the stream function for the steady state.If the vorticity distribution was exactly that of the Stuart vortex, it would obey Liouville's equation, and a plot of the log of the vorticity against the stream function should be a straight line.The authors of [42] claim that, at high Reynolds numbers, inside the core there is excellent agreement between their scatter-plots and a straight line, so that the vorticity distribution in the equilibrium states is the Stuart vortex.
Finally, Wang and Maxey [43] gave a kinematical description of the mixing process.They noted that, in the quasisteady state, the vorticity field develops into a well-defined structure similar to the Stuart vortex.In their Figure 1, they plot contours of vorticity, and their use of color in this figure enables the structure of the vorticity field to be clearly seen, with once again the vorticity developing into a round hill surrounded by the ridge-like braids.Wang and Maxey attributed this to the vorticity layer being stretched and wrapped into a spiral structure and then diffusion acting to smooth the vorticity in between the spirals, and they show a cartoon of this process in their Figure 4.
The simulations discussed above suggest that, rather than the vorticity homogenizing inside the vortex cores, an alternative scenario, presented by Rhines and Young [44] in a different context, occurs.Rhines and Young argued that homogenization of vorticity within closed streamlines will take place in two stages: a rapid phase dominated by shearaugmented diffusion in which the vorticity will tend to a smooth but not necessarily uniform distribution and a second much slower phase in which the process of homogenization will be completed.The simulations presented in [42] suggest that the first, rapid, phase postulated by Rhines and Young will take place during vortex rollup in a shear layer but that, as suggested by Stewartson, the vortices will decay before the second stage has occurred, so that the second stage would only occur if forcing were present.Thus, the vorticity inside the vortex cores will have a smooth but nonuniform distribution, and we believe, on entropy grounds, that the Stuart vortex is that distribution.These same arguments apply to the vorticity distribution inside the vortex cores of the plane jet considered here.
In our study, we have modeled a quasi-steady state with a steady state, and, as pointed out by a referee, a steady state may not in fact be achievable in practice.The viscosity needed to smooth the vorticity distribution inside the core would ultimately cause the cat's eye to decay, and we are modeling the quasi-steady state which develops before this decay sets in.The unattainability of a steady state is a longstanding concern about vortex solutions; for example, [45][46][47] would seem to indicate that the flow in the unsteady twodimensional critical layer analyses of [15,16] would never reach a steady state.
Caillol and Grimshaw have considered the unattainability of a steady state in some detail, and they invoked the Prandtl-Batchelor theorem [34] and the work of Rhines and Young [44] and extended it for the cases they studied.Caillol and Grimshaw considered several quasi-steady regimes with critical layers, including an internal solitary wave with a weakly stratified critical layer [48], a Rossby solitary wave [49], a Rossby elevation wave where no induced mean flow is created [50], a finite-period neutral shear wave in a vortex [30], and multipolar planar vortices [31].In their studies, Caillol and Grimshaw drew on earlier work by Troitskaya and collaborators in Nizhny Novgorod on the subtle role of viscosity in nonlinear critical layers [51,52].
In the present study, we will use a time-independent nonlinear critical layer analysis to model the vortices which develop during the rollup of a plane jet, using asymptotic techniques due to [13,14] but with a vorticity distribution inside the cores which is smooth but nonconstant, resembling Stuart vortices [19] rather than BBD vortices, and in that sense our solution may be thought of as Stuart vortices in a plane jet.This sort of approach has been used previously both by ourselves for rotating and stratified mixing layers [53][54][55] and by others [33,56].By contrast, compressible Stuart vortices were found numerically [57,58].The difference between the plane jet discussed here and the mixing layer is that the plane jet has two critical layers, with a separate expansion required for each of these critical layers, and we must include both the varicose and sinuous modes in these expansions.
The rest of the paper is as follows.In Section 2, we present our analysis and find the steady-state vortex solutions.Essentially, our analysis involves posing an outer expansion in the main body of the fluid away from the critical layer, at which the base velocity is equal to the phase velocity of the disturbance, and an inner expansion near to the critical layer, using rescaled variables, and matching the two expansions together.For the outer expansion, we use the inviscid incompressible Euler equations in the stream function formulation.We suppose that we have superimposed a small, spatially periodic disturbance on the base flow and expand in powers of the amplitude; this expansion becomes disordered at the critical layer where it is necessary to pose an inner expansion using stretched coordinates.In this inner solution, rather than assuming that vorticity is uniform (homogenized) inside the closed cat's eyes (as was done in the BBD theory and its extensions by Maslowe and Haberman), we use the Stuart vortex as our guide, for the reasons outlined above.Finally, in Section 3, we will make some concluding remarks.

Formulation and Outer Expansion.
The motion of an inviscid incompressible fluid in two dimensions is governed by the Euler equations in the vorticity-stream function formulation: where (, , ) is the stream function, with the velocity given by  = (  , −  , 0).The objective of the present study is to seek nonlinear waves of permanent form which are periodic in  = (−), where  is the wavenumber and  is the phase velocity.If we work in a frame traveling in the streamwise direction at the phase velocity, we can set  = 0 without affecting our analysis, so that  =  and our solutions are independent of time, and we can replace (3) by which has a solution of the functional form To find our solution, we will pose a weakly nonlinear expansion of the form  =  (0) +  [ (1)  1   +  (1)  2  2 + cc] +  2 [ (2)  0 + ( (2)  1   +  (2)  2  2 + (2)  3  3 +  (2)  4  4 + cc)] where cc means complex conjugate and  ∼ ( 1) is a scaling parameter, while  (0) = tanh  − 2/3 is the base flow which is related to the base velocity  (0) = sech 2 −2/3 by  (0) () =  (0) ().The "2/3" here is the well-known value of the phase velocity for singular neutral modes which have critical layers at the inflection points of the flow, that is, at  = ±  , where   = arccosh√3/2 with  (0) (±  ) =  (0) (±  ) = 0.
There are two such modes, a varicose mode with a streamwise wavenumber of 1 and a sinuous mode with a streamwise wavenumber of 2 [2,3].It has been shown [33] that there is no single-valued solution to (5) for the Bickley jet with this phase velocity; that is why we must seek a series solution.We will also expand the wavenumber  = 1 +  1 + ⋅ ⋅ ⋅ to allow for a small amount of detuning.By substituting expansion (6) into governing equation ( 4), we get a series of equations at successive powers of , and to abbreviate these equations, we introduce the operator ).At leading order, O(), we assume that the disturbance is comprised of two modes, which obey L 1  (1) 1 = 0 and L 2  (1) 2 = 0.There are well-known closed form solutions  (1)   1 =  1 sech tanh  and  (1)

2
=  2 sech 2  for these two modes [2,3], where the constants  1 and  2 are a measure of the relative magnitudes of the wavenumbers 1 and 2 modes.We need to expand these solutions as series about the critical layers at  = ±  , using the variables  ± =  ∓   , as follows: where we take the upper and lower elements in ± and ∓ for the expansion about  =   and  = −  , respectively.At the next order, O( 2 ), we need to include a mean flow distortion term,  (2)  0 ∼  (2)± 0 +  (2)± 0  ± + ( (2)± 0  2 ± /2!) + ⋅ ⋅ ⋅ , and we need to calculate four modes: which we can solve as series about the critical layers: where the 's and 's are constants which can be evaluated by enforcing the condition that the 's vanish as  → ±∞.
At O( 3 ), once again we need a mean flow distortion term, which at this order is singular at the critical layers,  (3)   0 2 ± /2!) + ⋅ ⋅ ⋅ , and we need to calculate six modes: which we can again solve as series: These solutions and their derivatives are singular in the limit  ± → 0, with these discontinuities being caused by the expansion becoming disordered near the critical layers, and to remedy this, it is necessary to use stretched variables in the critical layers.

Inner Expansion.
To remedy the singularities observed in the outer expansion, we will follow [13,14] and introduce the stretched variables  + =  −   =  1/2  + ,  − =  +   =  1/2  − , and  = Ψ in the critical layers.This scaling is chosen so that nonlinear terms enter at the same order as the singularities inside the critical layer.Once again, we will pose an expansion in , the form of which is suggested by the outer solution rewritten in the inner variables as follows: with the  3/2 (ln  1/2 ) terms caused by the logs in the outer solution.We will have slightly different solutions for the expansions inside the upper and lower critical layers.In what follows, we will use the superscripts "+" and "−" to denote the solutions at the upper and lower critical layers, respectively, and we will drop the superscript when the two expressions are the same.
At O( −1 ), we have a constant term, At O( 0 ), our solution must obey which has a general solution Ψ (0)±  ±  ± = F (0)± (Ψ (0)± ).In order to solve (13) and subsequent equations, we will use the outer solution as our guide.At this order, the outer solution written in inner variables becomes and since this is nonsingular and satisfies ( 13), we will use it as our inner solution.Solution ( 14) is of the same form as the leading order terms of both the Stuart vortex stream function [19] and the BBD solution but with of course two waves present rather than just one.Using ( 14), the O( 1/2 ) inner equation is and once again, our outer solution written in inner variables satisfies this and is nonsingular, so we may use it as our solution: Using ( 14) and ( 16), the O() inner equation is and once again, our outer solution written in inner variables satisfies this and is nonsingular, so we may use it as our solution as follows: At the first logarithmic order O( 3/2 ln( 1/2 )), the inner equation is and we may again use the outer expansion to construct our inner solution as follows:

Mathematical Problems in Engineering
At O( 3/2 ), the inner equation is It is at this order that the singularities appear in the outer solution written in inner variables, and as  ± → ±∞, we have To proceed, we must solve for the flow inside and outside the core separately.Inside the core, we can assume that the solution is regular and take Ψ (3) = Ψ (3) with which satisfies (21), while outside the core we must take which has a general solution Ψ (3)±  ±  ± = F (3)± (Ψ (0)± ).From the outer expansion, we can deduce where  = −1 above the critical layer and +1 below it.To find the functional form of F (3) , we need to follow [13,14] and reintroduce weak viscosity, which will allow us to obtain a viscous secularity condition.The viscous counterpart of ( 4) is where Re is the Reynolds number.To ensure that viscous terms enter into the critical layer at the correct order, we will write ( Re) −1 =  3/2 , where  is the so-called Benney-Bergeron parameter [13,14] measuring the relative importance of the viscous and nonlinear terms inside the critical layer.In our analysis, we will assume that viscosity is weak even inside the critical layer with  ≪ 1.If we denote the inner viscous solution by Ψ (V) , we may expand this as a series in : where the first term in the expansion is the inviscid solution and the second term is the viscous correction.We may also expand this viscous correction Ψ as a series in , just as we posed an expansion (12) of the inviscid solution Ψ.To determine the function F (3) , we need to arrive at an equation for the viscous correction to Ψ (3) , but to do this, we first need to calculate the viscous correction at lower powers of .At the first few orders, we find and the O( 0 ), O( 1/2 ), and O( 1 ) viscous corrections all vanish, Ψ(0) = Ψ(1) = Ψ(2) = 0.At O( 3/2 ), we find that Ψ(3) obeys If we use  =  and Ψ (0) as von Mises coordinates, we can write (29) in the form which can be integrated with respect to  to give If we integrate (30) between  = 0 and  = 2 and use periodicity, we are left with where  is a constant.If we integrate this with respect to Ψ (0) , we obtain where  is a second constant.This equation ( 33) holds for all values of Ψ (0) , and in particular, we can evaluate it for large values of Ψ (0) , which correspond to large values of ||.In this limit, we know F (3) from ( 25), while Using ( 34) and ( 25) in (33), we find that  ± = ±(32/ √ 3), while  = 0. Using these, we can write Comparing our two expressions for F (3) (Ψ (0) ) also tells us that  (3)± 0,−1 = ∓2 2 1  2 /3 1/2 , so that the singular part of the O( 3/2 ) mean flow distortion is only present when both waves are present.

We now need an expression for
.It is straightforward to show that we can write This expression simplifies if only one mode is present.When  2 = 0, it reduces to where  is the elliptic integral of the second kind [59], but when  1 = 0, ) .
Our results are shown graphically in Figure 1 where we plot the inner stream function.In Figures 1(a) and 1(b), we show the flow in the upper and lower critical layers, respectively, when only  1 is present, and the two critical layers are in phase with one vortex present in the box 0 ≤  ≤ 2.In Figures 1(c) and 1(d), we show the flow when only  2 is present, and the two critical layers are out of phase with two vortices present in 0 ≤  ≤ 2.In each of Figures 1(a vortices.In the remaining figures, both modes are present. In Figures 1(e) and 1(f),  1 ≫  2 and the flow resembles that in Figures 1(a) and 1(b).In Figures 1(g) and 1(h),  2 ≫  1 and the flow resembles that in Figures 1(c) and 1(d).In Figures 1(i) and 1(j), the two modes are of comparable size and the two critical layers differ from each other slightly.In the upper layer, one vortex is larger than the other, while in the lower layer, the two vortices are of comparable size but are starting to merge at the edge of the box.

Discussion
In the previous section, we have used nonlinear critical layer techniques to find nonlinear quasi-steady-state vortex solutions for the plane jet; these solutions represent the quasiequilibrium large scale vortices rings that emerge following vortex rollup in a plane jet.These solutions have a smooth but nonuniform vorticity distribution inside the cores, so they may be thought of as Stuart vortices in a plane jet, rather than BBD-style vorticity-homogenized cores [13,14].Our reasons for using the Stuart vortex as the solution inside the core were presented in the introduction and are essentially due to a combination of the Brown and Stewartson argument that the cat's eyes would decay before vorticity has had time to homogenize [35], which has been confirmed by numerical simulations and the claim by Sommeria et al. [42] that the Stuart vortex is a preferred state on equilibrium grounds.