Thesis (M.S.)--Oregon State University, 2008. / Printout. Includes bibliographical references (leaves 145-148). Also available on the World Wide Web.

Thesis (Ph. D., Mechanical Engineering)--University of Idaho, June 2009. / Major professor: Akira Tokuhiro. Includes bibliographical references. Also available online (PDF file) by subscription or by purchasing the individual file.

Thesis (Ph.D)--Aerospace Engineering, Georgia Institute of Technology, 2009. / Committee Chair: Mavris, Dimitri; Committee Member: Alley, Nicholas; Committee Member: Lakshmi, Sankar; Committee Member: Sriram, Rallabhandi; Committee Member: Stephen, Ruffin. Part of the SMARTech Electronic Thesis and Dissertation Collection.

Thesis (M.S.)--Brigham Young University. Dept. of Mechanical Engineering, 2008. / Includes bibliographical references (p. 79-82).

Thesis (Ph. D.)--Worcester Polytechnic Institute. / UMI no.: 99-37530. Includes bibliographical references (leaves 221-238).

Thesis (M.S. in Mechanical Engineering) Naval Postgraduate School, March 1997. / Thesis advisor, Ashok Gopinath. Includes bibliographical references (p. 47-48). Also available online.

Thesis supervisors: Larry J. Weber, Nathan C. Young. Includes bibliographic references (p. 152-155).

Thesis supervisors: Marcela Politano, Larry J. Weber. Includes bibliographic references (p. 95-97).

The marker and cell method is an efficient co-volume technique suitable for the solution of incompressible flows using a Cartesian mesh. For flows around complex geometries the use of an unstructured mesh is desirable. For geometric flexibility an unstructured mesh implementation is desirable. A co-volume technique requires a dual orthogonal mesh, in the triangular case the Delaunay-Voronoi dual provides the means for determining this dual orthogonal mesh in an unstructured mesh framework. Certain mesh criteria must be placed on the Delaunay-Voronoi to ensure it meets the dual orthogonal requirements. The two dimensional extension of the marker and cell method to an unstructured framework is presented. The requirements of the mesh are defined and methods in their production are discussed. Initially an explicit time stepping scheme is implemented which allows efficient simulation of incompressible fluid flow problems. Limitations of the explicit time stepping scheme that were discovered, mean that high Reynolds number flows that require the use of stretched meshes cannot produce solutions in a reasonable time period. A semi-implicit time stepping routine removes this limitation allowing these types of flows to be successfully modelled. To validate the solvers accuracy and demonstrate its performance, a number of test cases are presented. These include the lid driven cavity, flow over a backward facing step, inviscid flow around a circular cylinder, unsteady flow around a circular cylinder, flow around an SD7003 aerofoil, flow around a NACA0012 aerofoil and flow around a multi element aerofoil. The investigation although revealing a high dependence on the quality of the mesh still demonstrates that accurate results can be obtained efficiently. The efficiency is demonstrated by comparison to the in-house 2D incompressible finite volume solver for flow around a circular cylinder. For this case the unstructured MAC method produced a solution four times faster than the finite volume code.

In this research, the finite element method (FEM) was used to solve the nonlinear, incompressible, transient, three dimensional Navier-Stokes equations in their non-conservative form. Linear tetrahedron elements were employed with the elegant, equal order interpolation for both pressure and velocity. The characteristic based split scheme was formulated in a fully implicit manner to circumvent the time step restrictions of the classical explicit formulations. The monolithic (single step, fully coupled solution procedure for pressures and velocity) form of the CBS scheme was also derived and its suitability was positively demonstrated. Casting the CBS scheme in a monolithic framework, results in the generation of a pressure stabilization term in the mass conservation equation, thereby circumventing the LBB restriction by the elimination of the zero pressure block. An account of all the steps involved in discretizing the Navier-Stokes equations (both in split and monolithic frameworks) was presented in meticulous detail, which included the derivation of the convective and pressure stabilization terms, linearization of the non-linear terms and the consequent derivation of the highly efficient analytical jacobian matrix, along with the temporal and spatial discretizations of the corresponding terms. The monolithic and the split version of the CBS scheme were integrated into a parallel, scalable and extensible Fortran90 software called IFENs. The development of IFENs started during the course of this research and all of its components have been designed and implemented by the author of this thesis. Multi processor parallelism was achieved using the Intel implementation of the most widely used/preferred, Message Passing Interface (MPI) standard. The parallel support needed for the use of a variety of parallel, linear, iterative solvers belonging to the Krylov subspace family (e.g. GMRES and its variants, CG, BiCG, BiCG- stab, etc.), parallel non linear solvers belonging to the Newton-Krylov family (line search newton, trust region newton, nonlinear GMRES, etc.) and parallel preconditioners (incomplete LU, Additive Shwarz Method - ASM, algebraic multigrid, etc.), was provided by the incorporation of PETSc into IFENs. PETSc is a state of the art, non-trivial toolkit, which represents a collection of several parallel libraries useful in high performance scientific computing. Keeping in mind the specific requirements of IFENs, a custom mesh partitioner was implemented. It operated on meshes that were renumbered using bandwidth reducing algorithms like Revere Cuthill Mckee. The possibility of using established domain decomposition libraries like ParMETIS was explored and demonstrated to be counter productive for the demands of this research. After the preliminary testing and validation of the procedures adopted before and during the execution of IFENs, large, high definition domains representative of human arteries (specifically, carotid bifurcations, found in the neck) were considered and the complete incompressible set of Navier-Stokes equations were solved for pressure and velocity fields. During the tenure of this research more than 1000 recorded parallel test cases were executed to test various components of IFENs, as well as various simulations representative of a wide variety of problems. IFENs can easily handle meshes with tens of millions of elements. The largest mesh used for the purpose of this research contained 14.58 million tetrahedrons and 2.489 million nodes, which on average required just 7 minutes per timestep, while executing the classical split framework of the CBS scheme. Results from the simulation of 9 carotid meshes, representative of 4 carotid geometries were presented and found to be in good agreement with the available ultrasound data. The flow fields were analysed and post processed using different techniques for each case. The haemodynamic wall parameters like time averaged wall shear stress and oscillatory shear index were calculated and mapped onto the corresponding boundary nodes. The region in the carotid bifurcation susceptible to the deposition of plaques and consequent stenosis were pointed out and other anomalies were highlighted.