Study on small-strain behaviours of methane hydrate sandy sediments using discrete element method

Methane hydrate bearing soil has attracted increasing interest as a potential energy resource where methane gas can be extracted from dissociating hydrate-bearing sediments. Seismic testing techniques have been applied extensively and in various ways, to detect the presence of hydrates, due to the fact that hydrates increase the stiffness of hydrate-bearing sediments. With the recognition of the limitations of laboratory and field tests, wave propagation modelling using Discrete Element Method (DEM) was conducted in this study in order to provide some particle-scale insights on the hydrate-bearing sandy sediment models with pore-filling and cementation hydrate distributions. The relationship between shear wave velocity and hydrate saturation was established by both DEM simulations and analytical solutions. Obvious differences were observed in the dependence of wave velocity on hydrate saturation for these two cases. From the shear wave velocity measurement and particle-scale analysis, it was found that the small-strain mechanical properties of hydrate-bearing sandy sediments are governed by both the hydrate distribution patterns and hydrate saturation. © 2013 AIP Publishing LLC.

Modelling the entrainment and motion of particles in a gap: Application to abrasive wear

The microscale abrasion or ball-cratering test is being increasingly applied to a wide range of bulk materials and coatings. The response of materials to this test depends critically on the nature of the motion of the abrasive particles in the contact zone: whether they roll and produce multiple indentations in the coating, or slide causing grooving abrasion. Similar phenomena also occur when hard contaminant particles enter a lubricated contact. This paper presents simple quantitative two-dimensional models which describe two aspects of the interaction between a hard abrasive particle and two sliding surfaces. The first model treats the conditions under which a spherical abrasive particle of size d can be entrained into the gap between a rotating sphere of radius R and a plane surface. These conditions are determined by the coefficients of friction between the particle and the sphere, and the particle and the plane, denoted by μs and μp respectively. This model predicts that the values of (μs + μp) and 2μs should both exceed √2d/R for the particles to be entrained into the contact. If either is less than this value, the particle will slide against the sphere and never enter the contact. The second model describes the mechanisms of abrasive wear in a contact when an idealized rhombus-sectioned prismatic particle is located between two parallel plane surfaces separated by a certain distance, which can represent either the thickness of a fluid film or the spacing due to the presence of other particles. It is shown that both the ratio of particle size to the separation of the surfaces and the ratio of the hardnesses of the two surfaces have important influences on the particle motion and hence on the mechanism of the resulting abrasive wear. Results from this model are compared with experimental observations, and the model is shown to lead to realistic predictions. © IMechE 2003.

Numerical simulations of unsteady wet steam flows with non-equilibrium condensation in the nozzle and the steam turbine

Numerical techniques for non-equilibrium condensing flows are presented. Conservation equations for homogeneous gas-liquid two-phase compressible flows are solved by using a finite volume method based on an approximate Riemann solver. The phase change consists of the homogeneous nucleation and growth of existing droplets. Nucleation is computed with the classical Volmer-Frenkel model, corrected for the influence of the droplet temperature being higher than the steam temperature due to latent heat release. For droplet growth, two types of heat transfer model between droplets and the surrounding steam are used: a free molecular flow model and a semi-empirical two-layer model which is deemed to be valid over a wide range of Knudsen number. The computed pressure distribution and Sauter mean droplet diameters in a convergent-divergent (Laval) nozzle are compared with experimental data. Both droplet growth models capture qualitatively the pressure increases due to sudden heat release by the non-equilibrium condensation. However the agreement between computed and experimental pressure distributions is better for the two-layer model. The droplet diameter calculated by this model also agrees well with the experimental value, whereas that predicted by the free molecular model is too small. Condensing flows in a steam turbine cascade are calculated at different Mach numbers and inlet superheat conditions and are compared with experiments. Static pressure traverses downstream from the blade and pressure distributions on the blade surface agree well with experimental results in all cases. Once again, droplet diameters computed with the two-layer model give best agreement with the experiments. Droplet sizes are found to vary across the blade pitch due to the significant variation in expansion rate. Flow patterns including oblique shock waves and condensation-induced pressure increases are also presented and are similar to those shown in the experimental Schlieren photographs. Finally, calculations are presented for periodically unsteady condensing flows in a low expansion rate, convergent-divergent (Laval) nozzle. Depending on the inlet stagnation subcooling, two types of self-excited oscillations appear: a symmetric mode at lower inlet subcooling and an asymmetric mode at higher subcooling. Plots of oscillation frequency versus inlet sub-cooling exhibit a hysteresis loop, in accord with observations made by other researchers for moist air flow. Copyright © 2006 by ASME.

An analysis of the cooperative mechano-sensitive feedback between intracellular signaling, focal adhesion development, and stress fiber contractility

Cells communicate with their external environment via focal adhesions and generate activation signals that in turn trigger the activity of the intracellular contractile machinery. These signals can be triggered by mechanical loading that gives rise to a cooperative feedback loop among signaling, focal adhesion formation, and cytoskeletal contractility, which in turn equilibrates with the applied mechanical loads. We devise a signaling model that couples stress fiber contractility and mechano-sensitive focal adhesion models to complete this above mentioned feedback loop. The signaling model is based on a biochemical pathway where IP3 molecules are generated when focal adhesions grow. These IP3 molecules diffuse through the cytosol leading to the opening of ion channels that disgorge Ca2+ from the endoplasmic reticulum leading to the activation of the actin/myosin contractile machinery. A simple numerical example is presented where a one-dimensional cell adhered to a rigid substrate is pulled at one end, and the evolution of the stress fiber activation signal, stress fiber concentrations, and focal adhesion distributions are investigated. We demonstrate that while it is sufficient to approximate the activation signal as spatially uniform due to the rapid diffusion of the IP3 through the cytosol, the level of the activation signal is sensitive to the rate of application of the mechanical loads. This suggests that ad hoc signaling models may not be able to capture the mechanical response of cells to a wide range of mechanical loading events. © 2011 American Society of Mechanical Engineers.

Characterisation and simulation of asphalt road surfaces

Accurate simulation of rolling-tyre vibrations, and the associated noise, requires knowledge of road-surface topology. Full scans of the surface types in common use are, however, not widely available, and are likely to remain so. Ways of producing simulated surfaces from incomplete starting information are thus needed. In this paper, a simulation methodology based solely on line measurements is developed, and validated against a full two-dimensional height map of a real asphalt surface. First the tribological characteristics-asperity height, curvature and nearest-neighbour distributions-of the real surface are analysed. It is then shown that a standard simulation technique, which matches the (isotropic) spectrum and the probability distribution of the height measurements, is unable to reproduce these characteristics satisfactorily. A modification, whereby the inherent granularity of the surface is enforced at the initialisation stage, is introduced, and found to produce simulations whose tribological characteristics are in excellent agreement with the measurements. This method will thus make high-fidelity tyre-vibration calculations feasible for researchers with access to line-scan data only. In addition, the approach to surface tribological characterisation set out here provides a template for efficient cataloguing of road textures, as long as the resulting information can subsequently be used to produce sample realisations. A third simulation algorithm, which successfully addresses this requirement, is therefore also presented. © 2011 Elsevier B.V.

Validation of a lattice Boltzmann model for gas-solid reactions with experiments

A lattice Boltzmann method is used to model gas-solid reactions where the composition of both the gas and solid phase changes with time, while the boundary between phases remains fixed. The flow of the bulk gas phase is treated using a multiple relaxation time MRT D3Q19 model; the dilute reactant is treated as a passive scalar using a single relaxation time BGK D3Q7 model with distinct inter- and intraparticle diffusivities. A first-order reaction is incorporated by modifying the method of Sullivan et al. [13] to include the conversion of a solid reactant. The detailed computational model is able to capture the multiscale physics encountered in reactor systems. Specifically, the model reproduced steady state analytical solutions for the reaction of a porous catalyst sphere (pore scale) and empirical solutions for mass transfer to the surface of a sphere at Re=10 (particle scale). Excellent quantitative agreement between the model and experiments for the transient reduction of a single, porous sphere of Fe 2O 3 to Fe 3O 4 in CO at 1023K and 10 5Pa is demonstrated. Model solutions for the reduction of a packed bed of Fe 2O 3 (reactor scale) at identical conditions approached those of experiments after 25 s, but required prohibitively long processor times. The presented lattice Boltzmann model resolved successfully mass transport at the pore, particle and reactor scales and highlights the relevance of LB methods for modelling convection, diffusion and reaction physics. © 2012 Elsevier Inc.

A Physical Interpretation of Stagnation Pressure and Enthalpy Changes in Unsteady Flow

This paper provides a physical interpretation of the mechanism of stagnation enthalpy and stagnation pressure changes in turbomachines due to unsteady flow, the agency for all work transfer between a turbomachine and an inviscid fluid. Examples are first given to illustrate the direct link between the time variation of static pressure seen by a given fluid particle and the rate of change of stagnation enthalpy for that particle. These include absolute stagnation temperature rises in turbine rotor tip leakage flow, wake transport through downstream blade rows, and effects of wake phasing on compressor work input. Fluid dynamic situations are then constructed to explain the effect of unsteadiness, including a physical interpretation of how stagnation pressure variations are created by temporal variations in static pressure; in this it is shown that the unsteady static pressure plays the role of a time-dependent body force potential. It is further shown that when the unsteadiness is due to a spatial nonuniformity translating at constant speed, as in a turbomachine, the unsteady pressure variation can be viewed as a local power input per unit mass from this body force to the fluid particle instantaneously at that point. © 2012 American Society of Mechanical Engineers.

Heat source modelling and natural ventilation efficiency

We compare natural ventilation flows established by a range of heat source distributions at floor level. Both evenly distributed and highly localised line and point source distributions are considered. We demonstrate that modelling the ventilation flow driven by a uniformly distributed heat source is equivalent to the flow driven by a large number of localised sources. A model is developed for the transient flow development in a room with a uniform heat distribution and is compared with existing models for localised buoyancy inputs. For large vent areas the flow driven by localised heat sources reaches a steady state more rapidly than the uniformly distributed case. For small vent areas there is little difference in the transient development times. Our transient model is then extended to consider the time taken to flush a neutrally buoyant pollutant from a naturally ventilated room. Again comparisons are drawn between uniform and localised (point and line) heat source geometries. It is demonstrated that for large vent areas a uniform heat distribution provides the fastest flushing. However, for smaller vent areas, localised heat sources produce the fastest flushing. These results are used to suggest a definition for the term 'natural ventilation efficiency', and a model is developed to estimate this efficiency as a function of the room and heat source geometries. © 2006 Elsevier Ltd. All rights reserved.

The effect of cationic, non-ionic and amphiphilic surfactants on the intercalation of bentonite

Surfactant-clay interactions are key for the development of new clay applications and inorganic-organic nanocomposites. Bentonite, with montmorillonite as the principal clay mineral constituent, was modified with varying concentrations of hexadecethyltrimethylammonium chloride (HDTMA), as a reference cationic surfactant, polypropylene glycol (PPG) 1200 and 2000, as non-ionic surfactants, and lecithin and Topcithin®, as amphiphilic phospholipid surfactants, according to the cation exchange capacity (CEC). The modified bentonites were characterised by X-ray diffraction, thermogravimetric analysis (TGA), Fourier transform infrared (FTIR) spectrometry, specific surface area and pore volume. Three intercalation regions have been identified depending on the surfactant. The non-ionic surfactant caused only a crystalline expansion of bentonite interlayers, while the cationic surfactant induced an osmotic intercalation. The amphiphilic lecithin derivatives intercalated more extensively with the bentonite matrix. The TGA and the FTIR spectra showed that, at lower concentrations, the PPGs and HDTMA adopted a disordered conformation that required more energy to degrade, while at higher concentrations, the surfactants were ordered in the interlayer space of the bentonite. The lecithin derivative surfactant had a greater thermal and conformation stability. The specific surface area reduced with increasing surfactant concentrations. This study highlights the effect of surfactant type on the interlayer space of montmorillonite in the perspective of developing novel clay functions. © 2013.

Novel fluid grid and voidage calculation techniques for a discrete element model of a 3D cylindrical fluidized bed

A discrete element model (DEM) combined with computational fluid dynamics (CFD) was developed to model particle and fluid behaviour in 3D cylindrical fluidized beds. Novel techniques were developed to (1) keep fluid cells, defined in cylindrical coordinates, at a constant volume in order to ensure the conditions for validity of the volume-averaged fluid equations were satisfied and (2) smoothly and accurately measure voidage in arbitrarily shaped fluid cells. The new technique for calculating voidage was more stable than traditional techniques, also examined in the paper, whilst remaining computationally-effective. The model was validated by quantitative comparison with experimental results from the magnetic resonance imaging of a fluidised bed analysed to give time-averaged particle velocities. Comparisons were also made between theoretical determinations of slug rise velocity in a tall bed. It was concluded that the DEM-CFD model is able to investigate aspects of the underlying physics of fluidisation not readily investigated by experiment. © 2014 The Authors.

CO2-gasification of a lignite coal in the presence of an iron-based oxygen carrier for chemical-looping combustion

Chemical-looping combustion (CLC) has the inherent property of separating the product CO2 from flue gases. Instead of air, it uses an oxygen carrier, usually in the form of a metal oxide, to provide oxygen for combustion. All techniques so far proposed for chemical looping with solid fuels involve initially the gasification of the solid fuel in order for the gaseous products to react with the oxygen carrier. Here, the rates of gasification of coal were compared when gasification was undertaken in a fluidised bed of either (i) an active Fe-based oxygen carrier used for chemical looping or (ii) inert sand. This enabled an examination of the ability of chemical looping materials to enhance the rate of gasification of solid fuels. Batch gasification and chemical-looping combustion experiments with a German lignite and its char are reported, using an electrically-heated fluidised bed reactor at temperatures from 1073 to 1223 K. The fluidising gas was CO2 in nitrogen. The kinetics of the gasification were found to be significantly faster in the presence of the oxygen carrier, especially at temperatures above 1123 K. A numerical model was developed to account for external and internal mass transfer and for the effect of the looping agent. The model also included the effects of the evolution of the pore structure at different conversions. The presence of Fe2O3 led to an increase in the rate of gasification because of the rapid oxidation of CO by the oxygen carrier to CO2. This resulted in the removal of CO and maintained a higher mole fraction of CO2 in the mixture of gas around the particle of char, i.e. within the mass transfer boundary layer surrounding the particle. This effect was most prominent at about 20% conversion when (i) the surface area for reaction was at its maximum and (ii) because of the accompanying increase in porosity and pore size, intraparticle resistance to gas mass transfer within the particle of char had fallen, compared with that in the initial particle. Excellent agreement was observed between the rates predicted by the numerical model and those observed experimentally. ©2013 Elsevier Ltd. All rights reserved.

Elastodynamic erosion of thermal barrier coatings

The finite element method is used to analyze the elastodynamic response of a columnar thermal barrier coating due to normal impact and oblique impact by an erosive particle. An assessment is made of the erosion by crack growth from preexisting flaws at the edge of each column: it is demonstrated that particle impacts can be sufficiently severe to give rise to columnar cracking. First, the transient stress state induced by the normal impact of a circular cylinder or a sphere is calculated in order to assess whether a 2D calculation adequately captures the more realistic 3D behavior. It is found that the transient stress states for the plane strain and axisymmetric models are similar. The sensitivity of response to particle diameter and to impact velocity is determined for both the cylinder and the sphere. Second, the transient stress state is explored for 2D oblique impact by a circular cylindrical particle and by an angular cylindrical particle. The sensitivity of transient tensile stress within the columns to particle shape (circular and angular), impact angle, impact location, orientation of the angular particle, and to the level of friction is explored in turn. The paper concludes with an evaluation of the effect of inclining the thermal barrier coating columns upon their erosion resistance. © 2011 The American Ceramic Society.