Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffused interface model”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1999

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid interface shape and convection patterns during solidification of pure Al and Al-4 wt. % Cu alloy. Steady state calculations were performed for different gravity levels and orientations. Furthermore, transient modeling was undertaken to investigate effect of the solidification velocity. The furnace configuration used in this analysis is the proposed International Space Station Furnace. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. The Solidification model was adopted from previous work and was based on the finite element code FIDAP. Thermocouple data and quenched interface shape from a pure Al sample flown on the Life Sciences and Microgravity Spacelab (LMS) mission, July 1996 was used for model validation. Good agreement was obtained between the predicted interface shape and that measured from the quenched LMS sample. It was found that the imposed temperature boundary condition must be sufficiently smooth for the heat fluxes in the model to be self-consistent. The model predicted that alloy sample Al-Cu is more sensitive to variations in the gravity level compared to pure Al. For Al-Cu alloy, solute diffusivity is approximately 10(exp 4) times smaller than thermal diffusivity of pure Al. Hence, in a microgravity environment the weak convection has no measurable effect on the heat fluxes, but is still strong enough to affect concentration distribution. Since the alloy melting temperature is determined by concentration the interface shape depends on the level of convection. This was found not to be the case for pure metal.

Bune, Andris V.↗
DOE OSTI2020

Investigation of the discontinuous precipitation of U-Nb alloys via thermodynamic analysis and phase-field modeling

U-Nb’s discontinuous precipitation, $\gamma^{bcc}_{matrix}$, $\alpha^{orth}_{cellular}$, $\gamma\prime^{bcc}_{cellular}$, is intriguing in the sense that it allows formation and growth of the metastable $\gamma\prime$ phase during the course of its occurrence. Previous attempts to explain the thermodynamic origin of U-Nb’s discontinuous precipitation hypothesized that the energy of $\alpha$ forms an intermediate common tangent with the first potential of the double-well energy of $\gamma$ at the $\gamma\prime$ composition. In this work, we examine different possible mechanisms by which the discontinuous precipitation product in the U-Nb system can be stabilized. We put forward a mechanism by which the bulk free energy of the $\gamma$ can develop a non-monotonic dependency with composition. Additionally we posit that local contributions due to lattice mismatch between the $\alpha$ and $\gamma$ phases may be responsible for the generation of metastable states that may stabilized by thermodynamics as well as by kinetics. Our work suggests that local misfit strain tends to play a crucial role in the growth of the discontinuous precipitation product. Depending on the magnitude of strain developed at the $\alpha$/$\gamma\prime$ interfaces, either an increasing $\gamma\prime$ composition or a random distribution of $\gamma\prime$ compositions around the equiatomic value with respect to increasing temperature could be expected. Moreover, we show how it is possible to stabilize the discontinuous precipitation front through highly anisotropic and fast interface diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2002

Transport and Growth Kinetics in Microgravity Protein Crystal Growth

The dynamic coupling between mass transport and incorporation of growth units into the surface of a crystal growing from solution in microgravity is used to derive quantitative information on the crystal growth kinetics. To this end, new procedures for experiment preparation, interferometric data processing and model fitting have been developed. The use of experimental data from the bulk diffusive maw transport together with a model for steady state stagnant crystal growth allows the detailed quantitative understanding of the kinetics of both the concentration depletion zone around the crystal and the growth of the crystal interface. The protein crystal used in the experiment is shown to be growing in the mixed kinetic regime (0.2 x 10(exp -6) centimeters per second less than beta R/D less than 0.9 x 10(exp -6) centimeters per second).

Otalora, F.↗
NASA NTRS2003

Turbulent Mixing of Multiphase Flow

Thus we conduct numerical simulations of multiphase fluids stirred by two-dimensional turbulence to assess the possibility of self-similar drop size distribution in turbulence. In our turbulence simulations, we also explore the non-diffusive limit, where molecular mobility for the interface is vanishing. Special care is needed to transport the non-diffusive interface. Numerically, we use the particle level set method to evolve the interface. Instead of using the usual methods to calculate the surface tension force from the level set function, we reconstruct the interface based on phase- field modeling, and calculate the continuum surface tension forcing from the reconstructed interface.

Young, Y.-N.↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2020

A One‐Dimensional Model for Turbulent Mixing in the Benthic Biolayer of Stream and Coastal Sediments

Abstract In this paper, we develop and validate a rigorous modeling framework, based on Duhamel's Theorem, for the unsteady one‐dimensional vertical transport of a solute across a flat sediment‐water interface (SWI) and through the benthic biolayer of a turbulent stream. The modeling framework is novel in capturing the two‐way coupling between evolving solute concentrations above and below the SWI and in allowing for a depth‐varying diffusivity. Three diffusivity profiles within the sediment (constant, exponentially decaying, and a hybrid model) are evaluated against an extensive set of previously published laboratory measurements of turbulent mass transfer across the SWI. The exponential diffusivity profile best represents experimental observations and its reference diffusivity scales with the permeability Reynolds number, a dimensionless measure of turbulence at the SWI. The depth over which turbulence‐enhanced diffusivity decays is of the order of centimeters and comparable to the thickness of the benthic biolayer. Thus, turbulent mixing across the SWI may serve as a universal transport mechanism, supplying the nutrient and energy fluxes needed to sustain microbial growth, and nutrient processing, in the benthic biolayer of stream and coastal sediments.

Grant, Stanley B.↗
NASA NTRSDate not supplied

Delamination Depth Imaging for Quasi-Static and Fatigue Load Testing Using Passive Thermography

Passive thermography is used to monitor small increases in temperature resulting from delamination damage formation in composite hat-stiffened panels during both quasi-static and fatigue loading. It has been observed that rapid damage formation produces heating. This is composed of two heat generation components. The first component is an instantaneous response due to an irreversible thermoelastic strain release. The second component is mechanical heating, at the interface of failure, due to fracture damage. This produces a transient rise in temperature as the heat diffuses to the surface. The first component defines the thermal start time for the transient response. A one-dimensional thermal model is fitted to the data pixel by pixel, to produce imagery of the damage depth. The results are compared to ultrasonic and X-ray CT data.

Passive thermography↗
NASA NTRS1990

Process modelling for Space Station experiments

Examined here is the sensitivity of a variety of space experiments to residual accelerations. In all the cases discussed the sensitivity is related to the dynamic response of a fluid. In some cases the sensitivity can be defined by the magnitude of the response of the velocity field. This response may involve motion of the fluid associated with internal density gradients, or the motion of a free liquid surface. For fluids with internal density gradients, the type of acceleration to which the experiment is sensitive will depend on whether buoyancy driven convection must be small in comparison to other types of fluid motion, or fluid motion must be suppressed or eliminated. In the latter case, the experiments are sensitive to steady and low frequency accelerations. For experiments such as the directional solidification of melts with two or more components, determination of the velocity response alone is insufficient to assess the sensitivity. The effect of the velocity on the composition and temperature field must be considered, particularly in the vicinity of the melt-crystal interface. As far as the response to transient disturbances is concerned, the sensitivity is determined by both the magnitude and frequency of the acceleration and the characteristic momentum and solute diffusion times. The microgravity environment, a numerical analysis of low gravity tolerance of the Bridgman-Stockbarger technique, and modeling crystal growth by physical vapor transport in closed ampoules are discussed.

Alexander, J. Iwan D.↗
DOE OSTI2025

First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport

Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.

36 - MATERIALS SCIENCE↗
NASA NTRS1994

Axial and centrifugal pump meanline performance analysis

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) has been 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 is obtained from empirically derived correlations of loss to rotor specific speed. The rapid input setup and computer run time for the meanline pump flow code makes it an effective analysis and conceptual design tool. The map generation capabilities of the PUMPA code provide the information needed for interfacing with a rocket engine system modeling code.

Veres, Joseph P.↗
NASA NTRS1998

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗
NASA NTRS1997

Adsorption and solvation of ethanol at the water liquid-vapor interface: a molecular dynamics study

The free energy profiles of methanol and ethanol at the water liquid-vapor interface at 310K were calculated using molecular dynamics computer simulations. Both alcohols exhibit a pronounced free energy minimum at the interface and, therefore, have positive adsorption at this interface. The surface excess was computed from the Gibbs adsorption isotherm and was found to be in good agreement with experimental results. Neither compound exhibits a free energy barrier between the bulk and the surface adsorbed state. Scattering calculations of ethanol molecules from a gas phase thermal distribution indicate that the mass accommodation coefficient is 0.98, and the molecules become thermalized within 10 ps of striking the interface. It was determined that the formation of the solvation structure around the ethanol molecule at the interface is not the rate-determining step in its uptake into water droplets. The motion of an ethanol molecule in a water lamella was followed for 30 ns. The time evolution of the probability distribution of finding an ethanol molecule that was initially located at the interface is very well described by the diffusion equation on the free energy surface.

NASA Discipline Exobiology↗
NASA NTRS2005

Morphological Evolution of Directional Solidification Interfaces in Microgravity: An Analysis of Model Experiments Performed on the International Space Station

A series of experiments performed using the Pore Formation and Mobility Investigation (PFMI) apparatus within the glovebox facility (GBX) on board the International Space Station (ISS) has provided video images of the morphological evolution of a three-dimensional interface in a diffusion controlled regime. The experimental samples were prepared on ground by filling glass tubes, 1 cm ID and approximately 30 cm in length, with "alloys" of succinonitrile (SCN) and water in an atmosphere of nitrogen at 450 millibar pressure. The compositions of the samples processed and analyzed are 0.25,0.5 and 1.0 wt% water. Experimental processing parameters of temperature gradient and translation speed, as well as camera settings, were remotely monitored and manipulated from the ground Telescience Center (TSC) at the Marshall !3pace Flight Center. During the experiments, the sample was first subjected to a unidirectional melt back, generally at 10 microns per second, with a constant temperature gradient ahead of the melting interface. Following the melt back, the interface was allowed to stabilize before translation is initiated. The temperatures in the sample were monitored by six in situ thermocouples and the position is monitored by an optical linear encoder. For the experiments performed and analyzed, the gradients ranged from 2.5 - 3.3 K/mm and the initial pulling velocities ranged from 0.7 micron per second to 1 micron per second with subsequent transition velocities of up to 100 microns per second. The data provided by the PFMI for analysis includes near-real-time (NRT) video captured on the ground during the experiment runs, ISS Video Tape Recorder (VTR) data dumped from the VTR at the end of the experiment run and recorded on the ground, telemetry data including temperature and position measurements, and limited flight HI-8 tapes in 2 camera views of experiment runs for which tapes have been returned to the investigators from ISS. Because of limited down mass from the ISS, the majority of the initial analysis has been performed using the NRT and VTR video data but will be supplemented with the HI-8 video as it becomes available. hardware and procedures required to prepare samples for correlation to flight samples is described. Using this ground-based system, a series of experiments have been performed for direct comparison with the flight data. The results of these comparisons as well as implications to future microgravity experiments are presented and discussed. A ground-based thin-sample directional solidification system as well as all associated

Strutzenberg, Louise L.↗
NASA NTRS1993

Mixing and unmixedness in plasma jets 1: Near-field analysis

The flow characteristics in the near-field of a plasma jet are simulated with a two-fluid model. This model accounts for both gradient-diffusion mixing and uni-directional sifting motion resulting from pressure-gradient-body-force imbalance. This latter mechanism is believed to be responsible for the umixedness observed in plasma jets. The unmixedness is considered to be essentially a Rayleigh-Taylor kind instability. Transport equations are solved for the individual plasma and ambient gas velocities, temperatures and volume fractions. Empirical relations are employed for the interface transfers of mass, momentum and heat. The empirical coefficients are first established by comparison of predictions with available experimental data for shear flows. The model is then applied to an Argon plasma jet ejecting into stagnant air. The predicted results show the significant build-up of unmixed air within the plasma gas, even relatively far downstream of the torch. By adjusting the inlet condition, the model adequately reproduces the experimental data.

Ilegbusi, Olusegun J.↗
DOE OSTI2025

jaxhps: An elliptic PDE solver built with machine learning in mind

Elliptic partial differential equations (PDEs) can model many physical phenomena, such as electrostatics, acoustics, wave propagation, and diffusion. In scientific machine learning settings, a high-throughput PDE solver may be required to generate a training dataset, run in the inner loop of an iterative algorithm, or interface directly with a deep neural network. To provide value to machine learning users, such a PDE solver must be compatible with standard automatic differentiation frameworks, scale efficiently when run on graphics processing units (GPUs), and maintain high accuracy for a large range of input parameters. We have designed the jaxhps package with these use-cases in mind by implementing a highly efficient and accurate solver for elliptic problems with native hardware acceleration and automatic differentiation support.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Robust Dissimilar Metal Friction Welded Spool for Enhanced Capability for Steam Power Components (Final Program Report)

This project successfully demonstrated the feasibility and advantages of producing durable, friction-welded dissimilar-metal spools, both with and without transition pieces, as well as the use of advanced oxidation-protective coatings for steam power applications. The work encompassed the optimization and full-scale production of NFA tubing, overcoming manufacturing and processing challenges to achieve enhanced strength, creep, and fatigue properties validated by extensive microstructural characterization. Optimized welding techniques were developed and refined for joining dissimilar materials, resulting in defect-free welds optimized for both low- and high-pressure boiler circuits. Systematic evaluation of welding parameters, coupled with heat treatments, ensured reliable performance and minimized weld defects such as banding and cracking. Comprehensive mechanical and non-destructive testing, including tensile, bend, fatigue, and CT scanning, confirmed the superior integrity and durability of the friction-welded joints under representative service conditions. To address oxidation-related degradation, a range of protective alloy coatings—such as IN625, APMT, Duplex 2507, and FeCrAlY—were selected and evaluated for their ability to counteract oxide notch formation at dissimilar metal interfaces. These coatings demonstrated improved oxidation resistance and long-term stability under thermal cycling, thereby extending weld service life. Thermal cyclic testing and long-term studies captured the phase stability and diffusion behavior of coated and uncoated configurations. Validated models for low-cycle fatigue (LCF), creep, and thermal mechanical fatigue (TMF) were developed to predict component life and failure modes. These models confirmed that friction-welded configurations exhibited significantly enhanced performance compared to conventional gas-tungsten arc welding (GTAW), with up to 7× improvement in creep resistance and at least 2× enhancement in fatigue properties. The introduction of a tailored transition piece further reduced thermally induced strains, enhancing fatigue life. The results demonstrate a robust pathway for retrofitting existing steam fleet components and upgrading future equipment to deliver substantially improved cold-start cycling behavior and reduction of unplanned outages due to premature weld failures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗