Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pressure gradient”

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

When does turbulence spreading matter?

Few, if any, of the many papers on turbulence spreading address the key question of how turbulence spreading actually affects the profile structure. In this work, we are using a reduced model to answer that question. Turbulence spreading is most relevant near regions where the profiles support a strong intensity gradient ∇I. One such case is at the edge of an L mode ischarge, near a source of turbulence [i.e., either a localized source of edge turbulence or an influx of turbulence from the scrape-off layer (SOL)]. Another is in “No Man’s Land” (NML), which connects the pedestal to the stiff core in H mode. In the case of L mode, without an edge intensity source, the turbulence intensity profile is nearly flat and spreading has a weak effect. An edge localized source increases the edge ∇I, which then drives inward spreading. Invasion of turbulence from the SOL to the edge softens the edge pressure gradient. In H mode, the strong shear suppression of pedestal turbulence necessarily forces a sharp ∇I in NML. This sharp∇I drives a significant flux of turbulence from the core to the pedestal, where it is ultimately dissipated by shearing. Counter-intuitively, the results indicate that spreading actually increases the pedestal height and width and hence the energy content in H mode. This suggests that models of the pedestal structure should include NML turbulence spreading effects. The relation of avalanches to spreading is studied. Spreading weakly affects the avalanche distribution, but the spatiotemporal correlation of intensity increases with spreading.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Study of Boundary Layer Development in a Two-Stage Low-Pressure Turbine

Experimental data from jet-engine tests have indicated that unsteady blade row (wake) interactions and separation can have a significant impact on the efficiency of turbine stages. The effects of these interactions can be intensified in low-pressure turbine stages because of the low Reynolds number operating environment. Measured turbine efficiencies at takeoff can be as much as two points higher than those at cruise conditions. Thus, during the last decade a significant amount of effort has been put into determining the effects of transition and turbulence on the performance of low pressure turbine stages. Experimental investigations have been performed, for example, by Hodson et al. and Halstead et al. These investigations have helped identify/clarify the roles that factors such as the Reynolds number, free stream turbulence intensity, pressure gradient and curvature have in the generation of losses. In parallel to the experimental investigations, there have been significant analytical efforts to improve the modeling of transition. Examples of such efforts include the works of Mayle and Gostelow et al. These newer models show promise of providing accurate transition predictions over a wide range of flow conditions, although they have yet to be implemented into the numerical flow analyses used by the turbine design community. Some recent computational investigations of interest include the works of Chernobrovkin and Lakshminarayana and Eulitz and Engel. The focus of the current effort has been to -use a viscous, unsteady quasi-three-dimensional Navier-Stokes analysis to study boundary layer development in a two-stage low-pressure turbine. A two-layer algebraic turbulence model, along with a natural transition model and a bubble transition model, have been used, The geometry used in the simulations has been the subject of extensive experiments. The predicted results have been compared with experimental data, including airfoil loadings and time-averaged/unsteady integral boundary layer quantities.

Dorney, Daniel J.↗
NASA NTRS1992

Measurements in a leading-edge separation bubble due to a simulated airfoil ice accretion

The separation bubble formed on an airfoil at low Reynolds number behind a simulated leading-edge glaze ice accretion is studied experimentally. Surface pressure and split hot-film measurements as well as flow visualization studies of the bubble reattachment point are reported. The simulated ice generates an adverse pressure gradient that causes a laminar separation bubble of the long bubble type to form. The boundary layer separates at a location on the ice accretion that is independent of angle of attack and reattaches at a downstream location 5-40 percent chord behind the leading edge, depending on the angle of attack. Velocity profiles show a large region of reverse flow that extends up from the airfoil surface as much as 2.5 percent chord. After reattachment, a thick distorted turbulent boundary layer exists. The separation bubble growth and reattachment are clearly seen in the plots of boundary-layer momentum thickness vs surface distance. Local minima and maxima in the boundary-layer momentum thickness development compare well with the shear layer transition point as indicated by the surface pressures and the reattachment point as measured from surface oil flow, respectively.

Bragg, M. B.↗
NASA NTRS1994

Crossflow instability control on a swept-wing: Preliminary studies

The pressure distribution on a swept wing causes the streamlines at the edge of the boundary layer to be curved. This pressure gradient normal to the external streamline creates a velocity component normal to the external streamline within the boundary layer which is referred to as the crossflow velocity. Because the crossflow velocity profile perpendicular to the wing surface has an inflection point, the profile is unstable. The stationary instability mode takes the form of crossflow vortices. Under these conditions, the boundary layer on the wing is extremely unstable and transition to turbulent flow takes place much closer to the leading edge of the wing than it would on an unswept wing. Higher skin friction drag is associated with turbulent flow, and so better aircraft performance could be obtained if the crossflow could be eliminated One method of controlling crossflow that is being investigated is boundary-layer suction. An extensive airfoil suction experiment in the 8 feet Transonic Pressure Tunnel (TPT) at NASA Langley Research Center will begin late in 1994. Because of the size, complexity, and expense associated with this test, a number of 'risk-reduction' tests are currently being conducted. The 20 x 28 in. Shear Flow Control Tunnel at NASA Langley is being used for some of these tests. Prior to the summer of 1994, a flat plate with a swept leading edge was installed in the 20 x 28 in. tunnel, with a displacement body mounted on the tunnel ceiling that created a pressure distribution on the plate similar to the pressure distribution on a swept wing. The flow over the plate was investigated during the summer of 1994 using a laser Doppler velocimeter (LDV) system. The LDV measurements indicated the possible presence of multiple disturbance modes, a rarely-seen phenomena since, in most tests, one disturbance mode dominates. The possible existence of multiple disturbance modes in the flat plate boundary layer, however, means that the flow in the 20 x 28 in. tunnel is of interest itself, and will be investigated more thoroughly in the future. With a view to these investigations, the boundary layer traverse mechanism in the 20 x 28 in. tunnel was modified to improve its performance, and strain gauges were mounted on the traverse in order to monitor its deflection during a test. Other preliminary work conducted in the 20 x 28 in. tunnel included the use of an infrared camera system. Previous work with this system showed that transition indeed could be detected, but the signal produced by the crossflow vortices was too weak to be detected. It was hoped that spraying the flat plate with naphthalene would augment the heat transfer associated with the crossflow vortices so that they would show up in a IR image; however, experiments showed that this would not work. Another set of tests was conducted in the 20 x 28 in. tunnel to determine the tripping requirements for a set of airfoil-shaped struts that will be used in the 8 feet TPT experiment. Since the Reynolds number associated with these struts is small, a laminar boundary layer would separate early, causing large fluctuations in the flow field. A turbulent boundary layer would remain attached further back, but tripping from laminar to turbulent flow at low Reynolds number is very difficult. However, trip strip configurations were found that should effectively trip the boundary layer at the required conditions. Currently underway is an investigation of the data acquisition requirements for the 8 feet TPT experiment, with the purpose of the finding the minimum amount of data needed to characterize sufficiently the swept-wing boundary layer. This study is being conducted using a numerically generated data set.

Bridges, David H.↗
NASA NTRS2023

Understanding the Deflection of the “Cartwheel CME”: Data Analysis and Modeling

We study the low corona evolution of the “Cartwheel” coronal mass ejection (CME; 2008 April 9) by- reconstructing its three-dimensional path and modeling it with magnetohydrodynamic simulations. This event exhibited a double deflection that has been reported and analyzed in previous works but whose underlying cause remained unclear. The Cartwheel CME traveled toward a coronal hole (CH) and against the magnetic gradients. Using a high-cadence, full-trajectory reconstruction, we accurately determine the location of the magnetic flux rope (MFR) and, consequently, the magnetic environment in which it is immersed. We find a pseudostreamer (PS) structure whose null point may be responsible for the complex evolution of the MFR at the initial phase. From the preeruptive magnetic field reconstruction, we estimate the dynamic forces acting on the MFR and provide a new physical insight into the motion exhibited by the 2008 April 9 event. By setting up a similar magnetic configuration in a 2.5D numerical simulation we are able to reproduce the observed behavior, confirming the importance of the PS null point. We find that the magnetic forces directed toward the null point cause the first deflection, directing the MFR toward the CH. Later, the magnetic pressure gradient of the CH produces the reversal motion of the MFR.

Solar prominences↗
NASA NTRS2022

Semispan Test Results of an Active Flow Control Enabled High-Lift Common Research Model in Landing Configuration

A 10%-scale semispan, Active Flow Control (AFC) enabled, simplified high-lift version of the Common Research Model (CRM-SHL-AFC) was tested in the 14- by 22-Foot Subsonic Tunnel at the NASA Langley Research Center. The main objective of the test was to develop an AFC system that can provide the necessary lift recovery on a simple-hinged flap high-lift system while minimizing its pneumatic power requirement. Three new types of AFC approaches were examined: Double-Row Sweeping Jets (DRSWJ), Alternating Pulsed Jets (APJ), and High Efficiency Low Power (HELP) actuators. The DRSWJ and the APJ actuators used two rows of unsteady jets, whereas the HELP actuators used an upstream row of sweeping jets combined with a downstream row of steady jets to overcome strong adverse pressure gradients. The test was conducted mostly at a freestream Mach number of 0.20. For exploration purposes, a limited number of runs were made at lower Mach numbers or using vortex generators (VGs). Minimal sensitivity to Mach number or VGs, for the cases evaluated, were observed. The AFC-induced lift coefficient increment was maintained over the AFC-off case for most flow-control cases examined. The CRM-SHL-AFC configuration equipped with HELP actuation was the only actuator configuration able to match or exceed the targeted lift performance of a reference conventional high-lift configuration. The presented aerodynamic data include lift, drag, and pitching moment coefficients as a function of angle of attack, with and without the Transonic Wall Interference Correction System (TWICS) method applied. Lift increments as a function of AFC pneumatic power usage (i.e., nozzle pressure ratio, mass flow, momentum coefficient, and power coefficient) are also presented at a lower angle of attack (α = 8.9°) and at maximum lift (α= 17.1°). At the two angles of attack, the surface pressure distributions and autospectral densities from the unsteady pressure transducers for the AFC-off case and the best HELP actuation case are compared.

high-lift↗
NASA NTRS1976

Modifications of the law of the wall and algebraic turbulence modelling for separated boundary layers

Various modifications of the conventional algebraic eddy viscosity turbulence model are investigated for application to separated flows. Friction velocity is defined in a way that avoids singular behavior at separation and reattachment but reverts to the conventional definition for flows with small pressure gradients. This leads to a modified law of the wall for separated flows. The effect on the calculated flow field of changes in the model that affect the eddy viscosity at various distances from the wall are determined by (1) switching from Prandtl's form to an inner layer formula due to Clauser at various distances from the wall, (2) varying the constant in the Van Driest damping factor, (3) using Clauser's inner layer formula all the way to the wall, and (4) applying a relaxation procedure in the evaluation of the constant in Clauser's inner layer formula. Numerical solutions of the compressible Navier-Stokes equations are used to determine the effects of the modifications. Experimental results from shock-induced separated flows at Mach numbers 2.93 and 8.45 are used for comparison. For these cases improved predictions of wall pressure distribution and positions of separation and reattachment are obtained from the relaxation version of the Clauser inner layer eddy viscosity formula.

Baldwin, B. S.↗
NASA NTRS1998

Incompressible Turbulent Wing-Body Junction Flow

The overall objective of this study is to contribute to the optimized design of fan bypass systems in advanced turbofan engines. Increasing the engine bypass ratios have provided a major boost in engine performance improvement over the last fifty years. An engine with high bypass ratio (11-16:1) such as the Advanced Ducted Propulsion (ADP) is being developed and is expected to provide an additional 25% improvement in overall efficiency over the early turbofans. Such significant improvements in overall efficiency would reduce the cost per seat mile, which is a major government and Industry challenge for the 21th century. The research is part of the Advanced Subsonic Technology (AST) program that involves a NASA, U.S. Industry and FAA partnership with the goal of a safe and highly productive global air transportation system. The immediate objective of the study is to perform numerical simulation of duct-strut interactions to elucidate the loss mechanisms associated with this configuration that is typical of advanced turbofan engines such as ADP. However, at present experimental data for a duct-strut configuration are not available. Thus, as a first step a wing-body junction flow would be studied and is the specific objective of the present study. At the outset it is to be recognized that while duct-strut interaction flow is similar to that of wing-body junction flows, there are some differences owing to the presence of a wall at both ends of the strut. Likewise, some differences are due to the sheared inflow (as opposed to a uniform inflow) velocity profile. It is however expected that some features of a wing-body junction flow would persist. Next, some of the salient aspects of the complex flow near a wing-body junction, as revealed by various studies reported in the literature will be reviewed. One of the principle characteristics of the juncture flow, is the presence of the mean flow components in a plane perpendicular to the direction of the oncoming free-stream flow. The lateral curvature of the wing/strat causes the oncoming turbulent layer to skew about am axis (x-axis) parallel to the plane (xz-plane) of the mean shear. This is the principle mechanism for the generation of secondary flow. Such skew-induced secondary flows are slow to be attenuated by Reynolds stresses. Additional contribution to the generation of secondary flow comes from anisotropies in Reynolds stresses. Upstream of the strut, the mean-vorticity is directed span wise (along the y-direction). The presence of secondary flow in the vicinity of the strut causes the vorticity to stretch around the obstacle in a horse-shoe shape, with each leg having a vorticity of the opposite sense. The blockage effect of the strut imposes a severe adverse pressure gradient on the oncoming turbulent shear layer, causing boundary layer separation ahead of the leading edge, resulting in a vortex that rolls up and flows downstream into the juncture region. The separation vortices trailing in the wake of the wing can alter the lift or drag characteristics of the surfaces downstream of the wing-body juncture. Likewise, on submarines, the wake flow behind the appendage can degrade the performance of the propeller located downstream. The complex nature of this flow is caused by the presence of all six components of Reynolds stresses. Devenport and Simpson report that in the vicinity of the horse-shoe vortex there is intense recirculation with turbulent stresses being much larger than those normally observed in turbulent flows. These features contribute to making this flow a challenge to predict numerically. Some of the past studies provide useful insights into this flow that would guide our numerical efforts. In measurements reported by Shabaka and Bradshaw, the eddy viscosity tensor is seen to be non-isotropic and has negative components in certain regions. In an effort to evaluate the closure assumptions of various turbulence models, Devenport and Simpson used their own extensive measurements in juncture flows around the nose of a wing-body junction. Measured values of mean-velocity and/or turbulence kinetic energy was used to predict the magnitude of the shear stress vector. Algebraic stress models performed the best followed by Cebeci-Smith eddy viscosity model. The flow is reported to be dominated by a pressure field produced by the wing and the velocity field generated by the horseshoe vortex that is wrapped around the junction between the wing and wall. Kubendran et al. conclude from an experimental study that the shape of leading edge of the wing as characterized by its slenderness ratio is a major factor in determining the flow fields in the juncture region. The more thinner the leading edge of the juncture, the weaker the horseshoe vortex is. Also, with a slender leading edge, the secondary flow in the juncture would be due mainly to the cross-stream gradients of Reynolds stresses rather than due to a lateral skewing of the shear layer.

Krishnamurthy, R.↗
DOE OSTI2025

Algorithmic Advancements for High-Order Self-Gravitating Hydrodynamics

Self-gravity plays a key role in the formation and evolution of many astronomical objects. Though gravity is often dominant at large scales, other forces (e.g., gas pressure gradients, radiation, and/or magnetic fields) often compete. It is therefore essential for numerical simulations to evaluate their interplay accurately and robustly. Hanawa & Mullen derived a 4th-order accurate finite volume scheme to solve the equations of self-gravitating hydrodynamics on a uniform Cartesian grid. In this work, we supply improvements to the algorithm that (1) mitigate spurious gravitational circulation and (2) greatly simplify the evaluation of the high order corrections. The proposed algorithm provides the gravitational acceleration (ρg) and the gravitational energy release (ρv · g) as source terms for the hydrodynamic equations, all while preserving conservation of linear momentum. Spurious heating and/or cooling associated with truncation error in the numerical evaluation of the gravitational energy release decreases in proportion to the fourth power of the cell width. We demonstrate fourth order convergence on smooth problems (e.g., 3D inclined sound wave propagation and 3D equilibria). An application test tracks the spherical collapse of a polytrope by an imposed, sudden decrease of the central gas pressure; a bounce and second collapse (associated with a spherical accretion shock) are robustly captured by the high order algorithm.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

Boundary-Layer Separation Control under Low-Pressure Turbine Airfoil Conditions using Glow-Discharge Plasma Actuators

Modem low-pressure turbines, in general, utilize highly loaded airfoils in an effort to improve efficiency and to lower the number of airfoils needed. Typically, the airfoil boundary layers are turbulent and fully attached at takeoff conditions, whereas a substantial fraction of the boundary layers on the airfoils may be transitional at cruise conditions due to the change of density with altitude. The strong adverse pressure gradients on the suction side of these airfoils can lead to boundary-layer separation at the latter low Reynolds number conditions. Large separation bubbles, particularly those which fail to reattach, cause a significant degradation of engine efficiency. A component efficiency drop of the order 2% may occur between takeoff and cruise conditions for large commercial transport engines and could be as large as 7% for smaller engines at higher altitude. An efficient means of of separation elimination/reduction is, therefore, crucial to improved turbine design. Because the large change in the Reynolds number from takeoff to cruise leads to a distinct change in the airfoil flow physics, a separation control strategy intended for cruise conditions will need to be carefully constructed so as to incur minimum impact/penalty at takeoff. A complicating factor, but also a potential advantage in the quest for an efficient strategy, is the intricate interplay between separation and transition for the situation at hand. Volino gives a comprehensive discussion of several recent studies on transition and separation under low-pressure-turbine conditions, among them one in the present facility. Transition may begin before or after separation, depending on the Reynolds number and other flow conditions. If the transition occurs early in the boundary layer then separation may be reduced or completely eliminated. Transition in the shear layer of a separation bubble can lead to rapid reattachment. This suggests using control mechanisms to trigger and enhance early transition. Gad-el-Hak provides a review of various techniques for flow control in general and Volino discusses recent studies on separation control under low-pressure-turbine conditions utilizing passive as well as active devices. As pointed out by Volino, passive devices optimized for separation control at low Reynolds numbers tend to increase losses at high Reynolds numbers, Active devices have the attractive feature that they can be utilized only in operational regimes where they are needed and when turned off would not affect the flow. The focus in the present paper is an experimental Separation is induced on a flat plate installed in a closed-circuit wind tunnel by a shaped insert on the opposite wall. The flow conditions represent flow over the suction surface of a modem low-pressure-turbine airfoil ('Pak-B'). The Reynolds number, based on wetted plate length and nominal exit velocity, is varied from 50,000 to 300,000, covering cruise to takeoff conditions. Low (0.2%) and high (2.5%) Gee-stream turbulence intensities are set using passive grids. A spanwise-oriented phased-plasma-array actuator, fabricated on a printed circuit board, is surface- flush-mounted upstream of the separation point and can provide forcing in a wide frequency range. Static surface pressure measurements and hot-wire anemometry of the base and controlled flows are performed and indicate that the glow-discharge plasma actuator is an effective device for separation control. of active separation control using glow discharge plasma actuators.

Hultgren, Lennart S.↗
NASA NTRS2002

Application of the generalized vertical coordinate ocean model for better representing satellite data

It is found that two adaptive parametric functions can be introduced into the basic ocean equations for utilizing the optimal or hybrid features of commonly used z-level, terrain- following, isopycnal, and pressure coordinates in numerical ocean models. The two parametric functions are formulated by combining three techniques: the arbitrary vertical coordinate system of Kasahara (1 974), the Jacobian pressure gradient formulation of Song (1 998), and a newly developed metric factor that permits both compressible (non-Boussinesq) and incompressible (Boussinesq) approximations. Based on the new formulation, an adaptive modeling strategy is proposed and a staggered finite volume method is designed to ensure conservation of important physical properties and numerical accuracy. Implementation of the combined techniques to SCRUM (Song and Haidvogel1994) shows that the adaptive modeling strategy can be applied to any existing ocean model without incurring computational expense or altering the original numerical schemes. Such a generalized coordinate model is expected to benefit diverse ocean modelers for easily choosing optimal vertical structures and sharing modeling resources based on a common model platform. Several representing oceanographic problems with different scales and characteristics, such as coastal canyons, basin-scale circulation, and global ocean circulation, are used to demonstrate the model's capability for multiple applications. New results show that the model is capable of simultaneously resolving both Boussinesq and non-Boussinesq, and both small- and large-scale processes well. This talk will focus on its applications of multiple satellite sensing data in eddy-resolving simulations of Asian Marginal Sea and Kurosio. Attention will be given to how Topex/Poseidon SSH, TRMM SST; and GRACE ocean bottom pressure can be correctly represented in a non- Boussinesq model.

numerical ocean model generalized coordinate syste↗
NASA NTRS1996

The Evolution of Finite Amplitude Wavetrains in Plane Channel Flow

We consider a viscous incompressible fluid flow driven between two parallel plates by a constant pressure gradient. The flow is at a finite Reynolds number, with an 0(l) disturbance in the form of a traveling wave. A phase equation approach is used to discuss the evolution of slowly varying fully nonlinear two dimensional wavetrains. We consider uniform wavetrains in detail, showing that the development of a wavenumber perturbation is governed by Burgers equation in most cases. The wavenumber perturbation theory, constructed using the phase equation approach for a uniform wavetrain, is shown to be distinct from an amplitude perturbation expansion about the periodic flow. In fact we show that the amplitude equation contains only linear terms and is simply the heat equation. We review, briefly, the well known dynamics of Burgers equation, which imply that both shock structures and finite time singularities of the wavenumber perturbation can occur with respect to the slow scales. Numerical computations have been performed to identify areas of the (wavenumber, Reynolds number, energy) neutral surface for which each of these possibilities can occur. We note that the evolution equations will breakdown under certain circumstances, in particular for a weakly nonlinear secondary flow. Finally we extend the theory to three dimensions and discuss the limit of a weak spanwise dependence for uniform wavetrains, showing that two functions are required to describe the evolution. These unknowns are a phase and a pressure function which satisfy a pair of linearly coupled partial differential equations. The results obtained from applying the same analysis to the fully three dimensional problem are included as an appendix.

Hewitt, R. E.↗
DOE OSTI2025

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗
NASA NTRS2023

Exploring the Langtry-Menter Transition Model for High Speed Applications Using FUN3D

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed using the FUN3D flow solver to explore the capabilities of the Langtry-Menter Shear-Stress Transport (LM-SST) transition model for predicting transition for aircraft inlet applications. Two geometries were simulated: a zeropressure-gradient flat plate and an axisymmetric cone exposed to hypersonic flow. In addition to the transition-sensitized LM-SST model investigations, simulations were run with the one-equation SpalartAllmaras (SA) and the two-equation Menter Shear-Stress Transport (SST-V) RANS models in fully turbulent mode to identify the natural RANS model transition behavior as a function of Mach number when executed in fully turbulent mode. The flat plate simulations showed that (1) the transition model was able to predict rapid transition at a freestream Mach number of 0.2, which is expected but (2) the predicted transition location moved downstream as the freestream Mach number was increased for the simulations that used the SST-V turbulence model. The latter is significant as it is usually assumed that one- and two-equation turbulence models will produce fully turbulent flow very near the boundary layer origin. The flat plate simulation freestream Mach number trend was confirmed with simulations using the Wind-US code, which also saw a similar trend when employing the SA turbulence model. For the axisymmetric cone simulations, the transition location was highly sensitive to the inflow turbulence levels. This is significant as the prediction of the transition location is crucial when trying to predict inlet performance, especially for hypersonic vehicle applications. It was also noted that the predicted transition location for the cone when using the SST-V turbulence model agreed well with the predicted transition location from the equivalent zero-pressure-gradient cold wall flat plate case.

Transition Model↗
DOE OSTI2020

Impact of equilibrium radial electric field on energy loss process after pedestal collapse

An impact of the equilibrium radial electric field on energy loss processes after pedestal collapse is numerically investigated using the BOUT++ framework. Using linear stability analysis, the resistive ballooning mode is shown to be stabilized by the radial shear of the equilibrium radial electric field. On the other hand, the energy loss level after the pedestal collapse increases if the equilibrium radial electric field is taken into account. The spatio-temporal and phase diagram analyses show that the equilibrium radial electric field partially cancels the fluctuation-driven toroidally axisymmetric radial electric field and weakens the E × B shearing rate after pedestal collapse, weakening the turbulence suppression by vortex shearing. Finally, the equilibrium radial electric field therefore increases turbulence intensity in nonlinear cyclic oscillations among pressure gradient, E × B shearing rate, and turbulence intensity, which gives rise to subsequent bursts of turbulent transport and increases the energy loss level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hypersonic Jets of Detonation Products in the Hydrodynamic Collapse of Macroscopic Voids

Localizing the energetic output from detonation waves has been a long-standing challenge in applied detonation physics. Here, energy localization is achieved via machined millimeter scale voids in pressed samples of PBX 9501, an HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane)-based plastic bonded explosive. A main mechanism of energy localization in these systems, the formation of hydrodynamic jets of dense product gases, is characterized experimentally using a semicylindrical geometry in witness plate impact experiments and streak imaging of the jet propagating into the air. The distance at which the jet is optimally developed is identified and the supersonic flow structure in the vicinity of this feature is explored using hydrocode simulations. This analysis found that most of the kinetic energy of the hydrodynamic jet arises from pressure gradients induced by geometrically mediated squeeze flow lateral to the direction of detonation propagation. This work presents a new development in the control of energetic output from detonation waves and applications to detonation wave shaping are discussed.

42 ENGINEERING↗
DOE OSTI2021

On the interplay between magnetic field and anisotropy in holographic QCD

We investigate the combined effects of anisotropy and a magnetic field in strongly interacting gauge theories by the gauge/gravity correspondence. Our main motivation is the quark-gluon plasma produced in off-central heavy-ion collisions which exhibits large anisotropy in pressure gradients as well as large external magnetic fields. We explore two different configurations, with the anisotropy either parallel or perpendicular to the magnetic field, focusing on the competition and interplay between the two. A detailed study of the RG flow in the ground state reveals a rich structure where depending on which of the two, anisotropy or magnetic field, is stronger, intermediate geometries with approximate AdS 4 × $\mathbb{R}$ and AdS 3 × $\mathbb{R}$ 2 factors arise. This competition is also manifest in the phase structure at finite temperature, specifically in the dependence of the chiral transition temperature on anisotropy and magnetic field, from which we infer the presence of inverse magnetic and anisotropic catalyses of the chiral condensate. Finally, we consider other salient observables in the theory, including the quark-antiquark potential, shear viscosity, entanglement entropy and the butterfly velocity. We demonstrate that they serve as good probes of the theory, in particular, distinguishing between the effects of the magnetic field and anisotropy in the ground and plasma states. We also find that the butterfly velocity, which codifies how fast information propagates in the plasma, exhibits a rich structure as a function of temperature, anisotropy and magnetic field, exceeding the conformal value in certain regimes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Coupled Thermo-Hydro-Mechanical Processes in Fractured Rocks: Some Past Scientific Highlights and Future Research Directions

Abstract Coupled thermo-hydro-mechanical (THM) processes in fractured rocks have been a topic of intense scientific research for more than 30 years. The present paper takes a look into the past and highlights some scientific advances which are of an unusual “out-of-the-box” nature, and then looks forward and discusses possible directions of future research in this interesting field of study. Concerning future research directions, we see a trend from a focus on coupled THM processes in single fractures or a few interacting fractures, to the study of coupled THM behavior in complex fracture network systems where the fractures act collectively giving rise to local stress concentration points and points of large pressure gradients. Three examples of future research directions are presented. First is an effort towards identifying characterizing parameters of a fracture network that play a direct controlling role in major coupled THM phenomena (such as induced seismicity and flow channeling), rather than parameters of stochastic distributions of fractures in the network. The second example of research direction is accounting for the heterogeneity and hierarchy of fractures in a fault or fracture zone which has been associated with major THM events in a number of geo-energy projects. The third example is at the opposite end of the first; here it is recognized that in some cases, the coupled THM processes in fractured rocks may be controlled dominantly by only a few key bridges. Identification, characterization, and evaluation of these key bridges should be one of the important research directions in the coming days.

58 GEOSCIENCES↗