Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion current”

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

A new method for operating a continuous-flow diffusion chamber to investigate immersion freezing: assessment and performance study

Abstract. Glaciation in mixed-phase clouds predominantly occurs through the immersion-freezing mode where ice-nucleating particles (INPs) immersed within supercooled droplets induce the nucleation of ice. Model representations of this process currently are a large source of uncertainty in simulating cloud radiative properties, so to constrain these estimates, continuous-flow diffusion chamber (CFDC)-style INP devices are commonly used to assess the immersion-freezing efficiencies of INPs. This study explored a new approach to operating such an ice chamber that provides maximum activation of particles without droplet breakthrough and correction factor ambiguity to obtain high-quality INP measurements in a manner that previously had not been demonstrated to be possible. The conditioning section of the chamber was maintained at −20 ∘C and water relative humidity (RHw) conditions of 113 % to maximize the droplet activation, and the droplets were supercooled with an independently temperature-controlled nucleation section at a steady cooling rate (0.5 ∘C min−1) to induce the freezing of droplets and evaporation of unfrozen droplets. The performance of the modified compact ice chamber (MCIC) was evaluated using four INP species: K-feldspar, illite-NX, Argentinian soil dust, and airborne soil dusts from an arable region that had shown ice nucleation over a wide span of supercooled temperatures. Dry-dispersed and size-selected K-feldspar particles were generated in the laboratory. Illite-NX and soil dust particles were sampled during the second phase of the Fifth International Ice Nucleation Workshop (FIN-02) campaign, and airborne soil dust particles were sampled from an ambient aerosol inlet. The measured ice nucleation efficiencies of model aerosols that had a surface active site density (ns) metric were higher but mostly agreed within 1 order of magnitude compared to results reported in the literature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A new method for operating a continuous flow diffusion chamber to investigate immersion freezing: assessment and performance study

Glaciation in mixed-phase clouds predominately occurs through the immersion freezing mode where ice nucleating particles (INPs) immersed within supercooled droplets induce nucleation of ice. Currently, model representations of this process are a large source of uncertainty in simulating cloud radiative properties, and to constrain these estimates, continuous flow diffusion chamber (CFDC)-style INP devices are commonly used to assess the immersion freezing efficiencies of INPs. In this study, a new approach was explored to operating such an ice chamber that provides maximum activation of particles without droplet breakthrough and correction factor ambiguity to obtain high-quality INP measurements in a manner that has not been demonstrated as possible previously. The conditioning section of the chamber was maintained at ~ –20 °C and water relative humidity (RH w ) ~ 113 % conditions to maximize the droplet activation, and the droplets were supercooled with an independently temperature-controlled nucleation section at a steady cooling rate (0.5 °C min –1 ) to induce the freezing of droplets and evaporation of unfrozen droplets. The performance of the modified ice chamber was evaluated using four INP species: K-feldspar, illite-NX, Argentinian soil dust, and airborne arable dust that had shown ice nucleation over a wide span of supercooled temperatures. Dry dispersed and size-selected K-feldspar particles were generated in the laboratory. Illite-NX and soil dust particles were sampled during the second phase of the Fifth International Ice Nucleation Workshop (FIN-02) campaign, and airborne arable dust particles were sampled from the aerosol inlet located on the rooftop of the laboratory. The measured ice nucleation efficiencies of model aerosols with a surface active site density (n s ) metric were higher, but mostly agreed within one order of magnitude compared to literature results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Design Rules of the Mixing Phase and Impacts on Device Performance in High-Efficiency Organic Photovoltaics

In nonfullerene acceptor- (NFA-) based solar cells, the exciton splitting takes place at both domain interface and donor/acceptor mixture, which brings in the state of mixing phase into focus. The energetics and morphology are key parameters dictating the charge generation, diffusion, and recombination. It is revealed that tailoringthe electronic properties of the mixing region by doping with larger-bandgap components could reduce the density of state but elevate the filling state level, leading to improved open-circuit voltage ( V OC ) and reduced recombination. The monomolecular and bimolecular recombinations are shown to be intercorrelated, which show a Gaussian-like relationship with V OC and linear relationship with short-circuit current density ( J SC ) and fill factor (FF). The kinetics of hole transfer and exciton diffusion scale with J SC similarly, indicating the carrier generation in mixing region and crystalline domain are equally important. From the morphology perspective, the crystalline order could contribute to V OC improvement, and the fibrillar structure strongly affects the FF. These observations highlight the importance of the mixing region and its connection with crystalline domains and point out the design rules to optimize the mixing phase structure, which is an effective approach to further improve device performance.

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

Diffuse Reflectance Spectroscopy and Principal Component Analysis to Retrospectively Determine Production History of Plutonium Dioxide

Diffuse reflectance spectroscopy measurements in the shortwave infrared (930–1600 nm) spectral region were acquired for Pu 2 (C 2 O 4 ) 3 •9H 2 O and its thermal decomposition product, PuO 2 . We analyzed a total of eight PuO 2 samples that were produced at different calcination temperatures (300, 350, 450, 525, 600, 675, 750, and 900 °C). Our goal was to identify spectroscopic fingerprints that could be used to gain retrospective information regarding the production parameters of these important nuclear compounds. The diffuse reflectance spectrum of Pu 2 (C 2 O 4 ) 3 •9H 2 O features several broad bands that currently preclude detailed analysis. However, all PuO 2 samples produced relatively sharp spectral features that got sharper and more intense for samples that were produced at higher calcination temperatures. The electronic band observed at 1433 nm in the diffuse reflectance spectra of PuO 2 was found to be a sensitive indicator of crystallinity; a result that is corroborated by ancillary Raman spectroscopy measurements. Principal component analysis of diffuse reflectance spectra was able to clearly rank and categorize PuO 2 samples based on the calcination temperature that was employed during their production. Thus, we show herein that important retrospective information pertaining to the process history of PuO 2 can be gained through the relatively simplistic combination of diffuse reflectance spectroscopy and principal component analysis. This discovery presents a new method for determining the provenance and process history of PuO 2 and should have an impact in the fields of nuclear forensics and nuclear nonproliferation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Two-temperature effects in Hall-MHD simulations of the HIT-SI experiment

A two-temperature Hall-magnetohydrodynamic (Hall-MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet 3D MHD code and used to model plasma dynamics in the Helicity Injected Torus–Steady Inductive (HIT-SI) experiment. The two-temperature model is utilized for HIT-SI simulations in both the PSI-Tet and NIMROD codes at a number of different injector frequencies in the 14.5–68.5 kHz range. At all frequencies, the NIMROD two-temperature model results in increased toroidal current, lower chord-averaged density, higher average temperatures, outward radial shift of the current centroid, and axial symmetrization of the current centroid, relative to the single-temperature NIMROD simulations. The two-temperature PSI-Tet model illustrates similar trends, but at high frequency operation, it exhibits lower electron temperature, smaller toroidal current, and decreased axial symmetrization with respect to the single-temperature PSI-Tet model. With all models, average temperatures and toroidal currents increase with the injector frequency. Power balance and heat fluxes to the wall are calculated for the two-temperature PSI-Tet model and illustrate considerable viscous and compressive heating, particularly at high injector frequency. Parameter scans are also presented for artificial diffusivity, wall temperature, and density. Both artificial diffusivity and the density boundary condition significantly modify the plasma density profiles, leading to larger average temperatures, toroidal current, and relative density fluctuations at low densities. As a result, a low density simulation achieves sufficiently high current gain (G > 5) to generate significant volumes of closed flux lasting 1–2 injector periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Pulsed power accelerator surface Joule heating models

Understanding the effects of contaminant plasmas generated within the Z machine at Sandia is critical to understanding current loss mechanisms. The plasmas are generated at the accelerator electrode surfaces and include desorbed species found in the surface and substrate of the walls. These desorbed species can become ionized. The timing and location of contaminant species desorbed from the wall surface depend non-linearly on the local surface temperature. For accurate modeling, it is necessary to utilize wall heating models to estimate the amount and timing of material desorption. One of these heating mechanisms is Joule heating. In this paper, we propose several extended semi-analytic magnetic diffusion heating models for computing surface Joule heating and demonstrate their effects for several representative current histories. We quantitatively assess under what circumstances these extensions to classical formulas may provide a validatable improvement to the understanding of contaminant desorption timing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An effective description of charge diffusion and energy transport in a charged plasma from holography

We discuss the physics of sound propagation and charge diffusion in a plasma with non-vanishing charge density. Our analysis culminates the program initiated in to construct an open effective field theory of low-lying modes of the stress tensor and charge current in such plasmas. We model the plasma holographically as a Reissner-Nordström-AdS d+1 black hole, and study linearized fluctuations of longitudinally polarized scalar gravitons and photons in this background. We demonstrate that the perturbations can be decoupled and repackaged into the dynamics of two designer scalars, whose gravitational coupling is modulated by a non-trivial dilatonic factor. The holographic analysis allows us to isolate the phonon mode from the charge diffusion mode, and identify the combination of currents that corresponds to each of them. We use these results to obtain the real-time Gaussian effective action, which includes both the retarded response and the associated stochastic (Hawking) fluctuations, accurate to quartic order in gradients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

High-temperature compressive creep tests of U 3 Si 2 with spark plasma sintering: Experiments and finite element modeling

This paper reports a systematic high-temperature creep behavior of dense U 3 Si 2 pellets using a spark plasma sintering (SPS) apparatus at elevated temperatures and pressures under vacuum conditions. The stress exponent was subsequently derived to be 3.21 at 1173 K and 2.17 at 1223 K, respectively, indicating a grain boundary sliding creep mechanism. The creep activation energy was determined to be 203.6 ± 19.0 kJ/mol, which agrees well with the literature. Finite element modeling was performed using the creep parameters fitted from the strain-time plot. Additionally, the results suggest an excellent match with the experimental data, confirming the validity of the experiments. Microstructure characterizations indicate that the main phase of the specimens after creep tests remains to be U 3 Si 2 , with a 4 μm thick layer of nano-sized particles induced from the diffusion between U 3 Si 2 and alumina disc used to avoid electric current passing through the sample. The successful conduct of creep experiments demonstrates the great potential of SPS to perform high-temperature mechanical testing of nuclear fuels under vacuum conditions. The subsequent finite element modeling exhibits excellent capabilities for accurately predicting material performance in the creep tests and provides a practical tool in evaluating nuclear fuels’ performance for a much-extended time scale.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structure and dynamical properties of two-dimensional dusty plasmas on one-dimensional periodic substrates

In the presence of the substrate, various structural and dynamical properties of two-dimensional dusty plasma (2DDP) were investigated using Langevin dynamical simulations. This paper reviews a series of results of the structural and dynamical properties of 2DDP modified by one-dimensional periodic substrates (1DPSs) as follows. First, when the depth of the 1DPS increases gradually, it is found that the static structural order of 2DDP along each potential well of the 1DPS increases first, and then decreases gradually. When the width of the 1DPS increases gradually, the particle arrangement in each potential well of the 1DPS changes from one straight row gradually to two rows, including the stable zigzag structure. Second, when there are two rows in each potential well of the 1DPS, the phonon spectra are split into two branches, corresponding to the breathing and sloshing modes, respectively. In addition, due to the pure repulsion between dust particles, from the obtained phonon spectra, the sloshing wave propagates backward at small wave numbers. Third, the calculated mean-square displacement (MSD) shown that, at the intermediate timescale between the initial ballistic and final diffusive motions, the particle motion tends to be more subdiffusive while the depth of the 1DPS increases. While the width of the 1DPS increases gradually, the long-time diffusive motion first increases, then decreases, and finally increases again, exhibiting the oscillation-like diffusion, due to the stable zigzag structure. Finally, when an external direct-current driving force is applied on all particles of 2DDP, three different depinning phases are discovered, which are the pinned, disordered plastic flow, and moving ordered states, respectively, as the driving force increases from zero. In addition, the continuous/discontinuous property of the phase transition between these different depinning states is investigated, showing that the transition property is modulated by the depth of the applied 1DPS

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Inertial and anisotropic pressure effects on cross-field electron transport in low-temperature magnetized plasmas

In this paper, a one-dimensional (1D) particle-in-cell Monte Carlo collision (PIC-MCC) model is developed to investigate the effects of anisotropic pressure and inertial terms due to non-Maxwellian velocity distribution functions on cross-field electron transport. The conservation of momentum is evaluated by taking the moments of the first-principles gas-kinetic equation. A steady-state discharge is obtained without any low-frequency ionization oscillations by considering an anomalous electron scattering profile. Further, the results obtained from the 1D PIC-MCC model are compared with fluid models, including the quasi-neutral drift-diffusion (DD), non-neutral DD, and full fluid moment models. The discharge current obtained from the PIC-MCC model is in good agreement with the fluid models. The cross-field electron transport due to the inertial terms, i.e. the gradient of axial and azimuthal drift, is evaluated. Moreover, PIC-MCC simulation results show non-zero, anisotropic, off-diagonal pressure tensor terms due to asymmetric non-Maxwellian electron velocity distribution function, potentially contributing to cross-field electron transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Comparing Cetane Number Measurement Methods

Cetane number is one of the most important fuel performance metrics for mixing controlled compression-ignition “diesel” engines, quantifying a fuel’s propensity for autoignition when injected into end-of-compression-type temperature and pressure conditions. The historical default and referee method on a Cooperative Fuel Research (CFR) engine configured with indirect fuel injection and variable compression ratio is cetane number (CN) rating. A subject fuel is evaluated against primary reference fuel blends, with heptamethylnonane defining a low-reactivity endpoint of CN = 15 and hexadecane defining a high-reactivity endpoint of CN = 100. While the CN scale covers the range from zero (0) to 100, typical testing is in the range of 30 to 65 CN. Alternatively, several constant-volume combustion chamber (CVCC)-based cetane rating devices have been developed to rate fuels with an equivalent derived cetane number (DCN) or indicated cetane number (ICN). These devices measure ignition delay for fuel injected into a fixed volume of high-temperature and high-pressure air to simulate end-of-compression-type conditions. In this study, a range of novel fuel compounds are evaluated across three CVCC methods: the Ignition Quality Tester (IQT), Fuel Ignition Tester (FIT), and Advanced Fuel Ignition Delay Analyzer (AFIDA). Resulting DCNs and ICNs are compared for fuels within the normal diesel fuel range of reactivity, as well as very high (~100) and very low DCNs/ICNs (~5). Distinct differences between results from various devices are discussed. This is important to consider because some new, high-efficiency advanced compression-ignition (CI) engine combustion strategies operate with more kinetically controlled distributed combustion as opposed to mixing controlled diffusion flames. These advanced combustion strategies may benefit from new fuel chemistries, but current rating methods of CN, DCN, and ICN may not fully describe their performance. In addition, recent evidence suggests ignition delay in modern on-road diesel engines with high-pressure common rail fuel injection systems may no longer directly correlate to traditional CN fuel ratings. Simulated end-of-compression conditions are compared for CN, DCN, and ICN and discussed in the context of modern diesel engines to provide additional insight. Results highlight the potential need for revised and/or multiple fuel test conditions to measure fuel performance for advanced CI strategies.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Mathematical Model for Li-S Cell with Shuttling-Induced Capacity Loss Approximation

Lithium sulfur (Li-S) batteries have the potential to outperform the current lithium ion batteries and transform the technology of the future. However, dissolution, diffusion, and shuttling of the dissolved polysulfides result in parasitic reactions and substantial capacity loss. To provide a better understanding of the shuttling process, a 1D porous electrode mathematical model has been developed in this paper. An approximation method is used to account for the shuttling-induced capacity loss by adding an extra source/sink term in the material balance equations for the species involved in the parasitic reactions. Shuttling constants used in the source terms can be determined by fitting the model predictions to the experimental measurements. Here, the results showed that by including the approximation method, the model was able to predict the active material loss and the continuous decrease of volume fractions of Li 2 S on the cathode surface. The model sheds light on the capacity loss mechanism occurring inside the cell as a result of the shuttling of polysulfides.

25 ENERGY STORAGE↗
DOE OSTI2022

Fusion RF Modeling Machine Learning (FusionML_RF) v1.0

FusionML_RF consists of multiple codes and trained machine learning (ML) models that perform low-cost output modeling from the Genray-CQL3D. Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. For example, completing a single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time. On the other hand, these ML models achieve ~ms of inference time with high accuracy across the input parameter space. This software collection consists of multiple components. (1) codes that use ML methods and precomputed Genray-CQL3D simulation output to build regression models that enable approximate computations of Genray-CLQ3D outputs from arbitrary but physically meaningful input parameters (surrogate modeling); (2) three trained models created by the team, using a database of 16,000+ GENRAY/CQL3D simulations, to study the performance of ML models for surrogate modeling; (3) codes that load the trained models and simulation data, and then compute mean squared error between the models' predictions and the ground truth of simulation output data. This collection is being made available in conjunction with a scientific publication about the work to promote reusability and provide an artifact of the scientific work.

Bai, Zhe↗
DOE OSTI2020

UNR ETI Computational Investigation (Year 2)

The LANL ASC Arbitrary Lagrange/Eulerian (ALE) code, FLAG, is being used to model experiments conducted under this grant in order to understand the physics limiting high-current conduction and thermal plasma production. This includes the effects of nonlinear magnetic diffusion (NLMD), the magnetic Rayleigh-Taylor (MRT) instability, and the electro-thermal instability (ETI). FLAG is a multi-dimensional, unstructured-grid, multi-material, multi-temperature, shockphysics code with a single-fluid, resistive magnetohydrodynamics (MHD) package. Equations-ofstate (EOS), conductivity, material-strength, and melt can be treated with a variety of analytic or tabular models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Pre-Symptomatic COVID Screening

Temperature checks for fever are extensively used for preliminary COVID screenings but are ineffective during the incubation stage of infection when a person is asymptotic. Researchers at the European Centre for Disease Prevention and Control concluded that approximately 75% of passengers infected with COVID-19 and traveling from affected Chinese cities would not be detected by early screening. Core body temperature is normally kept within a narrow range and has the smallest relative standard deviation of all vital signs. Heat in the body is prioritized around internal organs at the expense of the periphery by controlling blood flow. In fact, blood flow to the skin may vary by a factor of 100 depending on thermal conditions. This adaptation causes rapid temperature fluctuations in different skin regions from changes in cardiac output, metabolism, and likely cytokine diffusion during inflammation that would not be seen in average core body temperature. Current IR and thermal scanners used for temperature checks are not necessarily reflective of core body temperatures and require cautious interpretation as they frequently result in false positive and false negative diagnosis. Hand held thermometers measure average skin temperatures and can get readings that differ from core body temperature by as much as 7°. Rather than focusing on a core body temperature threshold assessment we believe that variability of temperature patterns using a novel wearable transdermal microneedle sensor will be more sensitive to infections in the incubation stage and propose to develop a wearable transdermal temperature sensor using established Sandia microneedle technology for pre-symptomatic COVID screening that can additionally be used to monitor disease progression at later stages.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Pre Symptomatic COVID Screening

Temperature checks for fever are extensively used for preliminary COVID screenings but are ineffective during the incubation stage of infection when a person is asymptotic. Researchers at the European Centre for Disease Prevention and Control concluded that approximately 75% of passengers infected with COVID-19 and traveling from affected Chinese cities would not be detected by early screening. Core body temperature is normally kept within a narrow range and has the smallest relative standard deviation of all vital signs. Heat in the body is prioritized around internal organs at the expense of the periphery by controlling blood flow. In fact, blood flow to the skin may vary by a factor of 100 depending on thermal conditions. This adaptation causes rapid temperature fluctuations in different skin regions from changes in cardiac output, metabolism, and likely cytokine diffusion during inflammation that would not be seen in average core body temperature. Current IR and thermal scanners used for temperature checks are not necessarily reflective of core body temperatures and require cautious interpretation as they frequently result in false positive and false negative diagnosis. Hand held thermometers measure average skin temperatures and can get readings that differ from core body temperature by as much as 7°. Rather than focusing on a core body temperature threshold assessment we believe that variability of temperature patterns using a novel wearable transdermal microneedle sensor will be more sensitive to infections in the incubation stage and propose to develop a wearable transdermal temperature sensor using established Sandia microneedle technology for pre symptomatic COVID screening that can additionally be used to monitor disease progression at later stages.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Final-stage densification kinetics of direct current–sinteredZrB2

Final-stage sintering was analyzed for nominally phase pure zirconium diboride synthesized by borothermal reduction of high-purity ZrO2. Analysis was conducted on ZrB2 ceramics with relative densities greater than 90% using the Nabarro–Herring stress–directed vacancy diffusion model. Temperatures of 1900°C or above and an applied uniaxial pressure of 50 MPa were required to fully densify ZrB2 ceramics by direct current sintering. Ram travel data were collected and used to determine the relative density of the specimens during sintering. Specimens sintered between 1900 and 2100°C achieved relative densities greater than 97%, whereas specimens sintered below 1900°C failed to reach the final stage of sintering. The average grain size ranged from 1.0 to 14.7 μm. The activation energy was calculated from the slope of an Arrhenius plot that used the Kalish equation. The activation energy was 162 ± 34 kJ/mol, which is consistent with the activation energy for dislocation movement in ZrB2. The diffusion coefficients for dislocation motion that controls densification were 5.1 × 10−6 cm2/s at 1900°C and 5.1 × 10−5 cm2/s at 2100°C, as calculated from activation energy and average grain sizes. This study provides evidence that the dominant mechanism for final-stage sintering of ZrB2 ceramics is dislocation motion.

36 MATERIALS SCIENCE↗