Date of Award


Degree Type


Degree Name

Doctor of Philosophy


Mechanical Engineering

Major Professor

John Steinhoff

Committee Members

Gary A. Flandro, K. C. Reddy, Lloyd M. Davis


Vorticity Confinement (VC) is a computational technique used to compute fluid flows with thin vortical regions in an accurate and efficient manner. Recent results have shown that it accurately computes the turbulent wake behind blunt bodies at large Reynolds numbers. Physically, these flows are dominated by thin vortices that convect downstream with the flow, which accounts for the important higher-order velocity statistics. These thin vortices are artificially dissipated with traditional methods, whereas VC can accurately compute them.

VC consists of a set of discrete equations with the confinement appearing as a source term in the discretized momentum equations. The term is formulated to act solely on the small vortical regions of the flow which are artificially spread with traditional methods. The role of the confinement term is to contract the vortices to counteract the artificial spreading, while maintaining the conservation of properties.

Previous VC investigations into turbulence utilized a constant user-specified confinement strength which is analogous to the Boussinesq assumption. Here, two new models for free turbulence are developed and validated along with a near-wall turbulence model. The models replace the user-specified confinement strength with a confinement strength based upon properties of the flows. The models choose the most appropriate confinement strength to conserve energy, and in this way the new models are analogous to dynamic sub-grid models of LES. These models are developed such that energy behaves in a physically consistent manner for high Reynolds number turbulence-energy decays solely at the smallest scales.

The new models for free turbulence are validated against experiment, numerical and analytical solutions for the Taylor-Green vortex and the decay of forced homogenous turbulence. Results indicate that the inertial range can extend to the largest wavenumber resolvable by the computational grid. A new near-wall model is validated for the case of flat plate boundary layer flow. A final investigation for the flow over a back-facing step is also conducted and compared against experiment. VC results are found to be comparable to other LES computations but far more economical, since a coarser grid is used.

Files over 3MB may be slow to open. For best results, right-click and select "save as..."