Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “empirical”

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 343 records · Page 19

Inadequacy of current approaches for characterizing membrane transport properties at high salinities

Cost optimal design of osmotic membrane processes requires an accurate estimate of membrane transport parameters across their full operational range. However, standard approaches for estimating these parameters rely on empirical methods, the accuracy of which remains unquantified as a function of temperature, salinity, and measurement error. Herein, we present a systematic accuracy analysis of previously developed methods for estimation of membrane transport properties in reverse osmosis, high-pressure reverse osmosis, forward osmosis, pressure retarded reverse osmosis, and osmotically assisted reverse osmosis. We use a Monte Carlo approach to sample the full range of feasible membrane water permeabilities, salt permeabilities, structural parameters, and operating conditions for these processes. These material and process parameters are then incorporated into a physical transport model for each process. Our analysis shows that the statistical uncertainty of current empirical methods for estimating membrane parameters increases by 5 times from low-salinity to high-salinity conditions. Here, the result of this work demonstrates that empirical methods are inadequate for precisely estimating membrane transport properties at high salinity and highlight a critical need for the development of statistically validated higher accuracy methods.

42 ENGINEERING↗

Evaluation of neutron dosimetry capabilities with the MC-15 portable multiplicity counter

This work proposes a preliminary neutron dose rate estimation method for a neutron multiplicity detector through measurement- and simulation-based analyses. Uncharacterized neutron-emitting sources may be encountered in situations such as nuclear emergency response, safeguards, and treaty verification. These circumstances may present irradiation risk to personnel conducting field assay, search, and characterization measurements. It is therefore of interest to provide a field neutron dosimetry capability with the existing neutron multiplicity counting (NMC) capabilities. To date, no commercially-available neutron detection systems are capable of both accurate NMC and real-time neutron dosimetry. This work will focus on estimating dose rate using input from a single fielded NMC called the MC-15. The energy-dependent neutron detection efficiency response of the MC-15 was quantified in monoenergetic neutron simulations and evaluated in response to two neutron-emitting sources and to a polyethylene-moderated source. The results were compared to existing neutron dosimeters and established the proof of concept for further investigation of the MC-15 for dose estimation. Measurement results were also replicated in simulations; additional simulations were then conducted to expand upon the limited empirical data. The initial empirical results provided a conversion factor appropriate for use when measuring 252 Cf neutrons that is independent of polyethylene shielding presence and thickness. The simulated data sets were then used to evaluate a fit equation allowing estimation of the neutron dose rate for less restricted geometric configurations and dependent only on the distance between the source and the detector. Additionally, the energy dependence of the efficiency response indicates that further empirical evaluations could provide energy-dependent conversion factors for broader neutron dosimetry capabilities with a wider range of neutron-emitting sources.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Failure investigation of nuclear grade POCO graphite target in high energy neutrino physics through numerical simulation

In many high energy physics neutrino beamlines, targets are made up of isotropic POCO graphite grade for production of neutrinos for high energy particle physics research and are bombarded with highly energetic pulsed proton beam. The pulsed proton beam with a small beam size creates thermal stress waves as well as radiation damage such as displacement damage, void formation, swelling, and gas formation amongst others. As a result, after a few years of operation one of the longest operated targets appears to have undergone bulk swelling and fractured along the beam centerline. A complex interaction of dynamic loading due to beam, material degradation due to irradiation are thought to have caused such a failure. Here, we present a numerical simulation to capture the combined effects of swelling and the dynamic effect of thermal stress waves due to beam heating, to explain the failure of target. An empirical formula has been developed to take into account swelling as a function of temperature and proton fluence and has been implemented in a commercial finite element code. Extensive X-ray diffraction (XRD) were performed on graphite samples from the failed target to understand the lattice parameter changes due to irradiation and build a valid empirical model. Simulation results show that swelling plays a major role in elevating the stress state in the material exceeding the failure strength while the pulsed beam heating also introduces fatigue loading that would explain brittle fracture emanating from inside of the material and propagating outward. Such empirical formulations will help in identifying the threshold values of critical parameters to extend the service life of future multi-megawatt neutrino targets and targets in other accelerator environments.

43 PARTICLE ACCELERATORS↗

Validation of the single-event method for low-energy electron transport via stopping power calculations with $\mathrm{MCNP}$

Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Experimental and data-driven characterization of window-induced air leakage in residential buildings

Windows contributes up to 40% of envelope heat losses and around 9% of total building energy consumption due to air leakage. In the U.S., 48 million homes still use single-pane windows. Although the U.S. has an estimated 1.4 billion windows in its building stock and about 24 million windows are installed annually, only around 29 million individual window replacements (∼2%) occur each year. To address this gap, this study generates empirical evidence by (1) evaluating the contribution of windows to whole-building air leakage in 20 residential buildings using blower door tests before and after window replacement and (2) assessing whether building and window characteristics influence the measured change. Most simulation studies assume that replacing windows not only lowers the U-factor but also reduces air leakage by 10–20%. However, this assumption lacks empirical validation, highlighting the need for experimental analysis of air leakage specifically associated with windows. Using blower door tests in accordance with ASTM E779–19, the results indicated an average reduction in air infiltration of 6.1% within the range of 0.5–19.30% across all buildings and no significant correlations were found between air leakage improvements and any building/window characteristics. This research aims to help homeowners, and energy modelers to provide empirical data on importance of upgrading to more energy-efficient windows, supporting energy-efficient building standards.

Air leakage↗

Assessing the Accuracy of Property Model Predictions for Cost Optimization of Desalination Technologies

Accurate modeling of seawater thermophysical and thermodynamic properties is critical for optimizing desalination processes. This study compares three seawater property models, a Reaktoro multicomponent model, the thermophysical seawater properties library from the Massachusetts Institute of Technology, and a simplified sodium chloride model, in the context of levelized cost of water (LCOW) minimization for reverse osmosis (RO) and mechanical vapor compression systems. Process simulations and cost optimizations reveal that although all three models yield comparable LCOW and specific energy consumption (SEC) estimates under baseline conditions, deviations among their predictions increase with salinity. Relative differences in LCOW and SEC reach up to 6% and 8%, respectively. RO results show greater variability due to differences in osmotic pressure predictions, which affect pressure constraints at high recoveries. Computational performance varies substantially; specifically, Reaktoro simulations are up to 28 times slower than empirical models due to their detailed equilibrium calculations. These results suggest that empirical models offer acceptable accuracy for routine desalination process design, while Reaktoro provides advantages in scenarios requiring detailed speciation, such as scaling or pH adjustment studies. These findings underscore the importance of selecting appropriate property models based on the modeling objective of desalination applications and motivate future work integrating thermodynamic rigor with empirical efficiency.

Physical and chemical properties↗

Relative Controls of Vapor Pressure Deficit and Soil Water Stress on Canopy Conductance in Global Simulations by an Earth System Model

Climate change impacts drought and vapor pressure deficit (VPD) that critically influence land-atmosphere water vapor and CO 2 exchange. There is a lack of model-based global evaluation of the relative roles of the two types of water stresses in limiting canopy conductance (G c ). Using the Energy Exascale Earth System Model (E3SM), we conducted four land-atmosphere coupled simulations for the historical period with and without a plant hydraulics scheme (PHS) and under present-day and quadrupling of CO 2 concentrations. By fitting an empirical model of G c to the E3SM simulations, we evaluated the sensitivity of G c to water stresses caused by VPD and soil water (SW) in the growing season. The empirical model is based on a multiplicative algorithm that adjusts G c according to VPD, soil water saturation, and their interaction. Our results show positive trends of VPD over land during the historical period, with large differences between simulations with and without PHS. Despite the differences, a majority of land pixels that experience both SW and VPD stresses are in semi-arid and tropical savanna ecosystems. G c is more sensitive to VPD than SW in less than 20% of land pixels where the empirical model provides a good fit for the simulated G c . Ignoring SW and VPD interaction may underestimate the relative role of VPD on G c by as much as 28%. Including PHS accentuates SW stress limitation on G c by simulating wetter soil and higher ET that lower VPD. Elevated CO 2 reduces G c by increasing water use efficiency and ameliorates SW stress on G c .

54 ENVIRONMENTAL SCIENCES↗

Can a simple topological-constraints-based model predict the initial dissolution rate of borosilicate and aluminosilicate glasses?

Tuning glass composition to obtain targeted properties generally relies on empirical approaches. However, a deep understanding of the physical and chemical mechanisms linking glass composition to its structure and properties would enable developing reliable predictive models. Indeed, although empirical models are usually able to interpolate composition–property relationships within a given compositional envelope, they often fail at extrapolating predictions far from their training domain. Here, as an alternative route to empirical models, we show that a structural descriptor based on the number of topological constraints per atom can be used to predict the initial dissolution rate of aluminosilicate and borosilicate glasses after being parameterized on different families of glasses (specific series of borosilicate glasses). Sixteen glasses belonging to these families were studied and their initial dissolution rates were determined at 90 °C and pH 90 °C = 9, covering rates spanning over 5 orders of magnitude. The model based on topological constraints was trained based on seven select borosilicate glasses (R 2 = 0.997) and used to predict the dissolution rate of nine additional borosilicate and aluminosilicate glasses. We show that, provided that corrections are made for high alkali content glasses that dissolve incongruently (preferential release of Na), the model gives reasonable predictions, even far from its training domain.

36 MATERIALS SCIENCE↗

A Novel Relative Permeability Model for Gas and Water Flow in Hydrate-Bearing Sediments With Laboratory and Field-Scale Application

In a producing gas hydrate reservoir the effective porosity available for fluid flow constantly changes with dissociation of gas hydrate. Therefore, accurate prediction of relative permeability using legacy models (e.g. Brooks-Corey (B-C), van Genuchten, etc.) that were developed for conventional oil and gas reservoirs would require empirical parameters to be calibrated at various S h over its range of variation, but such calibrations are precluded because of lack of experimental relative permeability data. This study proposes a new relative permeability model for gas hydrate-bearing media that is a function of maximum capillary pressure, capillary entry pressure, pore size distribution index, residual saturations, hydrate saturation, and four other constants. The three novel features of the proposed model are: (i) requires fitting its six empirical parameters only once using experimental data from any single S h , and the same set of empirical parameters predict relative permeability at all S h , (ii) includes the effect of capillarity, and (iii) includes the effect of pore-size distribution. From practical standpoint, the model can be used to simulate multiphase flow in gas hydrate-bearing sediments where the proposed relative permeability can account for the evolving hydrate saturation. The proposed model is implemented in a numerical simulator and the wall time required to perform simulations using the proposed model is shown to be similar to the time it takes to run same simulations with the B-C model. The proposed model is a step forward towards achieving the goal of physically accurate modeling of multiphase flow for gas hydrate-bearing sediments that accounts for the effect of gas hydrate saturation change on relative permeability.

54 ENVIRONMENTAL SCIENCES↗

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↗

Dust and gas absorption in the high mass X-ray binary IGR J16318-4848

With an absorption column density on the order of 10 24 cm -2 , IGR J16318-4848 is one of the most extreme cases of a highly obscured high mass X-ray binary. In addition to the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. We aim at a self-consistent description of the continuum absorption, the absorption edges, and the fluorescence lines to constrain the properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescence emission from dust grains. We used XMM-Newton and NuSTAR spectra to first empirically constrain the incident continuum and fluorescence lines. Next we used XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In the third step we modeled the absorption and fluorescence in a dusty olivine absorber and employed both a simple analytical model for the fluorescence line emission and a Monte Carlo simulation of radiative transfer that generates line fluxes, which are very close to the observational data. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is able to produce the observed continuum absorption and edges. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318-4848. In particular, for low optical depth of individual grains, which is the case for typical volume densities and grain size distribution models, the dust will contribute significantly to the fluorescence emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Anomalies in the topology of the temperature fluctuations in the cosmic microwave background: An analysis of the NPIPE and FFP10 data releases

We present a topological analysis of the temperature fluctuation maps from the Planck 2020 Data Release 4 NPIPE dataset and the Planck 2018 Data Release 3 FFP10 dataset. We performed a multiscale analysis in terms of the homology characteristics of the maps, invoking relative homology to account for the analysis in the presence of masks. We performed our analysis for a range of smoothing scales spanning sub- and super-horizon scales corresponding to a full width at half maximum (FWHM) of 5',10',20',40',80',160',320', and 640', and employed simulations based on the standard model for comparison, which assumes the initial fluctuation field to be an isotropic and homogeneous Gaussian random field. Examining the behavior of topological components, represented by the 0D homology group, we find the observations to be approximately 2σ or less deviant from the simulations for all resolutions and scales for the NPIPE dataset. For the FFP10 dataset, we detect a 2.96σ deviation between the observations and simulations at N = 128, FWHM = 80'. For the topological loops, represented by the first homology group, the simulations and observations are consistent within 2σ for most resolutions and scales for both the datasets. However, for the NPIPE dataset, we observe a high deviation between the observation and simulations in the number of loops at FWHM = 320', but at a low dimensionless threshold ν = –2.5. Under a Gaussian assumption, this would amount to a deviation of ~4σ. However, the distribution in this bin is manifestly non-Gaussian and does not obey Poisson statistics either. In the absence of a true theoretical understanding, we simply note that the significance is higher than what may be resolved by 600 simulations, yielding an empirical p-value of at most 0.0016. Specifically in this case, our tests indicate that the numbers arise from a statistically stable regime, despite being based on small numbers. For the FFP10 dataset, the differences are not as strong as for the NPIPE dataset, indicating a 2.77σ deviation at this resolution and threshold. The Euler characteristic, which is the alternating sum of the ranks of relative homology groups, reflects the deviations in the components and loops. To assess the significance of combined levels for a given scale, we employed the empirical and theoretical versions of the χ 2 test as well as the nonparametric Tukey depth test. Although all statistics exhibit a stable distribution, we favor the empirical version of the χ 2 test in the final interpretation, as it indicates the most conservative differences. For the NPIPE dataset, we find that the components and loops differ at more than 95%, but agree within the 99% confidence level with respect to the base model at N = 32, FWHM = 320'. The Euler characteristic at this resolution displays a per mil deviation. In contrast, the FFP10 dataset shows that the observations are consistent with the base model within the 95% confidence level, at this and smaller scales. This is consistent with the observations of the Planck analysis pipeline via Minkowski functionals. For the largest smoothing scale, N = 16, FWHM = 640', both datasets exhibit an anomalous behavior of the loops, where FFP10 data exhibit a deviation that is larger by an order of magnitude than that of the NPIPE dataset. In contrast, the values for the topological components and the Euler characteristic agree between observations and model to within a confidence level of 99%. However, for the largest scales, the statistics are based on low numbers and may have to be regarded with caution. Even though both datasets exhibit mild to significant discrepancies, they also exhibit contrasting behaviors at various instances. Therefore, we do not find it feasible to convincingly accept or reject the null hypothesis. Disregarding the large-scale anomalies that persist at similar scales in WMAP and Planck, observations of the cosmic microwave background are largely consistent with the standard cosmological model within 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine Learned Hückel Theory: Interfacing Physics and Deep Neural Networks

The Hückel Hamiltonian is an incredibly simple tight-binding model known for its ability to capture qualitative physics phenomena arising from electron interactions in molecules and materials. Part of its simplicity arises from using only two types of empirically fit physics-motivated parameters: the first describes the orbital energies on each atom and the second describes electronic interactions and bonding between atoms. By replacing these empirical parameters with machine-learned dynamic values, we vastly increase the accuracy of the extended Hückel model. The dynamic values are generated with a deep neural network, which is trained to reproduce orbital energies and densities derived from density functional theory. The resulting model retains interpretability, while the deep neural network parameterization is smooth and accurate and reproduces insightful features of the original empirical parameterization. Altogether, this work shows the promise of utilizing machine learning to formulate simple, accurate, and dynamically parameterized physics models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The difference between molecules and materials: Reassessing the role of exact conditions in density functional theory

Exact conditions have long been used to guide the construction of density functional approximations. However, hundreds of empirical-based approximations tailored for chemistry are in use, of which many neglect these conditions in their design. We analyze well-known conditions and revive several obscure ones. Two crucial distinctions are drawn: that between necessary and sufficient conditions and that between all electronic densities and the subset of realistic Coulombic ground states. Simple search algorithms find that many empirical approximations satisfy many exact conditions for realistic densities and non-empirical approximations satisfy even more conditions than those enforced in their construction. The role of exact conditions in developing approximations is revisited.

Chemistry↗

Large-database cross-verification and validation of tokamak transport models using baselines for comparison

State-of-the-art 1D transport solvers ASTRA and TRANSP are verified, then validated across a large database of semi-randomly selected, time-dependent DIII-D discharges. Various empirical models are provided as baselines to contextualize the validation figures of merit using statistical hypothesis tests. For predicting plasma temperature profiles, no statistically significant advantage is found for the ASTRA and TRANSP simulators over a baseline empirical (two-parameter) model. For predicting stored energy, a significant advantage is found for the simulators over a baseline empirical model based on confinement time scaling. Uncertainty in the results due to diagnostic and profile fitting uncertainties is approximated and determined to be insignificant due in part to the large quantity of discharges employed in the study. Advantages are discussed for validation methodologies like this one that employ (1) large databases and (2) baselines for comparison that are specific to the intended use-case of the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dislocation theory of steady and transient creep of crystalline solids: Predictions for olivine

In applications critical to the geological, materials, and engineering sciences, deformation occurs at strain rates too small to be accessible experimentally. Instead, extrapolations of empirical relationships are used, leading to epistemic uncertainties in predictions. To address these problems, we construct a theory of the fundamental processes affecting dislocations: storage and recovery. We then validate our theory for olivine deformation. This model explains the empirical relationships among strain rate, applied stress, and dislocation density in disparate laboratory regimes. It predicts the previously unexplained dependence of dislocation density on applied stress in olivine. The predictions of our model for Earth conditions differ from extrapolated empirical relationships. For example, it predicts rapid, transient deformation in the upper mantle, consistent with recent measurements of postseismic creep.

36 MATERIALS SCIENCE↗

A Python Tool for Aqueous Plutonium Nitrate Density Law Input Preprocessing in MCNP6

Here, this work develops a predictive density tool in Python, named Plutonium Nitrate Solutions (PuNS), to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer method and an empirical method were implemented into the PuNS tool to generate atom densities for use in MCNP6 material cards. These material cards are directly prepared into an MCNP6 input text file and are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, and plutonium isotope weight percentages. The PuNS tool is validated and verified against the International Criticality Safety Benchmark Evaluation Project Handbook experiments and is observed to predict densities within a root mean square error of 0.89% for the Pitzer method and 1.82% for the empirical method. These errors in density lead to up to 1569 pcm difference in MCNP6 calculated k eff for the Pitzer method and up to a 1751 pcm difference for the empirical method when compared to experimental benchmarks. Simultaneous work is also being performed at Los Alamos National Laboratory and the University of New Mexico to create a similar tool for plutonium chloride solutions, named Plutonium Chloride Solution, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also facilitate an improved understanding of solution systems and a potential relaxation in the conservatism of current aqueous plutonium processing criticality safety limits.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Development of new baseline models for U.S. medium office buildings based on commercial buildings energy consumption survey data

Building energy estimation for the building sector under various scenarios are needed for building energy regulation and policy making. This often starts with representative baselines (either empirical baseline or modeled baseline). Commercial Buildings Energy Consumption Survey (CBECS) data is a widely used empirical baseline for U.S. commercial buildings, but none of the existing baseline model are developed to represent the CBECS data. This paper aims to develop new baseline models for the U.S. medium office buildings, which can produce modeled baselines consistent with the CBECS data. Here, we introduced the methodology to create baseline models and the criteria to evaluate the performance of baseline models. The methodology consists of three phases: (1) identification of model inputs, (2) model calibration, and (3) model validation with uncertainty analysis. The evaluation index is the coefficient of variation of the root-mean-square deviation (CV(RMSD)) of site energy use intensities (EUIs) between the modeled baseline and empirical baseline. Then 30 new baseline models for two vintages (pre- and post-1980) and 15 climate zones were created. The evaluation shows that the CV(RMSD) is lower than 0.05 for the modeled baselines produced by the new baseline models. As a comparison, the CV(RMSD) is higher than 0.1 for the existing modeled baselines generated by DOE Commercial Reference Building Models. Further analysis shows that the new baseline models are able to capture the uncertainties of the representative features of existing buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗