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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

[Two-phase Interfaces in Weak External Fields]

Our aim has been that of understanding from first principles the behavior of two-phase interfaces in the absence of gravitational constraints. This is fundamental to our ability to deal with the fluid structures that abound in the real biological, chemical, and physical world. A substantial effort was mounted to determine how familiar hydrodynamic concepts have to be modified and interpreted to make them appropriate to the multi-level structure alluded to above. This was primarily in the context of the microscopic symmetric pressure tensor, which was, for the first time, expressed in the invaluable density functional format, and the used to follow the predictions of popular microscopic models of the energetics of interfacial systems. In the course of these investigations, the previous murky relation between pressure tensor and thermodynamics was completely clarified. The process of extending thermodynamic information to interfacial dynamics was initiated along two paths. One was from the viewpoint of an inertialess lattice gas, resulting in the surprising conclusion that at this level, all transport is governed by precisely the thermodynamic free energy, albeit with a non-trivial effective particle mobility. The other aimed at understanding the fashion in which slow macroscopic motions, accounted for by a time-varying microscopic energy, generate effective hydrodynamic parameters. By examining a solvable model system, it was found that all current procedures for doing so are deficient, and suitable alleviation suggested. The major effect of this project was to set the stage for the analysis of the substantial dynamical regimes in which extensive equilibrium information provides the dominant background. This produces a smooth junction to the models of Araki and Munakata, Giacomin and Lebowitz, and Oxtoby. It is also crucial to our understanding of the complex interfacial equilibrium configurations required for intermediate stages of two-phase separation, for which the "phase-field" techniques we have been developing are uniquely effective. And it puts us within striking range of the physical experiments that can provide incisive tests of the theoretical framework.

Percus, J. K.↗
DOE OSTI2020

Thermodynamic Origin of Reaction Non-Uniformity in Battery Porous Electrodes and Its Mitigation

The development of non-uniform reaction current distribution within porous electrodes is a ubiquitous phenomenon during battery charging/discharging and frequently controls the rate performance of battery cells. Reaction inhomogeneity in porous electrodes is usually attributed to the kinetic limitation of mass transport within the electrolyte and/or solid electrode phase. In this work, however, we reveal that it is also strongly influenced by the intrinsic thermodynamic behavior of electrode materials, specifically the dependence of the equilibrium potential on the state of charge: the electrode reaction becomes increasingly non-uniform when the slope of the equilibrium potential curve is reduced. We employ numerical simulations and equivalent circuit model to elucidate such a correlation and show that the degree of reaction inhomogeneity and the resultant discharge capacity can be predicted by a dimensionless reaction uniformity number. For electrode materials that have equilibrium potentials insensitive to the state of charge and exhibit significant reaction non-uniformity, we demonstrate several approaches to spatially homogenizing the reaction current inside porous electrodes, including matching the electronic and ionic resistances, introducing graded electronic conductivity and reducing the surface reaction kinetics.

25 ENERGY STORAGE↗
DOE OSTI2021

Low-temperature domain-wall freezing and nonequilibrium dynamics in the transverse-field Ising model material CoNb 2 O 6

CoNb 2 O 6 is a rare realization of the transverse-field Ising model, making it a useful tool for studying both equilibrium and nonequilibrium many-body quantum physics. Despite a large body of work dedicated to characterizing this material, details of the ordered states in the presence of relatively weak transverse fields have not been discussed in detail. Here, we present a detailed study of CoNb 2 O 6 via ac susceptibility measurements in order to further characterize its low-temperature behavior in the presence of a transverse field. Specifically, we call attention to an unconventional freezing transition in zero field below T F = 1.2 K, occurring within the well-known commensurate antiferromagnetic (CAFM) state that onsets at T N2 = 1.9 K. We performed a series of transverse-field quenches into this frozen state, which resulted in a slowly relaxing susceptibility, χ' (t), that followed a logarithmic decay within the time range measured. Additionally, we discuss the frozen state in the context of the freezing of previously discussed “free” chains arising from domain walls between the four degenerate sublattices of the CAFM state. We also attempted to observe Kibble-Zurek scaling by quenching the transverse field into the frozen state at different rates. This produced a null result; the behavior can be fully explained by coarsening of domains over the time scale of the quenches. The absence of a clear Kibble-Zurek scaling is itself surprising, given the proposed ubiquity of the phenomenon for general second-order phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1994

Computation of iodine species concentrations in water

During an evaluation of the use of iodine as a water disinfectant and the development of methods for measuring various iodine species in water onboard Space Freedom, it became necessary to compute the concentration of the various species based on equilibrium principles alone. Of particular concern was the case when various amounts of iodine, iodide, strong acid, and strong base are added to water. Such solutions can be used to evaluate the performance of various monitoring methods being considered. The authors of this paper present an overview of aqueous iodine chemistry, a set of nonlinear equations which can be used to model the above case, and a computer program for solving this system of equations using the Newton-Raphson method. The program was validated by comparing results over a range of concentrations and pH values with those previously presented by Gottardi for a given pH. Use of this program indicated that there are multiple roots to many cases and selecting an appropriate initial guess is important. Comparison of program results with laboratory results for the case when only iodine is added to water indicates the program gives high pH values for the iodine concentrations normally used for water disinfection. Extending the model to include the effects of iodate formation results in the computer pH values being closer to those observed, but the model with iodate does not agree well for the case in which base is added in addition to iodine to raise the pH. Potential explanations include failure to obtain equilibrium conditions in the lab, inaccuracies in published values for the equilibrium constants, and inadequate model of iodine chemistry and/or the lack of adequate analytical methods for measuring the various iodine species in water.

Schultz, John R.↗
NASA NTRS1981

Effects of diffusion and mass flows on C IV and Si IV lines formed in the solar atmosphere

A model for the transition region is derived from an initial interpretation of EUV observations, assuming ionization equilibrium and constant elemental abundance with height. The effects of diffusion and mass flows are then included in the initial model and the emergent profiles of several C IV and Si IV lines are computed. It is found that diffusion and mass flows have a strong effect on both the emergent intensity and spectral shape of these lines. Diffusion acts to deplete the transition region of heavy ions to an extent which depends on the detailed temperature and density structure. The net effect is a weighting of the coronal emission relative to that in the transition. On the basis of the results, it is suggested that the downflows observed in the network in lines of C IV and Si IV could be due to gravitational settling of the ions following their injection, via diffusion, into the corona from spicules. When flows are superposed on the basic diffusion model, the rapid change in elemental abundance characteristic of the model is virtually eliminated and the coronal contribution to the emission in C IV and Si IV lines becomes negligible relative to that in the transition region. Flows have a strong effect on the computed line intensities and introduce large asymmetries into the line shapes.

Roussel-Dupre, R.↗
DOE OSTI2023

Gas-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

This report details a fully coupled neutronics thermal hydraulics reference plant model for a gas-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission's Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory high performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON MOOSE-based applications. Griffin provides the reactor physics capabilities including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient solutions. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions. BISON solves thermal conduction problems for the average pebble and TRISO in the pebble-bed core thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel and moderator temperatures. We present results for the coupled steady state equilibrium core and a protected loss of flow event. Although this model is to be considered prototypical with regard to the capabilities in BlueCRAB, its results are consistent with published work by the Institute of Nuclear and New Energy Technology (INET) in China and other research entities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1981

On the equilibrium distribution of the elements in the gas in the Coma cluster

The equilibrium distribution of the elements for the case of an isothermal plasma in a gravitational potential well is given governing equations and a rapidly converging iterative procedure for their solution. Results are applied to the intracluster gas in clusters of galaxies, and implications are identified for the controversy as to whether iron nuclei settle out in the cluster core, and with respect to the interpretation of distribution size and shape for continuum X-ray emission from clusters. A model for the Coma cluster is considered in which the intracluster gas is in equilibrium and partially processed. It is found that this model can account for the angular dependence of the continuum X-ray emission from Coma, as determined by Einstein Observatory observations, at least as well as the uniform composition model.

Abramopoulos, F.↗
DOE OSTI2021

A Consensus Equilibrium Approach for 3-D Land Seismic Shots Recovery

Physical and budget constraints often result in inadequate sampling for accurate subsurface imaging. Preprocessing approaches, such as missing trace interpolation, are typically employed to enhance seismic data in such cases. The compressed sensing (CS) framework has been applied for modeling missing seismic data, which is estimated by sparsity-based computational algorithms. While existing work mainly focuses on recovering missing traces resulting from receiver subsampling, source subsampling has greater economical advantages, as sources are more expensive than receivers. Moreover, stronger image models different from sparsity have not been explored for source recovery. This work presents a consensus equilibrium (CE) approach to recover missing seismic shots, which enables to incorporate various regularization operators modeling different data priors. Here, simulation results from a real 3-D land seismic dataset demonstrate that the CE approach provides more accurate estimations of the linear and hyperbolic events in the recovered shots, compared with pure sparsity-based reconstructions.

58 GEOSCIENCES↗
DOE OSTI2021

Role of Equilibrium Fluctuations in Light-Induced Order

Engineering novel states of matter with light is at the forefront of materials research. An intensely studied direction is to realize broken-symmetry phases that are “hidden” under equilibrium conditions but can be unleashed by an ultrashort laser pulse. Despite a plethora of experimental discoveries, the nature of these orders and how they transiently appear remain unclear. To this end, we investigate a nonequilibrium charge density wave (CDW) in rare-earth tritellurides, which is suppressed in equilibrium but emerges after photoexcitation. Using a pump-pump-probe protocol implemented in ultrafast electron diffraction, we demonstrate that the light-induced CDW consists solely of order parameter fluctuations, which bear striking similarities to critical fluctuations in equilibrium despite differences in the length scale. By calculating the dynamics of CDW fluctuations in a nonperturbative model, we further show that the strength of the light-induced order is governed by the amplitude of equilibrium fluctuations. Overall, these findings highlight photoinduced fluctuations as an important ingredient for the emergence of transient orders out of equilibrium. Our results further suggest that materials with strong fluctuations in equilibrium are promising platforms to host hidden orders after laser excitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Minimum conditions for accurate modeling of urea production via co-electrolysis

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

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

Thermodynamic Modeling of Point Defects in Triple Conducting Perovskite Ba 0.95 La 0.05 FeO 3- $_δ$ with Incorporation of the Hydride Defect Formation Reaction for Solid Oxide Cells

Distinct from the proton defect, the hydride defect species may be present in certain perovskite materials in reducing environments such as in fuel electrodes of solid oxide cells (SOCs) or on the reducing side of ceramic membranes. A generalized defect thermodynamic model was developed for the triple-conducting perovskites (La,Ba)Fe 1-x M x O 3-δ (M = Y and Zr) to allow inclusion of the hydride defect formation reaction in addition to the other three main defect reactions, namely, the oxygen vacancy formation, hydration, and charge disproportionation reactions. This comprehensive defect model also allows the incorporation of polynomial functional forms of oxygen nonstoichiometry δ to describe the defect reaction energies and entropies and to enable refinements of the defect reaction equilibrium constants in the defect thermodynamic analysis. As a first step, the developed model is applied to the Ba 0.95 La 0.05 FeO 3-δ material as an illustrative system to obtain its Brouwer diagrams with both the proton and hydride defects in relevant SOC conditions, particularly for more reducing environments. In conclusion, the results provide direct guidance on the influence of electronic and ionic defect concentrations upon thermodynamic properties and ultimately on the performance of Ba 0.95 La 0.05 FeO 3-δ and potentially other (La,Ba)Fe 1-x M x O 3-δ perovskite materials involved in SOC applications.

25 ENERGY STORAGE↗
NASA NTRS1996

Emission Spectral Measurements in the Plenum of an Arc Jet Wind Tunnel

Arc jet wind tunnel facilities are used to evaluate thermal protection system materials for re-entry vehicles. The high speed, high temperature flowfield generated by the arc jet can simulate the extreme aerodynamic heating environment experienced during re-entry so that the survivability of heat shield materials and performance of various designs options can be tested. Although the re-entry heating environment can be approximated in the arc jet facility, the flowfield only partially simulates the actual re-entry flight conditions. Reynolds numbers are not matched so that surface shear stress distributions and mass transfer rates due to ablation or other mechanisms are not modeled correctly. Unlike flight conditions the arc freestream air is in non-equilibrium because of the rapid expansion that occurs in the supersonic nozzle. To properly study the actual re-entry flow environment, computational fluid dynamics, computational chemistry and radiation models must be used. Arc jet tunnel tests serve to validate these models. To perform accurate simulations inlet and boundary-conditions are needed, which come from measurements of the flowfield. The present study is concerned with measurements in the plenum region of an arc heater. In the past, conditions in the arc heater flowfield have been predicted using simulations since conventional measurement techniques could not be used in the harsh extremely high temperature environment. The present study is part of a recent push to utilize optical techniques to help better characterize the arc jet flowfields. Emission measurements have been made in the shock layer and the constrictor section of the arc heater to determine temperatures and species number densities. LIF measurements have been made in the free stream to determine temperature and velocity.

Donohue, Jim↗
NASA NTRS1992

Fractal landscapes in biological systems: long-range correlations in DNA and interbeat heart intervals

Here we discuss recent advances in applying ideas of fractals and disordered systems to two topics of biological interest, both topics having common the appearance of scale-free phenomena, i.e., correlations that have no characteristic length scale, typically exhibited by physical systems near a critical point and dynamical systems far from equilibrium. (i) DNA nucleotide sequences have traditionally been analyzed using models which incorporate the possibility of short-range nucleotide correlations. We found, instead, a remarkably long-range power law correlation. We found such long-range correlations in intron-containing genes and in non-transcribed regulatory DNA sequences as well as intragenomic DNA, but not in cDNA sequences or intron-less genes. We also found that the myosin heavy chain family gene evolution increases the fractal complexity of the DNA landscapes, consistent with the intron-late hypothesis of gene evolution. (ii) The healthy heartbeat is traditionally thought to be regulated according to the classical principle of homeostasis, whereby physiologic systems operate to reduce variability and achieve an equilibrium-like state. We found, however, that under normal conditions, beat-to-beat fluctuations in heart rate display long-range power law correlations.

Non-NASA Center↗
DOE OSTI2021

High‐Efficiency Ion‐Exchange Doping of Conducting Polymers

Abstract Molecular doping—the use of redox‐active small molecules as dopants for organic semiconductors—has seen a surge in research interest driven by emerging applications in sensing, bioelectronics, and thermoelectrics. However, molecular doping carries with it several intrinsic problems stemming directly from the redox‐active character of these materials. A recent breakthrough was a doping technique based on ion‐exchange, which separates the redox and charge compensation steps of the doping process. Here, the equilibrium and kinetics of ion exchange doping in a model system, poly(2,5‐bis(3‐alkylthiophen‐2‐yl)thieno(3,2‐b)thiophene) (PBTTT) doped with FeCl 3 and an ionic liquid, is studied, reaching conductivities in excess of 1000 S cm −1 and ion exchange efficiencies above 99%. Several factors that enable such high performance, including the choice of acetonitrile as the doping solvent, which largely eliminates electrolyte association effects and dramatically increases the doping strength of FeCl 3 , are demonstrated. In this high ion exchange efficiency regime, a simple connection between electrochemical doping and ion exchange is illustrated, and it is shown that the performance and stability of highly doped PBTTT is ultimately limited by intrinsically poor stability at high redox potential.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Joint estimation of biogeochemical model parameters from multiple experiments: A bayesian approach applied to mercury methylation

Here, to characterize complex biogeochemical systems, results from multiple experiments, where each targets a specific subprocess, are commonly combined. The resulting datasets are interpreted through the calibration of biogeochemical models for process inference and predictions. Commonly used calibration approaches of fitting datasets from individual experiments to subprocess models one at a time is prone to missing information shared between datasets and incomplete uncertainty propagation. We propose a Bayesian joint-fitting scheme addressing the above-mentioned concerns by jointly fitting all the available datasets, thus calibrating the entire biogeochemical model in one go using Markov Chain Monte Carlo (MCMC). The identification of null spaces in the parameter distributions from MCMC guided the simplification of certain subprocess models. For example, fast kinetic sorption was replaced by equilibrium sorption, and Monod demethylation was replaced by first-order demethylation. Joint fitting of datasets resulted in complete uncertainty propagation with parameter estimates informed by all available data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Phenalenyl growth reactions and implications for prenucleation chemistry of aromatics in flames

The energetics and kinetics of phenalene and phenalenyl growth reactions were studied theoretically. Rate constants of phenalene and phenalenyl H-abstraction and C 2 H 2 addition to the formed radicals were evaluated through quantum-chemical and rate-theory calculations. The obtained values, assigned to all π radicals, were tested in deterministic and kinetic Monte Carlo simulations of aromatics growth under conditions of laminar premixed flames. Kekulé and non-Kekulé structures of the polycyclic aromatic hydrocarbons (PAHs) evolving in the stochastic simulations were identified by on-the-fly constrained optimization. The numerical results demonstrated an increased PAH growth and qualitatively reproduced experimental observations of Homann and co-workers of non-decaying PAH concentrations with nearly equal abundances of even and odd carbon-atom PAHs. The analysis revealed that the PAH growth proceeds via alternating and sterically diverse acetylene and methyl HACA additions. Furthermore, the rapid and diverse spreading in the PAH population supports a nucleation model as PAH dimerization, assisted by the non-equilibrium phenomena, forming planar aromatics first and then transitioning to the PAH–PAH stacking with size.

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

Novel intrinsic helical cores and MHD dynamo flux pumping evidence in DIII-D

Evidence of a MHD dynamo brought about by intrinsic helical core structures in negative triangularity (NT) and ITER baseline scenario (IBS) plasmas demonstrates paths to operating tokamaks without seeds for deleterious global MHD modes. A stationary sawtooth- and edge localized mode-free state maintains high performance for long durations in a strong NT shape in the DIII-D tokamak. Disparity between minimum safety factor values predicted by magnetic flux diffusion theory and derived from experiment suggests anomalous flux diffusion by an MHD dynamo, driven by flows associated with quasi-interchange modes. Stability calculations with the GATO code find NT and IBS experimental equilibria unstable to quasi-interchange modes, in agreement with observed saturated activity. IBS discharges with similar magnetic winding structure occupy robust helical core states following a transition from sawtooth instabilities. These states agree with prior modeling of helical core bifurcation thresholds with the VMEC equilibrium code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗