Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Error estimation”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Advanced Finite-Volume Numerics and Source Term Assumptions for Kernel and G-Equation Modelling of Propane/Air Flames

Here G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device. We tested the explicit (sub-cycled) vs. implicit formulation for the standard transport operators (advection, diffusion, compressibility). In addition to the accurate flame swept-volume method for chemistry and species source term, we developed a more accurate estimator for the burnt/unburnt split cell composition. Then, we developed a signed-distance-function (SDF) based method which provides a more stable reinitialization of the level-set field at every time-step. We found that simplifying assumptions common to several G-Equation implementations, for straightforward terms such as compressibility and advection, lead to large errors in predicting the propagation of even laminar flames, with deviations up to ~300% in simulated vs. formulated flame speed. Conversely, the enhanced numerics enabled through the SDF field reinitialization and improved chemistry source term improve simulation stability and smooth flame propagation even with significantly larger solver time-steps.

42 ENGINEERING↗
DOE OSTI2022

Large-scale scenarios of electric vehicle charging with a data-driven model of control

Transportation electrification is forecast to bring millions of new electric vehicles to roads worldwide this decade. Planning to support those vehicles depends on detailed scenarios of their electricity demand in both uncontrolled and controlled or smart charging scenarios. In this work, we present a novel modeling approach to enable rapid generation of demand estimates that represent the impact of controlled charging for large-scale scenarios with millions of individual drivers. To model the effect of load modulation control on aggregate charging profiles, we propose a novel machine learning approach that replaces traditional optimization approaches. We demonstrate its performance modeling workplace charging control under a range of electricity rate schedules, achieving small errors (2.5%–4.5%) while accelerating computations by more than 4000 times. To generate the uncontrolled charging demand for scenarios with residential, workplace, and public charging we use statistical representations of a large data set of real charging sessions. We demonstrate the methodology by generating diverse sets of scenarios for California's charging demand in 2030 which consider multiple charging segments and controls, each run locally in under 50 s. We further demonstrate support for rate design by modeling the large-scale impact of a new, custom rate schedule for workplace charging.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Simultaneous determination of electrical conductivity and thickness of nonmagnetic metallic foils using two-sided multifrequency eddy current techniques

Here, in this paper, we describe a simple method for performing multifrequency eddy current characterization of free-standing uniform-thickness metallic foils using a forked inductive coil arrangement. The method involves measuring the mutual inductance between two coils when a foil is present between the coils, and when it is not present; the ratio of these mutual inductances is compared with an analytical solution, and foil conductivity, thickness, and sheet resistance are simultaneously estimated using numerical inversion and least-squares fitting. This method was used to characterize 34 non-ferrous metallic samples with thicknesses between 50 and 640 μm and with conductivities between 0.8 × 10 7 and 5.8 × 10 7 S/m. The estimated thicknesses from eddy current characterization agreed well with those measured using confocal optical techniques; the two approaches agreed to within 1 μm for samples that were thinner than 200 μm, and to within 0.5% for samples that had a thickness of 200 μm or greater. The estimated conductivities from eddy current characterization were in close agreement with expected values, given knowledge of the materials used. A particular strength of this approach is that the instrumentation needed is broadly available in research and development laboratories and the associated fixturing is easy to manufacture and assemble. A calibration procedure is described that can be used to reduce errors from geometric uncertainties. This calibration requires a sample that has only a known conductivity or thickness; both do not need to be known. The method described herein is likely extensible to conductivities and thickness well outside the ranges measured as part of this work.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Towards real-time monitoring: data assimilated time-lapse full waveform inversion for seismic velocity and uncertainty estimation

SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.

58 GEOSCIENCES↗
DOE OSTI2023

zPerf: A Statistical Gray-Box Approach to Performance Modeling and Extrapolation for Scientific Lossy Compression

With the scaling up of simulation-based scientific discovery on high-performance computing systems, the disparity between compute and I/O has increased, forcing domain scientists to save only a small amount of simulation data to persistent storage. This can result in the loss of essential physics fields that are needed for data analysis. While error-bounded lossy compression has made tremendous progress in bridging the gap between compute and I/O, the lack of understanding of compression performance remains a key hurdle to its wide adoption. Here, in this work, we present zPerf, a statistical gray-box performance modeling approach for scientific lossy compression. Our contributions are threefold: 1) We develop zPerf to estimate the performance of lossy compression techniques, based on in-depth understanding and statistical modeling for data features and core compression metrics; 2) We demonstrate the in-detailed implementation of zPerf using two case studies, where we derive the performance modeling for SZ and ZFP, two leading lossy compressors; 3) We evaluate the effectiveness of zPerf on real-world datasets across various domains. Based on the evaluation, we demonstrate the efficacy of the zPerf performance model; 4) We further discuss three case studies where zPerf is applied to extrapolate the compression ratio of SZ and ZFP with alternative encoding schemes as well as ZFP with an alternative transform scheme. Through the case studies, we demonstrate the potential of zPerf for exploring the design space of lossy compression, which has hardly been studied in the literature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Simulation Based Inference with Domain Adaptation for Strong Gravitational Lensing

Simulation based inference leverages machine learning to carry out Bayesian inference for systems with intractable likelihoods. However, transitioning a network trained on simulated data to real data runs the risk of encountering domain shift, leading to performance losses. We attempt to implement domain adaptation into the sbi neural posterior estimation framework using the Maximum Mean Discrepancy as an additional network loss. We test two network architectures and use masked autoregressive flow for density estimation. We test the network on a set of 400,000 simulated strong gravitational lensing images generated using deeplenstronomy. The source domain is defined as low noise whereas the target domain has a noise profile sampled from experimentally derived DES survey conditions. We find that, while DA does lead to performance improvements, they are marginal at ~6% less inference error. We also find a similar marginal improvement in uncertainty calibration at around 8%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Synthetic method of analogues for emerging infectious disease forecasting

The Method of Analogues (MOA) has gained popularity in the past decade for infectious disease forecasting due to its non-parametric nature. In MOA, the local behavior observed in a time series is matched to the local behaviors of several historical time series. The known values that directly follow the historical time series that best match the observed time series are used to calculate a forecast. This non-parametric approach leverages historical trends to produce forecasts without extensive parameterization, making it highly adaptable. However, MOA is limited in scenarios where historical data is sparse. This limitation was particularly evident during the early stages of the COVID-19 pandemic, where the emerging global epidemic had little-to-no historical data. In this work, we propose a new method inspired by MOA, called the Synthetic Method of Analogues (sMOA). sMOA replaces historical disease data with a library of synthetic data that describe a broad range of possible disease trends. This model circumvents the need to estimate explicit parameter values by instead matching segments of ongoing time series data to a comprehensive library of synthetically generated segments of time series data. We demonstrate that sMOA has competitive performance with state-of-the-art infectious disease forecasting models, out-performing 78% of models from the COVID-19 Forecasting Hub in terms of averaged Mean Absolute Error and 76% of models from the COVID-19 Forecasting Hub in terms of averaged Weighted Interval Score. Additionally, we introduce a novel uncertainty quantification methodology designed for the onset of emerging epidemics. Developing versatile approaches that do not rely on historical data and can maintain high accuracy in the face of novel pandemics is critical for enhancing public health decision-making and strengthening preparedness for future outbreaks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Introducing the GeoRePORT Resource Size Tool: Reporting on Geothermal Resource Size Estimations Using the Geothermal Resource Portfolio Optimization and Reporting Technique (GeoRePORT): Preprint

The Geothermal Resource Portfolio Optimization and Reporting Technique (GeoRePORT) was developed with funding from the U.S. Department of Energy Geothermal Technologies Office to assist in identifying and pursuing long-term investment strategies through the development of a resource reporting protocol. The assessment protocols used in GeoRePORT allow for comparison of project attributes across locations and geological settings to understand the feasibility of geothermal development. This work introduces the Resource Size Tool, a new feature within the GeoRePORT package that compiles two independent methods for estimating geothermal resource size in terms of energy capacity in MW. Energy production potential for twenty-three case studies was estimated with the Resource Size Tool in order to 1) generate a reasonable range of resource size estimates for a particular geothermal field; 2) illustrate the advantages and limitations of each methodology (such as data input requirements, estimate accuracy and precision, and the appropriate circumstances of use); and 3) test the ability of the resource size tool to provide useful and accurate information for geothermal stakeholders. The tool employs two methods widely used in the geothermal industry: (1) USGS Volumetric and (2) Power Density. Results from our case studies show general overlap between these two methods in terms of resource size estimates; however, they also reveal key differences between the two approaches that should be considered when using such estimates to drive development. First, the two methods rely on different input parameters and therefore one method may be more appropriate and/or accurate for a given project than the other. Second, the Power Density method was found to generate wider ranges of resource size predictions, more consistently aligning with actual power production of the field but with larger scales of error; whereas the USGS Volumetric method predicts narrower ranges but tends to overestimate when compared to current MW production. Future work will refine variables used in the methods with input data from other sections of GeoRePORT and modify uncertainty levels based on the particular datasets used for a given project.

geological↗
DOE OSTI2021

A Novel Machine Learning Algorithm for Planetary Boundary Layer Height Estimation Using AERI Measurement Data

Accurately determining the height of the planetary boundary layer (PBL) is important since it can affect the climate, weather, and air quality. Ground-based infrared hyperspectral remote sensing is an effective way to obtain this parameter. Compared with radiosonde measurements, its temporal resolution is much higher. In this study, a method to retrieve the PBL height (PBLH) from the ground-based infrared hyperspectral radiance data is proposed based on machine learning. In this method, the channels that are sensitive to temperature and humidity profiles are selected as the feature vectors, and the PBLHs derived from radiosonde are taken as the true values. The support vector machine (SVM) is applied to train and test the data set, and the parameters are optimized in the process. The data set collected at the Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) from 2012 to 2015 is analyzed. The instruments used in this letter include Atmospheric Emitted Radiance Interferometer (AERI), Vaisala CL31 ceilometer, and radiosonde. It shows that the root mean square error (RMSE) between the PBLHs calculated by the proposed method using AERI data and those from radiosonde data can be within 370 m, and the square correlation coefficient (SCC) is greater than 0.7. Compared with the PBLHs derived from the ceilometer, it can be found that the new method is more stable and less affected by clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A machine learning estimator trained on synthetic data for real-time earthquake ground-shaking predictions in Southern California

Abstract After large-magnitude earthquakes, a crucial task for impact assessment is to rapidly and accurately estimate the ground shaking in the affected region. To satisfy real-time constraints, intensity measures are traditionally evaluated with empirical Ground Motion Models that can drastically limit the accuracy of the estimated values. As an alternative, here we present Machine Learning strategies trained on physics-based simulations that require similar evaluation times. We trained and validated the proposed Machine Learning-based Estimator for ground shaking maps with one of the largest existing datasets (<100M simulated seismograms) from CyberShake developed by the Southern California Earthquake Center covering the Los Angeles basin. For a well-tailored synthetic database, our predictions outperform empirical Ground Motion Models provided that the events considered are compatible with the training data. Using the proposed strategy we show significant error reductions not only for synthetic, but also for five real historical earthquakes, relative to empirical Ground Motion Models.

Environmental Sciences & Ecology↗
DOE OSTI2021

Dark Energy Survey Year 3 Results: Galaxy mock catalogs for BAO analysis

The calibration and validation of scientific analysis in simulations is a fundamental tool to ensure unbiased and robust results in observational cosmology. In particular, mock galaxy catalogs are a crucial resource to achieve these goals in the measurement of baryon acoustic oscillation (BAO) in the clustering of galaxies. Here we present a set of 1952 galaxy mock catalogs designed to mimic the Dark Energy Survey Year 3 BAO sample over its full photometric redshift range 0.6 < zphoto < 1.1. The mocks are based upon 488 ICE-COLA fast N-body simulations of full-sky light cones and were created by populating halos with galaxies, using a hybrid halo occupation distribution – halo abundance matching model. This model has ten free parameters, which were determined, for the first time, using an automatic likelihood minimization procedure. We also introduced a novel technique to assign photometric redshift for simulated galaxies, following a two-dimensional probability distribution with VIMOS Public Extragalactic Redshift Survey data. The calibration was designed to match the observed abundance of galaxies as a function of photometric redshift, the distribution of photometric redshift errors, and the clustering amplitude on scales smaller than those used for BAO measurements. An exhaustive analysis was done to ensure that the mocks reproduce the input properties. Finally, mocks were tested by comparing the angular correlation function w(θ), angular power spectrum Cℓ, and projected clustering ξp(r⊥) to theoretical predictions and data. The impact of volume replication in the estimate of the covariance is also investigated. The success in accurately reproducing the photometric redshift uncertainties and the galaxy clustering as a function of redshift render this mock creation pipeline as a benchmark for future analyses of photometric galaxy surveys.Key words: catalogs / large-scale structure of Universe / galaxies: distances and redshifts / Galaxy: halo / methods: numerical

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Kinetic equilibrium reconstructions of plasmas in the MAST database and preparation for reconstruction of the first plasmas in MAST upgrade

Reconstructions of plasma equilibria using magnetic sensors were routine during operation of the Mega Ampere Spherical Tokamak (MAST) device, but reconstructions using kinetic profiles were not. These are necessary for stability and disruption analysis of the MAST database, as well as for operation in the upgrade to the device, MAST-U. The three-dimensional (3D) code VALEN is used to determine eddy currents in the 3D vessel structures for vacuum coil test shots, which are then mapped to effective resistances in the two-dimensional vessel groupings in the EFIT equilibrium reconstruction code to be used in conjunction with nearby loop voltage measurements for estimated currents in the structures during reconstruction. Kinetic equilibrium reconstructions with EFIT, using all available magnetic sensors as well as Thomson scattering measurements of electron temperature and density, charge exchange recombination spectroscopy measurements of ion temperature, and internal magnetic field pitch angle measurements from a motional Stark effect (MSE) diagnostic are performed for a large database of MAST discharges. Excellent convergence errors are obtained for the portions of the discharges where the stored energy was not too low, and it is found that reconstructions performed with temperature and density measurements but without MSE data usually already match the pitch angle measurements well. A database of 275 kinetic equilibria is used to test the ideal MHD stability calculation capability for MAST. In conclusion, the necessary changes to conducting structure in VALEN, and diagnostic setup in EFIT have been completed for the upgrade from MAST to MAST-U, enabling kinetic reconstructions to commence from the first plasma discharges of the upgraded device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Velocity-dependent J-factors for Milky Way dwarf spheroidal analogues in cosmological simulations

We study the impact of the dark matter velocity distribution modelling on signals from velocity-dependent dark matter annihilation in Milky Way dwarf spheroidal galaxies. Using the high resolution APOSTLE simulations, we identify analogues corresponding to Milky Way dwarf spheroidal galaxies, and from these directly determine the dark matter pair-wise relative velocity distribution, and compare to best-fitting Maxwell-Boltzmann distribution models. For three velocity-dependent annihilation models, p-wave, d-wave, and the Sommerfeld model, we quantify the errors introduced when using the Maxwell-Boltzmann parameterization. We extract a simple power-law relation between the maximum circular velocity of the dwarf spheroidal analogue and the peak speed of the Maxwell-Boltzmann distribution. We show that this relation can be used to accurately calculate the dark matter relative velocity distribution, and find that it allows us to estimate the dark matter annihilation signal without the need to directly calculate the relative velocity distribution for each galaxy. The scatter in the J-factors calculated from the analogues dominates the uncertainty obtained when compared to the J-factor as determined from the observational data for each dwarf spheroidal, with the largest scatter from d-wave models and the smallest from Sommerfeld models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Sensor System and Observer Algorithm Co-Design For Modern Internal Combustion Engine Air Management Based on H2 Optimization

This paper outlines a novel sensor selection and observer design algorithm for linear time-invariant systems with both process and measurement noise based on H 2 optimization to optimize the tradeoff between the observer error and the number of required sensors. The optimization problem is relaxed to a sequence of convex optimization problems that minimize the cost function consisting of the H 2 norm of the observer error and the weighted l 1 norm of the observer gain. An LMI formulation allows for efficient solution via semi-definite programing. The approach is applied here, for the first time, to a turbo-charged spark-ignited engine using exhaust gas circulation to determine the optimal sensor sets for real-time intake manifold burnt gas mass fraction estimation. Simulation with the candidate estimator embedded in a high fidelity engine GT-Power model demonstrates that the optimal sensor sets selected using this algorithm have the best H 2 estimation performance. Sensor redundancy is also analyzed based on the algorithm results. This algorithm is applicable for any type of modern internal combustion engines to reduce system design time and experimental efforts typically required for selecting optimal sensor sets.

Zhang, Xu↗
DOE OSTI2022

Bayesian homodyne and heterodyne tomography

Continuous-variable (CV) photonic states are of increasing interest in quantum information science, bolstered by features such as deterministic resource state generation and error correction via bosonic codes. Data-efficient characterization methods will prove critical in the fine-tuning and maturation of such CV quantum technology. Although Bayesian inference offers appealing properties—including uncertainty quantification and optimality in mean-squared error—Bayesian methods have yet to be demonstrated for the tomography of arbitrary CV states. Here we introduce a complete Bayesian quantum state tomography workflow capable of inferring generic CV states measured by homodyne or heterodyne detection, with no assumption of Gaussianity. As examples, we demonstrate our approach on experimental coherent, thermal, and cat state data, obtaining excellent agreement between our Bayesian estimates and theoretical predictions. Our approach lays the groundwork for Bayesian estimation of highly complex CV quantum states in emerging quantum photonic platforms, such as quantum communications networks and sensors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Evidential Deep Learning: Enhancing Predictive Uncertainty Estimation for Earth System Science Applications

Abstract Robust quantification of predictive uncertainty is a critical addition needed for machine learning applied to weather and climate problems to improve the understanding of what is driving prediction sensitivity. Ensembles of machine learning models provide predictive uncertainty estimates in a conceptually simple way but require multiple models for training and prediction, increasing computational cost and latency. Parametric deep learning can estimate uncertainty with one model by predicting the parameters of a probability distribution but does not account for epistemic uncertainty. Evidential deep learning, a technique that extends parametric deep learning to higher-order distributions, can account for both aleatoric and epistemic uncertainties with one model. This study compares the uncertainty derived from evidential neural networks to that obtained from ensembles. Through applications of the classification of winter precipitation type and regression of surface-layer fluxes, we show evidential deep learning models attaining predictive accuracy rivaling standard methods while robustly quantifying both sources of uncertainty. We evaluate the uncertainty in terms of how well the predictions are calibrated and how well the uncertainty correlates with prediction error. Analyses of uncertainty in the context of the inputs reveal sensitivities to underlying meteorological processes, facilitating interpretation of the models. The conceptual simplicity, interpretability, and computational efficiency of evidential neural networks make them highly extensible, offering a promising approach for reliable and practical uncertainty quantification in Earth system science modeling. To encourage broader adoption of evidential deep learning, we have developed a new Python package, Machine Integration and Learning for Earth Systems (MILES) group Generalized Uncertainty for Earth System Science (GUESS) (MILES-GUESS) ( https://github.com/ai2es/miles-guess ), that enables users to train and evaluate both evidential and ensemble deep learning. Significance Statement This study demonstrates a new technique, evidential deep learning, for robust and computationally efficient uncertainty quantification in modeling the Earth system. The method integrates probabilistic principles into deep neural networks, enabling the estimation of both aleatoric uncertainty from noisy data and epistemic uncertainty from model limitations using a single model. Our analyses reveal how decomposing these uncertainties provides valuable insights into reliability, accuracy, and model shortcomings. We show that the approach can rival standard methods in classification and regression tasks within atmospheric science while offering practical advantages such as computational efficiency. With further advances, evidential networks have the potential to enhance risk assessment and decision-making across meteorology by improving uncertainty quantification, a longstanding challenge. This work establishes a strong foundation and motivation for the broader adoption of evidential learning, where properly quantifying uncertainties is critical yet lacking.

Schreck, John S.↗
DOE OSTI2020

Assessing the blockage effect of wind turbines and wind farms using an analytical vortex model

Wind farm blockage effects are currently neglected in the prediction of wind farm energy yield, typically leading to an overestimation of the production. This work presents a novel method to assess wind farm production, while accounting for blockage effects. We apply a vortex model, based on a cylindrical wake, to assess induction effects. We present variations of the model to account for finite wake length, finite tip-speed ratios, and the proximity to the ground. The results are applied to single rotors in aligned and yawed conditions and to different wind farm layouts. We provide far-field approximations for faster estimates of the velocity field. Further, this article includes a new methodology to couple the induction model to engineering wake models, such as the ones present in the FLOw Redirection and Induction in Steady State (FLORIS). We compare the results to actuator disk simulations for various operating conditions of a single turbine and different wind farm layouts. We found that the mean relative error of the model in the induction zone is typically around 0.2% compared with actuator disk simulations. The computational time of the velocity field using the analytical vortex model is three orders of magnitude less than the one obtained with the actuator disk simulation.

17 WIND ENERGY↗
DOE OSTI2020

Computational fluid dynamics analysis of char conversion in Sandia’s pressurized entrained flow reactor

Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia’s pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.

47 OTHER INSTRUMENTATION↗