Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Debye temperature”

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.

80 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The imaginary part of the heavy-quark potential from real-time Yang-Mills dynamics

We extract the imaginary part of the heavy-quark potential using classical-statistical simulations of real-time Yang-Mills dynamics in classical thermal equilibrium. The r-dependence of the imaginary part of the potential is extracted by measuring the temporal decay of Wilson loops of spatial length r. We compare our results to continuum expressions obtained using hard thermal loop theory and to semi-analytic lattice perturbation theory calculations using the hard classical loop formalism. We find that, when plotted as a function of m D r, where m D is the hard classical loop Debye mass, the imaginary part of the heavy-quark potential shows little sensitivity to the lattice spacing at small m D r ≲ 1 and agrees well with the semi-analytic hard classical loop result. For large quark-antiquark separations, we quantify the magnitude of the non-perturbative long-range corrections to the imaginary part of the heavy-quark potential. We present our results for a wide range of temperatures, lattice spacings, and lattice volumes. This work sets the stage for extracting the imaginary part of the heavy-quark potential in an expanding non-equilibrium Yang Mills plasma.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Non-thermal atmospheric pressure plasma–liquid synthesis of organic acids in aqueous solution from carbon monoxide

This work aims at understanding the conversion of CO to organic acids, namely oxalic acid and formic acid, using non-thermal atmospheric pressure plasma over aqueous solutions. CO exhibited significantly higher conversion to organic acids (more than 15×) compared to CO 2 under the same reaction conditions. The result bolsters a proposed two-step process for CO 2 fixation, whereby CO 2 is first converted to CO, and then CO is converted to organic acids. The organic acids produced from CO are intermediates in the water–gas shift (WGS) reaction of CO in the presence of an aqueous solution to dissolved CO 2 and hydrogen gas. Based on a simple thermodynamic analysis, the organic acid yield was increased by lowering the plasma–liquid reaction temperature using an ice bath to cool the reaction flask. The composition of the organic acids could be varied by changing the pH of the solution. Oxalate was formed in higher concentrations with increasing solution pH above the pK a of the radical species (CO 2 )˙ − . Below the pK a value, formate was the exclusive organic acid formed. The production of formate has a rather weak pH dependence but is enhanced slightly at a basic pH above 10. Furthermore, at basic pH, the effect of electrolyte concentration comes into play. Higher electrolyte concentrations, leading to shorter electrolyte Debye lengths, resulted in lowered organic acid yields. The highest yields of organic acids obtained in our system were 122 mg L −1 for oxalate and 77 mg L −1 for formate at an optimum 1 mM NaOH concentration in the starting solution. This work is a successful pioneering example of CO to organic acids conversion using non-thermal plasmas, which opens the pathway for a promising two-step conversion process of CO 2 to organic acids.

Sudagar, Alcina Johnson [Washington University in ↗
DOE OSTI2023

Fast correlation heating in moderately coupled electron–ion plasmas

If the electrons in a plasma are suddenly heated, the resulting change in Debye shielding causes the ion kinetic energy to quickly increase. For the first time, this correlation heating , which is much faster than collisional energy exchange, is rigorously derived for a moderately coupled, electron–ion plasma. The electron–ion mass ratio is taken to be the smallest parameter in the Bogoliubov–Born–Green–Kirkwood–Yvon hierarchy, smaller even than the reciprocal of the plasma parameter. This ordering differs from conventional kinetic theory by making the electron collision rates faster than the ion plasma frequency, which allows stronger coupling and makes the ion heating a function only of the total energy supplied to the electrons. The calculation uses known formulae for correlations in a two-temperature plasma, for which a new, elementary derivation is presented. Suprathermal ions may be created more rapidly by this mechanism than by ion–electron Coulomb collisions. This means that the use of a femtosecond laser pulse could potentially help to achieve ignition in certain fast ignition approaches to inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Fabrication of single color centers in sub-50 nm nanodiamonds using ion implantation

Diamond color centers have been widely studied in the field of quantum optics. The negatively charged silicon vacancy (SiV – ) center exhibits a narrow emission linewidth at the wavelength of 738 nm, a high Debye–Waller factor, and unique spin properties, making it a promising emitter for quantum information technologies, biological imaging, and sensing. In particular, nanodiamond (ND)-based SiV – centers can be heterogeneously integrated with plasmonic and photonic nanostructures and serve as in vivo biomarkers and intracellular thermometers. Out of all methods to produce NDs with SiV – centers, ion implantation offers the unique potential to create controllable numbers of color centers in preselected individual NDs. However, the formation of single color centers in NDs with this technique has not been realized. We report the creation of single SiV– centers featuring stable high-purity single-photon emission through Si implantation into NDs with an average size of ~20 nm. We observe room temperature emission, with zero-phonon line wavelengths in the range of 730–800 nm and linewidths below 10 nm. Our results offer new opportunities for the controlled production of group-IV diamond color centers with applications in quantum photonics, sensing, and biomedicine.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Artificial correlation heating in PIC simulations

The Particle-in-Cell (PIC) method, a cornerstone in plasma modeling, is widely employed for its ability to simulate kinetic phenomena in device-scale domains. Part of what makes this possible is that computational macroparticles represent many physical particles. It converges under certain constraints, including a grid spacing that resolves the Debye length and a time step small enough to respect the Courant–Friedrichs–Lewy condition and plasma frequency stability limit. Here, we introduce a new constraint necessary to avoid Artificial Correlation Heating (ACH). This requires that the macroparticle coupling strength be smaller than one, Γ ω < 1, where Γ ω ≡ Γω 2/3 ⁠, Γ = Z 2 e 2 /(4πε 0 ak B T) is the physical coupling strength, and w is the macroparticle weight. This is particularly relevant to 3D simulations of dense plasmas, which are becoming common with modern computing power. If this condition is violated, the finite macroparticle weight artificially enhances the coupling strength and causes the plasma to heat until the macroparticle coupling strength is near unity, depending on the grid resolution. A comprehensive model of ACH is developed that incorporates electron density, temperature, macroparticle weight, and grid resolution. It is then tested using PIC simulations, delineating the boundaries of the method's applicability and offering a predictive framework for ACH. Moreover, the research explores a runaway heating process induced by ACH in the presence of ionization, which can lead to numerical instability. A conclusion of this study is that the onset of ACH can impose a more stringent constraint on the macroparticle weight and average number of macroparticles per cell than what is typically expected, particularly in 3D simulations of dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Computational design of high-entropy rare earth aluminum garnets for advanced thermal and environmental barrier coatings

To enhance the protection of Ni-based superalloys in gas turbine engines’ high-temperature environments, it’s crucial to develop advanced thermal/environmental barrier coating (T/EBC) materials with a balanced combination of thermal and mechanical properties. This optimization is essential to safeguard against chemical and thermal challenges. Here, in this study, we harness the power of density functional theory (DFT) in conjunction with combinatorial chemistry methodologies to engineer high-performance high-entropy rare earth disilicates of the RE 3 Al 5 O 12 family (where RE denotes Y, Gd, Er, and Yb). These materials are meticulously designed to exhibit superior phase stability, a targeted coefficient of thermal expansion (CTE), low lattice thermal conductivity, and robust mechanical properties. The determination of CTE values is accomplished through phonon calculations at various volume settings within the quasi-harmonic approximation, while lattice thermal conductivities are rigorously assessed employing the Debye-Callaway model, accounting for three distinct phonon processes. Our findings highlight the remarkable attributes of the solid solution (Y 1/4 Gd 1/4 Er 1/4 Yb 1/4 ) 3 Al 5 O 12 , which displays a reduction in lattice thermal conductivity compared to its individual constituents while maintaining a favorable range of CTE values. The novel T/EBC material, distinguished by their multifaceted functionalities, are poised to use in substantial enhancements in the performance of engines.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Full 3D Kinetic Modeling and Quantification of Positive Streamer Evolution in an Azimuthally Swept Pin-to-Plane Wedge Geometry

Cathode-directed streamer evolution in near atmospheric air is modeled in 3D pin-to-plane geometries using a 3D kinetic Particle-In-Cell (PIC) code that simulates particle-particle collisions via the Direct Simulation Monte Carlo (DSMC) method. Due to the computational challenges associated with a complete 360° volumetric domain, a practical alternative was achieved using a wedge domain and a range of azimuthal angles was explored (5°, 15°, 30°, and 45°) to study possible effects on the streamer growth and propagation due to the finite wedge angle. A DC voltage of 6 kV is administered to a hemispherical anode of radius 100 μm, with a planar cathode held at ground potential, generating an over-volted state with an electric field of 4 MV/m across a 1500 μm gap. The domain is seeded with an initial ion and electron density of 10 18 m -3 at 1 eV temperature confined to a spherical region of radius 100 μm centered at the tip of the anode. The air chemistry model [1] includes standard Townsend breakdown mechanisms (electron-neutral elastic, excitation, ionization, attachment, and detachment collision chemistry and secondary electron emission) as well as streamer mechanisms (photoionization and ion-neutral collisions) via tracking excited state neutrals which can then either quench via collisions or spontaneously emit a photon based on specific Einstein-A coefficients [2, 3]. In this work, positive streamer dynamics are formally quantified for each wedge angle in terms of electron velocity and density as temporal functions of coordinates r, Φ, and z. Applying a random plasma seed for each simulation, particles of interest are tracked with near femtosecond temporal resolution out to 1.4 ns and spatially binned. This process is repeated six times and results are averaged. Prior 2D studies have shown that the reduced electric field, E/n, can significantly impact streamer evolution [4]. We extend the analysis to 3D wedge geometries, to limit computational costs, and examine the wedge angle’s effect on streamer branching, propagation, and velocity. Results indicate that the smallest wedge angle that produced an acceptably converged solution is 30°. The potential effects that a mesh, when under-resolved with respect to the Debye length, can impart on streamer dynamics and numerical heating were not investigated, and we explicitly state here that the smallest cell size was approximately 10 times the minimum λ D in the streamer channel at late times. This constraint on cell size was the result of computational limitations on total mesh count.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Thermal Conductivity of β-Phase Ga 2 O 3 and (Al x Ga 1– x ) 2 O 3 Heteroepitaxial Thin Films

Heteroepitaxy of β-phase gallium oxide (β-Ga 2 O 3 ) thin films on foreign substrates shows promise for the development of next-generation deep ultraviolet solar blind photodetectors and power electronic devices. In this work, the influences of the film thickness and crystallinity on the thermal conductivity of ($\bar{2}01$)-oriented β-Ga 2 O 3 heteroepitaxial thin films were investigated. Unintentionally doped β-Ga 2 O 3 thin films were grown on c-plane sapphire substrates with off-axis angles of 0° and 6° toward $\langle$$11\bar{2}0$$\rangle$ via metal–organic vapor phase epitaxy (MOVPE) and low-pressure chemical vapor deposition. The surface morphology and crystal quality of the β-Ga 2 O 3 thin films were characterized using scanning electron microscopy, X-ray diffraction, and Raman spectroscopy. The thermal conductivities of the β-Ga 2 O 3 films were measured via time-domain thermoreflectance. The interface quality was studied using scanning transmission electron microscopy. The measured thermal conductivities of the submicron-thick β-Ga 2 O 3 thin films were relatively low as compared to the intrinsic bulk value. The measured thin film thermal conductivities were compared with the Debye–Callaway model incorporating phononic parameters derived from first-principles calculations. The comparison suggests that the reduction in the thin film thermal conductivity can be partially attributed to the enhanced phonon-boundary scattering when the film thickness decreases. They were found to be a strong function of not only the layer thickness but also the film quality, resulting from growth on substrates with different offcut angles. Growth of β-Ga 2 O 3 films on 6° offcut sapphire substrates was found to result in higher crystallinity and thermal conductivity than films grown on on-axis c-plane sapphire. However, the β-Ga 2 O 3 films grown on 6° offcut sapphire exhibit a lower thermal boundary conductance at the β-Ga 2 O 3 /sapphire heterointerface. In addition, the thermal conductivity of MOVPE-grown ($\bar{2}01$)-oriented β-(Al x Ga 1–x ) 2 O 3 thin films with Al compositions ranging from 2% to 43% was characterized. Because of phonon-alloy disorder scattering, the β-(Al x Ga 1–x ) 2 O 3 films exhibit lower thermal conductivities (2.8–4.7 W/m∙K) than the β-Ga 2 O 3 thin films. The dominance of the alloy disorder scattering in β-(Al x Ga 1–x ) 2 O 3 is further evidenced by the weak temperature dependence of the thermal conductivity. This work provides fundamental insight into the physical interactions that govern phonon transport within heteroepitaxially grown β-phase Ga 2 O 3 and (Al x Ga 1–x ) 2 O 3 thin films and lays the groundwork for the thermal modeling and design of β-Ga 2 O 3 electronic and optoelectronic devices.

36 MATERIALS SCIENCE↗