A coupled fluid-structure interaction solver capability is developed and validated. A high fidelity fluids solver, Loci-Chem, is coupled with a finite-element structural dynamics toolkit, MAST. The coupled solver is validated for the prediction of several panel instability cases in uniform flows and in the presence of an impinging shock for a range of subsonic and supersonic Mach numbers, dynamic pressures, and pressure ratios. The panel deflections and limit-cycle oscillation amplitudes, frequencies, and bifurcation point predictions compare very well with benchmark results for 2D simulations. The same procedures outlined in the validation study have been applied to simulations of varying dynamic pressures at M = 2 for an impinging oblique shockwave. The influence of inviscid, laminar and turbulent boundary layer profiles on the development of flow field characteristics has been analyzed, and laminar predictions characterized by a large flow separation results in vastly different behavior than that of traditional flutter.

No description available.

No description available.

No description available.

No description available.

It has been shown that the classical Rosenzweig-MacArthur predator-prey model is sensitive to the functional form of the predator response. To see if this sensitivity remains in the highly controlled environment of the chemostat, we use a predator-prey model with three trophic levels and a Holling type II predator response function. We first focus on the analysis of the model using an Ivlev functional response. Local and global dynamics are studied, with global stability of the coexistence equilibrium point obtained under certain conditions. Bifurcation analysis reveals the existence of a stable periodic orbit that appears via a super-critical Hopf bifurcation. The uniqueness of this periodic orbit is explored. Finally, we make comparisons between the dynamics of the model with Ivlev response and Monod response, both of which have nearly identical graphs. The same sensitivity to functional form is observed in the chemostat as in the classical model. / Thesis / Master of Science (MSc)

This study focuses on the response of flat unsymmetric laminates to an inplane compressive loading that for symmetric laminates are of sufficient magnitude to cause bifurcation buckling, postbuckling, and secondary buckling behavior. In particular, the purpose of this study is to investigate whether or not the concept of bifurcation buckling is applicable to unsymmetric laminates. Past work by other researchers has suggested that such a concept is applicable for certain boundary conditions. The study also has as an objective the determination of the response of flat unsymmetric laminates if bifurcation buckling does not occur. The finite-element program ABAQUS is used to obtain results, and a portion of the study is devoted to becoming familiar with the way ABAQUS handles such highly geometrically nonlinear problems, particularly for composite materials and particularly when instabilities and dynamic behavior are involved. Familiarity with the problem, in general, and with the use of ABAQUS, in particular, is partially gained by considering semi-infinite unsymmetrically laminated cross- and angle-ply plates, a one-dimensional problem that can be solve in closed form and with ABAQUS by making the appropriate approximations for the infinite geometry. In this portion of the study it is found that semi-infinite cross-ply laminates with clamped boundary conditions and semi-infinite angle-ply plates with simple-support boundary conditions remain flat under a compressive load until the load magnitude reaches a certain level, at which time the out-of-plane deflection become indeterminate, essentially an eigenvalue problem as encountered with classic bifurcation buckling analyses. Obviously, a linear analysis of such problems would not reveal this behavior and, in fact, there are other revealed significant differences between the predictions of linear and nonlinear analyses. Transversely-loaded and inplane-loaded finite isotropic plates are studied by way of semi-closed form Rayleigh-Ritz-based solutions and ABAQUS in a step to approaching the problem with unsymmetric laminates. A method to investigate the unloading behavior of postbuckled finite isotropic plates is developed that reveal multiple plate configurations in the postbuckled region of the response, and this method is then extended to the study of finite inplane-loaded unsymmetric laminates. To that end, two specific laminates, a symmetric and an unsymmetric cross-ply laminates, and a variety of boundary conditions are used to study the response of inplane-loaded unsymmetric laminates. The symmetric laminate is included to provide a familiar baseline case and a means of comparison. Plates with all four edges clamped and a variety of inplane boundary conditions are studied. Of course the symmetric cross-ply laminate exhibits bifurcation behavior, and when the tangential displacement on the loaded edges and the normal displacement on the unloaded edges are restrained, secondary buckling behavior occurs. For the unsymmetric cross-ply laminate, bifurcation buckling behavior does not occur unless the tangential displacement on the loaded edges and the normal displacement on the unloaded edges are restrained, or the tangential displacement on the loaded edges and the normal displacement on the unloaded edges are free. If either of these conditions are not satisfied, the unsymmetric cross-ply laminate exhibits what could be termed 'near-bifurcation' behavior. In all cases rather complex behavior occurs for high levels of inplane load, including asymmetric postbuckling and secondary buckling behavior. For clamped loaded edges and simply-supported unloaded edges, bifurcation buckling behavior does not occur unless the tangential displacement on the loaded edges and the normal displacement on the unloaded edges are restrained. For this case, rather unusual asymmetric bifurcation and associated limit point behavior occur, as well as secondary buckling. This is a very interesting boundary condition case and is studied further for other unsymmetric cross-ply laminates, including the use of a Rayleigh-Ritz-based solution in attempt to quantify the problem parameters responsible for the asymmetric response. The overall results of the study have led to an increased understanding of the role of laminate asymmetry and boundary conditions on the potential for bifurcation behavior, on the response of the laminate for loads beyond that level. / Ph. D.

The instability modes and non-linear behavior of a cavitating flow have been studied based on the experimental data obtained from planar Particle Image Velocimetry (PIV). Three data-driven techniques, Proper Orthogonal Decomposition (POD), Dynamic Mode Decomposition (DMD), and Clustered-based Reduced Order Modeling (CROM), are applied to the snapshots of the fluctuating component of velocity to investigate instability modes of the cavitating flow. DMD and POD analysis yield multiple modes are corresponding to slow-varying drift flow, cloud-shedding, and Kelvin-Helmholtz (KH) instability for a fixed inlet flow condition. The high coherence measure obtained from the instabilities suggests a transfer of energy from the largest scales, fluctuating mean flow, to the smaller scales such as cloud cavitation and Kelvin-Helmholtz (KH) instability. It is demonstrated that the POD decorrelation of length scales yields inherently quasi-periodic time dynamics, e.g., incommensurate frequencies. Moreover, the eigenvalue obtained from DMD revealed multiple harmonic with different decay rates associated with the cloud cavitation. The above-mentioned intermittent transition between distinct cloud shedding regimes is investigated via Clustered-based Reduced Order Modeling (CROM). Four aperiodic shedding regimes are identified. 68% of the time, triplets of vortices are formed, while 28% of the time, a pair of vortices are formed in the near wake of the throat. Dominant mechanisms governing the momentum transport and the turbulence kinetic energy production, destruction, and redistribution in distinct regions of the flow field have been identified using Gaussian Mixture Models (GMMs). The preceding data-driven techniques and in-depth analysis of the results facilitated modeling of the cavitation inception and break-up. Accordingly, a phase transition field model is developed using the ultra-fast Time-Resolved Particle Image Velocimetry (TR-PIV) and vapor void fraction spatial and temporal data. The data acquisition is implemented in a Venturi-type test section. The approximate Reynolds number based upon the throat height is 10,000, and the approximate cavitation number is 1.95. The non-equilibrium cavitation model assumes that the phase production and destruction are correlated to the static pressure field, pressure spatial derivatives, void fraction, and divergence of the velocity field. Finally, the dependence of the model on the empirical constants has been investigated. / Doctor of Philosophy / A cavitation bubble occurs where the pressure field is below the saturation pressure of the liquid. Accumulation of the cavitation bubble forms a cavitating flow. This phenomenon is observed in pumps, propulsion systems, internal combustion engines, and rocket engines. Identifying the mechanisms leading to cavitation-induced damages is imperative in the design of the devices. In this regard, investigation of the cavitation bubble inception, deformation, collapse, and intermittent regime change is mandatory in learning the primary mechanisms of the stresses imposed on the device. Experiments and high-fidelity numerical and analytical methods can be employed to shed light on flow physics. The current study adopted joint experimental methods, data analysis techniques, and computational approaches to scrutinize the unsteady cavitating flow underlying physics as it occurs past the throat of a Venturi-type profile. Different mechanisms of instabilities are identified by applying the data-driven techniques to the raw images of the cavitating flow. The path of the transitions between physically different instabilities mechanisms is examined. The local dominant balance between stress terms in the conservation of momentum equation is identified, and the stress terms roles in cavitating flow instabilities and advective acceleration are determined. A new cavitation model is developed and validated against the experimental results. The new model improves the prediction of the void fraction in different regions of the flow field, making it feasible to simulate different regimes of cavitating flow. Finally, the dominant mechanism governing the liquid-vapor transition and the transport of the void fraction is described.

Hunting is a very common instability exhibited by rail vehicles operating at high speeds. The hunting phenomenon is a self excited lateral oscillation that is produced by the forward speed of the vehicle and the wheel-rail interactive forces that result from the conicity of the wheel-rail contours and the friction-creep characteristics of the wheel-rail contact geometry. Hunting can lead to severe ride discomfort and eventual physical damage to wheels and rails. A comprehensive study of the lateral stability of a single wheelset, a single truck, and the complete rail vehicle has been performed. This study investigates bifurcation phenomenon and limit cycles in rail vehicle dynamics. Sensitivity of the critical hunting velocity to various primary and secondary stiffness and damping parameters has been examined. This research assumes the rail vehicle to be moving on a smooth, level, and tangential track, and all parts of the rail vehicle to be rigid. Sources of nonlinearities in the rail vehicle model are the nonlinear wheel-rail profile, the friction-creep characteristics of the wheel-rail contact geometry, and the nonlinear vehicle suspension characteristics. This work takes both single-point and two-point wheel-rail contact conditions into account. The results of the lateral stability study indicate that the critical velocity of the rail vehicle is most sensitive to the primary longitudinal stiffness. A method has been developed to eliminate hunting behavior in rail vehicles by increasing the critical velocity of hunting beyond the operational speed range. This method involves the semi-active control of the primary longitudinal stiffness using the wheelset yaw displacement. This approach is seen to considerably increase the critical hunting velocity. / Master of Science

We study fast-slow versions of the SIR, SIRS and SIRWS epidemiological models, and of the SIRS epidemiological model on homogeneous graphs, obtained through the application of the moment closure method. The multiple time scale behavior is introduced to account for large differences between some of the rates of the epidemiological pathways. Our main purpose is to show that the fast-slow models, even though in nonstandard form, can be studied by means of Geometric Singular Perturbation Theory (GSPT). In particular, without using Lyapunov's method, we are able to not only analyze the stability of the endemic equilibria of the SIR and SIRS models, but also to show that in the remaining models limit cycles arise. We show that the proposed approach is particularly useful in more complicated (higher dimensional) models such as the SIRWS model and the SIRS on homogeneous graphs, for which we provide a detailed description of their dynamics by combining analytic and numerical techniques. In particular, for the latter we show that the model can give rise to periodic solutions, differently from the corresponding model based on homogeneous mixing.