The strand-Cartesian grid approach is a unique method of generating and computing fluid dynamic simulations. The strand-Cartesian approach provides highly desirable qual- ities of fully-automatic grid generation and high-order accuracy. This thesis focuses on the implementation of the Spalart-Allmaras turbulence model to the strand-Cartesian grid framework. Verification and validation is required to ensure correct implementation of the turbulence model.Mathematical code verification is used to ensure correct implementation of new algo- rithms within the code framework. The Spalart-Allmaras model is verified with the Method of Manufactured Solutions (MMS). MMS shows second-order convergence, which implies that the new algorithms are correctly implemented.Validation of the strand-Cartesian solver is completed by simulating certain cases for comparison against the results of two independent compressible codes; CFL3D and FUN3D. The NASA-Langley turbulence resource provided the inputs and conditions required to run the cases, as well as the case results for these two codes. The strand solver showed excellent agreement with both NASA resource codes for a zero-pressure gradient flat plate and bump- in-channel. The treatment of the sharp corner on a NACA 0012 airfoil is investigated, resulting in an optimal external sharp corner configuration of strand vector smoothing with a base Cartesian grid and telescoping Cartesian refinement around the trailing edge. Results from the case agree well with those from CFL3D and FUN3D. Additionally, a NACA 4412 airfoil case is examined, and shows good agreement with CFL3D and FUN3D, resulting in validation for this case.

Development of a novel high-order flux correction method on strand grids is presented. The method uses a combination of flux correction in the unstructured plane and summation-by-parts operators in the strand direction to achieve high-fidelity solutions. Low-order truncation errors are cancelled with accurate flux and solution gradients in the flux correction method, thereby achieving a formal order of accuracy of 3, although higher orders are often obtained, especially for highly viscous flows. In this work, the scheme is extended to high-Reynolds number computations in both two and three dimensions. Turbulence closure is achieved with a robust version of the Spalart-Allmaras turbulence model that accommodates negative values of the turbulence working variable, and the Menter SST turbulence model, which blends the k-ε and k-ω turbulence models for better accuracy. A major advantage of this high-order formulation is the ability to implement traditional finite volume-like limiters to cleanly capture shocked and discontinuous flow. In this work, this approach is explored via a symmetric limited positive (SLIP) limiter. Extensive verification and validation is conducted in two and three dimensions to determine the accuracy and fidelity of the scheme for a number of different cases. Verification studies show that the scheme achieves better than third order accuracy for low and high-Reynolds number flow. Cost studies show that in three-dimensions, the third-order flux correction scheme requires only 30% more walltime than a traditional second-order scheme on strand grids to achieve the same level of convergence. In order to overcome meshing issues at sharp corners and other small-scale features, a unique approach to traditional geometry, coined "asymptotic geometry," is explored. Asymptotic geometry is achieved by filtering out small-scale features in a level set domain through min/max flow. This approach is combined with a curvature based strand shortening strategy in order to qualitatively improve strand grid mesh quality.