Control of bulk melt crystal growth techniques is desirable for producing semiconductors with the highest purity and ternary alloys with tunable electrical properties. Because these molten materials are electrically conducting, external magnetic fields are often employed to regulate the flow in the melt. However, complicated by the coupled flow, thermal, electromagnetic and chemical physics, such magnetic control is typically empirical or even an educated guess. Two magnetic flow control mechanisms: flow damping by steady magnetic fields, and flow stirring by alternating magnetic fields, are investigated numerically. Magnetic damping during optically-heated float-zone crystal growth is modeled using a spectral collocation method. The Marangoni convection at the free melt-gas interface is suppressed by applying a steady magnetic field, measured by the Hartmann number Ha. Using normal mode linear stability analyses, suppression of detrimental flow instabilities is quantitatively determined in a range applicable to experiments (up to Ha = 300 for Pr = 0.02, and up to Ha = 500 for Pr = 0.001). The hydrodynamic flow instability for small Prandtl number P r float-zone is confirmed by energy analyses. Rotating magnetic field stirring during confined crystal growth in an ampoule is also modeled. Decoupled from the flow field at small magnetic Reynolds number, the electromagnetic field is solved in a finite element solver. At low AC frequencies, the force is only in the azimuthal direction but penetrates deep into the melt. In contrast, the magnetic shielding effect is observed at high alternating current (AC) frequencies, where the external magnetic field penetrates only by a skin depth into the electrically conducting media within the short AC cycle. As a result, the electromagnetic body force is primarily confined to the ampoule surface. At these high AC frequencies the magnetic flux lines are drastically distorted within the melt. The body force is fully three-dimensional and is much stronger than at low AC frequencies, but is confined to near the ampoule surface due to the magnetic shielding effect. These models promote fundamental understanding of flow dynamics regulated by electromagnetic body forces. They provide quantitative guidance for crystal growth to minimize trial and error experimentation that is slow and expensive.

Linearized stability analysis methodologies, system identification algorithms and optimal control approaches that are applicable to large scale, flexible multibody dynamic systems are presented in this thesis. For stability analysis, two classes of closely related algorithms based on a partial Floquet approach and on an autoregressive approach, respectively, are presented in a common framework that underlines their similarity and their relationship to other methods. The robustness of the proposed approach is improved by using optimized signals that are derived from the proper orthogonal modes of the system. Finally, a signal synthesis procedure based on the identified frequencies and damping rates is shown to be an important tool for assessing the accuracy of the identified parameters; furthermore, it provides a means of resolving the frequency indeterminacy associated with the eigenvalues of the transition matrix for periodic systems. For system identification, a robust algorithm is developed to construct subspace plant models. This algorithm uniquely combines the methods of minimum realization and subspace identification. It bypasses the computation of Markov parameters because the free impulse response of the system can be directly computed in the present computational environment. Minimum realization concepts were applied to identify the stability and output matrices. On the other hand, subspace identification algorithms construct a state space plant model of linear system by using computationally expensive oblique matrix projection operations. The proposed algorithm avoids this burden by computing the Kalman filter gain matrix and model dependency on external inputs in a small sized subspace. Balanced model truncation and similarity transformation form the theoretical foundation of proposed algorithm. Finally, a forward innovation model is constructed and estimates the input-output behavior of the system within a specified level of accuracy. The proposed system identification algorithms are computationally inexpensive and consist of purely post processing steps that can be used with any multi-physics computational tool or with experimental data. Optimal control methodologies that are applicable to comprehensive large-scale flexible multibody systems are presented. A classical linear quadratic Gaussian controller is designed, including subspace plant identification, the evaluation of linear quadratic regulator feedback gain and Kalman filter gain matrices and online control implementation.

The objective of this thesis is to investigate the transient response and optimal economic dispatch of an isolated power system with wind generators. Different types of wind turbines and the classification of Stability are introduced. Then, the process of Transient stability analysis and the concept of Genetic Algorithms are given for explanation. In this thesis, the practical power system of Kinmen is selected for case study. The disturbances introduced by gusting wind and N-1 system contingency are considered in the transient stability analysis. Furthermore, in order to obtain both accuracy and feasibility of the Optimal power dispatch by using Real-parameter Genetic Algorithms, the simulation results should be tested for the restrictions and requirements of the actual operation.

For linear system Ax = b, the traditional condition number is the worst case for all b¡¦s and often overestimated in many problems. For a specific b, the effective condition number is a better upper bound for the relative error of x. But, it is also possible that this effective condition number is overestimated. In this thesis, we study the true ratio of the relative error of x to the relative perturbation of b, called the true condition number. We obtain several new upper bounds and estimates for true condition number. We also explore to change the system to an equivalent one by shifting b to minimize its effective condition number. Finally we apply all our results to functional approximation.

Thermoacoustic tomography (TAT) is a newly emerging modality in biomedical imaging. It combines the good contrast of electromagnetic and good resolution of ultrasound imaging. The mathematical model of TAT is the observability problem for the wave equation: one observes the data on a hyper-surface and reconstructs the initial perturbation. In this dissertation, we consider several mathematical problems of TAT. The first problem is the inversion formulas. We provide a family of closed form inversion formulas to reconstruct the initial perturbation from the observed data. The second problem is the range description. We present the range description of the spherical mean Radon transform, which is an important transform in TAT. The next problem is the stability analysis for TAT. We prove that the reconstruction of the initial perturbation from observed data is not H¨older stable if some observability condition is violated. The last problem is the speed determination. The question is whether the observed data uniquely determines the ultrasound speed and initial perturbation. We provide some initial results on this issue. They include the unique determination of the unknown constant speed, a weak local uniqueness, a characterization of the non-uniqueness, and a characterization of the kernel of the linearized operator.

In the thesis, the error and stability analysis is made for the 2D Helmholtz equation by the method of fundamental solutions (MFS) using both Bessel and Neumann functions. The bounds of errors in bounded simply-connected domains are derived, while the bounds of condition number are derived only for disk domains. The MFS using Bessel functions is more efficient than the MFS using Neumann functions. Interestingly, for the MFS using Bessel functions, the radius R of the source points is not necessarily larger than the maximal radius r_max of the solution domain. This is against the traditional condition: r_max < R for MFS. Numerical experiments are carried out to support the analysis and conclusions made.

In this thesis, the analysis of the method of fundamental solution(MFS) is expanded for biharmonic equations. The bounds of errors are derived for the traditional and the Almansi's approaches in bounded simply-connected domains. The exponential and the polynomial convergence rates are obtained from highly and finite smooth solutions, respectively. Also the bounds of condition number are derived for the disk domains, to show the exponential growth rates. The analysis in this thesis is the first time to provide the rigor analysis of the CTM for biharmonic equations, and the intrinsic nature of accuracy and stability is similar to that of Laplace's equation. Numerical experiment are carried out for both smooth and singularity problems. The numerical results coincide with the theoretical analysis made. When the particular solutions satisfying the biharmonic equation can be found, the method of particular solutions(MPS) is always superior to MFS, supported by numerical examples. However, if such singular particular solutions near the singular points can not be found, the local refinement of collocation nodes and the greedy adaptive techniques can be used. It seems that the greedy adaptive techniques may provide a better solution for singularity problems. Beside, the numerical solutions by Almansi's approaches are slightly better in accuracy and stability than those by the traditional FS. Hence, the MFS with Almansi's approaches is recommended, due to the simple analysis, which can be obtained directly from the analysis of MFS for Laplace's equation.

In this thesis, we extend the spectral collocation methods(SCM) (i.e., pseudo-spectral method) in Quarteroni and Valli [27] for the semilinear, parameter-dependentproblems(PDP) in the square with the Dirichlet boundary condition. The optimal error bounds are derived in this thesis for both H1 and L2 norms. For the solutions sufficiently smooth, the very high convergence rates can be obtained. The algorithms of the SCM are simple and easy to carry out. Only a few of basis functions are needed so that not only can the high accuracy of the PDP solutions be achieved, but also a great deal of CPU time may be saved. Moreover, for PDP the stability analysis of SCM is also made, to have the same growth rates of condition number as those for Poisson¡¦s equation. Numerical experiments are carried out to verify the theoretical analysis made.

Dynamic stall deeply affects the response of helicopter rotor blades, making its modeling accuracy very important. Two commonly used dynamic stall models were implemented in a comprehensive code, validated, and contrasted to provide improved analysis accuracy and versatility. Next, computational fluid dynamics and computational structural dynamics loose coupling methodologies are reviewed, and a general tight coupling approach was implemented and tested. The tightly coupled computational fluid dynamics and computational structural dynamics methodology is then used to assess the stability characteristics of complex rotorcraft problems. An aeroelastic analysis of rotors must include an assessment of potential instabilities and the determination of damping ratios for all modes of interest. If the governing equations of motion of a system can be formulated as linear, ordinary differential equations with constant coefficients, classical stability evaluation methodologies based on the characteristic exponents of the system can rapidly and accurately provide the system's stability characteristics. For systems described by linear, ordinary differential equations with periodic coefficients, Floquet's theory is the preferred approach. While these methods provide excellent results for simplified linear models with a moderate number of degrees of freedom, they become quickly unwieldy as the number of degrees of freedom increases. Therefore, to accurately analyze rotorcraft aeroelastic periodic systems, a fully nonlinear, coupled simulation tool is used to determine the response of the system to perturbations about an equilibrium configuration and determine the presence of instabilities and damping ratios. The stability analysis is undertaken using an algorithm based on a Partial Floquet approach that has been successfully applied with computational structural dynamics tools on rotors and wind turbines. The stability analysis approach is computationally inexpensive and consists of post processing aeroelastic data, which can be used with any aeroelastic rotorcraft code or with experimental data.

This thesis consists of the research on linkage-based fingertips for prosthetic hands. These fingertips consists of small polycentric mechanisms attached to what would be the pulp in normal anatomical fingers. These mechanisms allow the prosthetic hand to conform to the shape of objects during grasp. The goal of these prosthetic fingertips is to maximize the functionality of the hand while minimizing the number of inputs that the user has to control. The stability of the fingertip mechanisms is analyzed using the principle of virtual work. From this analysis we are able to show that the fingertip mechanism is stable for a large range of rotation of the link and for a large range of directions on which the force is applied, and that the mechanism is indifferent to the magnitude of the force applied to it (assuming that the force does not damage/deform the mechanism). To assess if the four-bar mechanisms (fingertips) improve the grasping capabilities in robotics and prosthetics, tests were performed on prosthetic hands and robot grippers with and without the fingertips. Comparisons were made using the Southampton Hand Assessment Procedure (SHAP) protocol, which tests the differences and measures the functionality of particular types of grasp, such as power, spherical, lateral, tripod, tip and extension. In the human testing, the overall Index of Functionality (IOF) of the Hosmer hook is 66.65 and 66.21 for the hook with the fingertips. The hook with the fingertips had a better IOF in the spherical and power prehensile pattern. When the IOF is calculated for the tasks that the fingertips were used, in 10 of 11 of the tasks, the IOF is higher than using the Hosmer hook. In the robotic gripper testing, the Index of Functionality was not be calculated because the time to perform the tasks depended more on the robotic control system than on the physical characteristics of the gripper.