Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “water model”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

On Modeling Air/Space-Borne Radar Returns in the Melting Layer

The bright band is the enhanced radar echo associated with the melting of hydrometeors in stratiform rain where the melting process usually occurs below 0 C isotherm over a distance of about 500m. To simulate this radar signature, a scattering model of melting snow is proposed in which the fractional water content is prescribed as a function of the radius of a spherical mixed- phase particle consisting of air, ice and water. The model is based on the observation that melting starts at the surface of the particle and then gradually develops towards the center. To compute the scattering parameters of a non-uniform melting particle, the particle is modeled as a sphere represented by a collection of 64(exp 3) cubic cells of identical size where the probability of water at any cell is prescribed as a function of the radius. The internal field of the particle, used for deriving the effective dielectric constant, is computed by the Conjugate Gradient and Fast Fourier Transform (CGFFT) numerical methods. To make computations of the scattering parameters more efficient, a multi-layer stratified-sphere scattering model is introduced after demonstrating that the scattering parameters of the non-uniformly melting particle can be accurately reproduced by the stratified sphere. In conjunction with a melting layer model that describes the melting fractions and fall velocities of hydrometeors as a function of the distance from the 0 C isotherm, the stratified-sphere model is used to simulate the radar bright band profiles. These simulated profiles are shown to compare well with measurements from the Precipitation Radar (PR) aboard the Tropical Rainfall Measuring Mission (TRMM) satellite and a dual-wavelength airborne radar. The results suggest that the proposed model of a melting snow particle may be useful in studying the characteristics of the bright-band in particular and mixed- phase hydrometeors in general.

Liao, Liang↗
DOE OSTI2021

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree–Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Furthermore, our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

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

Evaluating Variable‐Resolution CESM Over China and Western United States for Use in Water‐Energy Nexus and Impacts Modeling

Abstract Climate models are critical tools to study earth system processes and are often further downscaled to provide refined‐resolution data sets for regional water‐energy impact analyses. In this study, the Variable‐Resolution Community Earth System Model (VR‐CESM), a new dynamical downscaling method, is employed to generate 1/8° simulations for China and the western United States over 1970–2006. Precipitation, temperature, snowpack, radiation, and wind speed are compared with two 1° global climate models, regional climate model ensemble, and reference data sets to evaluate their fidelity. Simulated precipitation skill is dependent upon the reference data sets used, a slight (<+10%) and spatially consistent wet bias in China during summer and the western United States during winter. Precipitation bias in VR‐CESM is the smallest among models and demonstrates the added value provided by resolution refinement. However, temperature biases are generally greater than precipitation due to more complicated interactions among processes. VR‐CESM simulates +2°C warm bias in the low‐elevation regions in China but −2°C cold bias in mountains. However, DJF cold bias (−2.23°C to−3.37°C) and JJA warm bias (+1.13°C to +1.87°C) exist in the western United States. The modeling skills of precipitation and temperature at various topographies are similar in the two evaluated regions. The cold bias was likely affected by downward shortwave radiation deficit (−15%) and slow bias of surface wind speeds (−16% to −29%). Despite the wet and cold bias, snow water equivalent and snow cover are underestimated in mountains. We conclude that VR‐CESM outperforms other models evaluated and provides the best estimation of downscaled climate model data for water‐energy studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Anomalies and Local Structure of Liquid Water from Boiling to the Supercooled Regime as Predicted by the Many-Body MB-pol Model

For the past 50 years, researchers have sought molecular models that can accurately reproduce water’s microscopic structure and thermophysical properties across broad ranges of its complex phase diagram. Herein, molecular dynamics simulations with the many-body MB-pol model are performed to monitor the thermodynamic response functions and local structure of liquid water from the boiling point down to deeply supercooled temperatures at ambient pressure. The isothermal compressibility and isobaric heat capacity show maxima near 223 K, in excellent agreement with recent experiments, and the liquid density exhibits a minimum at ~208 K. A local tetrahedral arrangement, where each water molecule accepts and donates two hydrogen bonds, is found to be the most probable hydrogen-bonding topology at all temperatures. Further, this work suggests that MB-pol may provide predictive capability for studies of liquid water’s physical properties across broad ranges of thermodynamic states, including the so-called water’s “no man’s land” which is difficult to probe experimentally.

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

Consequences of the failure of equipartition for the p – V behavior of liquid water and the hydration free energy components of a small protein

Earlier we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step δ t ≤ 0.5 fs to ensure equipartition between translational and rotational modes. Here we extend that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examine integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small δ t is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, and the dielectric constant, a measure of dipole moment fluctuations, are also seen to be sensitive to δ t . Using the mean volume estimated from the NpT simulation, we examine the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to δ t . In going from δ t = 2 fs to δ t = 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein.

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

Dataset for manuscript "Consequences of the failure of equipartition for the p-V behavior of liquid water and the hydration free energy components of a small protein"

Previously, we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step dt that is less than or equal to 0.5 fs to ensure equipartition between translational and rotational modes. We extended that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examined integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small time-step, dt, is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, is seen to be sensitive to dt. Using the mean volume estimated from the NpT simulation, we examined the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to dt. In going from a time-step of 2 fs to a time-step of 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein. The data-set contains the simulation metadata and log files that support the claims noted above.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Water and vapor transport in algal‐fungal lichen: Modeling constrained by laboratory experiments, an application for Flavoparmelia caperata

ABSTRACT Algal‐fungal symbionts share water, nutrients, and gases via an architecture unique to lichens. Because lichen activity is controlled by moisture dynamics, understanding water transport is prerequisite to understand their fundamental biology. We propose a model of water distributions within foliose lichens governed by laws of fluid motion. Our model differentiates between water stored in symbionts, on extracellular surfaces, and in distinct morphological layers. We parameterize our model with hydraulic properties inverted from laboratory measurements of Flavoparmelia caperata and validate for wetting and drying. We ask: (1) Where is the bottleneck to water transport? (2) How do hydration and dehydration dynamics differ? and (3) What causes these differences? Resistance to vapor flow is concentrated at thallus surfaces and acts as the bottleneck for equilibrium, while internal resistances are small. The model captures hysteresis in hydration and desiccation, which are shown to be controlled by nonlinearities in hydraulic capacitance. Muting existing nonlinearities slowed drying and accelerated wetting, while exaggerating nonlinearities accelerated drying and slowed wetting. The hydraulic nonlinearity of F. caperata is considerable, which may reflect its preference for humid and stable environments. The model establishes the physical foundation for future investigations of transport of water, gas, and sugar between symbionts.

Potkay, Aaron↗
NASA NTRS2002

Ocean Raman Scattering in Satellite Backscatter UV Measurements

Ocean Raman scattering significantly contributes to the filling-in of solar Fraunhofer lines measured by satellite backscatter ultraviolet (buy) instruments in the cloudless atmosphere over clear ocean waters. A model accounting for this effect in buy measurements is developed and compared with observations from the Global Ozone Monitoring Experiment (GONE). The model extends existing models for ocean Raman scattering to the UV spectral range. Ocean Raman scattering radiance is propagated through the atmosphere using a concept of the Lambert equivalent reflectively and an accurate radiative transfer model for Rayleigh scattering. The model and observations can be used to evaluate laboratory measurements of pure water absorption in the UV. The good agreement between model and observations suggests that buy instruments may be useful for estimating chlorophyll content.

Vasilkov, Alexander P.↗
NASA NTRS2005

Lunar and Planetary Science XXXVI, Part 6

Contents include the following: A Model for Multiple Populations of Presolar Diamonds. Characterization of Martian North Polar Geologic Units Using Mars Odyssey THEMIS Data. Effect of Flow on the Internal Structure of the Martian North Polar Layered Deposits. Elemental Abundance Distributions in Basalt Clays and Meteorites: Is It a Biosignature? Early Results on the Saturn System from the Composite Infrared Spectrometer. NanoSIMS D/H Imaging of Isotopically Primitive Interplanetary Dust Particles. Presolar (Circumstellar and Interstellar) Phases in Renazzo: The Effects of Parent Body Processing. Catastrophic Disruption of Hydrated Targets: Implications for the Hydrated Asteroids and for the Production of Interplanetary Dust Particles. Chemical and Mineralogical Analyses of Particles from the Stratospheric Collections Coinciding with the 2002 Leonid Storm and the 2003 Comet Grigg-Skjellerup Trail Passage. An Analysis of the Solvus in the CaS-MnS System. ESA s SMART-1 Mission at the Moon: First Results, Status and Next Steps. Europa Analog Ice-splitting Measurements and Experiments with Ice-Hunveyor on the Frozen Balaton-Lake, Hungary. Chromium on Eros: Further Evidence of Ordinary Chondrite Composition. Dust Devil Tracks on Mars: Observation and Analysis from Orbit and the Surface. Spatial Variation of Methane and Other Trace Gases Detected on Mars: Interpretation with a General Circulation Model. Mars Water Ice and Carbon Dioxide Seasonal Polar Caps: GCM Modeling and Comparison with Mars Express Omega Observations. Component Separation of OMEGA Spectra with ICA. Clathrate Formation in the Near-Surface Environment of Titan. Space Weathering: A Proposed Laboratory Approach to Explaining the Sulfur Depletion on Eros. Sample Collection from Small Airless Bodies: Examination of Temperature Constraints for the TGIP. Sample Collector for the Hera Near-Earth Asteroid Sample Return Mission. A Rugged Miniature Mass-Spectrometer for Measuring Aqueous Geochemistry on Mars. Martian and Lunar Pyroxene Microstructures Studied by Single-Crystal X-Ray Diffraction.

Source record↗
NASA NTRS1982

CONSERVB: A numerical method to compute soil water content and temperature profiles under a bare surface

A comprehensive, yet fairly simple model of water disposition in a bare soil profile under the sequential impact of rain storms and other atmospheric influences, as they occur from hour to hour is presented. This model is intended mostly to support field studies of soil moisture dynamics by our current team, to serve as a background for the microwave measurements, and, eventually, to serve as a point of departure for soil moisture predictions for estimates based in part upon airborne measurements. The main distinction of the current model is that it accounts not only for the moisture flow in the soil-atmosphere system, but also for the energy flow and, hence, calculates system temperatures. Also, the model is of a dynamic nature, capable of supporting any required degree of resolution in time and space. Much critical testing of the sample is needed before the complexities of the hydrology of a vegetated surface can be related meaningfully to microwave observations.

Vanbavel, C. H. M.↗
DOE OSTI2024

A Novel Standard Gibbs Energy of Formation Model for High-Enthalpy Water Systems

This work provides an advanced standard molar Gibbs energy of formation model, informed by molecular statistical thermodynamics (MST), and validated by mineral solubility and ion association reactions. This new model aligns with experimental data within uncertainties and highlights the role of specific MST interactions around the critical point. This work will provide a reliable tool for the optimization of industrial processes operating at the edges of the supercritical domain.

Hall, Derek↗
DOE OSTI2025

Extracted Water Induces Concentration Fluctuations in Model Ternary Liquid–Liquid Extraction System

Organic phase aggregation and phase splitting are two important (and connected) phenomena in liquid–liquid extraction of metal ions. Previously, we demonstrated how structure in binary mixtures of extractant and diluent originates from concentration fluctuations associated with the so-called “third-phase formation” liquid–liquid phase instability. For more complex systems, we sought to understand how extracted aqueous components affect these concentration fluctuations. Here, in this study, we use small-angle X-ray scattering and molecular dynamics simulations to investigate organic phase structural changes upon contact with water in a model ternary system of water/alkane/extractant. Critical exponents were determined from scattering with temperature variation and were consistent with mean field values. Then, by controlling the water content in the organic phase through direct dissolution of fixed quantities at constant temperature, we relate the observed increase in the fluctuation length scale to changes to the spinodal temperature of the third-phase formation phase boundary. We find a roughly linear relationship between water concentration and the spinodal temperature and a similar slope in that relationship for different linear alkane diluents. This suggests that the mechanism by which water extraction impacts the organic phase structure is to enhance concentration fluctuations through the introduction of additional polar extractant-water complexes, which increases the spinodal temperature and thereby decreases the reduced temperature, i.e., the distance to the instability. This explanation for how solutes induce aggregation in liquid–liquid extraction organic phases may extend to more complex systems, explaining the large increases in correlation lengths reported under high acid and metal loading.

Bonnett, Brittany L. [Argonne National Laboratory ↗
NASA NTRS2022

Instantaneous Photosynthetically Available Radiation (IPAR) prediction models based on Neural Network for ocean waters.

Instantaneous photosynthetically available radiation (IPAR) at the ocean surface and its vertical profile below the surface play a critical role in models to calculate net primary productivity of marine phytoplankton. In this work, we report two IPAR prediction models based on neural network (NN) approach, one for open ocean and the other for coastal waters. These models are trained, validated, and tested using a large volume of synthetic datasets for open ocean and coastal waters simulated by a radiative transfer model. Our NN models are designed to predict IPAR under a large range of atmospheric and oceanic conditions. The NN models can compute subsurface IPAR profile very accurately up to euphotic zone depth. The root mean square errors associated with the diffuse attenuation coefficient of IPAR are less than 0.011 𝑚−1 and 0.036 𝑚−1 for open ocean and coastal waters respectively. The performance of the NN models is better than presently available semi analytical models, with significant superiority in coastal waters.

PACE↗
DOE OSTI2021

Summary of CTF Modeling and Numerical Improvements for Boiling Water Reactor Simulation

This report documents geometry and numerical improvements made to CTF for the modeling of boiling water reactor (BWR) geometry and operating conditions. These activities are part of a larger program to extend the Virtual Environment for Reactor Applications (VERA) to better support BWR modeling and simulation. The activities documented in this report added features to CTF, including support for mixed-fuel cores, modeling of the upper plenum, and modeling of the lower tie plate form losses. A review of the spacer grid modeling approach was also performed, and a plan was discussed for future improvement. The parallelization of the model was improved, leading to a roughly 2× improvement in CTF runtime and a 1.6× improvement in total VERA runtime. An in-depth review of the governing equations and their linearization was performed and is documented in this report. Once implemented, this new linearization will allow CTF to take much larger timesteps, leading to more significant reductions in CTF and VERA runtimes.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1992

Distribution of water on Mars: Implications from SNC meteorites

There has been much speculation about the abundance of water and other volatiles on Mars. Attempts to calculate abundances of water on Mars indicate that Mars contains approx. 10-100 m of water. Numerous models have been put forth to determine the amount of water on Mars more closely. Some researchers infer that Chassigny parent magma contained greater than 1.5 percent water by weight and that the Martian mantle contained greater than 1000 parts per million water. This is too much water for a depleted region. Perhaps some of the water in Chassigny was assimilated at shallow depths, either in a crustal magma chamber or by interaction with superficial permafrost. Either is possible and provides an alternative to the dilemma of water-rich to depleted regions.

Jones, J. H.↗
NASA NTRS2004

Formation of Martian Gullies by the Flow of Simultaneously Freezing and Boiling Liquid Water

Geomorphic evidence suggests that recent gullies on Mars were formed by fluvial activity. The Martian gully features are significant because their existence implies the presence of liquid water near the surface on Mars in geologically recent times. Irrespective of the ultimate source of the fluid carving the gullies, we seek to understand the behavior of this fluid after it reaches the Martian surface. We find that, contrary to popular belief, the fluvially-carved Martian gullies require formation conditions such as now occur on Mars, outside of the temperature-pressure stability regime of liquid water. Mars Global Surveyor observations of gully length and our modeling of water stability are consistent with gully formation from the action of pure liquid water that is simultaneously boiling and freezing.

Heldmann, Jennifer L.↗
DOE OSTI2022

Modelling single atom catalysts for water splitting and fuel cells: A tutorial review

Here, in this tutorial review, we report the state-of-the-art of the modeling approaches of Single Atom Catalysts (SAC) for water splitting and fuel cells reactions. The discussion applies for Hydrogen Evolution Reaction (HER), Oxygen Reduction Reaction (OER), Hydrogen Oxidation Reaction (HOR), and Oxygen Reduction Reaction (ORR). The main scope of this work is to underline the relevant aspects of SACs modelling. On the one hand, the review could help computational chemists aiming to start the study of SACs. On the other hand, experimentalists could find the critical analysis of DFT results of interest to understand better the strengths and weaknesses of simulations, and how to interpret computational results. After an introductory section of SACs, we start by briefly presenting the state-of-the-art methodologies. Then, we analyze the critical aspects for a reliable prediction of the electronic properties, and we discuss the robustness of the structural models and ways to validate them. Furthermore, we discuss the main approaches to predict catalytic activity and selectivity, which is the final goal of the computational catalysis. We conclude this review with a critical analysis of the current challenges in the field, and the main limitations of the modeling approaches that are described.

25 ENERGY STORAGE↗