## A. Detailed derivations

This appendix gives some derivations mentioned earlier in the various chapters. First is the derivation of the redistribution equations (4.8) and following. We write the difference between the Fourier transform of the redistributed vorticity (4.5) and the exactly diffused vorticity (4.6) in terms of the viscous scale and the scaled relative positions (4.7):

 (A.1)

In the redistribution method, this error is made small of order , or for any finite wavenumber , by equating the Taylor series expansions of the two terms within the parentheses to that order. This produces, for any ,
 (A.2)

Expanding using the binomial theorem, the individual equations become, for ,
 (A.3)

These are the redistribution equations written out in (4.8) and following.

The remaining error in the Fourier transform, needed in chapter 5, is according to the Taylor series remainder theorem

 (A.4)

where the values of , and represent the undetermined midpoints in the remainder theorem, and is the even integer or .

To find the lower bound to the redistribution radius mentioned in section 6.2.3, we integrate (A.2) over the unit circle to produce

 (A.5)

Hence, in terms of the even integer or ,
 (A.6)

which implies the lower bound for the redistribution radius.

Up to fourth-order accuracy, this estimate for the minimum radius is precise. It may be verified by direct substitution into the redistribution equations that a positive second-order solution is obtained by spreading the fractions evenly over the circle with scaled radius . Similarly, a positive fourth-order solution is obtained by giving the vortex being redistributed a fraction and spreading the other half evenly over the circle .

Next we verify an assertion made in section 6.2.3: as long as all vortices are redistributed, the region containing the vortices must expand a finite scaled amount in each direction. To do so, we derive a lower bound to using the cases and , 1, and 2 of (A.3):

 (A.7)

The solution of the quadratic shows that
 (A.8)

Applied to the vortex at the largest value of , this value describes how much the vortex distribution needs to expand in the -direction in order for the redistribution equations to be solvable at that vortex. Since the redistribution problem is independent of the angular position of the coordinate system, this minimum expansion applies in any direction.

Next we verify an assertion made in section 6.2.3 and subsection 9.5: for third-order accuracy or higher, the scaled spacing between the vortices cannot be arbitrarily large. Defining and using (A.5) for and , it is seen that

 (A.9)

Hence the vortices cannot be spaced further apart than a scaled distance .

Finally we verify an assertion made in section 6.2.3: there are finite values and so that a positive solution to the redistribution equations exists within the circle with scaled radius provided that there are no square holes exceeding a scaled size in the distribution of the vortices. To do so, we first note that the diffusing delta function

 (A.10)

gives an exact continuous solution to our redistribution equations, replacing by . We now discretize this continuous solution using a cut-off at radius and a subdivision of the remaining domain into squares of size . We then select one vortex in each square and give it fraction
 (A.11)

where the are for polynomials of degree satisfying
 (A.12)

and zero elsewhere. In the above expression for , the first term provides a positive approximate solution to the redistribution equations, since the sums over become straightforward numerical approximations to the corresponding integrals of the continuous function . The second term gives corrections that make this approximation exact for suitable values of the constants . When is large enough and is small enough, these corrections do not change the positivity of the first term. This can be seen as follows: since the polynomials are bounded, there is a finite maximum value for the constants below which the correction terms cannot change the sign of the first term to . Further, since the redistribution equations give a system of equations for the that tends to a unit matrix, there is a value of below which the maximum can be bounded by a multiple of the maximum error due to the first term in (A.11). That error can be reduced to any finite amount by selecting a large enough and a small enough to make the numerical integrals sufficiently accurate. Hence the required positive total solution (A.11) can always be assured for some finite and .

At least for the case of first-order accuracy, , for any greater than the minimum value , a finite hole size exists that ensures a positive solution. This can be seen by selecting nine vortices to satisfy the redistribution equations. Eight of these are chosen as closely as possible to eight equally spaced points on the outside circle and given a nominal fraction , and the last point is chosen to be the vortex being redistributed, and given a nominal weight . This satisfies the redistribution equations approximately, and it is readily seen that for these nominal positions, the needed corrections in the weights to make the approximation exact can be bounded by the errors. Thus, similar as in the derivation above, the corrections do not change the sign of the weights when is small enough. The actual value of is unknown, but clearly must tend to zero when ; the allowed hole size must be small enough to ensure that there are vortices outside the circle within which no solution exists.