Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface diffusion”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling a ring magnet in ALEGRA

We show here that Sandia's ALEGRA software can be used to model a permanent magnet in 2D and 3D, with accuracy matching that of the open-source commercial software FEMM. This is done by conducting simulations and experimental measurements for a commercial-grade N42 neodymium alloy ring magnet with a measured magnetic field strength of approximately 0.4 T in its immediate vicinity. Transient simulations using ALEGRA and static simulations using FEMM are conducted. Comparisons are made between simulations and measurements, and amongst the simulations, for sample locations in the steady-state magnetic field. The comparisons show that all models capture the data to within 7%. The FEMM and ALEGRA results agree to within approximately 2%. The most accurate solutions in ALEGRA are obtained using quadrilateral or hexahedral elements. In the case where iron shielding disks are included in the magnetized space, ALEGRA simulations are considerably more expensive because of the increased magnetic diffusion time, but FEMM and ALEGRA results are still in agreement. The magnetic field data are portable to other software interfaces using the Exodus file format.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Elucidating the Link Between Alkali Metal Ions and Reaction-Transport Mechanisms in Cathode Electrodes for Alkali-ion Batteries

Our long-term goal is to improve the reliability of electrode materials and their ability to transport and store various metal ions for electrochemical energy storage applications. The main objective of this work was to investigate the intrinsic relationship between the role of alkali metal ions and electrochemically driven mechanical stability and kinetic properties of battery materials. The overall question was “What is the role of alkali metal ions on the electrochemical and mechanical behavior of cathode electrodes? Our guiding hypothesis was that intercalation of larger alkali metal ions (Na and K) inevitably alters the coupled transport-reaction processes during battery operation in organic electrolytes, leading to more intensive chemo-mechanical instabilities in cathode electrodes, resulting in rapid capacity fade. To validate the hypothesis, we experimentally characterized the reaction-transport processes and governing forces driving the instability of electrode materials in different alkali metal-ion environments. The project had three main tasks. The first one was to investigate intercalation-induced strains and associated stress generation, and their impact on structural deformations in composite cathode electrodes. The second task focused on identifying potential-dependent dynamic changes in the electrode-electrolyte interface in alkali metal ion batteries. The last task was focused on determining how larger alkali metal ions with slower diffusivity affect the transport-mechanics coupling at faster scan rates, compared to smaller ions with faster diffusivity in electrodes. We shortly provided the outcome of each task in the accomplishment section. This project produced 10 peer-reviewed publications (9 research papers and one review manuscript) and supported two Ph.D. students, who graduated from Oklahoma State University.

25 ENERGY STORAGE↗
DOE OSTI2020

Investigating the Linear Thermal Expansion of Additively Manufactured Multi-Material Joining between Invar and Steel

This work investigated the linear thermal expansion properties of a multi-material specimen fabricated with Invar M93 and A36 steel. A sequence of tests was performed to investigate the viability of additively manufactured Invar M93 for lowering the coefficient of thermal expansion (CTE) in multi-material part tooling. Invar beads were additively manufactured on a steel base plate using a fiber laser system, and samples were taken from the steel, Invar, and the interface between the two materials. The CTE of the samples was measured between 40 °C and 150 °C using a thermomechanical analyzer, and the elemental composition was studied with energy dispersive X-ray spectroscopy. The CTE of samples taken from the steel and the interface remained comparable to that of A36 steel; however, deviations between the thermal expansion values were prevalent due to element diffusion in and around the heat-affected zone. The CTEs measured from the Invar bead were lower than those from the other sections with the largest and smallest thermal expansion values being 10.40 μm/m-K and 2.09 μm/m-K. In each of the sections, the largest CTE was measured from samples taken from the end of the weld beads. An additional test was performed to measure the aggregate expansion of multi-material tools. Invar beads were welded on an A36 steel plate. The invar was machined, and the sample was heated in an oven from 40 °C and 160 °C. Strain gauges were placed on the surface of the part and were used to analyze how the combined thermal expansions of the invar and steel would affect the thermal expansion on the surface of a tool. There were small deviations between the expansion values measured by gauges placed in different orientations, and the elongation of the sample was greatest along the dimension containing a larger percentage of steel. On average, the expansion of the machined Invar surface was 42% less than the expansion of the steel surface. The results of this work demonstrate that additively manufactured Invar can be utilized to decrease the CTE for multi-material part tooling.

36 MATERIALS SCIENCE↗
NASA NTRS1981

Assessment of high temperature metallizations for I(2)L and CMOS technologies

High temperature barrier metallizations were assessed and tested for I(2)L and CMOS applications. Life tests were accelerated to 375 C in view of the -55 C to +300 C temperature range established for engine-located electronics without fuel cooling. The gold-refractory metallizations evaluated were Au-TiW-PtSi, Au-TiW/TiO2/TiW-PtSi and Au-TiW(N)-PtSi. These metallization systems were thermally annealed to at least 375 C for up to 250 hours. The critical requirement for stable diffusion barrier is the TiW grain size. Small grain (250A-500A) films were observed to be stable up to 375 C. Deposition to TiW diffusion barrier in the presence of oxygen and nitrogen also results in an effective diffusion barrier. Auger electron spectroscopy profiles of the PtSi indicates some penetration by the TiW. In the case of PtSi/TiW interface, the redistribution of oxygen further passivates the system by forming a TiO2 layer at the interface.

Christou, A.↗
NASA NTRS1995

Single-drop reactive extraction/extractive reaction with forced convective diffusion and interphase mass transfer

An algorithm has been developed for the forced convective diffusion-reaction problem for convection inside and outside a droplet by a recirculating flow field hydrodynamically coupled at the droplet interface with an external flow field that at infinity becomes a uniform streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet or reactions can take place in both phases. The algorithm has been implemented and results are shown here for the case of no reaction and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Kleinman, Leonid S.↗
DOE OSTI2021

Turbulent mixing by microscopic self-assembled spinners

A system for mixing particles that includes a liquid comprising inert particles and defining a liquid and air interface; magnetic microparticles suspended at the liquid and air interface; and a magnetic source configured to apply a uniaxial alternating magnetic field parallel to the liquid and air interface, wherein the uniaxial alternating magnetic field promotes a turbulent motion of the magnetic microparticles, which in turn promotes a diffusive motion of the inert particles.

Snezhko, Oleksiy↗
NASA NTRS2011

Test-Anchored Vibration Response Predictions for an Acoustically Energized Curved Orthogrid Panel with Mounted Components

A rich body of vibroacoustic test data was recently generated at Marshall Space Flight Center for a curved orthogrid panel typical of launch vehicle skin structures. Several test article configurations were produced by adding component equipment of differing weights to the flight-like vehicle panel. The test data were used to anchor computational predictions of a variety of spatially distributed responses including acceleration, strain and component interface force. Transfer functions relating the responses to the input pressure field were generated from finite element based modal solutions and test-derived damping estimates. A diffuse acoustic field model was employed to describe the assumed correlation of phased input sound pressures across the energized panel. This application demonstrates the ability to quickly and accurately predict a variety of responses to acoustically energized skin panels with mounted components. Favorable comparisons between the measured and predicted responses were established. The validated models were used to examine vibration response sensitivities to relevant modeling parameters such as pressure patch density, mesh density, weight of the mounted component and model form. Convergence metrics include spectral densities and cumulative root-mean squared (RMS) functions for acceleration, velocity, displacement, strain and interface force. Minimum frequencies for response convergence were established as well as recommendations for modeling techniques, particularly in the early stages of a component design when accurate structural vibration requirements are needed relatively quickly. The results were compared with long-established guidelines for modeling accuracy of component-loaded panels. A theoretical basis for the Response/Pressure Transfer Function (RPTF) approach provides insight into trends observed in the response predictions and confirmed in the test data. The software modules developed for the RPTF method can be easily adapted for quick replacement of the diffuse acoustic field with other pressure field models; for example a turbulent boundary layer (TBL) model suitable for vehicle ascent. Wind tunnel tests have been proposed to anchor the predictions and provide new insight into modeling approaches for this type of environment. Finally, component vibration environments for design were developed from the measured and predicted responses and compared with those derived from traditional techniques such as Barrett scaling methods for unloaded and component-loaded panels.

Frady, Gregory P.↗
DOE OSTI2025

Dissolved CO 2 Modulates the Electrochemical Capacitance on Gold Electrodes

The presence of CO 2 at an electrified interface between an aqueous electrolyte and a metal electrode is the prerequisite for many electrochemical CO 2 capture technologies. To understand the behavior of dissolved CO 2 at an aqueous electrified interface, we characterized the electrochemical interface of planar gold electrodes with cyclic voltammetry, electrochemical impedance spectroscopy (EIS), electrochemical surface plasmon resonance (EC–SPR), and attenuated total reflectance surface-enhanced infrared spectroscopy (ATR–SEIRAS). Under all investigated conditions, we observed a decrease in the electrochemical capacitance upon saturation of the electrolyte with CO 2 , as compared to an electrolyte saturated with Ar. EIS and EC–SPR showed that this capacitance reduction was also potential dependent: it reached a minimum near the point of zero charge and became more significant as the applied potential moved further away from the point of zero charge. Hybrid quantum–classical simulations of the gold/aqueous electrolyte interface indicate that bicarbonate decreases the capacitance and modifies the composition of the electric double layer. In addition to the binding of bicarbonate under positive bias, we propose that molecular CO 2 can be induced by applied potential to concentrate in the diffuse layer of the electric double layer, leading to a reduction in the electrochemical capacitance under both negative and positive bias. Furthermore, this work advances the understanding of non-Faradaic effects of dissolved CO 2 at aqueous electrified interfaces of relevance for electrochemical CO 2 capture.

25 ENERGY STORAGE↗
DOE OSTI2021

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Diffusion Driven Combustion Waves in Porous Media

Filtration of gas containing oxidizer, to the reaction zone in a porous medium, due, e.g., to a buoyancy force or to an external pressure gradient, leads to the propagation of Filtration combustion (FC) waves. The exothermic reaction occurs between the fuel component of the solid matrix and the oxidizer. In this paper, we analyze the ability of a reaction wave to propagate in a porous medium without the aid of filtration. We find that one possible mechanism of propagation is that the wave is driven by diffusion of oxidizer from the environment. The solution of the combustion problem describing diffusion driven waves is similar to the solution of the Stefan problem describing the propagation of phase transition waves, in that the temperature on the interface between the burned and unburned regions is constant, the combustion wave is described by a similarity solution which is a function of the similarity variable x/square root of(t) and the wave velocity decays as 1/square root of(t). The difference between the two problems is that in the combustion problem the temperature is not prescribed, but rather, is determined as part of the solution. We will show that the length of samples in which such self-sustained combustion waves can occur, must exceed a critical value which strongly depends on the combustion temperature T(sub b). Smaller values of T(sub b) require longer sample lengths for diffusion driven combustion waves to exist. Because of their relatively small velocity, diffusion driven waves are considered to be relevant for the case of low heat losses, which occur for large diameter samples or in microgravity conditions, Another possible mechanism of porous medium combustion describes waves which propagate by consuming the oxidizer initially stored in the pores of the sample. This occurs for abnormally high pressure and gas density. In this case, uniformly propagating planar waves, which are kinetically controlled, can propagate, Diffusion of oxidizer decreases the wave velocity. In addition to the reaction and diffusion layers, the uniformly propagating wave structure includes a layer with a pressure gradient, where the gas motion is induced by the production or consumption of the gas in the reaction as well as by thermal expansion of the gas. The width of this zone determines the scale of the combustion wave in the porous medium.

Aldushin, A. P.↗
DOE OSTI2025

Microstructure characteristics of LPBF&HIP fabricated graded composite transition joint between ferritic steel and austenitic stainless steel

Graded composite transition joints (GCTJs) offer a promising alternative to conventional dissimilar metal welds (DMWs) by enabling smooth compositional and microstructural transitions. However, GCTJs fabricated solely through additive manufacturing (AM) face challenges such as heat accumulation, complex parameter control, and elemental segregation. In this study, we propose a novel approach that relies on AM to design a spatially graded structure in one alloy and then employs hot isostatic pressing (HIP) as a diffusion bonding method to join it with a second alloy. Here, this method combines the flexibility of AM with the powder net-shaping advantage of HIP. Specifically, a series of closely packed austenitic stainless steel 304 conical structures were printed using laser powder bed fusion (LPBF) and then combined with ferritic steel P91 powder via HIP. By using electron backscatter diffraction (EBSD), electron probe microanalysis (EPMA), and transmission electron microscopy (TEM) techniques, the microstructure characteristics of the GCTJ of 304&P91, especially the interdiffusion zone (IDZ), have been systematically investigated. The microstructure at the interface transitions from austenite-ferrite (A+F) to austenite-martensite-ferrite (A+M+F), and finally to martensite-ferrite (M+F) due to diffusion. Additionally, the diffusion width between 304 and P91 increases with the volume fraction of P91. This unique design also ensures a gradual transition in both hardness and thermal expansion coefficient from 304 to P91, thereby enabling a smooth gradient in functional properties. Overall, this study proposes a novel approach for fabricating GCTJs and contributes to advancing design concepts in the field of dissimilar metal joining.

Additively manufacturing (AM)↗
NASA NTRS1989

Electrochemical cell for obtaining oxygen from carbon dioxide atmospheres

For manned missions to Mars to become a reality, an efficient and reliable means of obtaining oxygen from the carbon dioxide-rich atmosphere will be required. Otherwise, the high cost of transporting the oxygen needed to sustain the astronauts will severely restrict the expedition to the martian surface. Recently, the use of electrochemical devices has been explored as a means of obtaining oxygen from the carbon dioxide-rich atmosphere. In these devices, oxygen ions diffuse through solid oxide membranes, thus, separating oxygen from the other gases presented. This phenomenon has only recently been explored as a means of obtaining large quantities of oxygen from toxic atmospheres, although first observed by Walter nernst in 1899. Nernst observed that stabilized zirconia will conduct oxygen ions when an electrical potential is applied across metallic electrodes applied to the ceramic membrane. Diatomic oxygen molecules are dissociated at the positive electrode/electrolyte interface. The oxygen ions enter the ceramic body due to the ion density gradient which is produced by the electrical potential across the electrolytic membrane. Once the ions have diffused through the membrane, they reform diatomic oxygen molecules at the anode. The separation of oxygen from carbon dioxide is achieved by the combination of thermal and electrochemical processes. The thermal decomposition of carbon dioxide (at 1000 C) results in the production of carbon monoxide and oxygen by the reaction.

Hooker, M. W.↗
NASA NTRS1992

GRO/EGRET data analysis software: An integrated system of custom and commercial software using standard interfaces

The Energetic Gamma Ray Telescope Experiment (EGRET) on the Compton Gamma Ray Observatory has been in orbit for more than a year and is being used to map the full sky for gamma rays in a wide energy range from 30 to 20,000 MeV. Already these measurements have resulted in a wide range of exciting new information on quasars, pulsars, galactic sources, and diffuse gamma ray emission. The central part of the analysis is done with sky maps that typically cover an 80 x 80 degree section of the sky for an exposure time of several days. Specific software developed for this program generates the counts, exposure, and intensity maps. The analysis is done on a network of UNIX based workstations and takes full advantage of a custom-built user interface called X-dialog. The maps that are generated are stored in the FITS format for a collection of energies. These, along with similar diffuse emission background maps generated from a model calculation, serve as input to a maximum likelihood program that produces maps of likelihood with optional contours that are used to evaluate regions for sources. Likelihood also evaluates the background corrected intensity at each location for each energy interval from which spectra can be generated. Being in a standard FITS format permits all of the maps to be easily accessed by the full complement of tools available in several commercial astronomical analysis systems. In the EGRET case, IDL is used to produce graphics plots in two and three dimensions and to quickly implement any special evaluation that might be desired. Other custom-built software, such as the spectral and pulsar analyses, take advantage of the XView toolkit for display and Postscript output for the color hard copy. This poster paper outlines the data flow and provides examples of the user interfaces and output products. It stresses the advantages that are derived from the integration of the specific instrument-unique software and powerful commercial tools for graphics and statistical evaluation. This approach has several proven advantages including flexibility, a minimum of development effort, ease of use, and portability.

Laubenthal, N. A.↗
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

Understanding interfacial chemistry of positive bias high-voltage degradation in photovoltaic modules

Photovoltaic module degradation from a high system voltage is a prevalent degradation mode in the field, where the enabling degradation mechanisms are inherently dependent on the voltage bias polarity of the installed system. In this study, the effects of positive bias on module performance are confirmed and the underlying chemical degradation processes are more thoroughly investigated to reveal different degradation pathways from those previously reported in negative bias studies. When cells are under +1000 V stress, crystalline silicon mini-modules with poly(ethylene-co-vinyl acetate) (EVA) encapsulant demonstrated a significant photocurrent loss due to EVA discoloration and delamination from increased chemical reactivity at the front-side EVA/cell metallization interface. Brown discoloration of the EVA encapsulant near the cell gridlines is linked to an electrochemical reaction at the Ag gridlines under hot and humid conditions (85 °C, 85% relative humidity). Chemical compositional analysis using X-ray photoelectron spectroscopy (XPS) confirmed that the discoloration is attributed to the formation of silver sulfide (Ag 2 S) and/or silver oxide (Ag 2 O) species at the EVA/Ag gridline interface. The subsequent migration of Ag ions from the cell gridlines into the bulk of the EVA was evident from XPS depth profiling and optical microscopy. However, the Ag signal was not detected at the EVA/glass interface, inferring limited ionic transport through the nominally 0.45 mm thick encapsulant. For the samples studied herein, the sulfur is believed by the process of elimination to come from the ambient air, diffusing into the module through the permeable polymer backsheet.

14 SOLAR ENERGY↗
DOE OSTI2025

Compatibility of Zircaloy-4 in eutectic PbSn and PbBiSn liquid metals

Here, a novel nuclear fuel concept has been proposed, consisting of UO₂ particles suspended in a liquid metal alloy. To evaluate its feasibility, the compatibility of Zircaloy-4 cladding with candidate eutectic PbSn and PbBiSn alloys—as well as with pure Pb, Bi, and Sn—was investigated through 1000-hour exposures at 400 °C and 600 °C. Severe degradation of the Zircaloy-4 specimens was observed in eutectic PbSn at 600 °C, whereas oxidation and oxygen diffusion were evident after exposure at 400 °C. At 600 °C, the formation of Zr₃Fe precipitates at the oxide–alloy interface was identified in specimens exposed to Pb and PbBiSn. Thermodynamic analysis suggested that the observed intermetallic formation resulted from the destabilization of Laves phases due to oxygen ingress at elevated temperatures.

Bismuth↗
DOE OSTI2025

Correlating Solvation Free Energy to Electrolyte Properties for Lithium Metal Batteries

The electrolyte plays a critical role in lithium metal batteries. In particular, ion solvation profoundly impacts key electrolyte properties and battery performance. Here, in this study, we systematically investigate solvation-property relationships in a series of electrolytes with different solvent-diluent ratios. Through potentiometric techniques that measure the relative solvation free energies of electrolytes, we find that weaker solvation correlates with larger ion clusters, lower ionic conductivity and diffusion coefficient, and superior electrochemical stability. Weaker solvation leads to the formation of a small number of Li clusters with large hydrodynamic radii, which lowers the Li + diffusivity and ionic conductivity of the electrolyte. Concurrently, weaker solvation leads to improved electrochemical stability at both the cathode and anode interfaces. Understanding these solvation-property relationships and trade-offs is important to designing electrolytes for optimized lithium metal battery performance.

25 ENERGY STORAGE↗
DOE OSTI2021

Dynamic Tilting of Ferroelectric Domain Walls Caused by Optically Induced Electronic Screening

Optical excitation perturbs the balance of phenomena selecting the tilt orientation of domain walls within ferroelectric thin films. The high carrier density induced in a low-strain BaTiO3 thin film by an above-bandgap ultrafast optical pulse changes the tilt angle that 90° a/c domain walls form with respect to the substrate-film interface. The dynamics of the changes are apparent in time-resolved synchrotron x-ray scattering studies of the domain diffuse scattering. Tilting occurs at 298 K, a temperature at which the a/b and a/c domain phases coexist but is absent at 343 K in the better ordered single-phase a/c regime. Furthermore, phase coexistence at 298 K leads to increased domain-wall charge density, and thus a larger screening effect than in the single-phase regime. The screening mechanism points to new directions for the manipulation of nanoscale ferroelectricity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗