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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1998

3-D Modeling of Double-Diffusive Convection During Directional Solidification of a Non-Dilute Alloy with Application to the HgCdTe Growth Under Microgravity Conditions

A numerical calculation for a non-dilute alloy solidification was performed using the FIDAP finite element code. For low growth velocities plane front solidification occurs. The location and the shape of the interface was determined using melting temperatures from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion with temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors dimensional rather then non-dimensional modeling was performed. Estimates of convection contributions for various gravity conditions was performed parametrically. For gravity levels higher then 1 0 -7 of earth's gravity it was found that the maximum convection velocity is extremely sensitive to gravity vector orientation and can be reduced at least by factor of 50% for precise orientation of the ampoule in the microgravity environment. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D modeling are compared with previous 2-D finding. A video film featuring melt convection will be presented.

Bune, Andris V.↗

DOE OSTI2023

A new approach for simulating inhomogeneous chemical kinetics

Abstract In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling $$\alpha$$ α - and $$\gamma$$ γ -radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Dynamics of Electric Polarization and Relaxation of Ions at Humid Calcite Surfaces

Mobile ions at mineral surfaces can respond to an applied electric field, adopting a new distribution that effectively represents polarization of the electrical double layer. When the field is released, the ions relax to their equilibrium distribution. In both cases, the dynamics are characteristic of the interface. However, current models of electrokinetic phenomena are not sufficiently robust to accurately predict collective ion dynamics at structurally and chemically specific mineral–water interfaces. Here, in this study, we use electrostatic force microscopy (EFM) to investigate the dynamics of ion relaxation at hydrated calcite (104) surfaces at controlled relative humidity (RH). Electrically biased probes are used to polarize the distributions of calcium and carbonate ions that are intrinsic to this interface across a range of RH values. Polarization kinetics are tracked by monitoring the tip–sample force gradient during charging, and EFM imaging is used to characterize the spatial relaxation dynamics after the applied field is released. Electrostatic finite element modeling of the sample/probe system across length-scales from nanometers to millimeters reproduces the observed stretched exponential charging response. Together, these results allow us to estimate the ion diffusivities at the interface across a wide range of RH values. These diffusivities increase by roughly 5 orders of magnitude as the RH is increased from 5 to 90%, highlighting the critical role of adsorbed water for surface ion solvation that enables ion mobility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1999

Effects of Gravity on the Double-Diffusive Convection During Directional Solidification of a Non-Dilute Alloy with Application to HgCdTe

A general 2-D and 3-D finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe in terrestrial and microgravity conditions. Verification of the 3-D model was undertaken by comparison with previously published results on convection in an inclined cylinder. For low growth velocities, plane front solidification occurs. The location and the shape of the interface were determined using melting temperatures obtained from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes a thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion due to the combined effects of composition and temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors, dimensional rather than non-dimensional modeling was performed. the predicted interface shape is in agreement with one obtained experimentally by quenching.

Bune, Andris V.↗

DOE OSTI2022

Molecular Beam Epitaxy (MBE) Growth of Model Cathodes to Study Interfacial Ion Diffusion

Lithium-ion batteries (LIBs) have been the focus of research for decades owing to their superior energy storage potential and wide range of applications. However, long-term stability issues still persist in LIB cathodes as the result of the formation of a solid electrolyte interface (SEI), structural transformations, or the loss of active cathode materials. In nanoscale cathodes, it is often impossible to isolate, let alone mitigate the causes of the observed capacity loss. Herein, a novel approach is presented to synthesizing LIB cathode model systems using thin-film molecular beam epitaxy (MBE). Specifically, 100 nm thick LiMn 2 O 4 (LMO) thin-film cathodes are grown on SrRuO 3 /SrTiO 3 (100) single-crystal substrates and characterized in their pristine and cycled states. The pristine thin films exhibit the electro-chemical cycling behavior and structural evolution previously seen in bulk LMO cathodes. The formation of surface-layer rocksalt MnO and spinel Mn 3 O 4 along with the mass-loss associated with the corrosion of active Mn ions is also reported. It is, therefore, demonstrated that these MBE-grown model systems can be used for detailed studies of their surface structures and the electrode/electrolyte interfacial evolution in LIB cathodes, leading to the development of materials for rechargeable batteries with higher capacity and better stability.

25 ENERGY STORAGE↗

DOE OSTI2021

Multi-physics Peridynamic Modeling of Damage Processes in Protective Coatings

Protective coatings can revolutionize the aerospace industry by providing a way to build a fuel-efficient and low-emission gas turbine engine operating at elevated temperature. Currently, significant challenges persist in the multi-physics modeling of damage processes in protective coatings exposed to the hot corrosive environment of gas turbine engine. Examples of such detrimental processes include an interface damage between different material layers forming protective coating, a formation of the through-the-thickness vertical cracks, a diffusion of oxygen and moisture through those cracks, an oxidation, and a delamination. In this paper, a feasibility of a multi-physics peridynamic model to provide insight to these detrimental processes has been demonstrated. 3-Dimensional peridynamic equations which are valid everywhere, including dynamically evolving discontinuities, have been presented. In conclusion, validation and demonstration have been shown for conduction heat transfer, reactive oxidation, and delamination of protective coating with a composition borrowed from the open literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Modeling Interfaces to Support Low-Level Waste Disposal System Performance Assessments - 20366

In low-level waste (LLW) disposal sites, interfaces between cementitious materials, used as waste forms and/or engineered barriers, and the surrounding soil or backfill material are often encountered. Reactions across these interfaces may lead to chemical and structural alteration of the cementitious and backfill materials that may ultimately affect long-term performance. For example, the ingress of carbon dioxide from soil gas into a waste tank concrete shell at the Hanford Site or Savannah River Site (SRS) or a low-level waste disposal vault at SRS leads to carbonation of cement hydration products (e.g. portlandite, calcium-silicate hydrate or C-S-H, and ettringite). The result of carbonation is a decrease in pH in the cement paste portion of the tank wall or vault concrete structure that could ultimately lead to de-passivation of the embedded structural steel and, eventually, to cracking. Cracking can lead to increased ingress of water into the structure, corresponding increased release of constituents of concern, and to subsequent increased transport of these constituents to the surrounding environment. The mobility of trace constituents (e.g., radionuclides of concern) in waste forms may increase in response to changes in pH, pore structure, and mineralogical gradients. Depletion of portlandite results in subsequent decalcification of C-S-H that leads to changes in the cement strength and may lead to structural failure. At the boundary between waste forms and concrete vaults, the migration of sulfate ions from the salt waste form (e.g., the saltstone waste form used for LLW disposal at the Savannah River Site in Saltstone Disposal Units, or SDUs) into the barrier or vault concrete has been predicted to result in expansive mineral phase formation (or 'sulfate attack'). Expansive mineral formation could result in cracking and potential loss of structural integrity in the vault concrete. In the performance assessment (PA) for the SRS SDUs, prediction of both the carbonation and sulfate attack fronts are critical to assessing the long-term performance of these LLW disposal vaults. A general purpose geochemical reactive transport model has been developed in LeachXS/ORCHESTRA to evaluate typical interfaces for LLW disposal environments, including SDUs at SRS. The model may be used to predict the long-term performance of interfaces between cementitious materials and saltstone or backfill materials with respect to primary phases by considering mass transport, cement chemistry, geochemical speciation, and multi-ionic diffusion over interfaces between different materials. The present study compares the performance and interactions between vault concrete and saltstone waste form for a representative SDU scenario. The materials were characterized using pH-dependence leaching test (EPA 1313) and semi-dynamic transport tests (EPA 1315) with test results used to parameterize the interface model. The development of a credible inter model provides an important basis for supporting future PAs where carbonation and sulfate attack are important aging and degradation mechanisms. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗

DOE OSTI2023

Characterizing pulse reverse-pulse uranium recovery in molten LiCl-KCl with intermittent electrochemical impedance spectroscopy

In this work, EIS was applied to characterize the temporal recovery of uranium from UCl 4 onto a tungsten electrode in molten LiCl-KCl salt between sets of pulse reverse-pulse (PRP) deposition steps. Impedance measurements were fitted to an equivalent circuit (EC) model of the electrode interface using Simplex optimization to calculate the double layer capacitance, charge transfer resistance, Warburg impedance, and CPE-restricted diffusion. Diffusion of uranium towards the electrode surface dominated the surface interactions for the first 80 cycles, followed by a transition to kinetic control for the remaining 520 cycles. A decrease in the Warburg impedance showed depletion of uranium from the bulk melt, which was confirmed to be a 52.5% reduction through ICP-OES. Reduction in the phasance from CPE-restricted diffusion suggested the formation of pores in the uranium deposits, which was confirmed through SEM characterization.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2022

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗

NASA NTRS2006

A Theoretical Solid Oxide Fuel Cell Model for System Controls and Stability Design

As the aviation industry moves towards higher efficiency electrical power generation, all electric aircraft, or zero emissions and more quiet aircraft, fuel cells are sought as the technology that can deliver on these high expectations. The Hybrid Solid Oxide Fuel Cell system combines the fuel cell with a microturbine to obtain up to 70 percent cycle efficiency, and then distributes the electrical power to the loads via a power distribution system. The challenge is to understand the dynamics of this complex multi-discipline system, and design distributed controls that take the system through its operating conditions in a stable and safe manner while maintaining the system performance. This particular system is a power generation and distribution system and the fuel cell and microturbine model fidelity should be compatible with the dynamics of the power distribution system in order to allow proper stability and distributed controls design. A novel modeling approach is proposed for the fuel cell that will allow the fuel cell and the power system to be integrated and designed for stability, distributed controls, and other interface specifications. This investigation shows that for the fuel cell, the voltage characteristic should be modeled, but in addition, conservation equation dynamics, ion diffusion, charge transfer kinetics, and the electron flow inherent impedance should also be included.

Kopasakis, George↗

NASA NTRS2008

A Theoretical Solid Oxide Fuel Cell Model for Systems Controls and Stability Design

As the aviation industry moves toward higher efficiency electrical power generation, all electric aircraft, or zero emissions and more quiet aircraft, fuel cells are sought as the technology that can deliver on these high expectations. The hybrid solid oxide fuel cell system combines the fuel cell with a micro-turbine to obtain up to 70% cycle efficiency, and then distributes the electrical power to the loads via a power distribution system. The challenge is to understand the dynamics of this complex multidiscipline system and the design distributed controls that take the system through its operating conditions in a stable and safe manner while maintaining the system performance. This particular system is a power generation and a distribution system, and the fuel cell and micro-turbine model fidelity should be compatible with the dynamics of the power distribution system in order to allow proper stability and distributed controls design. The novelty in this paper is that, first, the case is made why a high fidelity fuel cell mode is needed for systems control and stability designs. Second, a novel modeling approach is proposed for the fuel cell that will allow the fuel cell and the power system to be integrated and designed for stability, distributed controls, and other interface specifications. This investigation shows that for the fuel cell, the voltage characteristic should be modeled but in addition, conservation equation dynamics, ion diffusion, charge transfer kinetics, and the electron flow inherent impedance should also be included.

Kopasakis, George↗

DOE OSTI2024

Morphology selection in dealloying: A phase field study of the coupling among kinetic mechanisms

A phase field model is used to investigate how the corrosion rate and morphology of a dealloying metal is controlled by the kinetic rates of bulk solid diffusion, interfacial diffusivity, and interface dissolution. A specific focus is investigating the manner by which variations in these kinetic parameters synergistically influence the onset of a morphological instability at the solid-liquid interface, which is known to precede rapid dealloying corrosion towards highly topologically complex structures such as bicontinuous porous solids. This work demonstrates that this instability can be suppressed by either an increase of interdiffusion in the solid alloy, or by a decreasing rate of dissolution at the solid-liquid interface. When the instability is suppressed, the result is a relatively planar, passivated dealloying front. Importantly, however, coupled changes in kinetic parameters can also promote the morphological instability and sustain fine-scale ligament and pore formation. Finally, the relevance of the solid-liquid interfacial diffusion for determining the rate of dealloying and ligament formation is highlighted. By elucidating the mechanistic influence of bulk and interfacial kinetic factors on the dealloyed morphology, this work helps to rationalize the contrasting morphologies observed in various dealloying systems.

Diffusion↗

DOE OSTI2023

Effective Lanthanide Diffusivity through U-Zr Metallic Fuel

Fuel-cladding chemical interaction (FCCI) is one of the main concerns for the performance of U-Zr metallic fuels. Although FCCI results in complex phase transformations and morphological changes in both the fuel and cladding sides of the interface, the formation of a brittle layer in the cladding side of the interface (sometimes referred to as wastage) has one of the most significant impacts on the cladding mechanical integrity. This brittle layer is associated with the formation of intermetallic compounds between the cladding constituents Fe, Cr and the lanthanide fission product species. The formation of this cladding wastage layer has been the focus of recent efforts to develop a mechanistic modeling framework, including production of lanthanides in the fuel, transport through the fuel to the fuel-cladding interface, and phase transformation to form intermetallic phases. To calculate the rate of transport of lanthanides through the fuel, a multi-scale computational approach has been used. Atomistic calculations have been performed to determine the diffusivities of Nd, the most prevalent lanthanide species, through the solid fuel matrix and along pore surfaces. These diffusivities have been used as input parameters for a mesoscale model to determine an effective diffusion coefficient that accounts for the formation of porosity, its interconnection, and infiltration with bond sodium from the fuel-cladding gap. The effective diffusivity has been used in engineering-scale simulations using the BISON fuel performance code, and BISON simulations have been compared to representative validation cases from the EBR-II reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Fundamental study of gas species transport in the oxygen electrode of solid oxide fuel and electrolysis cells

A fundamental analysis of multicomponent gas transport models was performed in application to the oxygen electrodes of solid oxide cells. It is common practice to neglect the effect of pressure gradients within oxygen electrodes, even though a net molar flux at the electrolyte surface implies that a pressure gradient must exist. The influence of both Darcy velocity and Knudsen flux are considered in the context of ordinary (Fickian) diffusion, the dusty gas model, and the binary friction model. Comparisons between the models and different sets of assumptions are made via parametric studies on operating load, oxygen partial pressure, microstructural properties, and electrode thickness. Results show that the pressure gradient will have a significant impact on the oxygen concentration distribution and therefore the concentration overpotential. In electrolysis mode, pressure increases up to 1 atm are predicted, indicating that pressure at the electrode/electrolyte interface could contribute to electrode delamination. Additionally, it is found that Darcy's law is insufficient for calculating the pressure distribution without accounting for the flux due to Knudsen diffusion. Additionally, it is found that for the range of properties typical of oxygen electrodes, there is a negligibly small difference between the dusty gas model and binary friction model from a practical standpoint.

08 HYDROGEN↗

DOE OSTI2021

Multifidelity computing for coupling full and reduced order models

Hybrid physics-machine learning models are increasingly being used in simulations of transport processes. Many complex multiphysics systems relevant to scientific and engineering applications include multiple spatiotemporal scales and comprise a multifidelity problem sharing an interface between various formulations or heterogeneous computational entities. To this end, we present a robust hybrid analysis and modeling approach combining a physics-based full order model (FOM) and a data-driven reduced order model (ROM) to form the building blocks of an integrated approach among mixed fidelity descriptions toward predictive digital twin technologies. At the interface, we introduce a long short-term memory network to bridge these high and low-fidelity models in various forms of interfacial error correction or prolongation. The proposed interface learning approaches are tested as a new way to address ROM-FOM coupling problems solving nonlinear advection-diffusion flow situations with a bifidelity setup that captures the essence of a broad class of transport processes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Spatially Mapping the CO 2 Alkaline Sorbent Diffuse Microenvironment Using Operando Raman Spectroscopy

When designing a chemical process, the local balance of transport and kinetics, collectively referred to as the diffuse microenvironment, plays a critical role in performance but is difficult to directly observe. This work demonstrates a method of two-dimensional spatial chemical mapping of the diffuse microenvironment in the context of alkali metal hydroxide direct air capture of carbon dioxide using a custom operando gas-absorption flow cell along with confocal Raman spectroscopy. Notably, we observe the concentration boundary layer near the gas–liquid interface and elucidate the interplay of carbonate and bicarbonate ions within it while inferring local hydroxide depletion through continuum modeling. These first of their kind observations provide a technique to compare the performance of direct air capture solvents based on diffuse microenvironment dynamics while also providing metrics important for air contactor design such as boundary layer thickness. Overall, this work showcases a new experimental platform to study interfacial diffuse microenvironments in and outside of the field of direct air capture of carbon dioxide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1980

Liquid-phase electroepitaxy - Dopant segregation

A theoretical model is presented which accounts for the dopant segregation in liquid-phase electroepitaxy in terms of dopant transport in the liquid phase (by electromigration and diffusion), the growth velocity, and the Peltier effect at the substrate-solution interface. The contribution of dopant electromigration to the magnitude of the effective segregation coefficient is dominant in the absence of convection; the contribution of the Peltier effect becomes significant only in the presence of pronounced convection. Quantitative expressions which relate the segregation coefficient to the growth parameters also permit the determination of the diffusion constant and electromigration mobility of the dopant in the liquid phase. The model was found to be in good agreement with the measured segregation characteristics of Sn in the electroepitaxial growth of GaAs from Ga-As solutions. For Sn in Ga-As solution at 900 C the diffusion constant was found to be 4 x 10 to the -5 sq cm/s and the electromigration velocity (toward the substrate with a positive polarity 2 x 10 to the -5 cm/s current density of 10 A/sq cm.

Lagowski, J.↗