CFD modelling of the fast pyrolysis of an in-flight cellulosic particle subjected to convective heat transfer

The pyrolysis of a freely moving cellulosic particle inside a 41.7mgs -1 source continuously fed fluid bed reactor subjected to convective heat transfer is modelled. The Lagrangian approach is adopted for the particle tracking inside the reactor, while the flow of the inert gas is treated with the standard Eulerian method for gases. The model incorporates the thermal degradation of cellulose to char with simultaneous evolution of gases and vapours from discrete cellulosic particles. The reaction kinetics is represented according to the Broido–Shafizadeh scheme. The convective heat transfer to the surface of the particle is solved by two means, namely the Ranz–Marshall correlation and the limit case of infinitely fast external heat transfer rates. The results from both approaches are compared and discussed. The effect of the different heat transfer rates on the discrete phase trajectory is also considered.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors. Part B: Heat, momentum and mass transport in bubbling fluidised beds

The fluid–particle interaction inside a 150 g/h fluidised bed reactor is modelled. The biomass particle is injected into the fluidised bed and the heat, momentum and mass transport from the fluidising gas and fluidised sand is modelled. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Heat transfer from the bubbling bed to the discrete biomass particle, as well as biomass reaction kinetics are modelled according to the literature. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase. FLUENT 6.2 has been used as the modelling framework of the simulations with the whole pyrolysis model incorporated in the form of user-defined function (UDF). The study completes the fast pyrolysis modelling in bubbling fluidised bed reactors.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors: modelling the impact of biomass shrinkage

The fluid–particle interaction and the impact of shrinkage on pyrolysis of biomass inside a 150 g/h fluidised bed reactor is modelled. Two 500 View the MathML sourcem in diameter biomass particles are injected into the fluidised bed with different shrinkage conditions. The two different conditions consist of (1) shrinkage equal to the volume left by the solid devolatilization, and (2) shrinkage parameters equal to approximately half of particle volume. The effect of shrinkage is analysed in terms of heat and momentum transfer as well as product yields, pyrolysis time and particle size considering spherical geometries. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Heat transfer from the bubbling bed to the discrete biomass particle, as well as biomass reaction kinetics are modelled according to the literature. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase. FLUENT 6.2 has been used as the modelling framework of the simulations with the whole pyrolysis model incorporated in the form of user defined function (UDF).

Numerical Comparison of the Drag Models of Granular Flows Applied to the Fast Pyrolysis of Biomass

The paper presents a comparison between the different drag models for granular flows developed in the literature and the effect of each one of them on the fast pyrolysis of wood. The process takes place on an 100 g/h lab scale bubbling fluidized bed reactor located at Aston University. FLUENT 6.3 is used as the modeling framework of the fluidized bed hydrodynamics, while the fast pyrolysis of the discrete wood particles is incorporated as an external user defined function (UDF) hooked to FLUENT’s main code structure. Three different drag models for granular flows are compared, namely the Gidaspow, Syamlal O’Brien, and Wen-Yu, already incorporated in FLUENT’s main code, and their impact on particle trajectory, heat transfer, degradation rate, product yields, and char residence time is quantified. The Eulerian approach is used to model the bubbling behavior of the sand, which is treated as a continuum. Biomass reaction kinetics is modeled according to the literature using a two-stage, semiglobal model that takes into account secondary reactions.

Fast x-ray spectroscopy study of ethene on clean and SO4 precovered Pt{111}

The adsorption and decomposition of ethylene over a Pt{111} single crystalsurface has been investigated by fast x-ray spectroscopy. At 100 K ethene displays precursor-mediated adsorption kinetics, adopting a single environment with a saturation C2H4 coverage of 0.25 ML and binding energy of 283.2 eV. Thermal decomposition proceeds above 240 K via dehydrogenation to ethylidyne with an activation barrier of 57±3 kJ mol−1 and preexponential factor ν=1×1010±0.5 s−1. Site-blocking by preadsorbed SO4 reduces the saturation ethene coverage but induces a new, less reactive π-bonded ethene species centered around 283.9 eV, which in turn decomposes to ethylidyne at 350 K.

Pseudoperiodic "Living" and/or controlled cationic ring-opening copolymerization of oxetane with tetrahydropyran:microstructure of polymers vs kinetics of chain growth

The "living" and/or controlled cationic ring-opening bulk copolymerization of oxetane (Ox) with tetrahydropyran (THP) (cyclic ether with no homopolymerizability) at 35°C was examined using ethoxymethyl-1 -oxoniacyclohexane hexafluoroantimonate (EMOA) and (BF3 · CH3OH)THP as fast and slow initiator, respectively, yielding living and nonliving polymers with pseudoperiodic sequences (i.e., each pentamethylene oxide fragment inserted into the polymer is flanked by two trimethylene oxide fragments). Good control over number-average molecular weight (Mn up to 150000 g mol-1) with molecular weight distribution (MWD ∼ 1.4-1, 5) broader than predicted by the Poison distribution (MWDs > 1 +1/DPn) was attained using EMOA as initiating system, i.e., C 2H5OCH2Cl with 1.1 equiv of AgSbF6 as a stable catalyst and 1.1 equiv of 2,6-di-tert-butylpyridine used as a non-nucleophilic proton trap. With (BF3 · CH 3OH)THP, a drift of the linear dependence M n(GPC) vs Mn(theory) to lower molecular weight was observed together with the production of cyclic oligomers, ∼3-5% of the Ox consumed in THP against ∼30% in dichloromethane. Structural and kinetics studies highlighted a mechanism of chains growth where the rate of mutual conversion between "strain ACE species" (chain terminated by a tertiary 1-oxoniacyclobutane ion, Al) and "strain-free ACE species" (chain terminated by a tertiary 1-oxoniacyclohexane ion, Tl) depends on the rate at which Ox converts the stable species T1 (kind of "dormant" species) into a living "propagating" center A1 (i.e., k aapp[Ox]). The role of the THP solvent associated with the suspension of irreversible and reversible transfer reactions to polymer, when the polymerization is initiated with EMOA, was predicted by our kinetic considerations. The activation -deactivation pseudoequilibrium coefficient (Qt) was then calculated in a pure theoretical basis. From the measured apparent rate constant of Ox (kOxapp) and THP (kTHPapp = ka(endo)app) consumption, Qt and reactivity ratio (kp/kd, k a(endo)/ka(exo), and ks/ka(endo) were calculated, which then allow the determination of the transition rate constant of elementary step reactions that governs the increase of Mu with conversion. © 2009 American Chemical Society.

Influence of reaction atmosphere (H2O, N2, H2, CO2, CO) on fluidized-bed fast pyrolysis of biomass using detailed tar vapor chemistry in computational fluid dynamics

Secondary pyrolysis in fluidized bed fast pyrolysis of biomass is the focus of this work. A novel computational fluid dynamics (CFD) model coupled with a comprehensive chemistry scheme (134 species and 4169 reactions, in CHEMKIN format) has been developed to investigate this complex phenomenon. Previous results from a transient three-dimensional model of primary pyrolysis were used for the source terms of primary products in this model. A parametric study of reaction atmospheres (H2O, N2, H2, CO2, CO) has been performed. For the N2 and H2O atmosphere, results of the model compared favorably to experimentally obtained yields after the temperature was adjusted to a value higher than that used in experiments. One notable deviation versus experiments is pyrolytic water yield and yield of higher hydrocarbons. The model suggests a not overly strong impact of the reaction atmosphere. However, both chemical and physical effects were observed. Most notably, effects could be seen on the yield of various compounds, temperature profile throughout the reactor system, residence time, radical concentration, and turbulent intensity. At the investigated temperature (873 K), turbulent intensity appeared to have the strongest influence on liquid yield. With the aid of acceleration techniques, most importantly dimension reduction, chemistry agglomeration, and in-situ tabulation, a converged solution could be obtained within a reasonable time (∼30 h). As such, a new potentially useful method has been suggested for numerical analysis of fast pyrolysis.

Fast-Food Franchises: An Alternative Menu for Hotel/Casinos

In their discussion - Fast-Food Franchises: An Alternative Menu for Hotel/Casinos - by Skip Swerdlow, Assistant Professor of Finance, Larry Strate, Assistant Professor of Business Law, and Francis X. Brown, Assistant Professor of Hotel Administration at the University of Nevada, Las Vegas, their preview reads: Hotel/casino food service operations are adding some non-traditional fare to their daily offerings in the form of fast-food franchises. The authors review aspects of franchising and cite some new Las Vegas food ideas.” The authors offer that the statewide food and beverage figures, according to the Nevada Gaming Abstract of 1985, exceeded $1.24 billion. Most of that figure was generated in traditional coffee shops, gourmet dining rooms, and buffets. With that kind of food and beverage figure solidly on the table, it was inevitable that fast-food franchises would move into casinos to garner a share of the proceeds. In a March 1986 review of franchising, Restaurant Business reported the following statistics: “Over 60 percent of all restaurants are franchisee owned. This relationship is also paralleled in dollar sales, which has exceeded $53 billion.” “Restaurant franchising expansion has grown at an annual rate of 12 percent per year for the past five years.” The beginning of the article is dedicated to describing, in general, the franchise phenomenon; growth has been spectacular the authors inform you. “The franchise concept has provided an easy method of going into business for the entrepreneur with minimal business experience, but a desire to work hard to make a profit,” say professors Swerdlow, Strate, and Brown. Lured by tourist traffic, and the floundering Chapter 11 afflicted, Riviera Hotel and Casino in Las Vegas, Burger King saw an attractive opportunity for an experiment in non-traditional outlet placement, say the authors. Although innately transient, the tourist numbers were way too significant to ignore. That tourist traffic, the authors say, is ‘round-the-clock. Added to that figure is the 2000-3000 average employee count for many of the casinos on the ‘Vegas strip. Not surprisingly, the project began to look very appealing to both Burger King and the Riviera Hotel/Casino, the authors report. In the final analysis, the project did work out well; very well indeed. So it is written, “The successful operation of the Burger King in the Riviera has sparked interest by other existing hotel/casino operations and fast-food restaurant chains. Burger King's operation, like so many other industry leadership decisions, provides impetus for healthy competition in a market that is burgeoning not only because of expansion that recognizes traditional population growth, but because of bold moves that search for customers in non-traditional areas.” The authors provide an Appendix listing Las Vegas hotel/casino properties and the restaurants they contain.

Effects of Natural Organic Matter on the Dissolution Kinetics and Bioavailability of Metal Oxide Nanoparticles

The rapid development of nanotechnology and wider applications of engineered nanomaterials (ENMs) in the last few decades have generated concerns regarding their environmental and health risks. After release into the environment, ENMs undergo aggregation, transformation, and, for metal-based nanomaterials, dissolution processes, which together determine their fate, bioavailability and toxicity to living organisms in the ecosystems. The rates of these processes are dependent on nanomaterial characteristics as well as complex environmental factors, including natural organic matter (NOM). As a ubiquitous component of aquatic systems, NOM plays a key role in the aggregation, dissolution and transformation of metal-based nanomaterials and colloids in aquatic environments.

The goal of this dissertation work is to investigate how NOM fractions with different chemical and molecular properties affect the dissolution kinetics of metal oxide ENMs, such as zinc oxide (ZnO) and copper oxide (CuO) nanoparticles (NPs), and consequently their bioavailability to aquatic vertebrate, with Gulf killifish (Fundulus grandis) embryos as model organisms.

ZnO NPs are known to dissolve at relatively fast rates, and the rate of dissolution is influenced by water chemistry, including the presence of Zn-chelating ligands. A challenge, however, remains in quantifying the dissolution of ZnO NPs, particularly for time scales that are short enough to determine rates. This dissertation assessed the application of anodic stripping voltammetry (ASV) with a hanging mercury drop electrode to directly measure the concentration of dissolved Zn in ZnO NP suspensions, without separation of the ZnO NPs from the aqueous phase. Dissolved zinc concentration measured by ASV ([Zn]ASV) was compared with that measured by inductively coupled plasma mass spectrometry (ICP-MS) after ultracentrifugation ([Zn]ICP-MS), for four types of ZnO NPs with different coatings and primary particle diameters. For small ZnO NPs (4-5 nm), [Zn]ASV was 20% higher than [Zn]ICP-MS, suggesting that these small NPs contributed to the voltammetric measurement. For larger ZnO NPs (approximately 20 nm), [Zn]ASV was (79±19)% of [Zn]ICP-MS, despite the high concentrations of ZnO NPs in suspension, suggesting that ASV can be used to accurately measure the dissolution kinetics of ZnO NPs of this primary particle size.

Using the ASV technique to directly measure dissolved zinc concentration, we examined the effects of 16 different NOM isolates on the dissolution kinetics of ZnO NPs in buffered potassium chloride solution. The observed dissolution rate constants (kobs) and dissolved zinc concentrations at equilibrium increased linearly with NOM concentration (from 0 to 40 mg-C L-1) for Suwannee River humic acid (SRHA), Suwannee River fulvic acid and Pony Lake fulvic acid. When dissolution rates were compared for the 16 NOM isolates, kobs was positively correlated with certain properties of NOM, including specific ultraviolet absorbance (SUVA), aromatic and carbonyl carbon contents, and molecular weight. Dissolution rate constants were negatively correlated to hydrogen/carbon ratio and aliphatic carbon content. The observed correlations indicate that aromatic carbon content is a key factor in determining the rate of NOM-promoted dissolution of ZnO NPs. NOM isolates with higher SUVA were also more effective at enhancing the colloidal stability of the NPs; however, the NOM-promoted dissolution was likely due to enhanced interactions between surface metal ions and NOM rather than smaller aggregate size.

Based on the above results, we designed experiments to quantitatively link the dissolution kinetics and bioavailability of CuO NPs to Gulf killifish embryos under the influence of NOM. The CuO NPs dissolved to varying degrees and at different rates in diluted 5‰ artificial seawater buffered to different pH (6.3-7.5), with or without selected NOM isolates at various concentrations (0.1-10 mg-C L-1). NOM isolates with higher SUVA and aromatic carbon content (such as SRHA) were more effective at promoting the dissolution of CuO NPs, as with ZnO NPs, especially at higher NOM concentrations. On the other hand, the presence of NOM decreased the bioavailability of dissolved Cu ions, with the uptake rate constant negatively correlated to dissolved organic carbon concentration ([DOC]) multiplied by SUVA, a combined parameter indicative of aromatic carbon concentration in the media. When the embryos were exposed to CuO NP suspension, changes in their Cu content were due to the uptake of both dissolved Cu ions and nanoparticulate CuO. The uptake rate constant of nanoparticulate CuO was also negatively correlated to [DOC]×SUVA, in a fashion roughly proportional to changes in dissolved Cu uptake rate constant. Thus, the ratio of uptake rate constants from dissolved Cu and nanoparticulate CuO (ranging from 12 to 22, on average 17±4) were insensitive to NOM type or concentration. Instead, the relative contributions of these two Cu forms were largely determined by the percentage of CuO NP that was dissolved.

Overall, this dissertation elucidated the important role that dissolved NOM plays in affecting the environmental fate and bioavailability of soluble metal-based nanomaterials. This dissertation work identified aromatic carbon content and its indicator SUVA as key NOM properties that influence the dissolution, aggregation and biouptake kinetics of metal oxide NPs and highlighted dissolution rate as a useful functional assay for assessing the relative contributions of dissolved and nanoparticulate forms to metal bioavailability. Findings of this dissertation work will be helpful for predicting the environmental risks of engineered nanomaterials.

Study of the kinetics and mechanism of rapid self-assembly in block copolymer thin films during solvo-microwave annealing

Microwave annealing is an emerging technique for achieving ordered patterns of block copolymer films on substrates. Little is understood about the mechanisms of microphase separation during the microwave annealing process and how it promotes the microphase separation of the blocks. Here, we use controlled power microwave irradiation in the presence of tetrahydrofuran (THF) solvent, to achieve lateral microphase separation in high- lamellar-forming poly(styrene-b-lactic acid) PS-b-PLA. A highly ordered line pattern was formed within seconds on silicon, germanium and silicon on insulator (SOI) substrates. In-situ temperature measurement of the silicon substrate coupled to condition changes during "solvo-microwave" annealing allowed understanding of the processes to be attained. Our results suggest that the substrate has little effect on the ordering process and is essentially microwave transparent but rather, it is direct heating of the polar THF molecules that causes microphase separation. It is postulated that the rapid interaction of THF with microwaves and the resultant temperature increase to 55 degrees C within seconds causes an increase of the vapor pressure of the solvent from 19.8 to 70 kPa. This enriched vapor environment increases the plasticity of both PS and PLA chains and leads to the fast self-assembly kinetics. Comparing the patterns formed on silicon, germanium and silicon on insulator (SOI) and also an in situ temperature measurement of silicon in the oven confirms the significance of the solvent over the role of substrate heating during "solvo-microwave" annealing. Besides the short annealing time which has technological importance, the coherence length is on a micron scale and dewetting is not observed after annealing. The etched pattern (PLA was removed by an Ar/O-2 reactive ion etch) was transferred to the underlying silicon substrate fabricating sub-20 nm silicon nanowires over large areas demonstrating that the morphology is consistent both across and through the film.

LES of flame acceleration and DDT in hydrogen-air mixture using artificially thickened flame approach and detailed chemical kinetics

A large eddy simulation is performed to study the deflagration to detonation transition phenomenon in an obstructed channel containing premixed stoichiometric hydrogen–air mixture. Two-dimensional filtered reactive Navier–Stokes equations are solved utilizing the artificially thickened flame approach (ATF) for modeling sub-grid scale combustion. To include the effect of induction time, a 27-step detailed mechanism is utilized along with an in situ adaptive tabulation (ISAT) method to reduce the computational cost due to the detailed chemistry. The results show that in the slow flame propagation regime, the flame–vortex interaction and the resulting flame folding and wrinkling are the main mechanisms for the increase of the flame surface and consequently acceleration of the flame. Furthermore, at high speed, the major mechanisms responsible for flame propagation are repeated reflected shock–flame interactions and the resulting baroclinic vorticity. These interactions intensify the rate of heat release and maintain the turbulence and flame speed at high level. During the flame acceleration, it is seen that the turbulent flame enters the ‘thickened reaction zones’ regime. Therefore, it is necessary to utilize the chemistry based combustion model with detailed chemical kinetics to properly capture the salient features of the fast deflagration propagation.

Analysis and control of transient spectra using time-dependent density functional theory.

215 p.

Transient scaling and resurgence of chimera states in networks of Boolean phase oscillators.

We study networks of nonlocally coupled electronic oscillators that can be described approximately by a Kuramoto-like model. The experimental networks show long complex transients from random initial conditions on the route to network synchronization. The transients display complex behaviors, including resurgence of chimera states, which are network dynamics where order and disorder coexists. The spatial domain of the chimera state moves around the network and alternates with desynchronized dynamics. The fast time scale of our oscillators (on the order of 100ns) allows us to study the scaling of the transient time of large networks of more than a hundred nodes, which has not yet been confirmed previously in an experiment and could potentially be important in many natural networks. We find that the average transient time increases exponentially with the network size and can be modeled as a Poisson process in experiment and simulation. This exponential scaling is a result of a synchronization rate that follows a power law of the phase-space volume.

Ultrafast transient absorption spectroscopy: principles and applications

The study of photophysical and photochemical processes crosses the interest of many fields of research in physics, chemistry and biology. In particular, the photophysical and photochemical reactions, after light absorption by a photosynthetic pigment-protein complex, are among the fastest events in biology, taking place on timescales ranging from tens of femtoseconds to a few nanoseconds. Among the experimental approaches developed for this purpose, the advent of ultrafast transient absorption spectroscopy has become a powerful and widely used technique.[1,2] Focusing on the process of photosynthesis, it relies upon the efficient absorption and conversion of the radiant energy from the Sun. Chlorophylls and carotenoids are the main players in the process. Photosynthetic pigments are typically arranged in a highly organized fashion to constitute antennas and reaction centers, supramolecular devices where light harvesting and charge separation take place. The very early steps in the photosynthetic process take place after the absorption of a photon by an antenna system, which harvests light and eventually delivers it to the reaction center. In order to compete with internal conversion, intersystem crossing, and fluorescence, which inevitably lead to energy loss, the energy and electron transfer processes that fix the excited-state energy in photosynthesis must be extremely fast. In order to investigate these events, ultrafast techniques down to a sub-100 fs resolution must be used. In this way, energy migration within the system as well as the formation of new chemical species such as charge-separated states can be tracked in real time. This can be achieved by making use of ultrafast transient absorption spectroscopy. The basic principles of this notable technique, instrumentation, and some recent applications to photosynthetic systems[3] will be described. Acknowledgements M. Moreno Oliva thanks the MINECO for a “Juan de la Cierva-Incorporación” research contract. References [1] U. Megerle, I. Pugliesi, C. Schriever, C.F. Sailer and E. Riedle, Appl. Phys. B, 96, 215 – 231 (2009). [2] R. Berera, R. van Grondelle and J.T.M. Kennis, Photosynth. Res., 101, 105 – 118 (2009). [3] T. Nikkonen, M. Moreno Oliva, A. Kahnt, M. Muuronen, J. Helaja and D.M. Guldi, Chem. Eur. J., 21, 590 – 600 (2015).

A SCR Model based on Reactor and Engine Experimental Studies for a Cu-zeolite Catalyst

A NOx reduction efficiency higher than 95% with NH3 slip less than 30 ppm is desirable for heavy-duty diesel (HDD) engines using selective catalytic reduction (SCR) systems to meet the US EPA 2010 NOx standard and the 2014-2018 fuel consumption regulation. The SCR performance needs to be improved through experimental and modeling studies. In this research, a high fidelity global kinetic 1-dimensional 2-site SCR model with mass transfer, heat transfer and global reaction mechanisms was developed for a Cu-zeolite catalyst. The model simulates the SCR performance for the engine exhaust conditions with NH3 maldistribution and aging effects, and the details are presented. SCR experimental data were collected for the model development, calibration and validation from a reactor at Oak Ridge National Laboratory (ORNL) and an engine experimental setup at Michigan Technological University (MTU) with a Cummins 2010 ISB engine. The model was calibrated separately to the reactor and engine data. The experimental setup, test procedures including a surrogate HD-FTP cycle developed for transient studies and the model calibration process are described. Differences in the model parameters were determined between the calibrations developed from the reactor and the engine data. It was determined that the SCR inlet NH3 maldistribution is one of the reasons causing the differences. The model calibrated to the engine data served as a basis for developing a reduced order SCR estimator model. The effect of the SCR inlet NO2/NOx ratio on the SCR performance was studied through simulations using the surrogate HD-FTP cycle. The cumulative outlet NOx and the overall NOx conversion efficiency of the cycle are highest with a NO2/NOx ratio of 0.5. The outlet NH3 is lowest for the NO2/NOx ratio greater than 0.6. A combined engine experimental and simulation study was performed to quantify the NH3 maldistribution at the SCR inlet and its effects on the SCR performance and kinetics. The uniformity index (UI) of the SCR inlet NH3 and NH3/NOx ratio (ANR) was determined to be below 0.8 for the production system. The UI was improved to 0.9 after installation of a swirl mixer into the SCR inlet cone. A multi-channel model was developed to simulate the maldistribution effects. The results showed that reducing the UI of the inlet ANR from 1.0 to 0.7 caused a 5-10% decrease in NOx reduction efficiency and 10-20 ppm increase in the NH3 slip. The simulations of the steady-state engine data with the multi-channel model showed that the NH3 maldistribution is a factor causing the differences in the calibrations developed from the engine and the reactor data. The Reactor experiments were performed at ORNL using a Spaci-IR technique to study the thermal aging effects. The test results showed that the thermal aging (at 800°C for 16 hours) caused a 30% reduction in the NH3 stored on the catalyst under NH3 saturation conditions and different axial concentration profiles under SCR reaction conditions. The kinetics analysis showed that the thermal aging caused a reduction in total NH3 storage capacity (94.6 compared to 138 gmol/m3), different NH3 adsorption/desorption properties and a decrease in activation energy and the pre-exponential factor for NH3 oxidation, standard and fast SCR reactions. Both reduction in the storage capability and the change in kinetics of the major reactions contributed to the change in the axial storage and concentration profiles observed from the experiments.