Vorticity confinement (VC), a physics-based computational fluid dynamics model analogous to shock capturing methods, was invented by Dr. John Steinhoff, professor at the University of Tennessee Space Institute, in the late 1980s[1] to solve vortex dominated flows. It was first formulated to capture concentrated vortices shed from the wings, and later became popular in a wide range of research areas.[2] During the 1990s and 2000s, it became widely used in the field of engineering.[3] [4]
VC has a basic familiarity to solitary wave approach which is extensively used in many condensed matter physics applications.[5] The effect of VC is to capture the small scale features over as few as 2 grid cells as they convect through the flow. The basic idea is similar to that of compression discontinuity in Eulerian shock capturing methods. The internal structure is maintained thin and so the details of the internal structure may not be important.
Consider 2D Euler equations, modified using the confinement term, F:
| \partialu | |
| \partialt |
+u ⋅ \nablau+\nabla
| P | |
| \rho |
=FD(u)-FC(u)
The discretized Euler equations with the extra term can be solved on fairly coarse grids, with simple low order accurate numerical methods, but still yield concentrated vortices which convect without spreading. VC has different forms, one of which is VC1. It involves an added dissipation,
FD
FC
FD
FC
To enforce no-slip boundary conditions on immersed surfaces, first, the surface is represented implicitly by a smooth “level set” function, “f”, defined at each grid point. This is the (signed) distance from each grid point to the nearest point on the surface of an object – positive outside, negative inside. Then, at each time step during the solution, velocities in the interior are set to zero. In a computation using VC, this results in a thin vortical region along the surface, which is smooth in the tangential direction, with no “staircase” effects.[6] The important point is that no special logic is required in the “cut” cells, unlike many conventional schemes: only the same VC equations are applied, as in the rest of the grid, but with a different form for F. Also, unlike many conventional immersed surface schemes, which are inviscid because of cell size constraints, there is effectively a no-slip boundary condition, which results in a boundary layer with well-defined total vorticity and which, because of VC, remains thin, even after separation. The method is especially effective for complex configurations with separation from sharp corners. Also, even with constant coefficients, it can approximately treat separation from smooth surfaces. General blunt bodies, which typically shed turbulent vorticity that induces a velocity around an upstream body. It is inconsistent to use body fitted grids as the vorticity convects through a non fitted grid.
VC is used in many applications including rotor wake computations, computation of wing tip vortices, drag computations for vehicles, flow around urban layouts, smoke/contaminant propagation and special effects. Also, it is used in wave computations for communication purposes.