Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vortex methods”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗
DOE OSTI2024

Wind turbine performance enhancement with minimal structural load penalty: A design philosophy

The performance benefits of using tip devices on wind turbines has been well-documented. However, previous studies show that adding blade tip devices such as winglets leads to a significant increase in blade root bending moment, potentially requiring structural reinforcement with cost and weight drawbacks. A new and unique design philosophy for retrofit blade tip devices for wind turbines is presented. By balancing generated aerodynamic and centrifugal loads, these devices offer an increase in power production without the need for structural reinforcement. Predicted performance and cost benefits of using retrofit blade tip devices on the National Renewable Energy Laboratory 5 MW reference wind turbine are shown. Finally, the addition of blade tip devices resulted in significant improvements in the coefficient of power (C p ) and annual energy production (AEP).

17 WIND ENERGY↗
DOE OSTI2022

Accelerated Wind-Turbine Wake Recovery Through Actuation of the Tip-Vortex Instability

Advances in wind-plant control have often focused on more effectively balancing power between neighboring turbines. Wake steering is one such method that provides control-based improvements in a quasi-static way, but this does little to fundamentally change the wake recovery process, and thus, it has limited potential. This study investigates use of another control paradigm known as dynamic wake control (DWC) to excite the mutual inductance instability between adjacent tip-vortex structures, thereby accelerating the breakdown of the structures. The current work carries this approach beyond the hypothetical by applying the excitation via turbine control vectors that already exist on all modern wind turbines: blade pitch and rotor speed control. The investigation leverages a free-vortex wake method (FVWM) that allows a thorough exploration of relevant frequencies and amplitudes of harmonic forcing for each control vector (as well as the phase difference between the vectors for a tandem configuration) while still capturing the essential tip-vortex dynamics. The FVWM output feeds into a Fourier stability analysis working to pinpoint candidate DWC strategies suggesting fastest wake recovery. Near-wake length reductions of >80% are demonstrated, although without considering inflow turbulence. Analysis is provided to interpret these predictions considering the presence of turbulence in a real atmospheric inflow.

17 WIND ENERGY↗
DOE OSTI2025

Unsteady aerodynamics of large-scale floating offshore wind turbines in surge motion

Unsteady aerodynamic loads significantly influence the design and wake flow field of floating offshore wind turbines, especially due to wave- or vibration-induced tower top motions triggering various unsteady phenomena. Recent studies show that increasing turbine sizes amplify unsteady aerodynamic effects, as their impact typically grows with rotor diameter. This work combines recent findings from experiments and simulations on model-scale FOWT aerodynamics with new numerical analyses of large-scale rotors, providing a comprehensive understanding of unsteady phenomena occurrence and impact. Numerical analyses of the IEA 15-MW and 22-MW rotors undergoing surge motions characterise the combined influence of motion-induced unsteady phenomena on rotor thrust. Results indicate that unsteady effects can reduce thrust force variations by up to 40% at realistic surge periods. These findings contrast with prior model-scale rotor investigations, attributed to the specialised design of the model rotors. Comparisons between numerical methods - the dynamic blade element momentum method in OpenFAST and free vortex wake modules in panMARE and OpenFAST - reveal persistent differences in thrust predictions under both idealised and realistic conditions. This highlights that classical blade element momentum approaches require enhancement to accurately capture unsteady loads on large-scale floating offshore wind turbines at low wind speeds.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Structure of a chaotic tangle of quantum vortices in turbulent superfluid liquids and in a Bose – Einstein condensate

Based on the theory of the thermodynamic equilibrium in a system of quantum vortices in superfluid liquids and in a Bose – Einstein condensate in the presence of a counterflow of normal and superfluid components, we investigate the structure of a chaotic tangle of quantum vortices in turbulent superfluid liquids. Using the characteristic functional method, the properties of hydrodynamic vortex filaments are examined. It is shown that the average curvature of the vortex lines is on the order of the intervortex distance, with the proportionality coefficient being independent of the counterflow velocity. It is found that the degree of anisotropy of the vortex loops does not depend on the applied counterflow velocity. The obtained results explain the origin of the anisotropy as well as the relationship between the curvature of the lines and the intervortex space and their dependence on the parameters of the problem. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A vortex sheet based analytical model of the curled wake behind yawed wind turbines

Motivated by the need for compact descriptions of the evolution of non-classical wakes behind yawed wind turbines, we develop an analytical model to predict the shape of curled wakes. Interest in such modelling arises due to the potential of wake steering as a strategy for mitigating power reduction and unsteady loading of downstream turbines in wind farms. We first estimate the distribution of the shed vorticity at the wake edge due to both yaw offset and rotating blades. By considering the wake edge as an ideally thin vortex sheet, we describe its evolution in time moving with the flow. Vortex sheet equations are solved using a power series expansion method, and an approximate solution for the wake shape is obtained. The vortex sheet time evolution is then mapped into a spatial evolution by using a convection velocity. Apart from the wake shape, the lateral deflection of the wake including ground effects is modelled. Our results show that there exists a universal solution for the shape of curled wakes if suitable dimensionless variables are employed. For the case of turbulent boundary layer inflow, the decay of vortex sheet circulation due to turbulent diffusion is included. Finally, we modify the Gaussian wake model by incorporating the predicted shape and deflection of the curled wake, so that we can calculate the wake profiles behind yawed turbines. Model predictions are validated against large-eddy simulations and laboratory experiments for turbines with various operating conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Identification of Kelvin-Helmholtz generated vortices in magnetised fluids

The Kelvin-Helmholtz Instability (KHI), arising from velocity shear across the magnetopause, plays a significant role in the viscous-like transfer of mass, momentum, and energy from the shocked solar wind into the magnetosphere. While the KHI leads to growth of surface waves and vortices, suitable detection methods for these applicable to magnetohydrodynamics (MHD) are currently lacking. A novel method is derived based on the well-established λ-family of hydrodynamic vortex identification techniques, which define a vortex as a local minimum in an adapted pressure field. The J × B Lorentz force is incorporated into this method by using an effective total pressure in MHD, including both magnetic pressure and a pressure-like part of the magnetic tension derived from a Helmholtz decomposition. The λ MHD method is shown to comprise of four physical effects: vortical momentum, density gradients, fluid compressibility, and the rotational part of the magnetic tension. A local three-dimensional MHD simulation representative of near-flank magnetopause conditions (plasma β’s 0.5 – 5 and convective Mach numbers M f ∼ 0.4) under northward interplanetary magnetic field (IMF) is used to validate λ MHD . Analysis shows it correlates well with hydrodynamic vortex definitions, though the level of correlation decreases with vortex evolution. Overall, vortical momentum dominates λ MHD at all times. During the linear growth phase, density gradients act to oppose vortex formation. By the highly nonlinear stage, the formation of small-scale structures leads to a rising importance of the magnetic tension. Compressibility was found to be insignificant throughout. Finally, a demonstration of this method adapted to tetrahedral spacecraft observations is performed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

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↗
DOE OSTI2022

Linear stability analysis via simulated annealing and accelerated relaxation

Simulated annealing (SA) is a kind of relaxation method for finding equilibria of Hamiltonian systems. A set of evolution equations is solved with SA, which is derived from the original Hamiltonian system so that the energy of the system changes monotonically while preserving Casimir invariants inherent to noncanonical Hamiltonian systems. The energy extremum reached by SA is an equilibrium. Since SA searches for an energy extremum, it can also be used for stability analysis when initiated from a state where a perturbation is added to an equilibrium. The procedure of the stability analysis is explained, and some examples are shown. Because the time evolution is computationally time consuming, efficient relaxation is necessary for SA to be practically useful. An acceleration method is developed by introducing time dependence in the symmetric kernel used in the double bracket, which is part of the SA formulation described here. An explicit formulation for low-beta reduced magnetohydrodynamics (MHD) in cylindrical geometry is presented. In conclusion, since SA for low-beta reduced MHD has two advection fields that relax, it is important to balance the orders of magnitude of these advection fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Quantitative and interpretable order parameters for phase transitions from persistent homology

Here, we apply modern methods in computational topology to the task of discovering and characterizing phase transitions. As illustrations, we apply our method to four two-dimensional lattice spin models: the Ising, square ice, XY, and fully frustrated XY models. In particular, we use persistent homology, which computes the births and deaths of individual topological features as a coarse-graining scale or sublevel threshold is increased, to summarize multiscale and high-point correlations in a spin configuration. We employ vector representations of this information called persistence images to formulate and perform the statistical task of distinguishing phases. For the models we consider, a simple logistic regression on these images is sufficient to identify the phase transition. Interpretable order parameters are then read from the weights of the regression. This method suffices to identify magnetization, frustration, and vortex-antivortex structure as relevant features for phase transitions in our models. We also define “persistence” critical exponents and study how they are related to those critical exponents usually considered.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Challenges and transformative opportunities in superconductor vortex physics

In superconductors, the motion of vortices introduces unwanted dissipation that is disruptive to applications. Fortunately, material defects can immobilize vortices, acting as vortex pinning centers, which engenders dramatic improvements in superconductor material properties and device operation. This has motivated decades of research into developing methods of tailoring the disorder landscape in superconductors to increase the strength of vortex pinning. Yet, efficacious materials engineering still eludes us. The electromagnetic properties of real (disordered) superconducting materials cannot yet be reliably predicted, such that designing superconductors for applications remains a largely inefficient process of trial and error. This is ultimately due to large gaps in our knowledge of vortex dynamics: the field is challenged by the extremely complex interplay between vortex elasticity, vortex-vortex interactions, and material disorder. Here, we review obstacles and recent successes in understanding and controlling vortex dynamics in superconducting materials and devices. We further identify major open questions and discuss opportunities for transformative research in the field. This includes improving our understanding of vortex creep, determining and reaching the ceiling for the critical current, advanced microscopy to garner accurate structure-property relationships, frontiers in predictive simulations and the benefits of artificial intelligence, as well as controlling and exploiting vortices in quantum information applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Effect of slip on vortex shedding from a circular cylinder in a gas flow

Most studies of vortex shedding from a circular cylinder in a gas flow have explicitly or implicitly assumed that the no-slip condition applies on the cylinder surface. To investigate the effect of slip, vortex shedding is simulated using molecular gas dynamics (the direct simulation Monte Carlo method) and computational fluid dynamics (the incompressible Navier-Stokes equations with a slip boundary condition). A Reynolds number of 100, a Mach number of 0.3, and a corresponding Knudsen number of 0.0048 are examined. For these conditions, compressibility effects are small, and periodic laminar vortex shedding is obtained. Slip on the cylinder is varied using combinations of diffuse and specular molecular reflections with accommodation coefficients from zero (maximum slip) to unity (minimum slip). Although unrealistic, bounce-back molecular reflections are also examined because they approximate the no-slip boundary condition (zero slip). The results from both methods are in reasonable agreement. The shedding frequency increases slightly as the accommodation coefficient is decreased, and shedding ceases at low accommodation coefficients (large slip). The streamwise and transverse forces decrease as the accommodation coefficient is decreased. Based on the good agreement between the two methods, computational fluid dynamics is used to determine the critical accommodation coefficient below which vortex shedding ceases for Reynolds numbers of 60–100 at a Mach number of 0.3. Finally, conditions to observe the effect of slip on vortex shedding appear to be experimentally realizable, although challenging.

42 ENGINEERING↗
DOE OSTI2020

Assessing the blockage effect of wind turbines and wind farms using an analytical vortex model

Wind farm blockage effects are currently neglected in the prediction of wind farm energy yield, typically leading to an overestimation of the production. This work presents a novel method to assess wind farm production, while accounting for blockage effects. We apply a vortex model, based on a cylindrical wake, to assess induction effects. We present variations of the model to account for finite wake length, finite tip-speed ratios, and the proximity to the ground. The results are applied to single rotors in aligned and yawed conditions and to different wind farm layouts. We provide far-field approximations for faster estimates of the velocity field. Further, this article includes a new methodology to couple the induction model to engineering wake models, such as the ones present in the FLOw Redirection and Induction in Steady State (FLORIS). We compare the results to actuator disk simulations for various operating conditions of a single turbine and different wind farm layouts. We found that the mean relative error of the model in the induction zone is typically around 0.2% compared with actuator disk simulations. The computational time of the velocity field using the analytical vortex model is three orders of magnitude less than the one obtained with the actuator disk simulation.

17 WIND ENERGY↗
DOE OSTI2020

Reconfigurable Pinwheel Artificial-Spin-Ice and Superconductor Hybrid Device

The ability to control the potential landscape in a medium of interacting particles could lead to intriguing collective behavior and innovative functionalities. In this work, we utilize spatially reconfigurable magnetic potentials of a pinwheel artificial-spin-ice (ASI) structure to tailor the motion of superconducting vortices. The reconstituted chain structures of the magnetic charges in the pinwheel ASI and the strong interaction between magnetic charges and superconducting vortices allow significant modification of the transport properties of the underlying superconducting thin film, resulting in a reprogrammable resistance state that enables a reversible and switchable vortex Hall effect. Our results highlight an effective and simple method of using ASI as an in situ reconfigurable nanoscale energy landscape to design reprogrammable superconducting electronics, which could also be applied to the in situ control of properties and functionalities in other magnetic particle systems, such as magnetic skyrmions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Mechanism for Anomalous Hall Ferromagnetism in Twisted Bilayer Graphene

Motivated by the recent observation of an anomalous Hall effect in twisted bilayer graphene, here we use a lowest Landau level model to understand the origin of the underlying symmetry-broken correlated state. This effective model is rooted in the occurrence of Chern bands which arise due to the coupling between the graphene device and its encapsulating substrate. Our model exhibits a phase transition from a spin-valley polarized insulator to a partial or fully valley unpolarized metal as the bandwidth is increased relative to the interaction strength, consistent with experimental observations. In sharp contrast to standard quantum Hall ferromagnetism, the Chern number structure of the flat bands precludes an instability to an intervalley coherent phase, but allows for an excitonic vortex lattice at large interaction anisotropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Time-dependent Ginzburg-Landau treatment of rf magnetic vortices in superconductors: Vortex semiloops in a spatially nonuniform magnetic field

We apply time-dependent Ginzburg-Landau (TDGL) numerical simulations to study the finite frequency electrodynamics of superconductors subjected to an intense rf magnetic field. Much recent TDGL work has focused on spatially uniform external magnetic fields and largely ignores the Meissner state screening response of the superconductor. In this paper, we solve the TGDL equations for a spatially nonuniform magnetic field created by a point magnetic dipole in the vicinity of a semi-infinite superconductor. A two-domain simulation is performed to accurately capture the effect of the inhomogeneous applied fields and the resulting screening currents. The creation and dynamics of vortex semiloops penetrating deep into the superconductor domain are observed and studied, and the resulting third-harmonic nonlinear response of the sample is calculated. Further, the effect of pointlike defects on vortex semi-loop behavior is also studied. This simulation method will assist our understanding of the limits of superconducting response to intense rf magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗