The direct and adjoint methods are to linearize the time-averaged solution of bounded dynamical systems about one or more design parameters. Hence, such methods are one way to obtain the gradient necessary in locally optimizing a dynamical system’s time-averaged behavior over those design parameters. However, when analyzing nonlinear systems whose solutions exhibit chaos, standard direct and adjoint sensitivity methods yield meaningless results due to time-local instability of the system. The present work proposes a new method of solving the direct and adjoint linear systems in time, then tests that method’s ability to solve instances of the Lorenz system that exhibit chaotic behavior. Promising results emerge and are presented in the form of a regression analysis across a parametric study of the Lorenz system.

This thesis presents the development towards a system that can capture and quantify motion for applications in biomechanical and medical fields demanding precision motion tracking using the lighthouse technology. Commercially known as SteamVR tracking, the lighthouse technology is a motion tracking system developed for virtual reality applications that makes use of patterned infrared light sources to highlight trackers (objects embedded with photodiodes) to obtain their pose or spatial position and orientation. Current motion capture systems such as the camera-based motion capture are expensive and not readily available outside of research labs. This thesis provides a case for low-cost motion capture systems. The technology is applied to quantify motion to draw inferences about biomechanics capture and analysis, quantification of gait, and prosthetic alignment. Possible shortcomings for data acquisition using this system for the stated applications have been addressed. The repeatability of the system has been established by determining the standard deviation error for multiple trials based on a motion trajectory using a seven degree-of-freedom robot arm. The accuracy testing for the system is based on cross-validation between the lighthouse technology data and transformations derived using joint angles by developing a forward kinematics model for the robot’s end-effector pose. The underlying principle for motion capture using this system is that multiple trackers placed on limb segments allow to record the position and orientation of the segments in relation to a set global frame. Joint angles between the segments can then be calculated from the recorded positions and orientations of each tracker using inverse kinematics. In this work, inverse kinematics for rigid bodies was based on calculating homogeneous transforms to the individual trackers in the model’s reference frame to find the respective Euler angles as well as using the analytical approach to solve for joint variables in terms of known geometric parameters. This work was carried out on a phantom prosthetic limb. A custom application-specific motion tracker was also developed using a hardware development kit which would be further optimized for subsequent studies involving biomechanics motion capture.

The author presents some generalized Burgers' equations for incompressible and isothermal flow of viscous non-Newtonian fluids based on the Cross model, the Carreau model, and the Power-Law model and some simple assumptions on the flows. The author numerically solves the traveling wave equations for the Cross model, the Carreau model, the Power-Law model by using industrial data. The author proves existence and uniqueness of solutions to the traveling wave equations of each of the three models. The author also provides numerical estimates of the shock thickness as well as maximum strain $\varepsilon_{11}$ for each of the fluids.

In this work, we design a linear, two step implicit finite difference method to approximate the solutions of a biological system that describes the interaction between a microbial colony and a surrounding substrate. Three separate models are analyzed, all of which can be described as systems of partial differential equations (PDE)s with nonlinear diffusion and reaction, where the biological colony grows and decays based on the substrate bioavailability. The systems under investigation are all complex models describing the dynamics of biological films. In view of the difficulties to calculate analytical solutions of the models, we design here a numerical technique to consistently approximate the system evolution dynamics, guaranteeing that nonnegative initial conditions will evolve uniquely into new, nonnegative approximations. This property of our technique is established using the theory of M-matrices, which are nonsingular matrices where all the entries of their inverses are positive numbers. We provide numerical simulations to evince the preservation of the nonnegative character of solutions under homogeneous Dirichlet and Neumann boundary conditions. The computational results suggest that the method proposed in this work is stable, and that it also preserves the bounded character of the discrete solutions.

Several targeting algorithms are developed and analyzed for possible future use onboard a spacecraft. Each targeter is designed to determine the appropriate propulsive burn for translunar injection to obtain desired orbital parameters upon arrival at the moon. Primary design objectives are to minimize the computational requirements for each algorithm but also to ensure reasonable accuracy, so that the algorithm’s errors do not force the craft to conduct large mid-course corrections. Several levels of accuracy for dynamical models are explored, the convergence range and speed of each algorithm are compared, and the possible benefits of the Broyden and trust-region targeters are evaluated. These targeters provide a proof of concept for the feasibility of a translunar injection targeting algorithm. Anticipating some future improvements, these algorithms could serve as a viable alternative to uploading ground-based targeting solutions and bypass the problems of delays and disruptions in communication, enabling the craft to conduct a translunar injection burn autonomously.

In this thesis the Ramberg-Osgood nonlinear model for describing the behavior of many different materials is investigated. A brief overview of the model as it is currently used in the literature is undertaken and several misunderstandings and possible pitfalls in its application is pointed out, especially as it pertains to more recent approaches to finding solutions involving the model. There is an investigation of the displacement of a cantilever beam under a combined loading consisting of a distributed load across the entire length of the beam and a point load at its end and new solutions to this problem are provided with a mixture of numerical techniques, which suggest strong mathematical consistency within the model for all theoretical assumptions made. A physical experiment was undertaken and the results prove to be inaccurate when using parameters derived from tensile tests, but when back calculating parameters from the beam test the model has a 14.40% error at its extreme against the experimental data suggesting the necessity for further testing.

This dissertation focused on developing a numerical model of spatial and temporal changes in bed morphology of ridge and slough features in wetlands with respect to hydrology and sediment transport when a sudden change in hydrologic condition occurs. The specific objectives of this research were: (1) developing a two-dimensional hydrology model to simulate the spatial distribution of flow depth and velocity over time when a pulsed flow condition is applied, (2) developing a process-based numerical model of sediment transport coupled with flow depth and velocity in wetland ecosystems, and (3) use the developed model to explore how sediment transport may affect the changes in bed elevation of ridge and slough landscape patterns observed in wetlands when a conditional pulsed flow was applied. The results revealed the areas within deep sloughs where flow velocities and directions change continuously. This caused enhanced mixing areas within the deep slough. These mixing areas may have had the potential to affect processes such as sediment redistribution and nutrient transport. The simulation results of solute/sediment transport model also supported the existence of areas within the domain where the mixing processes happened. These areas may have caused that nutrients and suspended particles stay longer time rather than entraining toward downstream and exiting the system. The results of bed simulation have shown very small magnitude of change in bed elevation inside deep slough and no changes on the ridge portion of the study area, when a conditional pulsed flow is applied. These findings may suggest that implementing pulsed flow condition did not increase suspended sediment concentration, which results in insignificant changes in bed morphology of a ridge and slough landscape. Therefore sediment transport may not play an important role in wetland bed morphology and ridge and slough stability. Results from the model development and numerical simulations from this research will provide an improved understanding of how wetland features such as ridge may have formed and degraded by changes in water management that resulted from increasing human activity in wetlands such as The Florida Everglades, over the past decades.

Congestion on roadways and high level of uncertainty of traffic conditions are major considerations for trip planning. The purpose of this research is to investigate the characteristics and patterns of spatial and temporal correlations and also to detect other variables that affect correlation in a freeway setting. 5-minute speed aggregates from the Performance Measurement System (PeMS) database are obtained for two directions of an urban freeway – I-10 between Santa Monica and Los Angeles, California. Observations are for all non-holiday weekdays between January 1st and June 30th, 2010. Other variables include traffic flow, ramp locations, number of lanes and the level of congestion at each detector station. A weighted least squares multilinear regression model is fitted to the data; the dependent variable is Fisher Z transform of correlation coefficient. Estimated coefficients of the general regression model indicate that increasing spatial and temporal distances reduces correlations. The positive parameters of spatial and temporal distance interaction term show that the reduction rate diminishes with spatial or temporal distance. Higher congestion tends to retain higher expected value of correlation; corrections to the model due to variations in road geometry tend to be minor. The general model provides a framework for building a family of more responsive and better-fitting models for a 6.5 mile segment of the freeway during three times of day: morning, midday, and afternoon. Each model is cross-validated on two locations: the opposite direction of the freeway, and a different location on the direction used for estimation. Cross-validation results show that models are able to retain 75% or more of their original predictive capability on independent samples. Incorporation of predictor variables that describe road geometry and traffic conditions into the model works beneficially in capturing a significant portion of variance of the response. The developed regression models are thus transferrable and are apt to predict correlation on other freeway locations.

<p>In this study, we concentrate on spin-1/2 systems. A series of tools using the Liouville space method have been developed for simulating of NMR of arbitrary pulse sequences.</p> <p>We have calculated one- and two-spin symbolically, and larger systems numerically of steady states. The one-spin calculations show how SSFP converges to continuous wave NMR. A general formula for two-spin systems has been derived for the creation of double-quantum signals as a function of irradiation strength, coupling constant, and chemical shift difference. The formalism is general and can be extended to more complex spin systems.</p> <p>Estimates of transverse relaxation, R₂, are affected by frequency offset and field inhomogeneity. We find that in the presence of expected B₀ inhomogeneity, off-resonance effects can be removed from R₂ measurements, when ||omega||<= 0.5 gamma\,B₁ in Hahn echo experiments, when ||omega||<=gamma\,B₁ in CPMG experiments with specific phase variations, by fitting exact solutions of the Bloch equations given in the Lagrange form.</p> <p>Approximate solutions of CPMG experiments show the specific phase variations can significantly smooth the dependence of measured intensities on frequency offset in the range of +/- 1/2 gamma\,B₁. The effective R₂ of CPMG experiments when using a phase variation scheme can be expressed as a second-order formula with respect to the ratio of offset to pi-pulse amplitude.</p> <p>Optimization problems using the exact or approximate solution of the Bloch equations are established for designing optimal broadband universal rotation (OBUR) pulses. OBUR pulses are independent of initial magnetization and can be applied to replace any pulse of the same flip angles in a pulse sequence. We demonstrate the process to exactly and efficiently calculate the first- and second-order derivatives with respect to pulses. Using these exact derivatives, a second-order optimization method is employed to design pulses. Experiments and simulations show that OBUR pulses can provide more uniform spectra in the designed offset range and come up with advantages in CPMG experiments.</p> / Doctor of Philosophy (PhD)

<p>In this thesis we study several inverse sampling procedures to test for homogeneity in a multivariate hypergeometric distribution. The procedures are finite population analogues of the procedures introduced in Panchapakesan et al. (1998) for the multinomial distribution. In order to develop some exact calculations for critical values not considered in Panchapakesan et al. we introduce some terminologies for target probabilities, transfer probabilities, potential target points, right intersection, and left union. Under the null and the alternative hypotheses, we give theorems to calculate the target and transfer probabilities, we then use these results to develop exact calculations for the critical values and powers of one of the procedures. We also propose a new approximate calculation. In order to speed up some of the calculations, we propose several fast algorithms for multiple summation.</p> <p>N >= 1680000, all the results are the same as those in the multinomial distribution.</p> <p>The computing results showed that the simulations agree closely with the exact results. For small population sizes the critical values and powers of the procedures are different from the corresponding multinomial procedures, but when</p> / Master of Science (MSc)