Explosion pressure prediction via polynomial mathematical correlation based on advanced CFD Modelling

Computational Fluid Dynamics CFD can be used as a powerful tool supporting engineers throughout the steps of the design. The combination of CFD with response surface methodology can play an important role in such cases. During the conceptual engineering design phase, a quick response is always a matter of urgency. During this phase even a sketch of the geometrical model is rare. Therefore, the utilisation of typical response surface developed for congested and confined environment rather than CFD can be an important tool to help the decision making process, when the geometrical model is not available, provided that similarities can be considered when taking into account the characteristic of the geometry in which the response surface was developed. The present work investigates how three different types of response surfaces behave when predicting overpressure in accidental scenarios based on CFD input. First order, partial second order and complete second order polynomial expressions are investigated. The predicted results are compared with CFD findings for a classical offshore experiment conducted by British Gas on behalf of Mobil and good agreement is observed for higher order response surfaces. The higher order response surface calculations are also compared with CFD calculations for a typical offshore module and good agreement is also observed. © 2011 Elsevier Ltd.

Discrete dislocation plasticity analysis of size effects in single crystals

The effect of size and slip system configuration on the tensile stress-strain response of micron-sized planar crystals as obtained from discrete dislocation plasticity simulations is presented. The crystals are oriented for either single or symmetric double slip. With the rotation of the tensile axis unconstrained, there is a strong size dependence, with the flow strength increasing with decreasing specimen size. Below a certain specimen size, the flow strength of the crystals is set by the nucleation strength of the initially present Frank-Read sources. The main features of the size dependence are the same for both the single and symmetric double slip configurations.

Discrete element calculations of the impact of a sand column against rigid structures

Discrete particle simulations of column of an aggregate of identical particles impacting a rigid, fixed target and a rigid, movable target are presented with the aim to understand the interaction of an aggregate of particles upon a structure. In most cases the column of particles is constrained against lateral expansion. The pressure exerted by the particles upon the fixed target (and the momentum transferred) is independent of the co-efficient of restitution and friction co-efficient between the particles but are strongly dependent upon the relative density of the particles in the column. There is a mild dependence on the contact stiffness between the particles which controls the elastic deformation of the densified aggregate of particles. In contrast, the momentum transfer to a movable target is strongly sensitive to the mass ratio of column to target. The impact event can be viewed as an inelastic collision between the sand column and the target with an effective co-efficient of restitution between 0 and 0.35 depending upon the relative density of the column. We present a foam analogy where impact of the aggregate of particles can be modelled by the impact of an equivalent foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally and yet give predictions to within 5% of the full discrete particle calculations. They also suggest that "model" materials can be used to simulate the loading by an aggregate of particles within a laboratory setting. © 2012 Elsevier Ltd. All rights reserved.

Stress-strain response of hydrate-bearing sands: Numerical study using discrete element method simulations

Gas hydrate is a crystalline solid found within marine and subpermafrost sediments. While the presence of hydrates can have a profound effect on sediment properties, the stress-strain behavior of hydrate-bearing sediments is poorly understood due to inherent limitations in laboratory testing. In this study, we use numerical simulations to improve our understanding of the mechanical behavior of hydrate-bearing sands. The hydrate mass is simulated as either small randomly distributed bonded grains or as "ripened hydrate" forming patchy saturation, whereby sediment clusters with 100% pore-filled hydrate saturation are distributed within a hydrate-free sediment. Simulation results reveal that reduced sand porosity and higher hydrate saturation cause an increase in stiffness, strength, and dilative tendency, and the critical state line shifts toward higher void ratio and higher shear strength. In particular, the critical state friction angle increases in sands with patchy saturation, while the apparent cohesion is affected the most when the hydrate mass is distributed in pores. Sediments with patchy hydrate distribution exhibit a slightly lower strength than sediments with randomly distributed hydrate. Finally, hydrate dissociation under drained conditions leads to volume contraction and/or stress relaxation, and pronounced shear strains can develop if the hydrate-bearing sand is subjected to deviatoric loading during dissociation.

Centrifuge and discrete element modelling of tunnelling effects on pipelines

The paper presents centrifuge test data of the problem of tunnelling effects on buried pipelines and compares them to predictions made using DEM simulations. The paper focuses on the examination of pipeline bending moments, their distribution along the pipe, and their development with tunnel volume loss. Centrifuge results are obtained by PIV analysis and compared to results obtained using the DEM model. The DEM model was built to replicate the centrifuge model as closely as possible and included numerical features formulated specially for this task, such as structural elements to replicate the tunnel and pipeline. Results are extremely encouraging, with deviations between DEM and centrifuge test bending moment results being very small. © 2010 Taylor & Francis Group, London.

Representations and algorithms for finite-state bisimulations of linear discrete-time control systems

While a large amount of research over the past two decades has focused on discrete abstractions of infinite-state dynamical systems, many structural and algorithmic details of these abstractions remain unknown. To clarify the computational resources needed to perform discrete abstractions, this paper examines the algorithmic properties of an existing method for deriving finite-state systems that are bisimilar to linear discrete-time control systems. We explicitly find the structure of the finite-state system, show that it can be enormous compared to the original linear system, and give conditions to guarantee that the finite-state system is reasonably sized and efficiently computable. Though constructing the finite-state system is generally impractical, we see that special cases could be amenable to satisfiability based verification techniques. ©2009 IEEE.

Modeling fatigue crack growth in crystalline solids with discrete dislocation plasticity

Analyses of crack growth under cyclic loading conditions are discussed where plastic flow arises from the motion of large numbers of discrete dislocations and the fracture properties are embedded in a cohesive surface constitutive relation. The formulation is the same as used to analyse crack growth under monotonic loading conditions, differing only in the remote loading being a cyclic function of time. Fatigue, i.e. crack growth in cyclic loading at a driving force for which the crack would have arrested under monotonic loading, emerges in the simulations as a consequence of the evolution of internal stresses associated with the irreversibility of the dislocation motion. A fatigue threshold, Paris law behaviour, striations, the accelerated growth of short cracks and the scaling with material properties are outcomes of the calculations. Results for single crystals and polycrystals will be discussed.

Robust network reconstruction in polynomial time

This paper presents an efficient algorithm for robust network reconstruction of Linear Time-Invariant (LTI) systems in the presence of noise, estimation errors and unmodelled nonlinearities. The method here builds on previous work [1] on robust reconstruction to provide a practical implementation with polynomial computational complexity. Following the same experimental protocol, the algorithm obtains a set of structurally-related candidate solutions spanning every level of sparsity. We prove the existence of a magnitude bound on the noise, which if satisfied, guarantees that one of these structures is the correct solution. A problem-specific model-selection procedure then selects a single solution from this set and provides a measure of confidence in that solution. Extensive simulations quantify the expected performance for different levels of noise and show that significantly more noise can be tolerated in comparison to the original method. © 2012 IEEE.