Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Benchmarking magnetised three-wave coupling for laser backscattering: analytic solutions and kinetic simulations

Understanding magnetised laser–plasma interactions is important for controlling magneto-inertial fusion experiments and developing magnetically assisted radiation and particle sources. For nanosecond pulses at non-relativistic intensities, interactions are dominated by coherent three-wave interactions, whose nonlinear coupling coefficients became known only recently when waves propagate at oblique angles with the magnetic field. In this paper, backscattering coupling coefficients predicted by warm-fluid theory are benchmarked using particle-in-cell simulations in one spatial dimension, and excellent agreements are found for a wide range of plasma temperatures, magnetic field strengths and laser propagation angles, when the interactions are mediated by electron-dominant hybrid waves. Systematic comparisons between theory and simulations are made possible by a rigorous protocol. On the theory side, the initial boundary value problem of linearised three-wave equations is solved, and the transient-time solutions allow the effects of growth and damping to be distinguished. On the simulation side, parameters are carefully chosen and calibration runs are performed to ensure that comparisons are well controlled. Fitting simulation data to analytical solutions yields numerical growth rates that match theory predictions within error bars. Although warm-fluid theory is found to be valid for a wide parameter range, genuine kinetic effects have also been observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

In situ observation of the magnetization configuration and reversal in cylindrical nanowires

Curvilinear magnetic structures often have unique magnetic behavior compared to their rectilinear counterparts. This is due to the unique curvilinear boundary conditions as well as the curvature induced Dzyaloshinskii–Moriya-like interaction and the curvature induced anisotropy. The effects of a curvilinear geometry are best studied in 3D structures, where the curvature can have a significant spatial extent. Of these 3D structures, the simplest structure to study is the cylindrical nanowire. Here, we have simulated the magnetization reversal in cylindrical NiFe nanowires and present in situ Lorentz TEM images to support the findings of the simulations. We studied the domain formation and reversal of nanowires with two distinct diameters that give rise to a different reversal behavior. We have, thus, found that the zero-field magnetization configuration in these wires can take on a double helix type of configuration. The reversal in these structures then proceeds through the winding and unwinding of these helical configurations rather than through domain wall propagation.

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

Diamond of infrared equivalences in abelian gauge theories

We demonstrate a tree-level equivalence between four distinct infrared objects in (d+2)-dimensional Abelian gauge theories. These are (i) the large gauge charge Q ϵ where the function ϵ on the sphere parametrizing large gauge transformations is identified with the Goldstone mode θ of spontaneously broken large gauge symmetry; (ii) the soft effective action that captures the dynamics of the soft and Goldstone modes; (iii) the edge mode action with Neumann boundary conditions; and (iv) the Wilson line dressing of a scattering amplitude, including a novel dressing for soft photons, which have local charge distributions despite having vanishing global charge. The promotion of the large gauge parameter to the dynamical Goldstone and the novel dressing of soft gauge particles give rise to intriguing possibilities for the future study of infrared dynamics of gauge theories and gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Accelerated Creep Test Qualification of Creep-Resistance Using the Wilshire–Cano–Stewart Constitutive Model and Stepped Isostress Method

In this paper, a qualification of accelerated creep-resistance of Inconel 718 is assessed using the novel Wilshire–Cano–Stewart (WCS) model and the stepped isostress method (SSM) and predictions are made to conventional creep data. Conventional creep testing is a long-term continuous process; in fact, the ASME B&PV III requires that 10,000+ h of experiments must be conducted to each heat for materials employed in boilers and/or pressure vessel components. This process is costly and not feasible for rapid development of new materials. As an alternative, accelerated creep testing techniques have been developed to reduce the time needed to characterize the creep resistance of materials. Most techniques are based upon the time-temperature-stress superposition principle that predicts minimum-creep-strain-rate (MCSR) and stress-rupture behaviors but lack the ability to predict creep deformation and consider deformation mechanisms that occur for experiments of longer duration. The SSM has been developed, which enables the prediction of creep deformation response as well as reduce the time needed for qualification of materials. The SSM approach has been successful for polymer, polymeric composites, and recently has been introduced for metals. In this study, the WCS constitutive model, calibrated to SSM test data, qualifies the creep resistance of Inconel 718 at 750 °C and predictions are compared to conventional creep testing data. The WCS model has proven to make long-term predictions for stress-rupture, MCSR, creep deformation, and damage in metallic materials. The SSM varies stress levels after time interval adding damage to the material, which can be tracked by the WCS model. The SSM data is calibrated into the model and the WCS model generates realistic predictions of stress-rupture, MSCR, damage, and creep deformation. The calibrated material constants are used to generate predictions of stress-rupture and are postaudit validated using the National Institute of Material Science database. Similarly, the MCSR predictions are compared from previous studies. Finally, the creep deformation predictions are compared with real data and is determined that the results are well in between the expected boundaries. Material characterization and mechanical properties can be determined at a faster rate and with a more cost-effective method. This is beneficial for multiple applications such as in additive manufacturing, composites, spacecraft, and industrial gas turbines.

42 ENGINEERING↗
DOE OSTI2023

Osmotic Pressure Gradient Effects on Water Diffusion in Porous Rock: Can They Pervert Permeability Tests?

Abstract Generation of a large network of hydraulic cracks is of key importance not only for the success of fracking of shale but also for the recent scheme of sequestration of CO2 in deep formations of basalt and peridotite, which are mafic and ultramafic rocks that combine chemically with CO2. In numerical simulation of the creation of a fracture network in porous rock, an important goal is to enhance the rock permeability. The objective of this article is to calculate the effect of osmotic pressure gradients caused by gradients of concentration of the ions of Ca, Mg, Na, etc. on the effective permeability of the rock. The basic differential equations are formulated, and their explicit solutions for appropriate initial and boundary conditions are obtained under certain plausible simplifications. The main result is explicit approximate formulas for the critical time before which no water permeation through a test specimen can be observed. Depending on various parameters, this time can be unacceptably long, which is manifested as a zero water outflow. The solution may also explain the unreasonably small permeability values reported for some shales.

Mechanics↗
DOE OSTI2022

Indications of a Decrease in the Depth of Deep Convective Cores with Increasing Aerosol Concentration during the CACTI Campaign

Abstract An aerosol indirect effect on deep convective cores (DCCs), by which increasing aerosol concentration increases cloud-top height via enhanced latent heating and updraft velocity, has been proposed in many studies. However, the magnitude of this effect remains uncertain due to aerosol measurement limitations, modulation of the effect by meteorological conditions, and difficulties untangling meteorological and aerosol effects on DCCs. The Cloud, Aerosol, and Complex Terrain Interactions (CACTI) campaign in 2018–19 produced concentrated aerosol and cloud observations in a location with frequent DCCs, providing an opportunity to examine the proposed aerosol indirect effect on DCC depth in a rigorous and robust manner. For periods throughout the campaign with well-mixed boundary layers, we analyze relationships that exist between aerosol variables (condensation nuclei concentration > 10 nm, 0.4% cloud condensation nuclei concentration, 55–1000-nm aerosol concentration, and aerosol optical depth) and meteorological variables [level of neutral buoyancy (LNB), convective available potential energy, midlevel relative humidity, and deep-layer vertical wind shear] with the maximum radar-echo-top height and cloud-top temperature (CTT) of DCCs. Meteorological variables such as LNB and deep-layer shear are strongly correlated with DCC depth. LNB is also highly correlated with three of the aerosol variables. After accounting for meteorological correlations, increasing values of the aerosol variables [with the exception of one formulation of aerosol optical depth (AOD)] are generally correlated at a statistically significant level with a warmer CTT of DCCs. Therefore, for the study region and period considered, increasing aerosol concentration is mostly associated with a decrease in DCC depth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

ACCELERATED CREEP TEST (ACT) QUALIFICATION OF CREEP RESISTANCE USING THE WCS CONSTITUTIVE MODEL AND STEPPED ISOSTRESS METHOD (SSM)

In this study, a qualification of accelerated creep-resistance of Inconel 718 is assessed using the novel Wilshire-Cano-Stewart (WCS) model and the stepped isostress method (SSM) and predictions are made to conventional creep data. Conventional creep testing (CCT) is a long-term continuous process, in fact, the ASME B&PV III requires that 10,000+ hours of experiments must be conducted to each heat for materials employed in boilers and/or pressure vessel components. This process is costly and not feasible for rapid development of new materials. As an alternative, accelerated creep testing techniques have been developed to reduce the time needed to characterize the creep resistance of materials. Most techniques are based upon the time-temperature-stress superposition principle (TTSSP) that predicts minimum-creep-strain-rate (MCSR) and stress-rupture behaviors but lack the ability to predict creep deformation and consider deformation mechanisms that occur for experiments of longer duration. The stepped isostress method (SSM) has been developed which enables the prediction of creep deformation response as well as reduce the time needed for qualification of materials. The SSM approach has been successful for polymer, polymeric composites, and recently has been introduced for metals. In this study, the WCS constitutive model, calibrated to SSM test data, qualifies the creep resistance of Inconel 718 at 750°C and predictions are compared to CCT data. The WCS model has proven to make long-term predictions for stress-rupture, minimum-creep-strain-rate (MCSR), creep deformation, and damage in metallic materials. The SSM varies stress levels after time interval adding damage to the material, which can be tracked by the WCS model. The SSM data is calibrated into the model and the WCS model generates realistic predictions of stress-rupture, MSCR, damage, and creep deformation. The calibrated material constants are used to generate predictions of stress-rupture and are post-audit validated using the National Institute of Material Science (NIMS) database. Similarly, the MCSR predictions are compared from previous studies. Finally the creep deformation predictions are compared with real data and is determined that the results are well in between the expected boundaries. Material characterization and mechanical properties can be determined at a faster rate and with a more cost-effective method. This is beneficial for multiple applications such as in additive manufacturing, composites, spacecraft, and Industrial Gas Turbines (IGT).

36 MATERIALS SCIENCE↗
DOE OSTI2024

An effective cosmological collider

Effective field theories (EFTs) of heavy particles coupled to the inflaton are rife with operator redundancies, frequently obscured by sensitivity to both boundary terms and field redefinitions. We initiate a systematic study of these redundancies by establishing a minimal operator basis for an archetypal example, the abelian gauge-Higgs-inflaton EFT. Working up to dimension 9, we show that certain low-dimensional operators are entirely redundant and identify new non-redundant operators with potentially interesting cosmological collider signals. Our methods generalize straightforwardly to other EFTs of heavy particles coupled to the inflaton.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Corrosion Resistance of 309L Stainless Steel Claddings on Carbon Steel Produced with Wire-Fed Directed Energy Deposition

Additive manufacturing (AM) tools are capable of applying overlay austenitic stainless steel (SS) claddings to carbon steel components. The benefits of this approach over arc welding include a smaller heat-affected zone, residual stress reduction, and material savings. In particular, wire-directed energy deposition is a suitable technique because of its low material cost and high rate of production compared to other AM methods. However, metallurgical variations in composition, phase fraction, and microsegregation can potentially influence the corrosion behavior of such claddings. In this work, 309L SS is clad on carbon steel substrates and electrochemical methods are used to measure their general and pitting corrosion resistance in simulated marine environments (3.5 wt% NaCl solutions). Two-layer claddings are fabricated with four laser powers to understand the effects of bulk chemical composition, austenite/δ-ferrite phase fractions, and individual phase compositions on corrosion behavior. The two-layer claddings are compared to a single-layer cladding, wrought 304 SS, and the carbon steel substrate for a comprehensive assessment of corrosion performance. The two-layer claddings are remarkably resistant to general corrosion in the 3.5 wt% NaCl environment because of their high Cr content (21.6 wt% to 23.3 wt% Cr). The single-layer cladding exhibits localized corrosion at unmixed Fe-rich peninsulas that originate at the dissimilar metal boundary and protrude into the first cladding layer. All two-layer claddings possess higher pitting corrosion resistance than wrought 304 SS, demonstrating their effectiveness as a corrosion-resistant barrier. Here, the pitting corrosion resistance is superior for claddings made with lower laser powers, due to low dilution and greater δ-ferrite contents.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Partnership Center for High-fidelity Boundary Plasma Simulation

The University of Colorado Boulder contributed to the DOE Partnership Center for High‑fidelity Boundary Plasma Simulation (award DE‑SC0018271) by advancing gyrokinetic modeling of core–edge plasma interactions in tokamaks. Using the GEM and XGC codes, the project investigated turbulence spreading from the plasma edge to the core, finding that while edge‑driven trapped electron mode turbulence enhances core ion heat flux in DIII‑D L‑mode plasmas, it does not fully explain experimentally observed transport shortfalls, indicating a need for flux‑driven models. The team also demonstrated that edge micro‑instabilities are highly sensitive to small variations in magnetic equilibrium, particularly near the pedestal and edge, highlighting the importance of accurately characterizing equilibrium uncertainties. Significant progress was made in spatially coupling the core GEM and edge XGC simulations, with successful validation against standalone XGC results in realistic tokamak geometries, laying the groundwork for future extensions to kinetic electrons and electromagnetic effects.

42 ENGINEERING↗
DOE OSTI2020

Determination of Elastic Coefficients for a Composite with Crossing Steel - 20086

For a composite material with an embedding of crossing reinforcement bars, the elastic coefficients are calculated by solving a boundary-value problem (elastostatic problem). The background material is concrete whereas the bars are made of steel. The bars cross each other horizontally at a right angle. The boundary-value problem is solved numerically by using 3-D quadratic basis functions in ABAQUS{sup R}. Under external shear loading conditions, the effective elastic coefficients on the macroscale are quite larger for a cross-section in which the bars are immediately below the ones in the cross-section normal to it. The shear modulus in the off-plane direction is about one and half times larger than that in the in-plane direction. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Tensile-stress effect on ferroelectric Barkhausen noise

Here, this study examines the effect of tensile stress on the ferroelectric properties of Pb(Zr 0.4 Ti 0.6 )O 3 thin film, with a focus on Barkhausen noise, observed for the first time under such conditions. Tensile stress significantly alters domain wall motions, affecting Barkhausen noise more than average polarization. Frequency analysis identifies grain boundaries as primary pinning sites, consistent across stress levels. A nonlinear relationship between stress, domain wall mobility, and polarization is found, where increased stress initially enhances pinning and polarization changes, but this effect diminishes at higher stress levels, indicating a shift in behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cloud adjustments from large-scale smoke–circulation interactions strongly modulate the southeastern Atlantic stratocumulus-to-cumulus transition

Smoke from southern Africa blankets the southeastern Atlantic Ocean from June to October, producing strong and competing aerosol radiative effects. Smoke effects on the transition between overcast stratocumulus and scattered cumulus clouds are investigated along a Lagrangian (air-mass-following) trajectory in regional climate and large eddy simulation models. Results are compared with observations from three recent field campaigns that took place in August 2017: ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES), CLouds and Aerosol Radiative Impacts and Forcing: Year 2017 (CLARIFY), and Layered Atlantic Smoke Interactions with Clouds (LASIC). The case study is set up around the joint ORACLES–CLARIFY flight that took place near Ascension Island on 18 August 2017. Smoke sampled upstream on an ORACLES flight on 15 August 2017 likely entrained into the marine boundary layer later sampled during the joint flight. The case is first simulated with the WRF-CAM5 regional climate model in three distinct setups: (1) FireOn, in which smoke emissions and any resulting smoke–cloud–radiation interactions are included; (2) FireOff, in which no smoke emissions are included; (3) RadOff, in which smoke emissions and their microphysical effects are included but aerosol does not interact directly with radiation. Over the course of the Lagrangian trajectory, differences in free tropospheric thermodynamic properties between FireOn and FireOff are nearly identical to those between FireOn and RadOff, showing that aerosol–radiation interactions are primarily responsible for the free tropospheric effects. These effects are non-intuitive: in addition to the expected heating within the core of the smoke plume, there is also a “banding” effect of cooler temperature (~1–2 K) and greatly enhanced moisture (>2 g kg -1 ) at the plume top. This banding effect is caused by a vertical displacement of the former continental boundary layer in the free troposphere in the FireOn simulation resulting from anomalous diabatic heating due to smoke absorption of sunlight that manifests primarily as a few hundred meters per day reduction in large-scale subsidence over the ocean. A large eddy simulation (LES) is then forced with free tropospheric fields taken from the outputs for the WRF-CAM5 FireOn and FireOff runs. Cases are run by selectively perturbing one variable (e.g., aerosol number concentration, temperature, moisture, vertical velocity) at a time to better understand the contributions from different indirect (microphysical), “large-scale” semi-direct (above-cloud thermodynamic and subsidence changes), and “local” semi-direct (below-cloud smoke absorption) effects. Despite a more than 5-fold increase in cloud droplet number concentration when including smoke aerosol concentrations, minimal differences in cloud fraction evolution are simulated by the LES when comparing the base case with a perturbed aerosol case with identical thermodynamic and dynamic forcings. A factor of 2 decrease in background free tropospheric aerosol concentrations from the FireOff simulation shifts the cloud evolution from a classical entrainment-driven “deepening–warming” transition to trade cumulus to a precipitation-driven “drizzle-depletion” transition to open cells, however. The thermodynamic and dynamic changes caused by the WRF-simulated large-scale adjustments to smoke diabatic heating strongly influence cloud evolution in terms of both the rate of deepening (especially for changes in the inversion temperature jump and in subsidence) and in cloud fraction on the final day of the simulation (especially for the moisture “banding” effect). Such large-scale semi-direct effects would not have been possible to simulate using a small-domain LES model alone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗
DOE OSTI2023

Effects of scattering on the field-induced T c enhancement in thin superconducting films in a parallel magnetic field

The problem of the normal-superconducting phase boundary for films in a parallel magnetic field, discussed in the classical paper by Ginzburg and Landau for temperatures close to the critical, is revisited with the help of the microscopic BCS theory for arbitrary temperatures taking pair-breaking and transport scattering into account. Although confirming experimental findings of the T c enhancement by the magnetic field, we find that the transport scattering pushes the phase transition curve to higher fields and higher temperatures for nearly all practical scattering rates. Still, the T c enhancement disappears in the dirty limit. Here, we also consider intriguing changes, such as reentrant superconductivity, caused to the phase boundary by pair-breaking magnetic ions spread on one of the film faces. These features await experimental verification.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Quantifying leakage fields at ionic grain boundaries using off-axis electron holography

The electrical properties of interfaces in semiconductors and ionic conductors are immensely important in a wide range of applications. Electron holography is ideally suited for the direct measurement of the electrostatic potential of such interfaces. A key challenge with this approach is the contribution of the leakage field from the sample to the observed electron phase shift. This leakage field cannot be a priori independently determined and can cause an overestimation of the phase shift. In this work, we use finite element simulations to compute the three-dimensional electrostatic potential in the vicinity of an interface associated with a given interfacial charge density distribution. We then evaluate the predicted phase shift and demonstrate that the leakage field strongly affects the recovery of the projected interface potential. From the difference between the true potential and uncorrected, recovered potential, we propose a method to correct for this effect. We then demonstrate the application of this methodology to the analysis of experimental off-axis electron holography data acquired from the grain boundaries in lightly doped ceria.

36 MATERIALS SCIENCE↗