Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

Insights on keto-hydroperoxide formation from O 2 addition to the beta-tetrahydrofuran radical

Cyclic ethers provide an interesting case study of low-temperature oxidation chemistry, especially in relevance to biofuels. A recent experimental study revealed new questions regarding the low-temperature oxidation mechanism of tetrahydrofuran (THF) concerning the formation of keto-hydroperoxides. In particular, keto-hydroperoxides originating from the THF-β radical were not captured accurately by current literature models, motivating this work to calculate the energetics of the first and second O 2 -addition pathways for THF radicals. Electronic structure calculations at the CCSD(T)/cc-pV∞Z//M06-2X/cc-pVTZ level of theory were used to generate potential energy surfaces for the α-C and β-C THF radicals and subsequent pathways to the formation of the keto-hydroperoxide isomers. These are the first theoretical calculations of the second O 2 -addition radical pathways for the THF-β radical. Results from the theorical work provided further insight into the low-temperature oxidation of THF. This included identifying the pathways most likely to form the keto-hydroperoxide isomers observed in prior experimental work; and detecting that the shortcomings in prior models are likely due to uncertainties in R + THF abstraction reaction rates. Here, these conclusions will motivate future work for accurate THF kinetic model development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Performance characterization of nanofiltration, reverse osmosis, and ion exchange technologies for acetic acid separation

A major obstacle to widespread implementation of bio-based fuels and chemicals is the lack of efficient and costeffective separation methods. To purify acetic acid produced by biochemical conversion of biomass via anaerobic digestion, this work employs two commonly used separation technologies: (1) ion-exchange (IX) resin and (2) pressure-driven membranes. This study tested five commercially available strong- and weak-base anion-exchange resins and five commercially available nanofiltration (NF) and reverse osmosis (RO) membranes. The pH of the feed solution significantly affected adsorption capacity. At pH 6.3, a strong-base IX resin (IRN-78) performed best (95.1% acetate removal). With strong-base IX resins, the Langmuir isotherm model fit well, whereas for weak-base IX resins, the Freundlich isotherm provided a better fit. A pseudo-second-order kinetic model fit well for both IRN-78 and IRA-67. Regarding membrane separation, RO (BW30XFR membrane) achieved the highest rejection (98.6% acetate rejection), whereas an NF membrane (NF*) achieved the best combination of permeate flux (105 L/(h∙m 2 )) and rejection (83.1% acetate rejection). For membrane performance, the experimental data were fit using the solution diffusion model. Increased pH in the feed solution lowered permeate flux but increased acetic acid rejection. Finally, then the acetic acid concentration in the feed solution increased, both permeate flux and acetic acid rejection decreased for membrane NF*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Environmentally relevant organophosphate triesters in herring gulls: In vitro biotransformation and kinetics and diester metabolite formation using a hepatic microsomal assay

The in vitro biotransformation and kinetics of six organophosphate triester (OPE) flame retardants were investigated in herring gulls (Larus argentatus) from the Great Lakes using a hepatic microsomal metabolism assay. Administration of each individual OPE (tri-n-butyl phosphate (TNBP), tris(2-butoxyethyl) phosphate (TBOEP), triphenyl phosphate (TPHP), triethyl phosphate (TEP), tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) and tris(2-chloroisopropyl) phosphate (TCIPP)) to the in vitro assay (concentration range 0.01 to 10 μM) resulted in rapid depletion with the exception of TEP. Following the Michaelis-Menten enzyme kinetics model, a preliminary 2-minute incubation period was used to estimate the V{sub max} (± SE) values (i.e., the maximal rate of reaction for a saturated enzyme system), which ranged from 5.0 ± 0.4 (TPHP) to 29 ± 18 pmol/min/mg protein (TBOEP), as well as the K{sub M} (± SE) values (i.e., the OPE concentration corresponding to one half of the V{sub max}), which ranged from 9.8 ± 1 (TPHP) to 189 ± 135 nM (TBOEP). Biotransformation assays over a 100-minute incubation period revealed that TNBP was metabolized most rapidly (with a depletion rate of 73 ± 4 pmol/min/mg protein), followed by TBOEP (53 ± 8 pmol/min/mg), TCIPP (27 ± 1 pmol/min/mg), TPHP (22 ± 2 pmol/min/mg) and TDCIPP (8 ± 1 pmol/min/mg). In vitro biotransformation of OP triesters was clearly structure-dependent where non-halogenated alkyl OP triesters were metabolized more rapidly than halogenated alkyl triesters. Halogenated OP triesters were transformed to their respective diesters more efficiently relative to non-halogenated OP triesters. To our knowledge, this is the first study to investigate OP triester metabolism and OP diester formation in an avian or wildlife model system, which is important to understand the fate and biological activity of OPEs in an exposed organism. - Highlights: • The metabolism and kinetics of 6 OPEs were examined in herring gull liver microsomes. • The metabolism of OPEs was rapid, with the exception of TEP (no metabolism observed). • Relative metabolism rates (pmol/min/mg protein): TNBP > TBOEP > TPHP > TDCIPP > TCIPP. • OPE metabolite (i.e., OP diesters) formation was shown for all 5 OP triesters. • Structure-dependent metabolism rates and OP diester formation were observed.

54 ENVIRONMENTAL SCIENCES↗

Catalysis Science Initiative: Catalyst Design by Discovery Informatics

Catalysts selectively enhance the rates of chemical reactions toward desired products. Such reactions provide great benefit to society in major commercial sectors such as energy production, protecting the environment, and polymer products and thereby contribute heavily to the country’s gross national product. Our premise is that the level of fundamental understanding of catalytic events at the atomic and molecular scale has reached the point that more predictive methods can be developed to shorten the cycle time to new processes. The field of catalysis can be divided into two regimes: heterogeneous and homogeneous. For the heterogeneous catalysis regime, we have used the water-gas shift (WGS) reaction (CO + H2O + CO2 + H2O) over supported metals as a test bed. Detailed analysis and strong coupling of theory with experiment have led to the following conclusions: • The sequence of elementary steps goes through a COOH intermediate • The CO binding energy is a strong function of coverage of CO adsorbed on the surface in many systems • In the case of Au catalysts, the CO adsorption is generally too weak on surface with close atomic packing, but the enhanced binding at corner atoms (which are missing bonding partners) of cubo-octahedral nanoparticles increases the energy to a near optimal value and produces very active catalysts. • Reaction on the metal alone cannot account for the experimental results. The reaction is dual functional with water activation occurring at the metal-support interface. It is clear from our work that the theory component is essential, not only for prediction of new systems, but also for reconciling data and testing hypotheses regarding potential descriptors. Particularly important is the finding that the interface between nano-sized metal particles and the oxides that are used to support them represent a new state of matter in the sense that the interfacial bonding perturbs the chemical state of both metals atoms and the support atoms in the interfacial region. Some of the first theoretical descriptions of this important chemistry and potential new source of control of catalyst properties are be in preparation for submission. On the homogeneous catalysis side, we have used single site olefin polymerization as the testbed. This system is important because changes in a single ligand bonded to the catalytically active metal site can alter the rates of individual steps in the polymerization sequence and thereby change the properties of the resulting polymer, potentially improving its value in a hundred million pound per year industry. We have made a major advance in understanding such systems by developing a population balance kinetic model that allows us to predict the molecular weight distribution (MWD) of the product. That, in turn, allows use of MWD data to fit kinetic parameters. By combining monomer loss data, MWD, measurement of the number of working active sites, and polymer end group analysis, we have a rich data set that is highly discriminating of kinetic mechanism. Thus, we have a robust tool for producing high quality, detailed kinetic parameters, which we have used to refine mechanisms presented in the literature and discover relationships between steric and electronic properties of group IV catalysts and individual rate constants in a number of systems. Our recent work on six-coordinate Zr, Ti, and Hf amine bis(phenolate) systems, we have shown that: • The sterics (bulkiness) of the ligands specifically affect the chain termination reaction • The electron density on the metal controls misinsertion (flipped orientation) of the olefin into the growing polymer • Steric effects related to the size of the ortho ligand on the catalyst have been shown to strongly affect its the degree of dormancy, i.e. tendency to stop reacting • Changes in the size of the amine pendent group on the catalyst can have such a strong effect on chain termination as to change the catalyst from one that produces only oligomers to one that incorporates oligomers in polymer chains to introduce chain branching. These effects control the length of molecular chains, the number of errors in the regularity of the chains, and the variation of chain lengths in the final product distribution. These characteristics, in turn, determine the properties of the resulting polymer material. Predictive modeling of polymerization process cannot be done without the quantitative rate constant determination that our unique and comprehensive approach has produced. This tool, combined with insights on how catalyst chemistry affects the rate constants, is an important step toward establishing the quantitative relationships between catalyst chemistry and polymer properties that will allow more efficient searches for new catalysts and speed the discovery process. Thus, in both the heterogeneous and homogeneous catalysis areas, this work has brought important new understanding to the field. The technical effectiveness of our approach is demonstrated by its success. For heterogeneous catalysis, experimental work alone showed the important of the metal-support interface, but the integration with theory has allowed deeper understanding of the fundamental source of the effects and shown the path for how to use this understanding for discovery. In the case of homogeneous catalysis, methodology for determination of which reactions are needed to account for the data and the extraction of quantitative values for rate constants for those reactions allows accurate modeling of reaction behavior. This allowed correlation of trends in catalyst properties with alterations of specific rates, thus contributing to a data base that will guide catalyst discovery. While not yet realized, the economic benefits of this approach will come in significantly shortened cycle times for discovery and ability to fill product demand at lower price.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

New insights into the low-temperature oxidation of 2-methylhexane

The low-temperature oxidation of a stoichiometric 2-methylhexane/O 2 /Ar mixture was examined in a jet-stirred reactor coupled with synchrotron vacuum UV photoionization molecular-beam MS. The initial gas mixture was composed of 2% 2-methyhexane, 22% O 2 and 76% Ar and the pressure of the reactor was kept at 780 Torr. Low-temperature oxidation intermediates with two to five oxygen atoms were observed. The detection of C 7 H 14 O 5 and C 7 H 12 O 4 species suggests that a third O 2 addition process occurs in 2-methylhexane low-temperature oxidation. A detailed kinetic model was developed that describes the third O 2 addition and subsequent reactions leading to C 7 H 14 O 5 (keto-dihydroperoxide and dihydroperoxy cyclic ether) and C 7 H 12 O 4 (diketo-hydroperoxide and keto-hydroperoxy cyclic ether) species. The kinetics of the third O 2 addition reactions are discussed and model calculations were performed that reveal that third O 2 addition reactions promote 2-methylhexane auto-ignition at low temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Refined Outgassing Model for the LiH/LiOH System

Previous models for H 2 and H 2 O outgassing from moisture exposed LiH specimens were based on temperature programmed desorption/decomposition (TPD) experiments and analyses that yielded kinetic descriptions which were specific to defined sets of initial conditions. Subsequently the inputs to the kinetic models were refined to correlate with infrared spectra of the adsorbed species. In this current enhancement, better agreement between prediction and experiment was obtained through: 1) refining the reaction mechanisms for some of the decomposition reactions; 2) developing the potential range of kinetic parameters by simultaneous fitting of all the TPD outgassing curves obtained at different heating rates; and, 3) limiting the scope of possible kinetic parameters to an optimized subset by matching the outgassing prediction curves at 3 different temperatures to the outcomes from the accurate but cumbersome, non-intuitive, and effort demanding isoconversional analysis.

36 MATERIALS SCIENCE↗

Thermal-Hydraulic Simulations of Single Pin and Assembly Sector for IVG- 1M Reactor

Thermal-hydraulic simulations have been performed using computational fluid dynamics (CFD) for the highly-enriched uranium (HEU) design of the IVG.1M reactor at the Institute of Atomic Energy (IAE) at the National Nuclear Center (NNC) in the Republic of Kazakhstan. Steady-state simulations were performed for both types of fuel assembly (FA), i.e. the FA in rows 1 & 2 and the FA in row 3, as well as for single pins in those FA (600 mm and 800 mm pins). Both single pin calculations and bundle sectors have been simulated for the most conservative operating conditions corresponding to the 10 MW output power, which corresponds to a pin unit cell Reynolds number of only about 7500. Simulations were performed using the commercial code STAR-CCM+ for the actual twisted pin geometry as well as a straight-pin approximation. Various Reynolds-Averaged Navier-Stokes (RANS) turbulence models gave different results, and so some validation runs with a higher-fidelity Large Eddy Simulation (LES) code were performed given the lack of experimental data. These singled out the Realizable Two-Layer k-ε as the most accurate turbulence model for estimating surface temperature. Single-pin results for the twisted case, based on the average flow rate per pin and peak pin power, were conservative for peak clad surface temperature compared to the bundle results. Also the straight-pin calculations were conservative as compared to the twisted pin simulations, as expected, but the single-pin straight case was not always conservative with regard to the straight-pin bundle. This was due to the straight-pin temperature distribution being strongly influenced by the pin orientation, particularly near the outer boundary. The straight-pin case also predicted the peak temperature to be in a different location than the twisted-pin case. This is a limitation of the straight-pin approach. The peak temperature pin was in a different location from the peak power pin in every case simulated, and occurred at an inner pin just before the enrichment change. The 600 mm case demonstrated a peak clad surface temperature of 370.4 K, while the 800 mm case had a temperature of 391.6 K. These temperatures are well below the necessary temperatures for boiling to occur at the rated pressure. Fuel temperatures are also well below the melting point. Future bundle work will include simulations of the proposed low-enriched uranium (LEU) design. Two transient scenarios were also investigated for the single-pin geometries. Both were “model” problems that were focused on pure thermal-hydraulic behavior, and as such were simple power changes that did not incorporate neutron kinetics modeling. The first scenario was a high-power, ramp increase, while the second scenario was a low-power, step increase. A cylindrical RELAP model was also constructed to investigate its accuracy as compared to the higher-fidelity CFD. Comparisons between the two codes showed good agreement for peak temperatures in the fuel and at the cladding surface for both cases. In the step transient, temperatures at four axial levels were also computed. These showed greater but reasonable discrepancy, with RELAP outputting higher temperatures. These results provide some evidence that RELAP can be used with confidence in modeling transients for IVG.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

In-Cylinder Diagnostics to Overcome Efficiency Barriers in Natural Gas Engines (LDRD 200166)

The high-level objective of this project is to solve national-security problems associated with petroleum use, cost, and environmental impacts by enabling more efficient use of natural-gas-fueled internal combustion engines. An improved science-base on end-gas autoignition, or “knock,” is required to support engineering of more efficient engine designs through predictive modeling. An existing optical diesel engine facility is retrofitted for natural gas fueling with laser-spark-ignition combustion to provide in-cylinder imaging and pressure data under knocking combustion. Zero-dimensional chemical-kinetic modeling of autoignition, adiabatically constrained by the measured cylinder pressure, isolates the role of autoignition chemistry. OH* chemiluminescence imaging reveals six different categories of knock onset that depend on proximity to engine surfaces and the in-cylinder deflagration. Modeling results show excellent prediction regardless of the knock category, thereby validating state-of-the-art kinetic mechanisms. The results also provide guidance for future work to build a science base on the factors that affect the deflagration rate.

33 ADVANCED PROPULSION SYSTEMS↗

Reaction kinetic analysis of the 3-hydroxypropionate/4-hydroxybutyrate CO 2 fixation cycle in extremely thermoacidophilic archaea

Here, the 3-hydroxypropionate/4-hydroxybutyrate (3HP/4HB) cycle fixes CO 2 in extremely thermoacidophilic archaea and holds promise for metabolic engineering because of its thermostability and potentially rapid pathway kinetics. A reaction kinetics model was developed to examine the biological and biotechnological attributes of the 3HP/4HB cycle as it operates in Metallosphaera sedula, based on previous information as well as on kinetic parameters determined here for recombinant versions of five of the cycle enzymes (malonyl-CoA/succinyl-CoA reductase, 3-hydroxypropionyl-CoA synthetase, 3-hydroxypropionyl-CoA dehydratase, acryloyl-CoA reductase, and succinic semialdehyde reductase). The model correctly predicted previously observed features of the cycle: the 35%–65% split of carbon flux through the acetyl-CoA and succinate branches, the high abundance and relative ratio of acetyl-CoA/propionyl-CoA carboxylase (ACC) and MCR, and the significance of ACC and hydroxybutyryl-CoA synthetase (HBCS) as regulated control points for the cycle. The model was then used to assess metabolic engineering strategies for incorporating CO 2 into chemical intermediates and products of biotechnological importance: acetyl-CoA, succinate, and 3-hydroxyproprionate.

3-hydroxypropionate↗

Elucidating the photodissociation fingerprint and quantifying the determination of organic hydroperoxides in gas-phase autoxidation

Hydroperoxides are formed in the atmospheric oxidation of volatile organic compounds, in the combustion autoxidation of fuel, in the cold environment of the interstellar medium, and also in some catalytic reactions. They play crucial roles in the formation and aging of secondary organic aerosols and in fuel autoignition. However, the concentration of organic hydroperoxides is seldom measured, and typical estimates have large uncertainties. In this work, we developed a mild and environmental-friendly method for the synthesis of alkyl hydroperoxides (ROOH) with various structures, and we systematically measured the absolute photoionization cross-sections (PICSs) of the ROOHs using synchrotron vacuum ultraviolet-photoionization mass spectrometry (SVUV-PIMS). A chemical titration method was combined with an SVUV-PIMS measurement to obtain the PICS of 4-hydroperoxy-2-pentanone, a typical molecule for combustion and atmospheric autoxidation ketohydroperoxides (KHPs). We found that organic hydroperoxide cations are largely dissociated by loss of OOH. This fingerprint was used for the identification and accurate quantification of the organic peroxides, and it can therefore be used to improve models for autoxidation chemistry. The synthesis method and photoionization dataset for organic hydroperoxides are useful for studying the chemistry of hydroperoxides and the reaction kinetics of the hydroperoxy radicals and for developing and evaluating kinetic models for the atmospheric autoxidation and combustion autoxidation of the organic compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗

Implementation of new mixture rules has a substantial impact on combustion predictions for H 2 and NH 3

Complex-forming reactions comprise a substantial fraction of all important combustion reactions and are central to combustion behavior. Despite being often called “pressure-dependent” reactions, their rate constants depend on not only the pressure but also the composition. While modern combustion codes allow arbitrarily high accuracy in treating pressure dependence, recent work has consistently demonstrated dramatic failures of essentially all available treatments of mixture dependence. In situations where mixture dependence is treated at all, it is inevitably treated through specification of pressure-dependent rate constants for a set of pure bath gases, which are then combined to estimate the rate constant in a mixture via a “mixture rule.” While there had been a generally unquestioning confidence in these mixture rules, they had, in reality, been scarcely tested until the last decade, when comparisons against master equation calculations revealed order-of-magnitude errors for important pressure-dependent reactions. New mixture rules, based on the reduced pressure, have recently been proposed and shown to reproduce master equation calculations for broad classes of complex-forming reactions very accurately. Here, in this work, we present an implementation of one such new mixture rule (“LMR-R”) in Cantera and then use it to enable simulations that use new high-accuracy ab initio data for individual bath gases (for the first time, since codes previously could not accommodate the complex bath gas dependence). Demonstrations focus on combustion of H 2 and NH 3 , where (1) high-accuracy ab initio data are available and (2) the impact is expected to be large due to the high fractions of efficient colliders (e.g., H 2 O and NH 3 ) in the burned and unburned gases. Indeed, we find the impact of this treatment to be substantial and may explain previous modeling difficulties for these important carbon-free fuels, particularly for NH 3 , whose extraordinarily high third-body efficiency (~20) is often omitted from kinetic models.

Ammonia↗

Morphological Evolution of Nanocluster Aggregates and Single Crystals in Alkaline Zinc Electrodeposition

The morphology of Zn electrodeposits is studied on carbon-coated transmission electron microscopy grids. At low over-potentials (eta = -50 mV), the morphology develops by aggregation at two distinct length scales: similar to 5 nm diameter monocrystalline nanoclusters form similar to 50 nm diameter polycrystalline aggregates, and the aggregates form a branched network. Epitaxial (00 (0) over bar2) growth above an overpotential of vertical bar eta(c)vertical bar > 125 mV leads to the formation of hexagonal single crystals up to 2 mu m in diameter. Potentiostatic current transients were used to calculate the nucleation rate from Scharifker et al.'s model. The exp(eta) dependence of the nucleation rates indicates that atomistic nucleation theory explains the nucleation process better than Volmer-Weber theory. A kinetic model is provided using the rate equations of vapor solidification to simulate the evolution of the different morphologies. On solving these equations, we show that aggregation is attributed to cluster impingement and cluster diffusion while single-crystal formation is attributed to direct attachment.

Desai, D↗

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗

Patient-specific modeling of individual sickle cell behavior under transient hypoxia

Sickle cell disease (SCD) is a highly complex genetic blood disorder in which red blood cells (RBC) exhibit heterogeneous morphology changes and decreased deformability. We employ a kinetic model for cell morphological sickling that invokes parameters derived from patient-specific data. This model is used to investigate the dynamics of individual sickle cells in a capillary-like microenvironment in order to address various mechanisms associated with SCD. We show that all RBCs, both hypoxia-unaffected and hypoxia-affected ones, regularly pass through microgates under oxygenated state. However, the hypoxia-affected cells undergo sickling which significantly alters cell dynamics. In particular, the dense and rigid sickle RBCs are obstructed thereby clogging blood flow while the less dense and deformable ones are capable of circumnavigating dead (trapped) cells ahead of them by choosing a serpentine path. Informed by recent experiments involving microfluidics that provide in vitro quantitative information on cell dynamics under transient hypoxia conditions, we have performed detailed computational simulations of alterations to cell behavior in response to morphological changes and membrane stiffening. Our model reveals that SCD exhibits substantial heterogeneity even within a particular density-fractionated subpopulation. These findings provide unique insights into how individual sickle cells move through capillaries under transient hypoxic conditions, and offer novel possibilities for designing effective therapeutic interventions for SCD.

59 BASIC BIOLOGICAL SCIENCES↗

Concepts Relevant for the Kinetic Analysis of Reversible Reaction Systems

The net rate of a reversible chemical reaction is the difference between unidirectional rates of traversal along forward and reverse reaction paths. In a multistep reaction sequence, the forward and reverse trajectories, in general, are not the microscopic reverse of one another; rather, each unidirectional route is comprised of distinct rate-controlling steps, intermediates, and transition states. Consequently, traditional descriptors of rate (e.g., reaction orders) do not reflect intrinsic kinetic information but instead conflate unidirectional contributions determined by (i) the microscopic occurrence of forward/reverse reactions (i.e., unidirectional kinetics) and (ii) the reversibility of reaction (i.e., nonequilibrium thermodynamics). This review aims to provide a comprehensive resource of analytical and conceptual tools which deconvolute the contributions of reaction kinetics and thermodynamics to disambiguate unidirectional reaction trajectories and precisely identify rate- and reversibility-controlling molecular species and steps in reversible reaction systems. The extrication of mechanistic and kinetic information from bidirectional reactions is accomplished through equation-based formalisms (e.g., De Donder relations) grounded in principles of thermodynamics and interpreted in the context of theories of chemical kinetics developed in the past 25 years. Here, the aggregate of mathematical formalisms detailed herein is general to thermochemical and electrochemical reactions and encapsulates a diverse body of scientific literature encompassing chemical physics, thermodynamics, chemical kinetics, catalysis, and kinetic modeling.

42 ENGINEERING↗

On the Statistical Calibration of Physical Models

ABSTRACT We introduce a novel statistical calibration framework for physical models, relying on probabilistic embedding of model discrepancy error within the model. For clarity of illustration, we take the measurement errors out of consideration, calibrating a chemical model of interest with respect to a more detailed model, considered as “truth” for the present purpose. We employ Bayesian statistical methods for such model‐to‐model calibration and demonstrate their capabilities on simple synthetic models, leading to a well‐defined parameter estimation problem that employs approximate Bayesian computation. The method is then demonstrated on two case studies for calibration of kinetic rate parameters for methane air chemistry, where ignition time information from a detailed elementary‐step kinetic model is used to estimate rate coefficients of a simple chemical mechanism. We show that the calibrated model predictions fit the data and that uncertainty in these predictions is consistent in a mean‐square sense with the discrepancy from the detailed model data.

Sargsyan, K.↗