Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “well data”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Scale effects on core design, fuel costs, and spent fuel volume of pressurized water reactors

The desire to improve the economic competitiveness and deployment pace of nuclear energy through modularization, manufacturing, and series production had led to the development of smaller size reactors. As the standard 17x17 fuel technology is mainly maintained in the pressurized water reactors (PWRs) category, this translates into a lower number of fuel assemblies in the core and sometimes a reduced fuel height. To assess the impact of such scale change in core design on fuel cycle cost and spent fuel volume, a scoping analysis tool is developed based on infinite lattice calculations, leakage, fuel management reduced models, and levelized unit cost of electricity (LCOE) estimate. As such, cost dynamics driven by fuel specific power, burnup, core leakage, feed, cycle length, fuel assembly height as well as uranium market data are captured with consistent set of assumptions and analysis methods. A selection of 5 reactor designs representative of leading PWR developers is assessed and compared. Pursuing higher specific powers and optimal burnups are highlighted as the main fuel cost reduction drivers, nevertheless, practical limitations and opportunities must be evaluated to establish the feasibility of such enhanced fuel operation. In consequence, a detailed core design is performed using SIMULATE3 code for 5 PWR variations including natural and forced coolant circulation modes, two reactor scales, power densities of 73, 112, and 123 kW/l and higher discharge burnups. Design and optimization are performed at the lattice level, for the reflector, and at the core loading level. Satisfactory steady-state operation including power distribution, coolant operating limits, and reactivity requirements are analyzed and reported in this paper. The fuel economics of the detailed core designs confirm the scoping analysis findings. Despite the unlocked power uprates in small PWRs, the achievable burnup for a given fuel specific power requires more enrichment and shorter fuel height results in higher fabrication costs per mass of fuel, which makes scaling down core size a more expensive endeavor on the fuel cycle front. Spent fuel volumes are reported for the PWRs designed in this paper. Furthermore, these volumes are driven by the core average discharge burnup regardless of the scale in consideration. Additional cost and core performance aspects related to heavy reflector gains, fuel-reflector substitution, and disposal cost policy in the U.S. are examined.

42 ENGINEERING↗
DOE OSTI2020

Online adaptive basis refinement and compression for reduced-order models via vector-space sieving

In many applications, projection-based reduced-order models (ROMs) have demonstrated the ability to provide rapid approximate solutions to high-fidelity full-order models (FOMs). However, there is no a priori assurance that these approximate solutions are accurate; their accuracy depends on the ability of the low-dimensional trial basis to represent the FOM solution. As a result, ROMs can generate inaccurate approximate solutions, e.g., when the FOM solution at the online prediction point is not well represented by training data used to construct the trial basis. To address this fundamental deficiency of standard model-reduction approaches, this work proposes a novel online-adaptive mechanism for efficiently enriching the trial basis in a manner that ensures convergence of the ROM to the FOM, yet does not incur any FOM solves. The mechanism is based on the previously proposed adaptive $h$-refinement method for ROMs, but improves upon this work in two crucial ways. First, the proposed method enables basis refinement with respect to any orthogonal basis (not just the Kronecker basis), thereby generalizing the refinement mechanism and enabling it to be tailored to the physics characterizing the problem at hand. Second, the proposed method provides a fast online algorithm for periodically compressing the enriched basis via an efficient proper orthogonal decomposition (POD) method, which does not incur any operations that scale with the FOM dimension. We report these two features allow the proposed method to serve as (1) a failsafe mechanism for ROMs, as the method enables the ROM to satisfy any prescribed error tolerance online (even in the case of inadequate training), and (2) an efficient online basis-adaptation mechanism, as the combination of basis enrichment and compression enables the basis to adapt online while controlling its dimension.

42 ENGINEERING↗
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Plausible indication of gamma-ray absorption by dark matter in NGC 1068

NGC 1068 is the brightest extragalactic source in high-energy neutrinos as seen by IceCube, yet the accompanying gamma-ray flux is orders of magnitude weaker. It has been argued that this indicates that the bulk of neutrinos and gamma rays are emitted in the innermost vicinity of the central supermassive black hole, which is transparent to neutrinos, but opaque to gamma rays. Even in such extreme scenarios for the acceleration of cosmic rays, astrophysical models typically overestimate the low-energy gamma-ray flux and/or require some fine-tuning in the physical parameters. Here, in this work, we suggest instead that the dark matter surrounding the supermassive black hole may absorb the gamma rays, inducing the observed deficit. We show that for a dark matter–photon scattering cross section in the range σ DM−γ /m DM ≃ 10 −28 −10 −30 cm 2 /GeV, Fermi-LAT measurements can be well reconciled with IceCube data. We also present some simple particle physics examples that achieve the correct spectral energy dependence while respecting complementary constraints.

Active Galactic Nuclei↗
DOE OSTI2020

Numerical study on spray collapse process of ECN spray G injector under flash boiling conditions

This work features a numerical study of the spray collapse process that occurs under flash-boiling, gasoline-direct-injection conditions. The Engine Combustion Network's Spray G injector was modeled, and iso-octane, n-hexane, and propane were used as fuels with the goal to investigate fuel property effects. The development of individual jets and jet-to-jet interaction were analyzed to reveal the physical mechanisms of spray collapse. The results agreed qualitatively well with the experimental data of propane available from the literature and the predictions revealed the model's ability to capture the near-nozzle shock structures witnessed in the experiments. Simulations showed that the flow development of flashing sprays was strongly tied with the fuel properties. Specifically, under-expanded jets were not observed with iso-octane due to the fuel's low vapor pressure. In the case of n-hexane, the jets were under-expanded, and a low-pressure core was observed inside the jets due to the flow expansion. With propane, the very strong expansion caused the flow to become transonic and generate a unique shock structure near the nozzle exit. The shock structure revealed the formation of primary shock cells within the jets and secondary shock cells induced by jet-to-jet interaction. Based on the results, the spray collapse under flash boiling conditions was mainly attributed to two effects, i.e., the decrease in pressure at the spray center induced by the low-pressure cores inside and between the jets, and the isolation effect due to the jet-to-jet interaction which contributed to maintain the pressure differential between the spray center and the ambient.

42 ENGINEERING↗
DOE OSTI2021

Fatigue behavior of fine-grained magnesium under tension-tension loading at 0 °C

Here, this study investigated low-temperature fatigue behavior of fine-grained magnesium through tension–tension fatigue testing at the temperature of 0 °C. The applied stress profiles of fatigue testing were selected to cover a wide range of loading conditions including low, medium, and high loading stresses. Even when the maximum applied stress of fatigue testing was as high as 88% of the ultimate tensile strength of the material, the corresponding fatigue life was long, i.e., 3973 cycles. Deformation accumulation per loading cycle reduced with the progression of fatigue testing for all studied loading conditions. Fatigue life analyses were performed using both Smith-Watson-Topper (SWT) model and Jahed-Varvani energy-based (JV) model. The fatigue life predictions through both models agree well with the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting ductility of Mg/SiCp nanocomposite under multiaxial loading conditions based on unit cell modeling

This article presents an investigation on the ductility of Mg/SiCp under nine different multiaxial loading conditions based on a unit cell (or called representative volumetric element, RVE) finite element model. The constitutive relationship of magnesium matrix is characterized by a model combined Lode angle dependent plasticity model with the CPB06 plasticity model, which is able to describe the tension-compression asymmetry in magnesium matrix. The debonding of the matrix-particle interface and matrix fracture are considered by the modified Mohr-Coulomb (MMC) fracture model. Simulation results show that the interface debonding and subsequent matrix fracture are responsible for the ductility and fracture of nanocomposite. The computational results with parametric studies are compared and well correlated with experimental data, and it finally gives a stress-based full fracture locus of Mg/SiCp nanocomposite.

42 ENGINEERING↗
DOE OSTI2025

Geographic source attribution of honey by strontium isotope analyses: Latvia and India measurements compared to model predictions in a feasibility study

Strontium isotope ratios ( 87 Sr/ 86 Sr) can be used to determine the geographic origin of agricultural products. Here, in this study, we measured the 87 Sr/ 86 Sr of honey samples from Richland, WA, USA, Latvia, and India, and compared the results to published 87 Sr/ 86 Sr from surrounding areas and values predicted based upon a random forest isoscape model (Bataille et al., 2020). While the 87 Sr/ 86 Sr of the honey samples compared well with previously published data, the model did not accurately predict the 87 Sr/ 86 Sr of the honey samples, demonstrating that further refinement of the model would be beneficial. Because Sr cycling in terrestrial ecosystems is complex, the 87 Sr/ 86 Sr of honey is likely to reflect several sources (e.g., nectar, surface water, dust) and may change over time. In order to accurately predict the 87 Sr/ 86 Sr of honey, models must consider local geology, ecology, and bee behavior. Additional research examining how Sr is incorporated into honey is needed to inform models.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Fission gas diffusion and release for Cr 2 O 3 -doped UO 2 : From the atomic to the engineering scale

Here, the anticipated benefits of large grains in Cr 2 O 3 -doped UO 2 pellets include improved mechanical and fission gas retention properties. To support the assessment of fission gas release (FGR) from doped pellets, the impact of doping on fission gas diffusivity for in-reactor conditions must be understood. In this work, we tackle this issue by informing the fission gas model within the BISON fuel performance code using material models developed at the atomic scale. The investigation of intra-granular fission gas diffusivity in Cr 2 O 3 -doped UO 2 is carried out by adapting a cluster dynamics model that, accounting for UO 2 thermochemistry, is capable of describing Xe diffusion under irradiation in undoped UO 2 as the starting point. Using a thermodynamic analysis, it is shown that in stoichiometric UO 2 with additions of Cr 2 O 3 the oxygen potential is defined by the Cr-Cr 2 O 3 two-phase equilibrium. Using the cluster dynamics model, the predicted Xe diffusivity in doped UO 2 was significantly increased in both the intrinsic and irradiation-enhanced regimes compared to undoped UO 2 as a result of higher concentrations of uranium and oxygen vacancies, respectively. This is a consequence of the more oxidizing conditions at high temperature, and more reducing conditions at low temperature, as a result of doping. Arrhenius functions have been fitted to the cluster dynamics results to enable implementation of the new diffusivities in the BISON fission gas behavior model. BISON simulations were carried out, showing the competing effects of the enlarged grains and the new fission gas diffusivity model, which act to suppress and enhance fission gas release, respectively. The new physics-informed model was validated against in-reactor experimental measurements under normal operation. Additionally, benchmarking was carried out for power ramp conditions. The predicted fission gas release agreed well with the experimental data, showing noticeable improvements over the standard UO 2 model.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling fission product diffusion in TRISO fuel particles with BISON

Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Hydrogen effects on thermal diffusivity and electrical resistivity of zircaloy cladding

Thermal diffusivity measurements on Zirconium-based cladding materials have historically been a challenge due to the difficulty to measure on specimens with curved geometries, including nuclear grade Zircaloy cladding materials. Here, in this work, we first used laser flash analysis method and four-probe configuration method to measure the thermal diffusivity and electrical resistivity of Zircaloy tubes respectively, which show good agreement with Zircaloy plates in this work, as well as previously published data. The consistent results proved the applicability of the laser flash analysis setup and four-probe configuration method for investigating thermal diffusivity and electrical resistivity of Zircaloy tubes. We further investigated the hydrogen effect on thermal diffusivity and electrical resistivity of Zircaloy. Hydrogen plays significant roles in the thermal diffusivity of Zircaloy, which depends on the hydrogen concentration. For higher hydrogen contents (1130 and 1820 wppm in this work), where phase transformation (α-Zr + δ-hydride → α-Zr + β-Zr) occurs at 567 °C, thermal diffusivity decreases as a function of temperature at the α-Zr + δ-hydride phase regime, while increase at higher temperature at the α-Zr + β-Zr phase regime. For low hydrogen concentration, where hydride dissolved into α-Zr matrix phase at higher temperature, the thermal diffusivity is lower than non-hydrided Zircaloy-4 with similar temperature-dependent trends. Such observation demonstrates the hydrogen effects on reducing thermal diffusivity of α-Zr phase. Hydrogen increases the electrical resistivity of Zircaloy. Similar to thermal diffusivity results, the phase transformation causes a reversal of temperature-dependent trends in electrical resistivity results in Zircaloy-4 with higher hydrogen contents.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamic recrystallization of a wrought magnesium alloy: Grain size and texture maps and their application for mechanical behavior predictions

Effect of dynamic recrystallization (DRX) on the grain refinement and texture modification was studied by conducting a series of hot compression on AZ31B Mg alloy. Processing-grain size-texture maps were established as a function of temperature, strain rate, strain, and the Zener-Hollomon parameter (Z). Moreover, influence of simultaneous changes in the grain size and texture on tensile yield strength and ductility of the hot-worked Mg alloy was studied using the grain size and texture maps established. The processing-grain size map showed that the DRX grain size decreased with the increase in Z, also revealing various characteristics ranging from a grain growth, grain refinement, to a bimodal distribution of ultrafine grains and partially recrystallized grains. The effect of twinning on the grain refinement was also evident at high Z conditions. The processing-texture map revealed that the initial fiber texture was altered significantly to a shear, off-normal, or extension-twin texture with the increase in Z. The Schmid-factor maps were calculated to identify a dominant deformation mode during a subsequent tensile deformation of the hot-worked samples. Finally, a corresponding Hall-Petch relationship for each dominant deformation mode was used to establish tensile yield-strength maps, which agree well with the measured data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Co-ion specific effect on sodium halides sorption and transport in a cross-linked poly(p-styrene sulfonate-co-divinylbenzene) for membrane applications

Here, sodium halide sorption and transport in a cation exchange membrane based on cross-linked poly(p-styrene sulfonate-co-divinylbenzene) (CR61) were investigated experimentally and theoretically. At fixed co-ion valence, sodium halide partitioning between aqueous electrolyte solutions and ion exchange polymers is only minimally sensitive to co-ion chemical identity. NaCl, NaF and NaI exhibit similar solubility coefficients in CR61, which agrees reasonably well with predictions of the Donnan-Manning model, especially at high external salt concentrations. Ion activity coefficients in the membrane were evaluated and ion sorption thermodynamics were investigated. Ion diffusion coefficients in CR61 were predicted using a recently developed framework based on the Manning diffusion model, which accounts for electrostatic effects on ion diffusion, and the Mackie-Meares model, which accounts for tortuosity/area effects on ion diffusion in swollen polymers. Model predictions agree reasonably well with the experimental data, especially at high external salt concentrations. The observed decrease in salt diffusion coefficients in CR61 relative to values in aqueous solutions is due primarily to tortuosity/area effects. Salt diffusion coefficients correlate with co-ion hydrated radius and membrane water content. Finally, salt permeability coefficients predicted using the solution-diffusion model exhibit good agreement with experimental data, especially at external salt concentration greater than 0.1 mol/L.

42 ENGINEERING↗
DOE OSTI2026

ROADRUNNER uranium nitride MiniFuel: Experimental design, fabrication and pre-irradiation baseline characterization for accelerated burnup testing

Uranium nitride (UN) is a promising fuel candidate for advanced reactor systems owing to its high uranium density and thermal conductivity; however, its qualification remains constrained by the scarcity of well-controlled irradiation performance data. Here, to address this limitation, the ROADRUNNER (Research On ADvancing the peRformance of UraNium Nitrides in Extreme enviRonments) campaign employs the MiniFuel platform in the High Flux Isotope Reactor (HFIR) to enable accelerated burnup irradiation testing under tightly controlled and largely isothermal conditions. This paper presents the experimental design, fuel fabrication, and pre-irradiation baseline characterization of the ROADRUNNER UN MiniFuel campaign. Thirty-six UN minidisc specimens were fabricated with systematically varied as-fabricated density (86–96% of theoretical density), carbon impurity content (961–5240 ppm), oxygen content (≤ ∼2000 ppm), and grain size (2.5–24 μm). The irradiation matrix spans nominal fuel temperatures of 873 K, 1173 K, and 1473 K and target burnups of 3.75%, 6.0%, and 7.5% fissions per initial metal atom (FIMA). Neutronic and thermal analyses were performed to define specimen-specific burnup accumulation and temperature histories, establishing the boundary conditions for subsequent in-pile behavior. Comprehensive pre-irradiation characterization—including dimensional metrology, density verification, impurity analysis, X-ray diffraction, Raman spectroscopy, scanning electron microscopy, X-ray computed tomography, and confocal profilometry—provides a detailed baseline for post-irradiation examination. Pre-irradiation data were further used to generate predictive estimates of fission gas release and swelling using existing empirical correlations. This quantitative comparison reveals substantial inter-model divergence at intermediate and elevated temperatures that exceeds propagated input uncertainties, highlighting structural gaps in the historical irradiation database. The ROADRUNNER irradiation campaign is currently underway in HFIR, with initial firs cycle completed in late 2025 and remaining targets scheduled through 2027. The experimental design and baseline dataset presented here establish the framework needed to interpret forthcoming post-irradiation measurements and to provide discriminating data for the validation and refinement of physics-based UN fuel performance models.

Lopes, Denise Adorno [Oak Ridge National Laborator↗
DOE OSTI2021

Kaon–proton strong interaction at low relative momentum via femtoscopy in Pb–Pb collisions at the LHC

In quantum scattering processes between two particles, aspects characterizing the strong and Coulomb forces can be observed in kinematic distributions of the particle pairs. The sensitivity to the interaction potential reaches a maximum at low relative momentum and vanishing distance between the two particles. Ultrarelativistic heavy-ion collisions at the LHC provide an abundant source of many hadron species and can be employed as a measurement method of scattering parameters that is complementary to scattering experiments. This study confirms that momentum correlations of particles produced in Pb–Pb collisions at the LHC provide an accurate measurement of kaon–proton scattering parameters at low relative momentum, allowing precise access to the K – p $\rightarrow$ K – p process. This work also validates the femtoscopic measurement in ultrarelativistic heavy-ion collisions as an alternative to scattering experiments and a complementary tool to the study of exotic atoms with comparable precision. In this work, the first femtoscopic measurement of momentum correlations of K – p (K + $\bar{p}$) and K + p (K – $\bar{p}$) pairs in Pb–Pb collisions at centre-of-mass energy per nucleon pair of $\sqrt{s_{NN}}$ =5.02 TeV registered by the ALICE experiment is reported. The components of the K – p complex scattering length are extracted and found to be $\mathfrak{R}$f 0 = -0.91±0.03(stat)$^{+0.17}_{-0.03}$(syst) and $\mathfrak{T}$f 0 =0.92±0.05(stat)$^{+0.12}_{-0.33}$(syst). The results are compared with chiral effective field theory predictions as well as with existing data from dedicated scattering and exotic kaonic atom experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Detailed modeling of hydrogen release and particle shrinkage during pyrolysis of inhomogeneous wood

Hydrogen release during pyrolysis of woody biomass is studied considering anisotropicity and inhomogeneity of wood structure. A new anisotropic shrinkage model is proposed based on the decomposition of main wood constituents, i.e., cellulose, hemicellulose, and lignin. The new shrinkage model can predict the temporal evolution of the wood structure, and the differences between axial and radial shrinkage during pyrolysis. The model agrees very well with several experimental data from the literature. Based on particle temperature during conversion, the pyrolysis is partitioned into four stages, and the hydrogen release and H 2 formation from each stage are investigated. Stage (IV) of pyrolysis, from 1000 to 1273 K, is found to be efficient for H 2 production owing to the production of considerable mass of H with a minimal amount of tar species. Furthermore, the char quality is found to be different at the end of stages (II), (III), and (IV), where around 67.7, 80.5, and 93.4% wt. of solid residue is made of carbon, respectively. The model is also used to explain how the heating rate affects the temperature distribution inside the particle and how it shifts the peak of hydrogen release. Finally, the pyrolysis of two inhomogeneous wood samples — a beech twig with bark and a beech dowel with growth rings — are investigated. The bark can affect the pyrolysis rate, products, and flow pattern inside the particle. The growth rings do not have a considerable effect on the pyrolysis rate and products, but they have a significant impact on the flow pattern. This has an important implication for char conversion studies where the internal surface area and porosity field distribution have a significant effect on the gasification and oxidation rates.

09 BIOMASS FUELS↗
DOE OSTI2022

Discrepancy quantification between experimental and simulated data of CO 2 adsorption isotherm using hierarchical Bayesian estimation

Here, to quantitatively analyze the inconsistencies commonly observed between experimental and simulated adsorption isotherms, parameter estimation of adsorption isotherm models was conducted by hierarchical Bayesian estimation with parameter uncertainties being quantified as probability distributions. The estimation method was implemented using Markov Chain Monte Carlo (MCMC) to analyze multiple data sets obtained from different sources, including a publicly available database. To describe the discrepancies of experimental and simulated adsorption data, the simulation data was set as the reference to which experimental measurements were compared. We applied the proposed approach to analyze CO 2 adsorption isotherms that are measured and simulated on zeolite 13X and MIL-101(Cr). In these case studies, the discrepancy of CO 2 adsorption isotherm was successfully quantified between experimental measurements and predictions given by molecular simulations using Grand Canonical Monte Carlo (GCMC), where uncertainties were quantified as probability distributions. Furthermore, experimental data sets that agree well with the GCMC simulation have been identified, providing insights into experimental and measurement methods as well as choosing the right assumptions in the molecular simulation.

42 ENGINEERING↗