Image Reconstruction from Undersampled Confocal Microscopy Data using Multiresolution Based Maximum Entropy Regularization

We consider the problem of reconstructing 2D images from randomly under-sampled confocal microscopy samples. The well known and widely celebrated total variation regularization, which is the L1 norm of derivatives, turns out to be unsuitable for this problem; it is unable to handle both noise and under-sampling together. This issue is linked with the notion of phase transition phenomenon observed in compressive sensing research, which is essentially the break-down of total variation methods, when sampling density gets lower than certain threshold. The severity of this breakdown is determined by the so-called mutual incoherence between the derivative operators and measurement operator. In our problem, the mutual incoherence is low, and hence the total variation regularization gives serious artifacts in the presence of noise even when the sampling density is not very low. There has been very few attempts in developing regularization methods that perform better than total variation regularization for this problem. We develop a multi-resolution based regularization method that is adaptive to image structure. In our approach, the desired reconstruction is formulated as a series of coarse-to-fine multi-resolution reconstructions; for reconstruction at each level, the regularization is constructed to be adaptive to the image structure, where the information for adaption is obtained from the reconstruction obtained at coarser resolution level. This adaptation is achieved by using maximum entropy principle, where the required adaptive regularization is determined as the maximizer of entropy subject to the information extracted from the coarse reconstruction as constraints. We demonstrate the superiority of the proposed regularization method over existing ones using several reconstruction examples

Investigation of flat capacitor discharge electromagnetic launchers

In this thesis, studies of flat or pancake type electromagnetic launcher systems are described. The studies involved the development of several numerical models, and are supported throughout by experimental investigation. The models were based on a coaxial filamentary division technique, and the results they provided were compared with those from a commercial electromagnetic finite element modelling package. They were used to investigate some of the many possible launcher structures and power supply arrangements, as part of a wide-ranging parametric study. The aim of this thesis was to gain an insight into the factors that affect the performance of the launchers. Several different techniques were implemented to reduce the computation time. Practical experimentation provided a clear demonstration of the launcher technology, and supplied valuable model validation data. To aid the experimental work new projectile speed and yaw measurement systems were developed, and these were supported by results from a high-speed camera. A novel dual projectile launcher was tested, and was shown to improve the launch efficiency and to operate at higher energies, due to the reduction in drive coil recoil. Projectile deformation was investigated in both solid discs and flat annular projectiles

CAD-Based Porous Scaffold Design of Intervertebral Discs in Tissue Engineering

With the development and maturity of three-dimensional (3D) printing technology over the past decade, 3D printing has been widely investigated and applied in the field of tissue engineering to repair damaged tissues or organs, such as muscles, skin, and bones, Although a number of automated fabrication methods have been developed to create superior bio-scaffolds with specific surface properties and porosity, the major challenges still focus on how to fabricate 3D natural biodegradable scaffolds that have tailor properties such as intricate architecture, porosity, and interconnectivity in order to provide the needed structural integrity, strength, transport, and ideal microenvironment for cell- and tissue-growth. In this dissertation, a robust pipeline of fabricating bio-functional porous scaffolds of intervertebral discs based on different innovative porous design methodologies is illustrated. Firstly, a triply periodic minimal surface (TPMS) based parameterization method, which has overcome the integrity problem of traditional TPMS method, is presented in Chapter 3. Then, an implicit surface modeling (ISM) approach using tetrahedral implicit surface (TIS) is demonstrated and compared with the TPMS method in Chapter 4. In Chapter 5, we present an advanced porous design method with higher flexibility using anisotropic radial basis function (ARBF) and volumetric meshes. Based on all these advanced porous design methods, the 3D model of a bio-functional porous intervertebral disc scaffold can be easily designed and its physical model can also be manufactured through 3D printing. However, due to the unique shape of each intervertebral disc and the intricate topological relationship between the intervertebral discs and the spine, the accurate localization and segmentation of dysfunctional discs are regarded as another obstacle to fabricating porous 3D disc models. To that end, we discuss in Chapter 6 a segmentation technique of intervertebral discs from CT-scanned medical images by using deep convolutional neural networks. Additionally, some examples of applying different porous designs on the segmented intervertebral disc models are demonstrated in Chapter 6

A high-performance boundary element method and its applications in engineering

As a semi-numerical and semi-analytical method, owing to the inherent advantage, of boundary-only discretisation, the boundary element method (BEM) has been widely applied to problems with complicated geometries, stress concentration problems, infinite domain problems, and many others. However, domain integrals and non-symmetrical and dense matrix systems are two obstacles for BEM which have hindered the its further development and application. This thesis is aimed at proposing a high-performance BEM to tackle the above two drawbacks and broaden the application scope of BEM. In this thesis, a detailed introduction to the traditional BEM is given and several popular algorithms are introduced or proposed to enhance the performance of BEM. Numerical examples in heat conduction analysis, thermoelastic analysis and thermoelastic fracture problems are performed to assess the efficiency and correction of the algorithms. In addition, necessary theoretical derivations are embraced for establishing novel boundary integral equations (BIEs) for specific engineering problems. The following three parts are the main content of this thesis. (1) The first part (Part II consisting of two chapters) is aimed at heat conduction analysis by BEM. The coefficient matrix of equations formed by BEM in solving problems is fully-populated which occupy large computer memory. To deal with that, the fast multipole method (FMM) is introduced to energize the line integration boundary element method (LIBEM) to performs better in efficiency. In addition, to compute domain integrals with known or unknown integrand functions which are caused by heat sources or heterogeneity, a novel BEM, the adaptive orthogonal interpolation moving least squares (AOIMLS) method enhanced LIBEM, which also inherits the advantage of boundary-only discretisation, is proposed. Unlike LIBEM, which is an accurate and stable method for computing domain integrals, but only works when the mathematical expression of integral function in domain integrals is known, the AOIMLS enhanced LIBEM can compute domain integrals with known or unknown integral functions, which ensures all the nonlinear and nonhomogeneous problems can be solved without domain discretisation. In addition, the AOIMLS can adaptively avoid singular or ill-conditioned moment matrices, thus ensuring the stability of the calculation results. (2) In the second part (Part III consisting of four chapters), the thermoelastic problems and fracture problems are the main objectives. Due to considering thermal loads, domain integrals appear in the BIEs of the thermoelastic problems, and the expression of integrand functions is known or not depending on the temperature distribution given or not, the AOIMLS enhanced LIBEM is introduced to conduct thermoelasticity analysis thereby. Besides, a series of novel unified boundary integral equations based on BEM and DDM are derived for solving fracture problems and thermoelastic fracture problems in finite and infinite domains. Two sets of unified BIEs are derived for fracture problems in finite and infinite domains based on the direct BEM and DDM respectively, which can provide accurate and stable results. Another two sets of BIEs are addressed by employing indirect BEM and DDM, which cannot ensure a stable result, thereby a modified indirect BEM is proposed which performs much more stable. Moreover, a set of novel BIEs based on the direct BEM and DDM for cracked domains under thermal stress is proposed. (3) In the third part (Part IV consisting of one chapter), a high-efficiency combined BEM and discrete element method (DEM) is proposed to compute the inner stress distribution and particle breakage of particle assemblies based on the solution mapping scheme. For the stress field computation of particles with similar geometry, a template particle is used as the representative particle, so that only the related coefficient matrices of one template particle in the local coordinate system are needed to be calculated, while the coefficient matrices of the other particles, can be obtained by mapping between the local and global coordinate systems. Thus, the combined BEM and DEM is much more effective when modelling a large-scale particle system with a small number of distinct possible particle shapes. Furthermore, with the help of the Hoek-Brown criterion, the possible cracks or breakage paths of a particle can be obtained

Markerless deformation capture of hoverfly wings using multiple calibrated cameras

This thesis introduces an algorithm for the automated deformation capture of hoverfly wings from multiple camera image sequences. The algorithm is capable of extracting dense surface measurements, without the aid of fiducial markers, over an arbitrary number of wingbeats of hovering flight and requires limited manual initialisation. A novel motion prediction method, called the ‘normalised stroke model’, makes use of the similarity of adjacent wing strokes to predict wing keypoint locations, which are then iteratively refined in a stereo image registration procedure. Outlier removal, wing fitting and further refinement using independently reconstructed boundary points complete the algorithm. It was tested on two hovering data sets, as well as a challenging flight manoeuvre. By comparing the 3-d positions of keypoints extracted from these surfaces with those resulting from manual identification, the accuracy of the algorithm is shown to approach that of a fully manual approach. In particular, half of the algorithm-extracted keypoints were within 0.17mm of manually identified keypoints, approximately equal to the error of the manual identification process. This algorithm is unique among purely image based flapping flight studies in the level of automation it achieves, and its generality would make it applicable to wing tracking of other insects

A high-order integral solver for scalar problems of diffraction by screens and apertures in three-dimensional space

We present a novel methodology for the numerical solution of problems of diffraction by infinitely thin screens in three-dimensional space. Our approach relies on new integral formulations as well as associated high-order quadrature rules. The new integral formulations involve weighted versions of the classical integral operators related to the thin-screen Dirichlet and Neumann problems as well as a generalization to the open-surface problem of the classical Calderón formulae. The high-order quadrature rules we introduce for these operators, in turn, resolve the multiple Green function and edge singularities (which occur at arbitrarily close distances from each other, and which include weakly singular as well as hypersingular kernels) and thus give rise to super-algebraically fast convergence as the discretization sizes are increased. When used in conjunction with Krylov-subspace linear algebra solvers such as GMRES, the resulting solvers produce results of high accuracy in small numbers of iterations for low and high frequencies alike. We demonstrate our methodology with a variety of numerical results for screen and aperture problems at high frequencies—including simulation of classical experiments such as the diffraction by a circular disc (featuring in particular the famous Poisson spot), evaluation of interference fringes resulting from diffraction across two nearby circular apertures, as well as solution of problems of scattering by more complex geometries consisting of multiple scatterers and cavities

Laser doppler vibrometer for efficient structural health monitoring

The research effort in this thesis is devoted to develop techniques to accurately and rapidly identify the location, orientation, and magnitude of the defects by using structural health monitoring concepts that use Laser Doppler Vibrometer as a non-contact sensor with multi-point sensing capability. The first research area addresses the formulation and validation of an innovative Damage Measure that is based on the ratios of the strain energy distributions of the damaged and undamaged structure. The innovations include use of a single set of actuator/sensor pair to excite and detect the responses of a structure for low frequency vibrations as well as guided wave propagation studies. A second new capability is the estimation of the Damage Measure without requiring any knowledge of the undamaged baseline structure. This method is made possible because of the development of these new technologies: Spatial Decimation and Wavenumber/Frequency filtering. The third contribution is to develop analytical models for the structural dynamics of damaged structure and seek solutions that use perturbation methods to detect damage in a plate structure. The fourth contribution is the development of a comprehensive damage detection technique over a wide frequency dynamic range. The fifth topic of research involves automation in Structural Health Monitoring based on the comprehensive Damage Measure formulation. Under the control of software the Scanning Laser Doppler Vibrometer is used to acquire the low frequency vibration mode data for a coarse identification of all the suspect regions of damage using a threshold criterion on the Damage Measure. Each suspect region of damage is further investigated using the high frequency elastic wave propagation to clearly identify the location, orientation, and extent of the damage. The computer control of the Laser Doppler Vibrometer and a quantitative assessment of the damage provide the enabling technologies for the automation proof of concept. Finally the developed techniques of damage detection are successfully demonstrated on practical structures such as a turbine blade in the laboratory and an F-15 vertical tail in field maintenance conditionsPh.D.Committee Chair: Hanagud, Sathya; Committee Member: Apetre, Nicole; Committee Member: Engelstad, Steve; Committee Member: Glass, Brian; Committee Member: Kardomateas, George; Committee Member: Ruzzene, Massim