Sensitivity analysis quantifies the effect that of perturbations of the model inputs have on the model's outputs. Some of the key insights gained using sensitivity analysis are to understand the robustness of the model with respect to perturbations, and to select the most important parameters for the model. MATLODE is a tool for sensitivity analysis of models described by ordinary differential equations (ODEs). MATLODE implements two distinct approaches for sensitivity analysis: direct (via the tangent linear model) and adjoint. Within each approach, four families of numerical methods are implemented, namely explicit Runge-Kutta, implicit Runge-Kutta, Rosenbrock, and single diagonally implicit Runge-Kutta. Each approach and family has its own strengths and weaknesses when applied to real world problems. MATLODE has a multitude of options that allows users to find the best approach for a wide range of initial value problems. In spite of the great importance of sensitivity analysis for models governed by differential equations, until this work there was no MATLAB ordinary differential equation sensitivity analysis toolbox publicly available. The two most popular sensitivity analysis packages, CVODES [8] and FATODE [10], are geared toward the high performance modeling space; however, no native MATLAB toolbox was available. MATLODE fills this need and offers sensitivity analysis capabilities in MATLAB, one of the most popular programming languages within scientific communities such as chemistry, biology, ecology, and oceanogra- phy. We expect that MATLODE will prove to be a useful tool for these communities to help facilitate their research and fill the gap between theory and practice. / Master of Science / Sensitivity analysis is the study of how small changes in a model?s input effect the model’s output. Sensitivity analysis provides tools to quantify the impact that small, discrete changes in input values have on the output. The objective of this research is to develop a MATLAB sensitivity analysis toolbox called MATLODE. This research is critical to a wide range of communities who need to optimize system behavior or predict outcomes based on a variety of initial conditions. For example, an analyst could build a model that reflects the performance of an automobile engine, where each part in the engine has a set of initial characteristics. The analyst can use sensitivity analysis to determine which part effects the engine’s overall performance the most (or the least), without physically building the engine and running a series of empirical tests. By employing sensitivity analysis, the analyst saves time and money, and since multiple tests can usually be run through the model in the time needed to run just one empirical test, the analyst is likely to gain deeper insight and design a better product. Prior to MATLODE, employing sensitivity analysis without significant knowledge of computational science was too cumbersome and essentially impractical for many of the communities who could benefit from its use. MATLODE bridges the gap between computational science and a variety of communities faced with understanding how small changes in a system’s input values effect the systems output; and by bridging that gap, MATLODE enables more large scale research initiatives than ever before.

Infantile nystagmus is a condition that causes involuntary, bilateral and conjugate oscillations of the eyes, which are predominately restricted to the horizontal plane. In order to investigate the cause of nystagmus, computational models and nonlinear dynamics techniques have been used to model and analyse the oculomotor system. Computational models are important in making predictions and creating a quantitative framework for the analysis of the oculomotor system. Parameter estimation is a critical step in the construction and analysis of these models. A preliminary parameter estimation of a nonlinear dynamics model proposed by Broomhead et al. [1] has been shown to be able to simulate both normal rapid eye movements (i.e. saccades) and nystagmus oscillations. The application of nonlinear analysis to experimental jerk nystagmus recordings, has shown that the local dimensions number of the oscillation varies across the phase angle of the nystagmus cycle. It has been hypothesised that this is due to the impact of signal dependent noise (SDN) on the neural commands in the oculomotor system. The main aims of this study were: (i) to develop parameter estimation methods for the Broomhead et al. [1] model in order to explore its predictive capacity by fitting it to experimental recordings of nystagmus waveforms and saccades; (ii) to develop a stochastic oculomotor model and examine the hypothesis that noise on the neural commands could be the cause of the behavioural characteristics measured from experimental nystagmus time series using nonlinear analysis techniques. In this work, two parameter estimation methods were developed, one for fitting the model to the experimental nystagmus waveforms and one to saccades. By using the former method, we successfully fitted the model to experimental nystagmus waveforms. This fit allowed to find the specific parameter values that set the model to generate these waveforms. The types of the waveforms that we successfully fitted were asymmetric pseudo-cycloid, jerk and jerk with extended foveation. The fit of other types of nystagmus waveforms were not examined in this work. Moreover, the results showed which waveforms the model can generate almost perfectly and the waveform characteristics of a number of jerk waveforms which it cannot exactly generate. These characteristics were on a specific type of jerk nystagmus waveforms with a very extreme fast phase. The latter parameter estimation method allowed us to explore whether the model can generate horizontal saccades of different amplitudes with the same behaviour as observed experimentally. The results suggest that the model can generate the experimental saccadic velocity profiles of different saccadic amplitudes. However, the results show that best fittings of the model to the experimental data are when different model parameter values were used for different saccadic amplitude. Our parameter estimation methods are based on multi-objective genetic algorithms (MOGA), which have the advantage of optimising biological models with a multi-objective, high-dimensional and complex search space. However, the integration of these models, for a wide range of parameter combinations, is very computationally intensive for a single central processing unit (CPU). To overcome this obstacle, we accelerated the parameter estimation method by utilising the parallel capabilities of a graphics processing unit (GPU). Depending of the GPU model, this could provide a speedup of 30 compared to a midrange CPU. The stochastic model that we developed is based on the Broomhead et al. [1] model, with signal dependent noise (SDN) and constant noise (CN) added to the neural commands. We fitted the stochastic model to saccades and jerk nystagmus waveforms. It was found that SDN and CN can cause similar variability to the local dimensions number of the oscillation as found in the experimental jerk nystagmus waveforms and in the case of saccade generation the saccadic variability recorded experimentally. However, there are small differences in the simulated behaviour compared to the nystagmus experimental data. We hypothesise that these could be caused by the inability of the model to simulate exactly key jerk waveform characteristics. Moreover, the differences between the simulations and the experimental nystagmus waveforms indicate that the proposed model requires further expansion, and this could include other oculomotor subsystem(s).

Conventional coupling paradigms used nowadays to couple various physics components in reactor analysis problems can be inconsistent in their treatment of the nonlinear terms. This leads to usage of smaller time steps to maintain stability and accuracy requirements thereby increasing the computational time. These inconsistencies can be overcome using better approximations to the nonlinear operator in a time stepping strategy to regain the lost accuracy. This research aims at finding remedies that provide consistent coupling and time stepping strategies with good stability properties and higher orders of accuracy. Consistent coupling strategies, namely predictive and accelerated methods, were introduced for several reactor transient accident problems and the performance was analyzed for a 0-D and 1-D model. The results indicate that consistent approximations can be made to enhance the overall accuracy in conventional codes with such simple nonintrusive techniques. A detailed analysis of a monoblock coupling strategy using time adaptation was also implemented for several higher order Implicit Runge-Kutta (IRK) schemes. The conclusion from the results indicate that adaptive time stepping provided better accuracy and reliability in the solution fields than constant stepping methods even during discontinuities in the transients. Also, the computational and the total memory requirements for such schemes make them attractive alternatives to be used for conventional coupling codes.