Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “alloy development”

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 793 records · Page 44

Ordering effects on deformation substructures and strain hardening behavior of a CrCoNi based medium entropy alloy

Here, a CrCoNi based medium entropy alloy with small additions of Ti, Al and Nb (denoted as (CrCoNi) 93 Al 4 Ti 2 Nb) in the as-quenched condition, exhibits tensile properties comparable to those of the equiatomic CrCoNi alloy at room temperature. Dark field transmission electron microscopy (TEM), atomic resolution high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) together with atom probe tomography (APT) show that spatially-localized long range ordering (LRO) L1 2 domains exist in this alloy. The evolution of deformation substructure with plastic deformation in this alloy was characterized using electron backscatter diffraction (EBSD), electron channeling contrast imaging (ECCI) and STEM based techniques including the recently developed weak beam dark field STEM imaging. Plastic deformation occurs by the slip of a/2<110>dislocations, which are narrowly dissociated into Shockley partial dislocations on {111} slip planes. Their dissociation distances in the (CrCoNi) 93 Al 4 Ti 2 Nb alloy are much smaller than the widths of the corresponding partials in the equiatomic CrCoNi alloy due to one or more of the minor alloying elements (Al, Ti, Nb). Dislocation slip in this alloy has a pronounced planar character. The leading dislocations in slip bands glide as pairs due to the existence of LRO domains. Multipoles were formed through the slip of dislocations with opposite signs on adjacent {111} slip planes. Those multipoles serve as building blocks for the formation of subgrain structures consisting of fine slip bands. The distances between slip bands were continuously refined during plastic deformation and dynamic refinement of slip bands plays a crucial role in strain hardening. The effects of LRO domains on planar dislocation slip, the deactivation of deformation twinning and strain hardening of this alloy are discussed.

36 MATERIALS SCIENCE↗

Preliminary results addressing material qualification using combined ion irradiation and modeling data

Additively-manufactured (AM) materials have attracted increasing attention in recent years as a new method to make novel and customized components. While AM and conventionally produced materials are compositionally similar, they do possess different microstructures, necessitating assessment of materials produced via AM for their behavior in reactor environments. Some microstructures unique to AM materials, such as compositional micro-inhomogeneity and dislocation cell structures, are of particular importance since they may lead to different radiation performance. The performance of AM materials for advanced nuclear reactor applications is of interest to the Advanced Materials and Manufacturing Technologies (AMMT) program under the Department of Energy Office of Nuclear Energy. The AMMT program aims to demonstrate its new accelerated development and qualification methods via laser powder bed fusion (LPBF) 316 stainless steel (SS). Focusing on material bearing both 316L and 316H specifications, we integrate ion irradiation and modeling. This year, we focus on answering foundational questions related to process variability, alloy chemistry variation, and microchemical segregation. Experimental results provide information and motivate questions to the modeling effort, which aims to develop the ability to model radiation-driven microstructural evolution in additively-manufactured 316 stainless steel under a variety of advanced reactor conditions, including different temperatures, neutron spectra, and fluxes, in a sort of "virtual experiment". We perform in-situ and ex-situ ion irradiations and microstructural characterizations to support the development of AM materials for reactor applications, develop a phase field model of radiation-induced segregation in additively manufactured material with high angle grain boundaries and dislocation cells, and investigate the effect of carbon and chromium content on point defect behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermocouple Testing in Support of the AGR-5/6/7 Experiment

This report documents thermocouple testing performed in IRC Lab C-15 over a period of seven years. This testing supported selection and characterization of the thermocouple set used in the AGR-5/6/7 experiment. The following summary was taken directly from the report. Temperature measurement is a challenging aspect of very high temperature irradiation experiments because commonly used high-temperature commercial thermocouples such as platinum-rhodium (Types S, R, and B) and tungsten-rhenium (Type C), suffer dramatic drift because of neutron-induced transmutation. As a result, these types of thermocouples, which are used routinely for industrial temperature measurements outside of reactors, are used only in very special circumstances for reactor experiments. Conversely, because of their low neutron cross-sections, Type N thermocouples are affected to only a limited extent by neutron irradiation. However, the use of these nickel-based thermocouples is limited when the temperature exceeds 1050°C due to drift arising from minor alloying elements migrating from the thermocouple's metal sheath to the thermoelements. This change in the composition of the thermo-elements results in significant decalibration of the signal. The issues described above were recognized during the early planning stages of the final AGR experiment (designated AGR-5/6/7), and a thermocouple furnace testing program was performed over a seven-year period (2014-2019, 2021) to first select and then characterize the best thermocouple set for the high temperature regions of the AGR-5/6/7 irradiation experiment. The calculated temperature range of the AGR-5/6/7 experiment was 600–1500°C. For temperatures below 1000°C standard Type N thermocouples were deemed adequate. The furnace testing campaign identified two thermocouple types suitable for measuring temperatures above 1000°C, a Mo/Nb thermocouple developed at INL called HTIR-TC, and a Type N thermocouple developed by Cambridge University (called herein Cambridge Type N), which featured a custom high nickel alloy sheath. One of the original goals of the furnace testing program was to identify a thermocouple capable of low drift operation near the peak temperature expected in AGR-5/6/7, i.e., about 1400°C. The HTIR-TC design appeared promising in this regard, however a manufacturing difficulty proved to be a barrier and instead the furnace testing focused on drift performance at 1250°C. The manufacturing difficulty was that the Nb sheaths of the HTIR-TCs experienced extreme embrittlement when heat treated at 1600°C or greater. Heat treatment is needed to stabilize the emf output of this TC type, and the higher the heat treatment temperature the higher the peak temperature of stable operation. Because of the sheath embrittlement the heat treatment temperature had to be lowered to 1450°C resulting in a stable operating temperature of about 1250°C. One of the successes of the furnace testing program was identification of a shortcoming in the heat treatment procedure that had been traditionally used in the production of HTIR-TCs. The shortcoming was that the entire heated length of the HTIR-TC sensor was not being heat treated, but rather only the part of the sensor expected to experience temperatures above 1000°C. The problem manifested itself when the thermocouples were removed from the heat treat furnace and placed in another furnace with a different geometry, their indicated temperatures would be widely scattered, but mostly in the negative direction. The solution was to heat treat the entire heated length of the sensor. Since the deepest immersion depth in the AGR-5/6/7 experiment was about 40 inches, a heat treatment length of 48 inches was used. After this change was implemented, thermocouples which were moved into a new environment with a different temperature profile (i.e., a different furnace), produced accurate temperature measurements. Although assembly of the AGR-5/6/7 experiment was completed in September of 2017 (and irradiation begun in 2018), furnace testing of thermocouples continued in 2018 and 2019. The main purpose of this testing was to establish very long-term drift characteristics of the HTIR and Cambridge Type N thermocouples installed in the experiment. Representative thermocouples from the same lots as those installed in the AGR-5/6/7 experiment were used. Additionally, thermocouples of different designs, (particularly variations on the HTIR-TC design) were "piggy-backed" on this testing program to provide insights for instrumenting future very high temperature irradiation experiments. This two-year testing program demonstrated that HTIR-TCs and Cambridge Type N TCs could operate at 1250°C for up to 10,000 hrs (and in some cases longer) while experiencing negative drifts on the order of 2-4°C/1000 hrs. This performance was considered acceptable given the extreme operating environment the sensors faced.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Microstructure-sensitive modeling of high temperature creep in grade-91 alloy

Predicting the effects of microstructure on high-temperature creep responses of steel components is critical to minimize the risks of failure and maximize economic viability in the energy sector. In this research, a recently developed advanced mechanistic constitutive model is employed to study the effect of microstructure on the creep responses and to rationalize experimentally reported variability in the performances of grade-91 alloy. Hundreds of experimental creep tests from the literature are used to assess the predictability of the model. The model proposed is shown to accurately predict the steady-state creep rates and also the less commonly considered, yet important, primary-to-secondary transients for a wide range of creep conditions. Further, the model is capable of extrapolating the creep response of grade 91 alloy even outside of the calibration regime. Using this model, the roles of initial microstructure described in terms of grain size, dislocation density, and precipitate content, on the creep responses are investigated. The effect of initial microstructure on the steady-state creep rates and creep-strain is found to vary with stress and temperature. Grain size plays a significant role in the low-stress creep, where the diffusional creep is dominant. On the other hand, the effects of dislocation densities and precipitate content are significant in the dislocation-plasticity-dominated high-stress and temperature regimes. Furthermore, by comparing model predictions against a large experimental creep database, it is found that variability in the creep responses can be rationalized on the basis of differences in the initial microstructure. Overall, this work provides a robust pathway for microstructure-aware engineering design, verification, and validation of metallic components for energy applications.

36 MATERIALS SCIENCE↗

High-Temperature Creep Test Suite for Grizzly

A recent development in the Grizzly has focused on improving its ability to model the high-temperature response of alloys, enabling its use for predictive modeling of material response for advanced reactor applications. This report documents an effort to expand the example and assessment test suites in Grizzly to include a set of problems that exercise the models for high-temperature response of components relevant to advanced reactors. These problems can be used to assess the accuracy of models, provide examples for users interested in modeling this class of problems, monitoring performance on large-scale problems and as test cases for future feature development.

36 MATERIALS SCIENCE↗

The effects of dose, dose rate, and irradiation type and their equivalence on radiation-induced segregation in binary alloy systems via phase-field simulations

Radiation-induced segregation is a phenomenon commonly observed in many alloys which consists of the redistribution of elements (solute or interstitial impurities) under irradiation. The onset and development of radiation-induced segregation can only occur when a sufficient flux of defects is sustained and defect sinks are present. Irradiation dose, dose rate, and particle types all affect defect flux. In this work, we employ a phase-field model to examine the effects of dose, dose rate, and type of incident particles on radiation-induced segregation behavior in a model binary alloy. The phase-field model takes into account the formation and evolution of point defects as well as defect clusters, the diffusion and clustering of alloy species, the presence of additional extrinsic defect sinks in the form of dislocations, and two different methods of radiation-damage insertion, which are intended to simulate either light-ion/electron irradiation via Frenkel pairs or heavy-ion irradiation in the form of cascades. Our results show a dose-rate and particle-type dependence on the amount of solute segregation. We show that the material systems exposed to higher dose rates are less subjected to solute segregation at equivalent doses. We also show that such dose-rate-dependence behavior is due to a delay of the incubation dose at which radiation-induced segregation effectively starts. Particle type and the presence of dislocations can accentuate this behavior. Our model predictions correlate with many experimental observations made over the years on radiation-induced segregation providing credence to the simulation results. The methodology presented in this study allows for a first-order prediction of the dose rate at which proxy irradiation experiments could be performed to approximate radiation-induced segregation behaviors seen in targeted irradiation conditions.

36 MATERIALS SCIENCE↗

Predicting Rupture Pressure and Leak Rate of Alloy 600 U-bend Tubes with PWSCC and Comparing with Test Results (Technical Letter Report)

This report summarizes structural integrity work performed for the US NRC that predicts the structural integrity of steam generator (SG) tubes with flaws in the U-bend. The work covers model development and validation for U-bend tubes containing lab-induced primary water stress corrosion cracks (PWSCC). Limit load models for plastic collapse predicting ligament rupture pressure and leak rate of cracked U-bend tubes were developed and validated with pressure tests and finite element analyses (FEAs). Failure pressure and leak rate tests were conducted at room temperature on Alloy 600 U-bend tubes containing flaws located at various locations around the U-bend. Existing models for predicting structural integrity of straight tubes with well defined, rectangular, electro-discharge-machined (EDM) flaws were modified and used to predict the ligament rupture and unstable burst pressures as well as leak rates for the PWSCC flaws in the U-bend section.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Engineering Multi-scale B2 Precipitation in a Heterogeneous FCC Based Microstructure to Enhance the Mechanical Properties of a Al0.5Co1.5CrFeNi1.5 High Entropy Alloy

While ordered L12 or gamma prime precipitates in face centered cubic (FCC) based microstructures have been extensively used for strengthening Nickel or Cobalt base superalloys, and more recently in high entropy alloys (HEAs) or complex concentrated alloys (CCAs), the possibility of exploiting ordered B2 precipitates in FCC-based systems has been relatively less investigated. The present study shows the propensity of developing a heterogeneous microstructure, consisting of two different distributions of FCC grain sizes, and two different size scales of B2 precipitates, within an FCC-based Al0.5Co1.5CrFeNi1.5 HEA/CCA. This alloy composition has been designed using solution thermodynamics-based modeling such that it has a high phase fraction and solvus temperature of the B2 phase. The resulting heterogenous microstructure exhibited an approximately 400% increase in yield strength with respect to the single-phase FCC solid solution condition of the same alloy while maintaining very good tensile ductility ~20 %.

Dasari, Sriswaroop↗

High-Temperature Oxidation Behavior of FeCoCrNi+(Cu/Al)-Based High-Entropy Alloys in Humid Air

Previous studies showed some transition metal high-entropy alloy (HEA) compositions can have good oxidation resistance in air up to 800 °C. Four equiatomic HEAs have been developed based on FeCoCrNi with additions of Mn, Cu, Al or Al+Cu. The oxidation behavior of these HEAs was compared in humid (10 vol.% H2O) air at 800 °C for 100–500 h to investigate the influence of water vapor on the oxidation mechanisms. The Cu- and Al-containing alloys exhibited improved oxidation resistance over the Mn composition. For the Cu-containing alloy, a local attack of the Cu-rich phase was observed, which formed an Fe/Ni/Co/Cr spinel that was surrounded by Cr2O3. This oxide was thicker for the humid air atmosphere when compared to dry air, and the transition of the Cu oxide to the spinel was accelerated. The Al-containing HEA formed a thin Al2O3 scale with humidity suppressing AlN formation and forming a smoother oxide layer. The Al+Cu composition had the highest overall oxidation resistance (minimal local attack, no nitridation) and also showed a smooth oxide scale topography under humid air oxidation as opposed to a plate-like, rougher scale under dry air.

Crystallography↗

Salt-Assisted 2H-to-1T' Phase Transformation of Transition Metal Dichalcogenides [plus supplemental information]

Phase engineering of nanomaterials (PEN) has demonstrated great potential in the fields of catalysis, electronics, energy storage and conversion, and condensed matter physics. Recently, transition metal dichalcogenides (TMDs) with unconventional metastable phases (e.g., 1T and 1T') have attracted increasing research interest due to their unique and appealing physicochemical properties. However, there is still a lack of a simple, universal, and controlled method for the preparation of large-scale and high-purity unconventional-phase TMD crystals, restricting their further fundamental study and practical applications. Here, a facile, one-step salt-assisted general strategy is reported for the controlled phase transformation of commercially available TMDs with conventional 2H phase, yielding a large amount of metastable 1T'-phase TMDs, including WS 2 , WSe 2 , MoS 2 , and MoSe 2 . It is found that the easily accessible metal salts, such as K 2 C 2 O 4 ·H 2 O, K 2 CO 3 , Na 2 CO 3 , Rb 2 CO 3 , Cs 2 CO 3 , KHCO 3 , NaHCO 3 , and NaC 2 O 4 , can be used to assist the 2H-to-1T' phase transformation, greatly simplifying the synthetic process for producing metastable 1T'-TMDs. Importantly, this method can also be used to prepare 1T'-TMD alloys, such as 1T'-WS 2x Se 2(1–x) . We report this newly developed strategy is robust and highly effective, which can also be used for the phase engineering of other materials with various polymorphs.

1T′ phase↗

The effect of Cr alloying on defect migration at Ni grain boundaries

Mass transport along grain boundaries in alloys depends not only on the atomic structure of the boundary, but also its chemical make-up. In this work, we use molecular dynamics to examine the effect of Cr alloying on interstitial and vacancy-mediated transport at a variety of grain boundaries in Ni. We find that, in general, Cr tends to reduce the rate of mass transport, an effect which is greatest for interstitials at pure tilt boundaries. However, there are special scenarios in which it can greatly enhance atomic mobility. Cr tends to migrate faster than Ni, though again this depends on the structure of the grain boundary. Further, grain boundary mobility, which is sometimes pronounced for pure Ni grain boundaries, is eliminated on the time scales of our simulations when Cr is present. We conclude that the enhanced transport and grain boundary mobility often seen in this system in experimental studies is the result of non-equilibrium effects and is not intrinsic to the alloyed grain boundary. These results provide new insight into the role of grain boundary alloying on transport that can help in the interpretation of experimental results and the development of predictive models of materials evolution.

36 MATERIALS SCIENCE↗

A phase-field study of stainless-steel oxidation from high-temperature carbon dioxide exposure

An electrochemical phase-field model has been developed to investigate the oxidation mechanisms of the 21-2N valve stainless steel alloy exposed to carbon dioxide at 973 K. Three oxide phases observed in oxidation experiments are included in the model: Mn 3 O 4 , Cr 2 O 3 , and MnCr 2 O 4 . Local charge neutrality and the conserved current condition are assumed to include the impact of the electric potential on oxidation, though it is found to be negligible when the system is electrically isolated. The sensitivity of the oxidation processes to the values of the diffusion mobilities is examined, and the oxidation rate is calibrated against experimental data by modifying the sensitive mobilities. We find that both inward oxygen and outward metal diffusion are important for oxidation. Furthermore, by investigating the impact of the order of the initial oxide layers, we find that Cr 2 O 3 serves as a barrier to outward Mn diffusion. Moreover, controlling the initial oxide layer order can be used to limit oxidation.

36 MATERIALS SCIENCE↗

Formation and characterization of hydride rim structures in Zircaloy-4 nuclear fuel cladding tubes

Zirconium alloy tubes used as nuclear fuel cladding are subject to oxidation and subsequent hydrogen pickup during their long service in commercial light water reactors. The hydrogen picked up in the cladding can precipitate as a brittle hydride rim feature on the cladding outer surface. Here, to better understand the effect of hydride rims on the fracture behavior of stress-relieved zirconium alloy cladding tubes, a procedure to produce these rim-like structures has been developed and is described herein. Extensive characterization of the ‘as hydrided’ tubes is performed. The hydrogen charging apparatus consists of a tube furnace with a quartz chamber that is connected to a vacuum pump as well as a bottle of pure hydrogen gas. The charging station uses a static charge of hydrogen as opposed to a flowing gas. The chief advantage of this approach is the ability to monitor the pressure drop in the hydriding chamber and correlate this pressure drop to a known rate of hydrogen pickup in the cladding tube. It was found that the hydrogen partial pressure, metal temperature, and surface treatment all clearly played a role in whether a hydride rim was formed. Extensive characterization of the hydride rims shows they consist of needle like platelets of δ phase hydrides (ZrH1.66) oriented in the circumferential direction with a radial spacing of several microns in a sandwich-like structure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

From classical thermodynamics to phase-field method

Phase-field method is a density-based computational method at the mesoscale for modeling and predicting the temporal microstructure and property evolution during materials processes. The focus of this article is on connecting the most common phase-field equations to the very basic first and second laws of classical thermodynamics through rudimentary irreversible thermodynamics. It briefly discusses the relations of the continuum phase-field equations to their counter parts at the microscopic and atomic levels. It attempts to clarify the contributions of long-range elastic, electrostatic, and magnetic interactions to domain structure evolution during structural, ferroelectric, and ferromagnetic phase transformations by separating order parameter changes due to the presence of quasi-static fields and those arising from phase transformations. A few examples are presented to demonstrate the possibility of employing the phase-field method to provide guidance to designing materials for optimum properties or discovering novel mesoscale phenomena or new materials functionalities. Here, the article ends with a brief perspective on a number of potential future directions on the development and applications of phase-field method beyond its traditional applications to structural alloys.

36 MATERIALS SCIENCE↗

Characterizing defect structures in AM steel using direct electron detection EBSD

The mechanical properties of additive and traditionally manufactured alloys are largely dependent on the characteristics and distribution of dislocation cell networks that develop during the fabrication process. This work demonstrates the ability to quantitatively characterize these dislocation structures by high angular resolution electron backscatter diffraction analysis using a direct electron detector. The defect structures are characterized in terms of the geometrically necessary dislocation density and the associated Burgers vector and line direction. Furthermore, the results are discussed in terms of potential defect formation mechanisms.

36 MATERIALS SCIENCE↗

A crystal plasticity finite element method modeling of Zircaloy with hydride phases based on scanning electron micrographs

This work developed a microstructure-based finite element model to predict the stress state of alloys with second phase inclusions. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. Generation of digital microstructure was achieved through two steps of tessellations using software Neper and Matlab. The process is demonstrated using an example of Zircaloy material with secondary phases of hydrides dispersed within and stress-strain response of Zircaloy containing hydrides was predicted. The constitutive material model for Zircaloy in this study was based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as isotropic elasto-plastic material. A parametric study had been conducted to understand the effect of volume fraction, and lamellae thickness of the hydride phase on the mechanical properties of the overall material. Results can help designers to alter the manufacturing process to obtain the enhanced mechanical properties for components used in nuclear applications made by Zircaloy material.

Kulkarni, Shank S.↗

Physics of wurtzite ferroelectrics

Ferroelectricity was long considered incompatible with the wurtzite structure, but the recent discovery of switchable polarization in wurtzite alloys has renewed interest in these materials for integrated electronic and memory applications. The development of wurtzite ferroelectrics faces significant technological challenges, which can be addressed through a fundamental physical understanding of their dielectric and ferroelectric properties. This article focuses on the physics that govern the polarization switching behavior, emphasizing the atomic- and meso-scale (domain) mechanisms involved in the transition between polarization states. A distinguishing feature of this article is a deep dive into the role of intrinsic and extrinsic defects—an area that has received limited attention in prior reviews, but is increasingly recognized as central to polarization switching, coercive fields, leakage, and fatigue. We highlight how defect behavior evolves during processing and electrical cycling, often contributing to long-term degradation. We also introduce powerful first-principles defect calculations, common in semiconductors but not yet widespread in ferroelectrics, as tools to understand and design materials. By integrating recent theoretical and experimental insights, we aim to provide a framework for advancing wurtzite ferroelectrics.

36 MATERIALS SCIENCE↗