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 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Thermovelocimetric characterization of liquid metal convection in a rotating slender cylinder

Rotating turbulent convection occurs ubiquitously in natural convective systems encompassing planetary cores, oceans and atmospheres, as well as in many industrial applications. While the global heat and mass transfer of water-like rotating Rayleigh–Bénard convection is well-documented, the dynamics in low-Prandtl-number liquid metals remain less understood. In this study, we experimentally investigate rotating Rayleigh-Bénard convection in liquid gallium (Prandtl number Pr ≈ 0.027) within a slender cylinder (diameter-to-height aspect ratio Γ = D/H =1/2) using novel thermovelocimetric diagnostic techniques that integrate simultaneous multi-point thermometry and ultrasonic Doppler velocity measurements. Our results reveal the formation of a stable, global-scale azimuthal wavenumber m = 2 quadrupolar vortex at low supercriticality. We propose that enhanced wall modes facilitated by the slender cylinder geometry interact with the bulk flow to create these large-scale axialized vortices. Furthermore, our findings imply a distinct scaling behavior for the wall-mode precession frequency in liquid metals, extending previous results obtained for moderate-Pr fluids. In conclusion, this provides new insights into wall-bulk coupling mechanisms of low-Pr rotating convective turbulence.

Magnetohydrodynamics↗

DOE OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗

DOE OSTI2025

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Dynamics of an inertially collapsing gas bubble between two parallel, rigid walls

The collapse of cavitation bubbles in channel flows can give rise to structural damage along neighbouring walls. Although the collapse of a bubble near a single wall has been studied extensively, less is known about bubble collapse between two walls, e.g. as in a channel. We conduct highly resolved, direct simulations of the Navier–Stokes equations to investigate the bubble dynamics and pressures produced by the collapse of a bubble between two parallel rigid walls. We examine the dependence of the dynamics and pressures on the initial bubble location, confinement and driving pressure. For a fixed initial stand-off distance, as the channel width increases the bubble volume, migration distance and re-entrant jet speed approach their single-wall counterparts. We obtain an expression for the minimum channel width at which the confinement does not affect the bubble dynamics depending on the driving pressure difference and initial stand-off distance. For a fixed channel width, varying stand-off distance reduced the maximum wall pressures in the channel relative to the single wall; the trend was consistent for three different driving pressures. Two different jetting behaviours are seen when the bubble is centred in the channel, depending on the channel width. Under significant confinement, wall-parallel re-entrant jets impinge upon each other and further intensify the collapse of the vortex ring.

Rodriguez, Jr, Mauro (ORCID:0000000305450265)↗

DOE OSTI2020

Dependence of Enstrophy Transport and Mixed Mass on Dimensionality and Initial Conditions in the Richtmyer–Meshkov Instability Induced Flows

Here we present a comparative study of the enstrophy budget and mixed mass between two- and three-dimensional flows induced by Richtmyer–Meshkov instability (RMI). Specifically, the individual contributions to the enstrophy budget due to the production from baroclinicity and from vortex stretching (which vanishes in two-dimensional (2D) flow) are delineated. This is enabled by a set of two- and three-dimensional computations at Atwood 0.5 having both narrow- and broad-band perturbations. A further three-dimensional (3D) computation is conducted at Atwood 0.9 using an identical narrowband perturbation to the Atwood 0.5 case to examine the sensitivity to density ratio. The mixed mass is also considered with the goal to obtain insight on how faithfully a simplified calculation performed in two dimensions can capture the mixed mass for an inertial confinement fusion (ICF) or other practical application. It is shown that the late time power law decay of variable density enstrophy is substantially different in two and three dimensions for the narrowband initial perturbation. The baroclinic production term is negligible in three dimensions (aside from the initial shock interaction), as vortex stretching is larger by two orders of magnitude. The lack of vortex stretching considerably reduces the decay rate in both narrowband and broadband perturbations in two dimensions. In terms of mixed mass, the lack of vortex stretching reduces the mixed mass in two dimensions compared to three in all cases. In the broadband cases, the spectral bandwidth in the 2D case is wider; hence, there is a longer time period of sustained linear growth which reduces the normalized mixed mass further.

42 ENGINEERING↗

DOE OSTI2022

A Large Eddy Simulation Study of Flow Turbulence, Alumina Transport, and Bath Temperature Evolution in Conventional Aluminum-Smelting Cell Using OpenFOAM

In this study, a Large Eddy Simulation (LES) of the aluminum-smelting process is performed using OpenFOAM. To understand the coupled behavior of heat transfer, mass transfer, and flow of the smelting process, a multi-physics computational fluid dynamics (CFD) model based on the Eulerian–Eulerian multi-fluid approach is adopted. The model accounts for CO 2 bubble and magnetohydrodynamics (MHD)-driven flow, along with alumina dissolution, transport, and bath temperature evolution. The simulation predictions show small-scale turbulent vortical structures in the anode–cathode space caused by combined effect of MHD and CO 2 bubble-bath interactions and relatively large-scale asymmetric vortices in the inter-anode space caused by the CO 2 bubble-bath interactions. The vortex formation at the edges of the anodes evidently aids in transporting alumina from the central channel to the bottom of the anodes and prevents accumulation of gas bubbles in the periphery of the anode bottom. Symmetric bath cold spots are observed in the vicinity of the feeder. Cold spots are also observed in the anode–cathode distance space below the anode bottom due to the transport of undissolved solid to this region by the flow. The findings from the work are useful in developing and designing alumina-feeding strategy leading to reduced anode effects and smooth operation of the cell. The work also highlights the important flow structures in conventional aluminum-smelting cell.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Analysis of the Aeroelastic Dynamics of Lightweight Flexible Variations of the SNL-NRT Turbine

Current trends show that wind turbines are growing in size to meet a rising demand for renewable energy generation, and their upscaled rotors have inherently become more flexible to maintain a proportionally lighter design. This is because larger rotors must be less massive relative to their diameter to minimize the levelized cost of energy (LCOE), which means that blades that are notably less stiff are produced as a result. These structural changes to blades are often reflected in their compromised aeroelastic stability and amplified deformation during operation, which has the potential to decrease the blade’s expected lifetime and the performance of the machine overall. Variations in blade flexibility are also known to influence vortex-wake structures downstream of the turbine, causing patterns of velocity deficit to evolve in ways that affect the performance of other turbines in the farm. This research explores how the increased flexibility of modern utility-scale wind turbine blades influences rotor aeroelastic behavior and interactions with farm flow. High-fidelity simulations of Sandia National Laboratories’ (SNL) National Rotor Testbed (NRT) wind turbine are presented. Flexible variations of the NRT baseline blade are simulated in a variety of realistic operational conditions typically expected at the SNL’s SWiFT facility in Lubbock, Texas. Solutions are then compared to investigate how specific changes to the structural properties of the NRT baseline blade’s design and construction can influence its aeroelastic response at the rotor and the evolution of the turbine’s wake.

Farrell, Alayna (ORCID:000000023555720X)↗

DOE OSTI2025

Metastability of stratified magnetohydrostatic equilibria and their relaxation

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Zinc oxide (ZnO) hybrid metasurfaces exhibiting broadly tunable topological properties

Extreme light confinement observed in periodic photonic structures, such as the vortex singularities in momentum (k) space, has been associated with their topological nature. Consequently, by exploiting and tuning their topological properties, optical metasurfaces have been demonstrated as an attractive platform for active photonics. However, given the fact that most active media under external excitations can only provide limited refractive index change, the potential advancements offered by the topological character of active metasurfaces have remained mostly unexplored. Zinc oxide (ZnO), which has recently exhibited optically-induced extraordinarily large permittivity modulations at visible and near-infrared frequencies, is an excellent active material for dynamic metasurfaces exhibiting strong tuning. This work demonstrates that a hybrid metasurface consisting of an array of ZnO nanodisks on a silver backplane displays broadly tunable topological properties. In particular, by performing k-space scattering simulations using measured pump-fluence-dependent material properties of ZnO, we study in detail the light reflection from the hybrid metasurface. Our results validate that the large k-space topology tuning of the metasurface can result in enormously strong polarization manipulation of near-infrared light in the vicinity of the topological features. The observed polarization switching effect is highly sensitive to the polarization and wavelength of an incident wave, owing to the symmetry and dispersion characteristics of the proposed system. Our study indicates that leveraging a combination of the extraordinary material properties and the k-space topology, hybrid metasurfaces based on ZnO may open new avenues for creating all-optical switchable metadevices.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Path Properties of Atmospheric Transitions: Illustration with a Low-Order Sudden Stratospheric Warming Model

Many rare weather events, including hurricanes, droughts, and floods, dramatically impact human life. To accurately forecast these events and characterize their climatology requires specialized mathematical techniques to fully leverage the limited data that are available. Here we describe transition path theory (TPT), a framework originally developed for molecular simulation, and argue that it is a useful paradigm for developing mechanistic understanding of rare climate events. TPT provides a method to calculate statistical properties of the paths into the event. As an initial demonstration of the utility of TPT, we analyze a low-order model of sudden stratospheric warming (SSW), a dramatic disturbance to the polar vortex that can induce extreme cold spells at the surface in the midlatitudes. SSW events pose a major challenge for seasonal weather prediction because of their rapid, complex onset and development. Climate models struggle to capture the long-term statistics of SSW, owing to their diversity and intermittent nature. We use a stochastically forced Holton–Mass-type model with two stable states, corresponding to radiative equilibrium and a vacillating SSW-like regime. In this stochastic bistable setting, from certain probabilistic forecasts TPT facilitates estimation of dominant transition pathways and return times of transitions. These “dynamical statistics” are obtained by solving partial differential equations in the model’s phase space. With future application to more complex models, TPT and its constituent quantities promise to improve the predictability of extreme weather events through both generation and principled evaluation of forecasts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Fundamental Superconductivity of Nb films for Quantum Computing Application

Niobium is a widely accepted material for quantum computing device as well as superconducting radio frequency (SRF) technology. Superconducting niobium is a marginal type II superconductor which has a very narrow gap (~20-30 mT) of the mixed state at 2K, even showing the intermediate state (IMS) at the early stage of magnetic vortex penetration. Tremendous progress has been made in understanding the impact of the Nb surface and bulk superconductivities on SRF resonator performance. However, the effect of thin superconducting Nb films for quantum computing applications requires further examination. In this study, we explore the fundamental superconducting properties of various thin Nb films and compare them with respect to the energy relaxation time, T1, measured from superconducting qubit fabricated with these films. The Nb films are fabricated both with or without a surface protective layer that prevents the formation of lossy Nb2O5. Electromagnetic properties are characterized by means of bulk magnetization, electromagnetic transport, and dynamics of surface superconductivity. In addition, analytical electron microscopy is implemented to further connect the superconducting properties to microstructure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Fundamental Superconductivity of Nb films for Quantum Computing Application

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Biasing topological charge injection in topological matter

Here we explore the interplay between topologies in the momentum and real spaces to formulate a thermodynamic description of nonequilibrium injection of topological charges under external bias. We show that the edge modes engendered by the momentum-space topology can play a functional role of connecting the external reservoirs to the bulk transport of topological charges in the real space. We illustrate our general results with two examples: the spin-torque injection of skyrmions in an electrically biased integer quantum Hall system, and the vortex injection in a topological $p + ip$ superconductor coupled to heat reservoirs. Based on the universal fractional entropy of the Majorana zero modes bound to the vortices, their controllable injection proposed in this work could provide a route for creating and manipulating Majorana fermions.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Quantum turbulence, superfluidity, non-Markovian dynamics, and wave function thermalization

While quantum turbulence has been addressed both experimentally (predominantly for superfluid He 4 and He 3 ) and theoretically, the dynamics of various ensembles of quantized vortices has been followed in time only until the vortices have decayed into phonons. How this “thermalization” is achieved is still an unaddressed and thus an unelucidated question. The unitary Fermi gas (UFG) is a unique quantum system, which has no classical counterpart and is of relevance to neutron stars, cold atoms, condensed-matter and nuclear many-body systems. The non-Markovian evolution of an isolated UFG is put in evidence and its entire nonequilibrium evolution can be studied theoretically within a unified theoretical framework. The initial lattice of quantum vortices and antivortices evolves through a couple of vortex tangles and excitation of Kelvin waves, where vortices cross and reconnect, until very slowly thermalization sets in. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Apparently superluminal superfluids

We consider the superfluid phase of a specific renormalizable relativistic quantum field theory. We prove that, within the regime of validity of perturbation theory and of the superfluid effective theory, there are consistent and regular vortex solutions where the superfluid’s velocity field as traditionally defined smoothly interpolates between zero and arbitrarily large superluminal values. We show that this solution is free of instabilities and of superluminal excitations. We show that, in contrast, a generic vortex solution for an ordinary fluid does develop an instability if the velocity field becomes superluminal. All this questions the characterization of a superfluid velocity field as the actual velocity of “something”.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Emergence of self-organized multivortex states in flocks of active rollers

Active matter, both synthetic and biological, demonstrates complex spatiotemporal self-organization and the emergence of collective behavior. A coherent rotational motion, the vortex phase, is of great interest because of its ability to orchestrate well-organized motion of self-propelled particles over large distances. However, its generation without geometrical confinement has been a challenge. Here, we show by experiments and computational modeling that concentrated magnetic rollers self-organize into multivortex states in an unconfined environment. We find that the neighboring vortices more likely occur with the opposite sense of rotation. Our studies provide insights into the mechanism for the emergence of coherent collective motion on the macroscale from the coupling between microscale rotation and translation of individual active elements. These results may stimulate design strategies for self-assembled dynamic materials and microrobotics.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Reversible to irreversible transitions for cyclically driven particles on periodic obstacle arrays

Here, we examine the collective dynamics of disks moving through a square array of obstacles under cyclic square wave driving. Below a critical density, we find that the system organizes into a reversible state in which the disks return to the same positions at the end of every drive cycle. Above this density, the dynamics are irreversible and the disks do not return to the same positions after each cycle. The critical density depends strongly on the angle θ between the driving direction and a symmetry axis of the obstacle array, with the highest critical densities appearing at commensurate angles such as θ = 0° and θ = 45° and the lowest critical densities falling at θ = arctan (0.618), the inverse of the golden ratio, where the flow is the most degenerate. As the density increases, the number of cycles required to reach a reversible state grows as a power law with an exponent near ν = 1.36, similar to what is found in periodically driven colloidal and superconducting vortex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Dynamic Mode Decomposition of Random Pressure Fields over Bluff Bodies

Fluctuating surface pressures on a bluff body exposed to a boundary layer flow generally are characterized as a spatiotemporally varying random field. In this paper, a dynamic mode decomposition (DMD) was applied to extract dominant features embedded in these random pressure fields. Utilizing an unsupervised machine learning algorithm, spatial modes and their temporal variations were grouped into different clusters at scales, e.g., macro, meso, and micro. A proper orthogonal decomposition (POD) of the experimental data was carried out to observe commonalities and distinctive perspectives each decomposition offers. Here, a comprehensive examination of the DMD/POD for their convergence criteria, data sufficiency, and modal components analysis was conducted. The physical interpretation of the spatiotemporal pressure field based on these decomposition schemes was discussed. At different scales, the DMD modes can capture the evolution of aerodynamic features, e.g., convection of vortices (or vortex tubes) and other structures. The distribution of energy among these three broad scales also reflects an energy cascade in pressure fluctuations akin to turbulence.

97 MATHEMATICS AND COMPUTING↗