SEARCH · Engineering Papers

Results for “Density gradient function”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Shortcomings of meta-GGA functionals when describing magnetism

Several recent studies have shown that SCAN, a functional belonging to the meta-generalized gradient approximation (MGGA) family, leads to significantly overestimated magnetic moments in itinerant ferromagnetic metals. However, this behavior is not inherent to the MGGA level of approximation since TPSS, for instance, does not lead to such severe overestimations. In order to provide a broader view of the accuracy of MGGAfunctionals for magnetism, we extend the assessment to more functionals but also to antiferromagnetic solids. The results show that to describe magnetism there is overall no real advantage in using a MGGA functional compared to GGAs. For both types of approximation, an improvement in ferromagnetic metals is necessarily accompanied by a deterioration (underestimation) in antiferromagnetic insulators, and vice versa. Furthermore, we also provide some analysis in order to understand in more detail the relation between the mathematical form of the functionals and the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Two-fluid and kinetic transport physics of Kelvin–Helmholtz instabilities in nonuniform low-beta plasmas

Hall-magnetohydrodynamic (Hall-MHD) theory, two-fluid simulations, and kinetic simulations are used in this work to investigate the cross-field transport properties of Kelvin–Helmholtz instabilities in nonuniform low-beta collisionless plasmas. Hall-MHD analysis shows how the linear properties of the instability are modified by density gradients and magnetization. High-order accurate two-fluid and kinetic simulations, with complete dynamics of finite-mass electrons and ions, are applied to a suite of parameter cases to systematically assess the effects of diamagnetic drift, magnetization, charge separation, and finite Larmor motion. Initialization of exact two-species kinetic equilibria facilitates the study of isolated physical effects and enables detailed cross-comparisons between two-fluid and kinetic simulations, including for cases where ion gyroradii are comparable to gradient scale lengths. For nonuniform plasmas with significant space charge, the results of two-fluid and kinetic simulations are found to disagree with Hall-MHD predictions. Kelvin–Helmholtz instability growth rates, per unit shear, are shown to be smaller when ion diamagnetic drift and $E \times B$ drift are parallel and larger when the two drifts are antiparallel. The effect is attributed to polarization drift in the shear layer, which leads to redistribution of charge, alters the electric field that drives plasma advection, and consequently modifies growth rates. Instability-induced mass transport for different parameters is characterized in terms of the flux across the shear layer and a simplified diffusion model. Distribution functions from kinetic simulations are shown to deviate substantially from Maxwellian reconstructions, indicating the importance of kinetic physics during the nonlinear phase of the instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Nonequilibrium steady-state thermoelectrics of Kondo-correlated quantum dots

The transport across a Kondo-correlated quantum dot coupled to two leads with independent temperatures and chemical potentials is studied using a controlled nonperturbative, and in this sense numerically exact, treatment based on a hybrid numerical renormalization group combined with time-dependent density matrix renormalization group (NRG-tDMRG). In the Kondo regime, for sufficiently large fixed voltage bias V ≳ T K , with T K the Kondo temperature, we find a peak in the conductance vs the temperature gradient Δ⁢T = T R - T L across left and right lead. Focusing then on zero voltage bias but finite ΔT far beyond linear response, we reveal the dependence of the characteristic zero-bias conductance on the individual lead temperatures. Here, we find that the finite-Δ⁢T data behaves quantitatively similar to linear response with an effective equilibrium temperature derived from the different lead temperatures. The regime of sign changes in the Seebeck coefficient, signaling the presence of Kondo correlations, and its dependence on the individual lead temperatures provide a complete picture of the Kondo regime in the presence of finite-temperature gradients. The results from the zero-bias conductance and Seebeck coefficient studies unveil an approximate “Kondo circle” in the T L /T R plane as the regime within which the Kondo correlations dominate. We also study the heat current and the corresponding heat conductance vs finite Δ⁢T. We provide a polynomial fit for our numerical results for the thermocurrent as a function of the individual lead temperatures, which may be used to fit experimental data in the Kondo regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Accessing the Subsurface Biosphere Within Rocks Undergoing Active Low-Temperature Serpentinization in the Samail Ophiolite (Oman Drilling Project)

The Oman Drilling Project established an “Active Alteration” multi-borehole observatory in peridotites undergoing low-temperature serpentinization in the Samail Ophiolite. The highly serpentinized rocks are in contact with strongly reducing fluids. Distinct hydrological regimes, governed by differences in rock porosity and fracture density, give rise to steep redox (Eh +200 to –750 mV) and pH (pH range 8.5–11.2) gradients within the 300–400 m deep boreholes. The serpentinites and fluids host an active subsurface ecosystem. Microbial cell abundances in serpentinite vary at least six orders of magnitude, from ≤3.5 × 10 1 to 2.9 × 10 7 cells/g. Low levels of biological sulfate reduction (2–1,000 fmol/cm 3 /day) can be detected in rock cores, particularly in rocks in contact with reduced groundwaters with pH < 10.5. Thermodesulfovibrio is the predominant sulfate reducer identified via metagenomic sequencing of adjacent groundwater communities. Here, we infer that transport and reaction of microbially generated sulfide with the serpentine and brucite assemblages gives rise to optical darkening and sulfide overprinting, including the formation of tochilinite-vallerite group minerals, potentially serving as an indicator that this system is inhabited by microbial life. Olivine mesh-cores replaced with ferroan brucite and minor awaruite, abundant veins containing hydroandradite garnet and polyhedral serpentine, and late-stage carbonate veins are suggested as targets for future spatially resolved life-detection investigations. The high-quality whole-round core samples that have been preserved can be further probed to define how life distributes itself and functions within a system where chemical disequilibria are sustained by low-temperature water/rock interaction, and how biosignatures of in situ microbial activity are generated.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Density Variations as a Function of Compostion in the Venusian Crust: Implications for Tesserae Formation

The Venusian surface plains are thought to be predominantly basaltic, while tesserae occupy ~ 7 % and are potentially the oldest units preserved [1]. They consist of layered, folded, and eroded rocks [2], potentially giving an insight into what tectonic regimes were present during Venus’ history. Studies of the tesserae have noted lower emissivity measurements compared to the surrounding plains [3]. Several explanations for these lower emissivity values have been proposed, including a more felsic composition for the tesserae [4], that the lower emissivity could be a result of a difference in mineralogy caused by chemical weathering facilitated between the rocks and atmosphere by the lower temperatures at the higher elevations of tesserae [5]. Alternatively, the morphology of some features on tesserae, such as the Ovda Fluctus, are more consistent with a basaltic composition [6]. The presence of significant amounts of felsic igneous compositions could either imply partial melting of hydrous material and hence a higher water content in early Venus [7], or differentiation of basaltic melts [8]. Further, tesserae composition influences the density of the crust with major implications for the stability of crustal roots. Here, we use thermodynamic modeling to compute phase diagrams and extract rock densities for two end-member compositions: a basalt and a granite. We discuss major phase and density changes for each composition and describe the effect of compositional variations along geothermal gradients with implications for the stability of crustal roots with the aim to provide constraints on the composition of tesserae.

A R Baker↗

NASA NTRS1996

Adaptive grid methods for RLV environment assessment and nozzle analysis

Rapid access to highly accurate data about complex configurations is needed for multi-disciplinary optimization and design. In order to efficiently meet these requirements a closer coupling between the analysis algorithms and the discretization process is needed. In some cases, such as free surface, temporally varying geometries, and fluid structure interaction, the need is unavoidable. In other cases the need is to rapidly generate and modify high quality grids. Techniques such as unstructured and/or solution-adaptive methods can be used to speed the grid generation process and to automatically cluster mesh points in regions of interest. Global features of the flow can be significantly affected by isolated regions of inadequately resolved flow. These regions may not exhibit high gradients and can be difficult to detect. Thus excessive resolution in certain regions does not necessarily increase the accuracy of the overall solution. Several approaches have been employed for both structured and unstructured grid adaption. The most widely used involve grid point redistribution, local grid point enrichment/derefinement or local modification of the actual flow solver. However, the success of any one of these methods ultimately depends on the feature detection algorithm used to determine solution domain regions which require a fine mesh for their accurate representation. Typically, weight functions are constructed to mimic the local truncation error and may require substantial user input. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. These weight functions can then be used to construct blending functions for algebraic redistribution, interpolation functions for unstructured grid generation, forcing functions to attract/repel points in an elliptic system, or to trigger local refinement, based upon application of an equidistribution principle. The popularity of solution-adaptive techniques is growing in tandem with unstructured methods. The difficultly of precisely controlling mesh densities and orientations with current unstructured grid generation systems has driven the use of solution-adaptive meshing. Use of derivatives of density or pressure are widely used for construction of such weight functions, and have been proven very successful for inviscid flows with shocks. However, less success has been realized for flowfields with viscous layers, vortices or shocks of disparate strength. It is difficult to maintain the appropriate mesh point spacing in the various regions which require a fine spacing for adequate resolution. Mesh points often migrate from important regions due to refinement of dominant features. An example of this is the well know tendency of adaptive methods to increase the resolution of shocks in the flowfield around airfoils, but in the incorrect location due to inadequate resolution of the stagnation region. This problem has been the motivation for this research.

Thornburg, Hugh J.↗

DOE OSTI2022

Comparing measurement of limiting current in block copolymer electrolytes as a function of salt concentration with theoretical predictions

Optimizing electrolyte performance is crucial for the widespread adoption of electrochemical energy storage. We demonstrate that limiting current provides a robust criterion for determining the optimum electrolyte. Experiments were conducted on rigid block copolymer electrolytes comprising mixtures of polystyrene-block-poly(ethylene oxide) copolymer (SEO) and lithium bis(trifluoromethanesulfonyl) imide salt (LiTFSI) over a salt concentration range from r av = 0.04 to r av = 0.20 (r av is the molar ratio of lithium ions to ethylene oxide). Here we show that the maximum limiting current density is 4.3 mA cm –2 at r av = 0.12. The dependence of limiting current on salt concentration is in good agreement with predictions from Newman's concentrated solution theory with no adjustable parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1990

Detection of a human intracisternal A-type retroviral particle antigenically related to HIV

Sjogren's syndrome is an autoimmune disease that is characterized by dryness of the mouth and eyes. The loss of salivary and lacrimal gland function is accompanied by lymphocytic infiltration. Because similar symptoms and glandular pathology are observed in certain persons infected with human immunodeficiency virus (HIV), a search was initiated for a possible retroviral etiology in this syndrome. A human intracisternal A-type retroviral particle that is antigenically related to HIV was detected in lymphoblastoid cells exposed to homogenates of salivary tissue from patients with Sjogren's syndrome. Comparison of this retroviral particle to HIV indicates that they are distinguishable by several ultrastructural, physical, and enzymatic criteria.

Non-NASA Center↗

DOE OSTI2026

Response of soil nutrient pools and microbiomes to recurrent wildfire disturbance and varying burn severities in a mixed conifer forest

Wildfire is a pervasive disturbance in mixed-conifer forests, yet the relative influence of fire recurrence versus burn severity on soil biogeochemistry and microbial communities remains poorly quantified. We examined a natural gradient of fire history (0–3 prior fires) and burn severity (low–high) spanning 50 yr in a mixed-conifer ecosystem to assess how repeated fire shapes soil carbon (C) and nitrogen (N) pools, their isotopic signatures, mineral and particulate fractions, microbial community composition, carbon-use, CO₂ fluxes, and vegetation cover. Successive fires produced progressively higher bare-ground percentages and lower tree cover, which were tightly linked to declines in microbial diversity and reductions bulk %C, and %N. δ 13 C increased with fire frequency, indicating preferential loss of labile C through combustion or enhanced microbial oxidation, thereby explaining the observed net soil-C decline. Conversely, δ 15 N decreased and pH increased as tree density declined, reflecting altered N cycling and reduced acidification in post-fire soils. Fire recurrence, more than severity, corresponded with a marked shift in the bacterial community: for example, Xanthobacteraceae—key N-fixers and C-cyclers—diminished, while N-fixing Bacillaceae increased, underscoring the tightly coupled nature of soil nutrient dynamics and microbiome composition after repeated burns. Our results demonstrate that fire recurrence appears to be a stronger driver of post-fire soil ecosystem responses in this mixed-conifer forest, influencing both abiotic nutrient pools and the functional potential of the soil microbiome. These findings provide a more enhanced assessment and understanding to date of the biogeochemical consequences of repeated wildfire disturbance that can be used to inform management strategies aimed at preserving soil health in fire-prone landscapes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Role of Intermolecular Interactions in Deep Eutectic Solvents for CO 2 Capture: Vibrational Spectroscopy and Quantum Chemical Studies

Recent research and reviews on CO 2 capture methods, along with advancements in industry, have highlighted high costs and energy-intensive nature as the primary limitations of conventional direct air capture and storage (DACS) methods. In response to these challenges, deep eutectic solvents (DESs) have emerged as promising absorbents due to their scalability, selectivity, and lower environmental impact compared to other absorbents. However, the molecular origins of their enhanced thermal stability and selectivity for DAC applications have not been explored before. Therefore, the current study focuses on a comprehensive investigation into the molecular interactions within an alkaline DES composed of potassium hydroxide (KOH) and ethylene glycol (EG). Combining Fourier transform infrared (FT-IR) and quantum chemical calculations, the study reports structural changes and intermolecular interactions induced in EG upon addition of KOH and its implications on CO 2 capture. Experimental and computational spectroscopic studies confirm the presence of noncovalent interactions (hydrogen bonds) within both EG and the KOH-EG system and point to the aggregation of ions at higher KOH concentrations. Additionally, molecular electrostatic potential (MESP) surface analysis, natural bond orbital (NBO) analysis, quantum theory of atoms-in-molecules (QTAIM) analysis, and reduced density gradient-noncovalent interaction (RDG-NCI) plot analysis elucidate changes in polarizability, charge distribution, hydrogen bond types, noncovalent interactions, and interaction strengths, respectively. Evaluation of explicit and hybrid models assesses their effectiveness in representing intermolecular interactions. This research enhances our understanding of molecular interactions in the KOH-EG system, which are essential for both the absorption and desorption of CO 2 . The study also aids in predicting and selecting DES components, optimizing their ratios with salts, and fine-tuning the properties of similar solvents and salts for enhanced CO 2 capture efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2012

Spatial Variability in Column CO2 Inferred from High Resolution GEOS-5 Global Model Simulations: Implications for Remote Sensing and Inversions

Column CO2 observations from current and future remote sensing missions represent a major advancement in our understanding of the carbon cycle and are expected to help constrain source and sink distributions. However, data assimilation and inversion methods are challenged by the difference in scale of models and observations. OCO-2 footprints represent an area of several square kilometers while NASA s future ASCENDS lidar mission is likely to have an even smaller footprint. In contrast, the resolution of models used in global inversions are typically hundreds of kilometers wide and often cover areas that include combinations of land, ocean and coastal areas and areas of significant topographic, land cover, and population density variations. To improve understanding of scales of atmospheric CO2 variability and representativeness of satellite observations, we will present results from a global, 10-km simulation of meteorology and atmospheric CO2 distributions performed using NASA s GEOS-5 general circulation model. This resolution, typical of mesoscale atmospheric models, represents an order of magnitude increase in resolution over typical global simulations of atmospheric composition allowing new insight into small scale CO2 variations across a wide range of surface flux and meteorological conditions. The simulation includes high resolution flux datasets provided by NASA s Carbon Monitoring System Flux Pilot Project at half degree resolution that have been down-scaled to 10-km using remote sensing datasets. Probability distribution functions are calculated over larger areas more typical of global models (100-400 km) to characterize subgrid-scale variability in these models. Particular emphasis is placed on coastal regions and regions containing megacities and fires to evaluate the ability of coarse resolution models to represent these small scale features. Additionally, model output are sampled using averaging kernels characteristic of OCO-2 and ASCENDS measurement concepts to create realistic pseudo-datasets. Pseudo-data are averaged over coarse model grid cell areas to better understand the ability of measurements to characterize CO2 distributions and spatial gradients on both short (daily to weekly) and long (monthly to seasonal) time scales

Ott, L.↗

DOE OSTI2022

Structural and dynamic properties of solvated hydroxide and hydronium ions in water from ab initio modeling

Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water from ab initio molecular dynamics (AIMD) has been a challenging task. The difficulty mainly comes from a lack of accurate and efficient exchange–correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange–correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice for describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. In conclusion, our AIMD simulations indicate that hydroxide diffuses more slowly than hydronium in water, which is consistent with the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

Geometrical evidence for dark matter: X-ray constraints on the mass of the elliptical galaxy NGC 720

We describe (1) a new test for dark matter and alternate theories of gravitation based on the relative geometries of the X-ray and optical surface brightness distributions and an assumed form for the potential, of the optical light, (2) a technique to measure the shapes of the total gravitating matter and dark matter of an ellipsoidal system which is insensitive to the precise value of the temperature of the gas and to modest temperature gradients, and (3) a new method to determine the ratio of dark mass to stellar mass that is dependent on the functional forms for the visible star, gas and dark matter distributions, but independent of the distance to the galaxy or the gas temperature. We apply these techniques to X-ray data from the ROSAT Position Sensitive Proportional Counter (PSPC) of the optically flattened elliptical galaxy NGC 720; the optical isophotes have ellipticity epsilon approximately 0.40 extending out to approximately 120 sec. The X-ray isophotes are significantly elongated, epsilon = 0.20-0.30 for semimajor axis a approximately 100 sec. The major axes of the optical and X-ray isophotes are misaligned by approximately 30 deg +/- 15 deg. Spectral analysis of the X-ray data reveals no evidence of temperature gradients or anisotropies and demonstrates that a single-temperature plasma (T approximately 0.6 keV) having subsolar heavy element abundances and a two-temperature model having solar abundances describe the spectrum equally well. Considering only the relative geometries of the X-ray and optical surface brightness distributions and an assumed functional form for the potential of the optical light, we conclude that matter distributed like the optical light cannot produce the observed ellipticities of the X-ray isophotes, independent of the gas pressure, the gas temperature, and the value of the stellar mass; this comparison assumes a state of quasi-hydrostatic equilibrium so that the three-dimensional surfaces of the gas emissivity trace the three-dimensional isopotential surfaces -- we discuss the viability of this assumption for NGC 720. Milgrom's Modification of Newtonian Dynamics (MOND) cannot dispel this manifestation of dark matter. Hence, geometrical considerations require, without mention of pressure or temperature, the presence of an extended, massive dark matter halo in NGC 720. Employing essentially the technique of Buote & Canizares (1992; Buote 1992) we use the shape of the X-ray surface brightness to constrain the shape of the total gravitating matter. The total matter is modeled as either an oblate or prolate spheriod of constant shape and orientation having either a Ferrers (rho approximately r(exp -n)) or Hernquist density. Assuming the X-ray gas is in hydrostatic equilibrium, we construct a model X-ray gas distribution for various temperature profiles. We determine the ellipticity of the total gravitating matter to be epsilon approximately 0.50-0.70. Using the single-temperature model we estimate a total mass approximately (0.41-1.4) x 10(exp 12) h(sub 80) solar mass interior to the ellipsoid of semimajor axis 43.6 h(sub 80) kpc. Ferrers densities as steep as r(exp -3) do not fit the data, but the r(exp -2) and Hernquist models yield excellent fits. We estimate the mass distributions of the stars and the gas and fit the dark matter directly. For a given gas equation of state and functional forms for the visible stars, gas, and dark matter, these models yield a distance-independent and temperature-independent measurement of the ratio of dark mass to stellar mass M(sub DM)/M(sub stars). We estimate a minimum M(sub DM)/M(sub stars) greater than or equal to 4 which corresponds to a total mass slightly greater than that derived from the single-temperature models for distance D = 20h(sub 80) Mpc.

Buote, David A.↗

NASA NTRS1992

Kinetic plasma processes occurring in the outer plasmasphere

One area of data analysis work that was begun under this contract is the fitting of the perpendicular velocity distributions of equatorially trapped ions with a Kappa function. This type of characterization of the trapped ions will be very useful for comparison with velocity distributions produced by the model. A second area of data analysis is to study data from consecutive passes when DE 1's apogee was near the magnetic equator and the spacecraft was often skimming along nearly the same L shell. In 1982 three such periods occurred in May, June, and July. For these consecutive events we have Kp histories, density measurements from a number of sources (Whistler data, DE SFR, ISEE SFR) and consecutive samples of ion pitch angle distributions along field lines. It is clear from this data how the pitch angle distributions evolve during a flux tube refilling event. Our modeling of the flow of plasma along closed field lines is following two basic tracks. The first is a study of the basic refilling process without the effect of wave-particle heating near the equator or the effect of large or abrupt field-aligned electric potential drops. This model includes the effects of Coulomb self-collisions and collisions with the O+ ions in the topside ionosphere. The second track is a study of the effects of wave produced pitch-angle scattering and perpendicular heating occurring near the magnetic equator, in connection with the development of large potential drops that result from electron heating and the development of density gradients.

Wilson, Gordon R.↗

NASA NTRS1998

The Effects of Ground and Space Processing on the Properties of Organic, Polymeric, and Colloidal Materials

In recent years, a great deal of interest has been directed toward the use of organic materials in the development of high-efficiency optoelectronic and phototonic devices. There is a myriad of possibilities among organic materials which allow flexibility in the design of unique structures with a variety of functional groups. The use of nonlinear optical (NLO) organic materials as thin film wave-guides allows full exploitation of their desirable qualifies by permitting long interaction lengths and large susceptibilities allowing modest power input. There are several methods in use to prepare thin films such as Langmuir-Blodgett (LB) and self-assembly techniques, vapor deposition, growth from sheared solution or melt, and melt growth between glass plates. Organic-based materials have many features that make them desirable for use in optical devices, such as high second-and third-order nonlinearity, flexibility of molecular design, and damage resistance to optical radiation. However, their use in devices has been hindered by processing difficulties for crystals and thin films. We discuss the potential role of microgravity processing of a few organic and polymeric materials. It is of interest to note how materials with second-and third-order NLO behavior may be improved in a diffusion-limited environment and ways in which convection may be detrimental to these materials. We focus our discussion on third-order materials for all-optical switching, and second-order materials for frequency conversion and electrooptics. The goal of minimizing optical loss obviously depends on processing methods. For solution-based processes, such as solution crystal growth and solution photopolymerization, it is well known that thermal and solutal density gradients can initiate buoyancy-driven convection. Resultant fluid flows can affect transport of material to and from growth interfaces and become manifest in the morphology and homogeneity of the growing film or crystal. Likewise, buoyancy-driven convection can hinder production of defect-free, high-quality crystals or films during crystal and film growth by vapor deposition.

Frazier, Donald O.↗

DOE OSTI2021

Electronic and optical properties of Zn-doped β-Ga 2 O 3 Czochralski single crystals

β-Ga 2 O 3 has several soluble deep acceptors that impart insulating behavior. In this study, we investigate Zn doping (0.25at.%) in bulk Czochralski and vertical gradient freeze β-Ga 2 O 3 . Representative crystals were assessed for orientation (electron backscatter diffraction and Raman spectroscopy), purity (glow discharge mass spectrometry and secondary ion mass spectrometry), optical properties (ultraviolet to near infrared absorption), and electrical properties (resistivity and current–voltage). Purity measurements indicate that Zn evaporation is insufficient to inhibit doping of Zn into β-Ga 2 O 3 . Hybrid functional calculations show Zn substitutes nearly equally on tetrahedral and octahedral sites, with less than ~0.1eV preference for the octahedral (Ga II ) site. Furthermore, calculations show that Zn Ga acts as a deep acceptor with trapping levels ~1.3 and ~0.9eV above the valence band for one and two holes, respectively. The solubility and electronic behavior of Zn dopants are consistent with measured concentrations >1×10 18 atoms/cm 3 and electrical measurements that show resistivity 10 11 –10 13 Ωcm, with no p-type conduction.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Investigating instabilities in magnetized low-pressure capacitively coupled RF plasma using particle-in-cell (PIC) simulations

The effect of a uniform magnetic field on particle transport in low-pressure radio frequency (RF) capacitively coupled plasma (CCP) has been studied using a particle-in-cell model. Three distinct regimes of plasma behavior can be identified as a function of the magnetic field. In the first regime at low magnetic fields, asymmetric plasma profiles are observed within the CCP chamber due to the effect of $\overrightarrow{E}$ × $\overrightarrow{B}$ drift. As the magnetic field increases, instabilities develop and form self-organized spoke-shaped structures that are distinctly seen within the bulk plasma closer to the sheath. In this second regime, the spoke-shaped coherent structures rotate inside the plasma chamber in the -$\overrightarrow{E}$ × $\overrightarrow{B}$ direction, where $\overrightarrow{E}$ × $\overrightarrow{B}$ are the DC electric and magnetic field vectors, respectively, and the DC electric field exists in the sheath and pre-sheath regions. The spoke rotation frequency is in the megahertz range. As the magnetic field strength increases further, the rotating coherent spokes continue to exist near the sheath. The coherent structures are, however, accompanied by new small-scale incoherent structures originating and moving within the bulk plasma region away from the sheath. This is the third regime of plasma behavior. The threshold values of the magnetic field between these regimes were found not to vary with changing plasma reactor geometry (e.g., area ratio between ground and powered electrodes) or the use of an external capacitor between the RF-powered electrode and the RF source. The threshold values of the magnetic field between these regimes shift toward higher values with increasing gas pressure. Analysis of the results indicates that the rotating structures are due to the lower hybrid instability driven by density gradients and electron-neutral collisions. This paper provides guidance on the upper limit of the magnetic field for instability-free operation in low-pressure CCP-based semiconductor deposition and etch systems that use the external magnetic field for plasma uniformity control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

Hydrodynamics of Two-Phase Flows through Porous Media in Microgravity: Packed Bed Reactor Experiment onboard of the International Space Station

The objective of the Packed Bed Reactor Experiment was to develop a fundamental understanding of the hydrodynamics of two-phase flow with no phase transition through porous media in microgravity. This work describes two experimental campaigns flown on the International Space Station (ISS) that were designed to achieve the objective. This work presents the results on flow patterns observed, two-phase flow pressure drop in porous media, and the impact of capillary effects on pressure drop at various gas and liquid flow rates. This work is the first to present predictive correlations of the two-phase friction factor for the different regimes identified based on the superficial liquid and gas velocity in microgravity, the first to hypothesize the different regimes based on the change in slope of the pressure gradient versus flow rate plots, the first to address the effects of the capillarity on the pressure gradient and the first to measure column holdup and assess the effective porosity based on the packing and the trapped gas bubbles. Experimental data on pressure drop in gas (N2) – liquid (water) flow show the presence of four different flow regimes in microgravity. Dispersed bubble flow (DB) and pulse (P) flow regimes are detected at high liquid flow rates, whereas at low liquid flow rates “large bubble” or elongated bubble (LB) and “gas channeling” (GC) regimes are observed. For these different flow regimes, different two-phase flow friction factor and pressure gradient correlations are presented for the first time as a function of the gas and liquid modified Reynolds numbers and Suratman number. Within the viscous-capillary (V-C) regime, it is found that the capillary contribution is the dominant force that contributes to the pressure drop for the wetting case (glass). However, for the non-wetting packing (Teflon), the viscous contribution dominates. It was found that the gas hold-up and pressure drop are functions of bed history at low liquid and gas flow rates with the magnitude of the hysteresis decreasing with increasing flow rates. PBRE-2 results show that the capillary force is a strong function of the superficial liquid velocity but is a much weaker function of the superficial gas velocity and varies inversely with the particle diameter. Within the Viscous-Capillary (V-C) regime, over 90% of the pressure gradient is attributed to the capillary contribution in the gas continuous regime. However, in the large bubble regime, the viscous and capillary contributions were comparable. After the completion of PBRE-2 with glass spherical beads, another bed packed with alumina was installed. Pressure gradient data for the alumina packed bed (PBRE-Water Recovery) were transferred to the PI and were not analyzed by the authors of this work. The Packed Bed Reactor Experiment concluded its on-orbit operation after two successful campaigns of testing in 2017 and in 2021. The flight hardware was brought back from orbit and is being reconditioned for a series of future experiments referred to as PBRE-WRS in support of water recovery (WRS). The WRS series consists of testing different two-phase fluid system components and packed beds to assess their pressure gradient characteristics, which will be used for designing packed bed reactors for various applications relevant to life support. Although there is so much relevance of the experimental results obtained from the two PBRE campaign in microgravity to life support, these findings also apply to other applications that involve two phase flow in porous medium such as fuel cells, transport of nutrient to plants in space and other chemical and materials processing that involve two phase flows. These systems operate differently in microgravity because, due to the lack of buoyancy, the density difference between the phases becomes irrelevant and no longer leads to phase separation.

Packed Bed↗