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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

The Implementation and Evaluation of the Emergency Response Dose Assessment System (ERDAS) at Cape Canaveral Air Station/Kennedy Space Center

The Emergency Response Dose Assessment System (ERDAS) is a system which combines the mesoscale meteorological prediction model RAMS with the diffusion models REEDM and HYPACT. Operators use a graphical user interface to run the models for emergency response and toxic hazard planning at CCAS/KCS. The Applied Meteorology Unit has been evaluating the ERDAS meteorological and diffusion models and obtained the following results: (1) RAMS adequately predicts the occurrence of the daily sea breeze during non-cloudy conditions for several cases. (2) RAMS shows a tendency to predict the sea breeze to occur slightly earlier and to move it further inland than observed. The sea breeze predictions could most likely be improved by better parameterizing the soil moisture and/or sea surface temperatures. (3) The HYPACT/REEDM/RAMS models accurately predict launch plume locations when RAMS winds are accurate and when the correct plume layer is modeled. (4) HYPACT does not adequately handle plume buoyancy for heated plumes since all plumes are presently treated as passive tracers. Enhancements should be incorporated into the ERDAS as it moves toward being a fully operational system and as computer workstations continue to increase in power and decrease in cost. These enhancements include the following: activate RAMS moisture physics; use finer RAMS grid resolution; add RAMS input parameters (e.g. soil moisture, radar, and/or satellite data); automate data quality control; implement four-dimensional data assimilation; modify HYPACT plume rise and deposition physics; and add cumulative dosage calculations in HYPACT.

Evans, Randolph J.↗
NASA NTRS1985

Pattern selection in a boundary-layer model of dendritic growth in the presence of impurities

Presently analyzed, in the context of a boundary-layer model, is the problem of pattern selection in dendritic growth in a situation where impurities are present in the undercooled liquid. It is found that the tip-velocity selection criterion that has been proposed recently for the geometrical model and the boundary-layer model of a pure substance can be extended, in a nontrivial way, to this more complex situation where two coupled diffusion fields (temperature and solute) determine the interface dynamics. This model predicts a sharp enhancement of tip velocity in good qualitative agreement with experiment. This agreement is consistent with the conjecture that a solvability condition can be used to determine the operating point of the dendrite in the full nonlocal problem.

Karma, A.↗
DOE OSTI2022

Self-mixing in microtubule-kinesin active fluid from nonuniform to uniform distribution of activity

Active fluids have applications in micromixing, but little is known about the mixing kinematics of systems with spatiotemporally-varying activity. To investigate, UV-activated caged ATP is used to activate controlled regions of microtubule-kinesin active fluid and the mixing process is observed with fluorescent tracers and molecular dyes. At low Péclet numbers (diffusive transport), the active-inactive interface progresses toward the inactive area in a diffusion-like manner that is described by a simple model combining diffusion with Michaelis-Menten kinetics. At high Péclet numbers (convective transport), the active-inactive interface progresses in a superdiffusion-like manner that is qualitatively captured by an active-fluid hydrodynamic model coupled to ATP transport. Results show that active fluid mixing involves complex coupling between distribution of active stress and active transport of ATP and reduces mixing time for suspended components with decreased impact of initial component distribution. This work will inform application of active fluids to promote micromixing in microfluidic devices.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Electric Field Effects on Water and Ion Structure and Diffusion at the Orthoclase (001)–Water Interface

Understanding the electrochemical properties of mineral–water interfaces tends to rely upon electrical double layer (EDL) models, but these models are based on the assumption that electrostatic equilibrium is constantly maintained. In reality, interfacial reactions, ion diffusion, and their electrochemical signatures are based in nonequilibrium conditions of locally or globally imbalanced electrical fields where current EDL models have limited purview. In this work, we performed molecular dynamics (MD) simulations of the orthoclase (001) surface in contact with a 1 M NaCl aqueous solution under various electric fields, to explore the interplay between EDL structure and dynamics when perturbed by electric fields of different direction and strength, by confinement, and by different distributions of structural surface charge. The simulations showed that confinement between two opposing (001) surfaces led to the development of an induced field when the applied field was perpendicular to the surfaces and, as a result, to ionic diffusion coefficients that were independent of electric field strength. In contrast, when the applied field was parallel to the surfaces, confinement resulted in ionic diffusion coefficients that were more strongly dependent on the magnitude of the electric field than in bulk water. Differences in the density and distribution of aluminol groups on the two surfaces had a significant impact on how the interfacial structure and dynamics varied in the presence of an electric field. Notably, these differences resulted in an electro-osmotic flow with opposite directions at the two surfaces under parallel applied electric field. Overall, the MD simulations highlighted the importance of considering atomic-level structure and heterogeneities when developing models of the electrochemical properties of mineral–water interfaces.

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

Insights into Na + Diffusion in Silicon Modules under Operating Conditions: Measuring Low Concentrations by D-SIMS

Sodium induced shunting under an electric field is a challenging reliability issue in crystalline Si solar modules. The source of this Potential-Induced Degradation of the Shunting type (PID-s) is well understood and its influence on power loss has been intensively studied based on phenomenological models on cell or module level relating the experimental power-loss to stressing parameters (time, temperature, voltage). However, little is known about the Na ion migration kinetics, responsible for PID on a microscopic level, and its quantitative relation to the efficiency degradation. In this paper we present our investigations of sodium ion migration in Ethylene-Vinyl Acetate (EVA) and silicon through Dynamic Secondary Ion Mass Spectroscopy (D-SIMS). Each sample was annealed at field relevant temperatures from 60-90 °C to address typical migration mechanisms of common PV installations. Analysis of the SIMS migration profiles revealed a diffusivity constant D 0,EVA = 0.09 ± 0.14 cm 2 /s and an activation energy E A,EVA = 0.85 ± .04 eV for Na in EVA and diffusivities higher than extrapolated literature values in silicon (D 0,Si = (3.03 ± 2.42)x10 -5 cm 2 /s, and E A,Si = 0.98 ± 0.02 eV). The new insight will be included in a drift-diffusion based degradation model accounting for the partition coefficient across all relevant interfaces. This model can assist in predicting PID-failure in the field based on the given model stack and the diffusion of Na + through each material. Furthermore, this tool can be used for process optimization as well as material selection significantly reducing the cost and time to validate a technology.

14 SOLAR ENERGY↗
DOE OSTI2021

Analytical modeling for redox flow battery design

Deeper market penetration of redox flow batteries requires optimization of the cell performance. Though important for optimization, detailed analytical solutions have not been developed for electrolyte flow, mass and charge transport, and reaction kinetics within redox flow batteries. To this end, here we present analytical solutions to active species concentration and over-potential based on advection-diffusion transport for ions and Bulter-Volmer model for interface reaction kinetics. The solutions were validated with results from a finite element model. These solutions were then applied to investigate the relationship between over-potential and state of charge, current density, reaction rate constant, flow velocity, diffusivity, total active species concentration, and electrode structure. Explicit formulas were identified for minimum activation over-potential and limiting current density as well as their dependence on electrolyte properties, operation conditions, and electrode structure. With our new mathematical formulas, this work provides a theoretical framework for flow battery design.

25 ENERGY STORAGE↗
DOE OSTI2025

First-Principles Studies of Tritium Species Diffusivity Across the Interfaces of Ni-NiZr Alloy-Zircaloy-4

A tritium producing assembly, also known as tritium producing burnable absorber rods (TPBARs), consists of a metal getter tube located between the cladding and γ-LiAlO2 pellets. The metal getter tube is composed of a nickel (Ni) layer coated on Zircaloy-4. The getter assembly is used to capture tritium (3H) species (mainly 3H2 and 3H2O) generated from γ-LiAlO2 pellets during irradiation. Exploring 3H species (3H2, 3H2O) dissociation on the Ni surface and diffusion in the Ni layer and across the interface of Ni-plated Zircaloy-4 getters can provide insights on tritium transport and retention in the pellets and the getter materials. . In FY25, based on the ideal interface model of Ni-Zircaloy-4 generated in FY24, we will further explore the diffusion pathways of 3H across such Ni-ZrNi-alloy layer-Zircalory-4 interfaces under different conditions, including with oxide/hydroxide clusters on the Ni layer, and impurities located in ZrNi-alloy and Zircalory-4 layers. Due to the size limitation of DFT simulation, we will separate the 3-layer system into 3 sub-systems: Ni-ZrNi interface, ZrNi alloy layer, ZrNi-alloy-Zircaloy-4 interface, and simulate the diffusion barriers for tritium.

diffusion barrier↗
NASA NTRS1976

Diffusion analysis for two-phase metal-matrix composite

Diffusion controlled filament matrix interaction in a metal matrix composite, where the filaments and matrix comprise a two phase binary alloy system, was mathematically modeled. The problem of a diffusion controlled, two phase moving interface by means of a one dimensional, variable grid, finite difference technique was analyzed. Concentration dependent diffusion coefficients and equilibrium solubility limits were used, and the change in filament diameter and compositional changes in the matrix were calculated as a function of exposure time at elevated temperatures. With the tungsten nickel (W-Ni) system as a model composite system, unidirectional composites containing from 0.06 to 0.44 initial filament volume fraction were modeled. Compositional changes in the matrix were calculated by superposition of the contributions from neighboring filaments. Alternate methods for determining compositional changes between first and second nearest neighbor filaments were also considered. The results show the relative importance of filament volume fraction, filament diameter, exposure temperature, and exposure time as they affect the rate and extent of filament matrix interaction.

Tenney, D. R.↗
DOE OSTI2023

Modeling the Charge and Exciton Distributions in Phosphorescent White Organic Light-Emitting Diodes

A comprehensive numerical model is developed to quantitatively analyze the charge and exciton distributions in phosphorescent white organic light emitting diodes (WOLEDs). The model encompasses bipolar charge drift-diffusion transport with doping-dependent charge mobilities, field-dependent charge hopping across multiple interfaces, exciton generation, diffusion, energy transfer, radiative decay, and non-radiative exciton annihilation. The model is validated by direct measurement of the exciton distribution in monochromatic phosphorescent OLEDs. It is then applied to complex, multilayer WOLEDs where it is used to accurately predict their intensity and current-dependent spectral shifts. The magnitudes of different exciton loss channels in WOLEDs are calculated and quantitatively compared. This model provides an alternative for the many experimental iterations commonly used to optimize the color balance of complex OLEDs. In this context, the model inspired the addition of a hole blocking layer in the WOLED to balance the charge conduction and reduce spectral shifts – a design that is verified by experiment.

42 ENGINEERING↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Development and application of sorption mass transfer models for fission product transport in TRISO fuel systems using BISON

Fission product mass transfer within and between TRISO particle and compact layers is an important phenomenon. It directly impacts fission product release predictions, which are used as source terms for safety and licensing calculations. Modeling fission product transport within layers and across bonded interfaces is relatively straightforward under the assumptions of isotropic Fickian diffusion, concentration continuity, and flux continuity. Modeling fission product transport across gaps and debonded layers is more difficult. Gaps are known to form between the buffer and inner PyC (IPyC), and debonding may occur between the IPyC and silicon carbide (SiC). A new mass transfer model was developed to provide a more accurate and robust fission product release calculation by accounting for interlayer sorptivity. One importance of the new model is to account for temperature and material changes across the gap based on sorption isotherm. This work presents the development of the sorption behavioral model and a fission product trapping model, and analyses the subsequent fission product diffusion behavior, specifically cesium, through TRISO particles and compact materials using the BISON fuel performance code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Three-dimensional modeling of hyphal fusion, branching, and nutrient transport in filamentous fungi

Fungi exhibit behaviors distinct from other microbes. Filamentous fungi grow by extending complex networks of branched filaments collectively referred to as the mycelium. These networks can expand over large distances and traverse low-nutrient areas by translocating nutrients through the filament network. This spatial characteristic makes filamentous fungi crucial for soil ecosystems, supporting stable microbial communities and promoting plant growth. However, simulating these behaviors is complex. The elongated nature of fungal compartments results in different mechanical interactions compared to the commonly modeled spherical bacteria. These detailed hyphal mechanics require specialized consideration and are often excluded from conventional fungal simulation packages. Additionally, the extensive fungal networks in nature demand computationally intensive simulations, necessitating high-performance algorithms. Therefore, realistic fungi simulations require specialized software. Here, we introduce a fungal modeling expansion to the high-performance biological modelling and interface exchange (bmx) software suite. bmx leverages adaptive mesh refinement in AMReX for chemical diffusion and incorporates a full mechanical model for bacterial cells, accelerated by GPUs. By extending bmx to model filamentous particles, we demonstrate the formation of complex filament networks through interactions like hyphal branching and fusion (anastomosis). We show that the networks produced match real-world fungal structures through various metrics. This work supports computational studies of fungal growth dynamics and can be adapted to investigate the growth of other filamentous structures in biology or materials science. The expanded-BMX package is open-sourced and is available online.

Cell mechanics↗
DOE OSTI2021

Precipitate growth kinetics under inhomogeneous concentration fields using a phase-field model

We investigate precipitation dynamics in the presence of a local solute gradient using phase-field simulations. During the homogenization heat treatment of the solidified Inconel 718 alloy, high Nb concentration within the Laves phases or at the core of the secondary arms results in Nb diffusion into the γ matrix. The volume fraction and spatial distribution of precipitation during subsequent annealing can be controlled by tailoring the Nb concentration gradient in the matrix during homogenization. We use a surrogate Ni-Fe-Nb alloy for Inconel 718 to explore the growth dynamics of δ precipitates related to the local Nb concentration levels. The simulations indicate that in the presence of a Nb concentration gradient the growth rate of δ precipitates is higher than in a matrix of uniform average Nb concentration. The higher growth rate is a result of the higher local thermodynamic driving force at the interface between the solute-rich matrix and the δ interface. Here, we propose a phenomenological model to describe the diffusion-controlled growth kinetics of the δ phase under a solute concentration gradient.

36 MATERIALS SCIENCE↗
NASA NTRS1997

Numerical Modeling of HgCdTe Solidification: Effects of Phase Diagram, Double-Diffusion Convection and Microgravity Level

Melt convection, along with species diffusion and segregation on the solidification interface are the primary factors responsible for species redistribution during HgCdTe crystal growth from the melt. As no direct information about convection velocity is available, numerical modeling is a logical approach to estimate convection. Furthermore influence of microgravity level, double-diffusion and material properties should be taken into account. In the present study, HgCdTe is considered as a binary alloy with melting temperature available from a phase diagram. The numerical model of convection and solidification of binary alloy is based on the general equations of heat and mass transfer in two-dimensional region. Mathematical modeling of binary alloy solidification is still a challenging numericial problem. A Rigorous mathematical approach to this problem is available only when convection is not considered at all. The proposed numerical model was developed using the finite element code FIDAP. In the present study, the numerical model is used to consider thermal, solutal convection and a double diffusion source of mass transport.

Bune, Andris V.↗
DOE OSTI2025

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Real-Time KMC Simulation of Vacancy-Mediated Intermixing in Au@Ag Octahedral Core–Cubic Shell Nanocrystals with Ab Initio-Guided Kinetics

Utilization of core–shell rather than monometallic nanocrystals (NCs) facilitates fine-tuning of NC properties for applications. However, compositional evolution via intermixing can degrade these properties prompting recent experimental studies. We develop an atomistic-level stochastic model for vacancy-mediated intermixing exploiting a formalism which allows incorporation at an ab initio density functional theory level of not just the thermodynamics of vacancy formation, but also relevant diffusion barriers for a vast number of possible local environments (in the core and in the shell, at the interface, and in the intermixed phase). This facilitates a predictive treatment and comprehensive understanding of intermixing on the relevant time scale (e.g., 10 1 –10 3 s). In contrast, previous modeling at the atomistic level utilized only unrealistic generic prescriptions of barriers or employed simplified continuum treatments. For Au@Ag octahedral core–cubic shell NCs, our modeling not only captures the experimentally observed rate or time scale for intermixing of ~100 s at 450 °C for 60 nm NCs, but also elucidates the underlying rate controlling processes and the effective intermixing barrier.

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

Non-isothermal effects on SO2 absorption by water droplets. I - Model development. II - Results and discussion

An analytic model of SO2 absorption in a falling water droplet is developed and a simulation of SO2 washout is performed. Nonisothermic effects on drop growth, droplet physical parameters, reaction rates, and multicomponent diffusion are treated in the model. The gas-liquid interface is assumed to be at equilibrium, and interfacial resistance is negligible. Raindrops are simulated as falling from a 2 km height through an atmospheric region containing SO2. The droplets decrease in size from evaporation and cooling, and their slightly basic pH aids SO2 absorption. The simulation indicates higher SO2 absorption at higher altitudes, and desorption may occur at ground level. Isothermal effects are concluded to be significant, and quantification of effects will depend on further modelling.

Reda, M.↗
DOE OSTI2022

Phase instabilities in austenitic steels during particle bombardment at high and low dose rates

Disruption of phase stability by energetic particle bombardment is a major challenge in designing advanced radiation-tolerant alloys and ion beam processing of nanocomposites. Particularly, ballistic dissolution susceptibility of different solute nanocluster species in alloys is poorly understood. Here, low dose rate neutron irradiations were conducted on a Fe-Cr-Ni based austenitic steel in the BOR-60 reactor (9.4 × 10 -7 dpa/s, 318 °C) followed by accelerated dose rate ion irradiations at multiple temperatures ( 10 -3 dpa/s, 380 – 420 °C). Using atom probe tomography, the stability of radiation-enhanced Cu-rich and radiation-induced Ni-Si-Mn-rich nanoclusters was evaluated. During neutron irradiation, Cu-rich clusters nucleated with their core concentrations progressively increasing with dose, while Ni-Si-Mn-rich clusters formed and evolved into G-phase precipitates. Ion irradiations dramatically altered the nanoclusters. Cu-rich clusters were ballistically dissolved, but Ni-Si-Mn-rich clusters remained stable and coarsened with dose at 400 and 420 °C, highlighting that different nanocluster species in a single microstructure can have innately distinct ballistic dissolution susceptibilities. Solute-specific recoil rates were incorporated into the Heinig precipitate stability model, which shows that in addition to radiation-enhanced diffusion, recovery from ballistic dissolution depends on solute concentration gradient near cluster interfaces. Finally, the combined experimental-modeling study quantified the critical temperatures and damage rates where ballistic dissolution dominates for each cluster species.

36 MATERIALS SCIENCE↗