CFD Simulation
CFD Simulation
CFD Simulation
SUMMARY
We propose a method to visualize vortex cores based on manipulation of the pressure field
produced by isolated vortices in incompressible flow. Under ideal conditions the function
D 2 p 2 p yields an approximate distance to vortex centerlines. As opposed to local methods
to identify coherent structures, isosurfaces of D produce a field of vortex tubes equidistant to the
vortex core center which, ideally, are independent of vortex intensity or size. In contrast to other
line-vortex identification methods, which typically rely on algorithms to detect vortex core lines
and frequently need complex implementations, the proposed method can be computed from the
local Eulerian velocity and pressure fields as easily as vortex identification methods such as the Q
and λ2 criteria. D 2 p 2 p results in the exact distance to the core center for a Rankine vortex
and is in general valid for the region of a vortex where there is pure rotation, yielding an
approximation to the distance farther from the core in other simple one-dimensional vortex models.
The methodology performs well in all tests we attempted, though limitations are presented and
discussed. The method is demonstrated for a canonical Burgers vortex, a Bodewadt vortex,
homogeneous isotropic turbulent flow, the wake of a propeller, a heaving plate, a turning
containership and the airwake of a surface combatant. The proposed method helps to better
visualize vortical flow fields by displaying vortex cores, complementing methods like Q and λ2
which display vortical volumes.
KEY WORDS: Vortices, Computational Fluid Dynamics.
A vortex filament is generally recognized as
a region where rigid body rotation has a
1. Introduction
maximum and constitutes the center of the
While the definition of a vortex is still being vortex. The definition of a vortex core is more
debated, the need to visualize flow fields has complicated since it involves identifying the
led to the development of a considerable boundary between the rotational and the
number of vortex identification methods. irrotational flow, becoming subjective (see
Though none of these methods is perfect or different points of view on this in [1-3]). This
universal, due to limitations and a natural definition of the interface between rotational
consequence of the lack of a vortex definition, and irrotational flow is already arbitrary for
many of them have been in use for a number of simple one-dimensional vortex considering the
years, and are frequently tested by users to find diffusion of vorticity by viscosity, and
the method that best fits their particular becomes even more so for complex three
problem and needs. dimensional flows where vortices may not be
symmetric, and are subject to interaction with
A simple approach for vortex core visualization 2
Qt will produce false positives, and overly the center line crosses the element. This
busy structures. A high Qt , on the other hand, algorithm is simple and easy to parallelize,
though it tends to fail for vortex core lines with
will display less and tighter structures, clearing high curvature [12]. This problem was studied
the picture and allowing easier observation of and solutions proposed by Roth [13] in terms
the stronger vortices, but weaker vortices are of parallel vector operators, which can be used
lost when they may be relevant. to express most definitions of vortex cores. By
Other Eulerian local vortex identification choosing appropriate pairs of vectors and
methods have been proposed with various Boolean operators, Roth [13] was able to group
degrees of success. There is no consensus on the methods of Sujudi & Haimes [11], Levy et
what method is ultimately superior; instead all al. [14], Banks & Singer [15], Kida & Miura
have strengths and weaknesses. Jeong & [16] and Strawn et al. [17] into a single
Hussain [3] proposed the 2 -criterion, which algorithm implementation framework.
uses a negative second eigenvalue of the To our knowledge, all current vortex core
second order tensor Ω2 S2 to define a vortex. line detection methods fall in the category of
The 2 -criterion, as the Q-criterion, is widely algorithmic methods, with different degrees of
difficulty. Most of the proposed schemes can
used and relatively simple to evaluate, though
be very complex to implement and rarely used.
computation of the eigenvalues of Ω2 S2 and The method to visualize vortex core proposed
ordering is more involved than the in this paper, though it cannot find the exact
computation of Q, which only requires center of the vortex, can be used to
derivatives of the velocity. Chakraborty et al. approximate the core lines using isosurfaces of
[1] introduced the ci / cr enhanced swirling small D as shown in subsequent sections. In
strength criterion, an improvement to the addition, it is Eulerian and does not require
original ci swirling strength criterion method specific implementations, and can be used with
of Zhou et al. [10]. the capabilities to process velocity and
pressure fields available in most commercial or
1.2. Vortex core line identification open source software to postprocess flow
methods solutions.
Vortex core line identification methods are
mostly based on local or global algorithms that 2. Distance to a vortex core line
detect the center of rotation of vortices. One of
the most widely used methods was developed In this section we present a methodology to
by Sujudi & Haimes [11], and involves produce vortex tubes at a given distance
detection of a center of rotation in individual around vortex center lines as isosurfaces of a
cells. The velocity field is interpolated on each function D obtained from the Eulerian local
tetrahedra and v is computed. If the velocity and pressure.
eigenvalues of v has one real R and two 2.1. The distance function D
complex conjugate C eigenvalues, then the In absence of strong stretching by external
velocity field is projected into the plane normal strain, the inner portion of the core of a vortex
to the eigenvector corresponding to R . The can be approximated by rigid body rotation,
process to find the center of rotation involves thus the azimuthal velocity is a linear function
detecting if the line of rotation intersects the with the radial distance to the vortex center line
tetrahedron faces, if two points are found then and the angular rotation velocity in an inertial
A simple approach for vortex core visualization 4
where P n n is the pressure gradient which in some cases yields better results than
Eq. (10), at the cost of a more complex
projected in the direction of a unit vector n
computation. Alternatively, 2 can be used to
normal to the vortex axis. n can be obtained
from the eigenvector corresponding to 2 as compute D, since Q = −λ2/2 in planar flows
defined in the 2 criterion [3]. If the pressure P n
D (12)
gradient is dominated by the pressure field max 2 / 2,
produced by the vortex, which as will be
which essentially computes Q in the plane of
discussed in §3 is in most situations the case
rotation, removing effects of pressure
since D is relevant in locations inside vortex
variations in direction of the axis of rotation.
cores, we can then simplify Eq. (7) to
Eq. (12), however, is not immune to external
P pressure fields in the plane of rotation.
D (8)
1 2
P 2.2. Simple vortex models
2
Gerz et al. [18] presents several one-
since in this case the pressure gradient and
dimensional vortex models in the context of a
Laplacian are dominated by the local vortex.
study on aircraft wing tip vortices. All models
Though Eq. (8) can be computed if the pressure
share two important properties. Since the
field is available from either CFD or
center of the core rotates as a rigid body, the
experiments, it involves second derivatives of
vorticity gradient normal to the vortex center
the pressure, and in some cases may result in
A simple approach for vortex core visualization 5
line is zero at the vortex axis, r 0 , [18], but that results in a deviation from the
implying that the azimuthal velocity increases rigid body rotational speed for r 0 .
linearly with r , v ~ r . The second is that the
vorticity vanishes far from the core, and thus
the azimuthal velocity decreases as v ~ 1/ r .
The simplest vortex model is the Rankine
vortex, where the vorticity is constant within
the vortex core, and zero outside,
r r 1
v (13)
1 / r r 1
where the velocity is nondimensionalized with
/ 2 a , the azimuthal velocity at a
characteristic vortex core size a as defined in
Leweke et al. [19]. a is used as reference
length for the equations in dimensionless form.
The Rankine vortex suffers from a
discontinuity in the velocity gradient at the
vortex core edge r 1 . A more realistic model
is the Lamb-Oseen vortex model,
v 1/ r 1 er
2
(14)
D
r 1 e r
2
(17)
incurred when the D function is displayed.
Adding straining to the Lamb-Oseen vortex
e r 2r 2 1 1
2
results in the Burgers vortex model, a good
example of a vortex affected by an external
The azimuthal velocity v and the function D pressure gradient
for Rankine and Lamb-Oseen vortices are
vr r , v 1/ r 1 e r , vz 2 z (18)
2
straining velocity field vr , vz respect to the (see Fig. 2, top). For 0.05 the difference
strength of the vortical structure represented by between the total pressure gradient and the
v . The resulting total dimensionless pressure pressure gradient due to the azimuthal velocity
is 0.03 at r 0 , and quickly decreases to less
gradient is shown in Fig. 2 for different and than 0.01 for larger r . A stronger straining
z values. For weak straining, 0.05 , the field causes a considerable pressure field,
external pressure gradient is negligible which causes Eq. (10) to predict a deformed
compared to the pressure gradient induced by distance function D, while Eq. (7) still
the flow rotation, resulting in no deformation produces good results (see Fig. 2, bottom). In
of the predicted tubes using the approximate most practical cases use of Eq. (7) produces
Eq. (10) respect to predictions with the more results that are very similar to those produced
accurate Eq. (7) (see Fig. 3). A stronger by Eq. (10). While more accurate and immune
straining flow with 0.2 causes a high to effects of external pressure fields on the
pressure gradient in the axial direction of the ability to display vortex tubes, Eq. (7) will in
vortex z , resulting in higher total pressure most cases result in noisier images caused by
the computation of the second derivatives of
gradient as z increases. The result is a pressure.
distortion of the isosurface of D further out
from the straining plane at higher z , deviating
from the desired vortex tube. Notice that using
only the pressure gradient normal to the axis,
Eq. (7), the vortex tube is recovered with little
effect on the distance to the center, as shown in
Fig. 3. This is explained by the little effect that
the straining field has in the radial pressure
gradient, depicted in Fig. 2.
P vr2 v2 P P vz2
, 0, (19)
r r r z z
In the Burgers vortex the radial inflow
produces a pressure gradient caused by the
straining velocity field, independent of the
vortical structure represented by v , resulting
in a simple example where a pressure field,
external to the vortex, affects the behaviour of
different methodologies to compute D. It must
be noted that Q is not affected by the straining
velocity field, with Q 0 isosurfaces yielding
a cylinder of radius R 1.1a . In the case of
low to moderate strain, the pressure gradient is
dominated by the azimuthal component of the
velocity and isosurfaces of D 0.4a from Eq.
Fig. 2: Burgers vortex: total pressure gradient (top) and
(10) produce a tube with actual distance to the radial pressure gradient and velocity (bottom) for
vortex axis of approximately R 0.355a , various values of and z .
indistinguishable from results using Eq. (7)
A simple approach for vortex core visualization 7
ACKNOWLEDGEMENTS
This work was partially supported by the US
Office of Naval Research grants N00014-17-1-
2082 and N00014-17-1-2293, Drs. Thomas Fu
and Ki-Han Kim program officers.
REFERENCES
Fig. 9: Propeller-rudder interaction for containership in
zigzag maneuver. Isosurfaces of D 0.001 (top) and
Q 15000 (bottom) coloured with dimensionless
[1] Chakraborty, P., Balachandar, S., and
vorticity magnitude.
Adrian, R. J., 2005, "On the relationships
between local vortex identification schemes,"
J. Fluid Mech., 535, pp. 189-214.
4. Conclusion [2] Haller, G., 2005, "An objective definition
of a vortex," J. Fluid Mech., 525, pp. 1-26.
We discussed a method to compute a
function D that provides the approximate [3] Jeong, J., and Hussain, F., 1995, "On the
distance to a vortex centerline for Identification of a Vortex," J. Fluid Mech.,
incompressible flows. The method assumes 285, pp. 69-94.
that the pressure gradient is mainly caused by [4] Epps, B. P., 2017, "Review of Vortex
the local vortex, as opposed to pressure Identification Methods," AIAA SciTec Forum.
gradients caused by other nearby vortices or
flow boundaries. This effect of external
A simple approach for vortex core visualization 14