SEARCH · Search NASA

Results for “Liquid structure simulation”

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 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Planar imaging of the mixing characteristics of a multi-element triplet injector

Planar light scattering has been used to investigate the turbulent mixing in a multi-element triplet injector at mixture ratio conditions typical of liquid rocket engine gas generators. Experiments were conducted at atmospheric pressure using air seeded with smoke particles to simulate the oxidizer flow and unseeded air to simulate the fuel. Data sets of 200 images (containing 122,000 simultaneous measurements each) of the instantaneous distribution of the scattered light were acquired at several positions throughout the flow. These images have been processed digitally to produce local distributions of the average and rms fluctuation of the mixture fraction. In addition, the spatial autocorrelation of the images has been calculated and analyzed to determine a mean length scale related to the large scale structures in the flow.

Makel, D. B.↗

NASA NTRS2009

Explosion/Blast Dynamics for Constellation Launch Vehicles Assessment

An assessment methodology is developed to guide quantitative predictions of adverse physical environments and the subsequent effects on the Ares-1 crew launch vehicle associated with the loss of containment of cryogenic liquid propellants from the upper stage during ascent. Development of the methodology is led by a team at Sandia National Laboratories (SNL) with guidance and support from a number of National Aeronautics and Space Administration (NASA) personnel. The methodology is based on the current Ares-1 design and feasible accident scenarios. These scenarios address containment failure from debris impact or structural response to pressure or blast loading from an external source. Once containment is breached, the envisioned assessment methodology includes predictions for the sequence of physical processes stemming from cryogenic tank failure. The investigative techniques, analysis paths, and numerical simulations that comprise the proposed methodology are summarized and appropriate simulation software is identified in this report.

Baer, Mel↗

DOE OSTI2022

Ion Transport in 2D Nanostructured $\pi$-Conjugated Thieno[3,2- b ]thiophene-Based Liquid Crystal

Leveraging the self-assembling behavior of liquid crystals designed for controlling ion transport is of both fundamental and technological significance. Here, we have designed and prepared a liquid crystal that contains (2,5-bis(5-(2,5,8,11-tetraoxatridecan-13-yl)thiophen-2yl)thieno[3,2-b]thiophene (BTTT) as mesogenic core and conjugated segment and symmetric tetra(ethylene oxide) (EO4) as polar side chains for ionic conducting regions. Driven by the crystallization of the BTTT cores, BTTT/dEO4 exhibits well-ordered smectic phases below 71.5 °C as confirmed by differential scanning calorimetry, polarized optical microscope, temperature dependent wide-angle X-ray scattering and grazing incidence wide-angle X-ray scattering (GIWAXS). We adopted a combination of experimental GIWAXS and molecular dynamics (MD) simulations to better understand the molecular packing of BTTT/dEO4 films, particularly when loaded with the ion conducting salt, lithium bis(trifluoromethanesulfonyl)imide (LiTFSI). Ionic conduction of BTTT/dEO4 is realized by the addition of LiTFSI, with the material able to maintain smectic phases up to r = [Li + ]/[EO] = 0.1. In this work, The highest ionic conductivity of 8 × 10 -3 S/cm was attained at an intermedium salt concentration of r = 0.05. It was also found that ion conduction in BTTT/dEO4 is enhanced by forming a smectic layered structure with irregular interfaces between the BTTT and EO4 layers and by the lateral film expansion upon salt addition. This can be explained by the enhancement of the misalignment and configurational entropy of the side chains, which increase their local mobility and that of the solvated ions. Our molecular design thus illustrates how, beyond the favorable energetic interactions that drive the assembly and ion solvating domains, modulation of entropic effects can also be favorably harnessed to improve ion conduction.

thin film electrochemical impedance spectroscopy↗

NASA NTRS2008

Active Co-Storage of Cryogenic Propellants for Lunar Explortation

Long-term storage of cryogenic propellants is a critical requirement for NASA's effort to return to the moon. Liquid hydrogen and liquid oxygen provide the highest specific impulse of any practical chemical propulsion system, and thus provides the greatest payload mass per unit of launch mass. Future manned missions will require vehicles with the flexibility to remain in orbit for months, necessitating long-term storage of these cryogenic liquids. For decades cryogenic scientific satellites have used dual cryogens with different temperatures to cool instruments. This technology utilizes a higher temperature cryogen to provide a stage that efficiently intercepts a large fraction of the heat that would otherwise be incident on the lower temperature cryogen. This interception reduces the boil-off of the lower temperature cryogen and increasing the overall life-time of the mission. The Active Co-Storage concept is implemented similarly; the 101 K liquid oxygen thermally shields the 24 K liquid hydrogen. A thermal radiation shield that is linked to the liquid oxygen tank shrouds the liquid hydrogen tank, thereby preventing the liquid hydrogen tank from being directly exposed to the 300 K external environment. Modern cryocooler technology can eliminate the liquid oxygen boil-off and also cool the thermal radiation shield thereby reducing the liquid hydrogen boil-off to a small fraction of the unshielded rate. The thermal radiation shield can be a simple conductive shroud or a more sophisticated but lighter Broad Area Cooling (BAC) shroud. The paper describes the design impact of an active co-storage system for the Altair Descent Vehicle. This paper also compares the spacecraft-level impacts of the conductive shroud and the BAC shroud active co-storage concepts with a passive storage option in the context of the different scales of spacecraft that will be used for the lunar exploration effort - the Altair Ascent and Descent Vehicles, the Orion, and the Ares V Earth Departure Stage. The paper also reports on a subscale test of this active co-storage configuration. The test tank is 0.7 m in diameter, approximately one-third the dimension of tanks that would be needed in a lunar ascent module. A thin-walled fiberglass skirt supports and isolates the tank from a 100 K stage. A similar thin-walled skirt supports the lOOK stage from the ambient temperature structure. An aluminum shield with a heavy MLI blanket surrounds the tank and is attached at the 100 K stage. In this initial phase of the project, there is no tank on the 100 K stage, but it is actively cooled by a single-stage cryocooler similar in design to the one used on the RHESSI mission. The test configuration includes a number of innovative elements, including a helical support heat exchanger and an external thermodynamic vent/heat interception system. To avoid the complexity of an explosive gas handling system, testing will be done with liquid helium and liquid neon as simulant fluids. The properties of these fluids bracket the properties of liquid hydrogen. Instrumentation allows tank temperature and shield temperature profiles, tank liquid levels, and pressure drops through the flow lines, to be measured.

Mustafi, S.↗

NASA NTRSDate not supplied

Coupled Reactor Multiphysics and Mass Scalability Assessment for Crewed Megawatt-Class NEP System Architectures

Nuclear Electric Propulsion (NEP) is an in-space propulsion technology capable of enabling opposition and conjunction class crewed Mars missions. NEP subsystems include the reactor for heat generation, a power conversion system (PCS), power management and distribution, electric propulsion system, and heat rejection system. Specific mass, or α (kg/kWe), is a key performance parameter (KPP) of the propulsion system which is directly scalable with the performance and mass predictions for each individual component. To inform technology maturation planning activities, full system and component level parametric modeling is ongoing to explore the design trade space and illustrate the effect of subsystem design choices on the system KPPs. In this study, scaling of high-assay, low-enriched uranium reactor designs is assessed through coupled reactor physics and thermal hydraulics analyses. Scaling analyses evaluate the impact of system performance parameters (power level, interface temperatures) on mass for direct gas cooled, pumped liquid metal, and passively-cooled heat pipe reactor concepts. Each concept requires specific geometries, fluids, and power conversion interface conditions (temperature, pressure, flow rate) to meet desired performance and mass. The reactor assembly includes the active core (fuel, moderator, cladding, working fluid), axial and radial neutron reflectors, control drums, structural support / pressure vessel, and external radiation shielding. Each of these components are parametrically sized based on performance parameters for a megawatt-class power cycle. Results of this scaling analysis increase NEP propulsion system modeling fidelity and ultimately aim to support concept down-selection along with related technology development planning. The reactor and shield α are a function of several PCS and heat rejection system design choices, and reactor scaling with these parameters must be considered to enable an informed decision on an optimal reactor geometry and working fluid combination.

Nuclear Electric Propulsion↗

DOE OSTI2025

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Membrane hydrophobicity determines the activation free energy of passive lipid transport

The collective behavior of lipids with diverse chemical and physical features determines a membrane's thermodynamic properties. Yet, the influence of lipid physicochemical properties on lipid dynamics, in particular interbilayer transport, remains underexplored. Here, we systematically investigate how the activation free energy of passive lipid transport depends on lipid chemistry and membrane phase. Through all-atom molecular dynamics simulations of 11 chemically distinct glycerophospholipids, we determine how lipid acyl chain length, unsaturation, and headgroup influence the free energy barriers for two elementary steps of lipid transport: lipid desorption, which is rate limiting, and lipid insertion into a membrane. Consistent with previous experimental measurements, we find that lipids with longer, saturated acyl chains have increased activation free energies compared to lipids with shorter, unsaturated chains. Lipids with different headgroups exhibit a range of activation free energies; however, no clear trend based solely on chemical structure can be identified, mirroring difficulties in the interpretation of previous experimental results. Compared to liquid-crystalline phase membranes, gel phase membranes exhibit substantially increased free energy barriers. Overall, we find that the activation free energy depends on a lipid's local hydrophobic environment in a membrane and that the free energy barrier for lipid insertion depends on a membrane's interfacial hydrophobicity. Both of these properties can be altered through changes in lipid acyl chain length, lipid headgroup, and membrane phase. Thus, the rate of lipid transport can be tuned through subtle changes in local membrane composition and order, suggesting an unappreciated role for nanoscale membrane domains in regulating cellular lipid dynamics.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Structural evolution of liquid silicates under conditions in Super-Earth interiors

Molten silicates at depth are crucial for planetary evolution, yet their local structure and physical properties under extreme conditions remain elusive due to experimental challenges. In this study, we utilize in situ X-ray diffraction (XRD) at the Matter in Extreme Conditions (MEC) end-station of the Linear Coherent Linac Source (LCLS) at SLAC National Accelerator Laboratory to investigate liquid silicates. Using an ultrabright X-ray source and a high-power optical laser, we probed the local atomic arrangement of shock-compressed liquid (Mg,Fe)SiO 3 with varying Fe content, at pressures from 81(9) to 385(40) GPa. We compared these findings to ab initio molecular dynamics simulations under similar conditions. Results indicate continuous densification of theO-O and Mg-Si networks beyond Earth’s interior pressure range, potentially altering melt properties at extreme conditions. This could have significant implications for early planetary evolution, leading to notable differences in differentiation processes between smaller rocky planets, such as Earth and Venus, and super-Earths, which are exoplanets withmasses nearly three times that of Earth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Experimental determination of the temperature-dependent Van Hove function in a Zr 80 Pt 20 liquid

Even though the viscosity is one of the most fundamental properties of liquids, the connection with the atomic structure of the liquid has proven elusive. By combining inelastic neutron scattering with the electrostatic levitation technique, the time-dependent pair-distribution function (i.e., the Van Hove function) has been determined for liquid Zr 80 Pt 20 . We show that the decay time of the first peak of the Van Hove function is directly related to the Maxwell relaxation time of the liquid, which is proportional to the shear viscosity. This result demonstrates that the local dynamics for increasing or decreasing the coordination number of local clusters by one determines the viscosity at high temperature, supporting earlier predictions from molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Temperature dependence of structural, dynamical, and electronic properties of amorphous Bi 2 Te 3 : an ab initio study

Bismuth telluride (Bi 2 Te 3 ) has garnered significant interest in thermoelectric applications and three-dimensional topological insulators due to its unique electronic, transport, and thermal properties. Bi 2 Te 3 and Sb 2 Te 3 chalcogenide compounds have the same crystal structure. While Sb 2 Te 3 has been shown to be a prototypical phase change memory (PCM) compound along the pseudobinary tie-line of Ge-Sb-Te alloys, whether Bi 2 Te 3 can also exhibit PCM functionality is still not well established. In this work, a systematic study on the structural, dynamical, and electronic properties of amorphous Bi 2 Te 3 during the quenching process has been performed by using ab initio molecular dynamics simulations. Pair correlation function, coordination number, bond-angle distribution functions, and a novel atomistic cluster alignment method are used to explore the structural characteristics of Bi 2 Te 3 as a function of temperature. Our study shows that there are many distorted octahedral clusters in amorphous Bi 2 Te 3 . In comparison with the local structures in Sb 2 Te 3 , we found that the degree of distortion of the octahedrons in the Bi 2 Te 3 system is smaller than that in Sb 2 Te 3 system. Moreover, the changes in the dynamical properties of Bi 2 Te 3 from liquid to glassy state are also explored. The approximate range of liquid-to-glass transition temperature is determined to be between 673 and 723 K. The electronic properties of Bi2Te 3 and Sb 2 Te 3 are also analysed by density-of-states and Bader charge calculations, both of them in glass state are semiconductors. Our studies provide useful insights into the local structure and dynamical properties of Bi 2 Te 3 at the atomistic level during the fast cooling process, and suggest that the compound can be a candidate for PCM materials.

36 MATERIALS SCIENCE↗

NASA NTRS1999

Solidification of II-VI Compounds in a Rotating Magnetic Field

This project is aimed at using a rotating magnetic field (RMF) to control fluid flow and transport during directional solidification of elemental and compound melts. Microgravity experiments have demonstrated that small amounts of residual acceleration of less than a micro-g can initiate and prolong fluid flow, particularly when there is a static component of the field perpendicular to the liquid solid interface. Thus a true diffusion boundary layer is not formed, and it becomes difficult to verify theories of solidification or to achieve diffusion controlled solidification. The RMF superimposes a stirring effect on an electrically conducting liquid, and with appropriate field strengths and frequencies, controlled transport of material through a liquid column can be obtained. As diffusion conditions are precluded and complete mixing conditions prevail, the technique is appropriate for traveling solvent zone or float zone growth methods in which the overall composition of the liquid can be maintained throughout the growth experiment. Crystals grown by RMF techniques in microgravity in previous, unrelated missions have shown exceptional properties. The objective of the project is two-fold, namely (1) using numerical modeling to simulate the behavior of a solvent zone with applied thermal boundary conditions and demonstrate the effects of decreasing gravity levels, or an increasing applied RMF, or both, and (2) to grow elements and II-VI compounds from traveling solvent zones both with and without applied RMFs, and to determine objectively how well the modeling predicts solidification parameters. Numerical modeling has demonstrated that, in the growth of CdTe from a tellurium solution, a rotating magnetic field can advantageously modify the shape of the liquid solid interface such that the interface is convex as seen from the liquid. Under such circumstances, the defect structure is reduced as any defects which are formed tend to grow out and not propagate. The flow of liquid, however, is complex due to the competing flow induced by the rotating magnetic field and the buoyancy driven convection. When the acceleration forces are reduced to one thousandth of gravity, the flow pattern is much simplified and well controlled material transport through the solvent zone can be readily achieved. Triple axis diffractometry and x-ray synchrotron topography have demonstrated that there is no significant improvement in crystal quality for HgCdTe grown on earth from a tellurium solution when a rotating magnetic field is applied. However, modeling shows that the flow in microgravity with a rotating magnetic field would produce a superior product.

Gillies, D. C.↗

NASA NTRSDate not supplied

Lightweight Thrust Chamber Assemblies using Multi-Alloy Additive Manufacturing and Composite Overwrap

Additive Manufacturing (AM) has brought significant design and fabrication opportunities for complex components with internal features such as liquid rocket engine thrust chambers not previously possible. This technology allows for significant cost savings and schedule reductions in addition to new performance optimization through weight reduction and increased margins. Specific to regeneratively-cooled combustion chambers and nozzles for liquid rocket engines, additive manufacturing offers the ability to form the complex internal coolant channels and the closeout of the channels to contain the high pressure liquid propellants with a single operation. Much of additive manufacturing development has focused on monolithic alloys using Laser Powder Bed Fusion (L-PBF), which do not allow for complete optimization of the structure. The National Aeronautics and Space Administration (NASA) completed feasibility of an AM bimetallic L-PBF GRCop-84 copper-alloy combustion chamber with an AM electron beam freeform Inconel 625 structural jacket under the Low Cost Upper Stage Propulsion (LCUSP) Project. A follow-on project called Rapid Analysis and Manufacturing Propulsion Technology (RAMPT) is under development to further expand large-scale multi-alloy thrust chambers while maturing composite overwrap technology for significant weight savings opportunities. The RAMPT project has three primary objectives: 1) Advancing blown powder Directed Energy Deposition (DED) to fabricate integral-channel large scale nozzles, 2) Develop composite overwrap technology to reduce weight and provide structural capability for thrust chamber assemblies, and 3) Develop bimetallic and multi-metallic additively manufactured radial and axial joints to optimize material performance. In addition to these primary manufacturing developments, analytical modeling efforts compliment the process development to simulate the AM processes to reduce build failures and distortions. The RAMPT project is also maturing the supply chain for various manufacturing processes described above in addition to L-PBF of GRCop-42. This paper will present an overview of the RAMPT project, the process development and hardware progress to date, material and hot-fire testing results, along with future developments.

Additive Manufacturing↗

DOE OSTI2021

Ion Dissociation Dynamics in an Aqueous Premelting Layer

Using molecular dynamics simulations and methods of importance sampling, we study the thermodynamics and dynamics of sodium chloride in the aqueous premelting layer formed spontaneously at the interface between ice and its vapor. We uncover a hierarchy of time scales that characterize the relaxation dynamics of this system, spanning the picoseconds of ionic motion to the tens or hundreds of nanoseconds associated with fluctuations of the liquid-crystal interface in their presence. We find that ions distort both local interfaces, incurring restoring forces that result in the ions preferentially residing in the middle of the layer. While ion pair dissociation is thermodynamically favorable, these structural and dynamic effects cause its rate to vary by over an order of magnitude through the layer, with a maximum rate significantly depressed from the corresponding bulk value. The solvation environment of ions in the premelting layer is distinct from that in a bulk liquid, being dominated by slow reorganization of water molecules and a water structure intermediate between ice and its melt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Thermal Stability of π-Conjugated n -Ethylene-Glycol-Terminated Quaterthiophene Oligomers: A Computational and Experimental Study

This work represents a joint computational and experimental study on a series of $\textit{n}$-ethylene glycol (PEO$\textit{n}$)-terminated quaterthiophene (4T) oligomers for 1 < $\textit{n}$ < 10 to elucidate their self-assembly behavior into a smectic-like lamellar phase. This study builds on an earlier study for $\textit{n}$ = 4 that showed that our model predictions were consistent with experimental data on the melting behavior and structure of the lamellar phase, with the latter consisting of crystal-like 4T domains and liquid-like PEO4 domains. The present study aims to understand how the length of the terminal PEO$\textit{n}$ chains modulates the disordering temperature of the lamellar phase and hence the relative stability of the ordered structure. In this work, a simplified bilayer model, where the 4T domains are not explicitly described, is put forward to efficiently estimate the disordering effect of the PEO domains with increasing n; this method is first validated by correctly predicting that layers of alkyl (PE)-capped 4T oligomers (for 1 < $\textit{n}$ < 10) stay ordered at room temperature. Both 4T-domain implicit and explicit model simulations reveal that the order-disorder temperature decreases with the length of the PEO capping chains, as the associated increase in conformational entropy drives a tendency toward disorder that overtakes the cohesive energy, keeping the ordered packing of the 4T domains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Lithium oxidation and electrolyte decomposition at Li-metal/liquid electrolyte interfaces

In this work, we examine the evolution of events occurring when a Li metal surface is in contact with a 2 M solution of a Li salt in a solvent or mixture of solvents, via classical molecular dynamics simulations with a reactive force field allowing bond breaking and bond forming. The main events include Li oxidation and electrolyte reduction along with expansion of the Li surface layers forming a porous phase that is the basis for the formation of the solid-electrolyte interphase (SEI) components. Nucleation of the main SEI components (LiF, Li oxides, and some organics) is characterized. The analysis clearly reveals the details of these physical–chemical events as a function of time, during 20 nanoseconds. The effects of the chemistry of the electrolyte on Li oxidation and dissolution in the liquid electrolyte, and SEI nucleation and structure are identified by testing two salts: LiPF 6 and LiCF 3 SO 3 , and various solvents including ethers and carbonates and mixtures of them. The kinetics and thermodynamics of Li 6 F, the core nuclei in the LiF crystal, are studied by analysis of the MD trajectories, and via density functional theory calculations respectively. The SEI formed in this computational experiment is the “native” film that would form upon contact of the Li foil with the liquid electrolyte. As such, this work is the first in a series of computational experiments that will help elucidate the intricate interphase layer formed during battery cycling using metal anodes.

25 ENERGY STORAGE↗

DOE OSTI2021

Stochastic parametric skeletal dosimetry model for humans: General approach and application to active marrow exposure from bone-seeking beta-particle emitters

The objective of this study is to develop a skeleton model for assessing active marrow dose from bone-seeking beta-emitting radionuclides. This article explains the modeling methodology which accounts for individual variability of the macro- and microstructure of bone tissue. Bone sites with active hematopoiesis are assessed by dividing them into small segments described by simple geometric shapes. Spongiosa, which fills the segments, is modeled as an isotropic three-dimensional grid (framework) of rod-like trabeculae that “run through” the bone marrow. Randomized multiple framework deformations are simulated by changing the positions of the grid nodes and the thickness of the rods. Model grid parameters are selected in accordance with the parameters of spongiosa microstructures taken from the published papers. Stochastic modeling of radiation transport in heterogeneous media simulating the distribution of bone tissue and marrow in each of the segments is performed by Monte Carlo methods. Model output for the human femur at different ages is provided as an example. The uncertainty of dosimetric characteristics associated with individual variability of bone structure was evaluated. An advantage of this methodology for the calculation of doses absorbed in the marrow from bone-seeking radionuclides is that it does not require additional studies of autopsy material. The biokinetic model results will be used in the future to calculate individual doses to members of a cohort exposed to 89,90 Sr from liquid radioactive waste discharged to the Techa River by the Mayak Production Association in 1949–1956. Further study of these unique cohorts provides an opportunity to gain more in-depth knowledge about the effects of chronic radiation on the hematopoietic system. In addition, the proposed model can be used to assess the doses to active marrow under any other scenarios of 90 Sr and 89 Sr intake to humans.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1985

Measurements in liquid fuel sprays

A ground test facility is being established at NASA Lewis Research Center to simulate the environmental and flight conditions needed to study adverse weather effects. One of the most important components is the water spray system which consists of many nozzles fitted on spray bars. Water is injected through air-assisted atomizers to generate uniform size drops to simulate icing in clouds. The primary objective is to provide experimental data on drop size distribution over a wide range of operating conditions. Correlation equations for mean drop size and initial injection parameters are being determined to assist in the design and modification of the Altitude Wind Tunnel. Special emphasis is being placed on the study of the aerodynamic structure of the air-assisted atomizer sprays. Detailed measurements of the variation of drop size distribution and velocity as a function of time and space are being made. Accurate initial and boundary conditions are being provided for computer model evaluation.

Chigier, N.↗

DOE OSTI2022

Transferability of data-driven, many-body models for CO 2 simulations in the vapor and liquid phases

Here, extending on the previous work by Riera et al. [J. Chem. Theory Comput. 16, 2246–2257 (2020)], we introduce a second generation family of data-driven many-body MB-nrg models for CO 2 and systematically assess how the strength and anisotropy of the CO 2 –CO 2 interactions affect the models’ ability to predict vapor, liquid, and vapor–liquid equilibrium properties. Building upon the many-body expansion formalism, we construct a series of MB-nrg models by fitting one-body and two-body reference energies calculated at the coupled cluster level of theory for large monomer and dimer training sets. Advancing from the first generation models, we employ the charge model 5 scheme to determine the atomic charges and systematically scale the two-body energies to obtain more accurate descriptions of vapor, liquid, and vapor–liquid equilibrium properties. Challenges in model construction arise due to the anisotropic nature and small magnitude of the interaction energies in CO 2 , calling for the necessity of highly accurate descriptions of the multidimensional energy landscape of liquid CO 2 . These findings emphasize the key role played by the training set quality in the development of transferable, data-driven models, which, accurately representing high-dimensional many-body effects, can enable predictive computer simulations of molecular fluids across the entire phase diagram.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗