Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Contact thermodynamics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Thermodynamic Analysis of Thermo-vacuum Clothes Drying Operation

Clothes drying accounts for a significant amount of energy consumed in residential and commercial sectors. A thermal vacuum clothes drying technology (TVCD) is proposed as an advanced clothes dryer which can significantly reduce the energy requirements by expediting the drying process. In the conventional convective clothes dryer, hot dry air is introduced into the drum which gets in direct contact to dry the clothes. This process is energy inefficient since the significant amount of heat and the water carried out with the exhaust stream are wasted. In contrast to the conventional convective drying technique, the drying mechanism of TVCD is through nucleate boiling at low temperature due to reduced vessel pressure. The process is not only efficient but also reduces the required time for drying. The objective of this paper is to develop a comprehensive thermodynamic model to predict the transient drying process of TVCD. The three-stage system model including the detail analysis of individual components is capable of simulating the variation of the water content in the textile under various operational conditions. The preliminary results show that the drying time of 3 lb textile from 70% to 2.5% in TVCD is approximately four times less than the time required in the conventional clothes dryer. Parametric studies are carried out to understand the effect of operation conditions and component geometry on the system performance. Energy consumption of the system is also analyzed.

Yang, Cheng-Min↗
DOE OSTI2026

Proximity Ferroelectricity in Compositionally Graded Structures

Proximity ferroelectricity is a novel paradigm for inducing ferroelectricity in a non-ferroelectric polar material, such as AlN or ZnO that are typically unswitchable with an external field below their dielectric breakdown field. When placed in direct contact with a thin switchable ferroelectric layer (such as Al 1-x Sc x N or Zn 1-x Mg x O), they become a practically switchable ferroelectric. Using the thermodynamic Landau-Ginzburg-Devonshire theory, in this work, we perform the finite element modeling of the polarization switching in the compositionally graded AlN-Al 1-x Sc x N, ZnO-Zn 1-x Mg x O, and MgO-Zn 1-x Mg x O structures sandwiched in both a parallel-plate capacitor geometry as well as in a sharp probe-planar electrode geometry. We reveal that the compositionally graded structure allows the simultaneous switching of spontaneous polarization in the whole system by a coercive field significantly lower than the electric breakdown field of unswitchable polar materials. The physical mechanism is the depolarization electric field determined by the gradient of chemical composition “x”. The field lowers the steepness of the switching barrier in the otherwise unswitchable parts of the compositionally graded AlN-Al 1-x Sc x N and ZnO-Zn 1-x Mg x O structures. In the MgO-like regions of the compositionally graded MgO-Zn 1-x Mg x O structure, a shallow double-well free energy potential emerges. Proximity ferroelectric switching of the compositionally graded structures placed in the probe-electrode geometry occurs due to nanodomain formation under the tip. We predict that a gradient of chemical composition “x” significantly lowers effective coercive fields of the compositionally graded AlN-Al 1-x Sc x N and ZnO-Zn 1-x Mg x O structures compared to the coercive fields of the corresponding multilayers with a uniform chemical composition in each layer. A tip-induced switching further lowers the coercive field, enabling control of ferroelectric domains in otherwise unswitchable compositionally graded structures, which can provide nanoscale domain control for memory, actuation, sensing, and optical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Li + /H + Exchange in Solid-State Oxide Li-Ion Conductors

Understanding the moisture stability of oxide Li-ion conductors is important for their practical applications in solid-state batteries. Unlike sulfide or halide conductors, oxide conductors generally better resist degradation when in contact with water but can still undergo topotactic Li + /H + exchange (LHX). Here, we combine density functional theory (DFT) calculations with a machine-learning interatomic potential model to investigate the thermodynamic driving force of the LHX reaction for two representative oxide Li-ion conductor families: garnets and NASICONs. Li-stuffed garnets exhibit a strong driving force for proton exchange due to their high Li chemical potential. In contrast, NASICONs demonstrate a higher resistance against proton exchange due to the lower Li chemical potential and the lower O–H bond covalency for polyanion-bonded oxygens. Our findings reveal a critical trade-off: Li stuffing enhances conductivity but increases moisture susceptibility. This study underscores the importance of designing Li-ion conductors that possess both high conductivity and high stability in practical environments.

Li, Zhuohan [Lawrence Berkeley National Laboratory↗
NASA NTRS2002

Bubble Detachment in Variable Gravity Under the Influence of Electric Fields

The objective of the research is to investigate the behavior of individual air bubbles injected through an orifice into an electrically insulating liquid under the influence of a static electric field. Situations were considered with both uniform and nonuniform electric fields. Bubble formation and detachment were visualized in terrestrial gravity as well as for several levels of reduced gravity (lunar, martian and microgravity) using a high-speed video camera. Bubble volume, dimensions and contact angles at detachment were measured. In addition to the experimental studies, a simple model, predicting bubble characteristics at detachment in an initially uniform electric field was developed. The model, based on thermodynamic considerations, accounts for the level of gravity as well as the magnitude of the uniform electric field. The results of the study indicate that the level of gravity and the electric field magnitude significantly affect bubble behavior as well as shape, volume and dimensions.

Herman, Cila↗
DOE OSTI2022

The “good,” the “bad,” and the “hidden” in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl 2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca 2+ –Cl - interactions were primarily “hidden” in the Ca-NDIS signal due to overlapping with Ca 2+ –O w and Ca 2+ –H w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models “on average,” the question “which model is the best” has not been answered but replaced by the question “which model is better for a given research.” An overall “good” model can still be inappropriate or, in some instances, “bad” and, unfortunately, produce erroneous results. Additionally, the accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

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

Increasing Ionic Conductivity of Poly(ethylene oxide) by Reaction with Metallic Li

Poly(ethylene oxide) (PEO) was the first lithium-ion conducting polymer developed 50 years ago and is still the most popular electrolyte matrix for solid-state lithium metal batteries. While many studies focus on increasing PEO ionic conductivity through doping with Li salts, little work has addressed using PEO and Li directly to generate Li + -conducting species in situ. Reaction between PEO and Li leads to ionic conductivity largely from Li + , in contrast to the case of added salts where the anion contribution dominates. Herein, electrochemical impedance spectroscopy shows the ionic conductivity of PEO thin films increases up to three orders of magnitude (from 10 –7 to 10 –4 S cm –1 ) when contacted with Li at elevated temperature. This is due to the reduction of ether bonds, which produces lithium alkoxides that are responsible for Li + transport. Density functional theory analysis confirms this mechanism as thermodynamically favorable. X-ray photoelectron spectroscopy also shows the presence of organolithium species and Li 2 O, which are responsible for propagating reactions with PEO and forming an electronically insulating layer at the PEO–Li interface that halts further reaction, respectively. The underlying mechanisms of Li–polymer electrolyte reactions is clarified and new pathways for in situ Li + doping of polymer electrolytes is presented.

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

Comparative Energetics of Various Membrane Distillation Configurations and Guidelines for Design and Operation

This paper presents a comparative performance study of single-stage desalination processes with major configurations of membrane distillation (MD) modules. MD modules covered in this study are (a) direct contact MD (DCMD), (b) vacuum MD (VMD), (c) sweeping gas MD (SGMD), and (d) air gap MD (AGMD). MD-based desalination processes are simulated with rigorous theoretical MD models supported by molecular thermodynamic property models for the accurate calculation of performance metrics. The performance metrics considered in MD systems are permeate flux and energy efficiency, i.e., gained output ratio (GOR). A general criterion is established to determine the critical length of these four MDs (at fixed width) for the feasible operation of desalination in a wide range of feed salinities. The length of DCMD and VMD is restricted by the feed salinity and permeate flux, respectively, while relatively large AGMD and SGMD are allowed. The sensitivity of GOR flux with respect to permeate conditions is investigated for different MD configurations. AGMD outperforms other configurations in terms of energy efficiency, while VMD reveals the highest permeate production. With larger MD modules, utilization of thermal energy supplied by the hot feed for evaporation is in the order of VMD > AGMD > SGMD > DCMD. Simulation results highlight that energy efficiency of the overall desalination process relies on the efficient recovery of spent for evaporation, suggesting potential improvement in energy efficiency for VMD-based desalination.

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

Magnetohydrodynamic models of coronal transients in the meridional plane. II - Simulation of the coronal transient of 1973 August 21

A two-dimensional planar MHD model of solar atmospheric transient events is used to simulate the flare-associated events of August 21, 1973. This event, observed in H-alpha, He II 304 A, soft X-ray, and coronal white light, provided sufficient information (especially in the latter diagnostic) for an assessment to be made of the model's ability to simulate major features of an actual solar event. It was found that a thermodynamic input pulse based on data provided by the NASA Marshall Space Flight Center-Aerospace Corporation X-ray telescope (S-056) on Skylab was sufficient to produce the global geometry, shock and contact-surface velocities, excess mass contours, and energy budget which were, for the most part, observed by the High Altitude Observatory white-light coronagraph (S-052) on Skylab in the form of a forerunner and coronal transient.

Dryer, M.↗
DOE OSTI2025

Colloidal Behavior of Plutonium Oxide in Concentrated Electrolyte Solutions

The Hanford Site in Washington State manages legacy high-level radioactive waste streams that display major chemistry and engineering challenges, including the high salt levels and pH values that correspond to conditions under which many classical concepts describing chemical reactivity cannot be applied. One particular challenge that needs to be tackled at Hanford is that Pu concentrations, [Pu], in the soluble phases of the tank wastes are higher than expected based on the solubility of crystalline PuO2, which is widely accepted to be caused by the formation of PuO2 colloid, consisting of nano- to submicron-sized particles (PuO2 NPs). Fundamental research underpinning the behavior of PuO2 NPs under conditions not only relevant to the Hanford tank waste but at high ionic strength in general is needed to reliably predict the chemical reactivity of PuO2 NPs and develop engineering solutions to safely and efficiently process high-level radioactive waste into forms suitable for long-term storage. In this work, we study the behavior of PuO2 NPs (particle size ~100 nm) under high ionic strength conditions by reacting it with highly concentrated (up to 5 M) salt solutions. We explore different electrolyte compositions to elucidate the impact of different anions (NO3-, Cl-, ClO4-, SO42-, C2O42-, CO32-) on the stability of PuO2 NP in the acidic and alkaline pH regime. PuO2 NP aggregation and precipitation as function electrolyte concentration is tracked by a combination of liquid scintillation counting, dynamic light scattering for determination of particle size distributions, and zeta potentials as a proxy for particle charge. At acidic pH, electrolytes containing non-coordinating anions, such as NaNO3, NaCl, and NaClO4 mostly stabilize PuO2 NPs over a large electrolyte concentration range, showing only subtle differences in their reactivity. Other electrolyte anions show a more pronounced effect on the PuO2 NP stability: SO42-, binds directly to the particles’ surface, reverses the particle charge, and precipitates the PuO2 NPs efficiently even at intermediate sulfate concentrations (>0.1 M). In contrast, C2O42- is found to lead to high [Pu] in solution, in the milli-molar range, even at mildly acidic pH (~4). Thermodynamic modeling of the dissolved Pu concentrations using PHREEQC is unable to predict the observed [Pu] in the acidic pH regime, supporting the influence of colloids in maintaining elevated [Pu]. It is noteworthy that the current thermodynamic databases do not include constants for colloidal Pu phases and cannot accurately predict many of the high ionic strength solutions relevant to this work. The mechanisms and models responsible for these observations will need further investigation in the future. At high pH values (~12), PuO2 NPs exhibits classical sol-gel chemistry, meaning that upon destabilization of the colloidal sol, for example by addition of concentrated NaOH, highly porous and viscid PuO2 coagulates are formed that consist of a three-dimensional network likely held together by physical interactions. The PuO2 NP coagulate shows no significant reversibility of the aggregation when contacted with concentrated brines; however, PuO2 NPs can be efficiently resuspended in solution by addition of diluted electrolytes, alkaline solutions containing high amounts of carbonate, or simple addition of water. Especially carbonate is shown to stabilize PuO2 NPs in solution at high pH, characterized by stable colloidal suspensions that are resistant against sedimentation during centrifugation. Thermodynamic modeling of the carbonate system was able to predict an increasing dissolved Pu concentration with increasing carbonate concentration. However, the model was profoundly sensitive to the fixed redox potential and does not include any thermodynamic constants for colloidal Pu species.

Neumann, Julia↗
DOE OSTI2022

Ion Transport in Batteries with Polymer Electrolytes

We discuss polymer electrolytes for use in rechargeable lithium batteries. Polymer electrolytes have the potential to enable batteries with lithium metal anodes. These batteries have significantly higher theoretical energy densities than current lithium-ion batteries. We consider binary mixtures of polymers and salts. We also cover more complex systems such as polymer electrolytes swollen with a solvent (gel polymer electrolytes) and microphase separated polymer electrolytes. By covalently attaching the anions to the chains in a polymer solid, one obtains a single-ion conductor. We mainly focus on experiments wherein the polymer electrolyte is placed between two lithium metal electrodes. These experiments enable the determination of three transport parameters, ionic conductivity, salt diffusion coefficient, and transference number, and the thermodynamic factor. The properties of dry polymer electrolytes are contrasted with those of gel polymer electrolytes. The gel systems exhibit higher conductivity while the dry systems exhibit superior mechanical properties. We discuss interfacial impedance when lithium metal is contacted with polymer electrolytes and the importance of coulombic efficiency.

Yu, Xiaopeng↗
DOE OSTI2023

High-fidelity imaging of a band insulator in a three-dimensional optical lattice clock

Here, we report on the observation of a high-density, band insulating state in a three-dimensional optical lattice clock. Filled with a nuclear-spin-polarized degenerate Fermi gas of 87 Sr, the three-dimensional (3D) lattice has one atom per site in the ground motional state, thus guarding against frequency shifts due to contact interactions. At this high density where the average distance between atoms is comparable to the probe wavelength, conventional imaging techniques at saturation intensity suffer from large systematic errors. To spatially probe frequency shifts in the clock and measure thermodynamic properties of this system, accurate imaging techniques at high optical depths are required. Using a combination of highly saturated fluorescence and absorption imaging, we confirm the density distribution in our 3D optical lattice in agreement with a single spin band insulating state. Combining our clock platform with this high filling fraction opens the door to studying new classes of long-lived, many-body states arising from dipolar interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Nickel hydroxide–nickel carbonate competitive growth on carbonate surfaces

The thermodynamic and kinetic factors controlling the competitive heterogeneous nucleation and growth of ubiquitous metal carbonate and hydroxide phases are poorly understood. In this work, calcite (CaCO 3 ) and magnesite (MgCO 3 ) powders were reacted with NiCl 2 (0–600 μM) for 7 days at 22 °C and 5 °C. The reacted powders were analyzed with X-ray photoelectron spectroscopy (XPS), scanning electron microscopy, and energy-dispersive X-ray spectroscopy to characterize the Ni surface precipitates formed. Evidence from these techniques pointed to the formation of mixed Ni carbonate-Ni hydroxide amorphous surface precipitates. On calcite, XPS detected primarily Ni(OH) 2 despite the initial solutions being more supersaturated with respect to gaspéite (NiCO 3 ) than to theophrastite (Ni(OH) 2 ) by a factor of 17–18. In contrast, NiCO 3 was the dominant component detected by XPS on magnesite in the same conditions. Decreasing the temperature had the effect of increasing the proportion of NiCO 3 to the detriment of Ni(OH) 2 . The experimental observations were consistent with low lattice/cation size mismatch favoring NiCO 3 nucleation and temperature most influencing Ni(OH) 2 nucleation. Comparison to previous work on Co-reacted powders indicated the differences in lattice/cation size mismatch and/or water exchange rate impacted the composition of the surface precipitates more than the relative thermodynamic stabilities of the competing minerals. This work explored the heterogeneous growth regime of Ni carbonate and hydroxide phases on carbonate surfaces and shed light on the factors that control the competition between surface precipitates when mineral surfaces are in contact with aqueous solutions supersaturated with respect to multiple mineral phases. In conclusion, these results contribute to geochemists’ efforts toward interpreting data from geochemical systems with elevated Ni concentrations, improving Ni environmental remediation and recovery strategies, and predicting the fate and transport of Ni in geochemical systems.

Heterogeneous nucleation↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
DOE OSTI2022

Hole Trapping in Halide Perovskites Induces Phase Segregation

Metal halide perovskites have garnered a great deal of attention for their applications in photovoltaics, LEDs, and radiation detection. The ease of solution processing high-quality perovskite semiconductors with large absorption coefficients and tolerance to native defects is decidedly attractive. Additionally, the ability to precisely tune the band gap of halide perovskites through compositional alloying of the halide ion is of particular interest for a range of applications, especially for tandem solar cells. However, under steady state light irradiation, an initially homogeneous mixed halide perovskite (MHP) will form local domains that are rich in one halide ion (e.g., Br or I). This light-induced phase segregation in MHPs forms iodide-rich domains that act as charge carrier traps and lowers the efficiency of perovskite-based devices. Thus, phase segregation poses a serious challenge to the implementation of MHPs in real-world device settings. Interestingly, when a phase segregated MHP film is placed in the dark, entropic driving forces become dominant and the segregated perovskite remixes and returns to its initially homogeneous state. Several key mechanistic details of phase segregation have been elucidated over the years. However, there are still aspects of halide segregation that are not clear, and there is ongoing debate in the literature as to what are the key factors that contribute to the mechanism. This Account discusses recent results that point to the specific role of hole trapping in phase segregation. Interestingly, generation of holes through above-band-gap excitation or through electrochemical injection increases ion migration and leads to phase segregation. The thermodynamic and redox properties of halide perovskites provide a strong driving force for hole trapping and oxidation of iodide species in MHPs. However, mobile halide species within the perovskite lattice take time to migrate and generate halide-rich domains. When in contact with a nonpolar solvent, the migration of iodine species is further extended to expulsion of iodine from the perovskite film. Thus, the mobility of halides and their susceptibility to hole-induced oxidation play a crucial role in determining the long-term stability of metal halide perovskites. Strategies to gain kinetic control over ion migration to slow phase segregation are needed to overcome these hurdles and achieve stable mixed halide perovskites. Modification of the perovskite composition through introduction of different cations or halide ions, or introduction of low-dimensional perovskite phases may suppress phase segregation. Furthermore, in achieving stability and improving the efficiency of perovskite solar cells and light emitting devices with minimal impacts, suppression of segregation remains the key factor.

36 MATERIALS SCIENCE↗
NASA NTRS1974

A statistical-chemical and thermodynamic approach to the study of lunar mineralogy

Principal components analysis is used to study the chemical compositions of pyroxenes of five Apollo 12 specimens. Important correlations are recognized in the variation of oxide weight per cent. These correlations indicating substitutional relationships can be interpreted as representative of stable and metastable trends of crystallization by using crystal-chemical and thermodynamic information. The per cent variance of pyroxene groups with characteristic trends in each specimen can be evaluated and interpreted in terms of history of crystallization. Distribution of Fe and Mg in certain pairs of olivine and pyroxene, which are found in contact in the rock and which may have crystallized simultaneously, is useful in recognizing the tendency towards chemical equilibrium in Fe-Mg distribution during a limited interval in the liquidus or subsolidus stages.

Saxena, S. K.↗
NASA NTRSDate not supplied

Development and Validation of a Slosh-Based Ullage Collapse Model

Ullage collapse modeling is an essential tool for the success of spacecraft powered by cryogenic propellants. The degree of ullage collapse can have a major impact on the design of cryogenic propellant tanks as well as the mission’s concept of operations. Additional commodities (e.g., helium gas for pressure control) and their quantities are dependent on the amount of heat transfer and phase change of propellants expected during flight. Currently, high-fidelity computational fluid dynamics (CFD) with heat transfer and phase change is the most accurate method for predicting ullage collapse. CFD simulations of this kind, however, can become very costly for problems requiring large domains, long simulation times, or significant liquid motion. Consequently, there is a need for a reduced order model to provide fast responses which can be used in engineering design decisions. The proposed numerical model is a nodal approach, augmented by outputs from CFD analyses that are generated rapidly and reliably. These outputs include transient values for liquid-vapor interface area, interface velocity, and liquid position using a volume of fluid (VoF) method. There are several benefits of combining a nodal approach with slosh predictions from CFD simulations: leveraging of the Loci/STREAM-VoF solver that has been validated in high-g and low-g environments, rapid modeling of VoF with isothermal propellant simplifications and relatively coarse meshes, and improved model predictions from the nodal tool. The nodal approach, informed by the slosh predictions from CFD, assumes uniform conditions in the liquid and gas domains. Modifications to account for temperature stratification in the liquid and gas may be implemented as a future improvement. Thermodynamic changes within the tank are solved at user-defined timesteps. First, impacts due to mass entering or leaving the control volume are captured. Second, heat transfer impacts are accounted for by tracking the volume of the gas which has been contacted by liquid in that timestep and rebalancing the internal energy. Finally, phase change is approximated using the area and velocity of the interface. The tank conditions are re-computed by conserving mass and energy. Calibration of the model to several test flights and experiments are required to baseline the predictive power and uncertainty. The final product will enable designers and engineers to create quick-turnaround predictions for ullage collapse and inform design feasibility.

Computational Fluid Dynamics↗