Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vortex flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗
DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

What Controls the Entrainment Rate of Dry Buoyant Thermals with Varying Initial Aspect Ratio?

Here, this study uses theory and numerical simulations to analyze the nondimensional spreading rate a (change in radius with height) of buoyant thermals as they rise and entrain surrounding environmental fluid. A focus is on how a varies with initial thermal aspect ratio A r , defined as height divided by width of the initial buoyancy perturbation. An analytic equation for thermal ascent rate w t that depends on a is derived from the thermal-volume-averaged momentum bud get equation. The thermal top height when w t is maximum, defining a critical height z c , is inversely proportional to a. The height z c also corresponds to the thermal top height when buoyant fluid along the thermal’s vertical axis is fully replaced by entrained nonbuoyant environmental fluid rising from below the thermal. The time scale for this process is controlled by the vertical velocity of parcels rising upward through the thermal’s core. This parcel vertical velocity is approximated from Hill’s analytic spherical vortex, yielding an analytic inverse relation between a and A r . Physically, this a–A r relation is connected to changes in circulation as A r is modified. Numerical simulations of thermals with A r varied from 0.5 to 2 give a values close to the analytic theoretical relation, with a factor of ~3 decrease in a as A r is increased from 0.5 to 2. The theory also explains why a of initially spherical thermals from past laboratory and modeling studies is about 0.15. Overall, this study provides a theoretical underpinning for understanding the entrainment behavior of thermals, relevant to buoyantly driven atmospheric flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Wintertime Arctic Oscillation and North Atlantic Oscillation and their impacts on the Northern Hemisphere climate in E3SM

The characteristics of the wintertime Arctic Oscillation (AO) and North Atlantic Oscillation (NAO) and their impacts on climate variability over the Northern Hemisphere are important metrics for evaluating a climate system model. Observational analyses reveal that the horizontal and vertical structures in the AO and NAO exhibit a meridional dipole and a large-scale barotropic pattern between the Arctic and mid-latitudes. Historical model simulations from the Energy Exascale Earth System Model (E3SM-HIST) are used to identify how well it captures these major climate modes. It is found that the simulated AO and NAO modes have spatial structures similar to the observed features. In addition, the observed frequency bands in the AO and NAO-related time variability are captured well in the E3SM-HIST simulation. Associated with the positive phase in wintertime AO and NAO, zonal flow and warm advection in mid-latitude continents are enhanced, along with stronger cold flow from enhanced northerly winds over high latitudes. These features are linked to the atmospheric circulation pattern reflected by lower SLP anomalies over the Arctic and higher SLP anomalies over the mid-latitudes. In E3SM-HIST, these spatial associations and main structural features are analogous to those in observations. In the time-height evolution related to winter AO and NAO modes, it can also be seen that the simulations reproduce the downward propagating patterns in observations. Nevertheless, the vertical structures associated with AO and NAO in E3SM-HIST exhibit substantial biases in the lower stratosphere. The cause of these stratospheric biases is investigated using the strength of climatological stratospheric polar vortex (SPV) and wave activity fluxes. The results herein suggest that E3SM-HIST has a reasonable skill in reproducing the observed characteristics related to the winter AO and NAO, although there exist systematic biases in the associated climate variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Lagrangian dynamics of particle transport in oral and nasal breathing

Here, we present a large-eddy simulation (LES) of saliva particle transport during normal human breathing through the nose and mouth. The flow of the air–saliva mixture is modeled using an Eulerian LES that is coupled with a Lagrangian particle tracking module to obtain trajectories of saliva particles in a room with stagnant air conditions. The coupled Eulerian–Lagrangian simulation yields novel insights into the intricate dynamics of Lagrangian coherent structures (LCS) and fundamental material lines that emerge from the saliva particles' trajectories during several breathing cycles. Further, we systematically compare the quantitative LCS diagnostics of mouth breathing with those of mouth and nose normal breathing. Analyzing the simulation results of human breathing from the mouth and nose, we show that, soon after the first breathing cycle, saliva particles form a series of roll-up vortex rings that propagate forward. The forward propagation of these vortex rings leads to the formation of an asymmetrical primary forefront vortex. The individual vortex rings continuously propagate forward, merging with the forefront vortex, and ascending along the limb of the leading vortex.

42 ENGINEERING↗
DOE OSTI2022

Experimental pub crawl from Rayleigh–Bénard to magnetostrophic convection

The interplay between convective, rotational and magnetic forces defines the dynamics within the electrically conducting regions of planets and stars. Yet their triadic effects are separated from one another in most studies, arguably due to the richness of each subset. In a single laboratory experiment, we apply a fixed heat flux, two different magnetic field strengths and one rotation rate, allowing us to chart a continuous path through Rayleigh–Bénard convection (RBC), two regimes of magnetoconvection, rotating convection and two regimes of rotating magnetoconvection, before finishing back at RBC. Dynamically rapid transitions are determined to exist between jump rope vortex states, thermoelectrically driven magnetoprecessional modes, mixed wall- and oscillatory-mode rotating convection and a novel magnetostrophic wall mode. Thus, our laboratory ‘pub crawl’ provides a coherent intercomparison of the broadly varying responses arising as a function of the magnetorotational forces imposed on a liquid-metal convection system.

42 ENGINEERING↗
DOE OSTI2023

Nonlinear mode coupling and energetics of driven magnetized shear-flow turbulence

To comprehensively understand the saturation of two-dimensional (2D) magnetized Kelvin–Helmholtz-instability-driven turbulence, energy transfer analysis is extended from the traditional interaction between scales to include eigenmode interactions, by using the nonlinear couplings of linear eigenmodes of the ideal instability. While both kinetic and magnetic energies cascade to small scales, a significant fraction of turbulent energy deposited by unstable modes in the fluctuation spectrum is shown to be re-routed to the conjugate-stable modes at the instability scale. They remove energy from the forward cascade at its inception. The remaining cascading energy flux is shown to attenuate exponentially at a small scale, dictated by the large-scale stable modes. Guided by a widely used instability-saturation assumption, a general quasi-linear model of instability is tested by retaining all nonlinear interactions except those that couple to the large-scale stable modes. These complex interactions are analytically removed from the magnetohydrodynamic equations using a novel technique. Observations are an explosive large-scale vortex separation instead of the well-known merger of 2D, a dramatic enhancement in turbulence level and spectral energy fluxes, and a reduced small-scale dissipation length scale. These show the critical role of the stable modes in instability saturation. Possible reduced-order turbulence models are proposed for fusion and astrophysical plasmas, based on eigenmode-expanded energy transfer analyses.

Physics↗
DOE OSTI2023

Data-Driven Transition Path Analysis Yields a Statistical Understanding of Sudden Stratospheric Warming Events in an Idealized Model

Abstract Atmospheric regime transitions are highly impactful as drivers of extreme weather events, but pose two formidable modeling challenges: predicting the next event (weather forecasting) and characterizing the statistics of events of a given severity (the risk climatology). Each event has a different duration and spatial structure, making it hard to define an objective “average event.” We argue here that transition path theory (TPT), a stochastic process framework, is an appropriate tool for the task. We demonstrate TPT’s capacities on a wave–mean flow model of sudden stratospheric warmings (SSWs) developed by Holton and Mass, which is idealized enough for transparent TPT analysis but complex enough to demonstrate computational scalability. Whereas a recent article (Finkel et al. 2021) studied near-term SSW predictability, the present article uses TPT to link predictability to long-term SSW frequency. This requires not only forecasting forward in time from an initial condition, but also backward in time to assess the probability of the initial conditions themselves. TPT enables one to condition the dynamics on the regime transition occurring, and thus visualize its physical drivers with a vector field called the reactive current . The reactive current shows that before an SSW, dissipation and stochastic forcing drive a slow decay of vortex strength at lower altitudes. The response of upper-level winds is late and sudden, occurring only after the transition is almost complete from a probabilistic point of view. This case study demonstrates that TPT quantities, visualized in a space of physically meaningful variables, can help one understand the dynamics of regime transitions.

Meteorology & Atmospheric Sciences↗
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↗
DOE OSTI2025

Dynamics of Downdrafts Around a Growing Convective Cloud: A Numerical Study

We examine the dynamics of cloud-edge downdrafts over the growth phase of isolated cumuli, combining Eulerian and Lagrangian analyses. As in previous studies, our results show that growing cumuli are surrounded by downdrafts linked to cloud-scale quasi-toroidal circulations at all times at middle and upper cloud levels consistent with the thermal chain description of convective clouds. These toroidal circulations are responsible for the most intense cloud-edge downdrafts in our simulations. In the upper cloud half, roughly 30%–50% of the upward mass flux is typically compensated within a radius of about twice the updraft radius in quasi-laminar simulations forced by a warm bubble in an initially quiescent flow. In a turbulent cloud forced by surface fluxes, this compensation fraction is around 10%–30% over the same region. In contrast to the buoyancy-centered view of subsiding shells, Eulerian and Lagrangian vertical momentum budget analyses show that the most intense cloud-edge downdrafts in the turbulent setup, and after spin-up of the toroidal circulation in the quasi-laminar experiments, are predominantly mechanically forced (i.e., driven by dynamic pressure accelerations). This is consistent throughout the entire growth phase of the cumulus clouds and across tests with varying assumptions, including drier and moister environments. Despite dynamic pressure perturbations being the main driver of toroidal downdrafts, the downdraft speed (relative to the corresponding updraft velocity) exceeds the prediction of the non-buoyant Hill's spherical vortex—a simple model frequently used for cloud-scale circulations—by more than 30%.

Pardo, Lianet Hernández [Goethe Univ., Frankfurt (↗
DOE OSTI2021

Influence of Heave Plate Topology on Reaction Force

In this study, we explore the role of heave plate topology on fluid reaction forces using three, representative shapes: a hexagonal flat plate, a hexagonal conic with an open top, and the same with a closed top that encloses a fluid mass. We force each test article sinusoidally in a quiescent tank and decompose the reaction force using forms of the Morison equation for phase-invariant and phase-dependent parameterizations. We find that a flat plate generates 5.3 % more fluid reaction force than the open conic topology, and 21.4% more than the enclosed conic. Asymmetric topologies generate asymmetric reaction forces, but the magnitude of asymmetry is limited by nearly symmetric fluid inertia forces, which dominate over drag for these test articles. Additionally, we observe asymmetric vortex dynamics for the flat plate when the Keulegan-Carpenter number is between 1 and 2, accompanied by a shift in the phase of the peak force by of the oscillation period. As a consequence of this shift, the hydrodynamic coefficients estimated from the phase-dependent Morison equation decomposition are asymmetric, suggesting that phase-dependent representations may not provide physical insight in some hydrodynamic regimes.

added mass↗
DOE OSTI2026

Emergence of vorticity and viscous stress in finite-scale quantum hydrodynamics

The Madelung equations offer a hydrodynamic description of quantum systems, from single particles to quantum fluids. In this formulation, the probability density is mapped onto the fluid density and the phase is treated as a scalar potential generating the velocity field. As examples of potential flows, quantum fluids described in this way are inherently irrotational, but quantum vortices may arise at discrete points where the phase is undefined. In this paper, starting from this irrotational description of a quantum fluid, a coarse-graining procedure is applied to arrive at a macroscopic description of the quantum fluid in terms of a hierarchy of moments in which the role of velocity is played by a Favre average of the microscopic velocity field. This hierarchy is truncated using an explicit closure derived from an expansion in a finite length scale. The resulting coarse-grained fields are shown to allow for finite vorticity at any point in the fluid. Additionally, it is shown that this vorticity obeys a similar equation to the vorticity equation in classical hydrodynamics and includes a vortex-stretching term. The particular closure employed here also gives rise to a novel stress term in the fluid equations, which in the appropriate limit appears analogous to an artificial viscous stress from computational fluid dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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

Experimental study of vorticity-strain interactions in turbulent premixed counterflow flames

Here, the effects of heat release on interactions between vorticity (ω) and strain rate (s) in turbulent premixed CH 4 /O 2 /N 2 counterflow flames are investigated using simultaneous OH laser-induced fluorescence (LIF) and tomographic particle image velocimetry (TPIV) measurements. A comparison between the flames and a corresponding turbulent non-reacting variable density N 2 -vs-products counterflow reveals the impact of heat release on vorticity-strain rate alignment statistics. Vorticity and strain rate statistics in the flames and non-reacting flow are conditioned on distance from the local flame front and gas mixing layer interface (GMLI) contours, respectively. The magnitude, alignment, and spatial distribution of the vorticity and principal strain rates (s 1 , s 2 , s 3 ) are rather different when heat release is present. Density variations without heat release enhance the ω-s 2 alignment while significantly reducing the ω-s 3 alignment and modestly reducing the ω-s 1 alignment. In contrast, heat release at the flame front further reduces the ω-s 1 alignment but increases the ω-s 3 alignment and suppresses the preferential ω-s 2 alignment. Furthermore, increasing turbulence diminishes the effect of heat release on this preferential alignment. In regions with the largest vorticities, both the reacting and non-reacting counterflows show an increase in the probability of ω-s 2 alignment. All counterflow cases have a net positive vortex-stretching contribution to the enstrophy production with a peak production rate at the flame front or GMLI, but the peak values depend on the density variation, heat release, and turbulence level. Elucidation of the complex interplay between these factors contributes to the understanding of the dynamics of turbulence-flame interactions.

42 ENGINEERING↗
DOE OSTI2021

The Effect of the Degree of Premixedness on Self-Excited Combustion Instability

The use of lean, premixed fuel and air mixtures is a common strategy to reduce NO x emissions in gas turbine combustors. However, this strategy causes an increased susceptibility to self-excited instability, which manifests as fluctuations in heat release rate, flow velocity, and combustor acoustics that oscillate in-phase in a feedback loop. Here, this study considers the effect of the level of premixedness on the self-excited instability in a single-nozzle combustor. In this system, the fuel can be injected inside the nozzle to create a partially-premixed mixture or far upstream to create a fully-premixed mixture, varying the level of premixedness of the fuel and air entering the combustor. When global equivalence ratio is held constant, the cases with higher levels of premixing have higher instability amplitudes. High-speed CH* chemiluminescence imaging shows that the flame for these cases is more compact and the distribution of the heat release rate oscillations is more concentrated near the corner of the combustor in the outer recirculation zones. Rayleigh index images, which are a metric for the relative phase of pressure and heat release rate oscillations, suggest that vortex rollup in the corner region is primarily responsible for determining instability characteristics when premixedness is varied. This finding is further supported through analysis of local flame edge dynamics.

42 ENGINEERING↗
DOE OSTI2021

Blood Plasma Self-Separation Technologies during the Self-Driven Flow in Microfluidic Platforms

Blood plasma is the most commonly used biofluid in disease diagnostic and biomedical analysis due to it contains various biomarkers. The majority of the blood plasma separation is still handled with centrifugation, which is off-chip and time-consuming. Therefore, in the Lab-on-a-chip (LOC) field, an effective microfluidic blood plasma separation platform attracts researchers’ attention globally. Blood plasma self-separation technologies are usually divided into two categories: active self-separation and passive self-separation. Passive self-separation technologies, in contrast with active self-separation, only rely on microchannel geometry, microfluidic phenomena and hydrodynamic forces. Passive self-separation devices are driven by the capillary flow, which is generated due to the characteristics of the surface of the channel and its interaction with the fluid. Comparing to the active plasma separation techniques, passive plasma separation methods are more considered in the microfluidic platform, owing to their ease of fabrication, portable, user-friendly features. We propose an extensive review of mechanisms of passive self-separation technologies and enumerate some experimental details and devices to exploit these effects. The performances, limitations and challenges of these technologies and devices are also compared and discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Conjugate heat transfer↗
DOE OSTI2022

Hydrokinetic energy conversion using flow induced oscillations of single-cylinder with large passive turbulence control

Various types of flow-induced oscillations (FIOs) have been implemented in development of marine hydrokinetic (MHK) energy converters. With passive turbulence control (PTC), energy harvesting starts at a flow speed of about 0.5 m/s. However, there is worldwide MHK energy available in even slower currents. In the present study, the effect of damping on FIO and power extraction is investigated for a converter with large turbulence stimulation (PTC) consisting of straight strips with a height of 15% of the cylinder diameter and placed symmetrically on the cylinder surface. The oscillating amplitude decreases, as the damping ratio increases, with unchanged sinusoidal pattern of the displacement time-history. The frequency ratio is also affected by damping especially in the VIV initial branch and transition region between VIV and galloping. An important flow characteristic of the large-PTC cylinder is that a recirculation region is formed behind the PTC, causing appreciable disturbance to the flow past the cylinder. Power can be harvested in the whole FIO range and the harnessed power maximum appears at the largest inflow velocity tested. However, the optimum of harnessing efficiency is located at the beginning of the VIV upper branch. The gap between VIV and galloping is bridged when large PTC is used, eliminating the drop in power and efficiency even at higher damping, which would be a weakness of regular-PTC cylinder for energy harvesting. Finally, the mechanism behind the variation of harnessing efficiency with inflow velocity and damping ratio is revealed, and the optimality criterion for the converter design is discussed.

16 TIDAL AND WAVE POWER↗