SEARCH · Engineering Papers

Results for “thermodynamic driving force”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Microstructure-Dependent Rate Theory Model of Radiation-Induced Segregation in Binary Alloys

Conventional rate theory often uses the mean field concept to describe the effect of inhomogeneous microstructures on the evolution of radiation induced defect and solute/fission product segregation. However, the spatial and temporal evolution of defects and solutes determines the formation and spatial distribution of radiation-induced second phase such as precipitates and gas bubbles/voids, especially in materials with complicated microstructures and subject to high dose radiation. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed to investigate the effect of inhomogeneous thermodynamic and kinetics properties of defects on diffusion and accumulations of solute A in AB binary alloys. Four independent concentrations: atom A, interstitial A, interstitial B, and vacancy on [A, B] sublattice are used as field variables to describe temporal and spatial distribution and evolution of defects and solute A. The independent concentrations of interstitial A and interstitial B allow to describe their different generation rates, thermodynamic and kinetic properties, and release the assumptions of interstitial generation and sink strength used in the conventional rate theory. Microstructure and concentration dependent chemical potentials of defects are used to calculate the driving forces of defect diffusions. With the model, the effects of defect chemical potentials and mobilities on the RIS in polycrystalline AB model alloys have been simulated. The results demonstrated the model capability in predicting defect evolution in materials with inhomogeneous thermodynamic and kinetic properties of defects. The model can be extended to materials with complicated microstructures such as a wide range of grain size distribution, coating structure and multiphases as well as radiation-induced precipitation subject to severe radiation damage.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Nucleation kinetics model for primary crystallization in Al–Y–Fe metallic glass

The high density of aluminum nanocrystals (>10 21 m -3 ) that develop during the primary crystallization in Al-based metallic glasses indicates a high nucleation rate (~10 18 m -3 s -1 ). Several studies have been advanced to account for the primary crystallization behavior, but none have been developed to completely describe the reaction kinetics. Recently, structural analysis by fluctuation electron microscopy has demonstrated the presence of the Al-like medium range order (MRO) regions as a spatial heterogeneity in as-spun Al 88 Y 7 Fe 5 metallic glass that is representative for the class of Al-based amorphous alloys that develop Al nanocrystals during primary crystallization. From the structural characterization, an MRO seeded nucleation configuration is established, whereby the Al nanocrystals are catalyzed by the MRO core to decrease the nucleation barrier. The MRO seeded nucleation model and the kinetic data from the delay time (τ) measurement provide a full accounting of the evolution of the Al nanocrystal density (N v ) during the primary crystallization under isothermal annealing treatments. Moreover, the calculated values of the steady state nucleation rates (J ss ) predicted by the nucleation model agree with the experimental results. Moreover, the model satisfies constraints on the structural, thermodynamic, and kinetic parameters, such as the critical nucleus size, the interface energy, and the volume-free energy driving force that are essential for a fully self-consistent nucleation kinetics analysis. In conclusion, the nucleation kinetics model can be applied more broadly to materials that are characterized by the presence of spatial heterogeneities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Global decline in ocean memory over the 21st century

Ocean memory, the persistence of ocean conditions, is a major source of predictability in the climate system beyond weather time scales. We show that ocean memory, as measured by the year-to-year persistence of sea surface temperature anomalies, is projected to steadily decline in the coming decades over much of the globe. This global decline in ocean memory is predominantly driven by shoaling of the upper-ocean mixed layer depth in response to global surface warming, while thermodynamic and dynamic feedbacks can contribute substantially regionally. As the mixed layer depth shoals, stochastic forcing becomes more effective in driving sea surface temperature anomalies, increasing high-frequency noise at the expense of persistent signals. Reduced ocean memory results in shorter lead times of skillful persistence-based predictions of sea surface thermal conditions, which may present previously unknown challenges for predicting climate extremes and managing marine biological resources under climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Neuromorphic heat transport effects in a molecular junction

Understanding energy transport at the nanoscale is an open and fundamental challenge in the molecular sciences with direct implications for the design of new electronics, computing devices, and materials. While nanoscale energy transport under steady-state conditions has been studied extensively, there is much less known about energy transport under time-dependent driving forces, particularly in the far-from-equilibrium regime. In this work, we use nonequilibrium molecular dynamics simulations and stochastic thermodynamics to investigate energy transport in a well-studied nanoscale system—a molecular junction—subjected to a time-periodic temperature gradient. The primary observation is that molecular junctions can exhibit heat transport hysteresis, a phenomenon in which the heat flux through a system depends not only on the instantaneous value of a time-dependent temperature bias but also on the temporal history of that bias. The presented findings illustrate that molecular junctions can exhibit the specific memory effect—heat transport hysteresis—that is essential for the design of thermal neuromorphic computers. This work elucidates a potential pathway toward the realization of such devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

The Importance of an Ensemble Approach for Modeling Aerosol‐Convection Interactions

Aerosol-convection interactions modulate cloud microphysics, thermodynamics, and updraft intensity, contributing to climate-scale aerosol-radiative forcing. However, quantifying aerosol indirect effects in mixed-phase deep convection remains challenging due to uncertainties in parameterized physics and initial conditions driving nonlinear evolution of convective processes. This study investigates the convective updraft sensitivity to random initial temperature perturbations using an idealized ensemble modeling framework informed by in situ thermodynamic and aerosol observations from the DOE TRACER field campaign. We analyze the impact of small-scale initial perturbations on updraft velocity and supersaturation, and determine the ensemble size required to minimize stochastic internal variability. Results show that minor thermodynamic perturbations can produce updraft variability comparable to aerosol-induced changes reported in prior work. An ensemble of 10 members sufficiently reduces variability, enabling robust investigation of aerosol-related updraft invigoration and informing the design of future ensemble-based aerosol-convection interaction studies in environments with significant spatiotemporal mesoscale thermodynamic and aerosol heterogeneity.

aerosol-convection interactions↗

DOE OSTI2025

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗

DOE OSTI2020

Computational thermodynamics and its applications

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable or un- stable, when interacting with the surroundings. In this overview, fundamentals of thermodynamics are briefly reviewed through the combination of first and second laws of thermodynamics for open and nonequilibrium systems to demonstrate that the reversible equilibrium and irreversible nonequilibrium thermodynamics can be integrated to enhance the power and utilities of thermodynamics. The recent progresses in computational thermodynamics, the remaining challenges, and potential impacts in broad scientific fields are discussed here. It is shown that computational thermodynamics enables the modeling of thermodynamics of a state as a function of both external and internal variables and quantitative calculations of a broad range of properties of a multicomponent system in terms of first and second derivatives of energy, including not only equilibrium states when there are no driving forces for any internal processes and but also non-equilibrium states with driving forces for internal processes. Consequently, external constraints such as fixed strain and internal degree of freedoms such as ordering and defects can be described in a coherent framework and applied to materials design. Furthermore, two important but largely overlooked aspects in thermodynamics will be discussed, i.e. the rigorous application of statistical thermodynamics with the probability of configurations and their contributions to system properties, and the applications of second derivatives of energy with respect to either two extensive variables or two potentials or a mixture of them in terms of understanding and predicting emergent behaviors, critical phenomena, kinetic coefficients, and mechanical properties.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Energetics of Salt-Bearing Sodalites, Na 8 Al 6 Si 6 O 24 X 2 (X = SO 4 , ReO 4 , Cl, I): A Treatment Option for Pertechnetate-Enriched Nuclear Waste Streams

An alternative option for treating anion-enriched reprocessed nuclear waste streams is to immobilize technetium-99 ( 99 Tc, β = 293.7 keV, t 1/2 = 2.1 × 10 5 years) and other anions in micro- and mesoporous materials. Here we determine the thermodynamic stability of anion bearing sodalites, Na 8 Al 6 Si 6 O 24 X 2 (X = SO 4 , ReO 4 , Cl, I), to improve our understanding of the driving forces that control framework assembly using high temperature oxide melt solution calorimetry. Raman and FTIR spectroscopy illustrate a strong dependence for vibrational features on anion size and enabled the development of a linear model that predicted the vibrational features for numerous anion bearing sodalites to within ±20 cm –1 (i.e., OH, F, Br, ClO 4 , NO 3 , and MnO 4 ). The largest negative enthalpy of formation from elements and the lack of structural water demonstrate that the perrhenate sodalite (Na 8 Al 6 Si 6 O 24 [ReO 4 ] 2 ), a chemical analogue for pertechnetate sodalite (Na 8 Al 6 Si 6 O 24 [TcO 4 ] 2 ), is more thermodynamically stable than all other anion bearing sodalites evaluated. The enthalpies of the reaction between nepheline and the sodium salt, which provides the guest anion species, was negative only for the ReO 4 and NO 3 bearing sodalites. Overall, we report for the first time the enthalpy of the ion exchange reactions for different anion bearing sodalites relative to the perrhenate sodalite, which is a key step in gaining the ability to tune sodalite material properties and structure during treatment and the immobilization of 99 Tc in the presence of competing anions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗

DOE OSTI2021

A Time-Dependent Directional Damage Theory for Brittle Rocks Considering the Kinetics of Microcrack Growth

In this work, a novel microcrack damage theory for describing the time-dependent behavior of brittle rocks is proposed. Instead of using scalar or tensorial variables to approximate the effect of microcracks, the concept of directional damage is introduced by upscaling the directional distribution density of microcracks through the noninteracting homogenization scheme. The contribution of an individual crack on global strain is first derived from the crack-opening displacements including sliding, closure, and dilation. The model is then cast in the free energy form, and the range of parameters ensuring compliance with the second law of thermodynamics is derived rigorously. Considering that the subcritical propagation of microcracks is the dominant mechanism behind the time-dependent deformation and failure of brittle rocks, the driving force for directional damage is defined as the corresponding energy release rate, and the effect of time is introduced through a damage evolution law inspired from the kinetics of subcritical crack growth. The model is applied to simulate the behavior of basalt rock under various loading conditions. In addition to predicting the stress–strain–time response, the model offers enhanced resolution in representing the anisotropic characteristics of microcrack-induced damage in brittle rocks and their time-dependent evolutions.

58 GEOSCIENCES↗

DOE OSTI2026

Electrolyte-Dependent, “Microscopically Irreversible” H-Atom Transfer Kinetics of Ce-Based Metal–Organic Framework, Ce-MOF-808

Redox reactions at the interface of metal oxides and protic electrolytes almost always involve protons and electrons in equal amounts. Given the stoichiometry, these proton-coupled electron transfer (PCET) reactions are thermochemically equivalent to net H-atom transfer (HAT) reactions. The correlation between the chemical nature of solid catalysts and HAT kinetics has been employed for decades as the design principle for energy-relevant reactions (e.g., reactions of 2H + /H 2 ). More recently, chemists have experimentally determined that a change in liquid electrolytes that alters the microenvironment at the redox-active sites has an equally profound impact on electrocatalysis involving PCET/HAT. Yet, precise correlations between the chemical nature of electrolytes and the PCET kinetics are, to date, rare in the literature. Herein, we report our findings using the Ce-based metal−organic framework, Ce-MOF-808, as a model system. Each Ce 6 (μ 3 −O) 4 (μ 3 − OH) 4 (OH) 6 (H 2 O) 6 node of this MOF undergoes a 1H + /1e − redox reaction. Using chronoamperometry and the Cottrell analysis, we have determined that the PCET hopping kinetics within the pores of Ce-MOF-808 can change by orders of magnitude by altering the buffer species and the proton activity of the electrolyte. Furthermore, in all buffers, reductive reactions were ∼3−10 times faster in kinetics than the reverse oxidative reaction with the same electrochemical driving force, suggesting that the system, at first glance, violates the principle of microscopic reversibility. Isothermal titration calorimetry (ITC) and computational simulations corroborated that the buffer-node binding thermodynamics are quite distinct, depending on the chemical nature of the buffer and the oxidation state of the node. Together, these results suggest that the substrate and the product during the oxidative vs reductive reaction of Ce-MOF-808 are chemically different species, which explains the apparent ‘microscopic irreversibility.’ Thus, the rational modulation of electrolytes can dramatically enhance PCET kinetics, even though the solid electrodes remain identical. Implications of these findings are contrasted with the electrochemical/electrocatalytic behavior of other redox-active MOFs, heterogeneous catalysts, and enzymatic systems at the solid−liquid interface.

Ce-based MOF↗

DOE OSTI2022

f‐ block MOFs: A Pathway to Heterometallic Transuranics

Abstract A novel series of heterometallic f ‐block‐frameworks including the first examples of transuranic heterometallic 238 U/ 239 Pu‐metal–organic frameworks (MOFs) and a novel monometallic 239 Pu‐analog are reported. In combination with theoretical calculations, we probed the kinetics and thermodynamics of heterometallic actinide(An)‐MOF formation and reported the first value of a U‐to‐Th transmetallation rate. We concluded that formation of uranyl species could be a driving force for solid‐state metathesis. Density of states near the Fermi edge, enthalpy of formation, band gap, proton affinity, and thermal/chemical stability were probed as a function of metal ratios. Furthermore, we achieved 97 % of the theoretical maximum capacity for An‐integration. These studies shed light on fundamental aspects of actinide chemistry and also foreshadow avenues for the development of emerging classes of An‐containing materials, including radioisotope thermoelectric generators or metalloradiopharmaceuticals.

Park, Kyoung Chul↗

DOE OSTI2022

f‐ block MOFs: A Pathway to Heterometallic Transuranics

Park, Kyoung Chul↗

DOE OSTI2024

Experimental investigation of the small-scale prototype of the dew point water recovery system

The article presents experimental investigation of the Dew Point Water Recovery (DPWR) system designed for clean water recovery using the dew point cooling phenomenon. The technology is a novel distillation process which allows for highly efficient water recovery. It is a first distillation technology which allows for water distillation through direct contact with atmospheric air. The authors designed, built, and tested an original prototype of a laboratory-scale DPWR system to verify its basic operating principles and perform performance analysis. It was shown that on average 10 × 10 −6 m 3 of clean water could be recovered from each cubic meter of the supply air with an average GOR of 1.9. Although the operation of the DPWR system requires electricity to drive fans, pumps and other associated equipment, the actual driving force of the system is thermal energy. The average thermal energy demand of the DPWR system was 370 kWh/m 3 , and more than half of this energy was based on the thermodynamic potential of the supply air. Based on experimental data, the mathematical model of the DPWR system was verified and it was found that the average relative discrepancy in predicting air temperature did not exceed 8%, and in predicting air humidity: 5%. In conclusion, the conducted experiments confirmed that the proposed DPWR system is feasible and has the potential to fill the technological gap in terms of a low-energy and easily accessible water recovery method.

42 ENGINEERING↗

DOE OSTI2025

Kinetic Understanding of Field-Induced Phase Transition from Tetragonal to Ferroelectric Orthorhombic Phase in Ferroelectric CeO 2 –HfO 2 –ZrO 2 Films

The ferroelectric properties and structural phase transition behaviors of fluorite-type CeO 2 −HfO 2 −ZrO 2 films were investigated. The epitaxial films on indium tin oxide (ITO) (111)/yttria-stabilized zirconia (YSZ) (111) substrates were grown through pulsed laser deposition at room temperature and subsequently heat-treated at 1000 °C under a N 2 gas flow. The crystalline phases and Curie temperatures of the films were investigated by X-ray diffraction. An increase in the Ce or Zr content in the films led to a higher crystallographic symmetry, such as orthorhombic or tetragonal. In addition, electrical characterization revealed that the orthorhombic films and some of the tetragonal films displayed ferroelectricity. This was due to the field-induced phase transition from the tetragonal to ferroelectric orthorhombic phase in the films, where the Curie temperatures were relatively low. The tetragonal metastable phase was kinetically frozen and could not change into the stable orthorhombic phase at such a low temperature. The critical electric field where the field-induced phase transition occurred was below 0.8 MV/cm, which was sufficiently small compared to the coercive field. These results evidence the kinetic driving force that causes a field-induced phase transition from the paraelectric tetragonal phase to the ferroelectric orthorhombic phase in HfO 2 -based ferroelectrics. They also enhance our understanding of the thermodynamic phase stabilities of HfO 2 -based material polymorphs.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Microscopic origin of tunable assembly forces in chiral active environments

Across a variety of spatial scales, from nanoscale biological systems to micron-scale colloidal systems, equilibrium self-assembly is entirely dictated by—and therefore limited by—the thermodynamic properties of the constituent materials. In contrast, nonequilibrium materials, such as self-propelled active matter, expand the possibilities for driving the assemblies that are inaccessible in equilibrium conditions. Recently, a number of works have suggested that active matter drives or accelerates self-organization, but the emergent interactions that arise between solutes immersed in actively driven environments are complex and poorly understood. Here, we analyze and resolve two crucial questions concerning actively driven self-assembly: (i) how, mechanistically, do active environments drive self-assembly of passive solutes? (ii) Under which conditions is this assembly robust? We employ the framework of odd hydrodynamics to theoretically explain numerical and experimental observations that chiral active matter, i.e., particles driven with a directional torque, produces robust and long-ranged assembly forces. Overall, these developments constitute an important step towards a comprehensive theoretical framework for controlling self-assembly in nonequilibrium environments.

36 MATERIALS SCIENCE↗

DOE OSTI2025

The Observed and Projected Changes of Global Monsoons: Current Status and Future Perspectives

The global monsoon system, encompassing the Asian-Australian, African, and American monsoons, sustains two-thirds of the world’s population by regulating water resources and agriculture. Monsoon anomalies pose severe risks, including floods and droughts. Recent research associated with the implementation of the Global Monsoons Model Intercomparison Project under the umbrella of CMIP6 has advanced our understanding of its historical variability and driving mechanisms. Observational data reveal a 20th-century shift: increased rainfall pre-1950s, followed by aridification and partial recovery post-1980s, driven by both internal variability (e.g., Atlantic Multidecadal Oscillation) and external forcings (greenhouse gases, aerosols), while ENSO drives interannual variability through ocean-atmosphere interactions. Future projections under greenhouse forcing suggest long-term monsoon intensification, though regional disparities and model uncertainties persist. Models indicate robust trends but struggle to quantify extremes, where thermodynamic effects (warming-induced moisture rise) uniformly boost heavy rainfall, while dynamical shifts (circulation changes) create spatial heterogeneity. Volcanic eruptions and proposed solar radiation modification (SRM) further complicate predictions: tropical eruptions suppress monsoons, whereas high-latitude events alter cross-equatorial flows, highlighting unresolved feedbacks. The emergent constraint approach is booming in terms of correcting future projections and reducing uncertainty with respect to the global monsoons. Critical challenges remain. Model biases and sparse 20th-century observational data hinder accurate attribution. The interplay between natural variability and anthropogenic forcings, along with nonlinear extreme precipitation risks under warming, demands deeper mechanistic insights. Additionally, SRM’s regional impacts and hemispheric monsoon interactions require systematic evaluation. Addressing these gaps necessitates enhanced observational networks, refined climate models, and interdisciplinary efforts to disentangle multiscale drivers, ultimately improving resilience strategies for monsoon-dependent regions.

climate extreme events↗

DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗