Material model of lung parenchyma based on living precision-cut lung slice testing.

We describe a novel constitutive model of lung parenchyma, which can be used for continuum mechanics based predictive simulations. To develop this model, we experimentally determined the nonlinear material behavior of rat lung parenchyma. This was achieved via uni-axial tension tests on living precision-cut rat lung slices. The resulting force-displacement curves were then used as inputs for an inverse analysis. The Levenberg-Marquardt algorithm was utilized to optimize the material parameters of combinations and recombinations of established strain-energy density functions (SEFs). Comparing the best-fits of the tested SEFs we found Wpar = 4.1 kPa(I1-3)2 + 20.7 kPa(I1 - 3)3 + 4.1 kPa(-2 ln J + J2 - 1) to be the optimal constitutive model. This SEF consists of three summands: the first can be interpreted as the contribution of the elastin fibers and the ground substance, the second as the contribution of the collagen fibers while the third controls the volumetric change. The presented approach will help to model the behavior of the pulmonary parenchyma and to quantify the strains and stresses during ventilation.

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.

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.

Initiation of immersed granular avalanches.

By means of coupled molecular dynamics-computational fluid dynamics simulations, we analyze the initiation of avalanches in a granular bed of spherical particles immersed in a viscous fluid and inclined above its angle of repose. In quantitative agreement with experiments, we find that the bed is unstable for a packing fraction below 0.59 but is stabilized above this packing fraction by negative excess pore pressure induced by the effect of dilatancy. From detailed numerical data, we explore the time evolution of shear strain, packing fraction, excess pore pressures, and granular microstructure in this creeplike pressure redistribution regime, and we show that they scale excellently with a characteristic time extracted from a model based on the balance of granular stresses in the presence of a negative excess pressure and its interplay with dilatancy. The cumulative shear strain at failure is found to be ≃ 0.2, in close agreement with the experiments, irrespective of the initial packing fraction and inclination angle. Remarkably, the avalanche is triggered when dilatancy vanishes instantly as a result of fluctuations while the average dilatancy is still positive (expanding bed) with a packing fraction that declines with the initial packing fraction. Another nontrivial feature of this creeplike regime is that, in contrast to dry granular materials, the internal friction angle of the bed at failure is independent of dilatancy but depends on the inclination angle, leading therefore to a nonlinear dependence of the excess pore pressure on the inclination angle. We show that this behavior may be described in terms of the contact network anisotropy, which increases with a nearly constant connectivity and levels off at a value (critical state) that increases with the inclination angle. These features suggest that the behavior of immersed granular materials is controlled not only directly by hydrodynamic forces acting on the particles but also by the influence of the fluid on the granular microstructure.

Non-catalytic multi-fuel reformer for syngas production for fuel cell applications

Rich combustion of n-heptane, diesel oil, jet A-1 kerosene, and bio-diesel (rapeseed-oil methyl ester) were studied to produce hydrogen enriched gas, ready for the cleanup stages for fuel cell applications. n-heptane was successfully reformed up to an equivalence ratio of 3:1, reaching a conversion efficiency up to 83% for a packed bed of alumina bead burner. Diesel, kerosene and bio-diesel were reformed to synthesis gas with conversion efficiency up to 65%. At equivalence ratio of 2:1 and P=7 kw, stability, low HC formation, high conversion efficiency, and low soot emission were achieved. A common synthesis gas composition around this condition was 15 and 13% H2, 15 and 17% CO, and 4 and 4.5% CO2 for n-heptane and diesel, jet A-1 and bio-diesel, respectively, for burner A. This is an abstract of a paper presented at the 2010 Spring National Meeting (San Antonio, TX 3/21-25/2010).

Metal-polymer composite sensors for volatile organic compounds: Part 1. Flow-through chemi-resistors

A new type of chemi-resistor based on a novel metal-polymer composite is described. The composite contains nickel particles with sharp nano-scale surface features, which are intimately coated by the polymer matrix so that they do not come into direct physical contact. No conductive chains of filler particles are formed even at loadings above the percolation threshold and the composite is intrinsically insulating. However, when subjected to compression the composite becomes conductive, with sample resistance falling from ≥ 1012 Ω to < 0.01 Ω. The composite can be formed into insulating granules, which display similar properties to the bulk form. A bed of granules compressed between permeable frits provides a porous structure with a start resistance set by the degree of compression while the granules are free to swell when exposed to volatile organic compounds (VOCs). The granular bed presents a large surface area for the adsorption of VOCs from the gas stream flowing through it. The response of this system to a variety of vapours has been studied for two different sizes of the granular bed and for different matrix polymers. Large responses, ΔR/R0 ≥ 10^7, are observed when saturated vapours are passed through the chemi-resistor. Rapid response allows real time sensing of VOCs and the initial state is recovered in a few seconds by purging with an inert gas stream. The variation in response as a function of VOC concentration is determined.

Interaction between Fe-based oxygen carriers and n-heptane during chemical looping combustion

Chemical looping combustion (CLC) uses a metal oxide (the oxygen carrier) to provide oxygen for the combustion of a fuel and gives an inherent separation of pure CO2 with minimal energy penalty. In solid-fuel CLC, volatile matter will interact with oxygen carriers. Here, the interaction between iron-based oxygen carriers and a volatile hydrocarbon (n-heptane) was investigated in both a laboratory-scale fluidised bed and a thermogravimetric analyser (TGA). Experiments were undertaken to characterise the thermal decomposition of the n-heptane occurring in the presence and in the absence of the oxygen carrier. In a bed of inert particles, carbon deposition increased with temperature and acetylene appeared as a possible precursor. For a bed of carrier consisting of pure Fe2O3, carbon deposition occurred once the Fe2O3 was fully reduced to Fe. When the Fe2O3 was doped with 10 mol % Al2O3 (Fe90Al), deposition started when the carrier was reduced to a mixture of Fe and FeAl2O4, the latter being very unreactive. Furthermore, when pure Fe2O3 was fully reduced to Fe, agglomeration of the fluidised bed occurred. However, Fe90Al did not give agglomeration even after extended reduction. The results suggest that Fe90Al is promising for the CLC of solid fuels. © 2012 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Seafloor-riser interaction model

Fatigue stresses associated with extreme storms, vessel movements, and vortex-induced vibrations are critical to the performance of steel catenary risers. The critical location for fatigue damage often occurs within the touchdown zone, where cyclic interaction of the riser with the seabed occurs. Developing a model for seabed stiffness requires characterization of a number of complex nonlinear processes including trench formation, nonlinear soil stiffness, soil suction, and breakaway of the riser from the seafloor. The analytical framework utilized in this research considers the riser-seafloor interaction problem in terms of a pipe resting on a bed of springs, the stiffness characteristics of which are described by nonlinear load-deflection (P-y) curves. The P-y model allows for first penetration and uplift, as well as repenetration and small range motions within the bounding loop defined by extreme loading. The backbone curve is constructed from knowledge of the soil strength, the rate of strength increase with depth, trench width, and two additional parameters, while three parameters are necessary for the cyclic response. © ASCE 2009.

Seafloor interaction with steel catenary risers

This paper describes the key features of a seafloor-riser interaction model. The soil is represented in terms of non-linear load-deflection (P- y) relationships, which are also able to account for soil stiffness degradation due to cyclic loading. The analytical framework considers the riser-seafloor interaction problem in terms of a pipe resting on a bed of springs, and requires the iterative solution of a fourth-order ordinary differential equation. A series of simulations is used to illustrate the capabilities of the model. Thanks to the non-linear soil springs with stiffness degradation it is possible to simulate the trench formation process and estimate moments in a riser. Copyright © 2008 by The International Society of Offshore and Polar Engineers (ISOPE).