Search NASA⌕ Search

SEARCH · Search NASA

Results for “zone refining”

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

Impacts of spatial heterogeneity of anthropogenic aerosol emissions in a regionally refined global aerosol–climate model

Abstract. Emissions of anthropogenic aerosol and their precursors are often prescribed in global aerosol models. Most of these emissions are spatially heterogeneous at model grid scales. When remapped from low-resolution data, the spatial heterogeneity in emissions can be lost, leading to large errors in the simulation. It can also cause the conservation problem if non-conservative remapping is used. The default anthropogenic emission treatment in the Energy Exascale Earth System Model (E3SM) is subject to both problems. In this study, we introduce a revised emission treatment for the E3SM Atmosphere Model (EAM) that ensures conservation of mass fluxes and preserves the original emission heterogeneity at the model-resolved grid scale. We assess the error estimates associated with the default emission treatment and the impact of improved heterogeneity and mass conservation in both globally uniform standard-resolution (∼ 165 km) and regionally refined high-resolution (∼ 42 km) simulations. The default treatment incurs significant errors near the surface, particularly over sharp emission gradient zones. Much larger errors are observed in high-resolution simulations. It substantially underestimates the aerosol burden, surface concentration, and aerosol sources over highly polluted regions, while it overestimates these quantities over less-polluted adjacent areas. Large errors can persist at higher elevation for daily mean estimates, which can affect aerosol extinction profiles and aerosol optical depth (AOD). We find that the revised treatment significantly improves the accuracy of the aerosol emissions from surface and elevated sources near sharp spatial gradient regions, with significant improvement in the spatial heterogeneity and variability of simulated surface concentration in high-resolution simulations. In the next-generation E3SM running at convection-permitting scales where the resolved spatial heterogeneity is significantly increased, the revised emission treatment is expected to better represent the aerosol emissions as well as their lifecycle and impacts on climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Area V: Waste Management and Subsurface Science (FY21-23 Final Report)

The objective of the Waste Management and Subsurface Science Topic Area V project is to evaluate the feasibility and safety of an intermediate-depth borehole for disposal of Israel Atomic Energy Commission (IAEC) nuclear waste. The IAEC nuclear waste is currently stored in above ground and/or in near-surface facilities that provide short-term radiological worker safety but do not provide long-term, full security and safeguards protection for the nuclear material. Disposal of the IAEC nuclear waste several hundred meters below the surface in a sealed borehole would provide long-term security and safeguards due to its depth underground, as well as long-term radiological safety, due to its isolation in a stable geological formation. This project supports the nuclear security and safeguards initiatives of the Science and Technology Working Group of the NNSA/IAEC Memorandum of Understanding (MOU) and addresses the Topic Area V objective and associated goals. The goals for Topic Area V are: (1) Assess the suitability of the vadose zone in the Yamin Plain of the northern Negev Desert, Israel as a host location for intermediate-depth borehole disposal of radioactive waste, (2) Define and refine the necessary characterization needs to better understand the regional climate, structure, and stratigraphy, geomechanics, geochemistry, hydrogeology, and seismicity of the stratigraphic sequence from the surface to the water table in the vicinity of the proposed host location for the disposal borehole(s), and (3) Develop a preliminary plan and schedule, including necessary resources and costs, for a borehole disposal concept and robust and transparent pre-closure and post-closure safety assessments. The safety assessments will require iteration with process level models at both site and regional scales. Upon completion, the information collected by this effort will allow the IAEC to make decisions regarding the viability and applicability of an intermediate depth (<500m) borehole for disposal in the vadose zone of the Yamin Plain. This FY21-23 final report is provided in lieu of a workshop report that was originally planned for delivery in January, 2024. Due to the conflict in Israel and associated delays, the workshop has been delayed to April 2024. The workshop will be held at LLNL but offsite and will switch to hybrid if travel restrictions remain in place. The associated workshop report will be provided as part of a new FY24-26 LCP.

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

Regional Geologic Cross Sections for Potential Storage and Containment Zones in the MRCSP Region [plus Plates 1-6]

The regional characterization work conducted by the Geoteams during the MRCSP Phase III project period (2010 –2019) focused on the following tasks: (1) refinement of geologic seals/reservoir systems; (2) assessment of Atlantic Coastal Plain and offshore opportunities; (3) expanded assessments of oil and gas fields, particularly as they relate to enhanced recovery opportunities; (4) regional support for implementation of carbon capture utilization and storage (CCUS) in the partnership area; and (5) communication and data sharing. This report, entitled Regional Geology, has been prepared in association with (4), regional support for implementation of CCUS in the partnership area..

58 GEOSCIENCES↗
DOE OSTI2024

Early experiences on the OLCF Frontier system with AthenaPK and Parthenon–Hydro

The Oak Ridge Leadership Computing Facility (OLCF) has been preparing the nation's first exascale system, Frontier, for production and end users. Frontier is based on HPE Cray's new EX architecture and Slingshot interconnect and features 74 cabinets of optimized 3rd Gen AMD EPYC CPUs for HPC and AI and AMD Instinct 250X accelerators. As a part of this preparation, “real-world” user codes have been selected to help assess the functionality, performance, and usability of the system. This article describes early experiences using the system in collaboration with the Hamburg Observatory for two selected codes, which have since been adopted in the OLCF test harness. Experiences discussed include efforts to resolve performance variability and per-cycle slowdowns. Results are shown for a performance portable astrophysical magnetohydronamics code, AthenaPK, and a mini-application stressing the core functionality of a performance portable block-structured adaptive mesh refinement framework, Parthenon-Hydro. Here, these results show good scaling characteristics to the full system. At the largest scale, the Parthenon-Hydro miniapp reaches a total of $1.7$ $\times$ $10^{13}$ zone-cycles/s on 9216 nodes (73,728 logical GPUs) at ≈92% weak scaling parallel efficiency (starting from a single node using a second-order, finite-volume method).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Technoeconomic Design Optimization for Fast Reactors. Part II: Impact of Technoeconomic Constraints on Optimal Design

There is a current drive toward optimizing reactors, particularly small/micro reactors to minimize cost and maximize performance. Previous work has investigated the development of technoeconomic workflows for the design optimization of pool-type fast reactors that aim to deploy into district energy grids. Initial scoping studies verified that the workflow was capable of capturing design trends throughout a variety of design configurations and problem formulations while remaining sufficiently flexible. In this paper, this methodology is applied to understand how cost functions and technoeconomic constraints can drive optimal reactor design. Specifically, the UPu10Zr-fueled fast reactor model from Part I is adapted to include changes in the fissile content limits, control rod worth limits, control rod drive cost, and assumed fuel form. In the case of constraint relaxation at fixed power (fissile content and control rod worth limits), cost sensitivities of 5% to 10% were uncovered. Multi-objective optimization at varying reactor power levels with individualized control rod drives for each assembly (as opposed to one operational and one safety drive) increased cost by $\$10$ to $\$25$ million and substantially altered the optimal core geometry, favoring geometries with substantially fewer control rod placements relative to baseline optimization. Finally, a multi-objective optimization was performed at varying power levels with the fuel form overhauled to metallic, high-assay low-enriched uranium–based U10Zr with more refined fuel cost models. In the case of uranium fueling, the costs increased by at least $50 million relative to the baseline case. Furthermore, economic fuel zoning and lower reactivity swing cores were recovered. Each case serves to demonstrate the value of applying technoeconomic workflows to initial reactor design scoping studies to better understand the trade-off for a proposed concept between different design options.

Argonne Reactor Computation (ARC) codes↗
DOE OSTI2020

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

High resolution ferromagnetic and antiferromagnetic correlations in La 0.7 Ca 0.3 MnO 3

La 0.7 Ca 0.3 MnO 3 (LCMO30) is a perovskite manganite with a phase transition (T C = 238 K) between a paramagnetic insulator and a ferromagnetic metal. We use single-crystal, inelastic neutron scattering to refine previous work on the low energy (ℏω < 10 meV) magnetic excitations. At 300 K the spectrum consists of a diffusive ferromagnetic component at the pseudo-cubic Brillouin zone center coexisting with a quasielastic excitation located at the pseudo-cubic Brillouin zone edge. Both components are anti-correlated with ferromagnetic long-range order, showing weakened intensity at 215 K in the ferromagnetic phase. Polarized neutron scattering shows both features to be purely magnetic in origin. Furthermore, this indicates the simultaneous presence of ferromagnetic and antiferromagnetic correlations (located, respectively, at the zone center and zone edge) in a three-dimensional perovskite, while this coexistence had been previously observed only in a bilayer compound.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of die design on microstructure and mechanical joint strength in friction self-piercing riveted AA7055-T76 and AA7055-T76

In the present study, a unique friction self-piercing riveting (F-SPR) was applied to join high-strength, low-ductility aluminum alloy (AA) 7055-T76 and AA7055-T76 to mitigate a cracking issue. Also, the effect of die design on both joint formation and the mechanical joint performance were investigated. A crack-free joint was achieved when joining AA7055-T76 and AA7055-T76 by F-SPR. As the bottom diameter of the die increased, the mechanical interlocking distance initially increased then became almost saturated. At the same time, the solid-state bonding gradually disappeared, and the lap-shear fracture changed from bottom aluminum fracture to rivet pull-out with partial bottom aluminum fracture. The average maximum lap-shear joint strength of 11.57 kN and cross tension of 5.65 kN were achieved from the D2 die design owing to the combined contributions from high mechanical interlocking distance and good solid-state bonding at the joint interface. Ultrafine grain refinement was observed at the region next to the rivet shank outer surface and along the materials flow line between the top and bottom aluminum sheet in the rivet cavity. Further, because of the same joint-process parameters, the hardness profile at the joints made by different dies did not show obvious differences. As the distance from the base metal to the rivet outer surface decreased, frictional heat caused the hardness to initially decrease from 190 to 150 HV in the heat-affected zone. Hardness increased up to 175 HV in the thermomechanical affected zone (TMAZ-I) owing to large plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structural investigation of the temperature-stable relaxor dielectric Ba 0.8 Ca 0.2 TiO 3 -Bi(Mg 0.5 Ti 0.5 )O 3

Aberration corrected scanning transmission electron microscopy (STEM) and electron diffraction have been used to disclose local structure and nano-chemistry in a Ca modified BaTiO 3 -Bi(Mg 0.5 Ti 0.5 )O 3 relax or dielectric ceramic which exhibits high and near-invariant relative permittivity over a wide temperature range. High resolution, synchrotron X-ray diffraction indicated a globally cubic structure ($Pm$$\overline{3}$$m$), but direct atomic-scale imaging by STEM revealed local tetragonal distortions. Nanopolar clusters were identified from B-site atomic displacement vectors measured relative to oxygen ion positions along < 100 > and < 110 > zone axes of integrated differential phase contrast (iDPC) STEM images, highlighting cluster sizes of 2–5 nm. Chemical analysis by STEM-energy dispersive X-ray spectroscopy and full pattern refinements of X-ray powder diffraction data each implied high levels of Bi vacancies within the matrix. The possibility that these A-site vacancies modulate the nanopolar structure and promote flattening of the permittivity-temperature response in this class of dielectric is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Texture and precipitation behavior of friction stir welded Inconel 825 alloy

The microstructure, texture pattern, and precipitation behavior of friction stir butt welded Inconel 825 alloy are highlighted here. The microstructure of bright regions and black-banded regions in the stir zone (SZ) is characterized by SEM, EBSD, and TEM with EDS. EBSD analysis justifies that the SZ exhibits a dynamic recrystallized structure. Grain boundary mapping corroborates the insignificant difference in the low- and high-angle grain boundaries fraction between the SZ, the thermo-mechanically affected zone, and the base metal. TEM analysis confirms the formation of nano twins and dispersion of nano-size precipitations in the grains, and at the grain boundaries as well. Grain refinement and nano-size precipitation formation in the SZ significantly enhanced both microhardness and tensile strength.

Das, Hrishikesh↗
DOE OSTI2024

Quantum mechanical model of crossing and anti-crossing points in 3D full-band Monte Carlo simulations

This work presents a 3D quantum mechanics based model to address the physics at band structure crossing/anti-crossing points in full band Monte Carlo (FBMC) simulations. The model solves the Krieger and Iafrate (KI) equations in real time using pre-computed coefficients at k-points spatially sampled within the first Brillouin zone. Solving the KI equations in real time makes this model applicable for all electric fields, which enables its use in FBMC device simulations. In this work, a two-level refinement scheme is used to aggressively sample regions in proximity to band crossings for accurate solutions to the KI equations and coarsely sample everywhere else to limit the number of k-points used. The presented sampling method is demonstrated on the band structure of silicon but is effective for the band structure of any semiconductor material. Next, the adaptation of the fully quantum KI model into the semi-classical FBMC method is discussed. Finally, FBMC simulations of hole transport in 4H silicon carbide with and without the KI model are performed. Results along different crystallographic directions for a wide range of electric fields are compared to previously published simulation and experimental values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A systematic decision-making methodology to formalize the selection of degree of realism in screening analysis of probabilistic risk assessment

In the nuclear power domain, Probabilistic Risk Assessment (PRA) is used to inform decision-making for Nuclear Power Plants (NPPs). Recently, there has been an increase in the utilization of modeling and simulation (M&S) to support the estimation of PRA inputs. Risk analysts should carefully select the PRA items that require M&S and their degree of realism (DoR) with consideration of the required resources. To support this selection, this article formulates a systematic decision-making approach for the DoR selection. The DoR selection is made based on two predictive decision-making attributes: the predicted differences in safety risk estimate (ΔSaRi) and the cost of analysis (ΔCAN). This research also develops and quantifies causal models to estimate ΔSaRi and ΔCAN. The causal model-based prediction of ΔSaRi and ΔCAN helps reduce the trial-and-error nature of the DoR selection in the PRA screening analysis and provides insights for DoR selection and the gradual refinements of PRA realism. This approach is demonstrated for a case study on fire PRA of NPPs, where an adequate DoR is selected from two fire models: an engineering correlation and a zone model.

Alkhatib, Sari [Department of Nuclear, Plasma, and↗
DOE OSTI2021

Impact of soiling on Si and CdTe PV modules: Case study in different Brazil climate zones

Soiling, particulate accumulation on photovoltaic (PV) module surfaces, reduces the available solar resource and the resulting generated device power. This case-study summarizes initial results of 5-year research on the contrasting soiling conditions in the tropical, subtropical, and semi-arid climates in Brazil. A major objective is to present a case study of the effects of soiling on PV module performance in different Brazil climate zones that represent the primary areas for the current and expanding Brazil solar installations. For this, the paper presents methodologies to quantify the soiling ratio (SRatio) and soiling rate (SRate) for two representative commercial technologies, polycrystalline or multicrystalline silicon (mc-Si) and thin-film cadmium telluride (CdTe) modules, through soiling monitoring stations deployed in the selected climate regions. An aim is to add to the growing soiling-research knowledge base through addressing these key factors and their relationships to critical electrical, solar resource, thermal, and local meteorological and environmental parameters. This paper presents, evaluates, and compares soiling rates and losses in Belo Horizonte, Minas Gerais (Equatorial Tropical: 19.92° S, 43.99° W), Porto Alegre, Rio Grande do Sul (Humid-Subtropical: 30.05° S, 51.17° W), and Brotas de Macaúbas, Bahia (Semi-Arid: 12.00° S, 42.63° W). The results show that soiling is moderate in all 3-regions, for example with 0.1%/day < SRate < 0.2%/day for Belo Horizonte. Precipitation dominates the cleaning of the modules in the summertime in this climate zone, while it is the major factor year-round in Rio Grande do Sul. Wind is the major issue mitigating the soiling accumulation for the Bahia installation. The methodology incorporates several key refinements, including the normalization and adjustment for the meteorological parameters (temperature, irradiance, wind, precipitation). The evaluations include the region-specific differing effects of non-uniform soiling, natural cleaning, and ambient temperatures.

14 SOLAR ENERGY↗
DOE OSTI2023

A Time-Gated, Time-Correlated Single-Photon-Counting Lidar to Observe Atmospheric Clouds at Submeter Resolution

Most lidars used for cloud observations have the range resolution of about 10 m, so they are incapable of resolving submeter-scale processes that are crucial to cloud evolution. This article describes a prototype of a ground-based, vertically pointing, time-gated, time-correlated single-photon-counting lidar (referred to as the T2 lidar) developed to explore atmospheric clouds at range resolution two orders of magnitude finer than traditional atmospheric lidars. The T2 lidar emits green-light pulses (532 nm) at a repetition rate of 20.6 kHz and a pulse width of ∼650 ps, which enables the observation of aerosol and cloud layers at heights from a few hundred meters to 7.28 km above the ground level at range resolution down to 10 cm. In addition, a digital delay pulse generator controls the detector to only receive photons for a short period after each laser pulse. This time-gated technique blocks photons arriving from regions outside the target zone, thus significantly reducing the noise level and allowing observation even inside clouds. Initial observations show that the T2 lidar can detect sharp cloud boundaries and fine structures near the cloud base. Such refined measurements of cloud structure could lead to an improved understanding of microphysical processes such as droplet activation, entrainment and mixing, and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Nitrous oxide inhibition of methanogenesis represents an underappreciated greenhouse gas emission feedback

Methane (CH 4 ) and nitrous oxide (N 2 O) are major greenhouse gases that are predominantly generated by microbial activities in anoxic environments. N 2 O inhibition of methanogenesis has been reported, but comprehensive efforts to obtain kinetic information are lacking. Using the model methanogen Methanosarcina barkeri strain Fusaro and digester sludge-derived methanogenic enrichment cultures, we conducted growth yield and kinetic measurements and showed that micromolar concentrations of N 2 O suppress the growth of methanogens and CH 4 production from major methanogenic substrate classes. Acetoclastic methanogenesis, estimated to account for two-thirds of the annual 1 billion metric tons of biogenic CH 4 , was most sensitive to N 2 O, with inhibitory constants (K I ) in the range of 18–25 μM, followed by hydrogenotrophic (K I , 60–90 μM) and methylotrophic (K I , 110–130 μM) methanogenesis. Dissolved N 2 O concentrations exceeding these K I values are not uncommon in managed (i.e. fertilized soils and wastewater treatment plants) and unmanaged ecosystems. Future greenhouse gas emissions remain uncertain, particularly from critical zone environments (e.g. thawing permafrost) with large amounts of stored nitrogenous and carbonaceous materials that are experiencing unprecedented warming. Incorporating relevant feedback effects, such as the significant N 2 O inhibition on methanogenesis, can refine climate models and improve predictive capabilities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Structural Characterization of Potential Carbon Dioxide Reservoirs and Adjacent Strata within the Llandovery Silurian to Middle Devonian Strata of Ohio

The Midwest Regional Carbon Sequestration Partnership (MRCSP) has incorporated the work of geologic research teams (Geoteams) in its regional characterization, project planning and carbon dioxide (CO 2 ) injection implementation work since the partnership was established by the U.S. Department of Energy (DOE) in 2003. Over this 16-year period, the cohort of Geoteams has grown from five to ten states and has contributed to the characterization of geologic sequestration opportunities, refinement of reservoir and seal data, and supported injection efforts through both predictive and post-injection assessments. This report details the work led by the Ohio Department of Natural Resources, Division of Geological Survey to evaluate the structure of nine formations in the Appalachian Basin of Ohio to assess seal integrity and whether carbon dioxide (CO 2 ) migration pathways could pose a potential risk. As part of Phase II MRCSP work, the Pennsylvania Geological Survey led the characterization of the Llandovery Silurian (“middle Silurian”) to Middle Devonian regional geologic characterization CO 2 utilization and storage. This preliminary assessment included generalized, regional-scale mapping, as well as the development of relatively coarse structure and isochore maps of stratigraphic intervals of interest for the Appalachian Basin. Because this stratigraphic interval offers several enhanced oil recovery targets and secondary targets for CO 2 storage, the Ohio Geological Survey undertook efforts in Phase III to characterize the Ohio’s portion of the Appalachian Basin more precisely and to resolve inconsistencies within Ohio’s formation tops. This study refined the subsurface correlations of the Silurian-through-Devonian stratigraphic section in Ohio. Higher-resolution isochore mapping of the CO 2 storage reservoirs and seals provides greater detail of potential high-porosity zones and potential structural influence on the interval.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Refinements to the Graves–Pitarka Kinematic Rupture Generator, Including a Dynamically Consistent Slip-Rate Function, Applied to the 2019 Mw 7.1 Ridgecrest Earthquake

The main objective of this study is to develop physics–based constraints on the spatiotemporal variation of the slip–rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth–dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip–rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface–rupturing, planar strike–slip fault that includes a weak (negative to low–stress–drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high–stress–drop patches designed to mirror the large–slip areas seen in various rupture model inversions of the event. The locations of the high–stress–drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip–rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth–dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth–dependent parametric variation of the slip–rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth–dependent variation in the shape of the slip–rate function improves the fit to the observed near–fault ground motions in the 0.5–3 s period range.

58 GEOSCIENCES↗
DOE OSTI2023

Tunable phase structure in NaNbO 3 ceramics by grain-size effect, electric field and heat treatment

Large polarization and strain change during antiferroelectric - ferroelectric phase transition under electric field is the foundation for realizing excellent electrical properties in antiferroelectric ceramics, therefore, the adjustment of antiferroelectricity and clarification of the corresponding mechanism is the foundation for controlling electrical properties. NaNbO3 is the most complex perovskite system showing multiple antiferroelectric phases in a wide temperature range, in which the antiferroelectricity shows obvious instability with changing external and internal conditions, namely the antiferroelectric phase can be adjusted by grain-size effect, electric field and heat treatment. According to the systematical study in terms of the Rietveld refinement of synchrotron XRD and Raman, NaNbO3 exhibits a ferrielectric P21ma structure at room temperature, the ferroelectric component of which increases with decreasing grain size. Two antiferroelectric tetragonal phases exist around Curie temperature TC before the entrance of antiferroelectric R phase zone, while an antiferroelectric monoclinic phase, which can be maintained to room temperature by annealing treatment, acts as the bridge for the depolarization of the poled NaNbO3 with ferroelectric Q phase. A detailed phase diagram mainly focused on the antiferroelectric phase zones of NaNbO3 is plotted, which gives a clear understanding about the polymorphic phase transitions under different conditions. Finally, the results concluded in this work would give a clear guidance for designing high-performance NaNbO3-based lead-free ceramics from the point of structure.

36 MATERIALS SCIENCE↗