Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium modeling”

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

Modeling swelling behavior of hydrogels in aqueous organic solvents

Hydrogel, a three-dimensional elastic structure of cross-linked polymer chains, swells or deswells when immersed in aqueous organic solutions. The swelling of hydrogels depends on temperature, solvent composition, and hydrogel structure as the thermodynamic driving forces for hydrogel swelling are due to solvent – hydrogel interactions and expansion of elastic network. Here in this work, polymer Nonrandom Two-Liquid activity coefficient model is used to account for the solvent – hydrogel interactions and a semi-empirical expression for rubber elasticity is used to describe the expansion of elastic network. The model successfully correlates data for both swelling and equilibrium phase compositions for n-isopropyl acrylamide hydrogels in water – ethanol, water – acetone, water – n-butanol, and water – methyl isobutyl ketone mixed solvents.

42 ENGINEERING↗
DOE OSTI2020

EMISApproximateEquilibrium.jl [SWR-19-56]

The Electricity Markets Investment Suite Approximate Equilibrium (EMIS-AE) package is developed at NREL and provides a solution capability to solve generation expansion equilibrium problems in the electricity markets. EMIS-AE is designed to capture the evolution of the electricity generation portfolio resulting from the interactions of heterogeneous investors under different policy and market designs. The investment problem for each generation company (GENCO) is a bi-level problem with the investment decision made in the upper level and market clearing condition in the lower level, which traditionally is represented as a Mathematical Program with Equilibrium Constraint (MPEC). EMIS-AE provides a predictive model to be trained for estimating the system-wide revenues for each technology type across energy, ancillary services and capacity markets given the amount of installed capacity on the grid. The profit maximization investment problem for each GENCO is solved using a global search algorithm, which uses the predictive model to evaluate the objective function. To solve for the strategic equilibrium, each GENCO’s problem is plugged into a diagonalization algorithm that is generally used in multi-leader, single-follower bi-level problems.

Dalvi, Sourabh↗
DOE OSTI2020

Physical effects on compact high-velocity clouds in the circumgalactic medium

ABSTRACT We numerically investigate the evolution of compact high-velocity clouds (CHVCs) passing through a hot, tenuous gas representing the highly ionized circumgalactic medium (CGM) by applying the adaptive-mesh refinement code flash. The model clouds start from both hydrostatic and thermal equilibrium and are in pressure balance with the CGM. Here, we present 14 models, divided into two mass categories and two metallicities each and different velocities. We allow for self-gravity and thermal conduction or not. All models experience mass diffusion, radiative cooling, and external heating leading to dissociation and ionization. Our main findings are (1) self-gravity stabilizes clouds against Rayleigh–Taylor instability, which is disrupted within 10 sound-crossing times without; (2) clouds can develop Jeans-instable regions internally even though they are initially below Jeans mass; (3) all clouds lose mass by ram pressure and Kelvin–Helmholtz instability; (4) thermal conduction substantially lowers mass-loss rates, by this, extending the clouds’ lifetimes, particularly, more than doubling the lifetime of low-mass clouds; (5) thermal conduction leads to continuous, filamentary stripping, while the removed gas is heated up quickly and mixes efficiently with the ambient CGM; (6) without thermal conduction the removed gas consists of dense, cool, clumpy fragments; (7) thermal conduction might prevent CHVCs from forming stars; and (8) clouds decelerated by means of drag from the ambient CGM form head-tail shapes and collapse after they reach velocities characteristic for intermediate-velocity clouds. Conclusively, only sophisticated modelling of CHVCs as non-homogeneous and non-isothermal clouds with thermal conduction and self-gravity explains observed morphologies and naturally leads to the suppression of star formation.

Sander, Bastian↗
DOE OSTI2024

Enhanced prediction of Cs removal by CST from Hanford tank waste with K accountability

Abstract The treatment of Hanford tank waste is one of the most challenging environmental cleanup activities to date. To expedite the processing of liquid waste stored in underground tanks in Washington State it is necessary to remove the significant dose contributor, 137 Cs. Crystalline silicotitanate ion exchanger is currently used to remove 137 Cs from the aqueous phase of Hanford tank wastes in preparation for vitrification at the Waste Treatment and Immobilization Plant (WTP). Improving the understanding of potassium impacts on ion exchange behavior of Cs will help in the operation of a critical component of one of the most complex treatment processes in the world today. Optimization of this process can result in significant cost savings and less waste production. Toward this effort, a series of batch contact tests varied in potassium concentration were conducted to look at the impact of potassium concentration on Cs distribution. Experimental distribution ratios ( K d ) were compared to the distribution ratios predicted using the ZAM model. A significant underprediction of Cs capacity in the presence of potassium was seen with the existing model. A revision of the equilibrium constants was determined and provided a statistically better fit for determining the Cs K d values in tank waste matrices.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Anomalous near-equilibrium capillary imbibition induced by nanoscale surface topography

Here, this work reports that the nanoscale surface topography induces a crossover from the classical Lucas–Washburn imbibition dynamics to a long-lived anomalously slow regime observed near thermodynamic equilibrium conditions. A compact analytical model considering the interplay between interfacial wetting, thermal motion, and the nanoscale physical topography accounts for experimental observations for the case of capillary rise of different nonvolatile liquids in simple glass tubes. Analytical predictions supported by topographic analysis via atomic force microscopy indicate that topographic features with a specific range of nanoscale dimensions determine the crossover condition and anomalous imbibition rate. Our findings have important implications for the scientific understanding and technical application of capillary imbibition and suggest strategies to control the adsorption of specific liquids in porous materials with complex surface topography.

Capillary driven imbibition↗
DOE OSTI2023

Part I: Predicting performance of Purolite A532E resins for remediation of comingled contaminants in groundwater

Ion exchange (IX) resins are used in pump-and-treat (P&T) facilities to remove soluble groundwater contaminants. However, natural anions present at concentrations orders of magnitude higher than contaminants can compete for IX sites and impact resin lifecycles. Here, the Hanford Site’s 200 West Area P&T facility (Washington State, USA) was selected as a case study because it currently uses two IX resins: Purolite® A532E (A532E) to remove pertechnetate (TcO 4 - ) and DOWEX 21K (DOWEX) to remove uranium from groundwater. Nitrate (NO 3 - ), sulfate (SO 4 2- ), chloride (Cl - ), and carbonate (CO 3 2- ) anions have been identified to potentially compete for A532E and DOWEX IX sites. Hanford-relevant anion groundwater concentrations were used to design a series of laboratory-scale batch experiments to evaluate the impact of competing anions on resin performance and potential kinetic effects. These data are then modeled to obtain Cl--normalized equilibrium exchange coefficients (K) needed to predict IX resin performance. The work is presented in two parts, with IX performance evaluated for A532E in Part I and DOWEX in Part II. Part I results demonstrate that TcO 4 - uptake is not impacted by NO 3 - , SO 4 2- , Cl - , CO 3 2- (as HCO 3 - ) and U(VI) carbonate anions, with K TcO4-/Cl- > 4,000, likely due to the high selectivity of A532E trihexylammonium sites for the large, weakly hydrated TcO 4 - anion. Other anion K values were K NO3-/Cl- = 20, K SO4--/Cl- = 0.2, K HCO3-/Cl- = 0.09, K U/Cl- = 370–1000. These K values provide conservative parameters for predicting A532E performance, and demonstrate that, under these test conditions, A532E will remove TcO 4 - from current and future influent streams to meet groundwater treatment objectives.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Part II: Predicting performance of $\mathrm{DOWEX 21K}$ resin for remediation of comingled contaminants in groundwater

The selectivity of ion exchange (IX) resins for aqueous contaminant removal can be impacted by changing concentrations of competing natural groundwater ions. In a two-part investigation, the Hanford Site 200 West Area pump-and-treat (P&T) facility in Washington State, USA is used as a case study to evaluate the performance of two IX resins for groundwater treatment: Purolite® A532E for pertechnetate (TcO 4 - ) removal, explored in Part I, and DOWEX 21K (DOWEX) for uranium (U) removal. In Part II, DOWEX selectivity for U, as uranyl carbonate species, and uptake kinetics is quantified in a series of laboratory-scale aqueous batch experiments containing Hanford-relevant concentrations of competing anions nitrate (NO 3 - ), sulfate (SO 4 2- ), chloride (Cl - ), and carbonate (CO 3 2- ), as well as co-mingled contaminant TcO 4 - . Here the results demonstrate that DOWEX trimethylammonium functional groups are highly selective for U carbonate species (85–100 % uptake) under all conditions investigated. Only NO 3 - concentrations of 100 mM were shown to decrease U removal, with the extent (85–99 %) depending on competing anion concentrations present in solution. However, at the highest NO 3 - concentrations reported for groundwaters treated at the P&T facility (25 mM), the effect on U uptake is minimal. The batch sorption results are modeled to obtain chloride normalized equilibrium exchange coefficients (K) for predicting DOWEX performance: K SO4--/Cl- = 2.0, K NO3-/Cl- = 5.0, K HCO3-/Cl- = 1.5, K TcO4-/Cl- = 2,000, and K U/Cl- = 50,000. These K values predict little effect of current and future influent chemistries on U removal by DOWEX, where both uranyl carbonate species and TcO 4 - are removed such that effluent concentrations meet groundwater treatment requirements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Defect cluster and nonequilibrium gas bubble associated growth in irradiated UMo fuels – a cluster dynamics and phase field model

Irradiation examination shows that gas bubble swelling kinetics is much faster after irradiation-induced recrystallization than that prior recrystallization in UMo fuels. It implies that gas bubbles in coarse grains and small recrystallized grains have different growth behavior. In this work, a phase-field model of gas bubble evolution integrating microstructure dependent cluster dynamics has been developed, for the first time, to study the gas bubble swelling behavior in the recrystallization zone of UMo fuels. Generation, diffusion, reaction, sink, emission and clustering of vacancies and interstitials are described by the cluster dynamics model while a phase-field model is used to describe the evolution of non-equilibrium gas bubbles including nucleation and growth. With the coupled model, the effect of defect generation rate, clustering rate, interstitial emission and sink rates on grain boundaries on the gas bubble evolution are systematically simulated. A set of model parameters (defect generation rate, clustering rate, interstitial emission and sink rates) is determined by comparing measured and simulated gas bubble swelling kinetics. The results demonstrate that interstitial clustering is one of the important physical mechanisms which results in a fast gas bubble swelling kinetics in the recrystallization zone. The developed model can also be extended to study the associated growth of defect and second phase precipitates often observed in irradiated materials.

Hu, Shenyang↗
DOE OSTI2024

A structural underpinning of the lower critical solution temperature (LCST) behavior behind temperature-switchable liquids

In this work, we use state-of-the-art X-ray scattering and molecular dynamics to analyze amine-water mixtures that show the unusual lower critical solution temperature (LCST) behavior. The goal is to provide direct experimental evidence for the entropy-lowering molecular cluster formation hypothesized as necessary for LCST behavior. Differential wide-angle X-ray scattering and pair distribution analysis and small-angle X-ray scattering measurements were combined with molecular modeling and liquid-liquid equilibrium measurements, revealing direct experimental evidence for the hypothesis. However, the response of the amine phase to accommodating water is even more subtle than the simple hypothesis suggests, with the formation of robust nanoscale reverse micelles. The techniques developed in this paper can be expected to yield insights in the use of temperature-switchable liquids in solvent extraction and other separations, and the stabilization of organelles in living cells that do not have physical membranes but do require compositional gradients to operate.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The impact of market design and clean energy incentives on strategic generation investments and resource adequacy in low-carbon electricity markets

Well-designed electricity markets play a crucial role in maintaining reliable electric power systems, which are critical in modern society. Here, this study examines the impact of different electricity market designs and clean energy incentive schemes on supporting renewable energy integration and achieving clean energy goals. To this end, we utilize a game-theoretical generation expansion planning model where generation companies make investment and retirement decisions to maximize their expected profit. The model is structured as an equilibrium problem with equilibrium constraints (EPEC) and solved using a diagonalization approach combined with progressive hedging. We analyze three types of electricity market designs: an energy-only market, a capacity market, and a clean energy market, and consider a wide range of market parameters resulting in 14 total scenarios. Wind and solar capacity comprise the majority of new investments in all considered scenarios, but the resultant system planning reserve margin (PRM) can differ significantly depending on market parameters. We also find that profit-driven investments lead to lower PRMs than a traditional system cost minimization approach. These individual scenario results further demonstrate how different market designs and clean energy incentive schemes may influence investor decision-making and impact resource adequacy throughout the clean energy transition.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Adsorption Equilibrium, Kinetics, and Column Breakthrough Data of Acetic Acid, Butyric Acid, and Lactic Acid on IRN-78 Ion-exchange Resin at Initial pH ~3 – 7 and Temperature 25 – 55 °C

This work presents systematic aqueous-phase adsorption equilibrium, kinetics, and column breakthrough measurements with three key biointermediates that are common compounds in many bioprocesses. Adsorption equilibrium experiments were carried out with acetic acid, butyric acid, and lactic acid on a commercial ion-exchange resin, Amberlite IRN-78, at wide ranges of acid concentration (8-500 mmol/L), initial pH (similar to 3-7), and temperature (25-55 degrees C), simulating the effluent characteristics from different fermenter operations. The kinetics and column breakthrough experiments were conducted at an initial pH of 6 and a concentration of 200 mmol/L. The equilibrium study shows a higher loading at the initial pH < pK(a) and a lower loading at the initial pH > pK(a). Overall removal varies between 16 and 99% depending on the initial pH, temperature, and organic acid concentration and type. The study further indicates monolayer adsorption at the equilibrium pH > 10 and multilayer adsorption at the equilibrium pH < 6. The thermodynamic modeling of adsorption isotherm data was carried out using Langmuir and Freundlich isotherms. IRN-78 presents fast adsorption kinetics as the maximum loading was attained in <= 10 min and nearly the same breakthrough time for all three organic acids involved in this study.

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

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ice Supersaturation Variability in Cirrus Clouds: Role of Vertical Wind Speeds and Deposition Coefficients

Abstract Aircraft measurements reveal ice supersaturation statistics in cirrus (ISSs) with broad maxima around ice saturation and pronounced variance. In this study, processes shaping ISSs in midlatitude and tropical upper tropospheric conditions are systematically investigated. Water vapor deposition and sublimation of size‐resolved ice crystal populations are simulated in an air parcel framework. Mesoscale temperature fluctuations (MTFs) due to gravity waves force the temporal evolution of supersaturation. Various levels of background wave forcing and cirrus thickness are distinguished in stochastic ensemble simulations. Kinetic limitations to ice mass growth are brought about by supersaturation‐dependent deposition coefficients that represent efficient and inefficient growth modes as a function of ice crystal size. The simulations identify a wide range of deposition coefficients in cirrus, but most values stay above 0.01 such that kinetic limitations to water uptake remain moderate. Supersaturation quenching times are long, typically 0.5–2 hr. The wave forcing thus causes a remarkably large variability in ISSs and cirrus microphysical properties except in the thickest cirrus, producing ensemble‐mean ISSs in line with in‐situ measurements. ISS variance is controlled by MTFs and increases with decreasing cirrus integral radii. In comparison, the impact of ice crystal growth rates on ISSs is small. These results contribute to efforts directed at identifying and solving issues associated with ice‐supersaturated areas and non‐equilibrium cirrus physics in global models.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Tax Credits for Clean Electricity: The Distributional Impacts of Supply-Push Policies in the Power Sector

We evaluate distributional and efficiency consequences of the bulk power clean electricity tax credits authorized by the 2022 Inflation Reduction Act. To do so, we link detailed electricity capacity expansion, computable general equilibrium, microsimulation, and air pollution models to estimate economic welfare and health incidence across demographic groups. We evaluate trade-offs between policy efficiency and income progressivity by comparing the tax credits to cap-and-trade policies. The tax credits encourage increased clean electricity investment, resulting in a reallocation of capital from elsewhere in the economy, higher prices for capital and other goods, lower power prices, and lower emissions. The tax credits yield progressive outcomes for economic welfare at the expense of efficiency while all modeled policies demonstrate progressivity in health impacts. The health benefits, absent climate benefits, exceed total policy costs and provide greater benefits for low-income and historically marginalized households given coincidence of household locations and emissions exposure intensity.

distributional impacts↗
DOE OSTI2022

Radiative pulsed L-mode operation in ARC-class reactors

Abstract A new ARC-class, highly-radiative, pulsed, L-mode, burning plasma scenario is developed and evaluated as a candidate for future tokamak reactors. Pulsed inductive operation alleviates the stringent current drive requirements of steady-state reactors, and operation in L-mode affords ELM-free access to ∼ 90 % core radiation fractions, significantly reducing the divertor power handling requirements. In this configuration the fusion power density can be maximized despite L-mode confinement by utilizing high-field to increase plasma densities and current. This allows us to obtain high gain in robust scenarios in compact devices with P fus > 1000 MW despite low confinement. We demonstrate the feasibility of such scenarios here; first by showing that they avoid violating 0D tokamak limits, and then by performing self-consistent integrated simulations of flattop operation including neoclassical and turbulent transport, magnetic equilibrium, and radiofrequency current drive models. Finally we examine the potential effect of introducing negative triangularity with a 0D model. Our results show high-field radiative pulsed L-mode scenarios are a promising alternative to the typical steady state advanced tokamak scenarios which have dominated tokamak reactor development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Solar Thermochemical Carbon Dioxide Splitting Using Ceria and Iron Aluminate Foam Devices and Simulation of a Plant System for Demonstration

An international research project has been undertaken to integrate a unique solar thermal processing reactor system with ceria and iron aluminate as active redox materials for CO2 splitting. Experimental investigations for CO2 splitting were conducted using a solar simulator and tube furnace at Niigata University, followed by demonstrations using a high-flux solar furnace (HFSF) at the National Renewable Energy Laboratory (NREL) in Golden, CO. Each experimental setup consisted of foam devices composed of reticulated porous ceramic (RPC). The RPC has a full ceria or iron aluminate body. It fabricated using the replica method and subjected to a two-step redox reaction, which iteratively separated a stream of CO2 into O2 and CO. Reactivity was evaluated using CO production per mass of the reactive material. The tubular furnace yielded a CO production of 6.41 mL/g at a reduction temperature of 1600degrees C, showcasing a higher CO production rate and total amount than those obtained from experiments conducted with solar simulators and solar furnace setups. For iron aluminate RPC, the productivity was measured as 3.57 mL/g using HFSF at a reduction temperature of 1450degrees C. These results are somewhat higher than those of the previous experiment at lower reduction temperatures of 1400degrees C-1500degrees C. Additionally, the production of CO in the case of ceria RPC was compared with the steady flow model simulation, which assumed chemical equilibrium at various levels of oxygen partial pressure during the reduction process. On the basis of these results, this study proposes a solar fuel system with an open receiver that uses a high-temperature heat transfer fluid.

carbon dioxide thermochemical splitting↗
DOE OSTI2024

Non-Thermal Solar Wind Electron Velocity Distribution Function

The quiet-time solar wind electrons feature non-thermal characteristics when viewed from the perspective of their velocity distribution functions. They typically have an appearance of being composed of a denser thermal “core” population plus a tenuous energetic “halo” population. At first, such a feature was empirically fitted with the kappa velocity space distribution function, but ever since the ground-breaking work by Tsallis, the space physics community has embraced the potential implication of the kappa distribution as reflecting the non-extensive nature of the space plasma. From the viewpoint of microscopic plasma theory, the formation of the non-thermal electron velocity distribution function can be interpreted in terms of the plasma being in a state of turbulent quasi-equilibrium. Such a finding brings forth the possible existence of a profound inter-relationship between the non-extensive statistical state and the turbulent quasi-equilibrium state. The present paper further develops the idea of solar wind electrons being in the turbulent equilibrium, but, unlike the previous model, which involves the electrostatic turbulence near the plasma oscillation frequency (i.e., Langmuir turbulence), the present paper considers the impact of transverse electromagnetic turbulence, particularly, the turbulence in the whistler-mode frequency range. It is found that the coupling of spontaneously emitted thermal fluctuations and the background turbulence leads to the formation of a non-thermal electron velocity distribution function of the type observed in the solar wind during quiet times. This demonstrates that the whistler-range turbulence represents an alternative mechanism for producing the kappa-like non-thermal distribution, especially close to the Sun and in the near-Earth space environment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Structure Formation Through Magnetohydrodynamical Instabilities in Primordial Disks

The shear flow instabilities under the presence of magnetic fields in the primordial disk can greatly facilitate the formation of density structures that serve as seeds prior to the onset of the gravitational Jeans instability. We evaluate the effects of the Parker, magnetorotational and kinematic dynamo instabilities by comparing the properties of these instabilities. We calculate the mass spectra of coagulated density structures by the above mechanism in the radial direction for an axisymmetric magnetohydrodynamic (MHD) torus equilibrium and power density profile models. Our local three-dimensional MHD simulation indicates that the coupling of the Parker and magnetorotational instability creates spiral arms and gas blobs in an accretion disk, reinforcing the theory and model. Such a mechanism for the early structure formation may be tested in a laboratory. The recent progress in experiments involving shear flows in rotating tokamak, field reversed configuration (FRC) and laser plasmas may become a key element to advance in nonlinear studies.

Noguchi, Koichi↗