A vortex-panel method for potential flow is used as a basis for modeling surface effects and creating a finite-element interface so that an arbitrary body can be analyzed. The basic model consists of triangular panels of linearly varying vorticity which represent the body, vortex cores on the lifting edges of the body, and vortex filaments representing the wake. The interface modification is made by using a finite-element application's output as the basis for an input file for the model, executing the main program, and writing body and wake output readable by the finite-element application. The surface-effect modification is made by including an image of the body below the real body to create a surface boundary condition through symmetry. / Master of Science

A unique multifrequencied transfer matrix method performs three-dimensional harmonic, steady-state response calculations on geared-rotor systems. The full six degrees-of-freedom method includes physical branching to accommodate multiple shafting and frequency branching to simultaneously accommodate multiple frequencies and their interdependence resulting from time-varying mesh stiffness. Areas of emphasis include development of a modified transfer matrix to handle multiple frequencies and shafting; description of the time-varying stiffness tensor representing the involute spur gear mesh based on bending, shear, compression, and local contact deformation; development of the mesh transfer matrix; development of an automatic system solver to allow the engineer to analyze systems of arbitrary construction; and the development of a matrix solver to efficiently handle large systems. A computer analysis demonstrates the significance of terms included in the stiffness evaluation as compared with less rigorous treatment in the literature. An analytical example problem illustrates the automated model generation through complete rotor system dynamic response analysis produced by the current work with special attention to the significance of parametric excitation due to the gear mesh. / Ph. D.

No description available.

In this thesis I study a variety of two-dimensional turbulent systems using a mixed analytical, phenomenological and numerical approach. The systems under consideration are governed by the two-dimensional Navier-Stokes (2DNS), surface quasigeostrophic (SQG), alpha-turbulence and magnetohydrodynamic (MHD) equations. The main analytical focus is on the number of degrees of freedom of a given system, defined as the least value $N$ such that all $n$-dimensional ($n$ ≥ $N$) volume elements along a given trajectory contract during the course of evolution. By equating $N$ with the number of active Fourier-space modes, that is the number of modes in the inertial range, and assuming power-law spectra in the inertial range, the scaling of $N$ with the Reynolds number $Re$ allows bounds to be put on the exponent of the spectrum. This allows the recovery of analytic results that have until now only been derived phenomenologically, such as the $k$[superscript(-5/3)] energy spectrum in the energy inertial range in SQG turbulence. Phenomenologically I study the modal interactions that control the transfer of various conserved quantities. Among other results I show that in MHD dynamo triads (those converting kinetic into magnetic energy) are associated with a direct magnetic energy flux while anti-dynamo triads (those converting magnetic into kinetic energy) are associated with an inverse magnetic energy flux. As both dynamo and anti-dynamo interacting triads are integral parts of the direct energy transfer, the anti-dynamo inverse flux partially neutralises the dynamo direct flux, arguably resulting in relatively weak direct energy transfer and giving rise to dynamo saturation. These theoretical results are backed up by high resolution numerical simulations, out of which have emerged some new results such as the suggestion that for alpha turbulence the generalised enstrophy spectra are not closely approximated by those that have been derived phenomenologically, and new theories may be needed in order to explain them.

The relative rates of detrital decomposition in four vegetation communities within the Everglades' ridge-slough microtopography were evaluated during two trials. Litterbags with community-specific detritus in proportion to each community's composition were put into the four communities; namely, submerged marsh, emergent marsh, short Cladium ridge, and tall Cladium ridge. These litterbags were paired with litterbags containing control leaf litter from Chrysobalanus icaco and Salix caroliniana during the wet and dry season trials, respectively. No regional differences in decomposition were shown, but there were significant differences across communities, attributed to the initial C:N ratio of the detritus, with the fastest decomposition occurring in the deepest submerged marsh followed by emergent marsh, and the shallower ridge communities had equally slower decomposition. Additionally, both controls followed the same pattern. Thus, decomposition contributes to an active self-maintenance mechanism within the vegetation communities which ultimately helps to conserve the ridges and sloughs. / by Sheryl R. van der Heiden. / Thesis (M.S.)--Florida Atlantic University, 2008. / Includes bibliography. / Electronic reproduction. Boca Raton, Fla., 2008. Mode of access: World Wide Web.

In this dissertation, a group of vehicle dynamics simulation tools is developed with two primary goals: to accurately represent vehicle behavior and to provide insight that improves the understanding of vehicle performance. Three tools are developed that focus on tire modeling, vehicle modeling and lap time simulation. Tire modeling is based on Nondimensional Tire Theory, which is extended to provide a flexible model structure that allows arbitrary inputs to be included. For example, rim width is incorporated as a continuous variable in addition to vertical load, inclination angle and inflation pressure. Model order is determined statistically and only significant effects are included. The fitting process is shown to provide satisfactory fits while fit parameters clearly demonstrate characteristic behavior of the tire. To represent the behavior of a complete vehicle, a Nondimensional Tire Model is used, along with a three degree of freedom vehicle model, to create Milliken Moment Diagrams (MMD) at different speeds, longitudinal accelerations, and under various yaw rate conditions. In addition to the normal utility of MMDs for understanding vehicle performance, they are used to develop Limit Acceleration Surfaces that represent the longitudinal, lateral and yaw acceleration limits of the vehicle. Quasi-transient lap time simulation is developed that simulates the performance of a vehicle on a predetermined path based on the Limit Acceleration Surfaces described above. The method improves on the quasi-static simulation method by representing yaw dynamics and indicating the vehicle's stability and controllability over the lap. These improvements are accomplished while maintaining the simplicity and computational efficiency of the two degree of freedom method. / Graduation date: 2013

In this study, numerical simulations of a low emission lab-scale non-premixed combustor are conducted and analyzed. The objectives are to provide new insight into the physical phenomena in the SPRF (Stagnation Point Reverse Flow) combustor built in the Georgia Tech Combustion Lab, and to compare three Large Eddy Simulation (LES) combustion models (Eddy Break-Up [EBU], Steady Flamelet [SF] and Linear Eddy Model [LEM]) for non-premixed combustion. The nominal operating condition of the SPRF combustor achieves very low NOx and CO emissions by combining turbulent mixing of exhaust gases with preheated reactants and chemical kinetics. The SPRF numerical simulation focuses on capturing the complex interaction between turbulent mixing and heat release. LES simulations have been carried out for a non-reactive case in order to analyze the turbulent mixing inside the combustor. The LES results have been compared to PIV experimental data and the code has been validated. The dominating features of the operational mode of the SPRF combustor (dilution of hot products into reactants, pre-heating and pre-mixing) have been analyzed, and results from the EBU-LES, SF-LES and LEM-LES simulations have been compared. Analysis shows that the LEM-LES simulation achieves the best agreement with the observed flame structure and is the only model that captures the stabilization processes observed in the experiments. EBU-LES and SF-LES do not predict the correct flow pattern because of the inaccurate modeling of sub-grid scale mixing and turbulence-combustion interaction. Limitations of these two models for this type of combustor are discussed.

Large Eddy Simulation (LES) of bluff body stabilized premixed and partially premixed combustion close to the flammability limit is carried out in this thesis. The LES algorithm has no ad-hoc adjustable model parameters and is able to respond automatically to variations in the inflow conditions. Algorithm validation is achieved by comparison with reactive and non-reactive experimental data. In the reactive flow, two scalar closure models, Eddy Break-Up (EBULES) and Linear Eddy Mixing (LEMLES), are used and compared. Over important regions, the flame lies in the Broken Reaction Zone regime. Here, the EBU model assumptions fail. The flame thickness predicted by LEMLES is smaller and the flame is faster to respond to turbulent fluctuations, resulting in a more significant wrinkling of the flame surface. As a result, LEMLES captures better the subtle effects of the flame-turbulence interaction. Three premixed (equivalence ratio = 0.6, 0.65, and 0.75) cases are simulated. For the leaner case, the flame temperature is lower, the heat release is reduced and vorticity is stronger. As a result, the flame in this case is found to be unstable. In the rich case, the flame temperature is higher, and the spreading rate of the wake is increased due to the higher amount of heat release Partially premixed combustion is simulated for cases where the transverse profile of the inflow equivalence ratio is variable. The simulations show that for mixtures leaner in the core the vortical pattern tends towards anti-symmetry and the heat release decreases, resulting also in instability of the flame. For mixtures richer in the core, the flame displays sinusoidal flapping resulting in larger wake spreading. More accurate predictions of flame stability will require the use of detailed chemistry, raising the computational cost of the simulation. To address this issue, a novel algorithm for training Artificial Neural Networks (ANN) for prediction of the chemical source terms has been implemented and tested. Compared to earlier methods, the main advantages of the ANN method are in CPU time and disk space and memory reduction.

Models of the fate and transport of river plumes and the bacteria they carry into lakes are developed. They are needed to enable informed decisions about beach closures to avoid economic losses, and to help design water intakes and operate combined sewer overflow schemes to obviate exposure of the public to potential pathogens. This study advances our understanding of river plumes dynamics in coastal waters by means of field studies and numerical techniques. Extensive field measurements were carried out in the swimming seasons of 2006 and 2007 on the Grand River plume as it enters Lake Michigan. They included simultaneous aerial photography, measurements of lake physical properties, the addition of artificial tracers to track the plume, and bacterial sampling. Our observed results show more flow classes than included in previous studies (e.g. CORMIX). Onshore wind can have a significant effect on the plume and whether it impacts the shoreline. A new classification scheme based on the relative magnitude of plume-crossflow length scale and Richardson number based on the wind speed is devised. Previous studies on lateral spreading are complemented with a new relationship in the near field. The plume thickness decreased rapidly with distance from the river mouth and a new non-dimensional relationship to predict thickness is developed. Empirical near field models for surface buoyant plumes are reviewed and a near field trajectory and dilution model for large aspect ratio surface discharge channels is devised. Bacterial reductions due to dilution were generally small (less than 10:1) up to 4.5 km from the river mouth. E. coli decay rates were significantly affected by solar radiation and ranged from 0.2 to 2.2 day-1 which were within the range of previous studies in Lake Michigan. Total coliform survived longer than E. coli suggesting different die-off mechanisms. Mathematical models of the bacterial transport are developed that employ a nested modeling scheme to represent the 3D hydrodynamic processes of surface river discharges in the Great Lakes. A particle tracking model is used that provides the capability to track a decaying tracer and better quantify mixing due to turbulent diffusion. Particle tracking models have considerable advantages over gradient diffusion models in simulating bacterial behavior nearshore that results in an improved representation of bacteria diffusion, decay and transport. Due to the complexity and wide variation of the time and length scale of the hydrodynamic and turbulent processes in the near field (where plume mixing is dominated by initial momentum and buoyancy) and far field (where plume mixing is dominated by ambient turbulence), a coupling technique is adapted. The far field random walk particle tracking model incorporates the empirical near field model. It simulates the transport, diffusion and decay of bacteria as discrete particles and employs the near field output as the source and transports the particles based on ambient currents predicted by the 3D hydrodynamic model. The coupled model improves dilution predictions in the near field. The new techniques advance our knowledge of the nearshore fate and transport of bacteria in the Great Lakes and can be ultimately applied to the NOAA Great Lakes Coastal Forecasting System to provide a reliable prediction tool for bacterial transport in recreational waters.

The aim of this study is to investigate means of efficiently assessing the effects of distributed structural modification on the dynamic properties of a complex structure. The helicopter structure is normally designed to avoid resonance at the main rotor rotational frequency. However, very often military helicopters have to be modified (such as to carry a different weapon system or an additional fuel tank) to fulfill operational requirements. Any modification to a helicopter structure has the potential of changing its resonance frequencies and mode shapes. The dynamic properties of the modified structure can be determined by experimental testing or numerical simulation, both of which are complex, expensive and time-consuming. Assuming that the original dynamic characteristics are already established and that the modification is a relatively simple attachment such as beam or plate modification, the modified dynamic properties may be determined numerically without solving the equations of motion of the full-modified structure. The frequency response functions (FRFs) of the modified structure can be computed by coupling the original FRFs and a delta dynamic stiffness matrix for the modification introduced. The validity of this approach is investigated by applying it to several cases, 1) 1D structure with structural modification but no change in the number of degree of freedom (DOFs). A simply supported beam with double thickness in the middle section is treated as an example for this case; 2) 1D structure with additional DOFs. A cantilever beam to which a smaller beam is attached is treated as an example for this case, 3) 2D structure with a reduction in DOFs. A four-edge-clamped plate with a cut-out in the centre is treated as an example for this case; and 4) 3D structure with additional DOFs. A box frame with a plate attached to it as structural modification with additional DOFs and combination of different structures. The original FRFs were obtained numerically and experimentally except for the first case. The delta dynamic stiffness matrix was determined numerically by modelling the part of the modified structure including the modifying structure and part of the original structure at the same location. The FRFs of the modified structure were then computed. Good agreement is obtained by comparing the results to the FRFs of the modified structure determined experimentally as well as by numerical modelling of the complete modified structure.