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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2023

New insights on wind turbine wakes from large‐eddy simulation: Wake contraction, dual nature, and temperature effects

Abstract Large‐eddy simulation (LES) has been adopted to study wind turbine wakes because it can capture fine‐scale details of turbulent wind flows and interactions with wind turbines. Here, we use the LES version of the Weather Research and Forecasting (WRF) model with an actuator disk model to gain insights on several wake effects that have been traditionally difficult to measure. The first finding is that the wake has a “dual nature,” meaning that the wind speed deficit behaves differently from the added turbulent kinetic energy (TKE) and the two are not co‐located in space. For example, the wind speed deficit peaks at hub height and reaches the ground within 8D (D is the rotor diameter), but added TKE peaks near the rotor tip and generally remains aloft. Second, temperature changes near the ground are driven by the added TKE in the rotor area and by atmospheric stability. The combination of these two factors determines the sign and intensity of the vertical heat flux divergence below the rotor, with convergence and warming associated with stable conditions and weak divergence and modest cooling with unstable conditions. Third, wakes do not expand indefinitely, as suggested by similarity theory applied to the wind speed deficit, but eventually stop expanding and actually contract, at different rates depending on atmospheric stability. The implication of these findings is that, in order to study wakes, it is not sufficient to focus on wind speed deficit alone, because TKE is also important and yet behaves differently from the wind speed deficit.

17 WIND ENERGY↗

NASA NTRS1977

COMOC 2: Two-dimensional aerodynamics sequence, computer program user's guide

The COMOC finite element fluid mechanics computer program system is applicable to diverse problem classes. The two dimensional aerodynamics sequence was established for solution of the potential and/or viscous and turbulent flowfields associated with subsonic flight of elementary two dimensional isolated airfoils. The sequence is constituted of three specific flowfield options in COMOC for two dimensional flows. These include the potential flow option, the boundary layer option, and the parabolic Navier-Stokes option. By sequencing through these options, it is possible to computationally construct a weak-interaction model of the aerodynamic flowfield. This report is the user's guide to operation of COMOC for the aerodynamics sequence.

Manhardt, P. D.↗

NASA NTRS1986

Plasma wave observations at comet Giacobini-Zinner

The plasma wave instrument on the International Cometary Explorer (ICE) detected strong ion acoustic waves together with electromagnetic whistlers and low-level electron plasma oscillations when the spacecraft was within two million km of the nucleus of comet Giacobini-Zinner. As ICE approached the anticipated bow-shock location, electromagnetic and electrostatic wave levels increased significantly, but even amidst this turbulence, the wave instrument detected structures with familiar bow shock characteristics that were correlated with observations of localized electron heating phenomena. Just beyond the visible coma, high-amplitude broadband waves were detected accounting for the significant electron heating observed in this region. Near closest approach, broadband electrostatic noise was detected together with a changing pattern of weak electron plasma oscillations that yielded a density profile for the outer layers of the cold plasma tail. Near the tail axis, the plasma wave instrument also detected a nonuniform flux of dust impacts, and a preliminary profile of the Giacobini-Zinner dust distribution for micrometer-sized particles is presented.

Scarf, F. L.↗

DOE OSTI2022

Hydrokinetic energy conversion using flow induced oscillations of single-cylinder with large passive turbulence control

Various types of flow-induced oscillations (FIOs) have been implemented in development of marine hydrokinetic (MHK) energy converters. With passive turbulence control (PTC), energy harvesting starts at a flow speed of about 0.5 m/s. However, there is worldwide MHK energy available in even slower currents. In the present study, the effect of damping on FIO and power extraction is investigated for a converter with large turbulence stimulation (PTC) consisting of straight strips with a height of 15% of the cylinder diameter and placed symmetrically on the cylinder surface. The oscillating amplitude decreases, as the damping ratio increases, with unchanged sinusoidal pattern of the displacement time-history. The frequency ratio is also affected by damping especially in the VIV initial branch and transition region between VIV and galloping. An important flow characteristic of the large-PTC cylinder is that a recirculation region is formed behind the PTC, causing appreciable disturbance to the flow past the cylinder. Power can be harvested in the whole FIO range and the harnessed power maximum appears at the largest inflow velocity tested. However, the optimum of harnessing efficiency is located at the beginning of the VIV upper branch. The gap between VIV and galloping is bridged when large PTC is used, eliminating the drop in power and efficiency even at higher damping, which would be a weakness of regular-PTC cylinder for energy harvesting. Finally, the mechanism behind the variation of harnessing efficiency with inflow velocity and damping ratio is revealed, and the optimality criterion for the converter design is discussed.

16 TIDAL AND WAVE POWER↗

DOE OSTI2021

Gravito-turbulence and dynamo in poorly ionized protostellar discs – I. Zero-net-flux case

In their early stages, protoplanetary discs are sufficiently massive to undergo gravitational instability (GI). This instability is thought to be involved in mass accretion, planet formation via gas fragmentation, the generation of spiral density waves, and outbursts. A key and very recent area of research is the interaction between the GI and magnetic fields in young protoplanetary discs, in particular whether this instability is able to sustain a magnetic field via a dynamo. We conduct 3D, stratified shearing-box simulations using two independent codes, PLUTO and Athena++, to characterize the GI dynamo in poorly ionized protostellar discs subject to ambipolar diffusion. We find that the dynamo operates across a large range of ambipolar Elssaser number Am (which characterizes the strength of ambipolar diffusion) and is particularly strong in the regime Am = 10–100, with typical magnetic to thermal energy ratios of order unity. The dynamo is only weakly dependent on resolution (at least for Am ≲ 100), box size, and cooling law. The magnetic field is produced by the combination of differential rotation and large-scale vertical roll motions associated with spiral density waves. Our results have direct implications for the dynamo process in young protoplanetary discs and possibly some regions of active galactic nucleus discs.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1981

Prediction of the turbulent wake with second-order closure

A turbulence was envisioned whose energy containing scales would be Gaussian in the absence of inhomogeneity, gravity, etc. An equation was constructed for a function equivalent to the probability density, the second moment of which corresponded to the accepted modeled form of the Reynolds stress equation. The third moment equations obtained from this were simplified by the assumption of weak inhomogeneity. Calculations are presented with this model as well as interpretations of the results.

Taulbee, D. B.↗

NASA NTRS2015

Gravity Waves and Vertical Mixing in the Tropical Tropopause Layer during ATTREX

Vertical mixing in the Tropical Tropopause Layer is thought to be weak outside regionsof convective injection. Nevertheless, some studies indicate a rolefor nonconvective vertical mixing processes in maintaining tracer distributions.The TTL also has a broad spectrum of high amplitude gravity waves, which havebeen shown to play an important role in the water vapor distribution. These gravitywaves, a result of tropical convection, also havesubstantial shears associated with them (which may produce turbulent mixing). Recent aircraft experiments, specificallythe Airborne Tropical TRopopause EXperiment (ATTREX), have yielded an unprecendentedmeteorological dataset in the TTL, including not only the temperatures and horizontalwinds (which can be obtained from radiosondes), but high frequency vertical winds. These highfrequency measurements yield an opportunity to study turbulence in the nonconvectiveenvironment. The numerous vertical profiles, essentially producing curtains of high resolutionmeasurements across broad regions of the tropics, also yield an opportunity to obtain bothhorizontal and vertical information about medium scale gravity waves.This presentation examines three phases of ATTREX (fall, 2011, winter 2013, and winter 2014) inthe eastern and western Pacific. The focus will be on high amplitude wave motions, theirstructure, and turbulence.

Pfister, Leonhard↗

DOE OSTI2024

Advancing Thermal-Hydraulic Modeling Capabilities for Molten Salt Reactors in Pronghorn

This is an M2 report for the NEAMS program. This work focuses on advancing the modeling and simulation capabilities for Molten Salt Reactors (MSRs) within Pronghorn, addressing key thermal-hydraulic challenges that must be addressed for the safe and efficient operation of these advanced reactors. This report summarizes the development and validation of the following new models in Pronghorn: (i) two-equation RANS turbulence models, (ii) near-wall turbulence corrections due to curvature and roughness, and (iii) multidimensional two-phase drift-flux models. For each of these cases, along with the experimental validation, an application study of interest to MSRs is presented. In addition, this report also presents recent improvements in the domain-overlapping coupling between Pronghorn and SAM to address the cases of weakly compressible flow and flow reversal. Finally, an MSR action that greatly simplifies the setup of MSR models using NEAMS tools is described.

42 ENGINEERING↗

DOE OSTI2024

Staircase resiliency in a fluctuating cellular array

Inhomogeneous mixing by stationary convective cells set in a fixed array is a particularly simple route to layering. Layered profile structures, or staircases, have been observed in many systems, including drift-wave turbulence in magnetic confinement devices. The simplest type of staircase occurs in passive-scalar advection, due to the existence and interplay of two disparate timescales, the cell turn-over (τ H ), and the cell diffusion (τ D ) time. In this simple system, we study the resiliency of the staircase structure in the presence of global transverse shear and weak vortex scattering. The fixed cellular array is then generalized to a fluctuating vortex array in a series of numerical experiments. The focus is on regimes of low-modest effective Reynolds numbers, as found in magnetic fusion devices. By systematically perturbing the elements of the vortex array, we learn that staircases form and are resilient (although steps become less regular, due to cell mergers) over a broad range of Reynolds numbers. The criteria for resiliency are (a) τ D >>τ H and (b) a sufficiently high profile curvature (κ ≥ 1.5). We learn that scalar concentration travels along regions of shear, thus staircase barriers form first, and scalar concentration "homogenizes" in vortices later. The scattering of vortices induces a lower effective speed of scalar concentration front propagation. The paths are those of the least time. We observe that if background diffusion is kept fixed, the cell geometric properties can be used to derive an approximation for the effective diffusivity of the scalar. Furthermore, the effective diffusivity of the fluctuating vortex array does not deviate significantly from that of the fixed cellular array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2020

Tuning Neural Network Models for Improved Prediction of Boundary Layer Transition

Boundary layer transition can strongly impact flight vehicle performance as it influences surface skin friction and aerodynamic heating, making accurate transition prediction a key to designing next generation aircraft. Artificial neural networks (ANNs) have shown promise toward predicting laminar-turbulent transition based on linear stability correlations. The computational efficiency of ANNs and the substantially reduced user involvement in relation to direct computations based on the linear stability theory (LST) makes them an attractive methodology for integrating the LST based correlations in computational fluid dynamics codes. Tollmien-Schlichting (TS) waves correspond to the dominant transition mechanism in 2D or weakly 3D subsonic boundary layers, such as those encountered in general aviation applications. Improvements to neural network model accuracy in predicting the amplification rates of TS instability waves have been investigated by leveraging recent machine learning developments in conjunction with surrogate optimization techniques and via suitable augmentation of the data used to train the networks. The optimized models trained on the modified dataset reduced the average transition location errors on different airfoils at several flow conditions by 51% of the original manually-tuned network’s errors on the same flow cases. The actual transition locations were derived from the Langley Stability and Transition Analysis Code (LASTRAC).

Machine Learning↗

NASA NTRS1980

A numerical study of a marine subtropical stratus cloud layer and its stability

A model for numerical simulation of stratus cloud layers is constructed by combining a second-order closure, turbulent transfer model with a thermal radiative transfer model. The turbulent transfer model allows water vapor saturation. The combined turbulence-radiation model is applied to both a horizontally uniform one-dimensional case and a horizontally nonuniform two-dimensional case. In the latter, the dynamics of mesoscale circulations are also incorporated. Results of the two-dimensional simulation show that the layer cloud instability occurs where the sea surface temperature is high and the large-scale subsidence is weak. The simulated instability is analyzed in view of an instability criterion, the eddy kinetic energy budget, and evaporative cooling near the cloud top.

Moeng, C.-H.↗

NASA NTRS1984

Interaction of aerodynamic noise with laminar boundary layers in supersonic wind tunnels

The interaction between incoming aerodynamic noise and the supersonic laminar boundary layer is studied. The noise field is modeled as a Mach wave radiation field consisting of discrete waves emanating from coherent turbulent entities moving downstream within the supersonic turbulent boundary layer. The individual disturbances are likened to miniature sonic booms and the laminar boundary layer is staffed by the waves as the sources move downstream. The mean, autocorrelation, and power spectral density of the field are expressed in terms of the wave shapes and their average arrival rates. Some consideration is given to the possible appreciable thickness of the weak shock fronts. The emphasis in the interaction analysis is on the behavior of the shocklets in the noise field. The shocklets are shown to be focused by the laminar boundary layer in its outer region. Borrowing wave propagation terminology, this region is termed the caustic region. Using scaling laws from sonic boom work, focus factors at the caustic are estimated to vary from 2 to 6 for incoming shocklet strengths of 1 to .01 percent of the free stream pressure level. The situation regarding experimental evidence of the caustic region is reviewed.

Schopper, M. R.↗

NASA NTRS1993

Supersonic shock wave/vortex interaction

Although shock wave/vortex interaction is a basic and important fluid dynamics problem, very little research has been conducted on this topic. Therefore, a detailed experimental study of the interaction between a supersonic streamwise turbulent vortex and a shock wave was carried out at the Penn State Gas Dynamics Laboratory. A vortex is produced by replaceable swirl vanes located upstream of the throat of various converging-diverging nozzles. The supersonic vortex is then injected into either a coflowing supersonic stream or ambient air. The structure of the isolated vortex is investigated in a supersonic wind tunnel using miniature, fast-response, five-hole and total temperature probes and in a free jet using laser Doppler velocimetry. The cases tested have unit Reynolds numbers in excess of 25 million per meter, axial Mach numbers ranging from 2.5 to 4.0, and peak tangential Mach numbers from 0 (i.e., a pure jet) to about 0.7. The results show that the typical supersonic wake-like vortex consists of a non-isentropic, rotational core, where the reduced circulation distribution is self similar, and an outer isentropic, irrotational region. The vortex core is also a region of significant turbulent fluctuations. Radial profiles of turbulent kinetic energy and axial-tangential Reynolds stress are presented. The interactions between the vortex and both oblique and normal shock waves are investigated using nonintrusive optical diagnostics (i.e. schlieren, planar laser scattering, and laser Doppler velocimetry). Of the various types, two Mach 2.5 overexpanded-nozzle Mach disc interactions are examined in detail. Below a certain vortex strength, a 'weak' interaction exists in which the normal shock is perturbed locally into an unsteady 'bubble' shock near the vortex axis, but vortex breakdown (i.e., a stagnation point) does not occur. For stronger vortices, a random unsteady 'strong' interaction results that causes vortex breakdown. The vortex core reforms downstream of the rear stagnation point, and the reduced circulation distribution once again becomes self-similar in this region. A-new model of this interaction is proposed. Finally, a curve defining the approximate limits of supersonic vortex breakdown is presented.

Settles, G. S.↗

NASA NTRS1992

Mean winds and momentum fluxes over Jicamarca, Peru, during June and August 1987

Data from the mesosphere-stratosphere-troposphere radar at Jicamarca, Peru, together with other available data, are used to diagnose the mean structure of winds and gravity-wave momentum fluxes from the surface to 90 km during two-ten-day campaigns in June and August of 1987. In the stratosphere a layer of maximum eastward flow associated with the quasi-biennial oscillation (QBO) was seen to strengthen and descend rapidly from June to August, overlying persistent westward flow. A layer of enhanced signal return, suggestive of a turbulent layer, was observed just above the descending QBO eastward maximum. Notable zonal asymmetries were present during this transition and the local meridional circulation departed from zonal-mean QBO theory. A substantial northeastward momentum flux was found below 25 km, which may be related to topographic gravity waves excited by southeastward flow across the Andes. In the lower mesosphere a relatively weak 'second' mesopause semiannual oscillation is confirmed. In both the lower stratosphere and mesosphere, body forces were comparable in magnitude to inferred Coriolis torques.

Hitchman, Matthew H.↗

NASA NTRS1992

Lidar measurements using large liquid mirror telescopes

It is a well accepted practice in the discussion of lidars to compare systems by a performance factor given by the product of the transmitter power and receiver aperture. This form of the performance factor reflects the two key parameters that determine the amount of photons recorded from a given range bin; the amount of transmitted photons and the probability of collecting a scattered photon. We are interested in studying the propagation and breaking of gravity waves typically generated in the lower atmosphere. To study these waves requires high temporal and spatial resolution, as they have frequencies as high as 500 mHz and generate turbulent structures with vertical scales of less than 25 m. In the 80 - 100 km region of the atmosphere these measurements are most efficiently obtained by using resonance scattering from sodium atoms deposited by meteors. However, in the 30 - 80 km region, the most efficient scattering arises from molecules (Rayleigh scattering). Since the atmospheric scale height is nominally 7 km, this scattering becomes extremely weak above 60 km. Furthermore, present laser technology limits the power available for Rayleigh-scattering experiments to about 0.5 - 50 W. The transmitter for the lidar system under development at the University of Western Ontario uses a two-beam transmitter to simultaneously measure density perturbations and temperature from both Rayleigh and resonance scattering.

Sica, R. J.↗

NASA NTRS2002

Spin-Up Instability of a Levitated Molten Drop in MHD-Flow Transition to Turbulence

When an alternating magnetic field interacts with induced eddy currents in a conducting body, there will be a repulsive force between the body and the driving coil system generating the field. This repulsive force is the basis of electromagnetic levitation, which allows containerless processing of different materials. The eddy currents in the conducting body also generate Joule heating. Axial rotation of electromagnetically levitated objects is a common observation in levitation systems and often an undesirable side effect of such experiments on 1-g and -g. There have been recent efforts to use magnetic damping and suppress this tendency of body rotation. The first report of rotation in EML drops was attributed to a slight asymmetry of the shape and location of the levitation coils could change the axis and speed of rotation. Other theories of sample rotation include a frequency difference in the traveling electromagnetic waves and a phase difference in two different applied fields of the same frequency. All of these different mechanisms share the following characteristics: the torque is small, constant for constant field strength, and very weakly dependent on the sample's temperature and phase (solid or liquid). During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop is inferred from motion of impurities on the surface. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal levitated drop that undergo secondary flow instabilities. As rise in the fluid temperature rises, the viscosity falls and the internal flow accelerates and becomes oscillatory; and beyond a point in the experiments, the surface impurities exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. In this work, a background of these set of observations will be given followed by a presentation of our results on the digital particle tracking analysis that has been performed on a number of available videos. The analysis indicates that the levitated drop attains a constant rotational speed during the melting phase and formation of the co-rotating axi-symmetric laminar toroidal structures. However, the rate of axial rotation increases dramatically during the deformation of the toroidal structures anti their breakup into chaotic entities. This new data suggests an interaction between the flow inside the levitated molten drop and the driving coils in the experiments. Possible mechanisms for this interaction will be reviewed. The data will also be used to make an assessment of existing theories on droplet rotation.

Abedian, B.↗

NASA NTRS1990

Energy and spectral dynamics in forced compressible turbulence

Numerical simulations are used to study compressible turbulence with microscale Reynolds numbers up to 40 and rms Mach numbers M up to 0.9. The flows are randomly forced, with energy supplied to either the rotational or compressive components of kinetic energy, which is then transferred to internal energy through the pressure-dilatation interaction and viscous dissipation terms. Coupling between the two components of kinetic energy by the advection term is relatively weak, and most energy introduced to either component by the external force is transferred, without passing through the other component, to internal energy. A statistically quasi-equilibrium of kinetic energy is realized while internal energy increases steadily. The spectral form of the rotational component of velocity, which hardly depends on M, is very close to that for incompressible flow. On the other hand, the compressive component depends strongly on M, especially at large wave numbers.

Kida, Shigeo↗

NASA NTRS1996

Numerical simulation of jet aerodynamics using the three-dimensional Navier-Stokes code PAB3D

This report presents a unified method for subsonic and supersonic jet analysis using the three-dimensional Navier-Stokes code PAB3D. The Navier-Stokes code was used to obtain solutions for axisymmetric jets with on-design operating conditions at Mach numbers ranging from 0.6 to 3.0, supersonic jets containing weak shocks and Mach disks, and supersonic jets with nonaxisymmetric nozzle exit geometries. This report discusses computational methods, code implementation, computed results, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions. The Navier-Stokes method using the standard Jones-Launder two-equation kappa-epsilon turbulence model can accurately predict jet flow, and such predictions are made without any modification to the published constants for the turbulence model.

Pao, S. Paul↗