Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak turbulence”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

HST Studies of the Chromospheres, Wind, and Mass-Loss Rates of Cool Giant and Supergiant Stars

UV spectra of K-M giant and supergiant stars and of carbon stars have been acquired with the Goddard High Resolution Spectrograph (GHRS) on the Hubble Space Telescope (HST). These spectra have been used to measure chromospheric flow and turbulent velocities, study the acceleration of their stellar winds, acquire constraints on their outer atmospheric structure, and enable estimates of their mass-loss rates. Results from our observations of the giant stars Gamma Dra (K5 III hybrid), Alpha Tau (K5 III), Gamma Cru (M3.4 III), Mu Gem (M3 IIIab), and 30 Her (MG III), the supergiants Alpha Ori (M2 Iab) and Lambda Vel (K5 Ib), and the carbon stars TX Psc (NO; C6,2) and TW Hor (NO; C7,2) will be summarized and compared. The high resolution and wavelength accuracy of these data have allowed the direct measurement of the acceleration of the stellar winds in the chromospheres of several of these stars (from initial velocities of 3-9 km/s to upper velocities of 15-25 km/s) and of the chromospheric macroturbulence (-25-35 km/s). The high signal-to-noise and large dynamic range of these spectra have allowed the detection and identification of numerous new emission features, including weak C IV emission indicative of hot transition-region plasma in the non-coronal giant Alpha Tau, many new fluorescent lines of Fe II, and the first detection of fluorescent molecular hydrogen emission and of Ca II recombination lines in the UV spectrum of a giant star. The UV spectrum of two carbon stars have been studied with unprecedented resolution and reveal extraordinarily complicated Mg II lines nearly smothered by circumstellar absorptions. Finally, comparison of synthetic UV emission line profiles computed with the Lamers et al. (1987) Sobolev with Exact Integration (SEI) code with observations of chromospheric emission lines overlain with wind absorption features provides estimates of the mass-loss rates for four of these stars.

Carpenter, Kenneth G.↗
NASA NTRS1994

Acoustic instability driven by cosmic-ray streaming

We study the linear stability of compressional waves in a medium through which cosmic rays stream at the Alfven speed due to strong coupling with Alfven waves. Acoustic waves can be driven unstable by the cosmic-ray drift, provided that the streaming speed is sufficiently large compared to the thermal sound speed. Two effects can cause instability: (1) the heating of the thermal gas due to the damping of Alfven waves driven unstable by cosmic-ray streaming; and (2) phase shifts in the cosmic-ray pressure perturbation caused by the combination of cosmic-ray streaming and diffusion. The instability does not depend on the magnitude of the background cosmic-ray pressure gradient, and occurs whether or not cosmic-ray diffusion is important relative to streaming. When the cosmic-ray pressure is small compared to the gas pressure, or cosmic-ray diffusion is strong, the instability manifests itself as a weak overstability of slow magnetosonic waves. Larger cosmic-ray pressure gives rise to new hybrid modes, which can be strongly unstable in the limits of both weak and strong cosmic-ray diffusion and in the presence of thermal conduction. Parts of our analysis parallel earlier work by McKenzie & Webb (which were brought to our attention after this paper was accepted for publication), but our treatment of diffusive effects, thermal conduction, and nonlinearities represent significant extensions. Although the linear growth rate of instability is independent of the background cosmic-ray pressure gradient, the onset of nonlinear eff ects does depend on absolute value of DEL (vector differential operator) P(sub c). At the onset of nonlinearity the fractional amplitude of cosmic-ray pressure perturbations is delta P(sub C)/P(sub C) approximately (kL) (exp -1) much less than 1, where k is the wavenumber and L is the pressure scale height of the unperturbed cosmic rays. We speculate that the instability may lead to a mode of cosmic-ray transport in which plateaus of uniform cosmic-ray pressure are separated by either laminar or turbulent jumps in which the thermal gas is subject to intense heating.

Begelman, Mitchell C.↗
DOE OSTI2025

Assessment of flamelet/progress variable methods for supersonic combustion

Tabulated chemistry models, including the flamelet/progress variable approach, have been successfully used for a variety of turbulent flame simulations. The progress variable describes the progress of reactions in a system and parameterizes a lookup table of thermochemical variables. This approach reduces the cost of simulations, transporting only one scalar (progress variable) instead of the many species mass fractions required for detailed chemistry. Originally developed for low Mach number flame simulations, recent works have focused on extensions of this approach to compressible flames, supersonic combustion, and detonations, with applications such as scramjet combustors and rotating detonation engines. Unlike low Mach simulations, compressible flow simulations require solving the energy transport equation, which is coupled to the equation of state. This leads to additional modeling challenges regarding the thermodynamics and its impact on the chemistry. The validity of modeling assumptions, for example the relationship between energy and temperature, also varies with the combustion regime. The present work provides a detailed assessment of the existing strategies for chemistry tabulation for compressible/supersonic combustion, including detonations. A priori analysis indicates that approximations which are reasonable for weakly compressible flames may break down for shock-induced combustion. Furthermore, the analysis identifies specific assumptions and approximations that do not hold for detonations, emphasizing that care must be taken when applying tabulated chemistry models outside their intended combustion regimes.

Detonations↗
NASA NTRS1999

Towards Understanding the Mechanism of Receptivity and Bypass Dynamics in Laminar Boundary Layers

Three problems concerning laminar-turbulent transition are addressed by solving a series of initial value problems. The first problem is the calculation of resonance within the continuous spectrum of the Blasius boundary layer. The second is calculation of the growth of Tollmien-Schlichting waves that are a direct result of disturbances that only lie outside of the boundary layer. And, the third problem is the calculation of non-parallel effects. Together, these problems represent a unified approach to the study of freestream disturbance effects that could lead to transition. Solutions to the temporal, initial-value problem with an inhomogeneous forcing term imposed upon the flow is sought. By solving a series of problems, it is shown that: A transient disturbance lying completely outside of the boundary layer can lead to the growth of an unstable Tollmien-Schlichting wave. A resonance with the continuous spectrum leads to strong amplification that may provide a mechanism for bypass transition once nonlinear effects are considered. A disturbance with a very weak unstable Tollmien-Schlichting wave can lead to a much stronger Tollmien-Schlichting wave downstream, if the original disturbance has a significant portion of its energy in the continuum modes.

Lasseigne, D. G.↗
NASA NTRS2016

Nonlinear Wakes Behind a Row of Elongated Roughness Elements

This paper is concerned with the high Reynolds number flow over a spanwise periodic array of roughness elements with inter-element spacing of the order of the local boundary-layer thickness. While earlier work by Goldstein, Sescu, Duck and Choudhari (2010) and Goldstein, Sescu, Duck and Choudhari (2011) was mainly concerned with smaller roughness heights that produced relatively weak distortions of the downstream flow, the focus here is on extending the analysis to larger roughness heights and streamwise elongated planform shapes that together produce a qualitatively different, nonlinear behavior of the downstream wakes. The roughness scale flow now has a novel triple-deck structure that is somewhat different from related studies that have previously appeared in the literature. The resulting flow is formally nonlinear in the intermediate wake region, where the streamwise distance is large compared to the roughness dimensions but small compared to the downstream distance from the leading edge, as well as in the far wake region where the streamwise length scale is of the order of the downstream distance from the leading edge. In contrast, the flow perturbations in both of these wake regions were strictly linear in the earlier work by Goldstein et al (2010, 2011). This is an important difference because the nonlinear wake flow in the present case provides an appropriate basic state for studying the secondary instability and eventual breakdown into turbulence.

Boundary layer control↗
NASA NTRS2014

The Role of Higher-Order Modes on the Electromagnetic Whistler-Cyclotron Wave Fluctuations of Thermal and Non-Thermal Plasmas

Two fundamental challenging problems of laboratory and astrophysical plasmas are the understanding of the relaxation of a collisionless plasmas with nearly isotropic velocity distribution functions and the resultant state of nearly equipartition energy density with electromagnetic plasma turbulence. Here, we present the results of a study which shows the role that higher-order-modes play in limiting the electromagnetic whistler-like fluctuations in a thermal and non-thermal plasma. Our main results show that for a thermal plasma the magnetic fluctuations are confined by regions that are bounded by the least-damped higher order modes. We further show that the zone where the whistler-cyclotron normal modes merges the electromagnetic fluctuations shifts to longer wavelengths as the beta(sub e) increases. This merging zone has been interpreted as the beginning of the region where the whistler-cyclotron waves losses their identity and become heavily damped while merging with the fluctuations. Our results further indicate that in the case of nonthermal plasmas, the higher-order modes do not confine the fluctuations due to the effective higher-temperature effects and the excess of suprathermal plasma particles. The analysis presented here considers the second-order theory of fluctuations and the dispersion relation of weakly transverse fluctuations, with wave vectors parallel to the uniform background magnetic field, in a finite temperature isotropic bi-Maxwellian and Tsallis-kappa-like magnetized electron-proton plasma. Our results indicate that the spontaneously emitted electromagnetic fluctuations are in fact enhanced over these quasi modes suggesting that such modes play an important role in the emission and absorption of electromagnetic fluctuations in thermal or quasi-thermal plasmas.

Vinas, Adolfo F.↗
DOE OSTI2026

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗
DOE OSTI2020

Multimodel validation of single wakes in neutral and stratified atmospheric conditions

Previous research has revealed the need for a validation study that considers several wake quantities and code types so that decisions on the trade-off between accuracy and computational cost can be well informed and appropriate to the intended application. In addition to guiding code choice and setup, rigorous model validation exercises are needed to identify weaknesses and strengths of specific models and guide future improvements. Here, we consider 13 approaches to simulating wakes observed with a nacelle-mounted lidar at the Scaled Wind Technology Facility (SWiFT) under varying atmospheric conditions. We find that some of the main challenges in wind turbine wake modeling are related to simulating the inflow. In the neutral benchmark, model performance tracked as expected with model fidelity, with large-eddy simulations performing the best. In the more challenging stable case, steady-state Reynolds-averaged Navier–Stokes simulations were found to outperform other model alternatives because they provide the ability to more easily prescribe noncanonical inflows and their low cost allows for simulations to be repeated as needed. Dynamic measurements were only available for the unstable benchmark at a single downstream distance. These dynamic analyses revealed that differences in the performance of time-stepping models come largely from differences in wake meandering. This highlights the need for more validation exercises that take into account wake dynamics and are able to identify where these differences come from: mesh setup, inflow, turbulence models, or wake-meandering parameterizations. In addition to model validation findings, we summarize lessons learned and provide recommendations for future benchmark exercises.

17 WIND ENERGY↗
DOE OSTI2024

The influence of the cloud virial parameter on the initial mass function

ABSTRACT Crucial for star formation is the interplay between gravity and turbulence. The observed cloud virial parameter, $\alpha _{\mathrm{vir}}$, which is the ratio of twice the turbulent kinetic energy to the gravitational energy, is found to vary significantly in different environments, where the scatter among individual star-forming clouds can exceed an order of magnitude. Therefore, a strong dependence of the initial mass function (IMF) on $\alpha _{\mathrm{vir}}$ may challenge the notion of a universal IMF. To determine the role of $\alpha _{\mathrm{vir}}$ on the IMF, we compare the star-particle mass functions obtained in high-resolution magnetohydrodynamical simulations including jet and heating feedback, with $\alpha _{\mathrm{vir}}=0.0625$, 0.125, and 0.5. We find that varying $\alpha _{\mathrm{vir}}$ from $\alpha _{\mathrm{vir}}\sim 0.5$ to $\alpha _{\mathrm{vir}}< 0.1$ shifts the peak of the IMF to lower masses by a factor of $\sim 2$ and increases the star formation rate by a similar factor. The dependence of the IMF and star formation rate on $\alpha _{\mathrm{vir}}$ is non-linear, with the dependence subsiding at $\alpha _{\mathrm{vir}}< 0.1$. Our study shows a systematic dependence of the IMF on $\alpha _{\mathrm{vir}}$. Yet, it may not be measurable easily in observations, considering the uncertainties, and the relatively weak dependence found in this study.

Mathew, Sajay Sunny (ORCID:0000000283818195)↗
NASA NTRS1993

Three-dimensional numerical simulation of near-surface flows over the Martian north polar cap

Measurements made by Viking Lander VL-2 (48 N) have shown that the near-surface wind and temperature regime on Mars displays striking similarities to terrestrial counterparts. The diurnal radiative cycle is responsible for establishment of a well-defined thermal circulation in which downslope (Katabatic) flows prevail during the nighttime hours and weak upslope (anabatic) conditions prevail during the daytime. Previous work has indicated that the slope flows are much like those found on Earth, particularly the Katabatic winds, which show striking similarities to drainage flows observed over Antarctica. The low-level wind regime appears to be an important factor in the scouring of the martian landscape. The north polar cap shows evidence of eolian features such as dunes, frost streaks, and grooves from Viking imagery. The direction of the prevailing wind can in cases be inferred from the eolian features. We examine the thermally induced flows that result from the radiative heating and cooling of the martian north polar region using a comprehensive three-dimensional atmospheric mesoscale numerical model. The same model has been used previously for simulation of Antarctic Katabatic winds. The model equations are written in terrain-following coordinates to allow for irregular terrain; prognostic equations include the flux forms of the horizontal momentum equations, temperature, continuity. A surface energy budget equation is also incorporated in which the surface temperature is determined. Explicit parameterization of both terrestrial (longwave) and solar (shortwave) radiation is included. Turbulent transfer of heat and momentum in the martian atmosphere is assumed to follow the similarity expressions in the surface boundary layer on Earth. The terrain heights for the martian north polar region have been obtained from the U.S. Geological Survey map and digitized onto a 57x57 grid with a spacing of 75 km. The resulting terrain map is shown in Fig. 1. The vertical grid consists of 15 levels (delta equals 0.998, 0.99, 0.98, 0.97, 0.96, 0.94, 0.92, 0.90, 0.85, 0.775, 0.70, 0.60, 0.50, 0.30, 0.10). The high resolution in the lower atmosphere is necessary to capture details of the boundary layer flows. The lowest level corresponds to a height of approximately 20 m above the ground, the seconds level 100 m.

Parish, Thomas R.↗
DOE OSTI2022

Dimensionless parameters for cloudy Rayleigh-Bénard convection: Supersaturation, Damköhler, and Nusselt numbers

In steady-state Rayleigh-Bénard convection, heat is transported by turbulent thermal convection from the bottom, hot surface to the top, cold surface, leading to a height-independent sensible heat flux. When water vapor is present and cloud formation occurs, there is also an additional latent heat flux. Heat transport in cloudy Rayleigh-Bénard convection depends on turbulent flow as well as the microphysical state of the clouds: specifically, whether substantial supersaturations exist and whether cloud liquid water is removed through sedimentation/precipitation. In this article we bridge between the Rayleigh-Bénard convection literature and the atmospheric literature. We express the governing equations for cloudy convection in dimensionless form, thereby explicitly identifying the governing parameters relevant to the cloudy case, including Schmidt, Damköhler, supersaturation, and sedimentation numbers. We further connect to the atmospheric literature by obtaining a Nusselt number (dimensionless heat flux) for a cloud-convection system, directly from the conservation equations for temperature and water vapor. This flux has the same form as that identified by Zhang et al. [L. Zhang, K. L. Chong, and K.-Q. Xia, J. Fluid Mech. 874, 1041 (2019)] for convection with water vapor, but is extended to the cloudy case. For equal thermal and water vapor diffusivities, the flux corresponds to the widely used atmospheric quantities equivalent temperature and moist static energy. Using large eddy simulation (LES) of an idealized cloudy Rayleigh-Bénard convection system with fixed boundary conditions, we find that the equivalent heat flux (Nusselt number) is only weakly dependent on the microphysical details of the system, such as liquid water mixing ratio and cloud droplet number concentration. Finally, from the results, we show the vertical profiles of sensible and latent heat fluxes depend on the liquid water content, whereas the equivalent heat flux remains a constant throughout the height of the chamber.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2022

Patches of Magnetic Switchbacks and Their Origins

Parker Solar Probe (PSP) has shown that the solar wind in the inner heliosphere is characterized by the quasi omnipresence of magnetic switchbacks ("switchback" hereinafter), local backward bends of magnetic field lines. Switchbacks also tend to come in patches, with a large-scale modulation that appears to have a spatial scale size comparable to supergranulation on the Sun. Here we inspect data from the first 10 encounters of PSP focusing on different time intervals when clear switchback patches were observed by PSP. We show that the switchbacks modulation, on a timescale of several hours, seems to be independent of whether PSP is near perihelion, when it rapidly traverses large swaths of longitude remaining at the same heliocentric distance, or near the radial-scan part of its orbit, when PSP hovers over the same longitude on the Sun while rapidly moving radially inwards or outwards. This implies that switchback patches must also have an intrinsically temporal modulation most probably originating at the Sun. Between two consecutive patches, the magnetic field is usually very quiescent with weak fluctuations. We compare various parameters between the quiescent intervals and the switchback intervals. The results show that the quiescent intervals are typically less Alfvénic than switchback intervals, and the magnetic power spectrum is usually shallower in quiescent intervals. We propose that the temporal modulation of switchback patches may be related to the "breathing" of emerging flux that appears in images as the formation of "bubbles" below prominences in the Hinode/SOT observations.

Solar wind↗
DOE OSTI2020

Planck 2018 results XII. Galactic astrophysics using polarized dust emission

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle ψ, and dispersion function of polarization angles S . The joint distribution (one-point statistics) of p and NH confirms that the mean and maximum polarization fractions decrease with increasing NH. The uncertainty on the maximum observed polarization fraction, pmax = 22.0 -1.4 +3.5% at 353 GHz and 80' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and S found earlier is seen to be present at high latitudes. This follows the S ∝ p -1 relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, ψ, and S for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map S × p, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between NH = 10 20 cm -2 and NH = 2 x 10 22 cm -2 , out of the Galactic plane, this product S × p only decreases by about 25%. Because S is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of S × p with the dust temperature T d , whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion S tend, on the contrary, to have colder dust than lines of sight with low p and high S . We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R P/p and R S/V , which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p V /E(B - V), of 13% at high Galactic latitude, compatible with the polarization fraction p ≈ 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

NIMROD Development AND Applications for Advanced Simulations of Tokamak Plasmas

Magnetic fusion experiments have the largest gradients in the world that are in steady state. These gradients lead to instabilities, and thus, most of plasma theory over the past four decades have been devoted to the development of instability theory and their use in interpreting plasma results. Instabilities in tokamaks lead to 3 different phenomenological outcomes: coherent saturation, turbulent saturation, and sudden relaxation (which may include disruptions). For long-wavelength instabilities, the NIMROD extended MHD code has emerged as an important tool for understanding tokamak instabilities. Because of the long history of NIMROD, a new version of it has been started at Tech-X to be able to address multi-species capabilities as well as exploit modern GPU systems. The new version builds on the previous version and improves the workflow by enabling new equilibria from experiments to be resolved. The fundamental formulation of the equations underlying NIMROD has also been explored in multiples ways. First, understanding instabilities in tokamaks through numerical methods is aided by being able to understand the rich history of analytic studies. A paper was published which aid in understanding the literature by simplifying some of the analytic machinery inherent in these studies. Extended MHD, like gyrokinetics, are quasineutral models. From a theoretical point of view, plasma quasineutrality approximation is best understood as the same as the magneto-quasistatic approximation of the Maxwell equations. Finally, the standard model of tokamak theory is that of instability theory. A simple dynamical systems model has been developed to better illustrate the strengths and weaknesses of this model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Initial Results of DC Electric Fields, Associated Plasma Drifts, Magnetic Fields, and Plasma Waves Observed on the C/NOFS Satellite

Initial results are presented from the Vector Electric Field Investigation (VEFI) on the Air Force Communication/Navigation Outage Forecasting System (C/NOFS) satellite, a mission designed to understand, model, and forecast the presence of equatorial ionospheric irregularities. The VEFI instrument includes a vector DC electric field detector, a fixed-bias Langmuir probe operating in the ion saturation regime, a flux gate magnetometer, an optical lightning detector, and associated electronics including a burst memory. Compared to data obtained during more active solar conditions, the ambient DC electric fields and their associated E x B drifts are variable and somewhat weak, typically < 1 mV/m. Although average drift directions show similarities to those previously reported, eastward/outward during day and westward/downward at night, this pattern varies significantly with longitude and is not always present. Daytime vertical drifts near the magnetic equator are largest after sunrise, with smaller average velocities after noon. Little or no pre-reversal enhancement in the vertical drift near sunset is observed, attributable to the solar minimum conditions creating a much reduced neutral dynamo at the satellite altitude. The nighttime ionosphere is characterized by larger amplitude, structured electric fields, even where the plasma density appears nearly quiescent. Data from successive orbits reveal that the vertical drifts and plasma density are both clearly organized with longitude. The spread-F density depletions and corresponding electric fields that have been detected thus far have displayed a preponderance to appear between midnight and dawn. Associated with the narrow plasma depletions that are detected are broad spectra of electric field and plasma density irregularities for which a full vector set of measurements is available for detailed study. Finally, the data set includes a wide range of ELF/VLF/HF oscillations corresponding to a variety of plasma waves, in particular banded ELF hiss, whistlers, and lower hybrid wave turbulence triggered by lightning-induced sferics. The VEFI data represents a new set of measurements that are germane to numerous fundamental aspects of the electrodynamics and irregularities inherent to the Earth's low latitude ionosphere.

Pfaff, R.↗
NASA NTRS1994

Modification of the MML turbulence model for adverse pressure gradient flows

Computational fluid dynamics is being used increasingly to predict flows for aerospace propulsion applications, yet there is still a need for an easy to use, computationally inexpensive turbulence model capable of accurately predicting a wide range of turbulent flows. The Baldwin-Lomax model is the most widely used algebraic model, even though it has known difficulties calculating flows with strong adverse pressure gradients and large regions of separation. The modified mixing length model (MML) was developed specifically to handle the separation which occurs on airfoils and has given significantly better results than the Baldwin-Lomax model. The success of these calculations warrants further evaluation and development of MML. The objective of this work was to evaluate the performance of MML for zero and adverse pressure gradient flows, and modify it as needed. The Proteus Navier-Stokes code was used for this study and all results were compared with experimental data and with calculations made using the Baldwin-Lomax algebraic model, which is currently available in Proteus. The MML model was first evaluated for zero pressure gradient flow over a flat plate, then modified to produce the proper boundary layer growth. Additional modifications, based on experimental data for three adverse pressure gradient flows, were also implemented. The adapted model, called MMLPG (modified mixing length model for pressure gradient flows), was then evaluated for a typical propulsion flow problem, flow through a transonic diffuser. Three cases were examined: flow with no shock, a weak shock and a strong shock. The results of these calculations indicate that the objectives of this study have been met. Overall, MMLPG is capable of accurately predicting the adverse pressure gradient flows examined in this study, giving generally better agreement with experimental data than the Baldwin-Lomax model.

Conley, Julianne M.↗
DOE OSTI2022

Finding evidence for inflation and the origin of galactic magnetic fields with CMB surveys

The origin of the μG magnetic fields observed in galaxies is unknown. One promising scenario is that magnetic fields generated during inflation, larger than 0.1 nG on Mpc scales, were adiabatically compressed to μG strengths in galaxies during structure formation. Thus, detecting a scale-invariant primordial magnetic field (PMF) above 0.1 nG on Mpc scales just after recombination would indicate an inflationary origin of galactic magnetic fields. Furthermore, this would also provide compelling evidence that inflation occurred since only an inflationary mechanism could generate such a strong, scale-invariant magnetic field on Mpc scales. In contrast, constraining the scale-invariant PMF strength to be below 0.1 nG would imply an inflationary scenario is not the primary origin, since such weak PMFs cannot be amplified enough via adiabatic compression to produce the strength of the galactic fields we observe today. We find that measurements of anisotropic birefringence by future CMB surveys will be able to improve the sensitivity to Mpc-scale inflationary PMFs by an order of magnitude, and, in particular, that CMB-HD would lower the upper bound to 0.072 nG at the 95% C.L., which is below the critical 0.1 nG threshold for ruling out a purely inflationary origin. If inflationary PMFs exist, we find that a CMB-HD survey would be able to detect them with about 3σ significance or higher, providing evidence for inflation itself.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Progress in extending high poloidal beta scenarios on DIII-D towards a steady-state fusion reactor and impact of energetic particles

To prepare for steady-state operation of future fusion reactors (e.g. the International Thermonuclear Experimental Reactor and China Fusion Engineering Test Reactor (CFETR)), experiments on DIII-D have extended the high poloidal beta ( β P ) scenario to reactor-relevant edge safety factor q 95 ~ 6.0, while maintaining a large-radius internal transport barrier (ITB) using negative magnetic shear. Excellent energy confinement quality ( H 98y2 > 1.5) is sustained at high normalized beta ( β N ~ 3.5). This high-performance ITB state with Greenwald density fraction near 100% and q min ≥ 3 is achieved with toroidal plasma rotation V tor ~ 0 at ρ ≥ 0.6. This is a key result for reactors expected to have low V tor . At high β P (≥1.9), large Shafranov shift can stabilize turbulence leading to a high confinement state with a low pedestal and an ITB. At lower β P (<1.9), negative magnetic shear in the plasma core contributes to turbulence suppression and can compensate for reduced Shafranov shift to continue to access a large-radius ITB and excellent confinement with low V tor , consistent with the results of gyrofluid transport simulations. These high- β P cases are characterized by weak/no Alfvén eigenmodes (a.e.) and classical fast-ion transport. At high density, the fast-ion deceleration time decreases and Δ β fast is lower; these reduce a.e. drive. The reverse-shear Alfvén eigenmodes are weaker or stable because the negative magnetic shear region is located at higher radius, away from the peaked fast-ion profile. Resistive wall modes can be a limitation at simultaneous high β N , low internal inductance, and low rotation. Analysis suggests that additional off-axis external current drive could provide a more stable path at reduced q 95 . Based on a DIII-D high- β P plasma with large-radius ITB, two scenarios are proposed for CFETR Q = 5 steady-state operation with ~1 GW fusion power: a lower-$l_i$($l_i$ ~ 0.66) and a higher-$l_i$($l_i$ ~ 0.75) case. Using a Landau closure model, multiple energetic particle (EP) effects on the a.e. stability are analyzed modifying the growth rate of the a.e.s triggered by the neutral-beam-injection EPs and alpha particles, although the stabilizing/destabilizing effect is weak for the cases analyzed. The stabilizing effects of the combined EP species β , energy, and density profile in CFETR need further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗