Generalized HelmholtzKirchhoff model for two dimensional distributed vortex motion
Abstract
The twodimensional NavierStokes equations are rewritten as a system of coupled nonlinear ordinary differential equations. These equations describe the evolution of the moments of an expansion of the vorticity with respect to Hermite functions and of the centers of vorticity concentrations. We prove the convergence of this expansion and show that in the zero viscosity and zero core size limit we formally recover the HelmholtzKirchhoff model for the evolution of pointvortices. The present expansion systematically incorporates the effects of both viscosity and finite vortex core size. We also show that a loworder truncation of our expansion leads to the representation of the flow as a system of interacting Gaussian (i.e. Oseen) vortices which previous experimental work has shown to be an accurate approximation to many important physical flows [9].
1 Introduction
In this paper we represent solutions of the twodimensional NavierStokes equations as a system of interacting vortices. This expansion, which generalizes the HelmholtzKirchhoff model of interacting point vortices in an inviscid fluid, systematically incorporates the effects of both vorticity and finite vortex core size. Furthermore, we give conditions which guarantee the convergence of our expansion. Incompressible viscous flow has two standard analytic representations: a formulation in terms of the primitive velocity and pressure variables, and a formulation in terms of the velocity and vorticity variables [7]. The velocityvorticity representation has particular advantages when boundaries are unimportant, since vorticity cannot be created or destroyed in the interior of a fluid. The vorticity field can also be directly related to physically observed flow structures such as line and ring vortices.
In two space dimensions, the vorticity field has the additional advantage of reducing to a scalar. An early representation of twodimensional flow in terms of moving point vortices was developed by HelmholtzKirchhoff [4] and by Helmholtz [15] . The point vortex model has been studied extensively  a thorough review of the model and recent developments are described in the monograph by Newton [11]. While the HelmholtzKirchhoff point vortex model captures many of the basic physical phenomena observed in twodimensional rotational flows, experiments with even relatively simple vortex configurations exhibit complications far beyond the point vortex predictions [9]. Additionally, the classical point vortex model neglects the effects of viscosity. However, these experiments also reinforce the idea that in many circumstances the fluid flow may be well approximated by a collection of interacting vortices  albeit vortices with finite core size, subject to the effects of viscosity. A few recent studies of the interaction of viscous vortices can be found in [1, 5, 6, 8, 9, 14]. The main focus of these papers is the merger of two like signed vortices. In [8] the authors use a spatial moment model for 2D Euler equations which later incorporates weak Newtonian viscosity to derive equations of motions for two like signed vortices. A metastable state is found before merger which consists of two rotating, nearcircular, vortices. More recently, it has been conjectured that in a twovortex system the profiles relax to a pair of gaussian vortices before merging [6] . Thus, it is of interest to extend and generalize the HelmholtzKirchhoff point vortex model to a model that incorporates nonzero vortex core size and viscous effects while retaining its basic form. Such an extension is the goal of this paper.
The governing equations for the velocity and pressure variables are
(1) 
(2) 
where is the fluid density and is the kinematic viscosity. Taking the curl of (1) and (2) gives
(3) 
(4) 
which are the governing equations for the velocityvorticity variables. For twodimensional flows, the vorticity vector is perpendicular to the plane of the flow, and the third term on the lefthand side of (3) vanishes. The condition is identically satisfied, and equation (3) then reduces to the single scalar equation:
(5) 
where is the single, nonzero component of the vorticity. A drawback of the formulation in (5) is that the velocity of the fluid is still present in the equation. However, assuming that the vorticity field is sufficiently localized, the velocity vector can be computed in terms of the vorticity by the BiotSavart law
(6) 
where for a twovector , .
In this paper we use equations (5) and (6) to develop a vorticity representation of twodimensional viscous flow. Our representation is based on a decomposition of the vorticity field into a set of moving distributed vortices. Differential equations are derived for the motion of the vortex centers and for the time evolution of the vortex distributions. The evolution of each individual vortex is represented as an expansion with respect to a sequence of Hermite functions. Such expansions have proven useful in theoretical studies of twodimensional fluid flows ([2], [3]) and the leading order term in this expansion is precisely the Gaussian vortex (i.e. Oseen vortex [13],[12]) whose utility as an approximation for vortex interaction was shown in [9]. We show that the coefficients in this expansion satisfy a system of ordinary differential equations whose coefficients can be explicitly represented in terms of a fixed, computable kernel function. We also prove the convergence of this expansion. It is shown that our representation reduces in the appropriate limit to the HelmholtzKirchhoff model, and allows at the same time arbitrarily complex evolution and interaction of the moving vortices.
In the present paper we concentrate on the mathematical formulation of the generalized HelmholtzKirchhoff model. In future work we will explore the predictions of this model both numerically and analytically in a number of different physical settings.
2 The “multivortex” expansion
In this section we separate the solution of the vorticity equation into components and derive separate evolution equations for each component. From a physical point of view this decomposition will be most useful when each of the components corresponds to a localized region of vorticity (e.g. a vortex) well separated from the other lumps but the mathematical development described below is well defined without regard to these physical considerations. However, with this application in mind, we will often refer to each of the components as a “vortex”.
Consider the initial value problem for the twodimensional vorticity equation
(7)  
where is the velocity field associated to the vorticity field . We begin by decomposing the initial vorticity distribution by writing
(8) 
Of course this decomposition is not unique  even the number of pieces, , into which we decompose the vorticity is up to us to choose. In general the choice we make will be motivated by physical considerations, however, for the development below, all we require of the decomposition is that the total vorticity of each vortex is nonzero, i.e.
(9) 
If (9) is satisfied we define by
(10) 
or equivalently
(11) 
We now write the vorticity for as
(12) 
and the velocity field as
(13) 
where is the velocity field associated to by the BiotSavart Law. Of course we still have to define the equations of motion for and .
The centers of the vorticity regions, , and the vorticity regions themselves evolve via a coupled system of ordinarypartial differential equations constructed so that in the limit of zero viscosity and when the different components of the vorticity happen to be point vortices (i.e. Diracdelta functions) we recover the HelmholtzKirchhoff point vortex equations. If we take the partial derivative of (12) and use the equation satisfied by the vorticity, we find:
Given this equation it is natural to define as the solution of the equation:
(15)  
To close this system of equations we must specify how the centers of vorticity evolve. We impose the condition that the first moment of each vorticity region must vanish at every time , i.e. we require that
(16) 
(Note that this equation really contains two conditions  one for each component of .) We impose this condition to fix the evolution of because Gallay and Wayne have recently shown [2] that if one considers the evolution of general solutions of (2) the solution will approach an Oseen vortex, and the rate of the approach will be faster if the vorticity distribution has first moment equal to zero. Solutions of (2) preserve the first moment, and hence if the initial conditions have first moment equal to zero the solution will have first moment zero for all time. The equations (15) no longer preserve the first moment and thus we impose this condition for all time, which then defines the motion of the center of vorticity.
Note that if we change variables in (16) to , we find
(17) 
Since this equation holds for all we can differentiate both sides with respect to to obtain
(18) 
Using (15) we can insert the formula for into this integral and we obtain:
(19) 
We first note that if we integrate twice by parts, we have
(20) 
Next if we take the component and integrate by parts we find
(21) 
where and .
Remark 2.1.
Finally, recalling that the velocity field is incompressible we can rewrite
so, again considering the component and integrating by parts we have
(22)  
Thus, if we combine (20) (21), and (22) we see that (2) reduces to the system of ordinary differential equations for the centers of the vorticity distributions:
(23) 
supplemented by the initial conditions (10), while the components of the vorticity evolve according to the partial differential equations (15) with initial conditions
(24) 
Remark 2.2.
Consider (23) in the limit in which the components are all point vortices, i.e. , with the Diracdelta function. Recall that the velocity field associated with such a point vortex is
where is the antisymmetric tensor with two indices. Then if we ignore the (singular) term with in the sum on the right hand side of (23) we find
(25) 
These of course are just HelmholtzKirchhoff equations for the inviscid motion of a system of point vortices. Thus, our expansion can be regarded as a generalization of the this approximation which allows for both nonzero viscosity and vortices of finite size. To justify omitting the term with on the right hand side of (25) we note that if we approximate the delta function with a narrow, Gaussian vorticity distribution, and by the corresponding velocity field, this term will vanish by symmetry.
3 The moment expansion; case of a single center
In this section we introduce another idea  an expansion of the vorticity in terms of Hermite functions. Then, in the next section we will combine the Hermite expansion with the multivortex expansion of the previous section.
The moment expansion is an expansion of the solution of the vorticity equation in terms of Hermite functions. Define
(26) 
where . Three simple facts that we will use repeatedly are



Finally, and crucially for what follows, the vorticity function is an exact solution (called the Oseen, or Lamb, vortex) of the two dimensional vorticity equation for all values of .
Note that we will often supress the dependence of on when there is no fear of confusion.
We now define the Hermite functions of order by
(27) 
and the corresponding moment expansion of a function by
(28) 
Note that if the function in (28) is the vorticity field of some fluid the linearity of the BiotSavart law implies that we can expand the associated velocity field as:
(29) 
where
(30) 
and is the velocity field associated with the Gaussian vorticity distribution  explicitly we have:
(31) 
We define the Hermite polynomials via their generating function:
(32) 
Note that the “standard” Hermite polynomials correspond to taking .
Then using the standard orthogonality relationship for the Hermite polynomials:
(33) 
we see that the coefficients in the expansion (28) are defined by the projection operators:
(34) 
3.1 Convergence of the moment expansion
In this subsection we derive a criterion for the convergence of the moment expansion derived above and we show that if this criterion is satisfied for then it is satisfied for all subsequent times .
Our convergence criterion is based on the observation that the Hermite functions are, for any value of , the eigenfunctions of the linear operator
(35) 
This fact can be verified by direct computation and is related to the fact that can be transformed into the Hamiltonian quantum mechanical harmonic oscillator.
The Gaussian function plays a crucial role in the convergence proof, and its dependence on the parameter is particularly important in this discussion, so for this subsection only we will define
to empasize this dependence.
If one now proceeds as in Lemma 4.7 of [3] one can prove
Proposition 3.1.
The operator is selfadjoint in the Hilbert space
with innerproduct .
An immediate corollary of this proposition and the general theory of selfadjoint operators is
Corollary 3.2.
The eigenfunctions of form a complete orthogonal set in the Hilbert space .
and as a corollary of this result and the observation that the eigenfunctions of are precisely our Hermite functions , we have finally
Proposition 3.3.
Suppose that
then the expansion
converges with respect to the norm on the Hilbert space .
Thus, the following criterion guarantees that the expansion (28) for the vorticity converges:
(36) 
The main result of this subsection is the following theorem which proves that if our initial vorticity distribution satisfies (36) for some , then the solution of the vorticity equation with that initial condition will satisfy (36) for all time with and hence as a corollary if the initial vorticity distribution satisfies (36) then our moment expansion converges for all times .
Theorem 3.4.
Define the weighted enstrophy function
If the initial vorticity distribution is such that for some , and is bounded (in the norm) then is finite for all times .
Proof: The idea of the proof is to derive a differential inequality for which guarantees that if is finite then will be finite for all . Differentiating we obtain
(37)  
(38)  
We now consider the last term in (38). First note that upon integration by parts we have:
(39) 
The right hand side of the expression in (39) can again be broken up into two pieces and the second can be bounded by:
(40) 
Finally, we bound the last term in (37), which comes from the nonlinear term in the vorticity equation. In this estimate we use the fact (see [2], Lemma 2.1) that the norm of the velocity field can be bounded by a constant times the sum of the and norms of the vorticity field  i.e by . This observation, combined with the fact that , which is a consequence of the maximum principle, implies that
and hence that
(41)  
3.2 Differential equations for the moments
Assuming that the function is a solution of (2), we can derive differential equations satisfied by the moments in (28). Surprisingly the expressions for the coefficients in these expansions are quite simple and explicit.
If we differentiate (28) and assume that is a solution of the twodimensional vorticity equation we obtain:
(43)  
From the first of the “simple facts” we stated about , we see that the last term on the first line cancels the middle line and hence if we apply the projection operators defined in (34), we are left with the system of ordinary differential equations for the moments
The somewhat surprising fact which, in our opinion, makes the preceding straightforward calculations interesting is that the projection on the right hand side of (3.2) can be computed explicitly in terms of the derivatives of an relatively simple function.
We now explain how this is done. First recall that:
(45) 
In order to avoid confusing the two sets of derivatives we will rewrite these formulas as
(46)  
Inserting these formulas into the right hand side of (3.2) and using the formula for the projection operator in terms of the integration against a Hermite polynomial we obtain
(47)  
The last equality in this expression results from rewriting the Hermite polynomial in terms of its generating function.
The last step in deriving the equations for the moments is to evaluate the integral in the last line of (47). The key step in this evaluation is to recall that for these incompressible flows the velocity field can be written in terms of the derivatives of the stream function and that the Laplacian of the stream function is minus the vorticity. Thus, we can write:
(48) 
where .
Inserting this into the integral in (47) we find
(49)  
Now note that all three factors in the integrand are Gaussians and thus the integral can be evaluated explicitly, and we find
(50)  
We next compute the expression
(51)  
Recall that given a vorticity field , is the associated stream function and the velocity field associated with . Since the inverse Laplacian and derivatives in (51) act only on the dependent parts of the expression we need to evaluate
(52)  
But
is just the velocity field associated with a Gaussian vorticity distribution (i.e. an Oseen vortex) centered at the point which we know explicitly. Hence, the expression on the right hand side of (49) has the explicit representation:
(53)  
If we now return to (47) we see that in order to compute the coefficients in the moment equations we need to evaluate the divergence of this last expression with respect to which gives