Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power coefficient”

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 505 records · Page 28

Heat Transfer Characteristics of Particle Flow Through Additively Manufactured (SS 316L) Lattice Frame Material Based on Octet-Shape Topology

Abstract This paper presents experimentally obtained heat transfer characteristics of additively manufactured lattice frame material based on Octet-shaped unit cell. Binder jetting technology was used to print lattices in Stainless Steel 316L material. Lattice porosities ranging from 0.75 to 0.9 were investigated where thermal transport characteristics were obtained for void space occupied by air and particles. Particle diameters were varied from 266–966 microns. Effective thermal conductivity and averaged heat transfer coefficient was calculated through steady-state experiments. It was found that presence of lattice enhances the effective thermal conductivity by 2–4 times when compared to packed bed of particles alone. Furthermore, for gravity-assisted particle flow through lattice panel, significantly high convective heat transfer coefficients ranging from 200–400 W/m2K were obtained for the range of particle diameters tested. The superior thermal transport properties of Octet-shape based lattice frame material for particle flow through them makes it a very promising candidate for particle-to-supercritical carbon dioxide (sCO2) heat exchanger in concentrating solar power (CSP) application.

Aider, Youssef↗

Uncertainty Quantification and Sensitivity Analysis of Non-Nuclear Advanced Controls Testbed Reactor Mockup

The research presented in this report describes our progress in applying stochastic methods and uncertainty quantification, parametric study, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated with Griffin, with the goal of developing a reduced order (surrogate) model which can be rapidly sampled while perturbing multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums, which are controlled by a hybrid proportional, integral and derivative (PID) controller, actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time and $\beta$ from a transient Griffin simulation executed with the improved quasi-static method to provide the kinetic parameters) as inputs to functions which control the CD rotation angle. Using a number of stochastic method approaches, we developed a dual purpose training-surrogate model of the NTP system using polynomial regression. The trained model can be rapidly sampled while simultaneously perturbing various input parameters of the model, such as coefficients on the PID control, or temperature (directly affect the neutron cross section). The surrogate model delivers accurate results orders-of-magnitude faster (minutes, not days) than the base model. Once the base model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs and their effect on the output. For example, coefficients used in the PID control system may vary due to some physical interference, or there may be uncertainty in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters and operational boundaries can be determined. The goal of this work is to support development of an advanced control system to operate CDs in a functioning NTP system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Acidity of surface-infiltrated binary oxides as a sensitive descriptor of oxygen exchange kinetics in mixed conducting oxides

Improving the kinetics of O 2 reduction on oxide surfaces is critical in many energy and fuel conversion technologies. Here we show that the acidity scale for binary oxides is a powerful descriptor for tuning and predicting oxygen surface exchange kinetics on mixed conducting oxides. By infiltrating a selection of binary oxides from strongly basic (Li 2 O) to strongly acidic (SiO 2 ) onto the surface of Pr 0.1 Ce 0.9 O 2-δ samples, it was possible to vary the chemical surface exchange coefficient k chem by 6 orders of magnitude, with basic oxides such as Li 2 O increasing k chem by nearly 1,000 times, with surface concentrations as low as 50 ppm impacting k chem . Strikingly, although the pre-exponential of k chem scales linearly with the acidity of the infiltrated binary oxide, there is nearly no change in the activation energy. Here, the origin of these dramatic changes is proposed to arise from the systematic increase in electron concentration at the Pr 0.1 Ce 0.9 O 2-δ surface with the decreasing acidity of the infiltrated binary oxide.

30 DIRECT ENERGY CONVERSION↗

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

Learning an Algebriac Multrigrid Interpolation Operator Using a Modified GraphNet Architecture

This work, building on previous efforts, develops a suite of new graph neural network machine learning architectures that generate data-driven prolongators for use in Algebraic Multigrid (AMG). Algebraic Multigrid is a powerful and common technique for solving large, sparse linear systems. Its effectiveness is problem dependent and heavily depends on the choice of the prolongation operator, which interpolates the coarse mesh results onto a finer mesh. Previous work has used recent developments in graph neural networks to learn a prolongation operator from a given coefficient matrix. In this paper, we expand on previous work by exploring architectural enhancements of graph neural networks. A new method for generating a training set is developed which more closely aligns to the test set. Asymptotic error reduction factors are compared on a test suite of 3-dimensional Poisson problems with varying degrees of element stretching. Results show modest improvements in asymptotic error factor over both commonly chosen baselines and learning methods from previous work.

97 MATHEMATICS AND COMPUTING↗

Influence on Structural Loading of a Wave Energy Converter by Controlling Variable-Geometry Components and the Power Take-Off

Oceans are harsh environments and can impose significant loads on deployed structures. The deployment of wave energy converters (WECs) faces a design challenge with apparently contradictory goals. A WEC should be designed to maximize the energy absorbed while ensuring the operating wave condition does not exceed the failure limits of the device itself. Therefore, the loads endured by the support structure are a design constraint for the system. Adaptability to different sea states is, therefore, highly desirable. This work uses a WEC-Sim model of a variable-geometry oscillating wave energy converter (VGOSWEC) mounted on a support structure simulated under different wave scenarios. A VGOSWEC resembles a paddle pitching about a fixed hinge perpendicular to the incoming wave fronts. Therefore, the hinge experiences loads perpendicular to its axis as it maintains its position. The geometry of the VGOSWEC is varied by opening a series of controllable flaps on the pitching paddle when the structure experiences threshold loads. Because opening the flaps lets the waves transmit through the paddle, it is hypothesized that opening the flaps should result in load shedding at the base of the support structure. The load shedding is achieved by reducing the moments about the hinge axis. This work compares the hydrodynamic coefficients, natural periods, and response amplitude operators from completely closed to completely open configurations of the controllable flaps. The comparisons quantify the effects of letting the waves transmit through the VGOSWEC. This work shows that the completely open configuration can reduce the pitch and surge loads on the base of the support structure by as much as 80%. It was observed that at the paddle's resonance frequency, the loads on the structure increased substantially. This increase in loads can be mitigated by a rotational power take-off damping about the hinge axis. Changing the rotational power take-off damping was identified as an additional design parameter that can be used to control the loads experienced by the WEC's support structure.

fixed and floating structures↗

Hydrogen Diffusion in Slit Pores: Role of Temperature, Pressure, Confinement, and Roughness

Diffusion of hydrogen (H 2 ) is important to understand the leakage risk and transport behavior for H 2 geologic storage. We applied molecular dynamics simulations to investigate the influencing factors of H2 diffusion in the slit pores of calcite, hematite, and quartz, owing to their abundance. It is revealed that the H2 self-diffusion coefficient increases with the temperature, regardless of the type of pore minerals. The diffusion of H 2 in the 20 nm slit pores falls into the bulk diffusion regime when the pressure is 10 MPa. The self-diffusion of H 2 decreases with pressure in all three types of slit pores, following a power law model with the exponents ranging from -0.825 to -0.964. Furthermore, the impact of confinement on H 2 diffusion is more pronounced for the slit pores with stronger interactions with H 2 -like calcite. The role of surface roughness in H 2 diffusion depends on the slit aperture. The rough surface enhances H 2 diffusion in the larger slit pores due to the enlarged effective pore space, whereas it weakens H 2 diffusion in the small slit pores due to stronger adsorption. These findings will fill the knowledge gap on the coupling effect of different factors influencing H 2 diffusion.

08 HYDROGEN↗

High Thermoelectric Performance in Chalcopyrite Cu 1-x Ag x GaTe 2 –ZnTe: Nontrivial Band Structure and Dynamic Doping Effect

The understanding of thermoelectric properties of ternary I–III–VI 2 type (I = Cu, Ag; III = Ga, In; and VI = Te) chalcopyrites is less well developed. Although their thermal transport properties are relatively well studied, the relationship between the electronic band structure and charge transport properties of chalcopyrites has been rarely discussed. In this study, we reveal the unusual electronic band structure and the dynamic doping effect that could underpin the promising thermoelectric properties of Cu 1–x Ag x GaTe 2 compounds. Density functional theory (DFT) calculations and electronic transport measurements suggest that the Cu 1–x Ag x GaTe 2 compounds possess an unusual non-parabolic band structure, which is important for obtaining a high Seebeck coefficient. Moreover, a mid-gap impurity level was also observed in Cu 1–x Ag x GaTe 2 , which leads to a strong temperature-dependent carrier concentration and is able to regulate the carrier density at the optimized value for a wide temperature region and thus is beneficial to obtaining the high power factor and high average ZT of Cu 1–x Ag x GaTe 2 compounds. We also demonstrate a great improvement in the thermoelectric performance of Cu 1–x Ag x GaTe 2 by introducing Cu vacancies and ZnTe alloying. The Cu vacancies are effective in increasing the hole density and the electrical conductivity, while ZnTe alloying reduces the thermal conductivity. As a result, a maximum ZT of 1.43 at 850 K and a record-high average ZT of 0.81 for the Cu 0.68 Ag 0.3 GaTe 2 –0.5%ZnTe compound are achieved.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electron ITB formation in EAST high poloidal beta plasmas under dominant electron heating

Abstract Plasma confinement and transport in tokamaks play a crucial role in the development of high poloidal beta steady-state operation scenarios. Therefore, it is very important to study the relevant mechanism of internal transport barriers (ITBs), which can help plasmas to obtain better confinement in order to achieve higher fusion gain. This paper mainly introduces the analysis of the characteristics of electron heat transport of discharges with ITB in high β P operation regime in Experimental Advanced Superconducting Tokamak (EAST). Based on the statistical analysis of stable discharges with β P > 1.5, it is found that there is an obvious bifurcation of the normalised electron temperature gradient (ETG) ( R / L T e ) in the range of β P = 2–2.2. Then the discharges of lower β P ( β P < 2, where the value of β P is below the bifurcation threshold) and of higher β P ( β P > 2.2, where the value of β P is above the threshold) were selected for analysis. The diagnostic data provided by Thomson scattering, x-ray crystal spectrometry and charge exchange recombination spectroscopy are used to provide reliable parameter profiles and then combined with the data of external magnetic probe measurements and the polarisation interferometer diagnosis system to fully reconstruct the balance. A relevant plasma current calculation model is used to calculate and analyse the current density profiles and power deposition, and then the transport analysis is carried out. Interestingly, in the higher β P discharges, it is found that the turbulence intensity provided by the CO 2 laser collective scattering system gradually decreases and the normalised ETG R / L T e gradually increases. It is also found that the electron heat transport coefficient decreases in the discharge with higher β P and the growth rate of the electron-scale turbulence calculated by transport gyro-Landau fluid (TGLF) is significantly reduced. Meanwhile, similar conclusions are also obtained in the discharges when the β P is further increased.

Wang, Z. H. (ORCID:0009000203431036)↗

Blade and Rim Seal Design of a First Stage High Pressure Turbine for a 300 MWe Supercritical CO2 Power Cycle

A first stage high-pressure turbine (HPT) blade is optimized for a 300 MWe supercritical CO2 (sCO2) power cycle using the surrogate-assisted genetic algorithm optimizer in Numeca FINE/Design 3D with objectives of increasing efficiency and decreasing heat load to the blade. The National Institute of Standards and Technology Reference Fluid Thermodynamic and Transport Properties Database (NIST REFPROP) [1] data for supercritical CO2 is formatted into tables of bicubic polynomial coefficients for use in condensable gas simulations in FINE/Turbo. Nearly 3000 unique shapes are evaluated via three-dimensional Reynolds Averaged Navier Stokes simulations, yielding increases in efficiency of up to 0.85 percentage points and decreases in heat load of 14%. A final blade, deemed the advanced blade, is chosen for future experimental analysis. Following this, a squealer tip optimization is performed on both the baseline and advanced blade designs. This optimization resulted in a performance gain of 1.25 points in efficiency and 15% reduction in tip heat load compared to the baseline flat tip design at the same clearance. In tandem, an optimization of the rotor-stator platform rim seal is performed using a parametrized geometry allowing for straight, meandering, and knife seal cavities. This multi-objective optimization focuses on decreasing the cooling mass flow and increasing the heat flux from the rotor and stator disk. The optimization resulted in cooling mass flow decreases of up to 26% while maintaining the average heat flux on the rim seal.

Tuite, Logan↗

Influence on Structural Loading of a Wave Energy Converter by Controlling Variable-Geometry Components and the Power Take-Off: Preprint

Oceans are a harsh environment and can impose significant loads on deployed structures. Deployment of wave energy converters (WECs) faces a design challenge with apparently contradictory goals. A WEC should be designed to maximize the energy absorbed while ensuring the operating wave condition does not exceed the failure limits of the device itself. Therefore, the loads endured by the support structure are a design constraint for the system. Adaptability to different sea states is, therefore, highly desirable. This work uses a WEC-Sim model of a Variable Geometry Oscillating Wave Energy Converter (VGOSWEC) mounted on a support structure simulated under different wave scenarios. A VGOSWEC resembles a paddle pitching about a fixed hinge perpendicular to the incoming waves-fronts. Therefore, the hinge experiences loads perpendicular to its axis as it maintains its position. The geometry of the VGOSWEC was varied by opening a series of controllable flaps on the pitching paddle when the structure experiences threshold loads. Since opening the flaps lets the waves transmit through the paddle, it is hypothesized that opening the flaps should result in load shedding at the base of the support structure. This load shedding is achieved by reducing the moments about the hinge axis. This work compares the hydrodynamic coefficients, natural periods, and Response Amplitude Operators (RAOs) from completely closed to completely open configurations of the controllable flaps. These comparisons quantify the effects of letting the waves transmit through the VGOSWEC. This work shows that the completely open configuration can reduce the pitch and surge loads on the base of the support structure by as much as 80%. It was observed that at the paddle's resonance frequency, the loads on the structure increased substantially. This increase in loads can be mitigated by a rotational Power Take-off (PTO) damping about the hinge axis. Changing the rotational PTO damping was identified as an additional design parameter that can be used to control the loads experienced by the WEC's support structure.

fixed and floating structures↗

Exploring mounting solutions for cryogenically cooled thin crystal optics in high power density x-ray free electron lasers

This study investigates three mounting methods—clamping, soldering, and a hybrid clamping–soldering approach—for cryogenically cooled thin diamond crystals crucial to stable operation of X-ray Free Electron Laser (XFEL) systems. While clamping methods exhibit temperature resilience and flexibility, meticulous design is required to prevent stress-induced warping and reduce thermal contact area. Soldering methods offer reliable mechanical and thermal bonding but encounter challenges due to the coefficient of thermal expansion mismatch at cryogenic temperatures. The hybrid method, integrating clamping and soldering with strain relief cuts, effectively mitigates overall distortion caused by mounting and XFEL thermal loads. These findings offer a novel mounting solution for high-performance x-ray optics in XFEL research and applications, ensuring stability and optimal functionality in cryogenic conditions.

Instruments & Instrumentation↗

Effects of Particle Size and Concentration of Magnesium Oxide on the Lubricating Performance of a Chloride Molten Salt for Concentrating Solar Power

One engineering challenge for the next-generation concentrating solar power (CSP) is to identify appropriate heat transfer and thermal energy storage media. A chloride salt mixture (20 wt % NaCl + 40 wt % MgCl 2 + 40 wt % KCl) in a molten phase is proposed as a potential fluid for carrying the heat as well as lubricating the CSP pump bearings. Due to the hygroscopic nature of MgCl 2 , the presence of magnesium oxide (MgO) particles is inevitable. In order to understand the impact of MgO particles, this study investigated the tribological behavior of the molten chloride salt mixture containing MgO of various particle sizes and concentrations. Tests were conducted at 750 °C in a dry Ar environment using a zirconia ball sliding against a Haynes 244 alloy disc lubricated by the molten salt. There was a clear trend of the increased friction coefficient and wear loss for a larger particle size and/or a higher concentration of MgO. Distinct wear mechanisms were identified for the ceramic ball and alloy disc. Results provide fundamental insights for this new research avenue of molten salt lubrication for CSP.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Protein conformational switch discerned via network centrality properties

Network analysis has emerged as a powerful tool for examining structural biology systems. The spatial organization of the components of a biomolecular structure has been rendered as a graph representation and analyses have been performed to deduce the biophysical and mechanistic properties of these components. For proteins, the analysis of protein structure networks (PSNs), especially via network centrality measurements and cluster coefficients, has led to identifying amino acid residues that play key functional roles and classifying amino acid residues in general. Whether these network properties examined in various studies are sensitive to subtle (yet biologically significant) conformational changes remained to be addressed. Here, we focused on four types of network centrality properties (betweenness, closeness, degree, and eigenvector centralities) for conformational changes upon ligand binding of a sensor protein (constitutive androstane receptor) and an allosteric enzyme (ribonucleotide reductase). We found that eigenvector centrality is sensitive and can distinguish salient structural features between protein conformational states while other centrality measures, especially closeness centrality, are less sensitive and rather generic with respect to the structural specificity. We also demonstrated that an ensembleinformed, modified PSN with static edges removed (which we term PSN*) has enhanced sensitivity at discerning structural changes.

59 BASIC BIOLOGICAL SCIENCES↗

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗

FWP FEAA149: “Next Generation Environmental Barrier Coatings”

Environmental barrier coatings (EBCs) are required coatings for utilization of SiC/SiC ceramic matrix composite (CMC) components in gas turbines, where the EBC represents the life-limiting factor for such components. EBC/CMC systems have shown success in aero-engine applications with increased turbine inlet temperatures and improved efficiencies, which are achieved through higher temperature stability, lower density, and decreased reliance on cooling air compared to traditional superalloys. While industrial gas turbines (IGTs) do not currently utilize SiC/SiC CMCs, the current shift towards low-carbon or carbon-free fuel sources for power generation could result in a need for EBC/CMC components with higher temperature capabilities. In this work, three tasks were outlined to improve understanding of EBC lifetimes to encourage use in IGTs with carbon-free fuel such as hydrogen: 1. Define the bond coating oxidation kinetics and EBC failure criteria, 2. Measure thermal expansion coefficients of each layered material, and 3. Perform advanced characterization and modeling to assess EBC lifetimes. Cyclic steam oxidation tests were conducted on various EBC/Si/SiC chemistries and EBC/SiC architectures to define substrate oxidation kinetics and EBC failure modes. An open-source code was developed to quantify the undulating thermally grown oxide thickness with thousands of measurements from specimen cross-section images. Bond coating oxidation kinetics were determined and used to develop a kinetic and thermodynamic model for predicting EBC lifetimes. High-temperature Raman spectroscopy was utilized for determining the SiO 2 thermally grown oxide phase transformation as the life-limiting feature for EBCs. Model efforts supported the claim that the SiO 2 phase transformation causes elevated stress during thermal cycling with associated cracking that decreases the adhesion strength of the EBC, eventually resulting in coating spallation. The finite element model subroutine will be made publicly available upon internal review. Further development of an EBC lifetime model for IGTs involves definition of a critical SiO 2 thickness for EBC spallation and must also consider both environmental (gas velocity, pressure, etc.) and specimen (EBC dopants, layer architectures, etc.) effects into predicted bond coating oxidation kinetics for long-lifetime components.

36 MATERIALS SCIENCE↗

Performance of reverse osmosis membrane with large feed pressure fluctuations from a wave-driven desalination system

Wave-driven desalination systems are proposed water treatment systems that involve reverse osmosis of seawater powered directly by wave motion. Such a configuration would result in drastic feed pressure fluctuations. For a technology conventionally operated with a constant feed condition, the effect of these variable pressures on membrane integrity and performance is unknown. Here, experiments were conducted with spiral wound membranes coupled to a system capable of producing feed pressure fluctuations of more than 400 psi. Feed composition included 5, 20, and 35 g/L NaCl, and a synthetic seawater at normal and 1.5x concentration. The variable feed conditions included sine-like pressure waves swings of 200-500 and 500-900 psi with frequencies of 1.25, 7.5, and 12 waves/min, and a model-generated random waveform. Between each wave experiment we performed membrane integrity tests at 650 psi and 25 g/L NaCl feed, which showed a 7.4% drop in the membrane's water permeability coefficient, an 18.4% flux decline, and more than 99% salt rejection over 1770 h of cumulative experimental time. Analysis of permeate samples showed high salt rejection. In general, variable feed pressure had no significant deleterious effect on membrane integrity or performance.

16 TIDAL AND WAVE POWER↗

Defect-mediated diffusion of implanted Mg in GaN: Suppressing dopant redistribution by sequential thermal and microwave annealing

The diffusion behavior of Mg in Mg/N co-implanted GaN is investigated in response to a set of annealing conditions and methodologies, namely, 1000 °C/30 min thermal anneal, by high-temperature pulsed gyrotron microwave annealing at 1420 or 1500 °C, or by thermal and microwave annealing, sequentially. After 1000 °C annealing, the diffusion of Mg in GaN is found to be negligible, as measured by secondary ion mass spectrometry. Annealing by gyrotron microwave annealing alone induces the diffusion of Mg at a rate on the order of 10 −12 cm 2 /s. However, the use of a thermal anneal before microwave gyrotron annealing reduces this rate by an order of magnitude to 10 −13 cm 2 /s. We find that a model that considers Mg diffusion from an inhomogeneous medium that contains a defect-rich implanted region near-surface to a relatively pristine region below the implant range better explains the observed diffusion behavior than a conventional model that assumes a homogeneous medium. By analyzing the diffusion behavior using the Boltzmann–Matano method, we present a discussion of reduction in [V Ga ] by thermal annealing at 1000 °C, leading to a suppressed diffusion coefficient during subsequent high-temperature annealing relative to diffusion after 1420/1500 °C annealing alone. This effect holds potential for improvement in the precision of selectively doped regions for future applications based on the (Al)GaN material system. An improved field profile control in real devices can increase the breakdown and current-handling capabilities in power electronic applications.

Meyers, V.↗