Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Predictive Modeling”

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 829 records · Page 46

Physics-informed machine learning modeling for predictive control using noisy data

Due to the occurrence of over-fitting at the learning phase, the modeling of chemical processes via artificial neural networks (ANN) by using corrupted data (i.e., noisy data) is an ongoing challenge. Therefore, this work investigates the effect of both Gaussian and non-Gaussian noise on the performance of process-structure based recurrent neural networks (RNN) models, which take the form of partially-connected RNN models in this work, that are used to approximate a class of multi-input-multi-outputs nonlinear systems. Furthermore, two different techniques, specifically Monte Carlo dropout and co-teaching, are utilized in the development of partially-connected RNN models. Here, these two techniques are employed to reduce the over-fitting in ANNs when noisy data is used in the training process and, hence, to improve the open-loop accuracy as well as the closed-loop performance under a Lyapunov-based model predictive controller (MPC). Aspen Plus Dynamics, a well-known high-fidelity process simulator, is used to simulate a large-scale chemical process application in order to demonstrate the anticipated improvements in both open-loop approximation and closed-loop controller performance in the presence of Gaussian and non-Gaussian noise in the data set using physics-informed RNNs.

97 MATHEMATICS AND COMPUTING↗

A multiscale recurrent neural network model for predicting energy production from geothermal reservoirs

Optimization of energy production from geothermal reservoirs requires reliable prediction of energy production performance under alternative operation and development scenarios. Traditionally, reservoir simulation models are used for the evaluation and screening of alternative production and development plans. However, simulation models require extensive data collection and modeling efforts and are time-consuming to build, run, and update. Data-driven predictive models, on the other hand, can serve as efficient prediction tools that can be used for decision support and management of daily operations and surveillance activities. Data-driven models become particularly attractive when a reservoir simulation model for a field does not exist and/or is difficult to build. Machine learning (ML)-based data-driven models that have recently become popular in several fields exploit statistical patterns and relations in training data to generate predictions. As such, they tend to perform better in interpolation problems (that is, prediction within the training data range) than when they are used to extrapolate beyond the training data. Production data from geothermal reservoirs tend to exhibit short-term variabilities as well as long-term trends, such as monotonically declining production temperatures. Capturing both short-term features and long-term trends with ML-based models is not trivial. We evaluate the use of recurrent neural networks (RNN) for the prediction of energy production from geothermal reservoirs. RNN is a class of ML architectures that are used to represent and predict sequential/dynamic data. Thus, it can be challenging to apply RNN to problems where long-term trends must be captured and extrapolation beyond the training data range is needed. We introduce the multiscale RNN architecture to extend the application of RNN to detect and predict both short-term variabilities and long-term trends in geothermal data. The developed architecture consists of a long-term component to only capture low-frequency data patterns, and a short-term component to detect features with higher frequency and more nonlinearity. The final prediction is obtained by combining the long-term and short-term predictions. Both synthetic and field data are used to evaluate the presented multiscale RNN model. The prediction performance of the multiscale RNN is compared against those obtained from the regular RNN and the autoregressive (AR) model. The results suggest that the multiscale architecture improves the long-term prediction performance of the regular RNN and enhances its robustness against noise.

15 GEOTHERMAL ENERGY↗

Development of a fourth generation predictive capability maturity model.

The Predictive Capability Maturity Model (PCMM) is an expert elicitation tool designed to characterize and communicate completeness of the approaches used for computational model definition, verification, validation, and uncertainty quantification associated for an intended application. The primary application of this tool at Sandia National Laboratories (SNL) has been for physics-based computational simulations in support of nuclear weapons applications. The two main goals of a PCMM evaluation are 1) the communication of computational simulation capability, accurately and transparently, and 2) the development of input for effective planning. As a result of the increasing importance of computational simulation to SNL s mission, the PCMM has evolved through multiple generations with the goal to provide more clarity, rigor, and completeness in its application. This report describes the approach used to develop the fourth generation of the PCMM.

Hills, Richard Guy↗

SU-F-BRD-01: A Logistic Regression Model to Predict Objective Function Weights in Prostate Cancer IMRT

Purpose: To develop a statistical model that predicts optimization objective function weights from patient geometry for intensity-modulation radiotherapy (IMRT) of prostate cancer. Methods: A previously developed inverse optimization method (IOM) is applied retrospectively to determine optimal weights for 51 treated patients. We use an overlap volume ratio (OVR) of bladder and rectum for different PTV expansions in order to quantify patient geometry in explanatory variables. Using the optimal weights as ground truth, we develop and train a logistic regression (LR) model to predict the rectum weight and thus the bladder weight. Post hoc, we fix the weights of the left femoral head, right femoral head, and an artificial structure that encourages conformity to the population average while normalizing the bladder and rectum weights accordingly. The population average of objective function weights is used for comparison. Results: The OVR at 0.7cm was found to be the most predictive of the rectum weights. The LR model performance is statistically significant when compared to the population average over a range of clinical metrics including bladder/rectum V53Gy, bladder/rectum V70Gy, and mean voxel dose to the bladder, rectum, CTV, and PTV. On average, the LR model predicted bladder and rectum weights that are both 63% closer to the optimal weights compared to the population average. The treatment plans resulting from the LR weights have, on average, a rectum V70Gy that is 35% closer to the clinical plan and a bladder V70Gy that is 43% closer. Similar results are seen for bladder V54Gy and rectum V54Gy. Conclusion: Statistical modelling from patient anatomy can be used to determine objective function weights in IMRT for prostate cancer. Our method allows the treatment planners to begin the personalization process from an informed starting point, which may lead to more consistent clinical plans and reduce overall planning time.

60 APPLIED LIFE SCIENCES↗

A Discrete-Time Average Model Based Predictive Control for Quasi-Z-Source Inverter

A discrete-time average model-based predictive control (DTA-MPC) is proposed for a quasi-Z-source inverter (qZSI). As a single-stage inverter topology, the qZSI regulates the dc-link voltage and the ac output voltage through the shoot-through (ST) duty cycle and the modulation index. Several feedback strategies have been dedicated to produce these two control variables, among which the most popular are the proportional–integral (PI)-based control and the conventional model-predictive control (MPC). However, in the former, there are tradeoffs between fast response and stability; the latter is robust, but at the cost of high calculation burden and variable switching frequency. Moreover, they require an elaborated design or fine tuning of controller parameters. The proposed DTA-MPC predicts future behaviors of the ST duty cycle and modulation signals, based on the established discrete-time average model of the quasi-Z-source (qZS) inductor current, the qZS capacitor voltage, and load currents. The prediction actions are applied to the qZSI modulator in the next sampling instant, without the need of other controller parameters’ design. A constant switching frequency and significantly reduced computations are achieved with high performance. Transient responses and steady-state accuracy of the qZSI system under the proposed DTA-MPC are investigated and compared with the PI-based control and the conventional MPC. Simulation and experimental results verify the effectiveness of the proposed approach for the qZSI.

42 ENGINEERING↗

Optimized Machine Learning Model for Predicting Groundwater Contamination

The use of physical models to predict groundwater contaminant movement remains technically challenging due to the complexity of the phenomena, the heterogeneity of key parameters in nature, and the presence of poorly defined interactive and feedback processes. New approaches to address these challenges are needed. In this study, we evaluate various Artificial Intelligence (AI)-based approaches to understand a hexavalent chromium (Cr(VI)) plumes located on the U.S. Department of Energy’s (DOE) Hanford Site in Richland, WA. The groundwater monitoring dataset used in this study included data from the 100 Area along the Columbia River and included data collected between 2010 to 2019. This study investigates the most prominent contaminant, Cr(VI), with the Extreme Gradient Boosting (XGBoost) machine learning model. The XGBoost models were compared with optimized versions using an Empirical Bayes Search Cross-Validation technique for better prediction. The optimized XGBoost model yielded an R^2 value of 0.99 on the training set and 0.85 on the testing set, whereas XGBoost without optimization yielded a value of 0.83 on the training set and 0.85 on the testing set. This paper provides an overview of a computational method for groundwater contamination modeling that shows promise for improving current remediation efforts.

Mazumdar, Hirak↗

Combination of searches for heavy resonances decaying into bosonic and leptonic final states using 36 fb - 1 of proton-proton collision data at s = 13 TeV with the ATLAS detector

Searches for new heavy resonances decaying into different pairings of W, Z, or Higgs bosons, as well as directly into leptons, are presented using a data sample corresponding to 36.1 fb -1 of pp collisions at $ \sqrt{s}=13 $ TeV collected during 2015 and 2016 with the ATLAS detector at the CERN Large Hadron Collider. Analyses selecting bosonic decay modes in the qqqq, ννqq, ℓνqq, ℓℓqq, ℓνℓν, ℓℓνν, ℓνℓℓ, ℓℓℓℓ, qqbb, ννbb, ℓνbb, and ℓℓbb final states are combined, searching for a narrow-width resonance. Likewise, analyses selecting the leptonic ℓν and ℓℓ final states are also combined. These two sets of analyses are then further combined. No significant deviation from the Standard Model predictions is observed. Three benchmark models are tested: a model predicting the existence of a new heavy scalar singlet, a simplified model predicting a heavy vector-boson triplet, and a bulk Randall-Sundrum model with a heavy spin-2 Kaluza-Klein excitation of the graviton. Cross section limits are set at the 95% confidence level using an asymptotic approximation and are compared with predictions for the benchmark models. These limits are also expressed in terms of constraints on couplings of the heavy vector-boson triplet to quarks, leptons, and the Higgs boson. The data exclude a heavy vector-boson triplet with mass below 5.5 TeV in a weakly coupled scenario and 4.5 TeV in a strongly coupled scenario, as well as a Kaluza-Klein graviton with mass below 2.3 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Searches for heavy diboson resonances in pp collisions at $ \sqrt{s}=13 $ TeV with the ATLAS detector

Searches for new heavy resonances decaying to WW, WZ, and ZZ bosons are presented, using a data sample corresponding to 3.2 fb -1 of pp collisions at s=13 TeV colle cted with the ATLAS detector at the CERN Large Hadron Collider. Analyses selecting ννqq, ℓνqq, ℓℓqq and qqqq final states are combined, searching for an arrow-width resonance with mass between 500 and 3000 GeV. The discriminating variable is either an invariant mass or a transverse mass. No significant deviations from the Standard Model predictions are observed. Three benchmark models are tested: a model predicting the existence of a new heavy scalar singlet, a simplified model predicting a heavy vector-boson triplet, and a bulk Randall-Sundrum model with a heavy spin-2 graviton. Cross-section limits are set at the 95% confidence level and are compared to theoretical cross-section predictions for a variety of models. The data exclude a scalar singlet with mass below 2650 GeV, a heavy vector-boson triplet with mass below 2600 GeV, and a graviton with mass below 1100 GeV. These results significantly extend the previous limits set using pp collisions at √s=8 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Aluminum cladding oxide growth prediction for high flux research reactors

Aluminum cladding oxidation of research-reactor fuel elements at high power conditions has a disadvantageous effect on fuel performance due to the lower thermal conductivity of the oxide. The oxide growth prediction models available in the literature were mostly developed for low power conditions. To examine the applicability of the models to high power and high temperature test conditions, the models were studied by coupling with the most frequently employed heat transfer coefficient (HTC) correlations including the Dittus-Boelter correlation, the Colburn correlation, the Sieder-Tate correlation, and KAERI-developed correlation. The Griess model over-predicted the oxide growth while the KAERI-Griess model under-predicted the oxide growth for high power tests. The Kim model, coupled with the Colburn correlation, gave most consistent results with the measured data from two BR2 experiments. However, the Kim model was found to be inapplicable to the EUHFRR conditions at the peak power locations if it was coupled with the Dittus-Boelter correlation. A revision of the prediction models to more closely agree with the measured data was recommended. Furthermore, AG3NE and AlFeNi cladding types were tested in the E-FUTURE experiment, and a noticeable (although small) reduction in oxide thickness on the AlFeNi cladding was observed. However, this difference was believed to be only a secondary effect considering other uncertainties in model predictions, so no attempt was made to model the alloying effect.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development and Evaluation of Occupancy-Aware HVAC Control for Residential Building Energy Efficiency and Occupant Comfort

Occupancy-aware heating, ventilation, and air conditioning (HVAC) control offers the opportunity to reduce energy use without sacrificing thermal comfort. Residential HVAC systems often use manually-adjusted or constant setpoint temperatures, which heat and cool the house regardless of whether it is needed. By incorporating occupancy-awareness into HVAC control, heating and cooling can be used for only those time periods it is needed. Yet, bringing this technology to fruition is dependent on accurately predicting occupancy. Non-probabilistic prediction models offer an opportunity to use collected occupancy data to predict future occupancy profiles. Smart devices, such as a connected thermostat, which already include occupancy sensors, can be used to provide a continually growing collection of data that can then be harnessed for short-term occupancy prediction by compiling and creating a binary occupancy prediction. Real occupancy data from six homes located in Colorado is analyzed and investigated using this occupancy prediction model. Results show that non-probabilistic occupancy models in combination with occupancy sensors can be combined to provide a hybrid HVAC control with savings on average of 5.0% and without degradation of thermal comfort. Model predictive control provides further opportunities, with the ability to adjust the relative importance between thermal comfort and energy savings to achieve savings between 1% and 13.3% depending on the relative weighting between thermal comfort and energy savings. In all cases, occupancy prediction allows the opportunity for a more intelligent and optimized strategy to residential HVAC control.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Liquid Core Detection and Strand Condition Monitoring in a Continuous Caster Using Optical Fiber

Real-time monitoring of the liquid core position during the continuous casting of steel has been demonstrated using low-cost distributed optical-fiber-based strain sensors. These sensors were installed on the containment roll support structures in the segments of a production continuous caster to detect the position of the solid–liquid interface and monitor the strand condition during the continuous casting. Distributed Fiber Bragg Grating sensors (FBGs) were used in this work to monitor strain at six roll positions in the caster. The sensor performance was first validated by comparing optical strain measurements with conventional strain gauge measurements in the lab. Next, optical strain measurements were performed on an isolated caster segment in a segment maintenance facility using hydraulic jacks to simulate the presence of a liquid core under the roll. Finally, the sensors were evaluated during caster operation. The sensors successfully detected the load increase associated with the presence of a liquid core under each instrumented roll location. Incidents of bulging and roll eccentricity were also detected using frequency analysis of the optical strain signal. The liquid core position measurements were compared using predictions from computer models (digital twins) in use at the production site. The measurements were in good agreement with the model predictions, with a few exceptions. Under certain transient caster operating conditions, such as spraying practice changes and SEN exchanges, the model predictions deviated slightly from the liquid core position determined from strain measurements.

47 OTHER INSTRUMENTATION↗

Mean Climate and Tropical Rainfall Variability in Aquaplanet Simulations Using the Model for Prediction Across Scales‐Atmosphere

Abstract Aquaplanet experiments are important tools for understanding and improving physical processes simulated by global models; yet, previous aquaplanet experiments largely differ in their representation of subseasonal tropical rainfall variability. This study presents results from aquaplanet experiments produced with the Model for Prediction Across Scales‐Atmosphere (MPAS‐A)—a community model specifically designed to study weather and climate in a common framework. The mean climate and tropical rainfall variability simulated by MPAS‐A with varying horizontal resolution were compared against results from a recent suite of aquaplanet experiments. This comparison shows that, regardless of horizontal resolution, MPAS‐A produces the expected mean climate of an aquaplanet framework with zonally symmetric but meridionally varying sea‐surface temperature. MPAS‐A, however, has a stronger signal of tropical rainfall variability driven by convectively coupled equatorial waves. Sensitivity experiments with different cumulus parameterizations, physics packages, and vertical grids consistently show the presence of those waves, especially equatorial Kelvin waves, in phase with lower‐tropospheric convergence. Other models do not capture such rainfall‐kinematics phasing. These results suggest that simulated tropical rainfall variability depends not only on the cumulus parameterization (as suggested by previous studies) but also on the coupling between physics and dynamics of climate and weather prediction models.

Rios‐Berrios, R.↗

Modeling the effectiveness of remediation efforts in contaminated urban areas: An EMRAS II Urban Areas Working Group exercise - 15631

The Urban Areas Working Group was part of the IAEA's EMRAS II (Environmental Modelling for Radiation Safety) Programme during 2009-2011. The Working Group's goal was to test and improve the capabilities of models used in assessment of radioactive contamination in urban settings. The 'contaminant transport and countermeasures' exercise was focused on modeling the effectiveness of potential countermeasures or remediation efforts with respect to reduction of human exposures and doses. The exercise was based on a hypothetical situation involving an initial concentration of Co-60 or Pu-239 in air in a specified urban area. Participants were given detailed input information, including descriptions of test locations in two types of urban areas, (1) a business area with buildings and asphalt, and (2) a park area near an apartment town. Several kinds of initial weather conditions (dry, light rain, and heavy rain) were explored. Modeling endpoints included contamination densities, dose rates (from each contributing surface and total), annual and cumulative doses (internal and external) for specified reference individuals, and effectiveness of countermeasures and remediation efforts in terms of dose rate reduction and dose reduction. External doses were emphasized for Co-60 and internal doses for Pu-239. Approaches and model predictions from five models are compared, and explanations for similarities and differences among model predictions are offered. Understanding the effects of different approaches, assumptions, and parameter values is essential to a comparison of results from different models and to understanding both the usefulness and the limitations of any given model. (authors)

54 ENVIRONMENTAL SCIENCES↗

Collaborative Research: Separating Forced and Unforced Decadal Predictability in Models and Observations

This report is a progress report of the accomplishments of the research grant “Collaborative Research: Separating Forced and Unforced Decadal Predictability in Models and Observations” during the period 1 May 2011- 31 August 2013. This project is a collaborative one between Columbia University and George Mason University. George Mason University will submit a final technical report at the conclusion of their no-cost extension. The purpose of the proposed research is to identify unforced predictable components on decadal time scales, distinguish these components from forced predictable components, and to assess the reliability of model predictions of these components. Components of unforced decadal predictability will be isolated by maximizing the Average Predictability Time (APT) in long, multimodel control runs from state-of-the-art climate models. Components with decadal predictability have large APT, so maximizing APT ensures that components with decadal predictability will be detected. Optimal fingerprinting techniques, as used in detection and attribution analysis, will be used to separate variations due to natural and anthropogenic forcing from those due to unforced decadal predictability. This methodology will be applied to the decadal hindcasts generated by the CMIP5 project to assess the reliability of model projections. The question of whether anthropogenic forcing changes decadal predictability, or gives rise to new forms of decadal predictability, also will be investigated.

54 ENVIRONMENTAL SCIENCES↗

PNNL Technical Support to The Implementation of EMTA and EMTA-NLA Models in Autodesk ® Moldflow ® Packages

Under the Predictive Engineering effort, PNNL developed linear and nonlinear property prediction models for long-fiber thermoplastics (LFTs). These models were implemented in PNNL’s EMTA and EMTA-NLA codes. While EMTA is a standalone software for the computation of the composites thermoelastic properties, EMTA-NLA presents a series of nonlinear models implemented in ABAQUS ® via user subroutines for structural analyses. In all these models, it is assumed that the fibers are linear elastic while the matrix material can exhibit a linear or typical nonlinear behavior depending on the loading prescribed to the composite. The key idea is to model the constitutive behavior of the matrix material and then to use an Eshelby-Mori-Tanaka approach (EMTA) combined with numerical techniques for fiber length and orientation distributions to determine the behavior of the as-formed composite. The basic property prediction models of EMTA and EMTA-NLA have been subject for implementation in the Autodesk ® Moldflow ® software packages. These models are the elastic stiffness model accounting for fiber length and orientation distributions, the fiber/matrix interface debonding model, and the elastic-plastic models. The PNNL elastic-plastic models for LFTs describes the composite nonlinear stress-strain response up to failure by an elastic-plastic formulation associated with either a micromechanical criterion to predict failure or a continuum damage mechanics formulation coupling damage to plasticity. All the models account for fiber length and orientation distributions as well as fiber/matrix debonding that can occur at any stage of loading. In an effort to transfer the technologies developed under the Predictive Engineering project to the American automotive and plastics industries, PNNL has obtained the approval of the DOE Office of Vehicle Technologies to provide Autodesk, Inc. with the technical support for the implementation of the basic property prediction models of EMTA and EMTA-NLA in the Autodesk ® Moldflow ® packages. This report summarizes the recent results from Autodesk Simulation Moldlow Insight (ASMI) analyses using the EMTA models and EMTA-NLA/ABAQUS ® analyses for further assessment of the EMTA-NLA models to support their implementation in Autodesk Moldflow Structural Alliance (AMSA). PNNL’s technical support to Autodesk, Inc. included (i) providing the theoretical property prediction models as described in published journal articles and reports, (ii) providing explanations of these models and computational procedure, (iii) providing the necessary LFT data for process simulations and property predictions, and (iv) performing ABAQUS/EMTA-NLA analyses to further assess and illustrate the models for selected LFT materials.

damage analysis↗

Validation of mathematical models for the prediction of organs-at-risk dosimetric metrics in high-dose-rate gynecologic interstitial brachytherapy

Purpose: Given the complicated nature of an interstitial gynecologic brachytherapy treatment plan, the use of a quantitative tool to evaluate the quality of the achieved metrics compared to clinical practice would be advantageous. For this purpose, predictive mathematical models to predict the D{sub 2cc} of rectum and bladder in interstitial gynecologic brachytherapy are discussed and validated.Methods: Previous plans were used to establish the relationship between D2cc and the overlapping volume of the organ at risk with the targeted area (C0) or a 1-cm expansion of the target area (C1). Three mathematical models were evaluated: D{sub 2cc}=α*C{sub 1}+β (LIN); D{sub 2cc}=α– exp(–β*C{sub 0}) (EXP); and a mixed approach (MIX), where both C{sub 0} and C{sub 1} were inputs of the model. The parameters of the models were optimized on a training set of patient data, and the predictive error of each model (predicted D{sub 2cc}− real D{sub 2cc}) was calculated on a validation set of patient data. The data of 20 patients were used to perform a K-fold cross validation analysis, with K = 2, 4, 6, 8, 10, and 20.Results: MIX was associated with the smallest mean prediction error <6.4% for an 18-patient training set; LIN had an error <8.5%; EXP had an error <8.3%. Best case scenario analysis shows that an error ≤5% can be achieved for a ten-patient training set with MIX, an error ≤7.4% for LIN, and an error ≤6.9% for EXP. The error decreases with the increase in training set size, with the most marked decrease observed for MIX.Conclusions: The MIX model can predict the D{sub 2cc} of the organs at risk with an error lower than 5% with a training set of ten patients or greater. The model can be used in the development of quality assurance tools to identify treatment plans with suboptimal sparing of the organs at risk. It can also be used to improve preplanning and in the development of real-time intraoperative planning tools.

61 RADIATION PROTECTION AND DOSIMETRY↗

Status of Chronic Oxidation Studies of Graphite

Graphite will undergo extremely slow, but continuous oxidation by traces of moisture that will be present, albeit at very low levels, in the helium coolant of HTGR. This chronic oxidation may cause degradation of mechanical strength and thermal properties of graphite components if a porous oxidation layer penetrates deep enough in the bulk of graphite components during the lifetime of the reactor. The current research on graphite chronic oxidation is motivated by the acute need to understand the behavior of each graphite grade during prolonged exposure to high temperature chemical attack by moisture. The goal is to provide the elements needed to develop predictive models for long-time oxidation behavior of graphite components in the cooling helium of HTGR. The tasks derived from this goal are: (1) Oxidation rate measurements in order to determine and validate a comprehensive kinetic model suitable for prediction of intrinsic oxidation rates as a function of temperature and oxidant gas composition; (2) Characterization of effective diffusivity of water vapor in the graphite pore system in order to account for the in-pore transport of moisture; and (3) Development and validation of a predictive model for the penetration depth of the oxidized layer, in order to assess the risk of oxidation caused damage of particular graphite grades after prolonged exposure to the environment of helium coolant in HTGR. The most important and most time consuming of these tasks is the measurement of oxidation rates in accelerated oxidation tests (but still under kinetic control) and the development of a reliable kinetic model. This report summarizes the status of chronic oxidation studies on graphite, and then focuses on model development activities, progress of kinetic measurements, validation of results, and improvement of the kinetic models. Analysis of current and past results obtained with three grades of showed that the classical Langmuir-Hinshelwood model cannot reproduce all data collected so far. Starting from here we propose a modification of the LH model to include temperature activation of graphite surface as a Boltzmann activation function. The enhanced Boltzmann-Langmuir-Hinshelwood model (BLH) was tested successfully on three grades of graphite. The model is a robust, comprehensive mathematical function that allows better fitting of experimental results spanning a wide range of temperature and partial pressures of water vapor and hydrogen. However, the model did not fit satisfactorily the data extracted from the old report on graphite H-451 oxidation by water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗