Engineering Papers⌕ Search

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.

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

At least 397 records · Page 22

Imaging of super-fast dynamics and flow instabilities of superconducting vortices

Quantized magnetic vortices driven by electric current determine key electromagnetic properties of superconductors. And while the dynamic behavior of slow vortices has been thoroughly investigated, the physics of ultrafast vortices under strong currents remains largely unexplored. Here, we use a nanoscale scanning superconducting quantum interference device to image vortices penetrating into a superconducting Pb film at rates of tens of GHz and moving with velocities of up to tens of km/s, which are not only much larger than the speed of sound but also exceed the pair-breaking speed limit of superconducting condensate. These experiments reveal formation of mesoscopic vortex channels which undergo cascades of bifurcations as the current and magnetic field increase. Our numerical simulations predict metamorphosis of fast Abrikosov vortices into mixed Abrikosov-Josephson vortices at even higher velocities. Our work offers an insight into the fundamental physics of dynamic vortex states of superconductors at high current densities, crucial for many applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimental and Numerical Investigation of a Miniature Additively Manufactured Vortex Tube

Ranque–Hilsch vortex tubes are simple devices that can produce a cooling effect using compressed air. A key advantage of vortex tubes is the lack of moving solid parts; however, their efficiencies are relatively low. The present study focuses on the development of a miniature variable-diameter tube using additive manufacturing. A metal-based 3D printing technique was utilized to fabricate this vortex tube monolithically. Computational fluid dynamics simulations employing software STAR-CCM + with a compressible Reynolds-Averaged Navier–Stokes (RANS) approach and the elliptic-blending lag k-epsilon turbulence model have been applied to model thermofluid processes inside the vortex tube, to good agreement with the experiment. In conclusion, a temperature decrease of 13.3 °C and a cooling power of approximately 4 W were experimentally achieved with a pressure ratio of 4 in the air at normal conditions. This result shows promise for the goal of utilizing additive manufacturing to design and build complex-geometry vortex tubes intended for use with cryogenic fluids.

08 HYDROGEN↗

The internal structure of a vortex in a two-dimensional superfluid with long healing length and its implications

We analyze the motion of quantum vortices in a two-dimensional spinless superfluid within Popov’s hydrodynamic description. In the long healing length limit (where a large number of particles are inside the vortex core) the superfluid dynamics is determined by saddle points of Popov’s action, which, in particular, allows for weak solutions of the Gross–Pitaevskii equation. We solve the resulting equations of motion for a vortex moving with respect to the superfluid and find the reconstruction of the vortex core to be a non-analytic function of the force applied on the vortex. This response produces an anomalously large dipole moment of the vortex and, as a result, the spectrum associated with the vortex motion exhibits narrow resonances lying within the phonon part of the spectrum, contrary to traditional view.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

In this work, we will discuss observations from images acquired from a time-resolved, high-speed self-aligned focusing Schlieren campaign that was performed at the CW-22 facility at NASA Glenn Research Center to understand the dynamic behavior of thick trailing-edge high-pressure turbine blades simulating a ceramic matrix composite (CMC) construction at high inlet turbulence conditions. For the CMC-9 blade, which has a trailing edge thickness of 9% of the axial chord, we identified a regime where an excessive total pressure loss (loss anomaly) is observed only for a narrow range of Reynolds numbers at a fixed exit Mach number of 0.74. The loss anomaly is qualitatively explained by our images, which were taken at 10 distinct blade Reynolds numbers spanning a factor of 6. The images show a significant increase in energy related to the oscillations due to transonic vortex shedding at the Reynolds numbers related to the high loss conditions. This increased energy leads to the formation of strong acoustic waves that turn into shock waves at the highest loss conditions. From our observations stemming from Spectral POD analysis of the high-speed images, we see the acoustic waves produced by the trailing edge vortex shedding exist in all conditions tested; but the shedding frequency has a very slight trend upwards as the Reynolds number is increased, varying about 6% in the range tested. Considering this variation of shedding frequency as a function of Reynolds number, which is well-established for other bluff-body flows, we stipulate there may be a potential feedback mechanism involving an acoustic information transfer path across neighboring blades that may explain why only a narrow range of Reynolds numbers displays strong, shock-forming vortex shedding. We consider a few feedback paths and examine the timing based on the mean flow field from a high-resolution LES simulation. It appears that all feedback mechanisms are viable, presenting an integer number of delay cycles with respect to disturbances generated at the trailing edge. Most noteworthy, however, is the acoustic/shock-boundary layer interaction mechanism at the blade geometric throat. Based on our analysis, this mechanism likely plays an important role in realistic, curved turbine blade passages.

Aerodynamics↗

Aerodynamics of vortex generators

An experimental and theoretical study was undertaken of the separation delay and dramatic boundary-layer thinning that can occur in vortex-generator installations. Wind tunnel measurements of the dynamic-pressure profile downstream of a vortex generator were found to compare under certain conditions with that downstream of a suction slit, while water-tunnel visualization studies of vortex-generator height and geometry suggested optimum configurations, and only a minor effect of base porosity. A series of progressively more complex inviscid flow models was developed to be applied to a 3-D integral boundary-layer code. This code predicted layer thinning downstream of the suction site of the vortex models, and other observed features. Thin-layer Navier-Stokes equations are now being used with the ultimate goal of clarifying the physical processes involved in vortex generator performance and developing calculational procedures capable of predicting it.

Breidenthal, Robert E., Jr.↗

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II

Here, in this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part I

In this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In conclusion, these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗

Lagrangian Transport Calculations Using UARS Data: Ozone - Part 2

Trajectory calculations are used to examine ozone transport in the polar winter stratosphere during periods of the Upper Atmosphere Research Satellite (UARS) observations. The value of these calculations for determining mass transport was demonstrated previously using UARS observations of long-lived tracers, In the middle stratosphere, the overall ozone behavior observed by the Microwave Limb Sounder in the polar vortex is reproduced by this purely dynamical model. Calculations show the evolution of ozone in the lower stratosphere during early winter to be dominated by dynamics in December 1992 in the Arctic. Calculations for June 1992 in the Antarctic show evidence of chemical ozone destruction and indicate that approx. 50% of the chemical destruction may be masked by dynamical effects, mainly diabatic descent, which bring higher ozone into the lower-stratospheric vortex. Estimating differences between calculated and observed fields suggests that dynamical changes masked approx. 20% - 35% of chemical ozone loss during late February and early March 1993 in the Arctic. In the Antarctic late winter, in late August and early September 1992, below approx. 520 K, the evolution of vortex-averaged ozone is entirely dominated by chemical effects; above this level, however, chemical ozone depletion can be partially or completely masked by dynamical effects. Our calculations for 1992 showed that chemical loss was nearly completely compensated by increases due to diabatic descent at 655 K.

Manney, Gloria L.↗

Activated vortex lattice transition in a superconductor with combined sixfold and twelvefold anisotropic interactions

Abstract Numerical simulations are used to examine the transition dynamics between metastable and ground state vortex lattice phases in a system with combined sixfold and twelvefold contributions to the vortex–vortex interactions. The system is initially annealed using a twelvefold anisotropy, yielding domains of two different orientations and separated by grain boundaries. The vortex–vortex interaction is then suddenly changed to a sixfold anisotropy, rendering the twelvefold state metastable. Applying a drive that mimics an oscillating magnetic field causes the metastable state to decay, indicated by the structure factor that evolves from twelve to six peaks. The results fit the behavior seen in recent small-angle neutron scattering studies of the vortex lattice in MgB 2 . At higher drive amplitudes, the decay exhibits a two step process in which the initial fast decrease is followed by a slower regime where avalanches or bursts are correlated with dislocation annihilation events. The results are compared to other types of metastable systems with quenched disorder that decay under a periodic external drive.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Local parametric instability near elliptic points in vortex flows under shear deformation

The dynamics of two point vortices embedded in an oscillatory external flow consisted of shear and rotational components is addressed. The region associated with steady-state elliptic points of the vortex motion is established to experience local parametric instability. The instability forces the point vortices with initial positions corresponding to the steady-state elliptic points to move in spiral-like divergent trajectories. This divergent motion continues until the nonlinear effects suppress their motion near the region associated with the steady-state separatrices. The local parametric instability is then demonstrated not to contribute considerably to enhancing the size of the chaotic motion regions. Instead, the size of the chaotic motion region mostly depends on overlaps of the nonlinear resonances emerging in the perturbed system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Metastability and hysteretic vortex pinning near the order-disorder transition in NbSe 2 : Interplay between plastic and elastic energy barriers

We studied thermal and dynamic history effects in the vortex lattice (VL) near the order-disorder transition in clean NbSe 2 single crystals. Comparing the evolution of the effective vortex pinning and the bulk VL structure, we observed metastable superheated and supercooled VL configurations that coexist with a hysteretic effective pinning response due to thermal cycling of the system. A novel scenario, governed by the interplay between (lower) elastic and (higher) plastic energy barriers, is proposed as an explanation for our observations: Plastic barriers, which prevent the annihilation or creation of topological defects, require dynamic assistance to be overcome and to achieve a stable VL at each temperature. Conversely, thermal hysteresis in the pining response is ascribed to low energy barriers, which inhibit rearrangement within a single VL correlation volume and are easily overcome as the relative strength of competing interactions changes with temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermal structure and dynamics of the Jovian atmosphere. I - The Great Red Spot

Temperatures and thermal winds over the Great Red Spot (GRS) are derived from Voyager IRIS data for atmospheric pressures between 3 and 500 mbar. A cold tropopause implies a decrease in anticyclonic vorticity with height above 500 mbar through the lower stratosphere. Observations at 5 microns indicate that there is little emission from the GRS itself, although there is enhanced emission from a ring around it. The behavior of the tropopause and 5 micron temperatures are considered to be results of circulation which rises within the GRS and subsides in the area around it. Results are discussed in the context of several theoretical models of the GRS, and latent heat release as an energy source is considered. Dynamical scalings based on an axisymmetric frictionally controlled vortex indicate that the large-scale dynamics of the GRS are linear, as opposed to those of a tropical cyclone which are nonlinear.

Flasar, F. M.↗

Prediction of Concentrated Vortex Aerodynamics: Current CFD Capability Survey

Concentrated vortex flows contribute to the aerodynamic performance of aircraft at elevated load conditions. For military interests, the vortex flows are exploited at maneuver conditions of combat aircraft and missiles. For transport interests, the vortex flows are exploited at takeoff and landing conditions as well as at select transonic conditions. Aircraft applications of these vortex flows are reviewed with a historical perspective followed by a discussion of the underlying physics of a concentrated vortex flow. A hierarchy of computational fluid dynamics simulation technology is then presented followed by findings from a capability survey for predicting concentrated vortex flows with computational fluid dynamics. Results are focused on military and civil aircraft; only limited results are included for missiles. Opportunities for predictive capability advancement are then reported with comments related to digital transformation interests. A hierarchical approach that merges a physics-based perspective of the concentrated vortex flows with a systems engineering viewpoint of the air vehicle is also used to frame much of the discussion.

vortex flow↗

Interplay of waves and eddies in rotating stratified turbulence and the link with kinetic-potential energy partition

The interplay between waves and eddies in stably stratified rotating flows is investigated by means of world-class direct numerical simulations using up to 3072 3 grid points. Strikingly, we find that the shift from vortex- to wave-dominated dynamics occurs at a wave number kR which does not depend on the Reynolds number, suggesting that the partition of energy between wave and vortical modes is not sensitive to the development of turbulence at the smaller scales. Additionally, we also show that k R is comparable to the wave number at which exchanges between kinetic and potential modes stabilize at close to equipartition, emphasizing the role of potential energy, as conjectured in the atmosphere and the oceans. Moreover, k R varies as the inverse of the Froude number as explained by the scaling prediction proposed, consistently with recent observations and modeling of the Mesosphere–Lower Thermosphere and of the ocean.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Rotational transition, domain formation, dislocations, and defects in vortex systems with combined sixfold and twelvefold anisotropic interactions

We introduce a phenomenological model for a pairwise repulsive interaction potential of vortices in a type-II superconductor, consisting of superimposed sixfold and twelvefold anisotropies. Using numerical simulations we study how the vortex lattice configuration varies as the magnitudes of the two anisotropic interaction terms change. A triangular lattice appears for all values, and rotates through 30 ° as the ratio of the sixfold and twlevefold anisotropy amplitudes is varied, in agreement with experimental results. The transition causes the vortex lattice to split into domains that have rotated clockwise or counterclockwise, with grain boundaries that are “decorated” by dislocations consisting of fivefold and sevenfold coordinated vortices. We also find intradomain dislocations and defects, and characterize them in terms of their energy cost. We discuss how this model could be generalized to other particle-based systems with anisotropic interactions, such as colloids, and consider the limit of very large anisotropy where it is possible to create cluster crystal states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Local imaging of diamagnetism in proximity-coupled niobium nanoisland arrays on gold thin films.

In this work, we study the effect of engineered disorder on the local magnetic response of proximity-coupled superconducting island arrays by comparing scanning superconducting quantum interference device (SQUID) susceptibility measurements to a model in which we treat the system as a network of one-dimensional (1D) superconductor-normal-metal-superconductor Josephson junctions, each with a Josephson coupling energy E-J determined by the junction length or distance between islands. We find that the disordered arrays exhibit a spatially inhomogeneous diamagnetic response which, for low local applied magnetic fields, is well described by this junction network model, and we discuss these results as they relate to inhomogeneous 2D superconductors. Our model of the static magnetic response of the arrays does not fully capture the onset of nonlinearity and dissipation with increasing applied field, as these effects are associated with vortex motion due to the dynamic nature of the scanning SQUID susceptometry measurement. This work demonstrates a model 2D superconducting system with engineered disorder, and it highlights the impact of dissipation on the local magnetic properties of 2D superconductors and Josephson junction arrays.

Bishop-Van Horn, Logan↗

The effect of entrance flow development on vortex formation and wall shear stress in a curved artery model

In this work, we numerically investigate the effect of entrance condition on the spatial and temporal evolution of multiple three-dimensional vortex pairs and the wall shear stress distribution in a curved artery model. We perform this study using a Newtonian blood-analog fluid subjected to a pulsatile flow with two inflow conditions. The first flow condition is fully developed while the second condition is undeveloped (i.e., uniform). We discuss the connection along the axial direction between regions of organized vorticity observed at various cross sections of the model and compare results between the different entrance conditions. We model a human artery with a simple, rigid 180° curved pipe with a circular cross section and constant curvature, neglecting the effects of taper, torsion, and elasticity. Numerical results are computed from a discontinuous high-order spectral element flow solver. The flow rate used in this study is physiological. We observe differences in secondary flow patterns, especially during the deceleration phase of the physiological waveform where multiple vortical structures of both Dean-type and Lyne-type coexist. The results indicate that decreased axial velocities under an undeveloped condition produce smaller secondary flows that ultimately inhibit growth of any interior flow vortices. We highlight the effect of the entrance condition on the formation of these structures and subsequent appearance of abnormal inner wall shear stresses, which suggest there may be a lower prevalence of cardiovascular disease in curved arteries where the flow is rather undeveloped—a potentially physiologically significant result to help understand the influence of blood flow development on disease.

60 APPLIED LIFE SCIENCES↗