Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffuse interface methods”

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 OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1980

Advanced methods for preparation and characterization of infrared detector materials

Mercury cadmium telluride crystals were prepared by the Bridgman method with a wide range of crystal growth rates and temperature gradients adequate to prevent constitutional supercooling under diffusion-limited, steady state, growth conditions. The longitudinal compositional gradients for different growth conditions and alloy compositions were calculated and compared with experimental data to develop a quantitative model of the crystal growth kinetics for the Hg(i-x)CdxTe alloys, and measurements were performed to ascertain the effect of growth conditions on radial compositional gradients. The pseudobinary HgTe-CdTe constitutional phase diagram was determined by precision differential thermal analysis measurements and used to calculate the segregation coefficient of Cd as a function of x and interface temperature. Computer algorithms specific to Hg(1-x)CdxTe were developed for calculations of the charge carrier concentrations, charge carrier mobilities, Hall coefficient, optical absorptance, and Fermi energy as functions of x, temperature, ionized donor and acceptor concentrations, and neutral defect concentrations.

Lehoczky, S. L.↗
NASA NTRS1981

A quantitative study of factors influencing lamellar eutectic morphology during solidification

The factors that influence the shape of the solid-liquid interface of a lamellar binary eutectic alloy are evaluated. Alloys of carbon tetrabromide and hexachloroethane which serve as a transparent analogue of lamellar metallic eutectics are used. The observed interface shapes are analyzed by computer-aided methods. The solid-liquid interfacial free energies of each of the individual phases comprising the eutectic system are measured as a function of composition using a 'grain boundary groove' technique. The solid-liquid interfacial free energy of the two phases are evaluated directly from the eutectic interface. The phase diagram for the system, the heat of fusion as a function of composition, and the density as a function of composition are measured. The shape of the eutectic interface is controlled mainly by the solid-liquid and solid-solid interfacial free energy relationships at the interface and by the temperature gradient present, rather than by interlamellar diffusion in the liquid at the interface, over the range of growth rates studied.

Kaukler, W. F. S.↗
DOE OSTI2022

Scalable Gas Diffusion Electrode Fabrication for Electrochemical CO 2 Reduction Using Physical Vapor Deposition Methods

Electrochemical CO 2 reduction (ECR) promises the replacement of fossil fuels as the source of feedstock chemicals and seasonal storage of renewable energy. While much progress has been made in catalyst development and electrochemical reactor design, few studies have addressed the effect of catalyst integration on device performance. Using a microfluidic gas diffusion electrolyzer, we systematically studied the effect of thickness and the morphology of electron beam (EB) and magnetron-sputtered (MS) Cu catalyst coatings on ECR performance. We observed that EB-Cu outperforms MS-Cu in current density, selectivity, and energy efficiency, with 400 nm thick catalyst coatings performing the best. Additionally, the superior performance of EB-Cu catalysts is assigned to their faceted surface morphology and sharper Cu/gas diffusion layer interface, which increases their hydrophobicity. Tests in a large-scale zero-gap electrolyzer yielded similar product selectivity distributions with an ethylene Faradaic efficiency of 39% at 200 mA/cm 2 , demonstrating the scalability for industrial ECR applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On a nonlocal Cahn–Hilliard model permitting sharp interfaces

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. In this work, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2004

An Investigation of Wave Propagations in Discontinuous Galerkin Method

Analysis of the discontinuous Galerkin method has been carried out for one- and two-dimensional system of hyperbolic equations. Analytical, as well as numerical, properties of wave propagation in a DGM scheme are derived and verified with direct numerical simulations. In addition to a systematic examination of the dissipation and dispersion errors, behaviours of a DG scheme at an interface of two different grid topologies are also studied. Under the same framework, a quantitative discrete analysis of various artificial boundary conditions is also conducted. Progress has been made in numerical boundary condition treatment that is closely related to the application of DGM in aeroacoustics problems. Finally, Fourier analysis of DGM for the Convective diffusion equation has also be studied in connection with the application of DG schemes for the Navier-Stokes equations. This research has resulted in five(5) publications, plus one additional manuscript in preparation, four(4) conference presentations, and three(3) departmental seminars, as summarized in part II. Abstracts of papers are given in part 111 of this report.

Hu, Fang Q.↗
NASA NTRS1995

Centrifugal and Axial Pump Design and Off-Design Performance Prediction

A meanline pump-flow modeling method has been developed to provide a fast capability for modeling pumps of cryogenic rocket engines. Based on this method, a meanline pump-flow code PUMPA was written that can predict the performance of pumps at off-design operating conditions, given the loss of the diffusion system at the design point. The design-point rotor efficiency and slip factors are obtained from empirical correlations to rotor-specific speed and geometry. The pump code can model axial, inducer, mixed-flow, and centrifugal pumps and can model multistage pumps in series. The rapid input setup and computer run time for this meanline pump flow code make it an effective analysis and conceptual design tool. The map-generation capabilities of the code provide the information needed for interfacing with a rocket engine system modeling code. The off-design and multistage modeling capabilities of PUMPA permit the user to do parametric design space exploration of candidate pump configurations and to provide head-flow maps for engine system evaluation.

Veres, Joseph P.↗
NASA NTRS1998

3-D Modeling of Directional Solidification of a Non-Dilute Alloy with Temperature and Concentration Fields Coupling via Materials Properties Dependence and via Double Diffusive Convection

Numerical simulation of the HgCdTe growth by the vertical Bridgman method was performed using FIDAP finite element code. Double-diffusive melt convection is analyzed, as the primary factor at controls inhomogeneity of the solidified material. Temperature and concentration fields in the model are also coupled via material properties, such as thermal and solutal expansion coefficients with the dependence on both temperature and concentration, and melting temperature evaluation from pseudobinary CdTe-HgTe phase diagram. Experimental measurements were used to obtain temperature boundary conditions. Parametric study of the melt convection dependence on the gravity conditions was undertaken. It was found, that the maximum convection velocity in the melt can be reduced under certain conditions. Optimal conditions to obtain a near flat solidified interface are discussed. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D calculations are compared with previous 2- D findings. A video film featuring 3-D melt convection will be presented.

Bune, Andris V.↗
DOE OSTI2022

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Simplifying computational workflows with the Multiscale Atomic Zeolite Simulation Environment (MAZE)

Zeolites, an important class of 3-dimensional nanoporous materials, have been widely explored for a variety of applications including gas storage, separations, and catalysis. As the properties of these aluminosilicate materials depend on a number of factors (e.g., framework topology, Si/Al ratio, extra-framework cations etc.), detailed experiments (e.g., catalytic properties, adsorption capacities etc.) are often limited to only a handful of materials. Computational methods have played an important role in (1) providing molecular level insights to rationalize experimental observations, and (2) screening large libraries of zeolites to identify promising candidates for experimental synthesis and validation. Different levels of theory and computational chemistry codes are necessary to describe the range of relevant phenomena such as adsorption (e.g., grand canonical Monte Carlo), diffusion (e.g., molecular dynamics), and chemical reactions (e.g., density functional theory). Manipulation of atomic structures, handling of input files, and developing robust workflows becomes quite cumbersome. To mitigate these challenges, we describe the development of the Multiscale Atomic Zeolite Simulation Environment (MAZE) – a Python package that simplifies zeolite-specific calculation workflows by providing a user-friendly interface for systematically manipulating zeolite structures

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

Structures Self-Assembled Through Directional Solidification

Nanotechnology has created a demand for new fabrication methods with an emphasis on simple, low-cost techniques. Directional solidification of eutectics (DSE) is an unconventional approach in comparison to low-temperature biomimetic approaches. A technical challenge for DSE is producing microstructural architectures on the nanometer scale. In both processes, the driving force is the minimization of Gibb's free energy. Selfassembly by biomimetic approaches depends on weak interaction forces between organic molecules to define the architectural structure. The architectural structure for solidification depends on strong chemical bonding between atoms. Constituents partition into atomic-level arrangements at the liquid-solid interface to form polyphase structures, and this atomic-level arrangement at the liquid-solid interface is controlled by atomic diffusion and total undercooling due to composition (diffusion), kinetics, and curvature of the boundary phases. Judicious selection of the materials system and control of the total undercooling are the keys to producing structures on the nanometer scale. The silicon-titanium silicide (Si-TiSi2) eutectic forms a rod structure under isothermal cooling conditions. At the NASA Glenn Research Center, directional solidification was employed along with a thermal gradient to promote uniform rods oriented with the thermal gradient. The preceding photomicrograph shows the typical transverse microstructure of a solidified Si-TiSi2 eutectic composition. The dark and light gray regions are Si and TiSi2, respectively. Preferred rod orientation along the thermal gradient was poor. The ordered TiSi2 rods have a narrow distribution in diameter of 2 to 3 m, as shown. The rod diameter showed a weak dependence on process conditions. Anisotropic etch behavior between different phases provides the opportunity to fabricate structures with high aspect ratios. The photomicrographs show the resulting microstructure after a wet chemical etch and a dry plasma etch. The wet chemical etches the silicon away, exposing the TiSi2 rods, whereas plasma etching preferentially etches the Si-TiSi2 interface to form a crater. The porous architectures are applicable to fabricating microdevices or creating templates for part fabrication. The porous rod structure can serve as a platform for fabricating microplasma devices for propulsion or microheat exchangers and for fabricating microfilters for miniatured chemical reactors. Although more work is required, self-assembly from DSE can have a role in microdevice fabrication.

Dynys, Frederick W.↗
DOE OSTI2026

Aluminum Based Solvent-Free Organic–Inorganic Hybrid Materials

In emerging materials, molecular hybrids are especially promising, as they have molecular level mixing of the organic and inorganic components, producing homogeneous materials without interfaces that can deteriorate properties. However, the current methods of manufacturing molecular hybrids are based on solution processing, which is impractical for bulk materials such as may be used for optically clear radiation and electromagnetic shielding components or photonics. Here we examine molecular hybrids composed of aluminum isopropoxide (AIP) and epoxy resins aiming to understand the molecular scale chemistry and manufacturability of these hybrid materials. DSCmonitored cure revealed the ideal cure temperature for these materials is 160−170 °C and demonstrated that an AIP concentration of 16.7 wt % maximizes the extent of reaction. Kinetic analysis of the curing reaction showed the Sestak−Berggren autocatalytic model is effective at temperatures over 140 °C but the reaction has diffusion limitations at a temperature of 120 °C. Mechanical testing with custom resin molds revealed a decrease in properties of the bulk samples with increasing AIP content due to an increase in defects but further testing with nanoindentation demonstrated comparable or improved mechanical properties of AIP-epoxy hybrids compared to epoxy resin with a standard hardener. Ultimately, this work lays the foundation for hardener-free epoxy-aluminum inorganic/organic hybrids and presents opportunities to expand on properties for specific applications such as thermal conductivity, optical clarity, and dielectric constant.

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

Utilization of 29Si MAS-NMR to Understand Solid State Diffusion in Energy Storage Materials

The properties of many solid-state materials arise from critical interfaces tied to the structure, morphology, and composition of the materials under study. For many materials, identifying components that may be invisible to diffraction techniques or other bulk sensitive techniques (i.e. inductively coupled plasma (ICP)), may cause important information to be overlooked. These can include grain boundary phases, nanoscale coatings, amorphous layers, or second phases that influence the materials environment. In this short review, the use of 29 Si MAS NMR as a local probe to detect silicon-containing phases in complex energy storage systems is explored with a focus is on silicon-containing materials and silicon electrodes. Examples highlighting the utility of 29 Si MAS NMR include 1) examining copper diffusion into silicon as a method to create 3 dimensional electrodes, 2) using Mg(II) electrolyte additives to create in-situ nanoscale silicide coatings to inhibit low voltage parasitic side reactions and extend calendar life, and 3) studying the lithiation reactions of passivated silicon on different time scales.

Dogan, Fulya↗
DOE OSTI2019

High Conductivity and Flexible Hybrid Solid State Electrolyte (Final Report: DOE-EE0007807)

The project objective is to develop flexible hybrid electrolyte with garnet nanofibers to achieve the following: (1) flexible, with greater mechanical strength (~ 10 MPa) and thermal stability than polymer electrolytes; (2) high room-temperature ionic conductivity, ~ 0.5 mS/cm; (3) stable interface with lithium metal and effective blocking of lithium dendrites at current densities up to 3 mA/cm 2 ; and (4) battery performance with Li S chemistry with an energy density of ≥ 450 Wh/kg (and ≥ 1000 Wh/L) and maintaining ≥ 80% of capacity up to 500 cycles. We synthesized garnet nanofibers, filled the porous region with polymer electrolyte, and characterized and tuned the flexible hybrid membrane properties. The flexible hybrid SSE microstructure was analyzed by leveraging the UMD AIMLab Focused Ion Beam Scanning Electron Microscopy (FIB-SEM) tools and integrated with electrochemical methods at the UMERC Battery Fabrication Lab to investigate the properties and stability with Li metal anode. We developed a fundamental understanding of the mechanism of Li ionic diffusion in garnet nanofibers and their mechanical properties, as well as these properties for hybrid garnet-fiber/polymer hybrids. Work progressed toward the study of the electrode assembly during electrochemical cycling of the anode.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Loci/STREAM Sharp Interface Model Validation of Cryogenic Propellant Tank Self-Pressurization

This paper documents a model validation effort for the prediction of cryogenic propellant tank selfpressurization physics using the Loci/STREAM sharp interface CFD model developed at NASA MSFC. Cryogenic fluid management (CFM) applications for long duration spaceflight missions require accurate modeling of heat transfer and thermodynamics. Therefore, detailed validation of those physics is important to ensure simulation accuracy. In the past, the MSFC Propulsion Fluid Dynamics Branch has validated and used the Loci/STREAM Volume of Fluid (VOF) module extensively for propellant slosh applications, which are dominated by fluid motion. For cases where the interface gas/liquid interface is essentially static, the Loci/STREAM sharp interface was developed for more rapid and efficient CFM analysis. In the present study, Loci/STREAM sharp interface simulations were compared with the NASA K-site self-pressurization experiment of a flight weight, partially full, LH2 tank in normal gravity. This study investigated several parametric effects and made a number of important observations on the performance of the Loci/STREAM sharp interface method: - First, a homogeneous model was derived to verify the CFD model. The homogeneous model assumes both gas and liquid phases are uniform and are at saturation temperature. The predicted pressure rise rates at 3.5 W/m2 and 2.0 W/m2 with a 49% fill level are all lower than the experimental measurement. This is expected due to the observed temperature stratification near the gas-liquid interface during testing. - The k-w and k-e turbulence models predict comparable pressurization rates, which are much lower than the experiment measurement. On the other hand, the pressure values and pressurization rates from both models are almost the same as those of the homogeneous model. The current turbulence model is too diffusive and is not capable of capturing the thermal stratification in the liquid and ullage. - When using the laminar flow model, the pressurization rate is found to be very close to the experimental value when a quasi-steady flow field is allowed to develop by continuously venting the tank prior to pressurization as was done during testing. On the other hand, without initial venting, the rate is not constant. Additionally, the development time to reach a quasi-steady flow field differed from the experiment due to the lack of turbulent mixing away from the gas-liquid interface that is expected in reality. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 15% of the experimental value for 49% fill level at 3.5 W/m 2 and is within 18% for 2.0 W/m 2 heat loading. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 18% of the experimental value for a 29% fill level and within 15% for an 83% fill level.

H. Q. Yang↗
DOE OSTI2021

In situ synchrotron diffraction and modeling of non-equilibrium solidification of a MnFeCoNiCu alloy

The solidification mechanism and segregation behavior of laser-melted Mn 35 Fe 5 Co 20 Ni 20 Cu 20 was firstly investigated via in situ synchrotron x-ray diffraction at millisecond temporal resolution. The transient composition evolution of the random solid solution during sequential solidification of dendritic and interdendritic regions complicates the analysis of synchrotron diffraction data via any single conventional tool, such as Rietveld refinement. Therefore, a novel approach combining a hard-sphere approximation model, thermodynamic simulation, thermal expansion measurement and microstructural characterization was developed to assist in a fundamental understanding of the evolution of local composition, lattice parameter, and dendrite volume fraction corresponding to the diffraction data. This methodology yields self-consistent results across different methods. Via this approach, four distinct stages were identified, including: (I) FCC dendrite solidification, (II) solidification of FCC interdendritic region, (III) solid-state interdiffusion and (IV) final cooling with marginal diffusion. It was found out that in Stage I, Cu and Mn were rejected into liquid as Mn 35 Fe 5 Co 20 Ni 20 Cu 20 solidified dendritically. During Stage II, the lattice parameter disparity between dendrite and interdendritic region escalated as Cu and Mn continued segregating into the interdendritic region. After complete solidification, during Stage III, the lattice parameter disparity gradually decreases, demonstrating a degree of composition homogenization. The volume fraction of dendrites slightly grew from 58.3 to 65.5%, based on the evolving composition profile across a dendrite/interdendritic interface in diffusion calculations. Postmortem metallography further confirmed that dendrites have a volume fraction of 64.7 ± 5.3% in the final microstructure.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Effect of Vacuum on Force Response of an Ultrasonic Penetrator

Introduction: The Apollo astronauts encountered higher than expected resistances when interacting with the lunar soil via the Apollo Lunar Surface Drill (ALSD) and the trenching tool. Reducing the force required to move tools or other mechanical components through regolith will impact many steps of the resource extraction process. Force reduction has been achieved in soil materials by imparting vibration to tooling interfaces such as a vibratory farming cultivator, a percussive scoop, and ultrasonically resonant penetrators. Vibration-assisted tools in granular media reduce interaction forces by fluidizing a volume around the tool, allowing the tool to progress through a dynamic (fluid) medium instead of a static (solid) medium. This work seeks to quantify ultrasonic vibration’s effect on the force response of a penetrator in lunar soil simulant in vacuum sufficient to be within the molecular flow regime of any disturbed gases. Methods: A custom vacuum chamber setup, CUBEvac, was designed and built to facilitate penetration testing in a high vacuum environment, for comparison to penetration behavior in ambient terrestrial environment. A two-stage pumping system (Agilent Triscroll 600 roughing pump, Agilent VHS-6 oil diffusion pump) reached chamber pressures of about 5x10-6 Torr with regolith simulant in place. Figure 1 is a schematic of the heart of the assembly (note the penetration drive mechanisms above the chamber feedthrough and the regolith simulant sample in the bottom are not shown). The penetration actuation stack was comprised of a stepper motor driving a lead screw to move the ultra-sonic probe vertically inside the chamber. Motion was coordinated with an Arduino Uno. GRC-3 lunar simulant was used for this set of experiments. Samples were prepared in a four-liter stain-less steel, cylindrical pot with an internal diameter of approximately 15.56 cm (6.125 in) and a depth of 19.37 cm (7 5/8 in) for testing. The soil was baked out prior to compaction preparation as a measure to reduce soil moisture which interfered with pump down capacity. The soil was not baked again if it was removed from the vacuum chamber, prepped, and immediately returned to the vacuum chamber for pump down. The soil was compacted using a 60 Hz vibration table with a surcharge of 34 kg place on top of the soil in the container. Prepared soil samples weighed approximately 6.5 kg (bulk density 1.895-1.934 g/cm3). Two probe end effectors were tested: A cone penetrometer (static only) with a nominal diameter of 12.7 mm (0.5 in) and a nominal height of 28.6 mm (1.125 in); and a vibrating cylindrical probe measuring 12.7 mm in diameter and 50.8 mm in effective length from the tip (Figure 2). The cylindrical probe vibrated resonantly at 20 kHz with an amplitude 23 μm. The cone penetration tests were conducted to assess potential soil behavior differences in vacuum. The cylinder probe tests were conducted as the primary subject of this investigation to assess force response in vacuum. For each test, a regolith simulant sample was load-ed and compacted in the chamber, which was then evacuated for roughly 18 hours to reach the lowest possible pressure (approximately 5x10-6 Torr for most tests). The probe was then moved to about 10 mm above the soil surface before being pushed to a depth of 50 mm for the cylinder probe tests and to a depth of 100 mm for the cone penetrometer tests, both at a speed of 2 mm/s. The simulant samples were prepared the same for all tests. Ideally, they would respond consistently to probe penetration under ambient and vacuum conditions. This was evaluated by measuring the resistance of representative prepared simulant beds with a standard cone penetrometer in both environments. Results and Discussion: The resistance of the simulant samples in the vacuum tests was consistently lower than in the ambient tests as determined by the cone penetration tests. Thus, the ambient and vacuum results cannot be compared directly; work is underway to de-confound and better correlate the data. Still, figures 3-6 show that probe penetration forces are lower overall in the vacuum environment. In both environments, resonant vibration of the probe provides two useful effects: It reduces the probe penetration force and smooths the force-depth curve, significantly reducing local maxima. These effects have implications for various potential applications, such as astronaut hand-tools, where benefits (reducing astronaut effort) outweigh the cost of the additional energy re-quired to generate vibration. These results demonstrate that resonantly vibrating tools can meaningfully reduce the penetration force required for excavation, probing, and drilling tools in simulated lunar regolith deposits under vacuum levels approaching those that will be experience on the Moon’s surface. Lunar-gravity, ambient environment tests are scheduled soon. The effects of realistic temperatures and temperature gradients and deeper vacuum remain to be tested.

E Rezich↗