Characterization of multiphase flows integrating X-ray imaging and virtual reality

Multiphase flows are used in a wide variety of industries, from energy production to pharmaceutical manufacturing. However, because of the complexity of the flows and difficulty measuring them, it is challenging to characterize the phenomena inside a multiphase flow. To help overcome this challenge, researchers have used numerous types of noninvasive measurement techniques to record the phenomena that occur inside the flow. One technique that has shown much success is X-ray imaging. While capable of high spatial resolutions, X-ray imaging generally has poor temporal resolution. This research improves the characterization of multiphase flows in three ways. First, an X-ray image intensifier is modified to use a high-speed camera to push the temporal limits of what is possible with current tube source X-ray imaging technology. Using this system, sample flows were imaged at 1000 frames per second without a reduction in spatial resolution. Next, the sensitivity of X-ray computed tomography (CT) measurements to changes in acquisition parameters is analyzed. While in theory CT measurements should be stable over a range of acquisition parameters, previous research has indicated otherwise. The analysis of this sensitivity shows that, while raw CT values are strongly affected by changes to acquisition parameters, if proper calibration techniques are used, acquisition parameters do not significantly influence the results for multiphase flow imaging. Finally, two algorithms are analyzed for their suitability to reconstruct an approximate tomographic slice from only two X-ray projections. These algorithms increase the spatial error in the measurement, as compared to traditional CT; however, they allow for very high temporal resolutions for 3D imaging. The only limit on the speed of this measurement technique is the image intensifier-camera setup, which was shown to be capable of imaging at a rate of at least 1000 FPS. While advances in measurement techniques for multiphase flows are one part of improving multiphase flow characterization, the challenge extends beyond measurement techniques. For improved measurement techniques to be useful, the data must be accessible to scientists in a way that maximizes the comprehension of the phenomena. To this end, this work also presents a system for using the Microsoft Kinect sensor to provide natural, non-contact interaction with multiphase flow data. Furthermore, this system is constructed so that it is trivial to add natural, non-contact interaction to immersive visualization applications. Therefore, multiple visualization applications can be built that are optimized to specific types of data, but all leverage the same natural interaction. Finally, the research is concluded by proposing a system that integrates the improved X-ray measurements, with the Kinect interaction system, and a CAVE automatic virtual environment (CAVE) to present scientists with the multiphase flow measurements in an intuitive and inherently three-dimensional manner

Accurate 3D-reconstruction and -navigation for high-precision minimal-invasive interventions

The current lateral skull base surgery is largely invasive since it requires wide exposure and direct visualization of anatomical landmarks to avoid damaging critical structures. A multi-port approach aiming to reduce such invasiveness has been recently investigated. Thereby three canals are drilled from the skull surface to the surgical region of interest: the first canal for the instrument, the second for the endoscope, and the third for material removal or an additional instrument. The transition to minimal invasive approaches in the lateral skull base surgery requires sub-millimeter accuracy and high outcome predictability, which results in high requirements for the image acquisition as well as for the navigation. Computed tomography (CT) is a non-invasive imaging technique allowing the visualization of the internal patient organs. Planning optimal drill channels based on patient-specific models requires high-accurate three-dimensional (3D) CT images. This thesis focuses on the reconstruction of high quality CT volumes. Therefore, two conventional imaging systems are investigated: spiral CT scanners and C-arm cone-beam CT (CBCT) systems. Spiral CT scanners acquire volumes with typically anisotropic resolution, i.e. the voxel spacing in the slice-selection-direction is larger than the in-the-plane spacing. A new super-resolution reconstruction approach is proposed to recover images with high isotropic resolution from two orthogonal low-resolution CT volumes. C-arm CBCT systems offers CT-like 3D imaging capabilities while being appropriate for interventional suites. A main drawback of these systems is the commonly encountered CT artifacts due to several limitations in the imaging system, such as the mechanical inaccuracies. This thesis contributes new methods to enhance the CBCT reconstruction quality by addressing two main reconstruction artifacts: the misalignment artifacts caused by mechanical inaccuracies, and the metal-artifacts caused by the presence of metal objects in the scanned region. CBCT scanners are appropriate for intra-operative image-guided navigation. For instance, they can be used to control the drill process based on intra-operatively acquired 2D fluoroscopic images. For a successful navigation, accurate estimate of C-arm pose relative to the patient anatomy and the associated surgical plan is required. A new algorithm has been developed to fulfill this task with high-precision. The performance of the introduced methods is demonstrated on simulated and real data

Modeling and Simulation in Engineering

This book provides an open platform to establish and share knowledge developed by scholars, scientists, and engineers from all over the world, about various applications of the modeling and simulation in the design process of products, in various engineering fields. The book consists of 12 chapters arranged in two sections (3D Modeling and Virtual Prototyping), reflecting the multidimensionality of applications related to modeling and simulation. Some of the most recent modeling and simulation techniques, as well as some of the most accurate and sophisticated software in treating complex systems, are applied. All the original contributions in this book are jointed by the basic principle of a successful modeling and simulation process: as complex as necessary, and as simple as possible. The idea is to manipulate the simplifying assumptions in a way that reduces the complexity of the model (in order to make a real-time simulation), but without altering the precision of the results

High Performance Optical Computed Tomography for Accurate Three-Dimensional Radiation Dosimetry

Optical computed tomography (CT) imaging of radiochromic gel dosimeters is a method for truly three-dimensional radiation dosimetry. Although optical CT dosimetry is not widely used currently due to previous concerns with speed and accuracy, the complexity of modern radiotherapy is increasing the need for a true 3D dosimeter. This thesis reports technical improvements that bring the performance of optical CT to a clinically useful level. New scanner designs and improved scanning and reconstruction techniques are described. First, we designed and implemented a new light source for a cone-beam optical CT system which reduced the scatter to primary contribution in CT projection images of gel dosimeters from approximately 25% to approximately 4%. This design, which has been commercially implemented, enables accurate and fast dosimetry. Second, we designed and constructed a new, single-ray, single-detector parallel-beam optical CT scanner. This system was able to very accurately image both absorbing and scattering objects in large volumes (15 cm diameter), agreeing within ∼1% with independent measurements. It has become a reference standard for evaluation of optical CT geometries and dosimeter formulations. Third, we implemented and characterized an iterative reconstruction algorithm for optical CT imaging of gel dosimeters. This improved image quality in optical CT by suppressing the effects of noise and artifacts by a factor of up to 5. Fourth, we applied a fiducial-based ray path measurement scheme, combined with an iterative reconstruction algorithm, to enable optical CT reconstruction in the case of refractive index mismatch between different media in the scanner’s imaged volume. This improved the practicality of optical CT, as time-consuming mixing of liquids can be avoided. Finally, we applied the new laser scanner to the difficult dosimetry task of small-field measurement. We were able to obtain beam profiles and depth dose curves for 4 fields (3x3 cm2 and below) using one 15 cm diameter dosimeter, within 2 hours. Our gel dosimetry depth-dose curves agreed within ∼1.5% with Monte Carlo simulations. In conclusion, the developments reported here have brought optical CT dosimetry to a clinically useful level. Our techniques will be used to assist future research in gel dosimetry and radiotherapy treatment techniques

CM Murchison : nebular formation of fine-grained chondrule rims followed by impact processing on the CM parent body

We examine the primitive carbonaceous chondrite, CM Murchison, to infer details concerning its formation and subsequent processing on the CM parent body. We use X-ray computed tomography (XCT) to measure the 3D morphology and spatial relationship of fine-grained rims (FGRs) of Type I chondrules and find that the relationship between FGR volume and interior core radius is well described by a power law function as proposed for FGR accretion in a turbulent nebula by Cuzzi (2004). We also find evidence that the rimmed chondrules were slightly larger than Kolmogorov-stopping-time nebular particles. Evidence against parent body FGR formation includes a positive correlation between rim thickness and chondrule size and no correlation between interior chondrule roughness (used as a proxy for degree of aqueous alteration) and FGR volume. We find that the chondrules are foliated and that the FGRs are compressed in the direction of maximum stress, resulting in rims that are consistently thicker in the plane of foliation. After accretion to the CM parent body, the material within Murchison experienced brittle deformation, porosity loss, and aqueous alteration. XCT reveals that partially altered chondrules define a prominent foliation and weak lineation. The presence of a lineation and evidence for a component of rotational, noncoaxial shear suggest that the deformation was caused by impact. Olivine optical extinction indicates that the sample is classified as shock stage S1 and electron microscopy reveals that plastic deformation was minimal and that brittle deformation was the dominant microstructural strain-accommodating mechanism. Evidence such as serpentine veins parallel to the foliation fabric and crosscutting alteration veins strongly suggest that some aqueous alteration post-dated or was contemporaneous with the deformation and that multiple episodes of fracturing and mineralization occurred. Finally, using the deformed shape of the chondrules we estimate the strain and infer that the original bulk porosity of Murchison before deformation was 32.2 – 53.4%. Our findings suggest that significant porosity loss, deformation, and compaction from impact can occur on chondrite parent bodies whose samples record only a low level of shock, and that significant chondrule deformation can result from brittle processes and does not require plastic deformation of grains.Geological Science

A Multi-Material Approach to Beam Hardening Correction and Calibration in X-Ray Microtomography

PhDX-ray microtomogaphy is a non-clinical, non-destructive, and quantitative technique for determining three-dimensional mineral concentration distributions in variably radiolucent samples with a spatial resolution on the micron scale. For reasons of practicality, particularly for longterm studies, it is often not possible or desirable to utilise a monochromatic X-ray source and so microtomography using a conventional impact-source X-ray generator to produce a polychromatic photon beam is carried out instead. The use of photons of multiple energies causes the production of projection artefacts arising from preferential absorption, which impair the greyscale accuracy of the resulting reconstruction and the material concentration measurements that are derived from the linear attenuation coefficients (LACs). The purpose of the project described in this thesis is to identify weaknesses in the current method of beam hardening correction and to develop and test a tomographic calibration and projection processing method which may demonstrably improve the quality of current beam hardening correction methods as used with the MuCAT microtomography equipment, which provides a worldclass impact-source microtomography research and production facility at Barts and The London School of Medicine and Dentistry. An overview of the physical basis of X-ray computed tomography and X-ray microtomography is given from first principles, and examples of quantitative applications of the techniques, which generally depend on accurate reconstruction of linear attenuation coefficient values, are discussed. The major sources of artefacts in X-ray microtomography are discussed, particularly those with a direct impact on reconstructed linear attenuation coefficient values. Beam hardening is identified as an error source of particular interest, with secondary research on the effects of any beam hardening correction method on the severity of Compton scatter artefacts, and a critical review is carried out of historical attempts to reduce or mitigate these artefacts, particularly the single-material parameter-optimisation approach in service at the beginning of the research project. A ‘carousel’ test piece comprising multiple attenuators of multiple materials along with attenuation optimisation software based on varying multiple system parameters in order to extend the functionality and usability of the existing correction model, and qualitative results have so far been gathered suggesting the use of this system over the pre-existing attenuation wedge and parameter-optimisation method. A study of the effects of tuning the photon energy to which calibrations are made is carried out, showing improved linear attenuation coefficient recovery at a higher energy than was previously believed to be optimal, and a significant effect arising from X-ray generator target evaporation leading to spatial changes and time-dependence of the target thickness parameter is measured, suggesting that automated calibration as a standard part of the measurement process is required. A stability experiment is carried out using this method in order to examine the possibility of inconsistency resulting from ageing of the filament cathode, which is found not to significantly impact the quality of results. An immersion tank is developed in order to ensure beam hardening correction validity in the case of dual-material specimens where a radiodensity-matching fluid can be provided and the sample is suitable for immersion. Experimental comparison using a commercial beam hardening calibration device as the specimen reveals significantly improved hydroxyapatite concentration measurement recovery. An in-scatter experiment was carried out on the immersion tank, and it was found that there was a significant scatter contribution when the tank was filled in the case where the sample thickness is much less than the tank thickness. Proposals are presented for further work to improve reconstruction quality through of scatter reduction techniques in impactsource microtomographic systems, and to utilise the immersion tank for in situ chemical erosion experiments. The effects of the improvements to the beam hardening process are demonstrated using a biological specimen to demonstrate qualitative changes in reconstruction, particularly in improved dark levels surrounding the specimen. A second experiment is carried out in order to test the reproducibility of results, which is found to be improved by approximately four times over the same dataset corrected using the pre-existingbeam-hardening calibration methodEngineering and Physical Sciences Research Council grant number EP/G007845/1