Search NASA⌕ Search

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.

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

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Bond-orientational order in liquid Si

Bond-orientational order in liquid Si via Monte Carlo simulation in conjuncation with empirical two- and three-body potentials of the form proposed by Stillinger and Weber are studied. Bond-orientational order (BOO) is described in terms of combinations of spherical harmonic functions. Liquid Si is found to have pronounced short-range BOO corresponding to l = 3, as expected for a structure with local tetrahedral order. No long-range BOO is found either in the equilibrium or the supercooled liquid. When the three-body potential is artificially removed, the tetrahedral bond-orientation order disappears and the liquid assumes a close-packed structure.

Wang, Z. Q.↗
NASA NTRS2002

Incompressible Navier-Stokes Solvers in Primative Variables and their Applications to Steady and Unsteady Flow Simulations

This paper reviews recent progress made in incompressible Navier-Stokes simulation procedures and their application to problems of engineering interest. Discussions are focused on the methods designed for complex geometry applications in three dimensions, and thus are limited to primitive variable formulation. A summary of efforts in flow solver development is given followed by numerical studies of a few example problems of current interest. Both steady and unsteady solution algorithms and their salient features are discussed. Solvers discussed here are based on a structured-grid approach using either a finite -difference or a finite-volume frame work. As a grand-challenge application of these solvers, an unsteady turbopump flow simulation procedure has been developed which utilizes high performance computing platforms. In the paper, the progress toward the complete simulation capability of the turbo-pump for a liquid rocket engine is reported. The Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for evaluation of two parallel computing algorithms that have been implemented in the INS3D code. The relative motion of the grid systems for the rotorstator interaction was obtained using overact grid techniques. Unsteady computations for the SSME turbo-pump, which contains 114 zones with 34.5 million grid points, are carried out on SCSI Origin 3000 systems at NASA Ames Research Center. The same procedure has been extended to the development of NASA-DeBakey Ventricular Assist Device (VAD) that is based on an axial blood pump. Computational, and clinical analysis of this device are presented.

Kiris, Cetin C.↗
NASA NTRS2004

A Study of the Response of Deep Tropical Clouds to Mesoscale Processes: Sensitivities to Microphysics, Radiation, and Surface Fluxes - Part 2

The Goddard Cumulus Ensemble (GCE) model is used to examine the sensitivities of surface fluxes, explicit radiation, and ice microphysical processes on multi-day simulations of deep tropical convection over the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE). The simulations incorporate large-scale advective temperature and moisture forcing, as well as large-scale momentum, that are updated every time step on a periodic lateral boundary grid. This study shows that when surface fluxes are eliminated, the mean atmosphere is much cooler and drier, convection and CAPE are much weaker, precipitation is less, and cloud coverage in stratiform regions much greater. Surface fluxes using the TOGA COARE flux algorithm are weaker than with the aerodynamic formulation, but closer to the observed fluxes. In addition, similar trends noted above for the case without surface fluxes are produced for the TOGA flux case, albeit to a much lesser extent. The elimination of explicit shortwave and longwave radiation is found to have only minimal effects on the mean thermodynamics, convection, and precipitation. However explicit radiation does have a significant impact on cloud temperatures and structure above 200 mb and on the overall mean vertical circulation. The removal of ice processes produces major changes in the structure of the cloud. Much of the liquid water is transported aloft and into anvils above the melting layer (600 mb), leaving narrow, but intense bands of rainfall in convective regions. The elimination of melting processes leads to greater hydrometeor mass below the melting layer, and produces a much warmer and moister boundary layer, leading to a greater mean CAPE. Finally, the elimination of the graupel species has only a small impact on mean total precipitation, thermodynamics, and dynamics of the simulation, but does produce much greater snow mass just above the melting layer. Some of these results differ from previous CRM studies of tropical systems, which is likely due to the type of simulated system, total time integration, and model setup.

Johnson, Daniel↗
DOE OSTI2018

Understanding Electric Double-Layer Gating Based on Ionic Liquids: from Nanoscale to Macroscale

In electric double-layer transistors (EDLTs), it is well known that the EDL formed by ionic liquids (ILs) can induce an ultrahigh carrier density at the semiconductor surface, compared to solid dielectric. However, the mechanism of device performance is still not fully understood, especially at a molecular level. In this study, we evaluate the gating performance of amorphous indium gallium zinc oxide (a-IGZO) transistor coupled with a series of imidazolium-based ILs, using an approach combining of molecular dynamics simulation and finite element modeling. Results reveal that the EDL with different ion structures could produce inhomogeneous electric fields at the solid–electrolyte interface, and the heterogeneity of electric field-induced charge distributions at semiconductor surface could reduce the electrical conductance of a-IGZO during gating process. Meanwhile, a resistance network analysis was adopted to bridge the nanoscopic data with the macroscopic transfer characteristics of IL-gated transistor, and showed that our theoretical results could well estimate the gating performance of practical devices. Finally, our findings could provide both new concepts and modeling techniques for IL-gated transistors.

36 MATERIALS SCIENCE↗
NASA NTRS2021

Application of Radon Transform to Multi-Angle Measurements Made by the Research Scanning Polarimeter: A New Approach to Cloud Tomography. Part I: Theory and Tests on Simulated Data.

The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP was a prototype for the Aerosol Polarimetry Sensor (APS) launched on-board the NASA Glory satellite. Currently the retrieval algorithms developed for the RSP are being adopted for the measurements of the space-borne polarimeters on the upcoming NASA’s Plankton, Aerosol, Cloud Ocean Ecosystem (PACE)satellite mission. The RSP’s uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid water cloud droplet sizes using the polarized rainbow structure. It also provides geometric constraints on the cumulus cloud’s 2D cross section yielding the cloud’s geometric shape estimates. In this study we further build on the latter technique to develop a new tomographic approach to retrieval of cloud internal structure from remote sensing measurements. While tomography in the strict definition is a technique based on active measurements yielding a tomogram (directional optical thickness as a function of angle and offset of the view ray), we developed a “semi-tomographic” approach in which tomogram of the cloud is estimated from passive observations instead of being measured directly. This tomogram is then converted into 2D spatial distribution of the extinction coefficient using inverse Radon transform (filtered back projection) which is the standard tomographic procedure used e.g., in medical CT scans. This algorithm is computationally inexpensive compared to techniques relying on highly-multi-dimensional least-square fitting; it does not require iterative 3D RT simulations. The resulting extinction distribution is defined up to an unknown constant factor, so we discuss the ways to calibrate it using additional independent measurements. In the next step we use the profile of the droplet size distribution parameters from the cloud’s side (derived by fitting the polarized rainbows) to convert the 2D extinction distribution into that of the droplet number concentration. We illustrate and validate the proposed technique using 3D-RT-simulatedRSP observations of a LES-generated Cu cloud. Quantitative comparisons between the retrieved and the original optical and microphysical parameters are presented.

clouds↗
DOE OSTI2023

Multiscale molecular simulations for the solvation of lignin in ionic liquids

Lignin, the second most abundant biopolymer found in nature, has emerged as a potential source of sustainable fuels, chemicals, and materials. Finding suitable solvents, as well as technologies for efficient and affordable lignin dissolution and depolymerization, are major obstacles in the conversion of lignin to value-added products. Certain ionic liquids (ILs) are capable of dissolving and depolymerizing lignin but designing and developing an effective IL for lignin dissolution remains quite challenging. To address this issue, the COnductor-like Screening MOdel for Real Solvents (COSMO-RS) model was used to screen 5670 ILs by computing logarithmic activity coefficients (ln(γ)) and excess enthalpies (H E ) of lignin, respectively. Based on the COSMO-RS computed thermodynamic properties (ln(γ) and H E ) of lignin, anions such as acetate, methyl carbonate, octanoate, glycinate, alaninate, and lysinate in combination with cations like tetraalkylammonium, tetraalkylphosphonium, and pyridinium are predicted to be suitable solvents for lignin dissolution. The dissolution properties such as interaction energy between anion and cation, viscosity, Hansen solubility parameters, dissociation constants, and Kamlet–Taft parameters of selected ILs were evaluated to assess their propensity for lignin dissolution. Furthermore, molecular dynamics (MD) simulations were performed to understand the structural and dynamic properties of tetrabutylammonium [TBA] + -based ILs and lignin mixtures and to shed light on the mechanisms involved in lignin dissolution. MD simulation results suggested [TBA] + -based ILs have the potential to dissolve lignin because of their higher contact probability and interaction energies with lignin when compared to cholinium lysinate.

09 BIOMASS FUELS↗
DOE OSTI2021

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

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

Characterization of Whole Biomasses in Pyridine Based Ionic Liquid at Low Temperature by 31 P NMR: An Approach to Quantitatively Measure Hydroxyl Groups in Biomass As Their Original Structures

In this paper, the dissolution of biomass components—cellulose, hemicellulose, and lignin, and two whole biomasses—switchgrass and poplar in a pyridine based ionic liquid at a low temperature—50°C has been examined, which will provide an opportunity to explore the original structures of biomass components. The following phosphitylation, and 31 P NMR measurement could provide quantitative results for various hydroxyl groups, including aliphatic, condensed phenolic, guaiacyl phenolic, p -hydroxyl phenyl and carboxylic hydroxyl groups in the biomass components, and whole biomass. By employing various biomass model compounds (glucose, cellotriose, and cellohexose), artificial mixtures of biomass components (cellulose, hemicellulose, and lignin), and computational simulation for the assignments by using density functional theory calculation in Gaussian, reliability and accuracy of this method have been examined as well, which indicated that this method is a reliable and accurate way to quantitatively characterize five different types of hydroxyl groups in biomass and its components.

09 BIOMASS FUELS↗
DOE OSTI2018

Unraveling Kinetically-Driven Mechanisms of Gold Nanocrystal Shape Transformations Using Graphene Liquid Cell Electron Microscopy

Mechanisms of kinetically driven nanocrystal shape transformations were elucidated by monitoring single particle etching of gold nanocrystals using in situ graphene liquid cell transmission electron microscopy (TEM). By systematically changing the chemical potential of the oxidative etching and then quantifying the facets of the nanocrystals, nonequilibrium processes of atom removal could be deduced. Etching at sufficiently high oxidation potentials, both cube and rhombic dodecahedra (RDD)-shaped gold nanocrystals transform into kinetically stable tetrahexahedra (THH)-shaped particles. Whereas {100}-faceted cubes adopt an { hko }-faceted THH intermediate where h/k depends on chemical potential, {110}-faceted RDD adopt a {210}-faceted THH intermediate regardless of driving force. For cube reactions, Monte Carlo simulations show that removing 6-coordinate edge atoms immediately reveals 7-coordinate interior atoms. The rate at which these 6- and 7-coordinate atoms are etched is sensitive to the chemical potential, resulting in different THH facet structures with varying driving force. Conversely, when RDD are etched to THH, removal of 6-coordinate edge atoms reveals 6-coordinate interior atoms. Thus, changing the driving force for oxidation does not change the probability of edge atom versus interior atom removal, leading to a negligible effect on the kinetically stabilized intermediate shape. These fundamental insights, facilitated by single-particle liquid-phase TEM imaging, provide important atomic-scale mechanistic details regarding the role of kinetics and chemical driving force in dictating shape transformations at the nanometer length scale.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2018

Energy Storage in Metallic Glasses via Flash Annealing

Abstract When the temperature rises in a kinetically frozen and out‐of‐equilibrium system, such as a structural glass, thermal activation usually drives a relaxation to a lower enthalpy state. This is the essence of the ubiquitous thermal ageing of glassy materials below the glass transition temperature. Here, it is shown that an ultrafast temperature rise without subsequent quenching can elevate the system to a higher enthalpy state, thereby erasing the relaxation history. The key of this temperature‐driven rejuvenation is a fast‐heating rate (flash annealing) that brings the glass into a temperature regime above the glass transition within less than a millisecond. Here, the structure is entropically driven towards the liquid state, lowering its free energy, but raising its enthalpy. Subsequent cooling freezes the increased enthalpy before the liquid state is reached, thereby leading to a global energy storage. This ultrafast rejuvenation process is demonstrated for a Zr‐based metallic glass, combining calorimetry, fast in‐situ X‐ray diffraction, and atomistic simulations, showing a tenfold increase in stored enthalpy via flash annealing.

Küchemann, Stefan↗
DOE OSTI2022

Miscibility of rock and ice in the interiors of water worlds

Abstract Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO 3 ), a major end-member silicate phase, and water (H 2 O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO 3 . Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

MHD Analysis of Dual-Coolant Lead-Lithium Blanket for Spherical Tokamak Advanced Reactor

The tritium breeding blanket is vital for future fusion power plants, with the Spherical Tokamak Advanced Reactor (STAR) project highlighting the dual-coolant lead-lithium (DCLL) design. The DCLL blanket performs shielding, energy exhaust, and tritium breeding using a lead-lithium alloy , with lithium as the breeder and lead as the neutron multiplier. It also serves as the primary coolant, with helium providing supplemental cooling. Reduced-activation ferritic/martensitic steel is used for the blanket structure. Magnetohydrodynamic (MHD) phenomena influence the liquid metal flow in a magnetic field, affecting heat transfer in the breeder affected by energetic neutrons. Understanding key flow parameters in such conditions is critical for efficient DCLL design. This study uses three-dimensional thermofluid MHD analysis with ANSYS CFX software, modified at Princeton Plasma Physics Laboratory, to simulate high Hartmann flows. The neutronics code MCNP, coupled with plasma equilibrium, provides heat source distribution. In conclusion, we examine electromagnetic interactions in adjacent fluid domains and analyze the magnetic field’s impact on flow distribution in the inboard and outboard blanket layout, using detailed mesh generation for accurate results.

DCLL↗
DOE OSTI2024

Validation Framework for Post-Combustion Carbon Capture CFD Simulations

First-principles based computational fluid dynamics (CFD) simulations are proposed as a fundamental tool for investigating solvent-based CO2 absorption in packed columns, due to the ability to accurately represent the underlying non-linear multiscale dynamics. In this work, we employ such models to investigate hydrodynamics of columns with structured by assessing the key hydrodynamic metrics, such as pressure drop and liquid holdup. Our models are validated with experimental data from a specifically designed column for this line of work. The test cases of gas and liquid flowrates and operating conditions were selected through a comprehensive sequential design of experiments approach offered by the CCSI2 toolset.

Panagakos, Grigorios↗
NASA NTRS1972

The design and operation of a planet emitted and albedo radiation simulator

A simulator is described, which provides a dynamic simulation of planetary albedo and planet-emitted thermal radiation. The simulator was designed for model tests of advanced thermal control concepts, the development of space hardware, and flight qualification of small payloads. The simulator was designed to be used in a vacuum chamber equipped with liquid nitrogen cold walls and a solar simulator. The simulator consists of two major components: a gimbaled model support system and an array of quartz infrared lamps. Tests indicate that the apparatus will satisfactorily simulate the secondary radiation on a 0.92-meter spherical body in near-earth orbits. Calculations indicate that other shapes can be accommodated without structural modifications.

Sweet, G. E.↗
NASA NTRS1978

Durability of foam insulation for LH2 fuel tanks of future subsonic transports

In connection with the potential short-supply of petroleum based fuels, NASA has initiated investigations concerning the feasibility of aircraft using as fuel hydrogen which is to be stored in liquid form. One of the problems to be solved for an operation of such aircraft is related to the possibility of a suitable storage of the liquid hydrogen. A description is presented of an experimental study regarding the suitability of commercially available organic foams as cryogenic insulation for liquid hydrogen tanks under extensive thermal cycling typical of subsonic airline type operation. Fourteen commercially available organic foam insulations were tested. The thermal performance of all insulations was found to deteriorate with increased simulated flight cycles. Two unreinforced polyurethane foams survived over 4200 thermal cycles (representative of approximately 15 years of airline service) without evidence of structural deterioration. The polyurethane foam insulations also exhibited excellent thermal performance.

Sharpe, E. L.↗
NASA NTRS2015

Testing and Analytical Modeling for Purging Process of a Cryogenic Line

The purging operations for cryogenic main propulsion systems of upper stage are usually carried out for the following cases: 1) Purging of the Fill/Drain line after completion of propellant loading. This operation allows the removal of residual propellant mass; and 2) Purging of the Feed/Drain line if the mission is scrubbed. The lines would be purged by connections to a ground high-pressure gas storage source. The flow-rate of purge gas should be regulated such that the pressure in the line will not exceed the required maximum allowable value. Exceeding the maximum allowable pressure may lead to structural damage in the line. To gain confidence in analytical models of the purge process, a test series was conducted. The test article, a 20-cm incline line, was filled with liquid hydrogen and then purged with gaseous helium (GHe). The influences of GHe flow-rates and initial temperatures were evaluated. The Generalized Fluid System Simulation Program, an in-house general-purpose computer program for flow network analysis, was utilized to model and simulate the testing. The test procedures, modeling descriptions, and the results will be presented in the final paper.

Hedayat, A.↗
DOE OSTI2020

Tensor network investigation of the double layer Kagome compound Ca 10 Cr 7 O 28

Quantum spin liquids are exotic quantum phases of matter that do not order even at zero temperature. While there are several toy models and simple Hamiltonians that could host a quantum spin liquid as their ground state, it is very rare to find actual, realistic materials that exhibit their properties. At the same time, the classical simulation of such instances of strongly correlated systems is intricate and reliable methods are scarce. Here, in this work, we investigate the quantum magnet Ca 10 Cr 7 O 28 that has recently been discovered to exhibit properties of a quantum spin liquid in inelastic neutron scattering experiments. This compound has a distorted bilayer Kagome lattice crystal structure consisting of Cr 5+ ions with spin- moments. Coincidentally, the lattice structure renders a tensor network algorithm in 2D applicable that can be seen as a new variant of a projected entangled simplex state algorithm in the thermodynamic limit. In this first numerical investigation of this material that takes into account genuine quantum correlations, good agreement with the experimental findings is found. Our study contributes to uplifting tensor networks from conceptual tools to methods to describe real two-dimensional quantum materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Experimentally validated high-fidelity simulations of a liquid jet in supersonic crossflow

Here, utilizing recent advancements in computational schemes for compressible, multiphase flows, this work features a parametric study of a pure liquid jet in supersonic crossflow that involves simulating the atomization process for four values of momentum-flux ratio. These simulations are validated against experimental results measured with high-speed X-ray imaging, which confirm the accuracy of the numerical approach. Also, the effect of numerical resolution on some flow behavior is investigated, revealing convergence of the jet shape and surface instability wavelength. Analysis of the resulting sprays includes statistical descriptions of the liquid distribution, liquid structures created through breakup, interfacial instabilities, and dominant flow features. As the flowrate increases, the spray penetrates further, it becomes more disperse, and less liquid impacts the wall, but the droplet size distribution changes little. The wavelength of instabilities on the windward side of the jet diverges from measured trends in subsonic crossflows. In a visualization of the time-averaged flow, counter-rotating vortices are observed along the jet core and in the wake, affecting the process of primary atomization and early droplet trajectories.

42 ENGINEERING↗