Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Carbonyl Identification and Quantification Uncertainties for Oxidative Polymer Degradation

The most revealing indicator for oxidative processes or state of degraded plastics is usually carbonyl formation, a key step in materials degradation as part of the carbon cycle for man-made materials. Hence, the identification and quantification of carbonyl species with infrared spectroscopy have been the method of choice for generations, thanks to their strong absorbance and being an essential intermediate in carbon oxidation pathways. Despite their importance, precise identification and quantification can be challenging and rigorous fully traceable data are surprisingly rare in the existing literature. An overview of the complexity of carbonyl quantification is presented by the screening of reference compounds in solution with transmission and polymer films with ATR IR spectroscopy, and systematic data analyses. Significant variances in existing data and their past use have been recognized. Guidance is offered how better measurements and data reporting could be accomplished. Experimental variances depend on the combination of uncertainty in exact carbonyl species, extinction coefficient, contributions from neighboring convoluting peaks, matrix interaction phenomena and instrumental variations in primary IR spectral acquisition (refractive index and penetration depth for ATR measurements). In addition, diverging sources for relevant extinction coefficients may exist, based on original spectral acquisition. For common polymer degradation challenges, a relative comparison of carbonyl yields for a material is easily accessible, but quantification for other purposes, such as degradation rates and spatially dependent interpretation, requires thorough experimental validation. All variables featured in this overview demonstrate the significant error margins in carbonyl quantification, with exact carbonyl species and extinction coefficients already being major contributors on their own.

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

Results from a wake-steering experiment at a commercial wind plant: investigating the wind speed dependence of wake-steering performance

Wake steering is a wind farm control strategy in which upstream wind turbines are misaligned with the wind to redirect their wakes away from downstream turbines, thereby increasing the net wind plant power production and reducing fatigue loads generated by wake turbulence. In this paper, we present results from a wake-steering experiment at a commercial wind plant involving two wind turbines spaced 3.7 rotor diameters apart. During the 3-month experiment period, we estimate that wake steering reduced wake losses by 5.6% for the wind direction sector investigated. After applying a long-term correction based on the site wind rose, the reduction in wake losses increases to 9.3%. As a function of wind speed, we find large energy improvements near cut-in wind speed, where wake steering can prevent the downstream wind turbine from shutting down. Yet for wind speeds between 6–8 m/s, we observe little change in performance with wake steering. However, wake steering was found to improve energy production significantly for below-rated wind speeds from 8–12 m/s. By measuring the relationship between yaw misalignment and power production using a nacelle lidar, we attribute much of the improvement in wake-steering performance at higher wind speeds to a significant reduction in the power loss of the upstream turbine as wind speed increases. Additionally, we find higher wind direction variability at lower wind speeds, which contributes to poor performance in the 6–8 m/s wind speed bin because of slow yaw controller dynamics. Further, we compare the measured performance of wake steering to predictions using the FLORIS (FLOw Redirection and Induction in Steady State) wind farm control tool coupled with a wind direction variability model. Although the achieved yaw offsets at the upstream wind turbine fall short of the intended yaw offsets, we find that they are predicted well by the wind direction variability model. When incorporating the expected yaw offsets, estimates of the energy improvement from wake steering using FLORIS closely match the experimental results.

17 WIND ENERGY↗
DOE OSTI2026

Corn grown on dredged sediments alters rat behavior in the elevated plus maze via the gut mycobiome

The use of marginal lands for growing food has increased worldwide. Many of these lands are supplemented with soil amendments to enhance their productivity, but they may also introduce contaminants depending on their origin. The use of sediments dredged from commercial waterways as a soil amendment has increased in recent years. However, there is significant concern regarding their potential to negatively affect public health by transferring contaminants, like metals, to the food supply. This study examined the behavior, microbiome and physiology of rats fed corn grown on dredged sediments with rats fed commercial feed corn and related those metrics to kernel metal content. Metal content did not differ between kernels grown on dredged sediments and feed corn for most metals except iron which was greater in commercial feed corn. Although there is some variability by animal response, animals fed corn grown on dredged sediments had a gut microbial community that generally promoted anxiety-like behavior while animals fed feed corn had a gut microbial community that promoted a more typical response. These findings suggest the effect of corn grown on dredged sediments on public health is complex and highlights the need for further investigation.

Rua, Megan [Wright State University, Dayton, OH]↗
DOE OSTI2026

Understanding and mitigating degradation in amine-based sorbents for CO 2 direct air capture

The success of direct air capture (DAC) of CO 2 depends on sorbents that combine high capacity, low energy requirements, and long-term durability. Amine-based sorbents-including solid-supported aminopolymers, grafted amines, and amine-functionalized resins-remain the leading candidates, but their limited lifetimes drive up costs and constrain deployment. In this review, we outline the current understanding of amine-based sorbent degradation with an emphasis on clearly identifying what is known about structure-property-performance relationships, as well as important knowledge gaps. More specifically, we discuss how polymer chemistry, sorbent design variables, and environmental and process conditions contribute to performance loss. In parallel, we outline how advances in spectroscopy, modeling, and accelerated testing are beginning to illuminate chemical and physical degradation mechanisms. Looking forward, we identify future research directions that will be critical for gaining a deeper understanding of degradation, as well as opportunities for developing innovative mitigation strategies for improving the lifetime of amine-based sorbents.

organic↗
DOE OSTI2025

Joint Modeling of Wind Speed and Wind Direction Through a Conditional Approach

Atmospheric near surface wind speed and wind direction play an important role in many applications, ranging from air quality modeling, building design, wind turbine placement to climate change research. It is therefore crucial to accurately estimate the joint probability distribution of wind speed and direction. In this work, we develop a conditional approach to model these two variables, where the joint distribution is decomposed into the product of the marginal distribution of wind direction and the conditional distribution of wind speed given wind direction. To accommodate the circular nature of wind direction, a von Mises mixture model is used; the conditional wind speed distribution is modeled as a directional dependent Weibull distribution via a two-stage estimation procedure, consisting of a directional binned Weibull parameter estimation, followed by a harmonic regression to estimate the dependence of the Weibull parameters on wind direction. A Monte Carlo simulation study indicates that our method outperforms two other approaches in estimation efficiency: one that utilizes periodic spline quantile regression and another that generates data from the commonly used Abe-Ley distribution for cylindrical data. We illustrate our method by using the output from a regional climate model to investigate how the joint distribution of wind speed and direction may change under some future climate scenarios. Our method indicates significant changes in the variation of wind speed with respect to some directions.

17 WIND ENERGY↗
DOE OSTI2023

Assessing sensitivities of climate model weighting to multiple methods, variables, and domains in the south-central United States

Given the increasing use of climate projections and multi-model ensemble weighting for a diverse array of applications, this project assesses the sensitivities of climate model weighting strategies and their resulting ensemble means to multiple components, such as the weighting schemes, climate variables, or spatial domains of interest. The purpose of this study is to assess the sensitivities associated with multi-model weighting strategies. The analysis makes use of global climate models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and their statistically downscaled counterparts created with the localized constructed analogs (LOCA) method. This work focuses on historical and projected future mean precipitation and daily high temperatures of the south-central United States. Results suggest that the model weights and the corresponding weighted model means can be sensitive to the weighting strategy that is applied. For instance, when estimating model weights based on Louisiana precipitation, the weighted projections show a wetter and cooler south-central domain in the future compared to other weighting strategies. Alternatively, for example, when estimating model weights based on New Mexico temperature, the weighted projections show a drier and warmer south-central domain in the future. However, when considering the entire south-central domain in estimating the model weights, the weighted future projections show a compromise in the precipitation and temperature estimates. As for uncertainty, our matrix of results provided a more certain picture of future climate compared to the spread in the original model ensemble. If future impact assessments utilize weighting strategies, then our findings suggest that how the specific weighting strategy is used with climate projections may depend on the needs of an impact assessment or adaptation plan.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Evaluating Atmospheric River Impacts on Energy and Moisture Transport in the Arctic Using Different Detection Algorithms

Atmospheric rivers (ARs) significantly impact the Arctic climate system by enhancing atmospheric heat and moisture transport and altering the local energy budget. Developing AR detection tools (ARDTs) is critical yet challenging. This study evaluates 12 ARDTs in the Arctic to assess their performance in representing atmospheric heat (represented by moist static energy) and moisture transport, as well as surface downward longwave radiation (LWD) and precipitation impacts, spanning 2000 to 2019 using ERA5 reanalysis. We find that AR occurrence frequency in the Arctic varies widely, from less than 1% to over 13%, depending on the ARDT. This variability leads to differences in contributions to poleward atmospheric heat (<1%–33%) and moisture (<1%–49%) transport. The highest AR frequency, and corresponding contributions to atmospheric heat and moisture transport, occurs over the Atlantic sector during non-summer seasons for most ARDTs. This region aligns with the primary poleward moisture pathway and the end of climatological mid-latitude storm tracks, highlighting strong connections between Arctic ARs and mid-latitude cyclones. ARs induce significant LWD anomalies, largest in winter, smallest in summer, and also substantially contribute to the seasonal precipitation. Global ARDTs detect fewer ARs with larger anomalies (>100 W m−2 in higher Arctic), but contribute <1% to seasonal climatological LWD and precipitation. In contrast, polar-specific ARDTs detect higher AR occurrences and account for up to 16% of seasonal LWD and 41% precipitation. This suggests that algorithms emphasizing extreme events with large anomalies do not necessarily indicate a large climate radiative and precipitation impact.

Zhang, Chen [University of Colorado, Boulder, CO (↗
DOE OSTI2025

Predictive dynamic wetting, fluid–structure interaction simulations for braze run-out

Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. Here, the quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of the space weather conditions on the Earth magnetosphere soft X-ray emissivity

The aim of the study is to model and characterize the soft X-ray emissivity on the Earth magnetosphere for different space weather conditions (SWC), providing information to interpret the soft X-ray measurements of the Solar wind Magnetosphere Ionosphere Link Explorer space mission. The MHD code pluto in spherical coordinates is used to perform parametric studies with respect to the solar wind (SW) dynamic pressure (considering density and velocity effects independently) as well as the IMF intensity and orientation, predicting the soft X-ray emissivity for different SWC. The integrated soft X-ray emissivity inside the magnetosheath is calculated as a proxy of the soft X-ray emission dependencies with the SWC independently of the satellite orbit and camera line of sight. The analysis indicates fluctuations of the interplanetary magnetic field (IMF) orientation and magnitude may significantly affect the measured soft X-ray emission although changes in the SW dynamic pressure should be the main source of variability. The southward IMF orientation leads to the configuration with the largest soft X-ray emissivity and northward to the lowest. Strongly distorted magnetospheres explored in configurations showing SW and IMF parameters comparable to the impact of interplanetary coronal mass ejections may show a decrease of the soft X-ray emissivity as the IMF magnitude increases, explained by the strong magnetosphere compression and constriction of the magnetosheath region where the soft X-ray emissivity maximum is located. The simulations also indicate large excursions of the soft X-ray emissivity maximum inside the magnetosheath as the IMF magnitude and SW dynamic pressure fluctuate particularly for radial and ecliptic IMF orientations.

Earth↗
DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A decade of insights: Delving into calendar aging trends and implications

Lithium-ion batteries remain at rest for extended periods and experience calendar aging. Although lithium-ion batteries are expected to perform for over 10 years at room temperature, long-term calendar aging data are seldom reported over such timescales. We present a dataset from 232 commercial cells across eight cell types and five manufacturers that underwent calendar aging across various temperatures and states of charge (SOCs) for up to 13 years. We analyze calendar aging across these conditions by tracking capacity loss and resistance growth as the cells degrade. This dataset is used to validate simple models, primarily the Arrhenius law and the power law, which explain the temperature and storage time on calendar aging. Certain applications of Arrhenius and power law fail to describe the dependence of capacity loss on temperature and resistance growth on storage time. Through this dataset, we demonstrate the complexity of calendar aging and the challenges in reducing trends into phenomenological models.

25 ENERGY STORAGE↗
DOE OSTI2021

The effects of storms and a transient sandy veneer on the interannual planform evolution of a low-relief coastal cliff and shore platform at Sargent Beach, Texas, USA

Coastal cliff erosion is alongshore-variable and episodic, with retreat rates that depend upon sediment as either tools of abrasion or protective cover. However, the feedbacks between coastal cliff planform morphology, retreat rate, and sediment cover are poorly quantified. This study investigates Sargent Beach, Texas, USA, at the annual to interannual scale to explore (1) the relationship between temporal and spatial variability in cliff retreat rate, roughness, and sinuosity and (2) the response of retreat rate and roughness to changes in sand and shell hash cover of the underlying mud substrate as well as the impact of major storms using field measurements of sediment cover, erosion, and aerial images to measure shore platform morphology and retreat. A storm event in 2009 increased the planform roughness and sinuosity of the coastal cliff at Sargent Beach. Following the storm, aerial-image-derived shorelines with annual resolution show a decrease in average alongshore erosion rates from 12 to 4 m yr –1 , coincident with a decrease in shoreline roughness and sinuosity (smoothing). Like the previous storm, a storm event in 2017 increased the planform roughness and sinuosity of the cliff. Over shorter timescales, monthly retreat of the sea cliff occurred only when the platform was sparsely covered with sediment cover on the shore platform, indicating that the tools and cover effects can significantly affect short-term erosion rates. The timescale to return to a smooth shoreline following a storm or roughening event, given a steady-state erosion rate, is approximately 24 years, with the long-term rate suggesting a maximum of ~107 years until Sargent Beach breaches, compromising the Gulf Intracoastal Waterway (GIWW) under current conditions and assuming no future storms or intervention. The observed retreat rate varies, both spatially and temporally, with cliff face morphology, demonstrating the importance of multi-scale measurements and analysis for interpretation of coastal processes and patterns of cliff retreat.

58 GEOSCIENCES↗
DOE OSTI2024

Not even 6 dB: Gaussian quantum illumination in thermal background

Abstract In analyses of target detection with Gaussian state transmitters in a thermal background, the thermal occupation is taken to depend on the target reflectivity in a way which simplifies the analysis of the symmetric quantum hypothesis testing problem. However, this assumption precludes comparison of target detection performance between an arbitrary transmitter and a vacuum state transmitter, i.e. ‘detection without illumination’, which is relevant in a bright thermal background because a target can be detected by its optical shadow or some other perturbation of the background. Using a target-agnostic thermal environment leads to the result that the oft-claimed 6 dB possible reduction in the quantum Chernoff exponent for a two-mode squeezed vacuum transmitter over a coherent state transmitter in high-occupation thermal background is an unachievable limiting value, only occurring in a limit in which the target detection problem is ill-posed. Further analyzing quantum illumination in a target-agnostic thermal environment shows that a weak single-mode squeezed transmitter performs worse than ‘no illumination’, which is explained by the noise-increasing property of reflected low-intensity squeezed light.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Arctic Impact Identification with Less Data Using Variable Relationships: An Exploratory Express LDRD project.

Regional impacts from sea ice loss can be challenging to separate from internal climate variability, potentially requiring thousands of ensemble members. East Asian wintertime cooling has been linked to sea ice loss from present day conditions in the Polar Amplification Model Intercomparison Project with these large ensemble counts. This cooling is theorized to arise from a strengthened Siberian High and East Asian Jet response. The strengthened Siberian High can be detected with one fifth the ensemble members needed for the East Asian wintertime cooling in a single model. We thus hypothesize that leveraging relationships between multiple variables in a conditional pathways-based approach would reduce the number of required ensemble members to conclusively attribute East Asian wintertime cooling to future sea ice concentrations. In all analyzed cases, confidence was increased when evaluating sea ice loss’s responsibility for the joint effects of East Asian cooling, East Asian Jet strengthening, and Siberian High strengthening over just East Asian cooling. However, we were not able to confidently attribute future East Asian wintertime cooling to sea ice loss in a single model. We found that significant intra-ensemble variability within single Earth System Models (ESMs) produced highly uncertain forcing response models upon which attribution results were undermined. We were able to show that ensemble mean seasonally averaged metrics from multiple ESMs greatly improved the accuracy of the forcing response linear models and exposed the necessity of all three steps in the pathway (sea ice area, Siberian High pressure, and East Asian Jet speed) for accurate prediction of East Asian wintertime cooling. Although all three steps were necessary, East Asian wintertime cooling possesses a large dependence on the Siberian High pressure, which weakens the confidence associated with overall strong joint-attribution comparing present day and future scenarios. We believe transitioning the pathway nodes to relative changes between the Siberian High and Aleutian Low as well as between the midlatitude westerlies and subtropical jet in the East Asianj Jet region may be able to produce significant attribution more fully dependent upon all three steps. Ultimately, this research demonstrates the simple extensibility of conditional pathways-based attribution to sea ice loss forcing on the Earth system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Direct sensitivity analysis on the parameterization of crystal plasticity models

Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.

Elementary Effects Test↗
DOE OSTI2025

The influence of laser power modulation on melt pool dynamics in laser powder bed fusion

While the majority of laser powder bed fusion (LPBF) metal additive manufacturing uses a continuous wave (CW) laser heat source, some commercial applications of LPBF additive manufacturing instead involve the modulation of the laser power on tens-of-microsecond timescales as an adjustable process variable. This article reports the use of in situ, high speed x-ray and optical imaging to probe melt pool fluid flow, defect formation, and nearby powder motion during LPBF with both modulated and CW laser heat sources. We observe melt pool dynamics unique to modulated laser melting even at very high duty cycles that are related to fluctuations in vapor depression depth, complex pore formation mechanisms, and changes to denudation physics when compared to CW melting. These behaviors are present in Ti–6Al–4V, 316L stainless steel, and AL1100 alloys but vary slightly as a function of material, indicating a substantial dependence on the viscosity and surface tension of the liquid metal. While high duty cycles produce weld tracks of comparable quality to CW melting, lower duty cycles introduce substantial defect concentrations. At intermediate duty cycles, careful control of modulation parameters can repeatably and precisely yield one pore per laser pulse, suggesting a method for intentionally inserting engineered porosity at specific sites during an LPBF build.

3D printing↗
DOE OSTI2023

A Systems Approach to Increasing Carbon Flux to Seed Oil for Biofuels and Bioproducts Production in Camelina sativa (Final Report)

To combat climate change and alleviate the dependency of the United States on fossil fuels, the transition to biofuel crops has long been proposed as a crucial part of the long-term solution. Camelina sativa has emerged as one of the leading commercially viable options for biofuel and bioproduct production for the U.S. Camelina has the advantages of low agronomic inputs and natural resistance to diverse biotic and abiotic stresses relative to other oilseed crops, and Camelina oil-based blends have been tested and approved as liquid transportation fuels. A major limitation in the widespread adoption of Camelina as a viable industrial oilseed crop is its modest oil yields. This project directly investigated possible paths towards increasing oil, by employing tissue-specific and whole plant systems approaches to identify major regulatory bottlenecks. We developed a new high-throughput method for the identification of multi-gene transformants in polyploid species like Camelina sativa and used this to determine that, contrary to what mathematical models suggested, introduction of the microbial Entner-Doudoroff (ED) pathway into Camelina sativa did not result in significant seed oil increases. We developed flux maps providing accurate and statistically robust models of central carbon metabolism in Camelina, including of cultured Camelina embryos. This knowledge was used in the capacitation of a number of junior researchers in gaining knowledge on approaches used for quantitative flux map analyses. We identified several novel Camelina transcription factor genes regulating fatty acid biosynthesis that showed variable effects on seed oil accumulation in transgenic plants. Among a number of community resources, we also developed a knowledge base web resource (CamRegBase) to integrate Camelina gene regulatory information. Together, research outcomes from this project have contributed to a much better understanding of the regulation and bottlenecks to engineer seed oil production in Camelina.

09 BIOMASS FUELS↗
DOE OSTI2020

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water (Rev. 1)

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗