Search NASA⌕ Search

SEARCH · Search NASA

Results for “Vortex flows”

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 1,243 records · Page 69

Graphical representation of three-dimensional vortical flows by means of helicity density and normalized helicity

The helicity density and normalized helicity are introduced as important tools for the graphical representation of three-dimensional flowfields that contain concentrated vortices. The use of these two quantities filters out the graphical representation of the flowfield regions of low vorticity, as well as regions of high vorticity but low speed where the angle between the velocity and vorticity vectors is large (such as in the boundary layer). Their use permits the researcher to identify and accentuate the concentrated vortices, differentiate between primary and secondary vortices and mark their separation and reattachment lines. The method allows locating the singular points in type I primary separations and tracing the vortex-core streamlines that emanate from them. It also allows locating the vortex-core axes in flows of other types.

Seginer, Arnan↗

High resolution numerical simulations of methane pool fires using adaptive mesh refinement

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. Here, in this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

42 ENGINEERING↗

Aerodynamic Inner Workings of Circumferential Grooves in a Transonic Axial Compressor

The current paper reports on investigations of the fundamental flow mechanisms of circumferential grooves applied to a transonic axial compressor. Experimental results show that the compressor stall margin is significantly improved with the current set of circumferential grooves. The primary focus of the current investigation is to advance understanding of basic flow mechanics behind the observed improvement of stall margin. Experimental data and numerical simulations of a circumferential groove were analyzed in detail to unlock the inner workings of the circumferential grooves in the current transonic compressor rotor. A short length scale stall inception occurs when a large flow blockage is built on the pressure side of the blade near the leading edge and incoming flow spills over to the adjacent blade passage due to this blockage. The current study reveals that a large portion of this blockage is created by the tip clearance flow originating from 20% to 50% chord of the blade from the leading edge. Tip clearance flows originating from the leading edge up to 20% chord form a tip clearance core vortex and this tip clearance core vortex travels radially inward. The tip clearance flows originating from 20% to 50% chord travels over this tip clearance core vortex and reaches to the pressure side. This part of tip clearance flow is of low momentum as it is coming from the casing boundary layer and the blade suction surface boundary layer. The circumferential grooves disturb this part of the tip clearance flow close to the casing. Consequently the buildup of the induced vortex and the blockage near the pressure side of the passage is reduced. This is the main mechanism of the circumferential grooves that delays the formation of blockage near the pressure side of the passage and delays the onset of short length scale stall inception. The primary effect of the circumferential grooves is preventing local blockage near the pressure side of the blade leading edge that directly determines flow spillage around the leading edge. The circumferential grooves do not necessarily reduce the over all blockage built up at the rotor tip section.

Hah, Chunill↗

Saturation of Vortex Rings Ejected from Shock-Accelerated Interfaces

Structures evoking vortex rings can be discerned in shock-accelerated flows ranging from astrophysics to inertial confinement fusion. By constructing an analogy between vortex rings produced in conventional propulsion systems and rings generated by a shock impinging upon a high-aspect-ratio protrusion along a material interface, we extend classical, constant-density vortex-ring theory to compressible multifluid flows. Here, we further demonstrate saturation of such vortex rings as the protrusion aspect ratio is increased, thus explaining morphological differences observed in practice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of the Use of Second Generation Wavelets in the Coherent Vortex Simulation Approach

The objective of this study is to investigate the use of the second generation bi-orthogonal wavelet transform for the field decomposition in the Coherent Vortex Simulation of turbulent flows. The performances of the bi-orthogonal second generation wavelet transform and the orthogonal wavelet transform using Daubechies wavelets with the same number of vanishing moments are compared in a priori tests using a spectral direct numerical simulation (DNS) database of isotropic turbulence fields: 256(exp 3) and 512(exp 3) DNS of forced homogeneous turbulence (Re(sub lambda) = 168) and 256(exp 3) and 512(exp 3) DNS of decaying homogeneous turbulence (Re(sub lambda) = 55). It is found that bi-orthogonal second generation wavelets can be used for coherent vortex extraction. The results of a priori tests indicate that second generation wavelets have better compression and the residual field is closer to Gaussian. However, it was found that the use of second generation wavelets results in an integral length scale for the incoherent part that is larger than that derived from orthogonal wavelets. A way of dealing with this difficulty is suggested.

Goldstein, D. E.↗

Vortex bursting and associated twist dynamics on helical vortex tubes and vortex rings

The interaction of opposite-signed twist waves on vortex tubes can lead to vortex bursting, a process where the core expands into a double ring-like structure with strong swirling flows. Previous works have studied vortex bursting on rectilinear vortices by axially perturbing the initial core size to generate the twist waves, and observed largely axisymmetric bursting dynamics. In this work, we numerically study bursting on vortical structures with curved centrelines, analysing the interaction between the centreline dynamics, twist wave generation and propagation, and vortex bursting. We focus on axially perturbed helical vortex tubes with small radius-to-pitch ratios up to$0.0625$, as well as vortex rings with a large radius-to-core size ratio$10$, both at a circulation-based Reynolds number$5000$. The results show that though the initial twist wave propagation speeds are relatively unaffected by the curvature and torsion of the centreline, the bursting process is altered significantly compared with rectilinear vortices. The self-induced rotation of the centreline of the helical tube induces a non-axisymmetric distortion of the bursting structure, which rapidly breaks up the vortex core into small-scale helical structures. A similar destabilization of the bursting structure also occurs on vortex rings. The enstrophy increase and accelerated energy decay associated with bursting are predominantly determined by the twist wave strength, rather than the curvature and torsion of the centreline. Combined, our findings imply that bursting could play an important role in transferring and dissipating energy of vortical structures in wakes, and turbulent flows in general.

Mechanics↗

Ignition and structure of a laminar diffusion flame in the field of a vortex

The distortion of flames in flows with vortical motion is examined via asymptotic analysis and numerical simulation. The model consists of a constant density, one step, irreversible Arrhenius reaction between initially unmixed species occupying adjacent half-planes which are then allowed to mix and react in the presence of a vortex. The evolution in time of the temperature and mass fraction fields is followed. Emphasis is placed on the ignition time and location as a function of vortex Reynolds number and initial temperature differences of the reacting species. The study brings out the influence of the vortex on the chemical reaction. In all phases, good agreement is observed between asymptotic analysis and the full numerical solution of the model equations.

Macaraeg, Michele G.↗

Geometric Perturbations and Asymmetric vortex shedding About Slender Pointed Bodies

The flow about slender, pointed bodies can be characterized by different states with angle of attack. At moderate-to-high angles of attack (alpha approximately equals 40deg), a steady, asymmetric vortex pattern develops along the body, leading to a net lateral force. At higher angles of attack (alpha approximately equals 60deg), the aft-end of the body develops an unsteady von Karman shedding. As the angle of attack approaches 90deg, the entire body length exhibits a time-dependent vortex shedding pattern. The current work uses three-dimensional, thin-layer Navier-Stokes simulations to investigate the physical mechanisms of asymmetric vortex shedding at alpha = 40deg and alpha = 60deg. The development of an asymmetric vortex pattern via a convective instability mechanism is investigated using tip bumps, surface roughness, and tip curvature. It's found that surface roughness simulations can incite an asymmetric vortex state at alpha = 60deg which is consistent with the application of a tip bumps, and the experimentally observed flowfield. The unsteady von Karman vortex shedding on the aft portion of the body is also well resolved. The use of surface roughness did not incite a flow asymmetry at alpha = 40deg, and it was necessary to simulate tip curvature at this angle of attack in order to generate an asymmetric vortex state.

Murman, Scott M.↗

Vortex/surface interaction

This paper considers the interaction of a vortex generated upstream in a flow field with a downstream aerodynamic surface that possesses a large chord. The flow is assumed to be steady, incompressible, inviscid and irrotational, and the surface to be semiinfinite. The vortex is considered to be a straight vortex filament. To lowest order the problem is modeled using potential theory, where the 3D Laplace's equation for the velocity potential on the surface is solved exactly. The closed-form equation for pressure distribution obtained from this theory is found to have a square root singularity at the leading-edge. It also converges, as x goes to infinity, to the solution of the 2D point-vortex/infinite plane problem. The pressure coefficient presents an anti-symmetric behavior, near the leading-edge and a symmetric behavior as x goes to infinity.

Bodstein, G. C. R.↗

Geometric Perturbations and Asymmetric Vortex Shedding about Slender Pointed Bodies

The flow about slender, pointed bodies can be characterized by different states with angle of attack. At moderate-to-high angles of attack (alpha approx. = 40 deg.), a steady, asymmetric vortex pattern develops along the body, leading to a net lateral force. At higher angles of attack (alpha approx. = 60 deg.), the aft-end of the body develops an unsteady von Karman shedding. As the angle of attack approaches 90 deg., the entire body length exhibits a time-dependent vortex shedding pattern. The current work uses three-dimensional, thin-layer Navier-Stokes simulations to investigate the physical mechanisms of asymmetric vortex shedding at alpha = 40 deg. and alpha = 60 deg. The development of an asymmetric vortex pattern via a convective instability mechanism is investigated using tip bumps, surface roughness, and tip curvature. It's found that surface roughness simulations can incite an asymmetric vortex state at alpha = 60 deg. which is consistent with the application of a tip bump, and the experimentally observed flowfield. The unsteady von Karman vortex shedding on the aft. portion of the body is also well resolved. The use of surface roughness did not incite a flow asymmetry at alpha = 40 deg., and it was necessary to simulate tip curvature at this angle of attack in order to generate an asymmetric vortex state.

Murman, Scott M.↗

Scattering of sound waves by a compressible vortex

Scattering of plane sound waves by a compressible vortex is investigated by direct computation of the two-dimensional Navier-Stokes equations. Nonreflecting boundary conditions are utilized, and their accuracy is established by comparing results on different sized domains. Scattered waves are directly measured from the computations. The resulting amplitude and directivity pattern of the scattered waves is discussed, and compared to various theoretical predictions. For compact vortices (zero circulation), the scattered waves directly computed are in good agreement with predictions based on an acoustic analogy. Strong scattering at about + or - 30 degrees from the direction of incident wave propagation is observed. Back scattering is an order of magnitude smaller than forward scattering. For vortices with finite circulation refraction of the sound by the mean flow field outside the vortex core is found to be important in determining the amplitude and directivity of the scattered wave field.

Colonius, Tim↗

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.↗

Shock-vortex interaction over a 65-degree delta wing in transonic flow

Transonic flow over a 65-deg swept-back, sharp-edged, cropped delta wing is investigated computationally using the time-accurate solution of the unsteady, compressible, full Navier-Stokes equations with an implicit, upwind, flux-difference splitting, finite-volume scheme. Coarse and fine O-H grids are used to obtain the solution. The grid consists of 125 x 85 x 84 points in the wrap-around, normal, and axial directions, respectively. The results are presented for an angle of attack of 20 deg Mach number of 0.85 and Reynolds number of 3.23 x 10 exp 6. With the fine grid, the results show that a system of shocks has been captured over the upper wing surface and that the leading-edge vortex core experiences an unsteady supersonic vortex breakdown after passing through a spanwise shock near the wing trailing edge. The computed results at a certain time are in good agreement with the experimental data. Topological aspects of the vortex breakdown flowfield are also presented and discussed.

Kandil, Osama A.↗

Reduced-order modeling on a near-term quantum computer

Quantum computing is an advancing area of research in which computer hardware and algorithms are developed to take advantage of quantum mechanical phenomena. In recent studies, quantum algorithms have shown promise in solving linear systems of equations as well as systems of linear ordinary differential equations (ODEs) and partial differential equations (PDEs). Reducedorder modeling (ROM) algorithms for studying fluid dynamics have shown success in identifying linear operators that can describe flowfields, where dynamic mode decomposition (DMD) is a particularly useful method in which a linear operator is identified from data. In this work, DMD is reformulated as an optimization problem to propagate the state of the linearized dynamical system on a quantum computer. This reformulation was chosen as a means of facilitating implementation on a near-term quantum computer. Quadratic unconstrained binary optimization (QUBO), a technique for optimizing quadratic polynomials in binary variables, allows for quantum annealing algorithms to be applied. A quantum circuit model (quantum approximation optimization algorithm, QAOA) is utilized to obtain predictions of the state trajectories. Results are shown for the quantum-ROM predictions for flow over a 2D cylinder at Re = 220 and flow over a NACA0009 airfoil at Re = 500 and α = 15°. The quantum-ROM predictions are found to depend on the number of bits utilized for a fixed point representation and the truncation level of the DMD model. Comparisons with DMD predictions from a classical computer algorithm are made, as well as an analysis of the computational complexity and prospects for future, more fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗

Flutter Analysis of the Transonic Truss-Braced Wing Aircraft Using Transonic Correction

This paper describes a flutter analysis method for the Transonic Truss-Braced Wing aircraft using a vortex lattice method coupled to an unsteady transonic correction method to account for unsteady aerodynamics in transonic flow. A steady-state vortex-lattice model of the Truss-Braced Wing aircraft is developed using vortex lattice code VSPAERO. A transonic and viscous flow correction method is implemented in the VSPAERO model to account for steady-state transonic and viscous flow effects using transonic small disturbance code TSFOIL coupled to an in-house integral boundary layer code. In addition, a wing-strut interference correction method is developed to account for the transonic interference aerodynamics in the strut juncture region using high-fidelity CFD code FUN3D. A structural dynamic finite-element model of the Truss-Braced Wing aircraft is developed using BEAM3D in-house finite-element code and is coupled to the VSPAERO. The BEAM3D model includes a geometric nonlinearity due to the tension in the strut which causes a deflection-dependent nonlinear stiffness. An unsteady transonic correction method is developed to better capture the unsteady aerodynamics in transonic flow. The unsteady transonic correction method makes use of the Theodorsen’s theory to account for the amplitude and phase shift of the unsteady lift coefficient in transonic flow. A preliminary flutter analysis of the Truss-Braced Wing aircraft is conduct to illustrate the unsteady transonic correction approach.

Nhan Nguyen↗

Describing the Mechanism of Instability Suppression Using a Central Pilot Flame With Coupled Experiments and Simulations

Pilot flames are commonly used to extend combustor operability limits and suppress combustion oscillations in low-emissions gas turbines. Combustion oscillations, a coupling between heat release rate oscillations and combustor acoustics, can arise at the operability limits of low-emissions combustors where the flame is more susceptible to perturbations. While the use of pilot flames is common in land-based gas turbine combustors, the mechanism by which they suppress instability is still unclear. In this study, we consider the impact of a central jet pilot on the stability of a swirl-stabilized flame in a variable-length, single-nozzle combustor. Previously, the pilot flame was found to suppress the instability for a range of equivalence ratios and combustor lengths. We hypothesize that combustion oscillation suppression by the pilot occurs because the pilot provides hot gases to the vortex breakdown region of the flow that recirculate and improve the static, and hence dynamic, stability of the main flame. This hypothesis is based on a series of experimental results that show that pilot efficacy is a strong function of pilot equivalence ratio but not pilot flow rate, which would indicate that the temperature of the pilot products as well as the combustion intensity of the pilot flame play more of a role in oscillation stabilization than the length of the pilot flame relative to the main flame. Further, the pilot-flame efficacy increases with pilot-flame equivalence ratio until it matches the main-flame equivalence ratio; at pilot equivalence ratios greater than the main equivalence ratio, the pilot-flame efficacy does not change significantly with pilot equivalence ratio. To understand these results, we use large-eddy simulation (LES) to provide a detailed analysis of the flow in the region of the pilot flame and the transport of radical species in the region between the main flame and pilot flame. Furthermore, the simulation, using a flamelet/progress variable-based chemistry tabulation approach and standard eddy viscosity/diffusivity turbulence closure models, provides detailed information that is inaccessible through experimental measurements.

42 ENGINEERING↗

Effect of delta tabs on mixing and axis switching in jets from asymmetric nozzles

The effect of delta tabs on mixing and the phenomenon of axis switching in free air jets from various asymmetric nozzles was studied experimentally. Flow visualization and Pitot probe surveys were carried out with a set of small nozzles (D = 1.47 cm) at a jet Mach number, Mj = 1.63. Hot wire measurements for streamwise vorticity were carried out with larger nozzles (D = 6.35 cm) at Mj = 0.31. Jet mixing with the asymmetric nozzles, as indicated by the mass fluxes downstream, was found to be higher than that produced by a circular nozzle. The circular nozzle with four delta tabs, however, produced fluxes much higher than that produced by a asymmetric nozzles themselves or by most of the tab configurations tried with them. Even higher fluxes could be obtained with only a few cases, e.g., with 3:1 rectangular nozzle with two large delta tabs placed on the narrow edges. In this case, the jet 'fanned out' at a large angle after going through one axis switch. The axis switching could be either stopped or augmented with suitable choice of the tab configurations. Two mechanisms are identified governing the phenomenon. One, as described in Ref. 12 and referred to here as the omega(sub Theta)-induced dynamics, is due to differential induced velocities of different segments of a rolled up azimuthal vortical structure. The other is the omega(sub x)-induced dynamics due to the induced velocities of streamwise vortex pairs in the flow. While the former dynamics are responsible for rapid axis switching in periodically forced jets, the effect of the tabs is governed mainly by the latter. It is inferred that both dynamics are active in a natural asymmetric jet issuing from a nozzle having an upstream contraction. The tendency for axis switching caused by the omega(sub Theta)-induced dynamics is resisted by the omega(sub x)-induced dynamics, leading to a delayed or no switch over in that case. In jets from orifices and in screeching jets, the omega(sub Theta)-induced dynamics dominate causing a faster switch over.

Zaman, K. B. M. Q.↗