The staged combustion aft-injected hybrid (SCAIH) rocket motor is a promising design for the future of hybrid rocket propulsion. Advances in computational fluid dynamics and scientific computing have made computational modelling an effective tool in design and development. The focus of this thesis is the numerical modelling of the SCAIH rocket motor in a turbulent combustion, high-speed, reactive flow accounting for solid soot transport and radiative heat transfer. The SCAIH motor has a shear coaxial injector with liquid oxygen injected centrally at sub-critical conditions: 150K, 150m/s (Mach≈0.9), and a gas-generator gas-solid mixture of one-third carbon soot by mass injected in the annual opening at 1175K, and 460m/s (Mach≈0.6). Flow conditions in the near injector region and the flame anchoring mechanism are of particular interest. Overall, the flow is shown to exhibit instabilities and the flame is shown to anchor directly on the injector faceplate with temperatures in excess of 2700K.

Microbubble ultrasound (US) contrast agents have been used to monitor the progression of anti-angiogenic chemotherapies. However, US backscatter measurements used in contrast imaging are inherently variable, given the presence of many microbubbles of random position and size. A model was developed to investigate the influence of US scanner and microbubble characteristics on these variable measurements. The Coefficient of Variation was used to measure variability. It was found that an optimum excitation frequency exists that minimizes this variability. In the case of DefinityTM, a 2.25 MHz centre-frequency pulse yielded a less variable measurement than at 5 MHz. Conversely, decreasing microbubbble concentration was found to significantly increase variability. Evidence suggests that microbubbles are no longer Rayleigh scatterers at sufficient low concentrations. Post-processing was found to aid in reducing measurement variability by averaging samples where microbubble positions are uncorrelated. As well, reduction can be achieved by averaging about a region-of-interest of uniform perfusion.

Many problems in fluid mechanics and material sciences deal with liquid-vapour flows. In these flows, the ideal gas assumption is not accurate and the van der Waals equation of state is usually used. This equation of state is non-convex and causes the solution domain to have two hyperbolic regions separated by an elliptic region. Therefore, the governing equations of these flows have a mixed elliptic-hyperbolic nature. Numerical oscillations usually appear with standard finite-difference space discretization schemes, and they persist when the order of accuracy of the semi-discrete scheme is increased. In this study, we propose to use a Chebyshev pseudospectral method for solving the governing equations. A comparison of the results of this method with very high-order (up to tenth-order accurate) finite difference schemes is presented, which shows that the proposed method leads to a lower level of numerical oscillations than other high-order finite difference schemes, and also does not exhibit fast-traveling packages of short waves which are usually observed in high-order finite difference methods. The proposed method can thus successfully capture various complex regimes of waves and phase transitions in both elliptic and hyperbolic regimes

Longwall Top Coal Caving (LTCC) technology has great potential for more efficient mining of Australian thick coal seams. LTCC can potentially double (or greater) the longwall recoverable tonnes, per metre of gateroad development and improve the safety standards in longwall operation. Accurate cavability assessment of the coal seam is the key pre-requisite for successful application of LTCC method. Although some cavability criteria were developed in previous studies, their shortcomings limit their application. Apart from the lack of suitable cavability criterion, the caving principles and mechanisms in longwall operation (and LTCC) is still not well understood. The main objectives of this research work were: 1. to improve the general understanding about the caving mechanisms involved in LTCC operation, 2. to develop a new cavability assessment criterion, and 3. to characterise and classify the caving behaviour in selected Australian longwall operations. For this study a combination of analytical, observational and empirical engineering methods were used to reduce the shortcoming of individual design methods. For analytical study, the Distinct Element Method (DEM) was selected for computational analysis as the most suitable technique for this type of study. For observational methods, the size distribution of particles in goaf was used as a monitoring measure for caving performance evaluation and this method was performed in three Australian longwall operations (Ulan, West Wallsend and Broadmeadow). The Empirical assessment was performed by using the results of observational/computational analyses along with basic geological/geotechnical data from selected mine sites to back analyse and re-evaluate the results. The most significant outcomes of this study include: 1. a new cavability assessment system (Top Coal Cavability Rating, TCCR) that was developed based on computational analysis and by back analysis of past LTCC experience in China, 2. Improved understanding of caving mechanism, deformation mechanism and drawing mechanism of top coal in LTCC method, 3. advanced numerical models (UDEC, 3DEC and PFC models) that not only can be used for future research but can also be employed for LTCC design purposes, 4. characterisation and classification of caving behaviour in three Australian longwall operations. The results of this study can be effectively used by Australian coal industry in feasibility and design studies of LTCC operation. The new cavability criterion can be used to identify and quantify the major risks that may be involved in LTCC application.

The tropical cyclone is a majestic, yet violent atmospheric weather system occurring over tropical waters. Their majesty evolves from the significant range of spatial scales they operate over: from the mesoscale, to the larger synoptic-scale. Their associated violent winds and seas, however, are often the cause of damage and destruction for settlements in their path. Between 10/11/07 and 16/11/07, tropical cyclone Sidr formed and intensified into a category 5 hurricane over the southeast tropical waters of the northern Indian Ocean. Sidr tracked west, then north, during the course of its life, and eventually made landfall on 15/11/07, as a category 4 cyclone near the settlement of Barguna, Bangladesh. The storm affected approximately 2.7 million people in Bangladesh, and of that number 4234 were killed. In this study, the dynamics of tropical cyclone Sidr are simulated using version 2.2.1 of Advanced Weather Research and Forecasting — a non-hydrostatic, two-way interactive, triply-nested-grid mesoscale model. Three experiments were developed examining model sensitivity to ocean-atmosphere interaction; initialisation time; and choice of convective parameterisation scheme. All experiments were verified against analysed synoptic data. The ocean-atmosphere experiment involved one simulation of a cold sea surface temperature, fixed at 10 °C; and simulated using a 15 km grid resolution. The initialisation experiment involved three simulations of different model start time: 108-, 72-, and 48-hours before landfall respectively. These were simulated using a 15 km grid resolution. The convective experiment consisted of four simulations, with three of these using a different implicit convective scheme. The three schemes used were, the Kain-Fritsch, Betts-Miller-Janjic, and Grell-Devenyi ensemble. The fourth case simulated convection explicitly. A nested domain of 5km grid spacing was used in the convective experiment, for high resolution modelling. In all experiments, the Eta-Ferrier microphysics scheme, and the Mellor-Yamada-Janjic planetary boundary layer scheme were used. As verified against available observations, the model showed considerable sensitivity in each of the experiments. The model was found to be well suited for combining ocean-atmosphere interactions: a cool sea surface caused cyclone Sidr to dissipate within 24 hours. The initialisation simulations indicated moderate model sensitivity to initialisation time: variations were found for both cyclone track and intensity. Of the three simulations, an initialisation time 108 hours prior to landfall, was found to most accurately represent cyclone Sidr’s track and intensity. Finally, the convective simulations showed that considerable differences were found in cyclone track, intensity, and structure, when using different convective schemes. The Kain-Fritsch scheme produced the most accurate cyclone track and structure, but the rainfall rate was spurious on the sub-grid-scale. The Betts-Miller-Janjic scheme resolved realistic rainfall on both domains, but cyclone intensity was poor. Of particular significance, was that explicit convection produced a similar result to the Grell-Devenyi ensemble for both model domain resolutions. Overall, the results suggest that the modelled cyclone is highly sensitive to changes in initial conditions. In particular, in the context of other studies, it appears that the combination of convective scheme, microphysics scheme, and boundary layer scheme, are most significant for accurate track and intensity prediction.

In the Asian and other monsoon regions of the world most of the severe weather observed is local or mesoscale in nature. Forecasting convective storms or mesoscale systems in the monsoon regions, especially in the tropics, has always been a challenging task to operational meteorologists. Maldives Islands, being situated in the tropical Indian Ocean, are affected by monsoon depressions and tropical cyclones. Thunderstorms and the passage of squall lines are well known sources of heavy rainfall. However, due to the lack of professional people and necessary equipment the weather systems around these islands are seldom studied. Therefore the aim of this thesis is to investigate whether the small islands can create sufficient perturbations in the mesoscale environment to result in the development of convective systems. In this regard, two numerical models, Weather Research and Forecasting model (WRF version 2.2.1) and Regional Atmospheric Modelling System (RAMS version 6.0) were used in this study. Two experiments were performed using the WRF model. In the first experiment, a case study was investigated where the selected day experienced heavy rainfall and thunderstorms. In the second experiment, the same case study was used but with the topographical and surface properties removed in order to investigate the influence of the island in modifying the mesoscale environment. All the experiments were initialized using the re-analysis data from NECP. WRF was able to predict the large scale synoptic features with reasonable accuracy when compared to the observations. Development of the boundary layer and the downstream advection of the temperature anomaly generated by the island were well represented. However, the magnitude of the effects was shown to be weak, probably due to the influence of large scale synoptic features. Even though the model was able to predict the large scale features and some of the mesoscale features, it did not predict any storm development and underestimated the precipitation. Therefore, it was decided to idealize the storm development using the RAMS model. RAMS model was used in a two-dimensional framework. The model was initialized horizontally homogenous using a single sounding and six simulations were performed. The simulation results clearly depicted that the small island can generate its own circulation and influence the mesoscale environment. The daytime heating of the island and the downstream advection of the temperature anomaly in a moist unstable atmosphere could trigger a thunderstorm later in the day. The storm becomes mature approximately 40-80 km offshore. This also suggests that triggering of a storm on one side of an atoll could influence the islands on the downstream side. Sensitivity of storm development to the thermodynamics showed that even with an unstable atmosphere, enough moisture in the lower and mid-troposphere is needed to trigger the storm. Sensitivity to the change of SST showed that convective development was suppressed with a drop of 1 oC. However, this needs further investigation. Assessment of sensitivity to the size of the island showed that the time of triggering of the storm was later and the scale of influence was smaller with a smaller island.

In the Kane silicon-based electron-mediated nuclear spin quantum computer architecture, phosphorus is doped at precise positions in a silicon lattice, and the P donor nuclear spins act as qubits. Logical operations on the nuclear spins are performed using externally applied magnetic and electric fields. There are two important interactions: the hyperfine and exchange interactions, crucial for logical qubit operations. Single qubit operations are performed by applying radio frequency magnetic fields resonant with targeted nuclear spin transition frequencies, tuned by the gate-controlled hyperfine interaction. Two qubit operations are mediated through the exchange interaction between adjacent donor electrons. It is important to examine how these two interactions vary as functions of experimental parameters. Here we provide such an investigation. First, we examine the effects of varying several experimental parameters: gate voltage, magnetic field strength, inter donor separation, donor depth below the silicon oxide interface and back gate depth, to explore how these variables affect the donor electron density. Second, we calculate the hyperfine interaction and the exchange coupling as a function of these parameters. These calculations were performed using various levels of effective mass theory. In the first part of this thesis we use a multi-valley effective mass approach where we incorporate the full Si crystal Bloch structure in calculating the donor electron energy in the bulk silicon. Including the detailed Bloch structure is very computationally intensive, thus when we considered the effect of the externally applied fields in the second and third part, we employ an approach where we focus on the smooth donor-modulated envelope function to determine the response of the donor electron to the applied electric and magnetic fields and qubit position in the lattice. The electric field potential was obtained using Technology Computer Aided Design software, and the interfaces were modelled as a barrier using a step function. One of the critical results of this theoretical study was finding that there exist two regimes for the behaviour of the donor electron in response to the applied gate voltage, dependent on donor distance from the gate. When the qubit is in close proximity to the gate the electron transfer to the gate is gradual. However if the qubit is located far enough from the gate, we found that the donor electron is ionised toward the gate for gate voltages above a certain threshold. Another significant development we have made is in our calculations of the exchange coupling between two adjacent donor electrons. We extended our original Heitler-London basis to describe the two-electron system, and adopted a molecular orbital method where we included a a basis of 78 singlet and 66 triplet two-electron states. In addition to calculating a more accurate exchange coupling, we also evaluated the energy spectrum of the two electron double donor system. We aim to provide relevant information for the experimental design of these devices and highlight the significance of environmental factors other than gate potential that affect the donor electron.

[Truncated abstract] The axial performance of piles in sand remains an area of great uncertainty in geotechnical engineering. Over the years, database studies have shown that the existing method for offshore piles (e.g. API 2000) is unreliable. There is therefore a clear need for an improved predictive method, which incorporates the state-ofthe- art understanding of the underlying controlling mechanisms. This Thesis is dedicated to address the factors influencing the end bearing performance of displacement piles in siliceous sand with a view to proposing and justifying an improved design formulation. Firstly, a database of displacement pile load tests in sand with CPT data was compiled in collaboration with James Schneider (Schneider 2007). It features the widest database with also the latest available pile load test data (e.g. Euripides, Ras Tanajib, Drammen etc) in electronic form. Evaluation of the three new CPTbased methods (Fugro-05, ICP-05 & NGI-05) against this database has revealed a broadly similar predictive performance despite their end bearing formulations being remarkably different. This anomaly promoted the author to extend the database to include additional tests with base capacity measurements to form new base capacity databases for driven and jacked piles, which resulted in the UWA- 05 method for end bearing of displacement piles in sand. This method accounts for the pile effective area ratio, differentiates between driven and jacked piles, and employs a rational qc averaging technique. ... Field tests were performed in Shenton Park, Perth to supplement the database study and, in particular, to examine the effect of the incremental filling ratio (IFR). 10 open-ended and 2 closed-ended piles were tested in compression followed by tension. The test results provide strong support for the UWA-05 method for base capacity evaluation employing the CPT qc values and the effective area ratio. A series of jacked pile tests was carried out on the UWA beam centrifuge, to further explore the factors affecting pile base response. In total, four uniform and four layered centrifuge samples were prepared and tested at various stress levels and relative densities using three separate pile diameters. The resistance ratio (qb0.1/qc,avg) is found to be independent of the absolute pile diameter, effective stress and soil relative density. The tests in layered soil enabled quantification of the reduction in penetration resistance when a pile/cone approaches a weak layer and revealed the significant influence on base stiffness of underlying soft clay layers. The stiffness decay curves (G/GIN vs. w/D, where GIN is initial operational shear stiffness) measured in static load tests were found to vary with ratios of GIN/qc, while there was a unique relationship between G/GIN and qb/qc. A detailed parametric study was carried out (using the FE code PLAXIS) by idealising pile penetration using a spherical cavity expansion analogue in layered soil. The numerical predictions compare well with the centrifuge results and their generalization enabled guidelines to be established for end bearing in layered soil.

The following report is the results of a master thesis project performed on the behalf of Boliden Tara Mines. The work includes a numerical modelling study with the main purpose to evaluate the room- and pillar stability for the mining method "drift-and-slash" mining currently used in the Tara Mine when mining in thin, low-rise lenses. Also, the potential of applying the "drift-and-fill" mining method was studied during the project. The numerical modelling work was performed both two-dimensionally and three-dimensionally, using the Itasca software FLAC and FLAC3D.  The study resulted in recommendations regarding the room- and pillar dimensions, mining sequence, rock support, mining method and further studies. It was concluded that the currently used room- and pillar dimensions could be altered to increase the extraction ratio and that the mining sequence and rock support should be determined depending on the room- and pillar dimensions. Also, that it could be possible to change from the drift-and-slash method to the drift-and-fill method from a stability perspective. Lastly, it was recommended to further study the rock strength and to monitor the pillar behaviour.

Alpine valleys are rarely closed systems, implying that the atmospheric boundary layer of a particular valley section is influenced by the surrounding terrain and large-scale flows. A detailed characterisation and quantification of these effects is required in order to design appropriate parameterisation schemes for complex terrains. The focus of this work is to improve the understanding of the effects of surrounding terrain (plains, valleys or tributaries) on the heat and mass budgets of the stable boundary layer of a valley section, under dry and weak large-scale wind conditions. Numerical simulations using idealised and real frameworks are performed to meet this goal. Several idealised terrains (configurations) were considered: an infinitely long valley (i.e. two-dimensional), and upstream valleys opening either on a plain (valley-plain), on a wider valley (draining) or on a narrower valley (pooling). In three-dimensional valleys, two main regimes can be identified for all configurations: a transient regime, before the down-valley flow develops, followed by a quasi-steady regime, when the down-valley flow is fully developed. The presence of a downstream valley reduces the along-valley temperature difference, therefore leading to weaker down-valley flows. As a result, the duration of the transient regime increases compared to the respective valley-plain configuration. Its duration is longest for the pooling configuration. For strong pooling the along-valley temperature difference can reverse, forcing up-valley flows from the narrower towards the wider valley. In this regime, the average cooling rate at the valley-scale is found to be a maximum and its magnitude is dependent on the configuration considered. Therefore pooling and draining induce colder and deeper boundary layers than the respective valley-plain configurations. In the quasisteady regime the cooling rate is smaller than during the transient regime, and almost independent of the configuration considered. Indeed, as the pooling character is more pronounced, the warming contribution from advection to the heat budget decreases because of weaker down-valley flows, and so does the cooling contribution from the surface sensible heat flux. The mass budget of the valley boundary layer was found to be controlled by a balance between the convergence of downslope flows at the top of the boundary layer and the divergence of the down-valley flow along the valley axis, with negligible contributions of subsidence far from the valley sidewalls. The mass budget highlighted the importance of the return current above the down-valley flow, which may contribute significantly to the inflow of air at the top of the boundary layer. A case-study of a persistent cold-air pool event which occurred in February 2015 in the Arve River Valley during the intensive observation period 1 (IOP1) of the PASSY- 2015 field campaign, allowed us to quantify the effects of neighbouring valleys on the heat and mass budgets of a real valley atmosphere. The cold-air pool persisted as a result of warm air advection at the valley top, associated with the passage of an upper-level ridge over Europe. The contributions from each tributary valley to the mass and heat budgets of the valley atmosphere were found to vary from day to day within the persistent stage of the cold-air pool, depending on the large-scale flow. Tributary flows had significant impact on the height of the inversion layer and the strength of the cold-air pool, transporting a significant amount of mass within the valley atmosphere throughout the night. The strong stratification of the near-surface atmosphere prevented the tributary flows from penetrating down to the valley floor. The evolution of the large-scale flow during the episode had a profound impact on the near-surface circulation of the valley. The channelling of the large-scale flow at night, can lead to the decrease of the horizontal temperature difference driving the near-surface down-valley flow, favouring the stagnation of the air close to the ground.