Search NASA⌕ Search

SEARCH · Search NASA

Results for “binary fluids”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Impact of particle arrays on phase separation composition patterns

We examine the symmetry-breaking effect of fixed constellations of particles on the surface-directed spinodal decomposition of binary blends in the presence of particles whose surfaces have a preferential affinity for one of the components. Our phase-field simulations indicate that the phase separation morphology in the presence of particle arrays can be tuned to have a continuous, droplet, lamellar, or hybrid morphology depending on the interparticle spacing, blend composition, and time. In particular, when the interparticle spacing is large compared to the spinodal wavelength, a transient target pattern composed of alternate rings of preferred and non-preferred phases emerges at early times, tending to adopt the symmetry of the particle configuration. We reveal that such target patterns stabilize for certain characteristic length, time, and composition scales characteristic of the pure phase-separating mixture. To illustrate the general range of phenomena exhibited by mixture-particle systems, we simulate the effects of single-particle, multi-particle, and cluster–particle systems having multiple geometrical configurations of the particle characteristic of pattern substrates on phase separation. Our simulations show that tailoring the particle configuration, or substrate pattern configuration, a relative fluid-particle composition should allow the desirable control of the phase separation morphology as in block copolymer materials, but where the scales accessible to this approach of organizing phase-separated fluids usually are significantly larger. Limited experiments confirm the trends observed in our simulations, which should provide some guidance in engineering patterned blend and other mixtures of technological interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2003

Simulating the Effect of Space Vehicle Environments on Directional Solidification of a Binary Alloy

Space microgravity missions are designed to provide a microgravity environment for scientific experiments, but these missions cannot provide a perfect environment, due to vibrations caused by crew activity, on-board experiments, support systems (pumps, fans, etc.), periodic orbital maneuvers, and water dumps. Therefore, it is necessary to predict the impact of these vibrations on space experiments, prior to performing them. Simulations were conducted to study the effect of the vibrations on the directional solidification of a dendritic alloy. Finite element ca!cu!attie?ls were dme with a simd2titcr based on a continuum model of dendritic solidification, using the Fractional Step Method (FSM). The FSM splits the solution of the momentum equation into two steps: the viscous intermediate step, which does not enforce continuity; and the inviscid projection step, which calculates the pressure and enforces continuity. The FSM provides significant computational benefits for predicting flows in a directionally solidified alloy, compared to other methods presently employed, because of the efficiency gains in the uncoupled solution of velocity and pressure. finite differences, arises when the interdendritic liquid reaches the eutectic temperature and concentration. When a node reaches eutectic temperature, it is assumed that the solidification of the eutectic liquid continues at constant temperature until all the eutectic is solidified. With this approach, solidification is not achieved continuously across an element; rather, the element is not considered solidified until the eutectic isotherm overtakes the top nodes. For microgravity simulations, where the convection is driven by shrinkage, it introduces large variations in the fluid velocity. When the eutectic isotherm reaches a node, all the eutectic must be solidified in a short period, causing an abrupt increase in velocity. To overcome this difficulty, we employed a scheme to numerically predict a more accurate value for the rate of change of fraction of liquid as the liquid in an element solidifies. The new method enables us to contrast results of simulations in which the alloy is subjected to no gravity or a steady-state acceleration versus simulations when the alloy is subjected to vibration disturbances; therefore, the effect of vibration disturbances can be assessed more accurately. To assess the impact of these vibration-perturbations, transient accelerometer data from a space shuttle mission are used as inputs for the simulation model. These on-orbit acceleration data were obtained from the Microgravity Science Division at Glenn Research Center (GRC- MSD) and are applied to the buoyancy term of the momentum equation in a simulation of a Pb-5.8 wt. % Sb alloy that solidifies in a thermal gradient of 4000 K/m and a translation velocity of 3 p d s . Figure 2 shows the vertical velocity of a node that begins in the all-liquid region and subsequently solidifies; the vibrations are applied at 5000 seconds in this simulation. An important difficulty, common to all solidification models based on finite elements or 2 The magnitudes of the velocity oscillations that are vibration-induced are very small and acceptable. The biggest concern is whether the concentration of the liquid near the dendrite tips is distorted because of the vibration-induced perturbations. Results for this case show no concentration oscillations present in the all-liquid region.

Westra, D. G.↗
DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗
NASA NTRS2001

Vapor Growth of Binary and Ternary Chalcogenides in Preparation for Microgravity Experiments

In the bulk crystal growth of some technologically important semiconducting chalcopyrites, such as ZnTe, CdS, ZnSe and ZnS, vapor growth techniques have significant advantages over melt growth techniques due to the high melting points of these materials. The realization of routine production of high-quality single crystals of these semiconductors requires a fundamental, systematic and in-depth study on the PVT growth process and crystal growth by vapor transport in low gravity offers a set of unique conditions for this study. Previously, two reasons have been put forward to account for this. The first is weight-related reductions in crystal strain and defects. These are thought to be caused by the weight of the crystals during processing at elevated temperatures and retained on cooling, particularly for materials with a low yield strength. The second, and more general, reason is related to the reduction in density-gradient driven convection. The PVT crystal growth process consists of essentially three processes: sublimation of the source material, transport of the vapor species and condensation of the vapor species to form the crystal. The latter two processes can be affected by the convection caused by gravitational accelerations on Earth. Reductions in such convection in low gravity is expected to yield a nearly diffusion-limited growth condition which results in more uniform growth rates (on the microscopic scale) and hence greater crystalline perfection and compositional homogeneity. The reduction of convective contamination by performing flight experiments in a reduced gravity environment will help to understand the relation between fluid phase processes (growth parameters) and defect and impurity incorporation in grown crystals.

Su, C.↗
NASA NTRS2001

Vapor Growth of Binary and Ternary Chalcogenides in Preparation for Microgravity Experiments

In the bulk crystal growth of some technologically important semiconducting chalcopyrites, such as ZnTe, CdS, ZnSe and ZnS, vapor growth techniques have significant advantages over melt growth techniques due to the high melting points of these materials. The realization of routine production of high-quality single crystals of these semiconductors requires a fundamental, systematic and in-depth study on the PVT growth process and crystal growth by vapor transport in low gravity offers a set of unique conditions for this study. Previously, two reasons have been put forward to account for this. The first is weight-related reductions in crystal strain and defects. These are thought to be caused by the weight of the crystals during processing at elevated temperatures and retained on cooling, particularly for materials with a low yield strength. The second, and more general, reason is related to the reduction in density-gradient driven convection. The PVT crystal growth process consists of essentially three processes: sublimation of the source material, transport of the vapor species and condensation of the vapor species to form the crystal. The latter two processes can be affected by the convection caused by gravitational accelerations on Earth. Reductions in such convection in low gravity is expected to yield a nearly diffusion-limited growth condition which results in more uniform growth rates (on the microscopic scale) and hence greater crystalline perfection and compositional homogeneity. The reduction of convective contamination by performing flight experiments in a reduced gravity environment will help to understand the relation between fluid phase processes (growth parameters) and defect and impurity incorporation in grown crystals.

Su, Ching-Hua↗
NASA NTRS2013

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant↗
DOE OSTI2020

Adsorption of furan, hexanoic acid, and alkanes in a hierarchical zeolite at reaction conditions: Insights from molecular simulations

Hierarchical zeolites containing both micropores and mesopores are valuable catalysts for facilitating reactions of large molecules. Furan acylation by fatty acids is a promising reaction for valorizing biomass, and the self-pillared pentasil (SPP) zeolite was found to perform particularly well for this reaction. To better understand the distribution of molecules in hierarchical zeolites at the elevated temperature (T = 523 K) and the elevated pressure (p >1 bar) associated with typical reaction conditions, unary and binary adsorption were predicted using Monte Carlo simulations in the isothermal–isobaric Gibbs ensemble. Adsorption of six species (furan, hexanoic acid, -hexane, -decane, -tetradecane, and 3,6-diethyloctane) was investigated from vapor, liquid, and supercritical phases, and loadings into the micropores, onto the mesopore surface, and in the mesopore interior of SPP were obtained. As pressure increases, -alkanes fill the micropores before loading the surface and then the interior of the mesopore, while furan and hexanoic acid adsorb strongly to the mesopore surface due to hydrogen bonding interactions with surface silanols. Hydrogen bonding interactions also draw hexanoic acid molecules in the micropore region toward the pore mouths, so their carboxylic acid group forms H-bonds with silanols, while the alkyl tails interact with the micropore walls. Mesopore condensation is observed for molecules below their critical point, and occurs when the Gibbs free energy of transfer into the mesopore interior and onto the mesopore surface converge. When hexanoic acid adsorption occurs in the presence of alkane solvents, then the selectivity and spatial distribution of hexanoic acid in the micropores and on the surface can be tuned by adjusting the fluid pressure and the alkane length and/or branching.

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

Mechanochemical Synthesis of High-Purity Anhydrous Binary Alkali and Alkaline Earth Chloride Mixtures

A direct synthesis route for high-purity, anhydrous binary salt mixtures has been developed. This atom efficient, solvent-free process is easily scalable, with the potential to produce salt mixtures that meet the purity standards required for industrial heat transfer and nuclear applications. The essence of the methodology lies in mechanochemical synthesis of carnallite precursors that can mitigate the hydrolysis of MgCl 2 ·6H 2 O under direct heating. Each dehydrated salt carnallite was then analyzed for purity and oxide content through subsequent powder X-ray diffraction, and strong acid titration. In conclusion, this process presents a more effective alternative route compared to previous methods for obtaining low-oxide, high-purity chloride salt mixtures.

Alkali↗
NASA NTRS2020

An Evolutionary System of Mineralogy. Part I: Stellar Mineralogy (>13 to 4.6 Ga)

Minerals preserve records of the physical, chemical, and biological histories of their origins and subsequent alteration, and thus provide a vivid narrative of the evolution of Earth and other worlds through billions of years of cosmic history. Mineral properties, including trace and minor elements, ratios of isotopes, solid and fluid inclusions, external morphologies, and other idiosyncratic attributes, represent information that points to specific modes of formation and subsequent environmental histories—information essential to understanding the co-evolving geosphere and biosphere. This perspective suggests an opportunity to amplify the existing system of mineral classification, by which minerals are defined solely on idealized end-member chemical compositions and crystal structures. Here we present the first in a series of contributions to explore a complementary evolutionary system of mineralogy—a classification scheme that links mineral species to their paragenetic modes. The earliest stage of mineral evolution commenced with the appearance of the first crystals in the universe at >13 Ga and continues today in the expanding, cooling atmospheres of countless evolved stars, which host the high-temperature (T > 1000 K), low-pressure (P < 10-2 atm) condensation of refractory minerals and amorphous phases. Most stardust is thought to originate in three distinct processes in carbon- and/or oxygen-rich mineral-forming stars: (1) condensation in the cooling, expanding atmospheres of asymptotic giant branch stars; (2) during the catastrophic explosions of supernovae, most commonly core collapse (Type II) supernovae; and (3) classical novae explosions, the consequence of runaway fusion reactions at the surface of a binary white dwarf star. Each stellar environment imparts distinctive isotopic and trace element signatures to the micro- and nanoscale stardust grains that are recovered from meteorites and micrometeorites collected on Earth’s surface, by atmospheric sampling, and from asteroids and comets. Although our understanding of the diverse mineral-forming environments of stars is as yet incomplete, we present a preliminary catalog of 41 distinct natural kinds of stellar minerals, representing 22 official International Mineralogical Association (IMA) mineral species, as well as 2 as yet unapproved crystalline phases and 3 kinds of non-crystalline condensed phases not codified by the IMA.

Robert M. Hazen↗
DOE OSTI2023

Accurate equation of state of H 2 – He binary mixtures up to 5.4 GPa

Brillouin scattering spectroscopy has been used to obtain an accurate (<1%) ρ-P equation of state (EOS) of 1:1 and 9:1 H 2 -He molar mixtures from 0.5 to 5.4 GPa at 296 K. Our calculated equations of state indicate close agreement with the experimental data right to the freezing pressure of hydrogen at 5.4 GPa. The measured velocities agree on average, within 0.5%, of an ideal mixing model. The ρ-P EOSs presented have a standard deviation of under 0.3% from the measured densities and under 1% deviation from ideal mixing. Furthermore, a detailed discussion of the accuracy, precision, and sources of error in the measurement and analyses of our equations of state is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗
NASA NTRS1988

Critical fluid light scattering

The objective is to measure the decay rates of critical density fluctuations in a simple fluid (xenon) very near its liquid-vapor critical point using laser light scattering and photon correlation spectroscopy. Such experiments were severely limited on Earth by the presence of gravity which causes large density gradients in the sample when the compressibility diverges approaching the critical point. The goal is to measure fluctuation decay rates at least two decades closer to the critical point than is possible on earth, with a resolution of 3 microK. This will require loading the sample to 0.1 percent of the critical density and taking data as close as 100 microK to the critical temperature. The minimum mission time of 100 hours will allow a complete range of temperature points to be covered, limited by the thermal response of the sample. Other technical problems have to be addressed such as multiple scattering and the effect of wetting layers. The experiment entails measurement of the scattering intensity fluctuation decay rate at two angles for each temperature and simultaneously recording the scattering intensities and sample turbidity (from the transmission). The analyzed intensity and turbidity data gives the correlation length at each temperature and locates the critical temperature. The fluctuation decay rate data from these measurements will provide a severe test of the generalized hydrodynamic theories of transport coefficients in the critical regions. When compared to equivalent data from binary liquid critical mixtures they will test the universality of critical dynamics.

Gammon, Robert W.↗
NASA NTRS2017

Vascular Patterning Analysis by VESGEN 2D/3D with Bioinformatics: Updates for Rodent Tissues

Fractally branching vascular systems are a complex physiological requirement shared by humans with all higher terrestrial life forms, including other vertebrates, insects, and higher land plants. Vascular trees, networks, and tree-network composites are therefore mapped and quantified by the VESsel GENeration Analysis (VESGEN) software according to weighted physiological vascular rules that include vessel connectivity, tapering and bifurcational branching. According to fluid dynamics, successful vascular transport depends upon a complex distributed system of highly regulated laminar flow. VESGEN has elucidated changes in vascular patterning resulting from inflammatory, developmental and other signaling pathways within numerous tissues of major model organisms important for Space Biology, especially for rodents. Important early stage regenerative opportunities have been identified by VESGEN vascular analysis for visual impairments in the human retina, and is currently being used for research into astronaut visual and ocular disorders associated with long duration missions. The VESGEN 2D software is a mature, automated, widely published capability for which beta testing and public release by NASA is planned for the upcoming year. Early-stage capabilities for VESGEN 3D analysis are under development for the rodent retina and intestine as prototype tissues. A prototype VESGEN 2D Bioinformatics software capability has also been developed to associate phenotypic changes in molecular expression with vascular structure and function. By new VESGEN bioinformatic innovations, expression patterns of the genetic, transcriptional, protein and other markers for regulatory molecules such as vascular endothelial growth factor (VEGF) and their receptors, often indicators of tissue oxygenation status, are co-localized with alterations in vascular pattern. Biomarkers are therefore mapped and quantified as information dimensions directly correlated with the spatial dimensions of a vascular pattern. Further important technology innovations by NASA include substantial image segmentation advances for more automated binary extraction of the grayscale vascular patterns, together with informative associated image quality assessments. Vascular mapping and quantification capabilities for the rodent retina and intestine are illustrated for VESGEN 2D, along with technology status reports on VESGEN 3D and Bioinformatic capabilities. Research partially supported by Ames Center Innovation Awards.

Parsons-Wingerter, P.↗
DOE OSTI2020

Operational Modes of a 2.0 MWth Chloride Molten-Salt Pilot-Scale System

The limit of traditional solar-salt thermal stability is around 600 °C with ambient air as the cover gas. Nitrate molten salt concentrating solar power (CSP) systems are currently deployed globally and are considered to be state-of the art heat transfer fluids (HTFs) for present day high-temperature operation. However, decomposition challenges occur with these salts for operation beyond 600 C. Although slightly higher limits may be possible with solar salt, to fully realize SunShot efficiency goals of $15/kWhth HTFs and an LCOE of 6¢/kWh, molten-salt technologies working at higher temperatures (e.g., 650 °C to 750 °C) will require an alternative salt chemistry composition, such as chlorides. In this investigation a 2.0MWth Pilot-scale CSP plant design is developed to assess thermodynamic performance potential for operation up to 720 C. Here, an Engineering Equation Solver (EES) model is developed with respect to 14 state-points from the base of a solar tower at the Sandia National Laboratories, National Solar Thermal Test Facility (NSTTF), to solar receiver mounted 120 ft. above the ground. The system design considers a ternary chloride ternary chloride (20%NaCl/40%MgCl/40%KCl by mol%) salt as the HTF, with 6 hrs. of storage and a 1 MWth primary salt to sCO2 heat exchanger. Preliminary system modelling results indicate a minimum non-dimensional Cv of 60 required for both cold and hot-side throttle recirculation valves for the operational pump operating between speeds of 1800 and 2400 RPM. Further receiver comparison study results suggest that the ternary salt requires an average 15.2% higher receiver flux with a slightly lower calculated receiver efficiency when compared to a binary carnelite salt to achieve a 2.0 MWth desired input power design.

associated liquids↗
DOE OSTI2023

Vapor and Liquid ( p –ρ– T–x ) Measurements of Binary Refrigerant Blends Containing R-32, R-152a, R-227ea, R-1234yf, and R-1234ze(E)

In this article, The pressure–density–temperature–composition (p–ρ–T–x) data of binary refrigerant mixtures containing R-32 (difluoromethane), R-152a (1,1,-difluoroethane), R-227ea (1,1,1,2,3,3,3-heptafluoropropane), R-1234yf (2,3,3,3-tetrafluoropropene), and R-1234ze(E) (trans-1,3,3,3-tetrafluoropropene) were measured in both the vapor and liquid phases using a two-sinker, magnetic suspension densimeter. The specific samples in this study comprised two compositions of approximately (0.3/0.7) and (0.7/0.3) mole fraction for each of the following four binary refrigerant blends: R-32 + R-1234yf, R-32 + R-1234ze(E), R-1234yf + R-152a, and R-1234ze(E) + R-227ea. Single-phase vapor densities were measured over a temperature range of approximately 253 to 293 K and pressures from 0.05 to 0.98 MPa. Single-phase liquid and supercritical densities were measured over a temperature range of approximately 230 to 400 K and pressures up to 22 MPa; for refrigerant blends containing R-1234yf, the maximum pressure was limited to 14 MPa. Overall relative combined, expanded (k = 2) uncertainties in density ranged from 0.025 to 0.191%, with an average uncertainty of approximately 0.05%. Here, we present measurement results, along with comparisons to available literature data and to default equations of state and mixture models included in REFPROP.

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

Development of a Gibbs Energy Minimiser for the MOOSE-based Corrosion Modelling App Yellowjacket and Validation of MSTDB

Nuclear materials are highly complex multiscale, multiphysics systems,and an effective prediction of nuclear reactor performance and safety requires simulation capabilities that tightly couple different physical phenomena. The Idaho National Laboratory’s Multiphysics Object Oriented Simulation Environment (MOOSE) provides the computational foundation for performing such simulations. With the move towards advanced reactors, such as the Molten Salt Reactor (MSR), that employ high temperature fluids compared to conventional reactors, corrosion has become a problem of great interest. A new application called Yellowjacket is currently under development to directly couple thermodynamic equilibrium and kinetics with phase field models in order to model corrosion in MSRs. As part of Yellowjacket, a Gibbs energy minimiser is being developed to perform thermochemical equilibrium calculations for a range of different materials, which is currently in its infancy. This report describes the further progress towards the development of Yellowjacket Gibbs energy minimiser. Ontario Tech University is developing a new Gibbs energy minimiser for Yellowjacket which is the primary contribution of this work. The aim to develop a thermochemistry solver for the MOOSE framework following the same development philosophy and using the same tools and libraries. A special focus is on performance, documentation and SQA. Furthermore through a scope extension partway through the fiscal year a thorough assessment of the MSTDB-TC v1.3 was performed at Ontario Tech University in the context of continuous improvement and quality assurance. The objective of the work was to have an arm’s length review of the database to give confidence that the database is performing as it was intended while assessing its current state to give recommendations to future developments. This assessment involved two parts: A) a quantitative assessment, and B) a qualitative assessment. Part A involved developing an automated test-suite that would compute values from the database using Thermochimica with comparisons to experimental measurements for validation purposes, which gives confidence to the database’s stakeholders that its functioning properly. Part B involved reviewing all binary systems in the database and making a qualitative assessment with two performance indicators: comprehensiveness and overall confidence. It is important to note that the models in the database are empirical, which is to say that the quality of any model is highly dependent on the experimental data used to inform its development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1994

Low-frequency modes and nonbarotropic effects in pseudo-Newtonian accretion disks

A recently developed formalism is used to reexamine the question of the existence of hydrodynamical modes that pulsate with very low frequencies in the inner regions of accretion disks. The formalism is valid in an exact sense for the adiabatic pulsations of rotating Newtonian fluids that are generally nonbarotropic (such as those with 'nonadiabatic temperature gradients,' for example), and hence its application in the present context represents an improvement over previous analyses that are more approximate. The formalism is applied to thin non-self-gravitating disks, with the gravitational potential of the central source modified in the usual way in order to simulate relativistic effects. In the barotropic limit, the analyses indicate that in many cases nearly Keplerian disks exhibit nonaxisymmetric modes of pulsation that are trapped in the inner disk regions, with pulsation periods much longer than the dynamical timescale. These results are similar to those of earlier calculations that assume disks pulsate without changing the temperature distribution. A method is developed for including lowest order nonbarotropic effects. Previous analyses have been incapable of accurately treating the nonbarotropic regime. The application of the present method to the low-frequency modes reveals that, due to unexpected cancellations among terms, the nonbarotropic correction to the pusation frequency omega is only of order tilde-omega(sub BV exp 2) omega, where tilde-omega(sub BV) is the appropriate dimensionless Brunt-Vaisala frequency. This correction is much smaller than the expected correction of order tilde-omega(sub BV) Omega, where Omega is the rotation angular velocity. The important conclusion drawn from this is that nonbarotropic corrections are generally small and hence that low-frequency modes persist into the nonbarotropic regime. For disk temperatures appropriate to X-ray emission, the adiabatic frequencies of trapped modes are of the same order as the frequencies associated with quasi-periodic oscillations (QPOs) of certain low-mass X-ray binaries and other X-ray sources. It is suggested that further consideration be given to the possiblity that modulation of emission processes by these modes can explain at least some of the QPO phenomena.

Ipser, James R.↗
NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗