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]
The method
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.
Example
Consider 2D Euler equations, modified using the confinement term, F:
\( {\displaystyle {\frac {\partial \mathbf {u} }{\partial t}}+\mathbf {u} \cdot \nabla \mathbf {u} +\nabla {\frac {P}{\rho }}=F_{D}(\mathbf {u} )-F_{C}(\mathbf {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, \( {\displaystyle F_{D}} \), to the partial differential equation, which when balanced with inward convection, \( {\displaystyle F_{C}} \), produce stable solutions. Another form is termed as VC2 in which dissipation is balanced with nonlinear anti-diffusion to produce stable solitary wave-like solutions.
\( {\displaystyle F_{D}} \) : Dissipation
\( {\displaystyle F_{C}} \) : Inward convection for VC1 and nonlinear anti-diffusion for VC2
The main difference between VC1 and VC2 is that in the latter the centroid of the vortex follows the local velocity moment weighted by vorticity. This should provide greater accuracy than VC1 in cases where the convecting field is weak compared to the self-induced velocity of the vortex. One drawback is that VC2 is not as robust as VC1 because while VC1 involves convection like inward propagation of vorticity balanced by an outward second order diffusion, VC2 involves a second order inward propagation of vorticity balanced by 4th order outward dissipation. This approach has been further extended to solve wave equation and is called Wave confinement (WC).
Immersed boundary
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.
Applications
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.
References
John Steinhoff (1994). "Vorticity Confinement: A New Technique for Computing Vortex Dominated Flows". Frontiers of Computational Fluid Dynamics. John Wiley & Sons. ISBN 978-0-471-95334-0.
Hu, Guangchu; Grossman, Bernard (2006-08-01). "The computation of massively separated flows using compressible vorticity confinement methods". Computers & Fluids. 35 (7): 781–789. doi:10.1016/j.compfluid.2006.03.001. ISSN 0045-7930.
Wenren, Y.; Fan, M.; Dietz, W.; Hu, G.; Braun, C.; Steinhoff, J.; Grossman, B. (2001-01-08). "Efficient Eulerian computation of realistic rotorcraft flows using Vorticity Confinement - A survey of recent results". 39th Aerospace Sciences Meeting and Exhibit. doi:10.2514/6.2001-996.
Murayama, Mitsuhiro; Nakahashi, Kazuhiro; Obayashi, Shigeru (2001-01-08). "Numerical simulation of vortical flows using vorticity confinement coupled with unstructured grid". 39th Aerospace Sciences Meeting and Exhibit. doi:10.2514/6.2001-606.
Bishop, A.R.; Krumhansl, J.A.; Trullinger, S.E. (1980). "Solitons in condensed matter: A paradigm". Physica D: Nonlinear Phenomena. 1 (1): 1–44. doi:10.1016/0167-2789(80)90003-2. ISSN 0167-2789.
Wenren, Y.; Fan, M.; Wang, L.; Xiao, M.; Steinhoff, J. (2003). "Application of Vorticity Confinement to Prediction of the Flow over Complex Bodies". AIAA Journal. 41 (5): 809–816. doi:10.2514/2.2042. ISSN 0001-1452.
vte
Numerical methods for partial differential equations
Finite difference
Parabolic
Forward-time central-space (FTCS) Crank–Nicolson
Hyperbolic
Lax–Friedrichs Lax–Wendroff MacCormack Upwind Method of characteristics
Others
Alternating direction-implicit (ADI) Finite-difference time-domain (FDTD)
Finite volume
Godunov High-resolution Monotonic upstream-centered (MUSCL) Advection upstream-splitting (AUSM) Riemann solver essentially non-oscillatory (ENO) weighted essentially non-oscillatory (WENO)
Finite element
hp-FEM Extended (XFEM) Discontinuous Galerkin (DG) Spectral element (SEM) Mortar Gradient discretisation (GDM) Loubignac iteration Smoothed (S-FEM)
Meshless/Meshfree
Smoothed-particle hydrodynamics (SPH) Moving particle semi-implicit method (MPS) Material point method (MPM) Particle-in-cell (PIC)
Domain decomposition
Schur complement Fictitious domain Schwarz alternating
additive abstract additive Neumann–Dirichlet Neumann–Neumann Poincaré–Steklov operator Balancing (BDD) Balancing by constraints (BDDC) Tearing and interconnect (FETI) FETI-DP
Others
Spectral Pseudospectral (DVR) Method of lines Multigrid Collocation Level-set Boundary element Immersed boundary Analytic element Isogeometric analysis Infinite difference method Infinite element method Galerkin method
Petrov–Galerkin method Validated numerics Computer-assisted proof Integrable algorithm Method of fundamental solutions
Hellenica World - Scientific Library
Retrieved from "http://en.wikipedia.org/"
All text is available under the terms of the GNU Free Documentation License
