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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multiscale Modeling of a Direct Nonoxidative Methane Dehydroaromatization Reactor with a Validated Model for Catalyst Deactivation

Due to the recent boom in shale gas production, aromatics production using direct nonoxidative methane dehydroaromatization (DHA) is being investigated extensively. However, due to rapid coke formation, catalysts in the nonoxidative methane DHA reactors get deactivated, which is one of the critical issues for the commercial success of the methane DHA process. In this paper, a model for catalyst deactivation is developed. Rate models for other DHA reactions are developed by considering the decrease in the catalyst activity with time. Due to the very fast coke formation rate on the fresh catalyst, there is coke formation immediately upon the introduction of the feed. Therefore, an algorithm is developed for estimation of the initial state of the reactor and the kinetic parameters by coupling an iterative direct substitution approach with an optimization approach. Transient experimental data from an in-house reactor are first reconciled and then used for developing the kinetic model including the coke formation model. Using the rate model, a dynamic, heterogeneous, multiscale reactor model with embedded heating is developed. Here, the model couples the catalyst pellet level model with a reactor level model. Impacts of temperature, L/D ratio, and scheduling of reactors on variability in conversion and yield with time are studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Origin of Multiple Peaks in the Potentiodynamic Oxidation of CO Adlayers on Pt and Ru-Modified Pt Electrodes

The research of the electrooxidation mechanism of CO ad on Pt based catalysts is very important for designing more effective CO-tolerant electrocatalysts for fuel cells. We have studied the origin of multiple peaks in the cyclic voltammograms of CO stripping from polycrystalline Pt and Ru modified polycrystalline Pt (Pt/Ru) surfaces in both acidic and alkaline media by differential electrochemical mass spectrometry (DEMS), DFT calculations, and kinetic Monte Carlo (KMC) simulations. A novel CO ad electrooxidation kinetic model on heterogeneous Pt and Pt/Ru catalysts is proposed to account for the multiple peaks experimentally observed. In this model, OH species prefer to adsorb at low-coordination sites or Ru sites and, thus, suppress CO repopulation from high-coordination sites onto these sites. Therefore, CO ad oxidation occurs on different facets or regions, leading to multiplicity of CO stripping peaks. This work provides a new insight into the CO electrooxidation mechanism and kinetics on heterogeneous catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

NMR spin-rotation relaxation and diffusion of methane

The translational diffusion-coefficient and the spin-rotation contribution to the 1 H NMR relaxation rate for methane (CH 4 ) are investigated using MD (molecular dynamics) simulations, over a wide range of densities and temperatures, spanning the liquid, supercritical, and gas phases. Here, the simulated diffusion-coefficients agree well with measurements, without any adjustable parameters in the interpretation of the simulations. A minimization technique is developed to compute the angular velocity for non-rigid spherical molecules, which is used to simulate the autocorrelation function for spin-rotation interactions. With increasing diffusivity, the autocorrelation function shows increasing deviations from the single-exponential decay predicted by the Langevin theory for rigid spheres, and the deviations are quantified using inverse Laplace transforms. The 1 H spin-rotation relaxation rate derived from the autocorrelation function using the “kinetic model” agrees well with measurements in the supercritical/gas phase, while the relaxation rate derived using the “diffusion model” agrees well with measurements in the liquid phase. 1 H spin-rotation relaxation is shown to dominate over the MD-simulated 1 H- 1 H dipole-dipole relaxation at high diffusivity, while the opposite is found at low diffusivity. At high diffusivity, the simulated spin-rotation correlation time agrees with the kinetic collision time for gases, which is used to derive a new expression for 1 H spin-rotation relaxation, without any adjustable parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Modeling Functional Organic Chemistry in Arctic Rivers: An Idealized Siberian System

Rivers of the Arctic will become ever more important for the global climate, since they carry a majority of continental dissolved organic carbon flux into the rapidly changing polar ocean. Aqueous organics comprise a wide array of functional groups, several of which are likely to impact coastal and open water biophysical properties. Light attenuation, interfacial films, aerosol formation, gas release and momentum exchange can all be cited. We performed Lagrangian kinetic modeling for the evolution of riverine organic chemistry as the molecules in question make their way from the highlands to Arctic outlets. Classes as diverse as the proteins, sugars, lipids, re-condensates, humics, bio-tracers and small volatiles are all included. Our reduced framework constitutes an idealized northward flow driving a major hydrological discharge rate and primarily representing the Russian Lena. Mountainous, high solute and tundra sources are all simulated, and they meet up at several points between soil and delta process reactors. Turnover rates are parameterized beginning with extrapolated coastal values imposed along a limited tributary network, with connections between different terrestrial sub-ecologies. Temporal variation of our total dissolved matter most closely resembles the observations when we focus on the restricted removal and low initial carbon loads, suggesting relatively slow transformation along the water course. Thus, channel combinations and mixing must play a dominant role. Nevertheless, microbial and photochemical losses help determine the final concentrations for most species. Chemical evolution is distinct for the various functionalities, with special contributions from pre- and post-reactivity in soil and delta waters. Several functions are combined linearly to represent the collective chromophoric dissolved matter, characterized here by its absorption. Tributaries carry the signature of lignin phenols to segregate tundra versus taiga sources, and special attention is paid to the early then marine behaviors of low molecular weight volatiles. Heteropolycondensates comprise the largest percentage of reactive carbon in our simulations due to recombination/accumulation, and they tend to be preeminent at the mouth. Outlet concentrations of individual structures such as amino acids and absorbers lie above threshold values for biophysical influence, on the monolayer and light attenuation. The extent of coastal spreading is examined through targeted regional box modeling, relying on salinity and color for calibration. In some cases, plumes reach the scale of peripheral arctic seas, and amplification is expected during upcoming decades. Conclusions are mapped from the Lena to other boreal discharges, and future research questions are outlined regarding the bonding type versus mass release as permafrost degrades. Dynamic aqueous organic coupling is recommended for polar system models, from headwaters to coastal diluent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Measurement of Kinetic Parameters by Oscillated Experiments

The availability of reliable values of the kinetic parameters (effective delayed neutron fraction {beta} and effective mean neutron generation time {Lambda}) is important for reactor kinetic analyses and for reactivity monitoring when adopting the point kinetic approach. The kinetic parameters may be obtained either by simulations (Monte Carlo or deterministic) of the system or by direct measurements. Several methods have been proposed for the direct experimental determination of the kinetic parameters, e.g. control rod drops and noise analysis approaches, e.g. In this work, the problem of the determination of kinetic parameters in subcritical systems is tackled. The system is subjected to an oscillated source and the value of the ratio {beta}/{Lambda} is related to a measurable quantity when the system response reaches an equilibrium condition. In the analysis presented, a pure frequency is considered: when analyzing more complicated source periodic behaviours, a Fourier analysis can allow to isolate the contributions of individual frequencies. The phase difference between the oscillating source and the system response is adopted as the experimental observable. The method for the interpretation of the measurements is based on the standard point kinetic model. The response of the system is obtained by localized neutron detectors and, therefore, spatial and spectral effects may also affect the results. This drawback needs some further discussion and will be considered in the future work. The phase shift between the source signal and the system response may be analytically retrieved by Laplace transforming the point kinetic equations in the presence of an oscillated source. An explicit formula relates the phase shift to the ratio {beta}/{Lambda} and to the frequency of the source. Based on this relationship, the ratio {beta}/{Lambda} can be determined. To validate the technique, the experimental data are simulated by a numerical code solving the time-dependent diffusion equations for a source-driven system. The work allows to assess the possibility of the technique and its limits of applicability. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Protein folding from heterogeneous unfolded state revealed by time-resolved X-ray solution scattering

One of the most challenging tasks in biological science is to understand how a protein folds. In theoretical studies, the hypothesis adopting a funnel-like free-energy landscape has been recognized as a prominent scheme for explaining protein folding in views of both internal energy and conformational heterogeneity of a protein. Despite numerous experimental efforts, however, comprehensively studying protein folding with respect to its global conformational changes in conjunction with the heterogeneity has been elusive. Here we investigate the redox-coupled folding dynamics of equine heart cytochrome c (cyt-c) induced by external electron injection by using time-resolved X-ray solution scattering. A systematic kinetic analysis unveils a kinetic model for its folding with a stretched exponential behavior during the transition toward the folded state. With the aid of the ensemble optimization method combined with molecular dynamics simulations, we found that during the folding the heterogeneously populated ensemble of the unfolded state is converted to a narrowly populated ensemble of folded conformations. These observations obtained from the kinetic and the structural analyses of X-ray scattering data reveal that the folding dynamics of cyt-c accompanies many parallel pathways associated with the heterogeneously populated ensemble of unfolded conformations, resulting in the stretched exponential kinetics at room temperature. This finding provides direct evidence with a view to microscopic protein conformations that the cyt-c folding initiates from a highly heterogeneous unfolded state, passes through still diverse intermediate structures, and reaches structural homogeneity by arriving at the folded state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Understanding fundamental effects of biofuel structure on ignition and physical fuel properties

The development of biochemicals and biofuels with advantaged properties of higher combustion efficiency at lower emission levels can be aided by a priori prediction of critical global characteristics based on molecular structure. Existing predictive techniques are incomplete for comprehensive ignition and physical fuel property prediction. Here this investigation focuses on filling this knowledge gap by developing a priori prediction techniques for predicting important ignition as well as physical fuel properties of promising advantaged biofuels. We demonstrate the utility of this a priori prediction technique on saturated (≤C 5 ) alcohols, both linear and branched in nature by correlating their molecular structure with both important ignition and physical properties of the fuels. This work utilizes a semi-automated method of predicting ignition characteristics via accurate fundamental-based calculations of the fuel radical cascades that govern low-temperature combustion. This numerical methodology is supported by the relevant fuel ignition delay time (IDT) and research octane number (RON) data experimentally measured by the Advanced Fuel Ignition Delay Analyzer (AFIDA). Phyiscal fuel property trends are also analyzed by correlating the molecular structure of the fuel to its physical characteristics. Based on their molecular structure, recommendations are made on an a priori method for selecting biofuels that have advantaged properties for the selected fuel application. Lastly, the utility of this investigation is extended beyond the subset of alcohols considered by analyzing the IDT trends of i-Pentanol and n-Hexanol, revealing the fidelity of the novel kinetic modelling approach and the accuracy of the predicted trends in this study.

09 BIOMASS FUELS↗
DOE OSTI2023

Establishing a Design Strategy for Corrosion Resistant Structural Materials in Molten Salt Technologies

To address the gap in material design strategies for structural components in molten salt-powered technologies, the current work aimed to identify key parameters governing corrosion induced degradation of materials during exposures in molten halide salts. The applicability of the data on dissolution rates of pure metals (e.g., Cr) in a given molten salt to enable predictions of corrosion behavior of multicomponent binary and ternary Fe-based model and Ni-based commercial alloys isothermally exposed in the purified KCl-MgCl 2 in quartz capsules at 700 and 800°C was evaluated. The influence of alloy chemistry on initial dissolution rates and the consequent impact on the time-dependent corrosion induced microstructural evolution in the alloy subsurface was predicted with a coupled thermodynamic-kinetic model. Finally, the model was well able to predict the Cr depletion in all investigated alloys demonstrating the ability of the proposed approach to enable a unique design strategy for structural materials in molten salts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic and cluster dynamics modeling of fission gas (Xe) diffusivity in TRISO fuel kernels

TRISO fuel particles are candidates for use in next generation reactors including gas reactors, fluoride salt-cooled high temperature reactors, and micro-reactors. The UCO fuel kernel consists of a uranium dioxide (UO) and uranium carbide mixture. The addition of UC helps suppress the formation of carbon monoxide gas, which led to failures during initial TRISO development. The addition of uranium carbide alters the chemistry of the UO kernel, which is known to influence performance parameters such as fission gas diffusivity, although the impact has not been quantified and no models exist that take the change in chemistry into account. Therefore, better understanding and more accurate models of the impact of chemistry on fuel performance are of high priority. In this paper, a first-principles density functional theory (DFT) and empirical potential based multi-scale study has been carried out to model the diffusivity of fission gas xenon (Xe) in UCO TRISO fuel kernels. The focus is on the UO component in the UCO fuel kernels, as that represents the largest volume fraction of the fuel kernels. The study relies on DFT and empirical potential calculations to determine Xe and point defect properties, which are then used in thermodynamic and kinetic models to predict diffusion for intrinsic conditions. In addition, the information is utilized in cluster dynamics simulations using the Centipede code to estimate the impact of irradiation on defect transport. Additionally, the presence of UC or UC in the UCO fuel kernels is shown to have a substantial impact on the UO non-stoichiometry by inducing oxygen vacancies and driving UO sub-stoichiometric, which causes much slower Xe diffusion in UCO compared to light water reactor UO fuel. The application of this model in fuel performance simulations using the Bison code is also demonstrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Second order kinetic theory of parallel momentum transport in collisionless drift wave turbulence

A second order kinetic model for turbulent ion parallel momentum transport is presented. A new nonresonant second order parallel momentum flux term is calculated. The resonant component of the ion parallel electrostatic force is the momentum source, while the nonresonant component of the ion parallel electrostatic force compensates for that of the nonresonant second order parallel momentum flux. The resonant component of the kinetic momentum flux can be divided into three parts, including the pinch term, the diffusive term, and the residual stress. By reassembling the pinch term and the residual stress, the residual stress can be considered as a pinch term of parallel wave-particle resonant velocity, and, therefore, may be called as “resonant velocity pinch” term. Considering the resonant component of the ion parallel electrostatic force is the transfer rate between resonant ions and waves (or, equivalently, nonresonant ions), a conservation equation of the parallel momentum of resonant ions and waves is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Mass Spectrometer Transient Analysis

This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.

Kristy, Stephen [Idaho National Laboratory (INL), ↗
DOE OSTI2015

HO x radical chemistry in oxidation flow reactors with low-pressure mercury lamps systematically examined by modeling

Oxidation flow reactors (OFRs) using OH produced from low-pressure Hg lamps at 254 nm (OFR254) or both 185 and 254 nm (OFR185) are commonly used in atmospheric chemistry and other fields. OFR254 requires the addition of externally formed O 3 since OH is formed from O 3 photolysis, while OFR185 does not since O 2 can be photolyzed to produce O 3 , and OH can also be formed from H 2 O photolysis. In this study, we use a plug-flow kinetic model to investigate OFR properties under a very wide range of conditions applicable to both field and laboratory studies. We show that the radical chemistry in OFRs can be characterized as a function of UV light intensity, H 2 O concentration, and total external OH reactivity (OHR ext , e.g., from volatile organic compounds (VOCs), NO x , and SO 2 ). OH exposure is decreased by added external OH reactivity. OFR185 is especially sensitive to this effect at low UV intensity due to low primary OH production. OFR254 can be more resilient against OH suppression at high injected O 3 (e.g., 70 ppm), as a larger primary OH source from O 3 , as well as enhanced recycling of HO 2 to OH, make external perturbations to the radical chemistry less significant. However if the external OH reactivity in OFR254 is much larger than OH reactivity from injected O 3 , OH suppression can reach 2 orders of magnitude. For a typical input of 7 ppm O 3 (OHR O 3 = 10 s −1 ), 10-fold OH suppression is observed at OHR ext ~ 100 s −1 , which is similar or lower than used in many laboratory studies. The range of modeled OH suppression for literature experiments is consistent with the measured values except for those with isoprene. The finding on OH suppression may have important implications for the interpretation of past laboratory studies, as applying OH exp measurements acquired under different conditions could lead to over a 1-order-of-magnitude error in the estimated OH exp . The uncertainties of key model outputs due to uncertainty in all rate constants and absorption cross-sections in the model are within ±25 % for OH exposure and within ±60 % for other parameters. These uncertainties are small relative to the dynamic range of outputs. Uncertainty analysis shows that most of the uncertainty is contributed by photolysis rates of O 3 , O 2 , and H 2 O and reactions of OH and HO 2 with themselves or with some abundant species, i.e., O 3 and H 2 O 2 . OH exp calculated from direct integration and estimated from SO 2 decay in the model with laminar and measured residence time distributions (RTDs) are generally within a factor of 2 from the plug-flow OH exp . However, in the models with RTDs, OH exp estimated from SO 2 is systematically lower than directly integrated OH exp in the case of significant SO 2 consumption. We thus recommended using OH exp estimated from the decay of the species under study when possible, to obtain the most appropriate information on photochemical aging in the OFR. Using HO x -recycling vs. destructive external OH reactivity only leads to small changes in OH exp under most conditions. Changing the identity (rate constant) of external OH reactants can result in substantial changes in OH exp due to different reductions in OH suppression as the reactant is consumed. We also report two equations for estimating OH exposure in OFR254. We find that the equation estimating OH exp from measured O 3 consumption performs better than an alternative equation that does not use it, and thus recommend measuring both input and output O 3 concentrations in OFR254 experiments. This study contributes to establishing a firm and systematic understanding of the gas-phase HO x and O x chemistry in these reactors, and enables better experiment planning and interpretation as well as improved design of future reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

A sub-pore model for multi-scale reaction–diffusion problems in porous media

Applications of reaction–diffusion systems in porous media pose a challenging problem for computational modeling approaches due to their multi-physics and multi-scale nature. The length scales usually span 3–4 orders of magnitude while physical phenomena involved include heat and mass transfer processes, and surface reactions. In this paper, a novel methodology that accounts for all the length scales and physical phenomena involved in a single framework is described. A length scale based dual approach is proposed – the larger pore channels (macro-pores) are resolved using conventional numerical techniques and a novel ‘sub-pore’ model is used to account for the unresolved pore channels (sub-pores) and the important physics therein. The porous network in the sub-pore system is composed of a fractal-like hierarchical system of straight cylindrical pores. Simplified governing equations for mass and energy transport are solved within the sub-pore system along with a reaction kinetics model to account for surface adsorption. An implicit coupling strategy is used to couple the macro-pore and the sub-pore systems so as to ensure conservation. The developed methodology is then applied to a few test cases and it is established that the proposed framework is necessary for problems where the adsorption time scale is much smaller than (diffusion-limited) or comparable to the diffusion time scale. Furthermore, it is also demonstrated that the framework can be potentially used to model the network of porous channels in its entirety thus significantly reducing computational costs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

LigninGraphs: lignin structure determination with multiscale graph modeling

Lignin is an aromatic biopolymer found in ubiquitous sources of woody biomass. Designing and optimizing lignin valorization processes requires a fundamental understanding of lignin structures. Experimental characterization techniques, such as 2D-heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectra, could elucidate the global properties of the polymer molecules. Computer models could extend the resolution of experiments by representing structures at the molecular and atomistic scales. We introduce a graph-based multiscale modeling framework for lignin structure generation and visualization. The framework employs accelerated rejection-free polymerization and hierarchical Metropolis Monte Carlo optimization algorithms. We obtain structure libraries for various lignin feedstocks based on literature and new experimental NMR data for poplar wood, pinewood, and herbaceous lignin. The framework could guide researchers towards feasible lignin structures, efficient space exploration, and future kinetics modeling. Its software implementation in Python, LigninGraphs, is open-source and available on GitHub.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

An Experimental Study of Uncertainty Considerations Associated with Predicting Auto-ignition Timing using Livengood-Wu Integral Method

The application of the Livengood-Wu (LW) integral method as a tool to estimate knock onset in spark ignited (SI) engines and combustion phasing in advanced compression ignition (ACI) engines has been demonstrated through simulations several times. In this study, the effect of uncertainties associated with parameters required for the LW integral method, when used as a tool for model based control of ignition timing in an ACI engine, were experimentally studied using five full boiling range gasoline fuels. As a first step, the method was applied to experimental data from a rapid compression machine and it was found that the ability of the LW integral method to predict ignition timing was very sensitive to the performance of the chemical kinetic model of each fuel. The method was subsequently applied to experimental data from a single-cylinder gasoline engine with simple approximations for the LW integral input parameters, and it was found that the predicted time of ignition was significantly different from the actual start of combustion. Systematic evaluation of various parametric uncertainties conducted thereafter showed that the uncertainty in cylinder charge temperature has the greatest influence. Improved methods of estimating cylinder charge temperature are proposed to account for the previously determined corrections, to enable the use of the LW integral method for model based control of ignition timing.

Compression Ignition↗
DOE OSTI2023

Mechanism and kinetic study of rare earth extraction from allanite by direct acid leaching

Direct acid leaching was conducted to recover rare earth elements (REEs) from an allanite ore assaying 9,861.4 ppm of REEs. The effects of operational variables such as acid type, sulfuric acid (H 2 SO 4 ) concentration, temperature, solid/liquid (S/L) ratio, and particle size on the total REE (TREE) recovery were systematically examined. The REE mineralogy and the involved leaching mechanisms were investigated by microscopic characterization and residue re-leaching tests. It was found that approximately 80% of REEs was extracted using 1 M H 2 SO 4 at 75 °C for 2 h. Increasing temperature significantly improved the REE extraction, while changing acid type, H 2 SO 4 concentration, and particle size did not. The TREE recovery also gradually decreased at higher S/L ratios. The low fraction of Si leached (<4.0%) suggested that the extraction did not significantly alter the crystal structure of the dominant mineral. The leaching kinetics of REEs were fast within the first 10 min, followed by a much slower stage (i.e., 10-120 min). The kinetic modeling results showed that each stage may be controlled by both chemical reaction and diffusion (i.e., mix control). The activation energies for REE leaching using 1 M H 2 SO 4 in the 0-10 min and 10-120 min were determined to be 20.3 and 10.8 kJ/mol, respectively. Finally, the residue re-leaching results suggested that the fast leaching kinetics at the beginning may be explained by the preferential decomposition of metamict allanite rather than the formation of silica gel, while the following slower leaching stage was due to the dissolution of difficult-to-leach REE-bearing minerals such as well-crystallized allanite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

HYDRATE v1.5 OPTION OF TOUGH+ v1.5

HYDRATE v1.5 is a numerical code that for the simulation of the behavior of hydrate-bearing geologic systems, and represents the third update of the code since its first release [Moridis et al., 2008]. It is an option of TOUGH+ v1.5 [Moridis and Pruess, 2014], a successor to the TOUGH2 [Pruess et al., 1999, 2012] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. HYDRATE v1.5 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. By solving the coupled equations of mass and heat balance, the fully operational TOUGH+HYDRATE code can model the non-isothermal gas release, phase behavior and flow of fluids and heat under conditions typical of common natural CH4-hydrate deposits (i.e., in the permafrost and in deep ocean sediments) in complex geological media at any scale (from laboratory to reservoir) at which Darcy's law is valid. TOUGH+HYDRATE v1.5 includes both an equilibrium and a kinetic model of hydrate formation and dissociation. The model accounts for heat and up to four mass components, i.e., water, CH4, hydrate, and water-soluble inhibitors such as salts or alcohols. These are partitioned among four possible phases (gas phase, liquid phase, ice phase and hydrate phase). Hydrate dissociation or formation, phase changes and the corresponding thermal effects are fully described, as are the effects of inhibitors. The model can describe all possible hydrate dissociation mechanisms, i.e., depressurization, thermal stimulation, salting-out effects and inhibitor-induced effects.

Moridis, GeorgeJ↗