Search NASA⌕ Search

SEARCH · Search NASA

Results for “ROTATIONAL FLOW”

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 55 records · Page 3

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aero-Thermal Characterization of Accelerating and Diffusing Passages Downstream of Rotating Detonation Combustors

Cycle benefits of rotating detonation engines show up to five percentage points of efficiency gain for low-pressure ratio engines. An optimal integration between the combustor and the turbine needs to be guaranteed to realize this potential gain. The rotating detonation combustor (RDC) exhausts transonic flow with shocks rotating at frequencies ranging from a few to tens of kilohertz depending on the number of present waves. Hence, the turbine design requires precise knowledge of the fluctuations and losses downstream of the combustor. This paper focuses on the quantification of fluctuations and losses for accelerating and diffusing passages. The analysis of the combustor is performed via reactive unsteady Reynolds Averaged Navier-Stokes (URANS) simulations. The unsteady RANS equations are solved via CFD++ from Metacomp with a one-step reaction mechanism for an H2-air mixture. The resolving of the boundary layer is achieved with a structured mesh of around 36 million cells. Inlet pressure of 10 bar and two different back pressures are applied to the combustor to model the interconnection with downstream turbines. Finally, we present and assess a methodology to reduce the computational time to model these passages ten times.

Braun, James↗

Structure Formation Through Magnetohydrodynamical Instabilities in Primordial Disks

The shear flow instabilities under the presence of magnetic fields in the primordial disk can greatly facilitate the formation of density structures that serve as seeds prior to the onset of the gravitational Jeans instability. We evaluate the effects of the Parker, magnetorotational and kinematic dynamo instabilities by comparing the properties of these instabilities. We calculate the mass spectra of coagulated density structures by the above mechanism in the radial direction for an axisymmetric magnetohydrodynamic (MHD) torus equilibrium and power density profile models. Our local three-dimensional MHD simulation indicates that the coupling of the Parker and magnetorotational instability creates spiral arms and gas blobs in an accretion disk, reinforcing the theory and model. Such a mechanism for the early structure formation may be tested in a laboratory. The recent progress in experiments involving shear flows in rotating tokamak, field reversed configuration (FRC) and laser plasmas may become a key element to advance in nonlinear studies.

Noguchi, Koichi↗

Figuring Out Gas & Galaxies in Enzo (FOGGIE). VI. The Circumgalactic Medium of L * Galaxies Is Supported in an Emergent, Nonhydrostatic Equilibrium

The circumgalactic medium (CGM) is often assumed to exist in or near hydrostatic equilibrium, with the regulation of accretion and the effects of feedback treated as perturbations to a stable balance between gravity and thermal pressure. We investigate global hydrostatic equilibrium in the CGM using four highly resolved L* galaxies from the Figuring Out Gas & Galaxies in Enzo (FOGGIE) project. The FOGGIE simulations were specifically targeted at fine spatial and mass resolution in the CGM (Δx ≲ 1 kpc h –1 and M ≃ 200M ⊙ ). We develop a new analysis framework that calculates the forces provided by thermal pressure gradients, turbulent pressure gradients, ram pressure gradients of large-scale radial bulk flows, centrifugal rotation, and gravity acting on the gas in the CGM. Thermal and turbulent pressure gradients vary strongly on scales of ≲5 kpc throughout the CGM. Thermal pressure gradients provide the main supporting force only beyond ~0.25R 200 , or ~50 kpc at z = 0. Within ~0.25R 200 , turbulent pressure gradients and rotational support provide stronger forces than thermal pressure. More generally, we find that global equilibrium models are neither appropriate nor predictive for the small scales probed by absorption line observations of the CGM. Local conditions generally cannot be derived by assuming a global equilibrium, but an emergent global equilibrium balancing radially inward and outward forces is obtained when averaging over the nonequilibrium local conditions on large scales in space and time. Approximate hydrostatic equilibrium holds only at large distances from galaxies, even when averaging out small-scale variations.

79 ASTRONOMY AND ASTROPHYSICS↗

Large-eddy simulations of turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines rotating at low tip speed ratios

Turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines (VAWTs) rotating at low tip speed ratios (TSRs) are studied numerically. The turbulent flows are simulated using the large-eddy simulation (LES) model, and the rotating turbine blades are modeled using the actuator line method. The helical VAWT has identical key parameters as the straight VAWT except for the 135° helical twist of the blades over the 0.3m vertical span. A set of LES runs are performed for two TSRs, 0.6 and 0.4, and the results are reported and analyzed. At these low TSRs, the wake behind the straight-bladed VAWT exhibits two-dimensional dominant flow motions (in the horizontal plane perpendicular to the straight blades) in the near-wake region that cause considerable spanwise expansion of the wake as it extends downstream. In contrast, the helical-bladed VAWT generates highly three-dimensional (3D) wake flow structures and upward/downward mean flow motions within the wake that cause the wake to expand mainly in the vertical direction. Turbulence statistical analyses also show that the 3D wake flow features induced by the helical blades accelerate the wake transition to turbulence and enhance the small-scale turbulent dissipation (as shown by the subgrid-scale turbulent dissipation in the LES), which leads to a more rapid decay of the wake turbulence intensity than that in the straight-bladed VAWT case at the same TSR. Finally, compared with the straight-bladed VAWT, the helical-bladed VAWT also exhibits much smaller temporal variations for the torque and power coefficients during the rotation cycle, which can be beneficial for wind power generation.

17 WIND ENERGY↗

Effect of outer-cylinder rotation on the radially heated Taylor–Couette flow

A Taylor–Couette setup with radial heating is considered where a Boussinesq fluid is sheared in the annular region between two concentric, independently rotating cylinders maintained at different temperatures. Linear stability analysis is performed to determine the Taylor number for the onset of instability. Two radius ratios corresponding to wide and thin gaps with several rotation rate ratios are considered. The rotation of the outer cylinder is found to have a general stabilizing effect on the stability threshold as compared to pure inner-cylinder rotation, with a few exceptions. The radial heating sets up an axial flow which breaks the reflection symmetry of isothermal Taylor–Couette flow in the axial coordinate. This symmetry breaking separates linear stability thresholds, and we find the fastest growing modes with both positive and negative azimuthal numbers for different parameters. Another important finding of the current study is the discovery of unstable modes in the Rayleigh-stable regime. Furthermore, closed disconnected neutral curves (CDNCs) are observed for both wide and thin gaps which can separate from or merge into open neutral–stability curves. Alternatively, CDNCs can also morph into open neutral stability curves as the rotation rate ratio is changed. CDNCs are observed to be sensitive to changes in control parameters, and their appearance/disappearance is shown to induce discontinuous jumps in the critical Taylor number. Finally, for both wide and thin gaps, the fastest-growing modes found in the pure corotation case are shown to have their origins in the instability islands at smaller values of rotation rate ratios.

97 MATHEMATICS AND COMPUTING↗

Flow Development in Radial Plane of Rotating Detonation Engine Integrated with Aerospike

In this study, axial and radial velocity components are measured downstream of a rotating detonation engine (RDE) integrated with an aerospike by using particle image velocimetry (PIV) at 100 kHz. The RDE is operated at high pressures obtained by restricting the RDE exit with a converging nozzle, which also directs the flow radially toward the aerospike. Reactants, methane fuel, and enriched air (67% oxygen and 33% nitrogen) are supplied from separate plenums at ambient temperature. PIV data are presented to investigate how pressure ratios within the RDE channel affect the flow evolution across the aerospike. The RDE is shown to operate stably and consistently in a single wave mode for all six test runs. The study highlights large temporal and spatial variations in both axial and radial flow velocities at the nozzle throat that persist downstream across the aerospike. Each point in the flowfield oscillates at the frequency recorded inside the RDE channel. Overall, the temporally and spatially varying aerospike flowfield is far from the ideal case of a uniform flow at its exit, and it points toward the need to condition the flow within the RDE channel to produce more uniform conditions at the nozzle throat.

Engineering↗

Review of rotating wing dynamic stall: Experiments and flow control

Dynamic stall has been a technical challenge and a fluid dynamical subject of interest for more than fifty years; but in the last decade significant advances have been made in the understanding, prediction, modeling, and control of dynamic stall on rotors. This paper provides a summary of the state of the art of dynamic stall experiments and future directions in the understanding of dynamic stall on rotors. Experimental data sets are discussed, as well the direction of future research for control of dynamic stall. Coordinated testing between airfoils and rotating blades, as well as close integration between computational and experimental studies were found to be productive approaches. Advanced analysis methods, including statistical methods, modal representations, and artificial intelligence methods have led to significant advances in the understanding of dynamic stall. Investigations of dynamic stall control devices have allowed many useful targeted investigations of the transition to separated flow, but have not yet resulted in a commercially implemented device.

42 ENGINEERING↗

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗

Large-Eddy Simulation of a wind turbine using a Filtered Actuator Line Model

When dealing with multirotor devices such as quadcopters or wind farms, the cost of blade-resolved large-eddy simulation (LES) becomes prohibitive. Combining LES with a family of lower-fidelity models, called actuator line models (ALMs), has grown in popularity in the past decade. ALM replaces full blade resolution with an array of actuator points or lines parameterized by aerodynamic lift/drag polar plots along the blades. Body forces computed based on these actuator points are then projected onto the LES flow mesh, mimicking the effect of rotating blades on the flow. However, the optimal projection radius and the associated LES grid size is often too restrictive for multirotor simulations. Recently, a new tip-correction-based filtered ALM (F-ALM) was proposed by Martinez-Tossas and Meneveau (2019), which allows coarser-than-optimal grids by avoiding the associated overprediction of thrust. In this work, F-ALM is implemented into a high-order, in-house LES code to simulate National Renewable Energy Laboratory Phase VI wind turbine. It is then followed by a comparison between the baseline ALM and the newly implemented F-ALM in terms of instantaneous and time-averaged flow fields and blade loads, revealing the advantages of F-ALM in preventing the overprediction of power on coarse grids. Finally, this encourages accurate and affordable simulations of multirotor devices in the future.

17 WIND ENERGY↗

Effects of Part-to-Part Flow Variations on Overall Effectiveness and Life of Rotating Turbine Blades

Abstract As firing temperatures in gas turbine engines continue to increase to achieve high efficiencies, components in the main gas path must be protected with cooling flows to ensure lifing targets are met. Manufacturing variations, however, influence the performance and life characteristics of components with the same nominal design. This study presents blade flow and overall cooling effectiveness measurements for nine true-scale, aero engine turbine blades with realistic manufacturing variations. Flow measurements were made through each blade at a fixed pressure ratio to determine flow variability between holes and between blades. Infrared thermography was used to capture spatially-resolved temperature measurements reported as overall effectiveness on the same nine blades under high-speed rotating conditions at the Steady Thermal Aero Research Turbine Laboratory. Thermal performance was correlated with blade flow performance indicating substantial blade-to-blade variations resulting from manufacturing differences. Measurements also indicated wide variations in cooling jet trajectories as well as overall cooling effectiveness. Finally, the observed blade-to-blade variations in effectiveness were scaled to engine conditions with lifing estimates showing some blades would be expected to last only half as long as others due to manufacturing variability.

Engineering↗

Polarization and Vorticity in the Quark–Gluon Plasma

The quark–gluon plasma (QGP) produced by collisions between ultrarelativistic heavy nuclei is well described in the language of hydrodynamics. Noncentral collisions are characterized by very large angular momentum, which in a fluid system manifests as flow vorticity. This rotational structure can lead to a spin polarization of the hadrons that eventually emerge from the plasma, and thus these collisions provide experimental access to flow substructure at unprecedented detail. Recently, the first observations of Λ hyperon polarization along the direction of collisional angular momentum were reported. These measurements are in broad agreement with hydrodynamic and transport-based calculations and reveal that the QGP is the most vortical fluid ever observed. However, there remain important tensions between theory and observation that might be fundamental in nature. In the relatively mature field of heavy-ion physics, the discovery of global hyperon polarization and 3D simulations of the collision have opened an entirely new direction of research. We discuss the current status of this rapidly developing area and directions for future research.

heavy-ion collisions↗

Laser Measurements and Modeling of Shock Tunnel Freestream Velocity and Multispecies Thermal Nonequilibrium

Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N 2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N 2 and O 2 , which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. Furthermore, this extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N 2 exhibits the most nonequilibrium, followed by O 2 and NO. The CFD model captures this trend, although it consistently overpredicts N 2 and O 2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.

Aerodynamics↗

Electrolyzers in focus: advances in CO 2 electrolyzer designs

Electrochemical CO 2 reduction (ECR) remains a viable method to reintegrate anthropogenic CO 2 into current energy infrastructures through its conversion into commodity chemicals. To facilitate the integration of ECR, electrochemical devices called electrolyzers must be implemented to overcome the inherent limitations that exist in current ECR experiments, namely kinetics and mass transport. In this review, we outline the current and advancing designs in ECR electrolyzers, with a focus on the following five electrochemical devices: membrane electrode assemblies (MEA), flow cell (FC), rotating disk electrode (RDE), rotating ring-disk electrode (RRDE), and rotating cylinder electrode (RCE). We highlight the tunable components of each electrolyzer with a forward outlook on the optimization and relevance of electrolyzer designs in upcoming ECR applications.

CO2 reduction↗

Taylor–Couette flow for astrophysical purposes

A concise review is given of astrophysically motivated experimental and theoretical research on Taylor–Couette flow. The flows of interest rotate differentially with the inner cylinder faster than the outer, but are linearly stable against Rayleigh’s inviscid centrifugal instability. At shear Reynolds numbers as large as 10 6 , hydrodynamic flows of this type (quasi-Keplerian) appear to be nonlinearly stable: no turbulence is seen that cannot be attributed to interaction with the axial boundaries, rather than the radial shear itself. Direct numerical simulations agree, although they cannot yet reach such high Reynolds numbers. This result indicates that accretion-disc turbulence is not purely hydrodynamic in origin, at least insofar as it is driven by radial shear. Theory, however, predicts linear magnetohydrodynamic (MHD) instabilities in astrophysical discs: in particular, the standard magnetorotational instability (SMRI). MHD Taylor–Couette experiments aimed at SMRI are challenged by the low magnetic Prandtl numbers of liquid metals. High fluid Reynolds numbers and careful control of the axial boundaries are required. The quest for laboratory SMRI has been rewarded with the discovery of some interesting inductionless cousins of SMRI, and with the recently reported success in demonstrating SMRI itself using conducting axial boundaries. Some outstanding questions and near-future prospects are discussed, especially in connection with astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 × 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E × B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E × B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimentally validated high-fidelity simulations of a liquid jet in supersonic crossflow

Here, utilizing recent advancements in computational schemes for compressible, multiphase flows, this work features a parametric study of a pure liquid jet in supersonic crossflow that involves simulating the atomization process for four values of momentum-flux ratio. These simulations are validated against experimental results measured with high-speed X-ray imaging, which confirm the accuracy of the numerical approach. Also, the effect of numerical resolution on some flow behavior is investigated, revealing convergence of the jet shape and surface instability wavelength. Analysis of the resulting sprays includes statistical descriptions of the liquid distribution, liquid structures created through breakup, interfacial instabilities, and dominant flow features. As the flowrate increases, the spray penetrates further, it becomes more disperse, and less liquid impacts the wall, but the droplet size distribution changes little. The wavelength of instabilities on the windward side of the jet diverges from measured trends in subsonic crossflows. In a visualization of the time-averaged flow, counter-rotating vortices are observed along the jet core and in the wake, affecting the process of primary atomization and early droplet trajectories.

42 ENGINEERING↗