Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 73 records · Page 4

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↗

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↗

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↗

Dean vortex-enhanced blood plasma separation in self-driven spiral microchannel flow with cross-flow microfilters

Point-of-care (POC) diagnostic devices have been developing rapidly in recent years, but they are mainly using saliva instead of blood as a test sample. A highly efficient self-separation during the self-driven flow without power systems is desired for expanding the point-of-care diagnostic devices. Microfiltration stands out as a promising technique for blood plasma separation but faces limitations due to blood cell clogging, resulting in reduced separation speed and efficiency. These limitations are mainly caused by the high viscosity and hematocrit in the blood flow. A small increment in the hematocrit of the blood significantly increases the pressure needed for the blood plasma separation in the micro-filters and decreases the separation speed and efficiency. Addressing this challenge, this study explores the feasibility of diluting whole blood within a microfluidic device without external power systems. This study implemented a spiral microchannel utilizing the inertial focusing and Dean vortex effects to focus the red blood cells and extract the blood with lower hematocrit. The inertial migration of the particles during the capillary flow was first investigated experimentally; a maximum of 88% of the particles migrated to the bottom and top equilibrium positions in the optimized 350 × 60 μm (cross-sectional area, 5.8 aspect ratio) microchannel. With the optimized dimension of the microchannel, the whole blood samples within the physiological hematocrit range were tested in the experiments, and more than 10% of the hematocrit reduction was compared between the outer branch outlet and inner branch outlet in the 350 × 60 μm microchannel.

Biochemistry & Molecular Biology↗

On closures for reduced order models - A spectrum of first-principle to machine-learned avenues

For over a century, reduced order models (ROMs) have been a fundamental discipline of theoretical fluid mechanics. Early examples include Galerkin models inspired by the Orr–Sommerfeld stability equation and numerous vortex models, of which the von Karman vortex street is one of the most prominent. Subsequent ROMs typically relied on first principles, like mathematical Galerkin models, weakly nonlinear stability theory, and two- and three-dimensional vortex models pioneered the data-driven proper orthogonal decomposition (POD) modeling. In early POD modeling, available data were used to build an optimal basis, which was then utilized in a classical Galerkin procedure to construct the ROM, but data have made a profound impact on ROMs beyond the Galerkin expansion. In this paper, we take a modest step and illustrate the impact of data-driven modeling on one significant ROM area. Specifically, we focus on ROM closures, which are correction terms that are added to the classical ROMs in order to model the effect of the discarded ROM modes in under-resolved simulations. Through simple examples, we illustrate the main modeling principles used to construct the classical ROMs, motivate and introduce modern ROM closures, and show how data-driven modeling, artificial intelligence, and machine learning have changed the standard ROM methodology over the last two decades. Finally, we outline our vision on how the state-of-the-art data-driven modeling can continue to reshape the field of reduced order modeling.

42 ENGINEERING↗

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗

Vortex rings in event-by-event relativistic heavy-ion collisions

We present event-by-event simulations for central asymmetric light+heavy and Au +Au collisions to investigate the formation and evolution of vortex-ring structures in the longitudinal flow velocity profile. The production-plane polarization of Λ hyperons, defined with respect to the Λ momentum and the beam, can track the “vortex-ring” feature in the event, a characteristic vortical structure generated by longitudinal flow gradients. We make comprehensive model predictions for the rapidity-dependent vortex-ring observables for different collision system sizes at √s NN = 200 and 72 GeV. Furthermore, our predictions at the latter energy can be explored in the future LHCb fixed-target experiment at the Large Hadron Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

LES simulations of a vacuum membrane distillation channel with geometric alterations

3D LES simulations were carried out to study the performance of a vacuum membrane distillation module. A wiggly wall profile with/without embedded stiffeners was considered to alleviate axial and radial temperature polarization, the cause of performance decrease in membrane distillation. Results of the flow field show the wiggles create unsteady vortex shedding inducing intense mixing in the channel. Vortex shedding intensity increases as the Reynolds number increases or stiffeners are added. 3D results show a 56% improvement in flux when moving from a flat sheet membrane to a wiggly membrane with stiffeners at a constant Reynolds number, corresponding to an improvement from 11.7 to 41.7 in the Nusselt number, showing that the alterations improved the flux performance by enhancing the heat transfer along the membrane surface and therefore alleviating temperature polarization, a critical bottleneck in membrane distillation systems. Further, a merit criterion was defined based on the Nusselt number and friction factor, and a 40% increase in merit was shown switching from a flat to a wiggly channel, while a 97% merit increase was seen going from flat to wiggly with stiffeners. A 42% enhancement in the flux was also seen moving from a straight channel to a wiggly channel at a higher Reynolds number which highlights the importance of strategically choosing the mass flow rate, as well as inducing flow separation and vortex shedding in the channel to promote mixing and dissipate the thermal boundary layer. The variation in results between the wiggly module with stiffeners and the wiggly module at high Reynolds numbers suggests enhancing mixing structures can be more impactful on flux for this geometry than increasing the flow rate to a turbulent/transitional regime. However, both are preferable for peak system flux performance. Furthermore, a 2D approximation was used to perform simulations on more extended channels to examine the length degradation. The modules with wiggly channels performed at the same flux level with a doubling of the length. On the other hand, the flat sheet modules experienced length degradation by temperature polarization and dropped in flux yield by around 13% of the short-channel value. This work illustrates that modeling the system and understanding how the performance decreases as the membrane surface area increases are critical for a larger module (scaling up from a lab to a prototype module) to maintain high flux performance.

42 ENGINEERING↗

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗

Validation of Calibrated k–ε Model Parameters for Jet-in-Crossflow

Previous efforts determined a set of calibrated, optimal model parameter values for Reynolds-averaged Navier–Stokes (RANS) simulations of a compressible jet in crossflow (JIC) using a $k–ε$ turbulence model. These parameters were derived by comparing simulation results to particle image velocimetry (PIV) data of a complementary JIC experiment under a limited set of flow conditions. Here, a $k–ε$ model using both nominal and calibrated parameters is validated against PIV data acquired from a much wider variety of JIC cases, including a realistic flight vehicle. The results from the simulations using the calibrated model parameters showed considerable improvements over those using the nominal values, even for cases that were not used in the calibration procedure that defined the optimal parameters. This improvement is demonstrated using a number of quality metrics that test the spatial alignment of the jet core, the magnitudes of multiple flow variables, and the location and strengths of vortices in the counter-rotating vortex cores on the PIV planes. These results suggest that the calibrated parameters have applicability well outside the specific flow case used in defining them and that with the right model parameters, RANS solutions for the JIC can be improved significantly over those obtained from the nominal model.

42 ENGINEERING↗

Turbulent flow characteristics in an 84-pin rod bundle for typical and damaged spacer grids

Hexagonal rod bundles arranged in a tightly packed triangular lattice are extensively used for heat transfer and energy generation applications. Staggered spacer grids are used to maintain the structural integrity of gas-cooled fast reactor (GFR) fuel assemblies, while inducing localized turbulence in flow. Damage to these spacer grids results in a disruption of flow fields within these hexagonal fuel bundles. Experimental flow visualizations are critical to identify the differences in local flow properties that the structural damage may cause. This experimental research investigates the flow-field characteristics at a near-wall and center plane in a prototypical 84-pin GFR fuel assembly. Newly installed typical spacers and spacers subject to naturally occurring damage due to material degradation over prolonged experimentation were investigated. Velocity fields were acquired by utilizing the matched-index-of-refraction method to obtain time-resolved particle image velocimetry measurements for a Reynolds number of 12 000. Reynolds decomposition statistical results divulged differences in the time-averaged velocity, velocity fluctuations, flow anisotropy, and Reynolds stress distributions. Galilean decomposition demarcated the influence of spacer grid damage on the velocity fields. To extract turbulent structures and elucidate mechanisms of flow instabilities, proper orthogonal decomposition analysis was employed. Reduced order flow reconstructions enabled the application of vortex identification algorithms to determine the spatial and statistical characteristics of vortices generated. This research work provides unique experimental data on the spacer grid condition-dependent flow. The results offer a deeper understanding of fluid dynamics behavior to support GFR rod bundle design efforts and computational fluid dynamics model validation.

36 MATERIALS SCIENCE↗

Laser powder bed fusion of Inconel 718 on 316 stainless steel

The joining of Inconel 718 and 316 stainless steel has notable industrial importance. However, traditional arc-welding has many issues. In particular, the formation of the detrimental Nb-rich Laves phase in the heat affected zone can promote crack initiation and propagation, which leads to poor mechanical properties. The small heat input and the rapid solidification of laser powder bed fusion (PBF) has the potential to suppress this problem. Despite the infancy of the multi-material PBF, it is of interest to investigate the relationship among the process parameter, the microstructure at the materials joining point, and the mechanical performance. In this study, we used high-speed x-ray imaging to in-situ monitor the printing of Inconel 718 tracks on a 316H SS substrate via laser powder bed fusion, followed by post-build characterization and testing with SEM, XRD and nanoindentation. We observed that a higher linear energy density resulted in a more diffused interface and a smaller compositional gradient, but as a trade-off, generated more keyhole pores and thermal cracks. Furthermore, we also showed that the built layer exhibited chemical inhomogeneity at two length scales: a tens-of-mu m-scale compositional vortex due to the melt flow, and a sub-mu m-scale Nb-rich dendritic structure due to the non-equilibrium solidification that led to the formation of the Laves phase. It was inferred that the chemical inhomogeneity may benefit the mechanical properties at the interface by providing more interlocking between the two materials and also by limiting the effect of the brittle Laves phase to potentially short range.

36 MATERIALS SCIENCE↗

Design and Testing of the Vortex Ring Facility

This document summarizes the efforts of building an experimental test facility to study the evolution of the cloud following a nuclear detonation above ground level and support the development of numerical models to describe it. The experimental facility allows nonintrusive flow visualization and measurements of buoyant vortex ring formation and evolution using background-oriented Schlieren and particle image velocimetry techniques. In addition, 3D unsteady computational fluid dynamics simulations using unsteady Reynolds-averaged Navier–Stokes and volume-of-fluid numerical methods were performed to support the experimental design as well as to provide insight about the formation and evolution of the vortex ring.

42 ENGINEERING↗

Applying Time-Parallelization to Turbulent Flows

Parallelization of the temporal domain is explored for the solution of turbulent flows. Multigrid reduction-in-time (MGRIT) is used to advance the large-scale fluid dynamics in time sequentially on the coarsest space-time grid but propagate the information in time parallel on all other levels. The goal of this process is to accurately and efficiently resolve the coarse-scale turbulence structure and use that to drive the fine-scales of the turbulent flow. The extra forcing from nonlinear multigrid facilitates the coupling and interaction between fine and coarse scales, through which the multiscale nonlinear physics is properly captured. Adaptive mesh refinement is employed to finely resolve only the regions with strong gradients, which provides further computational efficiency. The underlying computational fluid dynamics solver is a fourth-order finite-volume scheme with the standard 4-stage Runge-Kutta method. An advanced approach is devised and implemented to enable MGRIT to solve highly turbulent flows successfully. Furthermore, the method is applied to solve a Taylor-Green vortex problem and a doubleshear-layer turbulent mixing flow. Results are promising, validating that MGRIT with the filtering approach has the potential to efficiently solve general turbulent flows.

Computational Fluid Dynamics↗

Vortex Ordering and Dynamics of Santa Fe Artificial Ice Pinning Arrays

We numerically examine the ordering, pinning, and flow of superconducting vortices interacting with a Santa Fe artificial ice pinning array. We find that as a function of magnetic field and pinning density, a wide variety of vortex states occur, including ice rule obeying states and labyrinthine patterns. In contrast to square pinning arrays, we find no sharp peaks in the critical current due to the inherent frustration effect imposed by the Santa Fe ice geometry; however, there are some smoothed peaks when the number of vortices matches the number of pinning sites. For some fillings, the Santa Fe array exhibits stronger pinning than the square array due to the suppression of one-dimensional flow channels when the vortex motion in the Santa Fe lattice occurs through the formation of both longitudinal and transverse flow channels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effects of Downstream Vortex Generators on Film Cooling a Flat Plate Fed by Crossflow

Counter-rotating vortices, formed by the interaction of film-cooling jets and the hot gas flow, adversely affect the performance of conventional film-cooling designs. Downstream vortex generators have been shown to improve cooling effectiveness by mitigating the effects of the counter-rotating vortices and by deflecting the cooling jet laterally. In this study, computational and experimental methods were used to examine how cylindrical film-cooling holes (D = 3.2 mm, L/D = 6, p/D = 3, α = 30 deg) with and without downstream vortex generators perform when the coolant supply channel is perpendicular to the direction of the hot gas. For this study, the hot gas had a temperature of 650 K and an average Mach number of 0.23. Here, the hot-gas-to-coolant temperature ratio was 1.9, and two blowing ratios (0.75 and 1.0) were studied. Results from the computational fluid dynamics study show how crossflow affects the interaction between the film-cooling jet and hot gas flow with and without downstream vortex generators. The experimental measurements were based on infrared thermography in a conjugate heat transfer environment. Results were obtained for film-cooling performance in terms of overall effectiveness, film effectiveness, and local heat transfer coefficients. The downstream vortex generators can increase the laterally averaged effectiveness by a factor of 1.5 relative to cylindrical holes, but this higher performance is restricted to low crossflow velocities and higher blowing ratios.

17 WIND ENERGY↗

Vortex Induced Vibration Test of the WITT Device in a Tidal Channel

Fluid flow around a bluff body causes vibrations as the boundary layer separates and vortices are shed. This phenomenon, called vortex induced vibration (VIV), is minimized in most engineering applications, but can be utilized for power generation. Vortex induced vibration has potential for energy harvesting in low velocity water flows and therefore could open many locations to marine renewable energy. Here we present the results of a field test of a device that utilizes VIV to harvest energy from flow in a tidal channel. The pendulum-type device has been designed to generate electricity from motion around and along three axes. It is named the Whatever Input to Torsion Transfer (WITT) energy device and in this test is mounted on top of a pipe. During the field test, the WITT generated electricity at flow velocities from 0.27 to 2.82 m/s. The field test has proven that the WITT can utilize VIV to produce electricity and charge a battery in tidal flows with velocities as low as 0.27 m/s. In the future, the electricity could be used to charge a battery for ocean observation instruments such as temperature and nitrate sensors. In-situ power generation by the WITT would extend deployment times and increase the frequency of measurements.

Branch, Ruth A.↗

Hydrokinetic Energy Conversion by Flow-Induced Oscillation of Two Tandem Cylinders of Different Stiffness

Abstract The vortex-induced vibration for aquatic clean energy (VIVACE) converter harnesses hydrokinetic energy by enhancing flow-induced oscillations (FIOs) of elastically supported rigid cylinders in a river, tide, or ocean current. The harnessing power depends on the intensity of the oscillation, which is a consequence of the flow–structure interaction. The inflow condition for the downstream (second) cylinder is slowed down and perturbed by the upstream (first) cylinder, due to the shielding effect. Therefore, the optimal structural parameters, i.e., stiffness and damping ratio, for the second cylinder may be different from the first cylinder, in terms of energy harnessing. To improve the performance of the VIVACE converter, a series of experiments are conducted in a recirculating water channel, with various stiffness combinations of two cylinders in tandem. Results show that the stiffness of the second cylinder, K2, does not affect the energy harnessing power in vortex-induced vibration (VIV) occurring at low speeds, because the oscillation of the downstream cylinder in this velocity range is completely dominated by the wake of the upstream cylinder. K2 has a great influence on the harnessing power at higher velocities in the transition region from VIV to galloping and in galloping. Changing K2 onsets and enhances galloping at lower flow velocity and harnesses up to 110% more energy than the case of K1 = K2.

Engineering↗