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

Modeling discontinuous dynamic recrystallization containing second phase particles in magnesium alloys utilizing phase field method

Pre-existing second phase particles can induce particle stimulated nucleation (PSN) and pinning effect on grain boundaries during DRX process of magnesium alloys. The interaction among pinning effect, PSN mechanism and strain-induced grain boundary migration (SIBM) nucleation mechanism imposes challenge in quantitative study of microstructure evolution. A phase field model that describes discontinuous dynamic recrystallization process considering second phase particles (PF-DDRX-SP) which mainly distribute along grain boundaries is developed. The second phase order parameter with diffusion interface is introduced to the free energy function. With considering the pinning effect of second phase particles on grain boundaries, the grain boundary energy couples the interaction between second phase and grains. The effect of second phase on dislocation density evolution is considered in the model, and the PSN and SIBM mechanism are both executed by limiting the nucleation position to the grain boundary and the second phase interface. For the validation, the initial topology consistent with initial microstructure is applied to the simulation of DRX kinetics and flow behavior of AZ80 magnesium alloy, and great reliability and accuracy of developed model is indicated. Furthermore, through the model, the microstructure evolution with different particle size and with different volume fraction of round particle during DRX process were simulated. Finally, the results show that particles with higher dispersion reduce the number of DRX nuclei, which indicate that second phase particles suppress the nucleation at grain boundaries because of the pinning effect. More importantly, the proposed model could also quantitatively predict the relationships between the parameters of second phase particles and DRX behaviors, and enable to optimize the initial second phase structure in a uniform grain structure during thermomechanical process.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Improved analytic modeling of the linear Rayleigh–Taylor instability with plasma transport

The Rayleigh–Taylor (R–T) instability with transport effects has long been studied analytically in the neutral fluid regime. How this picture may differ in plasmas has been given less consideration, however. Leveraging previous numerical simulations with plasma viscosity and diffusion of R–T at a binary plasma interface, we show how plasma transport and kinetic effects alter the linear R–T dispersion relation. Our insight informs a much improved analytical dispersion relation, which better matches the simulation data than previously established models. Additionally, we demonstrate the importance of employing asymptotically correct viscosity coefficients in fluid codes for accurately capturing plasma R–T growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A dissolution model of alite coupling surface topography and ions transport under different hydrodynamics conditions at microscale

Portland cement is the most produced material in the world. The hydration process of cement consists of a group of complex chemical reactions. In order to investigate the mechanism of cement hydration, it is vital to study the hydration of each phase separately. An integrated model is proposed in this paper to simulate the dissolution of alite under different hydrodynamic conditions at microscale, coupling Kinetic Monte Carlo model (KMC), Lattice Boltzmann method (LBM) and diffusion boundary layer (DBL). The dissolution of alite is initialised with KMC. Two Multiple-relaxation-time (MRT) LB models are used to simulate the fluid flow and transport of ions, respectively. For solid-liquid interface, DBL is adapted to calculate the concentration gradient and dissolution flux. The model is validated with experiment from literature. The simulation results show good agreements with the results published in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Migration of inclusions in a matrix due to a spatially varying interface energy

An interfacial energy can be a function of a bulk field such as temperature or electric field. Here we find that in a system with a gradient in temperature or electric potential, the resulting variation in interfacial energy can induce a particle to migrate by either surface or bulk diffusion. For a circular particle under a constant unidirectional gradient in the bulk field, the field dependence of the interfacial energy induces a circular particle to move as a circle at a constant velocity along the direction from higher to lower interfacial energy. A linear stability analysis of this steady state migration suggests that perturbations will damp out as time evolves, and thus under these conditions a migrating circular particle is morphologically stable. Other spatial distributions of interface energy can lead to the distortion of an initially circular shaped particle during migration. A phase field model is developed that captures these distortions and verifies the theoretical results mentioned above.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thin-Film Paradigm to Probe Interfacial Diffusion during Solid-State Metathesis Reactions

We highlight a paradigm for studying complex reaction mechanisms that guide the synthesis of materials. Thin films of FeCl 2 and Na 2 S 2 were deposited to study the metathesis reaction to form FeS 2 and NaCl. In situ X-ray reflectivity was used to monitor the interface between materials, which revealed a slow, impeded reaction at high temperatures as compared to previous studies using powder samples. AC impedance and X-ray photoelectron spectroscopy provided insight into distribution of elements and conductivity of the phases present during the reaction, and phase-field modeling was used to elucidate the diffusion of ions throughout the thin-film bilayers. Here, the use of thin-film bilayers provides a simplified system to study solid-state metathesis reactions and highlights the complexity of diffusion at solid-state interfaces.

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

Simulation of Impedance Changes with Aging in Lithium Titanate-based Cells Using Physics-Based Dimensionless Modeling

Quantifying aging effects in lithium-ion cells with chemistries that have a flat open circuit potential is challenging. We implement a physics-based electrochemical model to track changes in the electrochemical impedance response of lithium titanate-based cells. Frequency domain equations of a pseudo two-dimensional model are made dimensionless, and the corresponding non-dimensional parameters are estimated using a Levenberg-Marquardt routine. The model weighs the relative contributions of changes in diffusion, ionic conduction within the electrolyte phase against solid phase electronic conduction towards cell aging. Solid-phase diffusion, charge transfer resistance and double layer capacitance at the solid-liquid interface are accounted for in the particle impedance. The estimation routine tracks dimensionless parameters using accelerated cycling data from full cells over 1000 cycles. The model can be deployed within a short time for state estimation using physics-based models without requiring prior knowledge of the battery chemistry, format, or capacity.

25 ENERGY STORAGE↗
DOE OSTI2023

Resolving experimental biases in the interpretation of diffusion experiments with a user-friendly numerical reactive transport approach

The reactive transport code CrunchClay was used to derive effective diffusion coefficients (D e ), clay porosities (ε), and adsorption distribution coefficients (K D ) from through-diffusion data while considering accurately the influence of unavoidable experimental biases on the estimation of these diffusion parameters. These effects include the presence of filters holding the solid sample in place, the variations in concentration gradients across the diffusion cell due to sampling events, the impact of tubing/dead volumes on the estimation of diffusive fluxes and sample porosity, and the effects of O-ring-filter setups on the delivery of solutions to the clay packing. Doing so, the direct modeling of the measurements of (radio)tracer concentrations in reservoirs is more accurate than that of data converted directly into diffusive fluxes. While the above-mentioned effects have already been described individually in the literature, a consistent modeling approach addressing all these issues at the same time has never been described nor made easily available to the community. A graphical user interface, CrunchEase, was created, which supports the user by automating the creation of input files, the running of simulations, and the extraction and comparison of data and simulation results. While a classical model considering an effective diffusion coefficient, a porosity and a solid/solution distribution coefficient (D e –ε–K D ) may be implemented in any reactive transport code, the development of CrunchEase makes it easy to apply by experimentalists without a background in reactive transport modeling. CrunchEase makes it also possible to transition more easily from a D e –ε–K D modeling approach to a state-of-the-art process-based understanding modeling approach using the full capabilities of CrunchClay, which include surface complexation modeling and a multi-porosity description of the clay packing with charged diffuse layers.

58 GEOSCIENCES↗
DOE OSTI2025

Multi-scale modeling of wastage layer formation in metallic fuel cladding

Fuel-cladding chemical interaction (FCCI) is a major concern for U-Zr metallic fuels' performance, primarily due to the formation of a brittle layer (wastage) in the cladding. This brittle layer, resulting from intermetallic compounds between cladding constituents Fe, Cr, and lanthanide fission products, significantly impacts the cladding's mechanical integrity. Recent efforts focus on developing a mechanistic modeling framework to understand lanthanide production, transport to the fuel-cladding interface, and phase transformation to intermetallic phases. A multi-scale computational approach has been used to calculate lanthanide transport rates, with atomistic calculations determining Nd diffusivities through the solid fuel matrix and along pore surfaces. These diffusivities inform a mesoscale model to determine an effective diffusion coefficient, accounting for porosity and infiltration with bond sodium. This effective diffusivity is used in engineering-scale simulations via the BISON fuel performance code, which has been validated against EBR-II and FFTF reactor experiments.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Decoupling Li out-diffusion and surface diffusion in the lithiation-assisted epitaxial growth of lithium tungstate

Lithiation-assisted epitaxy offers a flexible and robust approach for synthesizing high-quality Li-containing materials and interfaces with precise control. Here, in this study, we use lithium tungstate (Li x WO 3+x/2 , where x = 0 to 2) as a model system to investigate the intertwined effects of Li out-diffusion-induced compositional changes and surface-diffusion-induced morphological changes. By systematically varying synthesis and processing conditions, we uncover their impact on lithium tungstate film formation. Comprehensive characterizations, including X-ray diffraction, atomic force microscopy, X-ray photoemission spectroscopy and time-of-flight secondary ion mass spectrometry, reveal that low-temperature growth (< 300 °C) followed by high-temperature annealing yields continuous lithium tungstate films with significantly reduced surface roughness. In contrast, high-temperature deposition (≥ 300 °C) accelerates surface diffusion and Li out-diffusion, leading to island formation. Furthermore, in situ scanning transmission electron microscopy demonstrates the beam sensitivity of Li 2 WO 4 and reveals a phase transition from Li 2 WO 4 to LiWO 3.5 under prolonged electron beam exposure. These findings deepen our understanding of how to control composition and morphology of Li-containing films, providing valuable insights for the design and integration of energy materials.

Shi, Jueli [Pacific Northwest National Laboratory ↗
DOE OSTI2020

Interface learning in fluid dynamics: Statistical inference of closures within micro–macro-coupling models

Many complex multiphysics systems in fluid dynamics involve using solvers with varied levels of approximations in different regions of the computational domain to resolve multiple spatiotemporal scales present in the flow. The accuracy of the solution is governed by how the information is exchanged between these solvers at the interface and several methods have been devised for such coupling problems. In this article, we construct a data-driven model by spatially coupling a microscale lattice Boltzmann method (LBM) solver and macroscale finite difference method (FDM) solver for reaction-diffusion systems. The coupling between the micro-macro solvers has one to many mapping at the interface leading to the interface closure problem, and we propose a statistical inference method based on neural networks to learn this closure relation. The performance of the proposed framework in a bifidelity setting partitioned between the FDM and LBM domain shows its promise for complex systems where analytical relations between micro-macro solvers are not available.

42 ENGINEERING↗
DOE OSTI2022

Production and diffusion of H 2 O 2 during the interaction of a direct current pulsed atmospheric pressure plasma jet on a hydrogel

The interaction of cold atmospheric pressure plasma jets with hydrogels has been used as a model system to study the interaction of plasmas with tissues. In this study, we analyze the diffusion of reactive oxygen species (in particular H 2 O 2 ) and quantify the amount of plasma-produced H 2 O 2 species that penetrates into a gelatin hydrogel. We show that the diffusion constant of H 2 O 2 in 10% gelatin hydrogel is similar to its diffusion constant in water and that the production of H 2 O 2 in the hydrogel is significantly less than the production of H 2 O 2 in distilled water for the same plasma operation conditions suggesting that the scavenging of OH radicals at the plasma-gel interface significantly reduces the H 2 O 2 production.

60 APPLIED LIFE SCIENCES↗
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↗
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↗