Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction diffusion systems”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effects of diffusion barriers on reaction wave stability in Co/Al reactive multilayers

Bimetallic, reactive multilayers are uniformly structured materials composed of alternating sputter-deposited layers that may be ignited to produce self-propagating mixing and formation reactions. These nanolaminates are most commonly used as rapid-release heat sources. The specific chemical composition at each metal/metal interface determines the rate of mass transport in a mixing and formation reaction. The inclusion of engineered diffusion barriers at each interface will not only inhibit solid-state mixing but also may impede the self-propagating reactions by introducing instabilities to wavefront morphology. Here, this work examines the effect of adding diffusion barriers on the propagation of reaction waves in Co/Al multilayers. The Co/Al system has been shown to exhibit a reaction propagation instability that is dependent on the bilayer thickness, which allows for the occurrence of unstable modes in otherwise stable designs from the inclusion of diffusion barriers. Based on the known stability criteria in the Co/Al multilayer system, the way in which the inclusion of diffusion barriers changes a multilayer's heat of reaction, thermal conductivity, and material mixing mechanisms can be determined. These factors, in aggregate, lead to changes in the wavefront velocity and stability.

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

Mass transfer coefficients, kL, and air-CO2 ingassing rates in 3.4 m2 and 1-acre raceway ponds.

The US DOE 2016 Billion Ton Update (Langholtz 2016) projected that CO2 in flue-gas from power plants or similar sources would limit U.S. algal biofuel potential to under 5 billion gallons gasoline equivalent (gge) per year, due to lack of sufficient land, water and other requirements near such flue gas sources. 2nd generation carbon capture technologies are proposed to overcome this constraint (Davis 2017), and could expand this resource potential nearly 10-fold, assuming CO2 costs near the future, 2025 DOE NETL, goal of $40/metric ton CO2. An alternative pathway is the direct uptake of air CO2 into the algal ponds through an increase in the air-CO2 transfer rate resulting from the reaction of CO2 with hydroxide ions. However, such a ‘chemical enhancement’ in mass transfer depends on a high pH in the culture, and, by extension, the ability of microalgae to maintain high rates of carbon assimilation under such conditions. Abiotic experiments characterized the air-CO2 mass-transfer rate in 1-acre and 3.4 m2 raceway ponds as a function of pH, with rates approaching 10 g C/m2-day at pH 12 in the brackish water tested, equivalent to a biomass productivity potential of near 20 g AFDW/m2-day (@ 0.47 g C/g AFDW). The mass transfer rate was found to be nearly independent of the mass-transfer coefficient, indicating that at turbulence levels achievable in typical raceway ponds, the system is reaction rate, rather than diffusion limited. In biotic trials with an unknown green microalga, biomass productivity in ponds receiving CO2 from only direct surface air-CO2 exchange averaged 5 +/- 0.5 g AFDW/m2-day in August - September central California conditions, about half-that of experimental controls fertilized with supplemental CO2. pH in the air-only treatment reached a maximum of 10.5, supporting a model predicted air-CO2 ingassing rate between 2-3 g C/m2-day, consistent with the observed biomass carbon assimilation rate. Results suggest that unique alkaliphilic strains are required if reliant on direct in-pond air-CO2 transfer for algal inorganic carbon supply, and that reaching an economically viable biomass productivity will require strains that thrive at pH 11 and above.

09 BIOMASS FUELS↗
DOE OSTI2022

Long-time integration of parametric evolution equations with physics-informed DeepONets

Ordinary and partial differential equations (ODEs/PDEs) play a paramount role in analyzing and simulating complex dynamic processes across all corners of science and engineering. In recent years machine learning tools are aspiring to introduce new effective ways of simulating such equations, however existing approaches are not able to reliably return stable and accurate predictions across long temporal horizons. We aim to address this challenge by introducing an effective framework for learning evolution operators that map random initial conditions to associated ODE/PDE solutions within a short time interval. Such operators can be parametrized by deep neural networks that are trained in an entirely self-supervised manner without requiring one to generate any paired input-output observations. Global long-time predictions across a range of initial conditions can be then obtained by iteratively evaluating the trained model using each prediction as the initial condition for the next evaluation step. Here, this introduces a new approach to temporal domain decomposition that is shown to be effective in performing accurate long-time simulations for a wide range of parametric ODE and PDE systems, from wave propagation, to reaction-diffusion dynamics and stiff chemical kinetics, introducing a new way of rapidly emulating non-equilibrium processes in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Impact of cement composition, brine concentration, diffusion rate, reaction rate and boundary condition on self-sealing predictions for cement-CO 2 systems

Geological CO 2 storage (GCS) plays an important role in curbing CO 2 emissions by reducing the carbon footprint of difficult to decarbonize operations and achieve negative CO 2 emissions through activities like Bioenergy with Carbon Capture and Storage (BECCS) and Direct Air Carbon Capture and Storage (DACCS). Leakage of CO 2 through wells is an important concern when it comes to deployment of large-scale GCS. Existing wells in sites that are otherwise suitable for GCS can act as conduits for stored CO 2 to escape the reservoir. There is broad consensus that the main risk of leakage through wellbores is via fractures/damaged pathways. The results from several studies evaluating the permeability evolution of cement fractures in wells upon leakage of CO 2 agree that smaller fracture apertures, slower brine velocities and higher brine residence times promote self-sealing of fractures by mineral precipitation. Quantitatively, however, the differences in sealing conditions are significant and are typically attributed to differences in experimental conditions or model assumptions. Here we examine the sensitivity of our model, describing CO 2 leakage through wellbores, to cement composition, brine concentration, diffusion rates, and reaction rates. We also evaluate the impact of the boundary condition to allow comparisons between observations from experiments performed at constant flow rate and model predictions made at constant pressure conditions. Our results show that diffusion and reactions rates have the most impact on the self-sealing criteria for cement-CO 2 systems. In addition, conditions associated with self-sealing of fractures at constant flow rate require longer fractures, smaller fracture apertures, and slower velocities than under constant pressure.

58 GEOSCIENCES↗
DOE OSTI2021

Hyperbolic Reformulation Approach to Enable Efficient Simulation of Groundwater Flow and Reactive Transport

We apply Cattaneo's relaxation approach to the one-dimensional coupled Boussinesq groundwater flow and advection-diffusion-reaction equations, commonly used in engineering applications to simulate contaminant transport in the subsurface. The diffusion-type governing equations are reformulated as a hyperbolic system, augmented by an equation that can be interpreted as a momentum balance. The hyperbolization enables an efficient unified computation of the primary variable and its gradients, for example piezometric head and unit discharge in the Boussinesq equation. An augmented Roe scheme is used to solve the hyperbolic system. The hyperbolized system of equations is studied in a set of steady state and transient test cases with idealized geometry. These test cases confirm the equivalence of the hyperbolic system to its original formulation. The larger time step size of the hyperbolic equation is verified theoretically by means of a stability analysis and numerically in the test cases. Finally, a reach-scale application of flow and transport across a river meander is considered. This application case shows that the performance of the hyperbolic relaxation approach holds for more realistic groundwater flow and transport problems, relevant to water resources management.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

Excitable dynamics in high-Lewis number premixed gas combustion at normal and microgravity

Freely-propagating, premixed gas flames in high-Lewis (Le) number, quiescent mixtures are studied experimentally in tubes of various diameter at normal (lg) and microgravity (mu g). A premixture of lean butane and oxygen diluted with helium, argon, neon, nitrogen or a mixture of multiple diluents is examined such that the thermal diffusivity of the mixture (and to a lesser extent, the mass diffusivity of the rate-limiting component) is systematically varied. In effect, different diluents allow variation of the Le without changing the chemistry. The flames are recorded with high speed cinematography and their stability is visually assessed. Different modes of propagation were observed depending on the diameter of the tubes (different conductive heat loss), the composition of the mixture and the g-level. At 1g, four modes of propagation were observed in small and intermediate diameter tubes (large conductive heat loss): (1) steadily propagating flames, (2) radial and longitudinal pulsating flames, (3) 'wavering' flames, and (4) rotating spiral flames. As the diameter of the tube increases, the radial modes become more pronounced while the longitudinal modes systematically disappear. Also, multiple, simultaneous, spatially-separated 'pacemaker' sites are observed in intermediate and large diameter tubes. Each site starts as a small region of high luminosity and develops into a flamelet which assumes the form of one of the fore mentioned modes. These flamelets eventually interact, annihilate each other in their regions of intersection and merge at their newly created free-ends. For very large tubes, radially-propagating wave-trains (believed to be 'trigger waves') are observed. These are analogous to the radial pulsations observed in the smaller diameter tubes. At mu g, three modes of propagation have been observed: (1) steadily propagating flames, (2) radial and longitudinal pulsating flames, and (3) multi-armed, rotating flames. Since the pulsating mode exists at mu g and 1g, buoyant flicker is not the mechanism which drives the pulsations. Moreover, all of the instabilities at 1g and mu g have characteristic frequencies on the O(100Hz). This value is lower than the fundamental, longitudinal acoustic frequencies of the tubes which suggests that the instabilities are not acoustically driven. The patterns formed by this reaction bear remarkable similarities with the patterns formed in most excitable media when the behavior of the system is driven by couplings between chemical reaction and diffusion (e.g., Belousov-Zhabotinsky reaction, Patterns in slime molds, spiral waves in the retina of a bird's eye). While it is recognized that the chemical mechanism associated with this premixed gas reaction is exponentially sensitive to temperature and undoubtedly different from those which govern previously observed excitable media (most are isothermal, or weakly exothermic, liquid phase reactions), similar spatial and temporal patterns should not come as a complete surprise considering heat and mass diffusion are self similar. It is concluded that this premixed gas system is a definitive example of a diffusive-thermal, gas-phase oscillator based on these experimental results and their favorable comparison with theory.

Pearlman, Howard↗
DOE OSTI2022

Numerical Study of Mineral Nucleation and Growth on a Substrate

We report nucleation and growth processes of minerals and other crystals can significantly affect one another due to the transport limitations and local depletion of reactive ions in the solution. Most numerical models and experimental measurements have typically focused on either growth or nucleation, but not both. In this work, we incorporate a heterogeneous nucleation process based on classical nucleation theory into a microcontinuum model that implements the Darcy–Brinkman–Stokes approach to study the interplay between nucleation and crystal growth on a substrate in diffusive systems. We demonstrate how the Damköhler number (reaction rate) and nucleation rate prefactor change the effective nucleation rate on a substrate. Higher surface growth rates deplete the solute concentration around the nuclei that appear initially on the substrate, creating islands that screen against further nucleation. The model predicts that measured nucleation rates may be affected by the history of crystal nucleation on the substrate. In the extreme case of high growth rates relative to diffusion, it predicts that the rate of subsequent nucleation is limited by reactant depletion. We introduce a nondimensional number α to represent the relation between surface propagation rate during growth and the heterogeneous nucleation rate. We show that it is important to control Damköhler number and α to achieve similar precipitation regimes at different reaction and nucleation rates. We suggest that the observed universality can guide the interpretation of experimental results on nucleation rates, since matching experiment can be achieved by tuning transport, reaction, and nucleation parameters simultaneously. In addition, we show how the bulk solution concentration affects the structure and topology of precipitation on a substrate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2022

How the Hydrophobic Interface between a Perfluorosulfonic Acid Polymer and Water Vapor Controls Membrane Hydration

Stable hydration in perfluorinated polyelectrolyte membranes such as Nafion is essential to maintain good ion conductivity and manage permeation, especially in vapor-fed devices where water content depends on relative humidity in a gas stream. Extensive studies in the literature have shown that Nafion hydration in water vapor is controlled by its interfacial transport resistance. Nafion forms a fluorine-rich layer at the polymer-gas interface, and it has been proposed that this layer blocks water transport due to its hydrophobicity. To develop a molecular-level description of the physics underlying transport resistance in this system, we have performed a computational reaction-diffusion kinetics study of water evaporation from Nafion. Two distinct models are examined, one mimicking the blocking function proposed in the literature and the other assuming that there is no blocking, treating instead water evaporation as a dynamic balance between uptake from the gas and desorption from the polymer surface. Simulation results are compared to time-dependent infrared data over a range of 100-0% relative humidity from the literature. Only the dynamic model successfully reproduces experimental observations. This indicates that the physical nature of interfacial transport resistance is not slow diffusion across an interfacial layer; rather, it is due to the competition between dehydration and rehydration. The simulation data provide details on the accompanying water distributions throughout the membrane and on interfacial kinetics, showing that they are characterized by strong fluctuations.

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

Multi-resolution partial differential equations preserved learning framework for spatiotemporal dynamics

Traditional data-driven deep learning models often struggle with high training costs, error accumulation, and poor generalizability in complex physical processes. Physics-informed deep learning (PiDL) addresses these challenges by incorporating physical principles into the model. Most PiDL approaches regularize training by embedding governing equations into the loss function, yet this depends heavily on extensive hyperparameter tuning to weigh each loss term. To this end, we propose to leverage physics prior knowledge by “baking” the discretized governing equations into the neural network architecture via the connection between the partial differential equations (PDE) operators and network structures, resulting in a PDE-preserved neural network (PPNN). This method, embedding discretized PDEs through convolutional residual networks in a multi-resolution setting, largely improves the generalizability and long-term prediction accuracy, outperforming conventional black-box models. The effectiveness and merit of the proposed methods have been demonstrated across various spatiotemporal dynamical systems governed by spatiotemporal PDEs, including reaction-diffusion, Burgers’, and Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Reversible Reactions, Mesh Size, and Segmental Dynamics Control Penetrant Diffusion in Ethylene Vitrimers

The diffusion of two aromatic dyes with nearly identical sizes was measured in ethylene vitrimers with precise linker lengths and borate ester cross-links using fluorescence recovery after photobleaching (FRAP). One dye possessed a reactive hydroxyl group, while the second was inert. The reaction of the hydroxyl group with the network is slow relative to the hopping times of the dye, resulting in a large slowdown by a factor of 50 for a reactive probe molecule. A kinetic model was fit to the fluorescence intensity data to determine rate constants for the reversible reaction of the dye from the network, which confirms the role of slow reaction kinetics. A second network cross-linker was also investigated with a substituted boronic ester showing ∼10,000 times faster exchange kinetics. In this system, the two dyes show the same diffusion coefficient, as the reaction is no longer the rate-limiting step. The role of dense meshes on small and large dyes is also discussed in the context of the existing theories. Finally, these results highlight the potential of dynamic networks to control penetrant transport through synergistic effects of the mesh size, dynamic bond kinetics, and penetrant–network interactions.

confinement↗
DOE OSTI2023

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

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

Prediction of engine and near-field plume reacting flows in low-thrust chemical rockets

A computational model is employed to study the reacting flow within the engine and near-field plumes of several small gaseous hydrogen-oxygen thrusters. The model solves the full Navier-Stokes equations coupled with species diffusion equations for a hydrogen-oxygen reaction kinetics system and includes a two-equation q-omega model for turbulence. Predictions of global performance parameters and localized flowfield variables are compared with experimental data in order to assess the accuracy with which these flowfields are modeled and to identify aspects of the model which require improvement. Predicted axial and radial velocities 3 mm downstream of the exit plane show reasonable agreement with the measurements. The predicted peak in axial velocity in the hydrogen film coolant along the nozzle wall shows the best agreement; however, predictions within the core region are roughly 15 percent below measured values, indicating an underprediction of the extent to which the hydrogen diffuses and mixes with the core flow. There is evidence that this is due to three-dimensional mixing processes which are not included in the axisymmetric model.

Weiss, Jonathan M.↗
DOE OSTI2024

Controlled patterning of crystalline domains by frontal polymerization

Materials with hierarchical architectures that combine soft and hard material domains with coalesced interfaces possess superior properties compared with their homogeneous counterparts. These architectures in synthetic materials have been achieved through deterministic manufacturing strategies such as 3D printing, which require an a priori design and active intervention throughout the process to achieve architectures spanning multiple length scales. Here we harness frontal polymerization spin mode dynamics to autonomously fabricate patterned crystalline domains in poly(cyclooctadiene) with multiscale organization. This rapid, dissipative processing method leads to the formation of amorphous and semi-crystalline domains emerging from the internal interfaces generated between the solid polymer and the propagating cure front. The size, spacing and arrangement of the domains are controlled by the interplay between the reaction kinetics, thermochemistry and boundary conditions. Small perturbations in the fabrication conditions reproducibly lead to remarkable changes in the patterned microstructure and the resulting strength, elastic modulus and toughness of the polymer. Furthermore, this ability to control mechanical properties and performance solely through the initial conditions and the mode of front propagation represents a marked advancement in the design and manufacturing of advanced multiscale materials. Drawing inspiration from biological systems in which structural complexity develops through dissipative reaction–diffusion processes, this study explores a transformative synthetic manufacturing strategy aimed at harnessing the principles underpinning morphogenic growth, unlocking new avenues for advanced materials design and fabrication. Synthetic coupled reaction-transport processes offer a versatile yet relatively underexplored method to manipulate the spatial attributes of synthetic materials10. Here we introduce an innovative manufacturing approach based on frontal ring-opening metathesis polymerization (FROMP) that draws parallels with morphogenic growth and development, enabling the formation of patterned microstructures within polymeric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

From material properties to multiscale modeling to improve lithium-ion energy storage safety

Energy storage using lithium-ion cells dominates consumer electronics and is rapidly becoming predominant in electric vehicles and grid-scale energy storage, but the high energy densities attained lead to the potential for release of this stored chemical energy. This article introduces some of the paths by which this energy might be unintentionally released, relating cell material properties to the physical processes associated with this potential release. The selected paths focus on the anode–electrolyte and cathode–electrolyte interactions that are of typical concern for current and near-future systems. Relevant material processes include bulk phase transformations, bulk diffusion, surface reactions, transport limitations across insulating passivation layers, and the potential for more complex material structures to enhance safety. We also discuss the development, parameterization, and application of predictive models for this energy release and give examples of the application of these models to gain further insight into the development of safer energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2020

Simulating policy interventions in the interfirm diffusion of low-carbon technologies: An agent-based evolutionary game model

Policy interventions are frequently used by authorities around the world to mitigate carbon emissions. However, economic systems are essentially complex adaptive systems, which often exhibit unexpected responses to exogenous interventions and leave intervenors in a dilemma, even if the interventions are elaborately designed. To explore the behaviors of the diffusion system regarding low-carbon technologies, this study builds an agent-based model (ABM) to simulate enterprises' reactions to multiple policy interventions aimed at spurring low-carbon technology diffusion. The simulated enterprises are in a complex network where they play evolutionary games with their neighbors, which enables the model to possess two critical features of economic systems: adaptiveness and equilibrium. The model reveals a dilemma of policy interventions: intuitively, carbon taxes, asymmetric penalties (only imposed on un-low-carbon enterprises) and subsidies can improve the diffusion, which is consistent with previous research; counter-intuitively, all these policies turn out to be inefficient or even harmful to low-carbon enterprises because of the diffusion system's high adaptiveness. Specifically, when carbon taxes and penalties increase, both the low-carbon and un-low-carbon enterprises end up with the equilibrium of equivalent but lower profits. In contrast, all the enterprises earn equivalent but more profits even though subsidies are only given to the low-carbon enterprises, which implies that the un-low-carbon enterprises (indirectly) grab a portion of the subsidies and accordingly weaken the incentive effect of subsidies. These system behaviors are summarized as "equalizing effect" because the system tends to equalize the impact of both positive and negative interventions among all enterprises even though the policy interventions are asymmetrically imposed on one type of the enterprises. The findings also indicate that policies implemented to enlarge green market sizes can help policymakers bypass the dilemma.

Agent-based modeling↗
DOE OSTI2022

A Phase-Field Study on Internal to External Oxidation Transition in High-Temperature Structural Alloys

Structural alloys applied at high temperatures rely on an external dense layer of oxide scale for protection. As some reactive alloy components are selectively oxidized internally, understanding how these dispersed metal oxide particles reach the surface to form a protective scale—the so-called internal to external oxidation transition—is crucial for designing these alloys. While the literature is replete with experimental studies on oxidation of alloys, there is a lack of computational studies in this realm due to the complex nature of coupled reaction and diffusion processes in multicomponent multi-phase alloy systems. In this work, we apply a recently developed phase-field model to simulate the oxidation processes under different compositions and nucleation scenarios to gain insights into how a continuous oxide scale can be established. The results show that while alloy composition is critical for internal to external oxidation transition, the oxide nuclei size, shape and distribution also have significant impact on the transition kinetics.

36 MATERIALS SCIENCE↗