The paper by Uchiyama and Fukase (1) proposes a three-dimensional vortex method for particulate flows, which is applied to a circular jet of air laden with glass particles. The numerical method presented has the following characteristics: (1) it includes the two-way coupling of the particles and the fluid, (2) it discretizes the vorticity in the gas phase in the usual manner when vortex blob methods are used, and (3) it accounts for the effects of viscosity in the fluid by means of the less-than-usual core spreading method.

The method presented in the paper under discussion is claimed to be a 3D extension of the two-dimensional vortex method for particle-laden flows proposed by Uchiyama and Naruse (2). Upon reading this previous work, however, one notes important differences in the core spreading viscous scheme. There is an apparent inconsistency in the paper in relation to the implementation of the vortex method, which will be discussed below. There are three issues that should be made clear.

- 1
The core spreading formula presented is not the same as used in the previous work by Uchiyama and Naruse (2) (and several subsequent works by Uchiyama and others). The extension to 3D of (1) refers to the particle-fluid coupling only.

- 2
The core spreading scheme in (1) corresponds to that presented by Leonard (3); that scheme is derived for simple Gaussian vortex blobs, whereas in (1) a

*cubic*Gaussian is used. This is a severe inconsistency, and readers should be made aware that core spreading is derived for a*specific*blob distribution. - 3
The core spreading method used in (1) (Leonard’s) is a two-dimensional scheme, derived from an exact solution to the 2D Navier-Stokes equation. It is unclear how the authors are applying the scheme for 3D flows; in fact, the mere application of this scheme to 3D is questionable.

The 2D antecessor (2) of the method presented in (1) incorporated the two-way coupling for particles and fluid, by means of turning the force of particles acting on the fluid into a source term in the vorticity transport equation. This technique was well known at the time. The generation of vorticity due to the presence of the particles and by the action of the drag force was pointed out in the review by Crowe et al. (4). The concept was implemented successfully in a vortex-based method in (5).

The salient feature of the 2D method of Uchiyama and Naruse (2) is the use of the core spreading method for viscous effects in the gas phase. The core spreading method has not been used often, due to it falling out of favor with vortex methods workers when consistency objections were raised (6). The scheme was presented by Leonard (3), and consists of the idea that elementary vortex blobs can be made to grow in time to simulate viscosity.

The consistency issues with the core spreading vortex method stem from the fact that the vortex blobs are advected without deformation, and thus they cannot be allowed to grow larger and larger in time. Basically, the size of the blobs represents the smallest scales that the vortex blob method can resolve. If the size of the vortex blobs is $\sigma $, then the truncation error of the vortex method grows with $\sigma 2$, as proved in (3). Moreover, since the vorticity is advected with an average velocity and not with the actual local velocity, the convection of vorticity is incorrect even in the limit of infinitely many particles, as proved by (6). These issues were successfully addressed by Rossi (7), who proposed a vortex splitting technique to control the core sizes, and proved it to be convergent. This result prompted new recent work with the core spreading method (8,9). Most recently, a new formulation for the vortex method with core spreading has been developed, in which the core size control necessary for consistency and convergence is provided (without splitting) in a spatial adaption scheme using radial basis function interpolation (10).

In the paper under discussion, and the previous work of Uchiyama and co-workers, starting with (2), no correction is provided for the consistency problem of the core spreading method. As far as can be inferred from their publications, they allow the vortex blobs to grow uncontrolled throughout a time marching calculation. This is acceptable if the Reynolds number is large and the simulations short, but it is important to underscore that the algorithm is inconsistent if applied to another flow situation.

The cited work (11) contains a rather scant derivation, and is easy to question. They claim to obtain the above formula for the change in blob radius due to viscosity in three dimensions by *approximating it by two-dimensional viscous diffusion “for convenience.”* They seem to take the 3D vortex blobs, with spherical symmetry, and approximate these with a cylinder of equal volume, then use 2D core spreading to diffuse the (finite-length!) cylindrical vortex. Apparently they then switch back to a spherical blob of volume equal to the “diffused” cylindrical one.

*algebraic*cutoff, given by:

*different cutoff function*(Chorin’s). This is disconcerting, as the core spreading method is explicitly tied to the cutoff function. In subsequent papers (12,13,14), the authors continue to give the same core spreading formula as above, and continue to use Chorin’s cutoff function. This cutoff was introduced in (15) and is the following:

*3D method*, where an algebraic cutoff is used, and some transformations between spherical and cylindrical symmetry are applied to obtain the formula given by 1.

*exactly*by each blob when making $\sigma 2$ grow linearly according to 5. Since the diffusion equation is linear, the discretized vorticity field in 2D can be seen as a superposition of spreading line vortices of different strengths (circulation). This is the essence of the core spreading vortex method.

It is not explained how the authors use the two-dimensional core spreading formula with a three-dimensional cutoff function. It is also not clear why the standard core spreading formula is given, which applies only to simple Gaussian blobs, like 7.

The compilation of articles by Uchiyama and co-workers, from (2) to the paper under discussion, seems to indicate that the core spreading method has been misinterpreted. I tend to believe that the method that has been implemented is rather like a point vortex model, regularized. It seems quite probable that viscous effects are incorrectly modeled; the method seems to be mathematically inconsistent.

It is important, in my opinion, to discuss and acknowledge the problems just examined. The vortex method has sometimes suffered the reputation of being only a modeling approach for unsteady flows. It is a fact that the vortex method can produce direct numerical simulation results, and many researchers have shown accuracy in competition with finite difference methods and even spectral methods (see, for example, (18,19,20,21,22)). One should not, under the cover of considering the “application to engineering problems,” fail to have concern for numerical consistency and accuracy.

## Acknowledgment

For reading an early draft of this discussion, thanks to A. Leonard, C. Meneveau, and L. F. Rossi.