This thesis presents the implementation of several two-equation turbulence models into a finite difference, two- and three-dimensional Navier-Stokes Solver. Theories of turbulence modeling and the historical development of these theories are briefly investigated. Turbulence models that are defined by two partial differential equations, based on k-&amp / #969 / and k-&amp / #949 / models, having different correlations, constants and boundary conditions are selected to be adapted into the base solver. The basic equations regarding the base Navier-Stokes solver to which the turbulence models are implemented presented by briefly explaining the outputs obtained from the solver. Numerical work regarding the implementation of turbulence models into the base solver is given in steps of non-dimensionalization, transformation of equations into generalized coordinate system, numerical scheme, discretization, boundary and initial conditions and limitations. These sections of implementation are investigated and presented in detail with providing every steps of work accomplished. Certain trial problems are solved and outputs are compared with experimental data. Solutions for fluid flow over flat plate, in free shear, over cylinder and airfoil are demonstrated. Airfoil validation test cases are analyzed in detail. For three dimensional applications, computation of flow over a wing is accomplished and pressure distributions from certain sections are compared with experimental data.

Turbulent dynamics at two sites (C and D) in a hypoxic zone on the Texas- Louisiana continental shelf were studied by investigating turbulence quantities i.e. turbulence kinetic energy (TKE), dissipation rate of TKE (E), Reynolds stress (τ ), dissipation rate of temperature variance (χ), eddy diffusivity of temperature (ν't), and eddy diffusivity of density (ν'p). Numerical models were also applied to test their capability of simulating these turbulence quantities. At site D, TKE, E, and τ were calculated from velocity measurements in the bot- tom boundary layer (BBL), using the Kolmogorov’s -5/3 law in the inertial subrange of energy spectra of vertical velocity fluctuations in each burst measurement. Four second-moment turbulence closure models were applied for turbulence simulations, and modeled turbulence quantities were found to be consistent with those observed. It was found from inter-model comparisons that models with the stability functions of Schumann and Gerz predicted higher values of turbulence quantities than those of Cheng in the mid layer, which might be due to that the former stability functions are not sensitive to buoyancy. At site C, χ, E, v’t, and ν’p were calculated from profile measurements throughout the water column, and showed high turbulence level in the surface boundary layer and BBL, as well as in the mid layer where shear stress was induced by advected non-local water above a hypoxic layer. The relatively high dissolved oxygen in the non-local water resulted in upward and downward turbulent oxygen fluxes, and the bottom hypoxia will deform due to turbulence in 7.11 days. Two of the four models in the study at site D were implemented, and results showed that turbulence energy resulting from the non-local water was not well reproduced. We attribute this to the lack of high-resolution velocity measurements for simulations. Model results agreed with observations only for χ and E simulated from the model with the stability function of Cheng in the BBL. Discrepancies between model and observational results lead to the following conclusions: 1) the stability functions of Schumann and Gerz are too simple to represent the turbulent dynamics in stratified mid layers; 2) detailed velocity profiles measurements are required for models to accurately predict turbulence quantities. Missing such observations would result in underestimation,