Predicting active material utilisation in LiFePO4 electrodes using a multi-scale mathematical model

A mathematical model is developed to simulate the discharge of a LiFePO4 cathode. This model contains 3 size scales, which match with experimental observations present in the literature on the multi-scale nature of LiFePO4 material. A shrinking-core is used on the smallest scale to represent the phase-transition of LiFePO4 during discharge. The model is then validated against existing experimental data and this validated model is then used to investigate parameters that influence active material utilisation. Specifically the size and composition of agglomerates of LiFePO4 crystals is discussed, and we investigate and quantify the relative effects that the ionic and electronic conductivities within the oxide have on oxide utilisation. We find that agglomerates of crystals can be tolerated under low discharge rates. The role of the electrolyte in limiting (cathodic) discharge is also discussed, and we show that electrolyte transport does limit performance at high discharge rates, confirming the conclusions of recent literature.

A comparison of mathematical models for phase separation in high-rate LiFePO4 cathodes

We construct a two-scale mathematical model for modern, high-rate LiFePO4cathodes. We attempt to validate against experimental data using two forms of the phase-field model developed recently to represent the concentration of Li+ in nano-sized LiFePO4crystals. We also compare this with the shrinking-core based model we developed previously. Validating against high-rate experimental data, in which electronic and electrolytic resistances have been reduced is an excellent test of the validity of the crystal-scale model used to represent the phase-change that may occur in LiFePO4material. We obtain poor fits with the shrinking-core based model, even with fitting based on “effective” parameter values. Surprisingly, using the more sophisticated phase-field models on the crystal-scale results in poorer fits, though a significant parameter regime could not be investigated due to numerical difficulties. Separate to the fits obtained, using phase-field based models embedded in a two-scale cathodic model results in “many-particle” effects consistent with those reported recently.

Kinetics of the chemical looping oxidation of H2 by a co-precipitated mixture of CuO and Al2O3

Chemical-looping combustion (CLC) has the inherent property of separating 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. When used for the combustion of gaseous fuels, such as natural gas, or synthesis gas from the gasification of coal, the technique gives a stream of CO2 which, on an industrial scale, would be sufficiently pure for geological sequestration. An important issue is the form of the metal oxide, since it must retain its reactivity through many cycles of complete reduction and oxidation. Here, we report on the rates of oxidation of one constituent of synthesis gas, H2, by co-precipitated mixtures of CuO+Al2O3 using a laboratory-scale fluidised bed. To minimise the influence of external mass transfer, and also of errors in the measurement of [H2], particles sized to 355-500μm were used at low [H2], with the temperature ranging from 450 to 900°C. Under such conditions, the reaction was slow enough for meaningful measurements of the intrinsic kinetics to be made. The reaction was found to be first order with respect to H2. Above ∼800°C, the reaction of CuO was fast and conformed to the shrinking core mechanism, proceeding via the intermediate, Cu2O, in: 2CuO+H2→Cu2O+H2O, ΔH1073 K0=- 116.8 kJ/mol; Cu2O+H2→2Cu+H2O, ΔH1073 K0-80.9 kJ/mol. After oxidation of the products Cu and Cu2O back to CuO, the kinetics in subsequent cycles of chemical looping oxidation of H2 could be approximated by those in the first. Interestingly, the carrier was found to react at temperatures as low as 300°C. The influence of the number of cycles of reduction and oxidation is explored. Comparisons are drawn with previous work using reduction by CO. Finally, these results indicate that the kinetics of reaction of the oxygen carrier with gasifier synthesis gases is very much faster than rates of gasification of the original fuel. © 2010 The Institution of Chemical Engineers.

Study on the Dissociation Behavior of Gas Hydrate in Porous Sediment Based on Fractal Theory

The dissociation process of gas hydrate was regarded as a gas-solid reaction without solid production layer when the temperature was above the zero centigrade. Based on the shrinking core model and the fractal theory, a fractional dimension dynamical model for gas hydrate dissociation in porous sediment was established. The new approach of evaluating the fractal dimension of the porous media was also presented. The fractional dimension dynamical model for gas hydrate dissociation was examined with the previous experimental data of methane hydrate and carbon dioxide hydrate dissociations, respectively. The calculated results indicate that the fractal dimensions of porous media acquired with this method agree well with the previous study. With the absolute average deviation (AAD) below 10%, the present model provided satisfactory predictions for the dissociation process of methane hydrate and carbon dioxide hydrate.

Post consumer pet depolymerization by acid hydrolysis

The reaction of post-consumer poly(ethylene terephthalate) with aqueous solutions of sulfuric acid 7.5M was investigated in terms of temperature, time and particle size. The reaction extent reached 80% in four days at 100 degrees C and 90% in 5 hours at 135 degrees C. TPA obtained was purified and considered in the same level of quality of the commercial one after tests of elemental analysis, particle size and color. It was concluded that the hydrolysis occurred preferentially at the chain ends and superficially, having as controller mechanism the acid diffusion into the polymer structure. The shrinking-core model can explain the reaction kinetics.

Hydrometallurgical process modeling using advanced mathematical tools

Hydrometallurgical process modeling is the main objective of this Master’s thesis work. Three different leaching processes namely, high pressure pyrite oxidation, direct oxidation zinc concentrate (sphalerite) leaching and gold chloride leaching using rotating disc electrode (RDE) are modeled and simulated using gPROMS process simulation program in order to evaluate its model building capabilities. The leaching mechanism in each case is described in terms of a shrinking core model. The mathematical modeling carried out included process model development based on available literature, estimation of reaction kinetic parameters and assessment of the model reliability by checking the goodness fit and checking the cross correlation between the estimated parameters through the use of correlation matrices. The estimated parameter values in each case were compared with those obtained using the Modest simulation program. Further, based on the estimated reaction kinetic parameters, reactor simulation and modeling for direct oxidation zinc concentrate (sphalerite) leaching is carried out in Aspen Plus V8.6. The zinc leaching autoclave is based on Cominco reactor configuration and is modeled as a series of continuous stirred reactors (CSTRs). The sphalerite conversion is calculated and a sensitivity analysis is carried out so to determine the optimum reactor operation temperature and optimum oxygen mass flow rate. In this way, the implementation of reaction kinetic models into the process flowsheet simulation environment has been demonstrated.

Post consumer pet depolymerization by acid hydrolysis

The reaction of post-consumer poly(ethylene terephthalate) with aqueous solutions of sulfuric acid 7.5M was investigated in terms of temperature, time and particle size. The reaction extent reached 80% in four days at 100 degrees C and 90% in 5 hours at 135 degrees C. TPA obtained was purified and considered in the same level of quality of the commercial one after tests of elemental analysis, particle size and color. It was concluded that the hydrolysis occurred preferentially at the chain ends and superficially, having as controller mechanism the acid diffusion into the polymer structure. The shrinking-core model can explain the reaction kinetics.

Rock eval analysis of organic matter in sediment core 159-959D

The Cretaceous Equatorial Atlantic Gateway between the Central and South Atlantic basins is of interest not only for paleoceanographic and paleoclimatic studies, but also because it provided particularly favourable conditions for the accumulation and preservation of organic-rich sediments. Deposition of carbonaceous sediments along the Côte d'Ivoire-Ghana Transform Margin (Ocean Drilling Program Leg 159) was intimately linked to the plate tectonic and paleoceanographic evolution of this gateway. Notably, the formation of a marginal basement ridge on the southeastern border of the transform margin provided an efficient shelter of the landward Deep Ivorian Basin against erosive and potentially oxidizing currents. Different subsidence histories across the transform margin were responsible for the development of distinct depositional settings on the crest and on both sides of the basement ridge. Whereas the southern, oceanward flank of the basement ridge was characterized by rapid, continuous deepening since last Albian-early Cenomanian, marine sedimentation on the northern, landward flank was interrupted by a period of uplift and erosion in the late Albian, and rapid subsidence started after the early Coniacian. Organic-rich sediments occur throughout almost the entire Cretaceous section, but hydrogen-rich marine black shales were exclusively recovered from core sections above an uplift-related unconformity. These black shales formed when separation of Africa and South America was sufficient to allow permanent oceanic midwater exchange after the late Albian. Four periods of black shale accumulation are recovered, some of them are correlated with the global oceanic anoxic events: in the last Albian-earliest Cenomanian, at the Cenomanian-Turronian boundary, during the middle Coniacian-early Campanian, and in the mid-Maastrichtian. These periods were characterized by increasing carbon flux to the seafloor, induced by enhanced palaeoproductivity and intensified supply of terrestrial organic matter. Black shale depostion appears to be intimately linked to periods of rising or maximum eustatic sea level and to the expansion of the oxygen minimum zone, as indicated by foraminiferal biofacies. Intervals between black shales units, in contrast, indicate a shrinking oxygen minimum zone and enhanced detrital flux rates, probably related to lowering sea level. Upper Cretaceous detritral limestones with high porosities may provide excellent hydrocarbon reservoirs, alsthough their areal extent appears to be limited. Palaeogene porcellanites, capped by Neogene pelagic marls and clays, extend over a wider area and max provide another target for hydrocarbon exploration.

Progenitor-Dependent Explosion Dynamics in Self-Consistent, Axisymmetric Simulations of Neutrino-Driven Core-Collapse Supernovae

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11–28 M ⊙, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si–O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si–O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection timescales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.