Centrifuge modelling of submarine landslide flows

Landslides occur both onshore and offshore, however little attention has been given to offshore landslides (submarine landslides). The unique characteristics of submarine landslides include large mass movements and long travel distances at very gentle slopes. Submarine landslides have significant impacts and consequences on offshore and coastal facilities. This paper presents data from a series of centrifuge tests simulating submarine landslide flows on a very gentle slope. Experiments were conducted at different gravity levels to understand the scaling laws involved in simulating submarine landslide flows through centrifuge modelling. The slope was instrumented with miniature sensors for measurements of pore pressure beneath the flow. A series of digital cameras were used to capture the flow in flight. The results provide a better understanding of the scaling laws that needs to be adopted for centrifuge experiments involving submarine landslide flows and gives an insight into the flow mechanisms. © 2010 Taylor & Francis Group, London.

Progress towards the design and numerical analysis of a 3D microchannel biochip separator

This paper reports the design and numerical analysis of a three-dimensional biochip plasma blood separator using computational fluid dynamics techniques. Based on the initial configuration of a two-dimensional (2D) separator, five three-dimensional (3D) microchannel biochip designs are categorically developed through axial and plenary symmetrical expansions. These include the geometric variations of three types of the branch side channels (circular, rectangular, disc) and two types of the main channel (solid and concentric). Ignoring the initial transient behaviour and assuming that steady-state flow has been established, the behaviour of the blood fluid in the devices is algebraically analysed and numerically modelled. The roles of the relevant microchannel mechanisms, i.e. bifurcation, constriction and bending channel, on promoting the separation process are analysed based on modelling results. The differences among the different 3D implementations are compared and discussed. The advantages of 3D over 2D separator in increasing separation volume and effectively depleting cell-free layer fluid from the whole cross section circumference are addressed and illustrated. © 2011 John Wiley & Sons, Ltd.

Biofluid behaviour in 3D microchannel systems: Numerical analysis and design development of 3D microchannel biochip separators

This paper describes the design and development cycle of a 3D biochip separator and the modelling analysis of flow behaviour in the biochip microchannel features. The focus is on identifying the difference between 2D and 3D implementations as well as developing basic forms of 3D microfluidic separators. Five variants, based around the device are proposed and analysed. These include three variations of the branch channels (circular, rectangular, disc) and two variations of the main channel (solid and concentric). Ignoring the initial transient behaviour and assuming steady state flow has been established, the efficiencies of the flow between the main and side channels for the different designs are analysed and compared with regard to relevant biomicrofluidic laws or effects (bifurcation law, Fahraeus effect, cell-free phenomenon, bending channel effect and laminar flow behaviour). The modelling results identify flow features in microchannels, a constriction and bifurcations and show detailed differences in flow fields between the various designs. The manufacturing process using injection moulding for the initial base case design is also presented and discussed. The work reported here is supported as part of the UK funded 3D-MINTEGRATION project. © 2010 IEEE.

Finite element modelling of the sismic behaviour of water front structures

Water front structures have suffered significant damage in many of the recent earthquakes. These include gravity type quay walls, vertically composite walls, cantilever retaining walls, anchored bulkheads and similar structures. One of the primary causes for the poor performance of these classes of structures is the liquefaction of the foundation soil and in some instances liquefaction of the backfill soil. The liquefaction of the soil in-front of the quay wall tends to cause large lateral displacements and rotation of the wall. Often such gravity walls are placed on rubble mound deposited onto the sea bed.This paper presents finite element analyses of such a problem in which strength degradation of the foundation soil and the backfill material will be modelled using PZ mark III constitutive model. The performance of the wall in terms of its lateral displacement, vertical settlement and/or the rotation suffered by the wall will be presented. In addition, the contours of the horizontal and vertical effective stresses and the excess pore pressure ratio will be presented at different time instants together with hyrdraulic gradients. Immediately after the earthquake, the hydraulic gradients indicate migration of pore water into the region below the wall, suggesting further softening of the foundation soil below the wall.

Modelling and Fragility Analysis of Non-Ductile Reinforced Concrete Buildings in Low-to-Moderate Seismic Zones

Reinforced concrete buildings in low-to-moderate seismic zones are often designed only for gravity loads in accordance with the non-seismic detailing provisions. Deficient detailing of columns and beam-column joints can lead to unpredictable brittle failures even under moderate earthquakes. Therefore, a reliable estimate of structural response is required for the seismic evaluation of these structures. For this purpose, analytical models for both interior and exterior slab-beam-column subassemblages and for a 1/3 scale model frame were implemented into the nonlinear finite element platform OpenSees. Comparison between the analytical results and experimental data available in the literature is carried out using nonlinear pushover analyses and nonlinear time history analysis for the subassemblages and the model frame, respectively. Furthermore, the seismic fragility assessment of reinforced concrete buildings is performed on a set of non-ductile frames using nonlinear time history analyses. The fragility curves, which are developed for various damage states for the maximum interstory drift ratio are characterized in terms of peak ground acceleration and spectral acceleration using a suite of ground motions representative of the seismic hazard in the region.

Initial phase modelling in numerical explosion applied to process safety

The utilisation of computational fluid dynamics (CFD) in process safety has increased significantly in recent years. The modelling of accidental explosion via CFD has in many cases replaced the classical Multi Energy and Brake Strehlow methods. The benefits obtained with CFD modelling can be diminished if proper modelling of the initial phase of explosion is neglected. In the early stages of an explosion, the flame propagates in a quasi-laminar regime. Proper modelling of the initial laminar phase is a key aspect in order to predict the peak pressure and the time to peak pressure. The present work suggests a modelling approach for the initial laminar phase in explosion scenarios. Findings are compared with experimental data for two classical explosion test cases which resemble the common features in chemical process areas (confinement and congestion). A detailed analysis of the threshold for the transition from laminar to turbulent regime is also carried out. The modelling is implemented in a fully 3D Navier-Stokes compressible formulation. Combustion is treated using a laminar flamelet approach based on the Bray, Moss and Libby (BML) formulation. A novel modified porosity approach developed for the unstructured solver is also considered. Results agree satisfactorily with experiments and the modelling is found to be robust. © 2013 The Institution of Chemical Engineers.

Behaviour of cast-iron tunnel segmental joint from the 3D FE analyses and development of a new bolt-spring model

The behaviour of cast-iron tunnel segments used in London Underground tunnels was investigated using the 3-D finite element (FE) method. A numerical model of the structural details of cast-iron segmental joints such as bolts, panel and flanges was developed and its performance was validated against a set of full-scale tests. Using the verified model, the influence of structural features such as caulking groove and bolt pretension was examined for both rotational and shear loading conditions. Since such detailed modelling of bolts increases the computational time when a full scale segmental tunnel is analysed, it is proposed to replace the bolt model to a set of spring models. The parameters for the bolt-spring models, which consider the geometry and material properties of the bolt, are proposed. The performance of the combined bolt-spring and solid segmental models are evaluated against a more conventional shell-spring model. © 2014 Elsevier Ltd.

A practical system for modelling body shapes from single view measurements

This paper describes an interactive system for quickly modelling 3D body shapes from a single image. It provides the user with a convenient way to obtain their 3D body shapes so as to try on virtual garments online. For the ease of use, we first introduce a novel interface for users to conveniently extract anthropometric measurements from a single photo, while using readily available scene cues for automatic image rectification. Then, we propose a unified probabilistic framework using Gaussian processes, which predict the body parameters from input measurements while correcting the aspect ratio ambiguity resulting from photo rectification. Extensive experiments and user studies have supported the efficacy of our system. This system is now being exploited commercially online1. © 2011. The copyright of this document resides with its authors.

Modelling of settlement induced building damage

This thesis focuses on the modelling of settlement induced damage to masonry buildings. In densely populated areas, the need for new space is nowadays producing a rapid increment of underground excavations. Due to the construction of new metro lines, tunnelling activity in urban areas is growing. One of the consequences is a greater attention to the risk of damage on existing structures. Thus, the assessment of potential damage of surface buildings has become an essential stage in the excavation projects in urban areas (Chapter 1). The current damage risk assessment procedure is based on strong simplifications, which not always lead to conservative results. Object of this thesis is the development of an improved damage classification system, which takes into account the parameters influencing the structural response to settlement, like the non-linear behaviour of masonry and the soil-structure interaction. The methodology used in this research is based on experimental and numerical modelling. The design and execution of an experimental benchmark test representative of the problem allows to identify the principal factors and mechanisms involved. The numerical simulations enable to generalize the results to a broader range of physical scenarios. The methodological choice is based on a critical review of the currently available procedures for the assessment of settlement-induced building damage (Chapter 2). A new experimental test on a 1/10th masonry façade with a rubber base interface is specifically designed to investigate the effect of soil-structure interaction on the tunnelling-induced damage (Chapter 3). The experimental results are used to validate a 2D semi-coupled finite element model for the simulation of the structural response (Chapter 4). The numerical approach, which includes a continuum cracking model for the masonry and a non-linear interface to simulate the soil-structure interaction, is then used to perform a sensitivity study on the effect of openings, material properties, initial damage, initial conditions, normal and shear behaviour of the base interface and applied settlement profile (Chapter 5). The results assess quantitatively the major role played by the normal stiffness of the soil-structure interaction and by the material parameters defining the quasi-brittle masonry behaviour. The limitation of the 2D modelling approach in simulating the progressive 3D displacement field induced by the excavation and the consequent torsional response of the building are overcome by the development of a 3D coupled model of building, foundation, soil and tunnel (Chapter 6). Following the same method applied to the 2D semi-coupled approach, the 3D model is validated through comparison with the monitoring data of a literature case study. The model is then used to carry out a series of parametric analyses on geometrical factors: the aspect ratio of horizontal building dimensions with respect to the tunnel axis direction, the presence of adjacent structures and the position and alignment of the building with respect to the excavation (Chapter 7). The results show the governing effect of the 3D building response, proving the relevance of 3D modelling. Finally, the results from the 2D and 3D parametric analyses are used to set the framework of an overall damage model which correlates the analysed structural features with the risk for the building of being damaged by a certain settlement (Chapter 8). This research therefore provides an increased experimental and numerical understanding of the building response to excavation-induced settlements, and sets the basis for an operational tool for the risk assessment of structural damage (Chapter 9).

Modelling and comparison of trapped fields in (RE)BCO bulk superconductors for activation using pulsed field magnetization

The ability to generate a permanent, stable magnetic field unsupported by an electromotive force is fundamental to a variety of engineering applications. Bulk high temperature superconducting (HTS) materials can trap magnetic fields of magnitude over ten times higher than the maximum field produced by conventional magnets, which is limited practically to rather less than 2 T. In this paper, two large c-axis oriented, single-grain YBCO and GdBCO bulk superconductors are magnetized by the pulsed field magnetization (PFM) technique at temperatures of 40 and 65 K and the characteristics of the resulting trapped field profile are investigated with a view of magnetizing such samples as trapped field magnets (TFMs) in situ inside a trapped flux-type superconducting electric machine. A comparison is made between the temperatures at which the pulsed magnetic field is applied and the results have strong implications for the optimum operating temperature for TFMs in trapped flux-type superconducting electric machines. The effects of inhomogeneities, which occur during the growth process of single-grain bulk superconductors, on the trapped field and maximum temperature rise in the sample are modelled numerically using a 3D finite-element model based on the H-formulation and implemented in Comsol Multiphysics 4.3a. The results agree qualitatively with the observed experimental results, in that inhomogeneities act to distort the trapped field profile and reduce the magnitude of the trapped field due to localized heating within the sample and preferential movement and pinning of flux lines around the growth section regions (GSRs) and growth sector boundaries (GSBs), respectively. The modelling framework will allow further investigation of various inhomogeneities that arise during the processing of (RE)BCO bulk superconductors, including inhomogeneous Jc distributions and the presence of current-limiting grain boundaries and cracks, and it can be used to assist optimization of processing and PFM techniques for practical bulk superconductor applications. © 2014 IOP Publishing Ltd.

Automatic 3d free form shape matching using the graduated assignment algorithm.

Liu, Yonghuai. Automatic 3d free form shape matching using the graduated assignment algorithm. Pattern Recognition, vol. 38, no. 10, pp. 1615-1631, 2005.

Multi-scale modelling of atomic layer deposition

High-permittivity ("high-k") dielectric materials are used in the transistor gate stack in integrated circuits. As the thickness of silicon oxide dielectric reduces below 2 nm with continued downscaling, the leakage current because of tunnelling increases, leading to high power consumption and reduced device reliability. Hence, research concentrates on finding materials with high dielectric constant that can be easily integrated into a manufacturing process and show the desired properties as a thin film. Atomic layer deposition (ALD) is used practically to deposit high-k materials like HfO2, ZrO2, and Al2O3 as gate oxides. ALD is a technique for producing conformal layers of material with nanometer-scale thickness, used commercially in non-planar electronics and increasingly in other areas of science and technology. ALD is a type of chemical vapor deposition that depends on self-limiting surface chemistry. In ALD, gaseous precursors are allowed individually into the reactor chamber in alternating pulses. Between each pulse, inert gas is admitted to prevent gas phase reactions. This thesis provides a profound understanding of the ALD of oxides such as HfO2, showing how the chemistry affects the properties of the deposited film. Using multi-scale modelling of ALD, the kinetics of reactions at the growing surface is connected to experimental data. In this thesis, we use density functional theory (DFT) method to simulate more realistic models for the growth of HfO2 from Hf(N(CH3)2)4/H2O and HfCl4/H2O and for Al2O3 from Al(CH3)3/H2O.Three major breakthroughs are discovered. First, a new reaction pathway, ’multiple proton diffusion’, is proposed for the growth of HfO2 from Hf(N(CH3)2)4/H2O.1 As a second major breakthrough, a ’cooperative’ action between adsorbed precursors is shown to play an important role in ALD. By this we mean that previously-inert fragments can become reactive once sufficient molecules adsorb in their neighbourhood during either precursor pulse. As a third breakthrough, the ALD of HfO2 from Hf(N(CH3)2)4 and H2O is implemented for the first time into 3D on-lattice kinetic Monte-Carlo (KMC).2 In this integrated approach (DFT+KMC), retaining the accuracy of the atomistic model in the higher-scale model leads to remarkable breakthroughs in our understanding. The resulting atomistic model allows direct comparison with experimental techniques such as X-ray photoelectron spectroscopy and quartz crystal microbalance.

Modelling the fluid dynamics and coupled phenomena in arc weld pools

A 3D model of melt pool created by a moving arc type heat sources has been developed. The model solves the equations of turbulent fluid flow, heat transfer and electromagnetic field to demonstrate the flow behaviour phase-change in the pool. The coupled effects of buoyancy, capillary (Marangoni) and electromagnetic (Lorentz) forces are included within an unstructured finite volume mesh environment. The movement of the welding arc along the workpiece is accomplished via a moving co-ordinator system. Additionally a method enabling movement of the weld pool surface by fluid convection is presented whereby the mesh in the liquid region is allowed to move through a free surface. The surface grid lines move to restore equilibrium at the end of each computational time step and interior grid points then adjust following the solution of a Laplace equation.

Modelling electromagnetically levitated liquid droplet oscillations

This work comprises accurate computational analysis of levitated liquid droplet oscillations in AC and DC magnetic fields. The AC magnetic field interacting with the induced electric current within the liquid metal droplet generates intense fluid flow and the coupled free surface oscillations. The pseudo-spectral technique is used to solve the turbulent fluid flow equations for the continuously dynamically transformed axisymmetric fluid volume. The volume electromagnetic force distribution is updated with the shape and position change. We start with the ideal fluid test case for undamped Rayleigh frequency oscillations in the absence of gravity, and then add the viscous and the DC magnetic field damping. The oscillation frequency spectra are further analysed for droplets levitated against gravity in AC and DC magnetic fields at various combinations. In the extreme case electrically poorly conducting, diamagnetic droplet (water) levitation dynamics are simulated. Applications are aimed at pure electromagnetic material processing techniques and the material properties measurements in uncontaminated conditions.

AC & DC magnetic levitation and semi-levitation modelling

This work presents computation analysis of levitated liquid thermal and flow fields with free surface oscillations in AC and DC magnetic fields. The volume electromagnetic force distribution is continuously updated with the shape and position change. The oscillation frequency spectra are analysed for droplets levitation against gravity in AC and DC magnetic fields at various combinations. For larger volume liquid metal confinement and melting the semi-levitation induction skull melting process is simulated with the same numerical model. Applications are aimed at pure electromagnetic material processing techniques and the material properties measurements in uncontaminated conditions.