SEARCH · Engineering Papers

Results for “Vortex dynamics”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2021

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↗

DOE OSTI2024

Analyzing the Effects of Atmospheric Turbulent Fluctuations on the Wake Structure of Wind Turbines and Their Blade Vibrational Dynamics

In recent trends, a rising demand for renewable energy has driven wind turbines to larger proportions, where lighter blade designs are often adopted to reduce the costs associated with logistics and production. This causes modern utility-scale wind turbine blades to be inherently more flexible, and their amplified aeroelastic sensitivity results in complex multi-physics reactions to variant atmospheric conditions, including dynamic patterns of aerodynamic loading at the rotor and vortex structure evolutions within the wake. In this paper, we analyze the influence of inflow variance for wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex wake interactions for a range of flexible NRT blade variations, operating in differing conditions of variant atmospheric flow. CODEF solutions of turbine operation in Steady-In-The-Average (SITA) wind conditions are compared to SITA wind conditions featuring a controlled gust-like pulse overimposed, to isolate the effects of typical wind fluctuations. Finally, simulations of realistic time-varying wind conditions from SWiFT meteorological tower measurements are compared to the solutions of SITA wind conditions. These increasingly complex atmospheric inflow variations are tested to show the differing effects evoked by various patterns of spatiotemporal atmospheric flow fluctuations. An analysis is presented for solutions of wind turbine aeroelastic response and vortex wake evolution, to elucidate the consequences of variant inflow, which pertain to wind turbine dynamics at an individual and farm-collective scale. The comparisons of simulated farm flow for SITA and measured fluctuating wind conditions show that certain regions of the wake contain up to a 12% difference in normalized axial velocity, due to the introduction of wind fluctuations. The findings of this study prove valuable for practical applications in wind farm control and optimization strategies, with particular significance for modern utility-scale wind power plants operating in variant atmospheric conditions.

Farrell, Alayna (ORCID:000000023555720X)↗

DOE OSTI2024

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Gravitational vortex disruption in mesofluidic separators

Typically, bump array devices separate submicrometer particles where gravitational effects remain small. However, as particle size increases from a few micrometers to hundreds of micrometers, gravitational effects become important and settling may induce vortex formation. Here, this article shows under laminar flow that these vortices suffice to disrupt flow fields, frustrating particle separation. The limiting diameter, above which vortices do form, is determined, and the dependence of the limiting diameter on viscosity, particle properties, and separator orientation is quantified.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

On the interaction of a wind turbine wake with a conventionally neutral atmospheric boundary layer

In this work, we investigate the dynamics of wind turbine tip-vortex breakdown in a conventionally neutral atmospheric boundary layer (ABL). To this end, high-resolution data are collected from large-eddy simulations of a wind turbine operating within a neutral ABL and studied by means of proper orthogonal decomposition (POD) and Fourier analysis. The high resolution of the generated data in both space and time allows us to gain insight into the tip-vortex breakdown mechanisms by (i) capturing the energy modes of the coherent structures, (ii) studying their contribution to the tip-vortex breakdown through their power spectra functions and mean kinetic energy (MKE) flux, and (iii) analysing the growth rate of each contributing perturbation frequency along tip vortices. Our analysis shows that under a fully turbulent scenario, the growth rate of perturbations along the tip vortices is largest for low wave numbers, i.e. long-wave perturbations. Additionally, the MKE flux reaches its highest value at two diameters downstream of the rotor plane, a behaviour that can be attributed to the coexistence of multiple interacting POD modes, with the streamwise vortex roller mode being the primary contributor to the total MKE flux budget, contributing approximately . Finally, comparisons with a laminar, uniform flow scenario subject to a single-frequency perturbation highlight the differences between the two ambient flow conditions. In the non-turbulent, uniform flow scenario, the growth rate attains its maximum value at a wave number corresponding to the out-of-phase mutual-inductance mechanism, whereas the MKE flux exhibits local minima and maxima along the wake and at different downstream locations depending on the perturbation frequency. Our analyses suggest that the breakdown of the wind turbine tip vortices under a fully turbulent neutral ABL inflow is due to complex interactions across a range of excitation frequencies, in which the mutual-inductance instability may not be the dominant one.

17 WIND ENERGY↗

DOE OSTI2022

Fluid structure interaction of a subaqueous pendulum: Analyzing the effect of wake correction via large eddy simulations

The dynamic behavior of a subaqueous cylindrical pendulum and corresponding flow dynamics are investigated. The objectives were twofold: (i) to examine whether the two-dimensional model equations sufficiently capture the three-dimensional dynamics and (ii) to investigate the emerging three-dimensional vortical flow structures. Large eddy simulations with two-way coupling fluid structure interaction were carried out using the immersed boundary method to simulate the motion of the pendulum and its interactions with the initially stagnant water. The resulting pendulum motion is compared against measured data obtained in a series of experimental tests to validate the simulation results and the model equations with and without wake corrections. An analysis of the flow vorticity revealed the development of a vortex ring during the first swing and the formation of tip vortices. The evolution of the vortex rings emerging from the motion of the subaqueous cylindrical pendulum was visualized using Q-criteria showing a reasonable agreement with vortical structures observed in the experiment using particle imaging velocimetry. The hydrodynamic moments acting on the simulated pendulum and the moments calculated from the model equations are analyzed. Using the insights from these numerical simulations, a modification of the wake correction is proposed to enhance the accuracy of the rate of decay and period. The transient effect of coherent flow on pendulum dynamics, especially the added mass effect, is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Numerical Simulations of Cryogenic Hydrogen Cooling in Vortex Tubes with Smooth Transitions

Improving efficiency of hydrogen cooling in cryogenic conditions is important for the wider applications of hydrogen energy systems. The approach investigated in this study is based on a Ranque-Hilsch vortex tube (RHVT) that generates temperature separation in a working fluid. The simplicity of RHVT is also a valuable characteristic for cryogenic systems. In the present work, novel shapes of RHVT are computationally investigated with the goal to raise efficiency of the cooling process. Specifically, a smooth transition is arranged between a vortex chamber, where compressed gas is injected, and the main tube with two exit ports at the tube ends. Flow simulations have been carried out using STAR-CCM+ software with the real-gas Redlich-Kwong model for hydrogen at temperatures near 70 K. It is determined that a vortex tube with a smooth transition of moderate size manifests about 7% improvement of the cooling efficiency when compared vortex tubes that use traditional vortex chambers with stepped transitions and a no-chamber setup with direct gas injection.

08 HYDROGEN↗

DOE OSTI2022

Australian wildfire smoke in the stratosphere: the decay phase in 2020/2021 and impact on ozone depletion

Abstract. Record-breaking wildfires raged in southeastern Australia in late December 2019 and early January 2020. Rather strong pyrocumulonimbus (pyroCb) convection developed over the fire areas and lofted enormous amounts of biomass burning smoke into the tropopause region and caused the strongest wildfire-related stratospheric aerosol perturbation ever observed around the globe. We discuss the geometrical, optical, and microphysical properties of the stratospheric smoke layers and the decay of this major stratospheric perturbation. A multiwavelength polarization Raman lidar at Punta Arenas (53.2∘ S, 70.9∘ W), southern Chile, and an elastic backscatter Raman lidar at Río Grande (53.8∘ S, 67.7∘ W) in southern Argentina, were operated to monitor the major record-breaking event until the end of 2021. These lidar measurements can be regarded as representative for mid to high latitudes in the Southern Hemisphere. A unique dynamical feature, an anticyclonic, smoke-filled vortex with 1000 km horizontal width and 5 km vertical extent, which ascended by about 500 m d−1, was observed over the full last week of January 2020. The key results of the long-term study are as follows. The smoke layers extended, on average, from 9 to 24 km in height. The smoke partly ascended to more than 30 km height as a result of self-lofting processes. Clear signs of a smoke impact on the record-breaking ozone hole over Antarctica in September–November 2020 were found. A slow decay of the stratospheric perturbation detected by means of the 532 nm aerosol optical thickness (AOT) yielded an e-folding decay time of 19–20 months. The maximum smoke AOT was around 1.0 over Punta Arenas in January 2020 and thus 2 to 3 orders of magnitude above the stratospheric aerosol background of 0.005. After 2 months with strongly varying smoke conditions, the 532 nm AOT decreased to 0.03-0.06 from March–December 2020 and to 0.015–0.03 throughout 2021. The particle extinction coefficients at 532 nm were in the range of 10–75 Mm−1 in January 2020 and, later on, mostly between 1 and 5 Mm−1. Combined lidar–photometer retrievals revealed typical smoke extinction-to-backscatter ratios of 69 ± 19 sr (at 355 nm), 91 ± 17 sr (at 532 nm), and 120 ± 22 sr (at 1064 nm). An ozone reduction of 20 %–25 % in the 15–22 km height range was observed over Antarctica and New Zealand ozonesonde stations in the smoke-polluted air, with particle surface area concentrations of 1–5 µm2 cm−3.

Ohneiser, Kevin↗

DOE OSTI2026

Josephson dynamics in two-dimensional ring-shaped condensates

We investigate Josephson transport in a fully closed, two-dimensional superfluid circuit formed by a ring-shaped 87 Rb Bose-Einstein condensate that contains two optical barriers acting as movable weak links. Translating these barriers at controlled speeds imposes a steady bias current, enabling direct mapping of the current-chemical-potential (𝐼−Δ⁢𝜇) characteristics. For narrow junctions (𝑤 ≈ 1µ⁢m), the circuit exhibits a pronounced dc branch that terminates at a critical current 𝐼 𝑐 = 9⁢(1) × 10 3 s −1 ; above this threshold, the system switches to an ac, resistive regime. Classical-field simulations that include the moving barriers quantitatively reproduce both the nonlinear 𝐼−Δ⁢𝜇 curve and the measured 𝐼 𝑐 , validating the underlying microscopic picture. Analysis of the ensuing phase dynamics shows that dissipation is mediated by the nucleation and traversal of vortex-antivortex pairs through the junctions, while the bulk condensate remains globally phase locked—direct evidence of the ring's topological constraint enforcing quantized circulation. These results establish a cold-atom analog of a superconducting quantum interference device in which Josephson dynamics can be resolved at the single-vortex level, providing a versatile platform for atomtronic circuit elements, nonreciprocal Josephson devices, and on-chip Sagnac interferometers for multiaxis rotation sensing.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

The kinetic origin of the fluid helicity—A symmetry in the kinetic phase space

Helicity, a topological degree that measures the winding and linking of vortex lines, is preserved by ideal (barotropic) fluid dynamics. In the context of the Hamiltonian description, the helicity is a Casimir invariant characterizing a foliation of the associated Poisson manifold. Casimir invariants are special invariants that depend on the Poisson bracket, not on the particular choice of the Hamiltonian. The total mass (or particle number) is another Casimir invariant, whose invariance guarantees the mass (particle) conservation (independent of any specific choice of the Hamiltonian). In a kinetic description (e.g., that of the Vlasov equation), the helicity is no longer an invariant (although the total mass remains a Casimir of the Vlasov’s Poisson algebra). The implication is that some “kinetic effect” can violate the constancy of the helicity. To elucidate how the helicity constraint emerges or submerges, we examine the fluid reduction of the Vlasov system; the fluid (macroscopic) system is a “sub-algebra” of the kinetic (microscopic) Vlasov system. In the Vlasov system, the helicity can be conserved if a special helicity symmetry condition holds. To put it another way, breaking helicity symmetry induces a change in the helicity. We delineate the geometrical meaning of helicity symmetry and show that for a special class of flows (the so-called epi-two-dimensional flows), the helicity symmetry is written as ∂γ = 0 for a coordinate γ of the configuration space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Improvements to the Blade Element Momentum Formulation of OpenFAST for Skewed Inflows

In this work, we modify the blade element momentum algorithm of OpenFAST to improve its predictions under large skewed inflow conditions. We use the well-known Glauert's skew correction and introduce continuous extension of the model for high-thrust conditions. We present the rationale behind Glauert's empirical model and discuss the different conventions possible for the axial induction factor. We verify the model against the higher-fidelity lifting-line vortex method and blade-resolved computational fluid dynamics, and we observe that the new implementation enhances the accuracy and reliability of OpenFAST's aerodynamic modeling capabilities in conditions involving large skew angles. For the parametric studies run using the different codes, we find that the power changes with the skew angle as cos 1.7 (θ skew ) and the thrust as cos 0.65 (θ skew ). An analysis of the azimuthal variation of the induced velocities in the rotor plane reveals that current redistribution models used in blade element momentum codes may need to be refined.

17 WIND ENERGY↗

DOE OSTI2020

Phase-matching analysis in high-order harmonic generation with nonzero orbital angular momentum Laguerre-Gaussian beams

Through high-order harmonic generation driven by intense ultrashort vortex infrared or midinfrared lasers, a nonzero orbital angular momentum can be imprinted onto extreme ultraviolet (XUV) or soft-x-ray (SXR) light pulses. Here we simulate the generation of vortex XUV harmonics in the gas medium as well as their propagation in vacuum till reaching the far field. In this work, we find that the intensity and phase of generated high harmonics are very sensitive to the position of gas jet with respect to the laser focus. The topological charge of the $\textit{q}$th harmonic is found to be $\textit{q}$ times that of the driving Laguerre-Gaussian beam. Each harmonic in the far field appears as a single ring in the transverse plane with an invariant diameter which is scalable with the fundamental topological charge only when the gas jet is placed after the laser focus. The underlying phase-matching mechanism is analyzed by examining the spatial map of the coherence length and by calculating the evolution of harmonic emission in the medium. We anticipate this work to stimulate interest in generating intense vortex XUV or SXR attosecond pulses for probing dynamics of molecules where special molecular features are difficult to be detected with linear or circular XUV or SXR pulses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Tidal Turbine Benchmarking Project: Stage I - Steady Flow Blind Predictions: Preprint

This paper presents the first blind prediction stage of the Tidal Turbine Benchmarking Project being conducted and funded by the UK's EPSRC and Supergen ORE Hub. In this first stage, only steady flow conditions, at low and elevated turbulence (3.1%) levels, were considered. Prior to the blind prediction stage, a large laboratory scale experiment was conducted in which a highly instrumented 1.6m diameter tidal rotor was towed through a large towing tank in well-defined flow conditions with and without an upstream turbulence grid. Details of the test campaign and rotor design were released as part of this community blind prediction exercise. Participants were invited to use a range of engineering modelling approaches to simulate the performance and loads of the turbine. 26 submissions were received from 12 groups from across academia and industry using solution techniques ranging from blade resolved computational fluid dynamics through actuator line, boundary integral element methods, vortex methods to engineering Blade Element Momentum methods. The comparisons between experiments and blind predictions were extremely positive helping to provide validation and uncertainty estimates for the models, but also validating the experimental tests themselves. The exercise demonstrated that the experimental turbine data provides a robust data set against which researchers and design engineers can test their models and implementations to ensure robustness in their processes, helping to reduce uncertainty and provide increased confidence in engineering processes. Furthermore, the data set provides the basis by which modellers can evaluate and refine approaches.

benchmarking↗

DOE OSTI2023

Rate of Para-Orthohydrogen Conversion in Cryogenic Vortex Tube

Raising efficiency of the cooling process for cryogenic hydrogen and minimizing hydrogen boil-off during storage and transportation of liquid hydrogen are critically important factors for widening implementation of hydrogen-based clean energy systems. One novel approach involves utilization of relatively simple catalyzed vortex tubes, where hydrogen in the primarily para- nuclear spin isomer form is converted into the ortho- form while consuming a significant amount of heat. To design such systems effectively, better understanding is required about catalyst-assisted para-orthohydrogen conversion rates in high-speed vortical flows. In this study, a computational fluid dynamics simulation has been set up to model an experimental system with cryogenic hydrogen vortex tubes. Mesh-verification and a validation study has been conducted first for a non-catalyzed tube. Then, the rate coefficient of the para-orthohydrogen conversion of cryogenic hydrogen has been determined by matching numerical results with experimental data available for a catalyzed vortex tube. Furthermore, the presented information can help design and optimize cryogenic hydrogen cooling within vortex tubes.

08 HYDROGEN↗

DOE OSTI2021

Relaxation of a single defect spin by the low-frequency gyrotropic mode of a magnetic vortex

We excite the gyrotropic mode of a magnetic vortex and observe the resulting effect on the spin state of a nearby nitrogen-vacancy (NV) defect in diamond. Thin permalloy disks fabricated on a diamond sample are magnetized in a vortex state in which the magnetization curls around a central core. The magnetization dynamics of this configuration are described by a discrete spectrum of confined magnon modes as well as a low-frequency gyrotropic mode in which the vortex core precesses about its equilibrium position. Despite the spin transition frequencies being far-detuned from the modes of the ferromagnet, we observe enhanced relaxation of the NV spin when driving the gyrotropic mode. Moreover, we map the spatial dependence of the interaction between the vortex and the spin by translating the vortex core within the disk with an applied magnetic field, resulting in steplike motion as the vortex is pinned and de-pinned. Strong spin relaxation is observed when the vortex core is within approximately 250 nm of the NV center defect. We attribute this effect to the higher frequencies in the spectrum of the magnetic fringe field arising from the soliton-like nature of the gyrotropic mode when driven with sufficiently large amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Application of passive vortex generators to enhance vertical mixing in an open raceway pond

A novel use of vortex generators in a raceway pond is described here, with the flow field quantitatively simulated using computational fluid dynamics using the large eddy simulation turbulence model. Persistence lengths of the swirling motion generated by the vortex generators indicate that significant vertical mixing can be achieved by placing vortex generators in the straight section opposite the paddle wheel, downstream of the first hairpin bend. Relatively simple vortex generators are capable of creating stronger swirling motions that persist for a longer distance than those created by the paddle wheel. For optimal performance, vortex generators are positioned side by side but in opposite directions, and their diameters should be equal to or slightly less than the desired liquid depth. The optimal length of a 0.18 m diameter vortex generator in a 0.2 m deep pond was determined to be 0.3 m. As a result, it has been demonstrated that a longer persistence length is achieved by inducing a swirling motion with its rotational axis parallel to the primary flow direction.

09 BIOMASS FUELS↗

DOE OSTI2022

Dynamics and bifurcations of laminar annular swirling and non-swirling jets

This paper presents bifurcation analyses characterising the nonlinear dynamics of fully developed laminar annular jets with respect to the centrebody diameter ${ {d}}$ , Reynolds number ${ {Re}}$ , and swirl ratio ${ {S}}$ . Similar flows appear in numerous applications and feature a vibrant range of topological and dynamical characteristics associated with phenomena including shear layer separation and vortex breakdown. Our results begin by describing the non-monotonic evolution of the axisymmetric jet's steady topology under varying ${ {S}}$ . In accord with earlier reports, the jet progresses through a sequence of wake, breakdown and wall jet regimes in a qualitatively similar manner across a wide span of ${ {d}}$ and ${ {Re}}$ values. In the wake regime, the non-swirling jet bifurcates to a plane-symmetric, but not axisymmetric, steady flow pattern beyond a ${ {d}}$ -dependent critical ${ {Re}}$ value. With further increase in ${ {Re}}$ , the steady non-swirling jet destabilises subsequently via multiple distinct Hopf bifurcations. Introducing ${ {S}}>0$ to the jet also induces unsteadiness by twisting the singly azimuthally periodic ( $|m|=1$ ) asymmetric wake structure and causing it to precess periodically in time about the central axis. Intermediate swirl stabilises this unsteady dynamics and restores the jet's axisymmetry. This stabilising effect is then reversed in the breakdown regime at higher ${ {S}}$ , where a variety of different $|m|=1$ and $|m|=2$ instabilities bifurcate from the steady flow as ${ {S}}$ is increased. Several instances of hysteresis and subcritical behaviour are reported and discussed, including one that manifests precessing vortex core oscillations.

Douglas, Christopher M. (ORCID:0000000259683315)↗

DOE OSTI2025

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING↗