Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Implementation of the Turn Function Method in a three-dimensional, parallelized hydrodynamics code

Conference ·
OSTI ID:6792634
We describe the implementation of the Turn Function Method in KIVA-F90, a version of the KIVA computer program written in the FORTRAN 90 programming language that is used on some massively parallel computers. The Turn Function Method solves both linear momentum and vorticity equations in numerical calculations of compressible fluid flow. Solving a vorticity equation allows vorticity to be both conserved and transported more accurately than in traditional methods for computing compressible flow. This first implementation of the Turn Function Method in a three-dimensional hydrodynamics code involved some modification of the original method and some novel numerical difference approximations. In particular, a penalty method is used to keep the divergence of the computed vorticity field close to zero. Difference operators are also defined in such a way that the finite difference analog of [nabla] [center dot] ([nabla] [times] u) = 0 is exactly satisfied. Three example problems show the increased computational cost and the accuracy to be gained by using the Turn Function Method in calculations of flows with rotational motion. Use of the Method can increase by 60% the computational times of the Euler equation solver in KIVA-F90, but we conclude this increased cost is justified by the increased accuracy.
Research Organization:
Los Alamos National Lab., NM (United States)
Sponsoring Organization:
DOE; USDOE, Washington, DC (United States)
DOE Contract Number:
W-7405-ENG-36
OSTI ID:
6792634
Report Number(s):
LA-UR-92-3468; CONF-930372--1; ON: DE93003744
Country of Publication:
United States
Language:
English