Development of realistic liquid equations of state for multi-phase CFD modelling

Several equations of state (EOS) have been incorporated into a novel algorithm to solve a system of multi-phase equations in which all phases are assumed to be compressible to varying degrees. The EOSs are used to both supply functional relationships to couple the conservative variables to the primitive variables and to calculate accurately thermodynamic quantities of interest, such as the speed of sound. Each EOS has a defined balance of accuracy, robustness and computational speed; selection of an appropriate EOS is generally problem-dependent. This work employs an AUSM+-up method for accurate discretisation of the convective flux terms with modified low-Mach number dissipation for added robustness of the solver. In this paper we show a newly-developed time-marching formulation for temporal discretisation of the governing equations with incorporated time-dependent source terms, as well as considering the system of eigenvalues that render the governing equations hyperbolic.

A Force-Balanced Control Volume Finite Element Method for Multi-Phase Porous Media Flow Modelling

Peer reviewed

Multi-scale modelling of Suzuki segregation in gamma' precipitates in Ni and Co-base superalloys

The high temperature strength of alloys with (gamma +gamma') microstructure is primarily due to the resistance of the ordered precipitate to cutting by matrix dislocations. Such shearing requires higher stresses since it involves the creation of a planar fault. Planar fault energy is known to be dependent on composition. This implies that the composition on the fault may be different from that in the bulk for energetic reasons. Such segregation (or desegregation) of specific alloying elements to the fault may result in Suzuki strengthening which has not been explored extensively in these systems. In this work, segregation (or desegregation) of alloying elements to planar faults was studied computationally in Ni-3(Al, Ti) and Co-3(W, Al) type gamma' precipitates. The composition dependence of APB energy and heat of mixing were evaluated from first principle electronic structure calculations. A phase field model incorporating the first principles results, was used to simulate the motion of an extended superdislocation under stress concurrently with composition evolution. Results reveal that in both systems, significant (de) segregation occurs on equilibration. On application of stress, solutes were dragged along with the APB in some cases. Additionally, it was also noted the velocity of the superdislocation under an applied stress is strongly dependent on atomic mobility (i. e. diffusivity).

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.

Multi-scale modelling of granular avalanches

Avalanches, debris flows, and landslides are geophysical hazards, which involve rapid mass movement of granular solids, water and air as a single-phase system. The dynamics of a granular flow involve at least three distinct scales: the micro-scale, meso-scale, and the macro-scale. This study aims to understand the ability of continuum models to capture the micro-mechanics of dry granular collapse. Material Point Method (MPM), a hybrid Lagrangian and Eulerian approach, with Mohr-Coulomb failure criterion is used to describe the continuum behaviour of granular column collapse, while the micromechanics is captured using Discrete Element Method (DEM) with tangential contact force model. The run-out profile predicted by the continuum simulations matches with DEM simulations for columns with small aspect ratios ('h/r' < 2), however MPM predicts larger run-out distances for columns with higher aspect ratios ('h/r' > 2). Energy evolution studies in DEM simulations reveal higher collisional dissipation in the initial free-fall regime for tall columns. The lack of a collisional energy dissipation mechanism in MPM simulations results in larger run-out distances. Micro-structural effects, such as shear band formations, were observed both in DEM and MPM simulations. A sliding flow regime is observed above the distinct passive zone at the core of the column. Velocity profiles obtained from both the scales are compared to understand the reason for a slow flow run-out mobilization in MPM simulations. © 2013 AIP Publishing LLC.

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.

Multi-scale modelling of DP590 steel

Dual Phase (DP) steel one of the Advanced High Strength Steels (AHSS) has a two phase microstructure where soft and hard phase acts together to offer a high strength composite effect. The high strength, however, must be balanced with ductility so that complex parts and designs can be manufactured from AHSS sheets. However, during forming certain grades of DP steel a sudden crack can occur without any intimation of necking. Thus, due to this abnormal forming behaviour, is difficult to accurately predict because most classical modelling approaches are not designed for such micro-structurally heterogeneous materials. These modelling approaches are generally based on an average representation of the material behaviour in a continuum mechanics formulation. This works for materials that are homogenous, or at least could be assumed to be homogenous at scales lower than the naked eye can see. However, for a material like AHSS, the microstructure plays a significant role in dictating the mechanical behaviour at the macro-scale. This paper studies the multi-scale modelling ofDP590 steel. It is found that the sufficient accuracy can be achieved from multi-scale modelling while comparing with the experiments.

Dynamic load balancing of distributed memory parallel computational mechanics using unstructured meshes for multi-physical modelling

As the complexity of parallel applications increase, the performance limitations resulting from computational load imbalance become dominant. Mapping the problem space to the processors in a parallel machine in a manner that balances the workload of each processors will typically reduce the run-time. In many cases the computation time required for a given calculation cannot be predetermined even at run-time and so static partition of the problem returns poor performance. For problems in which the computational load across the discretisation is dynamic and inhomogeneous, for example multi-physics problems involving fluid and solid mechanics with phase changes, the workload for a static subdomain will change over the course of a computation and cannot be estimated beforehand. For such applications the mapping of loads to process is required to change dynamically, at run-time in order to maintain reasonable efficiency. The issue of dynamic load balancing are examined in the context of PHYSICA, a three dimensional unstructured mesh multi-physics continuum mechanics computational modelling code.

Multi-body modelling of wagon trains for crashworthiness analysis

Three-dimensional wagon train models have been developed for the crashworthiness analysis using multi-body dynamics approach. The contributions of the train size (number of wagon) to the frontal crash forces can be identified through the simulations. The effects of crash energy management (CEM) design and crash speed on train crashworthiness are examined. The CEM design can significantly improve the train crashworthiness and the consequential vehicle stability performance - reducing derailment risks.

Destinations, disasters and public relations: Stakeholder engagement in multi-phase disaster management

Research about disasters in tourism has emerged in earnest since the 1990s covering insights for preparedness and response. However, recently, authors have called for more systematic and holistic approaches to tourism disaster management research. To address this gap, this study adopted a public relations perspective to refocus attention to relationships and stakeholder expectations of destination communities across multiple phases of disaster management. The authors used a mixed method approach and developed a battery of disaster management attributes by conducting interviews and analyzing industry documents and the extant literature. These attributes formed part of a survey of tourism businesses. Exploratory factor analysis resulted in a two factor solution: - i) business disaster preparedness, and; - ii) destination disaster response and recovery. Findings also show that participants reported a gap between the importance and destination performance of these attributes. In particular, tourism businesses perceived destinations did not adequately engage in disaster preparedness activities, which had implications for disaster response and recovery.

Multi power port converter for hybrid electric vehicles using multi phase bidirectional fly-back topology

For hybrid electric vehicles the batteries and the drive dc-link may be at different voltages. The batteries are at low voltage to obtain higher volumetric efficiencies and the dc-link is at higher voltage to have higher efficiency on the motor side. Therefore a power interface between the batteries and the drive's dc-link is essential. This power interface should handle power flow from battery to motor, motor to battery, external genset to battery and grid to battery. This paper proposes a multi power port topology which is capable of handling multiple power sources and still maintains simplicity and features like obtaining any gain, wide load variations, lower output current ripple and capability of parallel battery energy due to the modular structure. The development and testing of a bi-directional fly-back DC-DC converter for hybrid electric vehicle is described in this paper. Simple hysteresis voltage control is used for DC link voltage regulation. The experimental results are presented to show the working of the proposed converter.

Differential Multi-phase DLL for Reconfigurable Radio Frequency Synthesizer

A multi phase, delay-locked loop (DLL) based frequency synthesizer is designed for harmonic rejection mixing in reconfigurable radios. This frequency synthesizer uses a 1 GHz input reference frequency, and achieves <= 20ns settling time by utilizing a wide loop bandwidth. The circuit has been designed in 0.13-mu m CMOS technology. It is designed for a frequency range of 500 MHz to 3 GHz with stuck/harmonic lock removal assist. Index Terms-stuck lock, harmonic lock, delay-locked loops, multi phase, phase detector, frequency synthesis

A Multi-phase Approach for Improving Information Diffusion in Social Networks

For maximizing influence spread in a social network, given a certain budget on the number of seed nodes, we investigate the effects of selecting and activating the seed nodes in multiple phases. In particular, we formulate an appropriate objective function for two-phase influence maximization under the independent cascade model, investigate its properties, and propose algorithms for determining the seed nodes in the two phases. We also study the problem of determining an optimal budget-split and delay between the two phases.

Large Scale Numerical Simulations For Multi-Phase Fluid Dynamics With Moving Interfaces

This short communication presents our recent studies to implement numerical simulations for multi-phase flows on top-ranked supercomputer systems with distributed memory architecture. The numerical model is designed so as to make full use of the capacity of the hardware. Satisfactory scalability in terms of both the parallel speed-up rate and the size of the problem has been obtained on two high rank systems with massively parallel processors, the Earth Simulator (Earth simulator research center, Yokohama Kanagawa, Japan) and the TSUBAME (Tokyo Institute of Technology, Tokyo, Japan) supercomputers.