Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic inverse modeling”

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

Kinetic beaming in radiative relativistic magnetic reconnection: a mechanism for rapid gamma-ray flares in jets

ABSTRACT Rapid gamma-ray flares pose an astrophysical puzzle, requiring mechanisms both to accelerate energetic particles and to produce fast observed variability. These dual requirements may be satisfied by collisionless relativistic magnetic reconnection. On the one hand, relativistic reconnection can energize gamma-ray emitting electrons. On the other hand, as previous kinetic simulations have shown, the reconnection acceleration mechanism preferentially focuses high energy particles – and their emitted photons – into beams, which may create rapid blips in flux as they cross a telescope’s line of sight. Using a series of 2D pair-plasma particle-in-cell simulations, we explicitly demonstrate the critical role played by radiative (specifically inverse Compton) cooling in mediating the observable signatures of this ‘kinetic beaming’ effect. Only in our efficiently cooled simulations do we measure kinetic beaming beyond one light crossing time of the reconnection layer. We find a correlation between the cooling strength and the photon energy range across which persistent kinetic beaming occurs: stronger cooling coincides with a wider range of beamed photon energies. We also apply our results to rapid gamma-ray flares in flat-spectrum radio quasars, suggesting that a paradigm of radiatively efficient kinetic beaming constrains relevant emission models. In particular, beaming-produced variability may be more easily realized in two-zone (e.g. spine-sheath) set-ups, with Compton seed photons originating in the jet itself, rather than in one-zone external Compton scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

An Assessment of Molten Metal Detachment Hazards During Electron Beam Welding in Space

The safety issue has been raised with regards to potential molten metal detachments from the weld pool and cold filler wire during electron beam welding in space. This investigation was undertaken to evaluate if molten metal could detach and come in contact with astronauts and burn through the fabric of the astronauts' Extravehicular Mobility Unit (EMU) during electron beam welding in space. Molten metal detachments from either the weld/cut substrate or weld wire could present harm to a astronaut if the detachment was to burn through the fabric of the EMU. Theoretical models were developed to predict the possibility and size of the molten metal detachment hazards during the electron beam welding exercises at Low Earth Orbit (LEO). The primary molten metal detachment concerns were those cases of molten metal separation from the metal surface due to metal cutting, weld pool splashing, entrainment and release of molten metal due to filler wire snap-out from the weld puddle, and molten metal accumulation and release from the end of the weld wire. Some possible ways of obtaining molten metal drop detachments would include an impulse force, or bump, to the weld sample, cut surface, or filler wire. Theoretical models were developed for these detachment concerns from principles of impact and kinetic energies, surface tension, drop geometry, surface energies, and particle dynamics. The surface tension represents the force opposing the liquid metal drop from detaching whereas the weight of the liquid metal droplet represents a force that is tending to detach the molten metal drop. Theoretical calculations have indicated that only a small amount of energy is required to detach a liquid metal drop; however, much of the energy of an impact is absorbed in the sample or weld plate before it reaches the metal drop on the cut edge or surface. The tendency for detachment is directly proportional to the weld pool radius and metal density and inversely proportional to the surface tension of the liquid metal. For a detachment the initial kinetic energy of the weld pool with respect to the plate has to exceed the energy to form the extra surface required for the detachment of the pool. The difficulty is in transferring the energy from the point of impact through the plate and sample to the cut edge. It is likely that not all of the kinetic energy is available for detaching the pool; some may be sequestered in weld pool oscillations. The coefficient of restitution for the collision will be lower than one if irreversible deformation, for example plastic flow deformation, takes place during the collision. Thus determining the amount of energy from an impact that actually reaches the molten metal droplet is critical. Various molten metal detachment scenarios were tested experimentally in an enclosed vacuum chamber using the Ukrainian Universal Hand Tool, an electron beam welder designed for space welding. The experimental testing was performed in a 4 ft. X 4 ft. vacuum chamber at Marshall Space Flight Center, evacuated to vacuum levels of at least 50 microTorr, and also some welding garment material was utilized to observe the effect of the molten metal detachments on the material. A "carillon" apparatus consisting of four pendulum hammer strikers, each weighing approximately 3.65 lbs, raised to predetermined specific heights was used to apply an impact force to the weld sample/plate during electron beam welding and cutting exercises. The strikers were released by switching on an electric motor to rotate a pin holding wires retaining the strikers at desired heights. The specimens were suspended so as to be free to respond to the blows with a sudden velocity increment. The specimens were mounted on a hinged plate for minimizing effective mass with the option to fasten it down so as to raise its effective mass closer to that anticipated for an actual space welding scenario. Measurements were made of the impact energy and the horizontal fling distances of the detached metal drops. It was not particularly easy to generate the detachments for this experiment. This document presents the details of the theoretical modeling effort and a summary of the experimental effort to measure molten metal drop detachments from terrestrial electron beam welding in the enclosed vacuum chamber. The results of the experimental effort have shown that molten metal detachments can occur from the sample/weld plate only if a sufficiently large impact force is applied to the weld plate. A "weld pool detachment parameter" was determined to indicate whether detachment would occur. Detachment can be either full or partial (dripping), Partial detachment means that the weld pool detached from one side of the liquid-solid boundary so as to leave a hole at the puddle site but remained attached over part of the liquid-solid boundary and dripped down the plate with no fully detached material detected. Full detachment, however, does not necessarily mean that the whole pool fully detached; in some cases only a smaller portion of the pool detached, the remainder dripping down the plate. The weld pool detachment parameter according to theory and according to the empirical data allows a determination of whether full detachments might occur. Theoretical calculations indicated titanium alloy would be the most difficult from which to detach molten metal droplets followed by stainless steel and then by aluminum. The experimental results were for the most part consistent with the theoretical analysis and predictions. The above theory is applicable to other situations as desired for assessing the potential for molten metal detachments.

Fragomeni, James M.↗
DOE OSTI2025

Combined TDLAS and chemiluminescence imaging in a flat flame burner operated with NH3/H2 blends

Ammonia is viewed as a viable hydrogen carrier due to favorable storage and transport characteristics. While it can be re-converted to hydrogen at point-of-use via thermal catalytic cracking, direct utilization in combustion systems can result in reduced costs and improved efficiency. A major barrier to this approach is the low flammability and potential for high nitrogen oxide emissions, driven by fuel-bound nitrogen and complex kinetic pathways. While a number of kinetic mechanisms currently exist for simulating ammonia combustion, a major need continues to be direct information about species profiles in easy-to-model systems capable of isolating chemical kinetics from multi-dimensional fluid dynamic effects. This paper reports on recent species measurements made in a flat flame burner using a combined tunable-diode-laser-absorption-spectroscopy (TDLAS) and chemiluminescence imaging approach. Three flame conditions were included representative of NH3/H2 blends, partially cracked NH3 (inc. N2), and 100% NH3 with enhanced air (30% oxygen). Two NIR distributed feedback (DFB) laser diodes were used to determine H2O concentration and temperature (via ratio thermometry) at various radial and axial positions, after which an inverse Abel transform was used to infer centerline values. Companion images were collected using a Princeton Instrument PI-MAX intensified camera equipped with a 105mm UV lens and multiple filter sets targeting OH*, NH*, and NH2* emission. Results were compared to companion Cantera burner-stabilized flame simulations using various kinetic mechanisms. A recent mechanism including excited species chemistry was also included, to investigate whether excited and ground state profiles exhibit significant differences.

ammonia combustion↗
DOE OSTI2025

In Situ Characterization of Surface Recombination in p-Si/SiO x Based Photoelectrochemical Cells

Time-resolved infrared (TRIR) and electrochemical impedance spectroscopy (EIS) were utilized to quantify surface states present at silicon oxide (SiO x ) protected crystalline p-Si electrolyte interfaces. The primary goal was to identify p-Si/SiO x photoelectrodes with both low surface recombination rates and efficient multi-electron transfer to an acceptor present in the external electrolyte. Three SiO x layers were investigated: native oxide (nOx), chemical oxide (cOx), and rapid thermal annealed (RTA) thermal oxide (tOx). Comparative study with [Ru(bpy) 3 ](PF 6 ) 2 as the electron acceptor indicated that tOx was most optimal with a small effective recombination rate, multi-electron transfer capability, and photovoltage of 500 ± 50 mV. A secondary goal was to analyze the surface recombination rates with the Shockley–Read–Hall (SRH) kinetic model. Two surface states were identified from this analysis, one closer to the CB edge (V t,1 ) and the other near the midgap (V t,2 ). EIS and SRH analyses revealed that a forming gas (5% H 2 /N 2 ) anneal (FGA) decreased surface recombination for tOx and nOx through a lower density of surface states. In the case of tOx, the infrared data indicated that V t,2 was completely removed. Here, the energetic positions of the band edges were correlated with the surface state density; low densities corresponded to more favorable potentials for inversion layer formation, which is expected to be most optimal for photocatalysis. Collectively this study indicates that the free carrier dynamics provided by TRIR represent a powerful in situ probe of the band edge and the surface state energetics in silicon based photoelectrochemical cells.

Electrochemical Impedance Spectroscopy↗
DOE OSTI2024

Impact of Momentum Perturbation on Convective Boundary Layer Turbulence

Mesoscale-to-microscale coupling is an important tool for conducting turbulence-resolving multiscale simulations of realistic atmospheric flows, which are crucial for applications ranging from wind energy to wildfire spread studies. Different techniques are used to facilitate the development of realistic turbulence in the large-eddy simulation (LES) domain while minimizing computational cost. Here, we explore the impact of a simple and computationally efficient Stochastic Cell Perturbation method using momentum perturbation (SCPM-M) to accelerate turbulence generation in boundary-coupled LES simulations using the Weather Research and Forecasting model. We simulate a convective boundary layer (CBL) to characterize the production and dissipation of turbulent kinetic energy (TKE) and the variation of TKE budget terms. Furthermore, we evaluate the impact of applying momentum perturbations of three magnitudes below, up to, and above the CBL on the TKE budget terms. Momentum perturbations greatly reduce the fetch associated with turbulence generation. When applied to half the vertical extent of the boundary layer, momentum perturbations produce an adequate amount of turbulence. However, when applied above the CBL, additional structures are generated at the top of the CBL, near the inversion layer. The magnitudes of the TKE budgets produced by SCPM-M when applied at varying heights and with different perturbation amplitudes are always higher near the surface and inversion layer than those produced by No-SCPM, as are their contributions to the TKE. This study provides a better understanding of how SCPM-M reduces computational costs and how different budget terms contribute to TKE in a boundary-coupled LES simulation.

17 WIND ENERGY↗
DOE OSTI2023

Fission with Exotic Nuclei (Abbreviated Report)

Nuclear fission is a key mechanism involved in the synthesis of heavy elements in the Cosmos and is the primary explanation for the stability of superheavy elements. Nevertheless, our knowledge of fission remains extremely fragmented. Most experiments have been conducted only on a tiny number of stable actinide nuclei and are often incomplete, leading to gaps in our basic understanding of the process. For many radioactive isotopes, basic fission data such as the charge or mass distribution of the fragments is unknown. These gaps cannot always be filled by simulation alone. Common fission models contain too many free parameters and lack predictive power. In contrast, the fundamental theory of fission under development at LLNL is much more predictive, but its current computational cost is too high to be used extensively for data evaluations. A unique window of opportunity to resolve these limitations has recently opened: the U.S. nuclear science community is ramping up major experimental programs at the Facility for Rare Isotope Beams (FRIB, the DOE flagship facility in low-energy nuclear science), and techniques from machine learning have shown great potential to simplify the use of a fundamental, quantum-mechanical theory of fission. This project has two components. On the experimental side, we acquired and deployed at the HIGS facility a new dual Frisch-Grid ionization chamber to measure correlated fragment-mass, kinetic energy, and angular distributions of fission fragments from induced fission. This new device was used to perform measurements of charge, mass and total kinetic energy of fission fragments in the photofission of 238 U and eight gamma-ray beam energies between 6.2 and 13 MeV, which allowed extracting high-precision independent yields for this reaction. The device was also used to perform measurements of the same quantities in the neutron-induced fission of 234 U with monoenergetic beams of energy between 5 and 8 MeV. In parallel, we collaborated with a team at Commissariat à l’énergie atomique et aux énergies alternatives (CEA) to perform a series of measurements of fission yields in inverse kinematics for the two isotopes of 236 U and 240 Pu. The experiment took place at the Grand Accélérateur National d’Ions Lourds in France in June 2023. The deployment of the VAMOS spectrometer with a new array called PISTA allowed determining the excitation energy of the fissioning system within 1 Mega-electronvolts. The second component of the project involved using deep neural networks to build fast and reliable emulators of our current fission models. In an invited paper published in Frontier in Physics, we showed that autoencoders could successfully compress nuclear wavefunctions in nuclear density functional theory. We achieved a dimensionality reduction of the order of two orders of magnitude while keeping the error in the total energy to less than 0.01%. In a second paper submitted to Physical Review Letters in June 2023 with our collaborators at CEA, we showed that variational autoencoders can learn the collective degrees of freedom driving the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1977

A numerical study of the unsteady leading edge separation bubble on an oscillating airfoil

A preliminary study of the unsteady viscous flow in the region of an airfoil leading edge was performed, in which the interaction between the viscous and inviscid flow fields is neglected. The solution method uses the finite difference form of the governing equations throughout the separated flow field and incorporates a transition model based on the integral turbulence kinetic energy equation. The validity of the numerical procedure is verified by making comparisons with analytical solutions to several test problems, including unsteady flow over a plate oscillating in its own plane. The method was then applied to the problem of unsteady viscous flow over a NACA 0012 airfoil oscillating sinusoidally in pitch. The flow field characteristics were in qualitative agreement with experimental results. The bubble moved forward on the airfoil and decreased in size as incidence was increased. Viscous flow in the leading edge region was found to be quasi-steady, while bubble height varied inversely with Reynolds number.

Kreskovsky, J. P.↗
NASA NTRS1991

A review of the kinetic detail required for accurate predictions of normal shock waves

Several aspects of the kinetic models used in the collision phase of Monte Carlo direct simulations have been studied. Accurate molecular velocity distribution function predictions require a significantly increased number of computational cells in one maximum slope shock thickness, compared to predictions of macroscopic properties. The shape of the highly repulsive portion of the interatomic potential for argon is not well modeled by conventional interatomic potentials; this portion of the potential controls high Mach number shock thickness predictions, indicating that the specification of the energetic repulsive portion of interatomic or intermolecular potentials must be chosen with care for correct modeling of nonequilibrium flows at high temperatures. It has been shown for inverse power potentials that the assumption of variable hard sphere scattering provides accurate predictions of the macroscopic properties in shock waves, by comparison with simulations in which differential scattering is employed in the collision phase. On the other hand, velocity distribution functions are not well predicted by the variable hard sphere scattering model for softer potentials at higher Mach numbers.

Muntz, E. P.↗
DOE OSTI2024

In situ aerobic methane oxidation rates in a stratified lake

Abstract Microbial aerobic methane oxidation is an important sink for aquatic methane worldwide. Despite its importance to global methane fluxes, few aerobic methane oxidation rates have been obtained in freshwater or marine environments without imposing changes to the microbial community through use of ex situ methods. A novel in situ incubation method for continuous time‐series measurements was used in Jordan Lake, North Carolina, during 2020–2021, to determine reaction kinetics for aerobic methane oxidation rates across a wide range of naturally varying methane (55–1833 nM) and dissolved oxygen (DO; 28–366 μM) concentrations and temperatures (17–30°C). Methane oxidation began immediately at the start of each of 21 incubations and methane oxidation rates were 1 st order with respect to methane. The data density allowed for accurate calculation of 1 st ‐order rate constants,k, that ranged from 0.018 to 0.462 h −1 (R 2 > 0.967). Addition of ammonium (20–45 μM) to natural concentrations ranging from 0.057 to 2.4 μM did not change aerobic methane oxidation rate kinetics, suggesting that the natural population of aerobic methane oxidizers in this eutrophic lake was not nitrogen limited. Values ofkinversely correlated most strongly with initial DO concentrations (R 2 = 0.82) rather than temperature. Values forkincreased with Julian day throughout our sampling period, suggesting seasonal influences on methane oxidation via responses to geochemical changes or shifts in microbial community abundance and composition. These experiments demonstrate a high variability in the enzymatic capacity for 1 st ‐order methane oxidation rates in this eutrophic lake that is tightly and inversely coupled to oxygen concentrations. Measurements of in situ aerobic methane oxidation rate constants allow for the direct quantification and modeling of the microbial community's capacity for methane oxidation over a wide range of natural methane concentrations.

Marine & Freshwater Biology↗
DOE OSTI2021

Self-consistent core-pedestal ITER scenario modeling

The purpose of these integrated ITER simulations is to identify dependencies that can impact the performance of ITER. The dependence of fusion power production, temperature and density pedestals, and core profiles on varying magnetic-q, edge density fueling strength, neoclassical transport, magnetic field strength, and alpha heating, are examined. It is noted that ion modes usually dominate the interior, while electron modes dominate the barrier. The electron mode without flow shear is substantially stronger. The q increase tends to cause more kinetic ballooning modes at the edge and the q decrease means that the peeling is getting worse, but the ballooning is getting better. The large edge q tends to provide hollow density while small edge q gives normal density profile. This is consistent with the observation that the density of particles is inversely proportional to the density of the plasma current. The rise in the source of the edge particle increases the density of the edge and the reaction rate close to the edge increases temperature fluxes, resulting in weak pedestal barriers. When neoclassical transport is turned off or the strength of B-field is increased, the pronounced edge barriers are identified. In addition, a distinct effect on the magnetohydrodynamic (MHD) modes is observed due to the increased B-field. The interesting aspect is that the barrier collapses about q = 2 as the stability of MHD has decreased to replace the barrier with a continuous slope but with roughly the same fusion power. This would actually be good for a reactor. The slope of the H-mode pedestal is found to be reduced due to the alpha heating. It is expected that the problem of edge-localized modes causing damage to the first wall will be eased. The effects of circular, elongated, and general geometry on plasma profiles are also compared and contrasted. The circular geometry has a hollow density profile and no edge or internal thermal transport barrier (ETB or ITB). However, both ETBs and ITBs can be found in the elongated geometry. In general geometry, there is a weaker ETB, comparable ITBs, but a higher edge particle transport barrier than in elongated geometry. Higher density near the edge, on the other hand, causes more wall erosion, so elongated geometry might be the best option.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Production and Persistence of Extreme Two-temperature Plasmas in Radiative Relativistic Turbulence

Turbulence is a predominant process for energizing electrons and ions in collisionless astrophysical plasmas, and thus is responsible for shaping their radiative signatures (luminosity, spectra, and variability). To better understand the kinetic properties of a collisionless radiative plasma subject to externally driven turbulence, we investigate particle-in-cell simulations of relativistic plasma turbulence with external inverse Compton cooling acting on the electrons. We find that ions continuously heat up while electrons gradually cool down (due to the net effect of radiation), and hence the ion-to-electron temperature ratio T i /T e grows in time. Furthermore, we show that Ti/Te is limited only by the size and duration of the simulations (reaching ${T}_{i}/{T}_{e}\sim {10}^{3}$), indicating that there are no efficient collisionless mechanisms of electron–ion thermal coupling. This result has implications for models of radiatively inefficient accretion flows, such as observed in the Galactic center and in M87, for which so-called two-temperature plasmas with ${T}_{i}/{T}_{e}\gg 1$ have been invoked to explain their low luminosity. Additionally, we find that electrons acquire a quasi-thermal distribution (dictated by the competition of turbulent particle energization and radiative cooling), while ions undergo efficient nonthermal acceleration (acquiring a harder distribution than in equivalent nonradiative simulations). There is a modest nonthermal population of high-energy electrons that are beamed intermittently in space, time, and direction; these beamed electrons may explain rapid flares in certain high-energy astrophysical systems (e.g., in the Galactic center). These numerical results demonstrate that extreme two-temperature plasmas can be produced and maintained by relativistic radiative turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

ALESEP: A computer program for the analysis of airfoil leading edge separation bubbles

The ALESEP program for the analysis of the inviscid/viscous interaction which occurs due to the presence of a closed laminar transitional separation bubble on an airflow is presented. The ALESEP code provides a iterative solution of the boundary layer equations expressed in an inverse formulation coupled to a Cauchy integral representation of the inviscid flow. This interaction analysis is treated as a local perturbation to a known solution obtained from a global airfoil analysis. Part of the required input to the ALESEP code are the reference displacement thickness and tangential velocity distributions. Special windward differencing may be used in the reversed flow regions of the separation bubble to accurately account for the flow direction in the discretization of the streamwise convection of momentum. The ALESEP code contains a forced transition model based on a streamwise intermittency function and a natural transition model based on a solution of the integral form of the turbulent kinetic energy equation. Instructions for the input/output, and program usage are presented.

Davis, R. L.↗
NASA NTRS1986

ALESEP. Part 2: A computer program for the analysis of leading edge separation bubbles on infinite swept wings

A program called ALESEP is presented for the analysis of the inviscid-viscous interaction which occurs due to the presence of a closed laminar-transitional separation bubble on an airfoil or infinite swept wing. The ALESEP code provides an iterative solution of the boundary layer equations expressed in an inverse formulation coupled to a Cauchy integral representation of the inviscid flow. This interaction analysis is treated as a local perturbation to a known solution obtained from a global airfoil analysis; hence, part of the required input to the ALESEP code are the reference displacement thickness and tangential velocity distributions. Special windward differencing may be used in the reversed flow regions of the separation bubble to accurately account for the flow direction in the discretization of the streamwise convection of momentum. The ALESEP code contains a forced transition model based on a streamwise intermittency function, a natural transition model based on a solution of the integral form of the turbulent kinetic energy equation, and an empirical natural transition model.

Davis, R. L.↗
DOE OSTI2023

An Analytic Formula for Entraining CAPE in Midlatitude Storm Environments

Abstract This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from an environmental sounding, rather than prescribed by the formula user. Entrainment is connected to the background environment using an eddy diffusivity approximation for lateral mixing, updraft geometry assumptions, and mass continuity. These approximations result in a direct correspondence between the storm-relative flow and the updraft radius and an inverse scaling between the updraft radius squared and entrainment rate. The aforementioned concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, the ECAPE formula is modified to account for updraft enhancement via kinetic energy extracted from the cloud’s background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity w max in an updraft. Hence, these formulas are evaluated using w max from past numerical modeling studies. Both of the new formulas improve predictions of w max substantially over commonly used diagnostic parameters, including undiluted CAPE and ECAPE with a constant prescribed entrainment rate. The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance of by w max , and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies. Significance Statement Substantial mixing occurs between the upward-moving air currents in thunderstorms (updrafts) and the surrounding comparatively dry environmental air, through a process called entrainment. Entrainment controls thunderstorm intensity via its diluting effect on the buoyancy of air within updrafts. A challenge to representing entrainment in forecasting and predictions of the intensity of updrafts in future climates is to determine how much entrainment will occur in a given thunderstorm environment without a computationally expensive high-resolution simulation. To address this gap, this article derives a new formula that computes entrainment from the properties of a single environmental profile. This formula is shown to predict updraft vertical velocity more accurately than past diagnostics, and can be used in forecasting and climate prediction to improve predictions of thunderstorm behavior and impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability

Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1979

Possible inverse cascade behavior for drift-wave turbulence

The turbulent spectral properties of the dynamical equation of Hasegawa and Mima (1978) governing the evolution of the electrostatic potential in drift-wave turbulence is investigated for two formulations of the problem: (1) as a nondissipative initial value problem, with the potential represented by a truncated Fourier series with large number of terms, and (2) as a dissipative problem with a small viscous dissipation at very short spatial scales, and a long wavelength forcing term at longer wavelengths. It is found that Hasegawa and Mima's prediction for the nondissipative, truncated initial value modal problem is accurate, but substantial differences exist for the forced dissipative case between computer results and analytical predictions based on a wave kinetic equation of Kadomtsev. Much better agreement is found with a simple dual-cascade model based on Kraichnan's generalization of Kolmogorov's cascade arguments.

Fyfe, D.↗
NASA NTRS1981

Hot-gas cold-dust pumping for water masers associated with H II regions

A collisional pump with an internal sink is proposed for the water masers associated with H II regions, where the population inversion occurs due to the absorption by cold ice-mantle grains in a highly dusty cloud of the far-infrared line radiation of hot water vapor. A new escape probability method is developed to calculate the transfer of line radiation in dusty medium. The pump mechanism explains the power of usual maser sources associated with H II regions and the enormous power of the sources associated with W49 N and external galaxies. Models of maser clouds have a radius of 5 x 10 to the 15th-10 to the 16th cm, an H2 number density of 4 x 10 to the 9th/cu cm, an expansion velocity of 10-30 km/s, a kinetic temperature of 350 K, and a grain temperature of 100 K. Giant maser sources require grains of the size about 1 micron. The apparent size of the emission spots (approximately 10 to the 13th cm) observed by VLBI is interpreted as due to a fluctuation in the cloud, and the assembly of the spots is spread within a size of 10 to the 16th cm. The temperature difference between the dust and gas is due to a relaxation process after an infrared burst accompanying protostar formation.

Deguchi, S.↗
DOE OSTI2021

Multi-scale transport in the DIII-D ITER baseline scenario with direct electron heating and projection to ITER

Multi-scale fluctuations measured by turbulence diagnostics spanning long and short wavelength spatial scales impact energy confinement and the scale-lengths of plasma kinetic profiles in the DIII-D ITER baseline scenario with direct electron heating. Contrasting discharge phases with ECH + neutral beam injection (NBI) and NBI only at similar rotation reveal higher energy confinement and lower fluctuations when only NBI heating is used. Modeling of the core transport with TGYRO using the TGLF turbulent transport model and NEO neoclassical transport reproduces the experimental profile changes upon application of direct electron heating and indicates that multi-scale transport mechanisms are responsible for changes in the temperature and density profiles. Intermediate and high-k fluctuations appear responsible for the enhanced electron thermal flux, and intermediate-k electron modes produce an inward particle pinch that increases the inverse density scale length. Projection to ITER is performed with TGLF and indicates a density profile that has a finite scale length due to intermediate-k electron modes at low collisionality and increases the fusion gain. For a range of E×B shear, the dominant mechanism that increases fusion performance is suppression of outward low-k particle flux and increased density peaking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗