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 649 records · Page 36

The 67 Hz Feature in the Black Hole Candidate GRS 1915+105 as a Possible Diskoseismic Mode

The Rossi X-Ray Timing Explorer has made feasible for the first time the search for high-frequency (greater than or equal to 100 Hz) periodic features in Black Hole Candidate (BHC) systems. Such a feature, with a 67 Hz frequency, recently has been discovered in the BHC GRS 1915+105 (Morgan, Remillard, & Greiner). This feature is weak (rms variability approx. 0.3%-1.6%), stable in frequency (to within approx. 2 Hz) despite appreciable luminosity fluctuations, and narrow (quality factor Q approx. 20). Several of these properties are what one expects for a 'diskoseismic' g-mode in an accretion disk about a 10.6 M(solar mass) (nonrotating) to 36.3 M(solar mass) (maximally rotating) black hole (if we are observing the fundamental-mode frequency). We explore this possibility by considering the expected luminosity modulation, as well as possible excitation and growth mechanisms-including turbulent excitation, damping, and 'negative' radiation damping. We conclude that a diskoseismic interpretation of the observations is viable.

Nowak, Michael A.↗

A Scalable Real-Time Data Assimilation Framework for Predicting Turbulent Atmosphere Dynamics

AI-based foundation models like FourCastNet, GraphCast are revolutionizing weather and climate predictions but are not yet ready for operational use. Their limitation lies in the absence of a data assimilation system to incorporate real-time Earth system observations, crucial for accurately forecasting events like tropical cyclones. To overcome these obstacles, we introduce a generic real-time data assimilation framework and demonstrate its end-to-end performance on the Frontier supercomputer. This framework comprises two primary modules: an ensemble score filter (EnSF), which significantly outperforms the state-of-the-art data assimilation method, and a vision transformer-based surrogate capable of real-time adaptation through the integration of observational data. We demonstrate both the strong and weak scaling of our framework up to 1024 GPUs on the Exascale supercomputer, Frontier. Our results not only illustrate the framework's exceptional scalability on high-performance computing systems, but also demonstrate the importance of supercomputers in real-time data assimilation for weather and climate predictions.

Lu, Dan↗

Improving GCM Predictability of Mixed-Phase Clouds and Aerosol Interactions at High Latitudes with ARM Observations

The goal of this project is to improve the predictability of mixed-phase clouds and aerosol interactions in the Community Atmosphere Model version 6 (CAM6) through comparison with the ARM observations. There are three main objectives of the proposed study: (1) Improve the representation of ice microphysical processes in mixed-phase clouds; (2) Test the performance of ice microphysics in CESM-CAM6 with the ARM observations in northern and southern high latitudes; and (3) Examine mixed-phase cloud microphysics-aerosol-turbulence-radiation interactions in CESM-CAM6. In this project, we have (1) Improved the representation of ice microphysical processes in mixed-phase clouds in CESM-CAM6 by implementing the marine organic aerosol (MOA) and treating the ice nucleating particles (INPs) from MOA and its impacts on mixed-phase clouds. We made a first attempt to represent different secondary ice production (SIP) mechanisms in a GCM (CESM2-CAM6). We found that misrepresentation of these ice formation processes in global climate models (GCMs) leads to too weak negative cloud feedback over the Southern Ocean (SO) and too high climate sensitivity in the models. In addition to ice formation processes, we improved the treatment of ice depositional growth through the Wegener–Bergeron–Findeisen (WBF) process by considering the subgrid heterogeneous distributions between liquid droplets and ice crystals in mixed-phase clouds; (2) Tested the performance of our improved representations of ice microphysics with the ARM observations at high latitudes.

54 ENVIRONMENTAL SCIENCES↗

Fluctuating wall pressures near separation in highly swept turbulent interactions

Fluctuating wall pressures have been measured in Mach 5 interactions generated by sharp, upswept fins at angles of attack of 16 to 28 deg. The results show that rms pressure distributions, like the mean, can be collapsed in conical coordinates. The wall pressure signal near separation is intermittent and is qualitatively similar to that measured in unswept interactions and other swept flows. However, the dominant separation shock frequencies in the swept flows are up to an order of magnitude higher than those in unswept interactions under identical incoming flow conditions. In light of the present observations, it appears that an earlier remark, by Gilson and Dolling, that separation is characterized by a shuddering compression system, in contrast to a translating separation shock is erroneous due to a combination of weak shock strength and inadequate spatial resolution. Furthermore, although the maximum rms pressure near separation increases with increasing interaction strength, as in swept compression ramp flows, comparison of data from the two flow types indicates that the appropriate correlating parameter is the interaction sweepback angle.

Schmisseur, J. D.↗

Theory of two-point correlations of jet noise

A large body of careful experimental measurements of two-point correlations of far field jet noise was carried out. The model of jet-noise generation is an approximate version of an earlier work of Ribner, based on the foundations of Lighthill. The model incorporates isotropic turbulence superimposed on a specified mean shear flow, with assumed space-time velocity correlations, but with source convection neglected. The particular vehicle is the Proudman format, and the previous work (mean-square pressure) is extended to display the two-point space-time correlations of pressure. The shape of polar plots of correlation is found to derive from two main factors: (1) the noncompactness of the source region, which allows differences in travel times to the two microphones - the dominant effect; (2) the directivities of the constituent quadrupoles - a weak effect. The noncompactness effect causes the directional lobes in a polar plot to have pointed tips (cusps) and to be especially narrow in the plane of the jet axis. In these respects, and in the quantitative shapes of the normalized correlation curves, results of the theory show generally good agreement with Maestrello's experimental measurements.

Ribner, H. S.↗

r -process Nucleosynthesis from Hyperaccreting Neutron Stars in Common Envelopes

We investigate nuclear reactions and feedback in hyperaccreting neutron star environments, considering accretion rates in the range 0.3–3 × 10 4 M ⊙ yr −1 , typical of short-period compact-object binaries in common envelopes. Our models account for weak reactions, neutrino energy loss, nuclear energy release, pair production, degenerate equations of state, and general relativistic hydrodynamics. Depending on the accretion rates, these systems can develop both proton- and neutron-rich atmospheres with strong convective instabilities linking the neutrino sphere to the outgoing accretion shock inside the radiation trapping zone. Convection drives nucleons through multiple heating and cooling cycles, with photodisintegration dominating during the heating phase and heavy element synthesis during the cooling phase, ejecting material with abundances that depend on the accretion rate and depth of the final decompression trajectory. The turbulent nature of convective currents is conducive to creating a wide range of nuclear products through a variety of effects, including nuclear statistical equilibrium freeze-out and the r-, p-, and γ-processes. We also observe a novel multistep process in reheated trajectories, consisting of proton-capture and photodissociation reactions operating on r-process seeds, producing overall neutron-deficient isotopes. A significant amount of infalling gas experiences high entropy and short (millisecond) freeze-out timescales capable of making r-process elements with high overabundances through a disequilibrium effect between neutrons and α-particles that does not require an excess of neutrons.

79 ASTRONOMY AND ASTROPHYSICS↗

Environmental sensitivities of shallow-cumulus dilution – Part 1: Selected thermodynamic conditions

Cumulus entrainment, and its consequent dilution of buoyant cloud cores, strongly regulates the life cycle of shallow cumuli yet remains poorly understood. Herein, new insights into this problem are obtained through large-eddy simulations that systematically investigate the sensitivity of shallow-cumulus dilution to cloud-layer relative humidity (RH), cloud- and subcloud-layer depths, and continentality (i.e., the land–ocean contrast). The simulated cloud-core dilution is found to be strongly sensitive to continentality, with fractional dilution rates twice as large over the ocean as over land. Using a similarity theory based on the turbulent-kinetic-energy (TKE) budget, the reduced cloud-core dilution over land is attributed to larger cloud-base mass flux (mb), driven by stronger surface heating and subcloud turbulence. As mb increases, the fractional dilution rate must decrease to maintain energetic equilibrium. A positive sensitivity is also found to cloud-layer RH, with the core dilution increasing by 25 %–50 % for a 10 % enhancement in RH. This sensitivity is interpreted using the buoyancy-sorting hypothesis, in that mixtures of cloud and environmental air are more likely to become negatively buoyant and detrain (rather than diluting the cloud core) in drier cloud layers. By contrast, the sensitivities of (marine) shallow-cumulus dilution to cloud- and subcloud-layer depths are weak, with a 3 % decrease for a doubling for the former and a 4 % reduction in dilution for a 50 % deeper subcloud layer. These surprisingly weak sensitivities are readily explained by offsetting effects in the TKE similarity theory. Altogether, these experimental findings provide useful, though still incomplete, guidance for flow-dependent shallow-cumulus entrainment formulations in large-scale models.

54 ENVIRONMENTAL SCIENCES↗

Numerical simulations of laser-driven experiments of ion acceleration in stochastic magnetic fields

We present numerical simulations used to interpret laser-driven plasma experiments at the GSI Helmholtz Centre for Heavy Ion Research. The mechanisms by which non-thermal particles are accelerated in astrophysical environments, e.g., the solar wind, supernova remnants, and gamma ray bursts, is a topic of intense study. When shocks are present, the primary acceleration mechanism is believed to be first-order Fermi, which accelerates particles as they cross a shock. Second-order Fermi acceleration can also contribute, utilizing magnetic mirrors for particle energization. Despite this mechanism being less efficient, the ubiquity of magnetized turbulence in the universe necessitates its consideration. Another acceleration mechanism is the lower-hybrid drift instability, arising from gradients of both density and magnetic field, which produce lower-hybrid waves with an electric field that energizes particles as they cross these waves. With the combination of high-powered laser systems and particle accelerators, it is possible to study the mechanisms behind cosmic-ray acceleration in the laboratory. In this work, we combine experimental results and high-fidelity three-dimensional simulations to estimate the efficiency of ion acceleration in a weakly magnetized interaction region. We validate the FLASH magneto-hydrodynamic code with experimental results and use OSIRIS particle-in-cell code to verify the initial formation of the interaction region, showing good agreement between codes and experimental results. We find that the plasma conditions in the experiment are conducive to the lower-hybrid drift instability, yielding an increase in energy ΔE of ~ 264 keV for 242 MeV calcium ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Slowly braked, rotating neutron stars

A slowly braked, rotating neutron star is believed to be a star which rapidly rotates, has no nebula, is nonpulsing, and has a long initial braking time of ten thousand to a million years because of a low magnetic field. Such an object might be observable as an extended weak source of infrared or radio wave radiation due to the scattering of low-frequency strong-wave photons by accelerated electrons. If these objects exist abundantly in the Galaxy, they would act as sources of relatively low-energy cosmic rays. Pulsars (rapidly braked neutron stars) are shown to have difficulties in providing an adequate amount of cosmic-ray matter, making these new sources seem necessary. The possibility that the acceleration mechanism around a slowly braked star may be not a direct acceleration by the strong wave but an acceleration due to plasma turbulence excited by the strong wave is briefly explored. It is shown that white dwarfs may also be slowly braked stars with braking times longer than 3.15 million years.

Sato, H.↗

Magnetic reconnection in the era of exascale computing and multiscale experiments

Astrophysical plasmas have the remarkable ability to preserve magnetic topology, which inevitably gives rise to the accumulation of magnetic energy within stressed regions including current sheets. This stored energy is often released explosively through the process of magnetic reconnection, which produces a reconfiguration of the magnetic field, along with high-speed flows, thermal heating and nonthermal particle acceleration. Either collisional or kinetic dissipation mechanisms are required to overcome the topological constraints, both of which have been predicted by theory and validated with in situ spacecraft observations or laboratory experiments. However, major challenges remain in understanding magnetic reconnection in large systems, such as the solar corona, where the collisionality is weak and the kinetic scales are vanishingly small in comparison with macroscopic scales. The plasmoid instability or formation of multiple plasmoids in long, reconnecting current sheets is one possible multiscale solution for bridging this vast range of scales, and new laboratory experiments are poised to study these regimes. In conjunction with these efforts, we anticipate that the coming era of exascale computing, together with the next generation of observational capabilities, will enable new progress on a range of challenging problems, including the energy build-up and onset of reconnection, partially ionized regimes, the influence of magnetic turbulence and particle acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

Scattering theory in noncanonical phase space: A Drift-Kinetic collision operator for weakly collisional plasmas

After developing a scattering theory for grazing collisions in general noncanonical phase spaces, we introduce a guiding center collision operator in five-dimensional phase space designed for plasma regimes characterized by long wavelengths (relative to the Larmor radius), low frequencies (relative to the cyclotron frequency), and weak collisionality (where repeated Coulomb collisions induce cumulatively small changes in particle magnetic moment). The collision operator is fully determined by the noncanonical Hamiltonian structure of guiding center dynamics and exhibits a metriplectic structure, ensuring the conservation of particle number, momentum, energy, and interior Casimir invariants. It also satisfies an H-theorem, allowing for deviations from an equilibrium Maxwellian distribution due to the nontrivial kernel of the noncanonical guiding center Poisson tensor, spanned by the magnetic moment. We propose that this collision operator and its underlying mathematical structure may offer valuable insight into the study of turbulence, transport, and self-organizing phenomena in both laboratory and astrophysical plasmas.

Hamiltonian mechanics↗

Convection and magnetic fields in an atmosphere with constant temperature gradient. II - The supergranules and magnetic fields

The hydrodynamic properties of convective motions in a stratified atmosphere suggest that the supergranules in the solar photosphere extend to depths in excess of their horizontal diameters. The calculations presented here show that the weak fields which appear at the photosphere are as strong as, or stronger than, the horizontal fields anywhere in the supergranules below the surface. Hence the general 10 to 100 gauss azimuthal field of the sun, whose presence is indicated by the emerging bipolar magnetic regions, must lie below the base of the supergranules, i.e., below a depth probably of the order of 20,000 to 50,000 km. Presumably then the azimuthal field is generated by nonuniform rotation below the supergranule convection and does not appear in the supergranule layer because of the high rate of turbulent diffusion loss from that layer.

Parker, E. N.↗

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E×B drift to cause a non-diffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should in principle relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

intrinsic rotation↗

Fully-Explicit and Self-Consistent Algebraic Reynolds Stress Models

A fully-explicit, self-consistent algebraic expression for the Reynolds stress, which is the exact solution to the Reynolds stress transport equation in the 'weak equilibrium' limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the 'weak equilibrium' assumption treat the production-to-dissipation (P/epsilon) ratio implicitly, resulting in an effective viscosity that can be singular away from the equilibrium limit. In the present paper, the set of simultaneous algebraic Reynolds stress equations are solved in the full non-linear form and the eddy viscosity is found to be non-singular. Preliminary tests indicate that the model performs adequately, even for three dimensional mean flow cases. Due to the explicit and non-singular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.

Girimaji, Sharath S.↗

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗

Effects of Cavities and Protuberances on Transition over Hypersonic Vehicles

Surface protuberances and cavities on a hypersonic vehicle are known to cause several aerodynamic or aerothermodynamic issues. Most important of all, premature transition due to these surface irregularities can lead to a significant rise in surface heating. To help understand laminar-turbulent transition induced by protuberances or cavities on a Crew Exploration Vehicle (CEV) surface, high-fidelity numerical simulations are carried out for both types of trips on a CEV wind tunnel model. Due to the large bluntness, these surface irregularities reside in an accelerating subsonic boundary layer. For the Mach 6 wind tunnel conditions with a roughness Reynolds number Re(sub kk) of 800, it was found that a protuberance with a height to boundary layer thickness ratio of 0.73 leads to strong wake instability and spontaneous vortex shedding, while a cavity with identical geometry only causes a rather weak flow unsteadiness. The same cavity with a larger Reynolds number also leads to similar spontaneous vortex shedding and wake instability. The wake development and the formation of hairpin vortices for both protuberance and cavity were found to be qualitatively similar to that observed for an isolated hemisphere submerged in a subsonic, low speed flat-plate boundary layer. However, the shed vortices and their accompanying instability waves were found to be slightly stabilized downstream by the accelerating boundary layer along the CEV surface. Despite this stabilizing influence, it was found that the wake instability spreads substantially in both wall-normal and azimuthal directions as the flow is evolving towards a transitional state. Similarities and differences between the wake instability behind a protuberance and a cavity are investigated. Computations for the Mach 6 boundary layer over a slender cylindrical roughness element with a height to the boundary layer thickness of about 1.1 also shows spontaneous vortex shedding and strong wake instability. Comparisons of detailed flow structures associated with protuberances at subsonic and supersonic edge Mach numbers indicate distinctively different instability mechanisms.

Chang, Chau-Lyan↗

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E × B drift to cause a nondiffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine f L , the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, f L scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should, in principle, relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear $\hat{s}$ was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of f L with the electron collisionality. Still, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Vertical Structure and Ice Production Processes of Shallow Convective Postfrontal Clouds over the Southern Ocean in MARCUS. Part I: Observational Study

A study of the vertical structure of postfrontal shallow clouds in the marine boundary layer over the Southern Ocean is presented here. The central question of this two-part study regards cloud phase (liquid/ice) of precipitation, and the associated growth mechanisms. In this first part, data from the Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) field campaign are analyzed, starting with a 75-h case with continuous sea surface-based thermal instability, modest surface heat fluxes, an open-cellular mesoscale organization, and very few ice nucleating particles (INPs). The clouds are mostly precipitating and shallow (tops mostly around 2 km above sea level), with weak up- and downdrafts, and with cloud-top temperatures generally around –18° to –10°C. Further, the case study is extended to three other periods of postfrontal shallow clouds in MARCUS. While abundant supercooled liquid water is commonly present, an experimental cloud-phase algorithm classifies nearly two-thirds of clouds in the 0° to –5°C layer as containing ice (cloud ice, snow, or mixed phase), implying that much of the precipitation grows through cold-cloud processes. The best predictors of ice presence are cloud-top temperature, cloud depth, and INP concentration. Measures of convective activity and turbulence are found to be poor indicators of ice presence in the studied environment. The water-phase distribution in this cloud regime is explored through numerical simulations in Part II. Significance Statement: Climate models generally predict a lower albedo than observed over the Southern Ocean, and this is largely attributed to a lack of cloudiness, especially in the postfrontal cold sector of midlatitude cyclones. This in turn may be due to an excess of ice in these simulated clouds, resulting in rapid precipitation fallout and an overly brief cloud lifespan. The objective of this study is to examine whether shallow postfrontal clouds over the Southern Ocean are dominated by supercooled drops, or by snow and ice, using data collected by a U.S. Department of Energy Atmospheric Radiation Measurement Mobile Facility deployed aboard an Australian Antarctic supply vessel. We find that these clouds contain much supercooled liquid, even though cloud-top temperatures generally are around –18° to –8°C, and that about two-thirds of the clouds just above the freezing level contain ice. Much of the precipitation appears to grow through cold-cloud processes above the freezing level, rather than drizzle/rain. Updrafts and/or turbulence in convection or in cloud-top generating cells do not initiate much ice, compared to observations elsewhere in a similar temperature range. This may be attributable to the extremely low concentration of ice nucleating particles in this environment. Ultimately, the deepest clouds with the coldest cloud tops are most likely to be ice dominated.

54 ENVIRONMENTAL SCIENCES↗