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

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗
DOE OSTI2023

Heterostructure engineering in electrode materials for sodium-ion batteries: Recent progress and perspectives

Sodium-ion batteries (SIBs) have stepped into the spotlight as a promising alternative to lithium-ion batteries for large-scale energy storage systems. However, SIB electrode materials, in general, have inferior performance than their lithium counterparts because Na + is larger and heavier than Li + . Heterostructure engineering is a promising strategy to overcome this intrinsic limitation and achieve practical SIBs. We provide a brief review of recent progress in heterostructure engineering of electrode materials and research on how the phase interface influences Na + storage and transport properties. Efficient strategies for the design and fabrication of heterostructures (in situ methods) are discussed, with a focus on the heterostructure formation mechanism. The heterostructure's influence on Na + storage and transport properties arises primarily from local distortions of the structure and chemomechanical coupling at the phase interface, which may accelerate ion/electron diffusion, create additional active sites, and bolster structural stability. Finally, we offer our perspectives on the existing challenges, knowledge gaps, and opportunities for the advancement of heterostructure engineering as a means to develop practical, high-performance sodium-ion batteries.

25 ENERGY STORAGE↗
NASA NTRS1996

MGLab: An Interactive Multigrid Environment

MGLab is a set of Matlab functions that defines an interactive environment for experimenting with multigrid algorithms. The package solves two-dimensional elliptic partial differential equations discretized using either finite differences or finite volumes, depending on the problem. Built-in problems include the Poisson equation, the Helmholtz equation, a convection-diffusion problem, and a discontinuous coefficient problem. A number of parameters controlling the multigrid V-cycle can be set using a point-and-click mechanism. The menu-based user interface also allows a choice of several Krylov subspace methods, including CG, GMRES(k), and Bi-CGSTAB, which can be used either as stand-alone solvers or as multigrid acceleration schemes. The package exploits Matlab's visualization and sparse matrix features and has been structured to be easily extensible.

Bordner, James↗
NASA NTRS1982

On the cause of the flat-spot phenomenon observed in silicon solar cells at low temperatures and low intensities

A model that explains the flat-spot power loss phenomenon is presented. Evidence suggests that the effect is due to localized metallurgical interactions between the silicon substrate and the contact metallization. These reactions are shown to result in localized regions in which the P-N junction is destroyed and replaced with a metal semiconductor-like interface. The effects of thermal treatment, crystallographic orientation, junction depth, and metallization are shown along with a method of preventing the effect through the suppression of vacancy formation at the free surface of the contact metallization. Data indicating the effectiveness of a TiN diffusion barrier in preventing the effect are also given.

Weizer, V. G.↗
NASA NTRS1988

Microlaminate composites as thermal barrier coatings

Thick multiple-layered Ni/NiCoCrAlY and Ti/CoCrAlY microlaminate composites are explored as thermal barrier coatings. The method of fabrication of these laminates and the measurement technique used to determine thermal diffusivity and thermal conductivity of these coatings are discussed. Results indicate that the thermal conductivity of the laminate composite perpendicular to the laminate plane decreases with increasing number of interfaces; the drop in thermal conductivity is likely to be associated with interfaces which act as a barrier to the transfer of heat across them; the variation in thermal conductivity with the number of interfaces is not linear.

Radhakrishna, M. C.↗
DOE OSTI2021

Microstructural and infiltration properties of woven preforms during chemical vapor infiltration

Interface-resolved direct numerical simulations (DNSs) of chemical vapor infiltration (CVI) have been performed over a range of furnace-operating conditions (Thiele moduli) and for practical woven preform geometries. A level-set method is used to resolve the geometry of the initial preform at tow scale. The interface between the vapor and solid phase is then evolved in time through the entire CVI densification cycle, fully resolving the time-varying topology between the two phases. In contrast to previous level-set methods for CVI simulation, the physical reaction and diffusion processes govern the level-set movement in the current approach. The surface deposition kinetics is described by the usual one-step model. In this paper, the DNS data are used to study the evolving porosity, surface-to-volume ratio, and flow infiltration properties (permeability and effective diffusivities). Additionally, comparisons are made to popularly-assumed structure functions and the standard, Kozeny–Carmen porous media model commonly employed in modeled CFD simulations of CVI. The virtual DNS experiments reveal a Thiele modulus and preform geometry (fabric layup) dependence which the existing microstructural and infiltration models are not able to describe throughout the entire densification process. The DNS-based, woven geometry-specific correlations can be applied directly to mean-field, furnace-scale CFD simulations.

36 MATERIALS SCIENCE↗
NASA NTRS1982

On the cause of the flat-spot phenomenon observed in silicon solar cells at low temperatures and low intensities

A model is presented that explains the "flat-spot" (FS) power loss phenomenon observed in silicon solar cells operating deep space (low temperature, low intensity) conditions. Evidence is presented suggesting that the effect is due to localized metallurgical interactions between the silicon substrate and the contact metallization. These reactions are shown to result in localized regions in which the PN junction is destroyed and replaced with a metal-semiconductor-like interface. The effects of thermal treatment, crystallographic orientation, junction depth, and metallurization are presented along with a method of preventing the effect through the suppression of vacancy formation at the free surface of the contact metallization. Preliminary data indicating the effectiveness of a TiN diffusion barrier in preventing the effect are also given.

Weizer, V. G.↗
NASA NTRS1982

On the cause of the flat-spot phenomenon observed in silicon solar cells at low temperatures and low intensities

A model is presented that explains the 'flat-spot' power-loss phenomenon observed in silicon solar cells operating under deep space (low temperature, low intensity) conditions. Evidence is presented suggesting that the effect is due to localized metallurgical interactions between the silicon substrate and the contact metallization. These reactions are shown to result in localized regions in which the PN junction is destroyed and replaced with a metal-semiconductor-like interface. The effects of thermal treatment, crystallographic orientation, junction depth, and metallization are presented along with a method of preventing the effect through the suppression of vacancy formation at the free surface of the contact metallization. Preliminary data indicating the effectiveness of a TiN diffusion barrier in preventing the effect are also given.

Weizer, V. G.↗
NASA NTRS1987

A survey of nested grid techniques and their potential for use within the MASS weather prediction model

A survey of various one- and two-way interactive nested grid techniques used in hydrostatic numerical weather prediction models is presented and the advantages and disadvantages of each method are discussed. The techniques for specifying the lateral boundary conditions for each nested grid scheme are described in detail. Averaging and interpolation techniques used when applying the coarse mesh grid (CMG) and fine mesh grid (FMG) interface conditions during two-way nesting are discussed separately. The survey shows that errors are commonly generated at the boundary between the CMG and FMG due to boundary formulation or specification discrepancies. Methods used to control this noise include application of smoothers, enhanced diffusion, or damping-type time integration schemes to model variables. The results from this survey provide the information needed to decide which one-way and two-way nested grid schemes merit future testing with the Mesoscale Atmospheric Simulation System (MASS) model. An analytically specified baroclinic wave will be used to conduct systematic tests of the chosen schemes since this will allow for objective determination of the interfacial noise in the kind of meteorological setting for which MASS is designed. Sample diagnostic plots from initial tests using the analytic wave are presented to illustrate how the model-generated noise is ascertained. These plots will be used to compare the accuracy of the various nesting schemes when incorporated into the MASS model.

Koch, Steven E.↗
DOE OSTI2022

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗
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

Scaling Laws of Exciton Recombination Kinetics in Low Dimensional Halide Perovskite Nanostructures

Carrier recombination is a crucial process governing the optical properties of a semiconductor. Although various theoretical approaches have been utilized to describe carrier behaviors, a quantitative understanding of the impact of defects and interfaces in low dimensional semiconductor systems is still elusive. Here, we develop a model system consisting of chemically tunable, highly luminescent halide perovskite nanocrystals to illustrate the role of carrier diffusion and material dimensionality on the carrier recombination kinetics and luminescence efficiency. Our advanced synthetic methods provide a well-controlled colloidal system consisting of nanocrystals with different aspect ratios, halide compositions, and surface conditions. Using this system, we reveal the scaling laws of photoluminescence quantum yield and radiative lifetime with respect to the aspect ratio of nanocrystals. The scaling laws derived herein are not only a phenomenological observation but proved a powerful tool disentangling the carrier dynamics of microscopic systems in a quantitative and interpretable manner. The investigation of our model system and theoretical formulation bring to light the dimensionality, as a hidden constraint on carrier dynamics, and identify the diffusion length as an important parameter that distinguishes nanoscale and macroscale carrier behaviors. Furthermore, the conceptual distinction in carrier dynamics in different dimensionality regimes informs new design rules for optical devices where complex microstructures are involved.

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

Time-dependent MOS breakdown

A general model for time-dependent breakdown in metal-oxide-silicon (MOS) structures is developed and related to experimental measurements on samples deliberately contaminated with Na. A statistical method is used for measuring the breakdown probability as a function of log time and applied field. It is shown that three time regions of breakdown can be explained respectively in terms of silicon surface defects, ion emission from the metal interface, and lateral ion diffusion at the silicon interface.

Li, S. P.↗
DOE OSTI2021

Magnetohydrodynamic Control of Interfacial Degradation in Lithium-Ion Batteries for Fast Charging Applications

Interfacial anodic degradation in graphitic materials under fast charging conditions causes severe performance loss and safety hazard in lithium ion batteries. We present a novel method for minimizing the growth of these aging mechanism by application of an external magnetic field. Under magnetic field, paramagnetic lithium ions experience a magnetohydrodynamic force, which rotates the perpendicularly diffusing species and homogenizes the ionic transport. This phenomenon minimizes the overpotential hotspots at the anode/separator interface, consequently reducing SEI growth, lithium plating, and interfacial fracture. In situ electrochemical measurements indicate an improvement in capacity for lithium cobalt oxide/graphite pouch cell (20 mAh) charged from 1–5 C under an applied field of 1.8 kG, with a maximum capacity gain of 22% at 5C. Post-mortem FE-SEM and EDS mapping shows that samples charged with magnetic field have a reduced lithium deposition at 3C and a complete suppression of interfacial fracture at 5C. At 5C, a 24% reduction in the lithium content is observed by performing XPS on the anodic interfacial film. Finally, fast charging performance under variable magnetic field strengths indicate a saturation behavior in capacity at high fields (>2 kG), thereby limiting the field and consequent energy requirements to obtain maximum capacity gain under extreme conditions.

electromagnetism↗
NASA NTRS2002

Local Discontinuous Galerkin Methods for Partial Differential Equations with Higher Order Derivatives

In this paper we review the existing and develop new continuous Galerkin methods for solving time dependent partial differential equations with higher order derivatives in one and multiple space dimensions. We review local discontinuous Galerkin methods for convection diffusion equations involving second derivatives and for KdV type equations involving third derivatives. We then develop new local discontinuous Galerkin methods for the time dependent bi-harmonic type equations involving fourth derivatives, and partial differential equations involving fifth derivatives. For these new methods we present correct interface numerical fluxes and prove L(exp 2) stability for general nonlinear problems. Preliminary numerical examples are shown to illustrate these methods. Finally, we present new results on a post-processing technique, originally designed for methods with good negative-order error estimates, on the local discontinuous Galerkin methods applied to equations with higher derivatives. Numerical experiments show that this technique works as well for the new higher derivative cases, in effectively doubling the rate of convergence with negligible additional computational cost, for linear as well as some nonlinear problems, with a local uniform mesh.

Yan, Jue↗
NASA NTRS1987

Oxygen-18 tracer study of the passive thermal oxidation of silicon

This work focuses on the thermal oxidation of silicon near 1273 K using the double-tracer oxidation method. The results confirm that oxidation occurs by the transport of electrically neutral nonnetwork oxygen through the interstitial space of the vitreous silica (v-SiO2) scale. Simultaneously, self- (or isotopic-) diffusion occurs in the network, resulting in characteristic isotopic fraction distributions near the gas-scale interface. The self-diffusion coefficients calculated from these profiles agree with those reported for tracer diffusion in v-SiO2, and the diffusion coefficient calculated from the scale growth is consistent with reported O2 permeation data. An important parameter that describes the double-oxidation behavior is a ratio relating the scale thickness grown during the second oxidation, the network self-diffusion coefficient for oxygen, and the time of the second oxidation.

Cawley, J. D.↗
NASA NTRS1992

Nature of fluid flows in differentially heated cylindrical container filled with a stratified solution

Semiconductor crystals such as Hg(1-x)Cd(x)Te grown by unidirectional solidification Bridgmann method have shown compositional segregations in both the axial and radial directions. Due to the wide separation between the liquidus and the solidus of its pseudobinary phase diagram, there is a diffusion layer of higher HgTe content built up in the melt near the melt-solid interface which gives a solute concentration gradient in the axial direction. Because of the higher thermal conductivity in the melt than that in the crystal there is a thermal leakage through the fused silica crucible wall near the melt-solid interface. This gives a thermal gradient in the radial direction. Hart (1971), Thorpe, Hutt and Soulsby (1969) have shown that under such condition a fluid will become convectively unstable as a result of different diffusivities of temperature and solute. It is quite important to understand the effects of this thermosolute convection on the compositional segregation in the unidirectionally solidified crystals. To reach this goal, we start with a simplified problem. We study the nature of fluid flows of a stratified solution in a cylindrical container with a radial temperature gradient. The cylindrical container wall is considered to be maintained at a higher temperature than that at the center of the solution and the solution in the lower gravitational direction has higher solute concentration which decrease linearly to a lower concentration and then remain constant to the top of the solution. The sample solution is taken to be salt water.

Wang, Jai-Ching↗