Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic thermal performance”

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 289 records · Page 16

FY24 Progress Report on Viscosity and Thermal Conductivity Measurements of Nuclear Industry Relevant Chloride Salts: An Experimental and Computational Study

As presented in this report, experimental and computational techniques were performed to assess the viscosity and thermal conductivity of key alkali and actinide chloride mixtures for molten salt reactor developers. These mixtures were pure LiCl, NaCl-KCl, LiCl-NaCl, LiCl-KCl, LiCl-NaCl-KCl, and NaCl-UCl 3 . Experimental measurements of viscosity were performed with a rolling ball viscometer, whereas experimental measurements of thermal conductivity were performed with a variable gap apparatus. Additional benchmarking work was performed using both property measurement systems to prepare for x-ray radiography in stainless-steel crucibles for viscosity and to ensure that calibration methods were accurate for thermal conductivity before assessing the NaCl-UCl 3 system. Validation data for the NaCl-UCl 3 in literature are minimal. Details on the calibration methods, salt measurement processes, and sources of error and uncertainty are discussed in detail for both property measurements. The computational methods described herein involved ab-initio molecular dynamics (AIMD) calculations using CP2K. The calculations were performed for the LiCl-KCl-NaCl and NaCl-UCl 3 systems. These calculations not only provided thermophysical property estimations for comparison to experimental data, but they also allowed for the determination of diffusion coefficients, coordination numbers, and radial distribution functions to provide insight into ion mobility and local coordination environments, which is linked to macroscopic property trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Uncertainty quantification of a physics-informed model based on sparse identification of a Thermal Energy Distribution System

Integrated energy systems (IES)s are crucial for enhancing the economy and efficiency of power generation sources (e.g., nuclear energy) necessary to unleash American energy dominance. These systems can be integrated with thermal energy storage (TES) and intermittent renewable energies to optimize overall energy use, peak-load regulation, and demand-side responses. However, the stabilization of energy generation, transport, and utilization introduces operational complexities that exceed the challenges of managing each sub-component individually. Currently, though IESs rely on human operators for efficiency and stability, reducing human error risk and enhancing performance through automation is highly desirable. Recent advances at Idaho National Laboratory have demonstrated successful control of the Thermal Energy Distributed System (TEDS). However, the automatic control system depends on a deterministic Sparse Identification of Nonlinear Dynamics with Control (SINDyC) model, which are trained based on simulation data from physics-based simulations. Because of uncertainties in physics-based simulation, SINDyC model results in large discrepancies against experimental data and cannot be reliably used in automatic control. In this paper, we present an innovative approach to address these discrepancies by quantifying uncertainties and developing a more robust model. We first generated trajectories by using first-principles physics codes to encapsulate the experiment. Next, we trained thousands of models by randomly sampling these trajectories. We then collapsed all those models into one probabilistic SINDyC by fitting a multivariate Gaussian distribution onto the resulting coefficient’s distribution. Despite its simplicity, our approach successfully produced 95% confidence intervals that captured the experimental trajectories. It even did so with a higher probability and better U-pooling score across six of the seven relevant quantities of interest (QoIs), as compared to other classical approaches. In conclusion, ongoing research is focusing on generating new experimental trajectories to validate this approach, and on employing Bayesian calibration to refine parametric uncertainties and guide future model development efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Weak Electron–Phonon Coupling and Enhanced Thermoelectric Performance in n-type PbTe–Cu 2 Se via Dynamic Phase Conversion

This study investigates Ga-doped n-type PbTe thermoelectric materials and the dynamic phase conversion process of the second phases via Cu 2 Se alloying. Introducing Cu 2 Se enhances its electrical transport properties while reducing its lattice thermal conductivity (κlat) via weak electron–phonon coupling. Additionally, Cu 2 Te and CuGa(Te/Se) 2 (tetragonal phase) nanocrystals precipitate during the alloying process, resulting in Te vacancies and interstitial Cu in the PbTe matrix. At room temperature, Te vacancies and interstitial Cu atoms serve as n-type dopants, increasing the carrier concentration and electrical conductivity from ≈1.18 × 10 19 cm –3 and ≈1870 S cm –1 to ≈2.26 × 10 19 cm –3 and ≈3029 S cm –1 , respectively. With increasing temperature, the sample exhibits a dynamic change in Cu 2 Te content and the generation of a new phase of CuGa(Te/Se) 2 (cubic phase), strengthening the phonon scattering and obtaining an ultralow k lat . Pb 0.975 Ga 0.025 Te-3%CuSe exhibits a maximum figure of merit of ≈1.63 at 823 K, making it promising for intermediate-temperature device applications.

36 MATERIALS SCIENCE↗

Dynamic System Scaling Application to Accelerated Nuclear Fuel Testing

The development of nuclear fuel and materials requires a continuous effort to investigate the acute and prolonged effects of irradiation, thermal-material stress, chemical change, or other conceivable damage mechanics acquired during normal operation or accident scenarios throughout its lifetime. As in-core fuel property measurement techniques advance to support in the real-time, non-invasive, and enhanced accuracy realm, it is the future of fuel development to pursue to higher degrees of control, predictability of integral test behavior via separate effects test (SET), and shorter test time intervals. The fuel development life cycle from initial concept to commercial licensing is approximated to be 20 years and current literature suggests by optimizing fuel performance codes with SETs, the process could possibly be compressed to 5 to 10 years. Recently, a research group in the Idaho National Laboratory (INL) is testing reduced scale fuel rods and increased power density to accelerate evolution of fuel phenomena in metallic fuels. In support of nuclear fuel rod development, compressing fuel test process, and accelerating fuel phenomena, it was the purpose of this study to investigate nuclear fuel performance phenomena via literature review and effectively scale the initial conditions, boundary conditions, and geometric properties to describe the time-dependent response including to fuel burnup, thermal and mechanical stress, transmutations and inter-diffusion, and other relevant observed phenomena. The study was based on BISON simulations of historic EBR-II metallic fuel experiments and Dynamical System Scaling (DSS) method are utilized to assess effects of scaling activity including fuel phenomena acceleration and calculations of time-dependent distortions. The research successfully scaled metallic fuel phenomena, accelerated fuel testing, and assessed the distortions for each scaled case derived.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Detailed Simulations of a Three-Stage Supercritical Carbon Dioxide Axial Compressor with a Focus on the Shrouded Stator Cavities Flow

This paper describes the findings of detailed simulations performed to investigate the impact of seal teeth cavity leakage flow on the aerodynamic and thermal performance of a three-stage supercritical CO 2 axial compressor. The study compares a shrouded stator configuration (with cavities) to a cantilevered stator configuration (without cavities) to highlight their differences. High-fidelity computational fluid dynamics simulations were performed using non-linear harmonic (NLH) and mixing plane assumptions, considering various possible rotor/stator interface configurations for mixing plane calculations. The key performance parameters for each case were compared, and the best-performing configuration selected for further analysis. The individual stage performance parameters are also examined and compared between the cantilevered and shrouded configurations. It was observed that in the shrouded case, the leakage flow enters the cavity downstream of the stator trailing edge and gets entrained into the primary flow upstream of the stator, leading to boundary layer changes at the hub and degradation of stator and downstream rotor performance. Vortical flow structures were also observed in the stator wells, which tended to change the flow angles around the region, thereby affecting mixing and velocity distribution, resulting in a slight deterioration of compressor performance. Additionally, the study examines windage heating due to shear work from rotating walls, including the seal teeth surface. The amount of shear work done on the leakage flow and the corresponding rise in fluid temperature were quantified, tabulated, and further compared with a simple analytical model, showing good agreement between them and, hence, validating the numerical approach used.

42 ENGINEERING↗

Analysis of silica fouling on nonwetting surfaces

Ground water sources used as coolant fluids in a variety of thermal systems such as heat exchangers and power plant condensers contain silica particles that accrete on heat transfer surfaces over time leading to reduction in thermal performance, a problem that is particularly exacerbated with temperature. Nonwetting superhydrophobic, lubricant-infused, and a new class of solid-infused surfaces introduced in this work are candidates for fouling mitigation, by virtue of their water repellency, but little is known about fouling of silica on the surfaces, especially under dynamic flow conditions and as a function of temperature. This article presents, for the first time, a systematic study of dynamic flow fouling of silica on nonwetting surfaces vis-à-vis conventional copper surface over a temperature range 20–50°C. The mechanism of silica aggregate formation and its adherence to the different surfaces is elucidated by scanning electron microscope (SEM) imaging. Sigmoidal growth model is used to describe the time evolution of fouling thermal resistance and an Arrhenius model is presented for the temperature-dependent increase in the asymptotic fouling resistance on nonwetting and conventional surfaces alike. Lubricant-infused and solid-infused surfaces are shown to reduce fouling thermal resistance by up to 25% and 13%, respectively, compared to a conventional surface, whereas superhydrophobic surfaces lose their non-wettability under flow conditions, leading to an adverse increase in the fouling resistance by up to 13%. Considering the possible lubricant depletion in lubricant-infused surfaces over prolonged exposure to a flowing fluid, solid-infused surfaces present a robust alternative.

42 ENGINEERING↗

Model Predictive Control for a Grid-interactive Efficient Thermal Storage-integrated Heat Pump System

Building heating and cooling systems can be used to overcome the mismatch between the intermittent supply of renewable power and the fluctuating demand for electricity. A novel underground thermal energy storage integrated with a dual-source heat pump has been proposed to mitigate the mismatch while meeting the thermal demand of buildings efficiently. Conventional thermostat control with heuristic rules cannot provide intelligent decisions to maximize the thermal efficiency and flexibility of the proposed system. Advanced control strategies like model predictive control (MPC) have provided a new paradigm for grid-interactive efficient building operation with the advancement of computation and sensing. This study developed an MPC for the proposed system to provide grid service for Demand Side Management and minimize the operating cost of building owners. A control-oriented dynamic model of the proposed system has been developed. Given an objective function and proper constraints, an optimization problem is formulated to determine the optimal control strategy of the system. Dynamic Programming is adopted to solve the optimization problem. A rule-based control (RBC) is also developed to achieve similar goals. Short-term simulations are conducted to compare the system performance resulting from the two controls. The simulation results indicate that the MPC performs more intelligently than the RBC in charging thermal energy storage and selecting heat pump sources by taking advantage of the predicted cooling demands of the building and the performance of the integrated system. As a result, the MPC could save energy and reduce operating costs compared with the RBC. A case study shows that, for a 3-day operation, the MPC saves 36.9% energy and reduces 38.5% operating cost compared with the RBC.

Shi, Liang↗

Nitrogen: A promising doping strategy for high-performance ovonic threshold switching selectors

The Ovonic Threshold Switching (OTS) selector serves as an essential component in the development of three-dimensional high-density memory integration technology. Nevertheless, the state-of-the-art high-performance OTS materials usually contain toxic elements such as arsenic (As), posing significant risks to both environmental and human health. Nitrogen (N), which belongs to the same group as arsenic (As), has emerged as a highly promising alternative for As doping. However, the underlying mechanisms that govern N-based OTS materials have not yet been extensively investigated. In this study, we delve into the effects of N doping on the structural, bonding, and electronic properties of amorphous GeSe (a-GeNSe) by ab initio molecular dynamics simulations to bridge the knowledge gap. Our findings indicate that upon N doping in a-GeSe, the formation of robust Ge-N bonds, along with N-centered tetrahedral and triangular structures, resulting in the sluggish atomic movement that enhances the thermal stability and endurance of a-GeNSe. The OTS characteristics are significantly influenced by the material’s electronic band structure, and thus the relatively slow performance drift can be attributed to the stabilization of mid-gap states, a result of N doping which effectively slows down the aging process of chalcogenide glass. Moreover, the increased mobility gap in a-GeNSe raises the threshold voltage (V th ), making it more compatible with commercially available phase-change memory materials. Furthermore, our findings reveal the extensive impact of the N element on a typical OTS material and offer valuable perspectives for alternative doping strategies that could potentially supplant As practices.

36 MATERIALS SCIENCE↗

Materials Selection Principles for Designing Electro‐Thermal Neurons

Artificial neurons exhibiting volatile threshold switching and action potential‐like oscillations are crucial for brain‐inspired computing. While Complimentary Metal‐Oxide‐Semiconductor (CMOS)‐based strategies require hundreds of transistors to simulate each neuron, neuronal oscillations arise spontaneously in individual electro‐thermal devices due to nonlinearities like the Mott transition in VO 2 . Despite improved understanding of the physics, quantitative connections between neuronal performance and material properties remain under‐explored, preventing predictive neuron design and rational materials selection. In this work, a physics‐aware forward design methodology is developed for interrogating a wide palette of materials with properties varying by orders of magnitude, and their performance (high frequency, high dynamical reconfigurability and low power) under external circuit and device geometry constraints is assessed. The space of viable materials is identified to be much larger than previously recognized, with candidates from a range of materials classes, including Ge, GaP and MoS 2 . CMOS‐compatible performance (such as 100 GHz oscillating frequencies) can be achieved with CMOS‐compatible node sizes (≈10 nm). Finally, combinations of material properties yielding desired neuronal performance under uncertain design constraints are considered. This work solidifies forward design principles for electro‐thermal neuron devices, a necessary pre‐condition for inverse design from desired neuronal performance to required materials properties.

compact model↗

An Enhanced Computational Framework for Energy Storage Participation in Transmission Planning with Electricity Market Participation: Dynamic Market Participation Restriction Estimation

In the previous study and report, we developed a computational framework to explore the technical and economic viability of the dual use of pumped storage hydropower (PSH) for both transmission assets and market participation. The rationale for the dual-use application is to improve the resource utilization rate and reduce the cost recovery burden. A PSH-based transmission solution generally has a large capacity, which corresponds to a higher investment cost compared to a traditional line solution. Market participation may allow the PSH asset to offset part of its high investment cost. A critical question is whether it is economically and technically viable to allow dual use of a PSH-based transmission project without jeopardizing transmission system reliability. And a more specific follow-up question is how to improve market participation performance while maintaining transmission reliability. In this study, we extend the computational framework by updating the method to determine transmission service obligation requirements dynamically based on the possibility of the need for transmission service in operation. Specifically, we correlate the need for transmission service with the probability of thermal overloading on a monitored transmission line, which is essentially based on the probability distribution of power flow and its probability of exceeding the line rating based on system dispatch and operational uncertainties. In addition, the need depends on the level of risk that a system operator is willing to take. The higher the risk tolerance, the lower the requirement that the system operator would apply to the PSH project as a transmission asset. The requirements are updated every day before the PSH project determines its market participation strategy. With the updated transmission service requirements, it is expected that a PSH project can serve as a transmission asset when the transmission system is at risk and really needs it, while it can participate in a market with sufficient information and certain restrictions. In this way, the transmission service quality and market participation performance can both be improved.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Synthetic Electricity Market Data Generation and HERON Use Case Setup of Advanced Nuclear Reactors Coupled with Thermal Energy Storage Systems

This study evaluates and optimizes advanced nuclear reactors coupled with thermal energy storage (TES) systems in an Integrated Energy System (IES) architecture to enable advanced nuclear power plants (A NPP) to participate in multi-commodity markets, thus enhancing their economic competitiveness. Nuclear-TES coupling scenarios studied herein are designed attenuate the nuclear heat dynamics and defer energy delivery to a later time, enabling the nuclear reactor to continue operating at or near steady-state design conditions as usual while also enabling flexible generation. Three A-NPPs, namely, an advanced light-water reactor (A LWR), a high temperature gas-cooled reactor (HTGR) and a liquid-metal fast reactor (LMFR) were selected as the initial use cases for demonstrating the technoeconomic of thermally balanced energy storage coupling design for thermal power extraction. Each of the reactor technologies were evaluated in two different electricity markets. Stochastic optimization approach was adopted which included the evaluation of price signals from the Pennsylvania-New Jersey-Maryland (PJM) market, and Electric Reliability Council of Texas (ERCOT), using an autoregressive moving average (ARMA) model. Risk Analysis Virtual Environment (RAVEN) tool and its dispatch optimization plugin, the Holistic Energy Resource Optimization Network (HERON), were used to perform dispatch and capacity optimization, using the price data provided by the ARMA models. The results from the Nuclear-TES use cases will be used to design and characterize dynamic integrated system behavior and feedback.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Infrared triggered dwell and active cooling thermal control effects on microstructural uniformity in DED

Directed Energy Deposition (DED) offers rapid large scale fabrication, but difficulty in delivering consistent microstructures and properties hinders the use of DED fabricated components in safety or performance critical applications. Variability stems from the complex thermal cycles generated by the toolpath used to print the required geometry. Several practical methods have become established in DED to regulate overheating, such as active cooling of the baseplate structure or the use of an infrared camera to inject interlayer pauses to ensure the top layer of the component cools to a set temperature, which have been shown to affect microstructure. However, no critical assessment has been performed as to how effective these controls are in promoting microstructural uniformity in the context of complex layer timing commonly generated by non-prismatic geometries. Here we show how controls influence the thermal field, phase transformations, and dynamic annealing of a low-temperature transformation steel using infrared imaging and operando neutron diffraction. Counterintuitively, common thermal homogenization process controls can reduce microstructural uniformity because these approaches stabilize peak temperature while overlooking temperatures near the solid-state phase transformation fronts. Instead, the cyclic reheating induces spatially-variant dynamically annealed regions which can be modulated via control parameters. We show that these controls have spatially linked effects centimeters away from the active weld, which implies that microstructure control must co-optimize thermal input across many subsequent layers. In conclusion, our results demonstrate the pressing need for higher order controls that integrate predictive elements of simulation data to stabilize printed properties for future qualification of DED components.

Directed energy deposition↗

Online Dynamic Mode Decomposition Based System Identification of Multi-Zone Building HVAC Systems

Many works have recently been conducted to reduce the electricity consumption of smart buildings and allow them to support various grid services. Most of these works require accurate system models for the various appliances in the building including heating, ventilation, and air conditioning (HVAC) units. In this paper, we investigate a recursive data-driven system identification strategy to construct the thermal model for a time-varying building with a multi-zone HVAC unit. The online dynamic mode decomposition (DMD)-based strategy is employed to identify the multi-zone thermal building dynamics, where a simple information update (rank-1) is selected to avoid computational complexity. The DMD-based identification strategy is validated using a real gymnasium building equipped with a 4-zone HVAC unit, and its performance is compared with that of the traditional nuclear-norm subspace identification (N2SID) strategy.

Wu, Tumin [University of Tennessee, Knoxville (UTK↗

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are the most important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a direct comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. Here, we used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a millimetre VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects, compared to the addition of thermal noise only. Using synthetic data based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M 87, we performed case studies to assess the image quality that can be obtained with the current and future EHT array for different weather conditions. Our synthetic observations show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of our GRMHD source models can be recovered robustly with the EHT2017 array after performing calibration steps, which include fringe fitting, a priori amplitude and network calibration, and self-calibration. With the planned addition of new stations to the EHT array in the coming years, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

79 ASTRONOMY AND ASTROPHYSICS↗

Flow and thermal modelling of the argon volume in the DarkSide-20k TPC

The DarkSide-20k dark matter experiment, currently under construction at LNGS, features a dual-phase time projection chamber (TPC) with a ∼ 50 t argon target from an underground well. At this scale, it is crucial to optimise the argon flow pattern for efficient target purification and for fast distribution of internal gaseous calibration sources with lifetimes of the order of hours. To this end, we have performed computational fluid dynamics simulations and heat transfer calculations. The residence time distribution shows that the detector is well-mixed on time-scales of the turnover time (∼ 40 d). Notably, simulations show that despite a two-order-of-magnitude difference between the turnover time and the half-life of 83m Kr of 1.83 h, source atoms have the highest probability to reach the centre of the TPC 13 min after their injection, allowing for a homogeneous distribution before undergoing radioactive decay. We further analyse the thermal aspects of dual-phase operation and define the requirements for the formation of a stable gas pocket on top of the liquid. We find a best-estimate value for the heat transfer rate at the liquid-gas interface of 62 W with an upper limit of 144 W and a minimum gas pocket inlet temperature of 89 K to avoid condensation on the acrylic anode. This study also informs the placement of liquid inlets and outlets in the TPC. The presented techniques are widely applicable to other large-scale, noble-liquid detectors.

47 OTHER INSTRUMENTATION↗

Electronic Structure of the Au Nanoparticle-TiO 2 Heterojunction: Influence of Nanoparticle Size, Shape, Oxygen Vacancies, and Temperature

The electronic structure of the gold nanoparticle-titanium dioxide (AuNP-TiO 2 ) heterojunction plays a critical role for charge transfer and recombination dynamics that underpin its photocatalytic function. However, building a representative model to capture the key physics remains a significant challenge. Here, we investigate the influence of the AuNP size and shape, as well as oxygen vacancy (V O ) defects at the anatase-phase TiO 2 (101) surface and the temperature of the heterojunction, on its interfacial electronic properties. Using density functional theory (DFT), we compare the closed-shell Au 20 and open-shell Au 19 clusters interfaced with pristine and V O defect TiO 2 surfaces. We find that the presence of a V O defect transforms pure TiO 2 from a p-type to an n-type semiconductor, reversing the interfacial band bending from downward to upward. For the heterosystem, density of states (DOS) analysis shows that V O minimally affects the open-shell Au 19 -TiO 2 system, but it significantly alters the closed-shell Au 20 -TiO 2 system, converting it from p-type to n-type at the Γ-point at 0 K. Furthermore, ab initio molecular dynamics (AIMD) simulations at 300 K reveal significant thermal fluctuations in AuNP positions relative to the TiO 2 surface. These fluctuations result in dynamic variations in the gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) in all systems studied: Au 19 NP-Pristine TiO 2 , Au 20 NP-Pristine TiO 2 , Au 19 NP-V O defect TiO 2 , and Au 20 NP-V O defect TiO 2 . Complementary Bader charge analysis performed at both 0 K and finite-temperature AIMD snapshots supports the emergence of an upward band bending and the formation of a Schottky barrier at the V O containing heterojunctions. Notably, we find that at finite temperature, an Au atom can dynamically passivate the V O , leading to a pronounced widening of the HOMO-LUMO gap in the AuNP-V O defect TiO 2 heterojunctions. Our computational findings underscore the pivotal role of the V O defect and thermal effects in modulating interfacial band alignment, electronic states, and HOMO-LUMO gaps, providing insights for designing Au-TiO 2 heterojunctions with tailored electronic properties.

Mora Perez, Carlos [Lawrence Berkeley National Lab↗

Force field for calcium sulfate minerals to predict structural, hydration, and interfacial properties

Calcium sulfates such as anhydrite, hemihydrate, and gypsum find widespread use in building materials, implants, and tissue healing. We introduce a simple and compatible atomistic force field for all calcium sulfate phases that reproduces a wide range of experimental data including lattice parameters, surface, hydration, mechanical, and thermal properties in 1% to 5% accuracy relative to experiments. The performance is several times better than prior force fields and DFT methods, which lead to errors in structures and energies up to 100%. We explain (hkl) cleavage energies, the dynamics of (hkl) water interfaces, and new insights into molecular origins of crystal-facet specific hydration and solubility. Impressive agreement of computed and experimentally measured hydration energies is shown. The models add to the Interface force field (IFF) and are compatible with multiple force fields (CHARMM, AMBER, GROMOS, CVFF, PCFF, OPLS-AA) for property predictions of sulfate-containing materials from atoms to the large nanometer scale.

36 MATERIALS SCIENCE↗