Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow velocity”

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 649 records · Page 36

The impact of swirl and wake strength on turbulent axisymmetric wake evolution

Here, an experimental investigation of swirl and wake strength influence on axisymmetric turbulent wake evolution was conducted. A novel wake generator design wire mounted in a wind tunnel test section with low free-stream turbulence produced wake Reynolds numbers based on momentum thickness and free-stream velocity in excess of 14000 and swirl numbers up to 0.4 with minimal blockage. Steady-state blade element momentum simulations of reference wind turbine designs indicated that wind turbines operate in the flow regimes studied, indicating the practical aspects of this work. Stereoscopic particle image velocimetry was used to acquire three components of velocity in the swirling wake at locations up to approximately ten diameters downstream. Quantitative measures of wake growth and decay were deduced using available equilibrium similarity scaling for the swirling wake. The results show an increase above 50% in growth and axial velocity decay rate constants over the range of swirl strength studied compared to those of the non-swirling wake. Tangential velocity decay constants were shown to decrease with swirl strength over the range of conditions studied. Notably, changes in wake strength have little influence on growth and decay rates when compared to changes in swirl strength for the flow regimes studied in this work.

17 WIND ENERGY↗

A Mixed Length Scale Model for Migrating Fluvial Bedforms

With the expansion of hydropower, in-stream converters, flood-protection infrastructures, and growing concerns on deltas fragile ecosystems, there is a pressing need to evaluate and monitor bedform sediment mass flux. It is critical to estimate real-time bedform size and migration velocity and provide a theoretical framework to convert easily accessible time histories of bed elevations into spatially evolving patterns. In this study, we collected spatiotemporally resolved bathymetries from laboratory flumes and the Colorado River in statistically steady, homogeneous, subcritical flow conditions. Wave number and frequency spectra of bed elevations show compelling evidence of scale-dependent velocity for the hierarchy of migrating bedforms observed in the laboratory and field. New scaling laws were applied to describe the full range of migration velocities as function of two dimensionless groups based on the bed shear velocity, sediment diameter, and water depth. Further simplification resulted in a mixed length scale model estimating scale-dependent migration velocities, without requiring bedform classification or identification.

58 GEOSCIENCES↗

Self-generated electrokinetic flows from active-charged boundary patterns

We develop a hydrodynamic description of self-generated electrolyte flow in capillaries whose bounding walls feature nonuniform distributions of charge nonuniform active ionic fluxes. The hydrodynamic velocity arising in such a system has components that are forbidden by symmetry in the absence of charge and fluxes. However, when these two boundary mechanisms are simultaneously present, they can lead to a symmetry broken state where steady flows with both unidirectional and circulatory components emerge. We show that these flow states arise when modulated boundary patterns of charge and fluxes are offset by a flux-charge phase difference, which is associated with the separation between sites of their peak densities on the wall. Mismatch in diffusivity of cationic and anionic species can modify the flow states and becomes an enhancing factor when fluxes of both ion species are being produced together at the same site. We demonstrate that this mechanism can be realized with a microfluidic generator that is powered by enzyme-coated patches that catalyze reactants in the solution to produce fluxes of ions. The local ionic elevation or depletion, which disrupts a nonuniform double layer, promotes self-induced gradients yielding persistent body forces to generate bulk fluid motion. Our work quantifies a boundary-driven mechanism behind self-sustained electrolyte flow in confined environments that exists without any external bulk-imposed fields or gradients. It provides a theoretical framework for understanding the combined effect of active and charged boundaries that are relevant in biological or soft matter systems, and can be utilized in electrofluidic and iontronic applications.

active matter↗

Experimental Measurements of the Wake of a Sphere at Subcritical Reynolds Numbers

This work experimentally investigated the flow phenomena and vortex structures in the wake of a sphere located in a water loop at Reynolds numbers of Re = 850, 1,250, and 1,700. Velocity fields in the wake region were obtained by applying the time-resolved stereoscopic particle image velocimetry (TR-SPIV) technique. From the acquired TR-SPIV velocity vector fields, the statistical values of mean and fluctuating velocities were computed. Spectral analysis, two-point velocity–velocity cross-correlation, proper orthogonal decomposition (POD) and vortex identification analyses were also performed. The velocity fields show a recirculation region that decreases in length with an increase of Reynolds numbers. The power spectra from the spectral analysis had peaks corresponding to a Strouhal number of St = 0.2, which is a value commonly found in the literature studies of flow over a sphere. The two-point cross-correlation analysis revealed elliptical structures in the wake, with estimated integral length scales ranging between 12% and 63% of the sphere diameter. The POD analysis revealed the statistically dominant flow structures that captured the most flow kinetic energy. It is seen that the flow kinetic energy captured in the smaller scale flow structures increased as Reynolds number increased. The POD modes contained smaller structure as the Reynolds number increased and as mode order increased. In addition, spectral analysis performed on the POD temporal coefficients revealed peaks corresponding to St = 0.2, similar to the spectral analysis on the fluctuating velocity. The ability of POD to produce low-order reconstructions of the flow was also utilized to facilitate vortex identification analysis, which identified average vortex sizes of 0.41D for Re1, 0.33D for Re2, and 0.32D for Re3.

Engineering↗

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗

A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity

In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.

42 ENGINEERING↗

Wind farm blockage effects: comparison of different engineering models

The work presents four engineering methods to estimate the induction zone in front of a wind turbine and account for the wind farm blockage effect. The methods comprise the vortex cylinder model, vortex dipole model, self-similar model, and wake projection model. The majority of the models presented account for yaw misalignments and ground effect. Actuator disk simulations are used to verify the individual models. The performance of each model is evaluated both in terms of precision and computational time. The induction models are coupled to wake models within the FLOw Redirection and Induction in Steady State framework to provide the full velocity field within a wind farm. Sample wind farm computations are presented, and the impact of including induction effects into wind farm performance predictions is reported. The different codes are publicly available online.

17 WIND ENERGY↗

Casimir interactions between two parallel graphene sheets carrying steady-state drift currents

Here, we investigate the fluctuation-induced Casimir interactions between two parallel graphene sheets carrying steady-state drift currents. The graphene properties are modeled based on the shifted Fermi disk model to capture the nonequilibrium optical response of the system. We find that the drift current introduces a repulsive correction to the perpendicular to the layers Casimir interaction, thereby reducing the overall attractive force. Although the correction is repulsive, it does not overcome the underlying attraction between the layers. It also generates a lateral force that opposes the carrier flow direction. Both contributions are studied in terms of distance and drift velocity functionalities showing pathways for Casimir force control.

Ke, Modi [University of South Florida, Tampa, FL (↗

Thermal chains and entrainment in cumulus updrafts, Part 1: Theoretical description

Recent studies have shown that cumulus updrafts often consist of a succession of discrete rising thermals with spherical vortex-like circulations. Herein, a theory is developed for why this “thermal chain” structure occurs. Theoretical expressions are obtained for a passive tracer, buoyancy, and vertical velocity in axisymmetric moist updrafts. Analysis of these expressions suggests that the thermal chain structure arises from enhanced lateral mixing associated with intrusions of dry environmental air below an updraft’s vertical velocity maximum. This dry air entrainment reduces buoyancy locally. Consequently, the updraft flow above levels of locally reduced buoyancy separates from below, leading to a breakdown of the updraft into successive discrete thermals. The range of conditions in which thermal chains exist is also analyzed from the theoretical expressions. A transition in updraft structure from isolated rising thermal, to thermal chain, to starting plume occurs with increases in updraft width, environmental relative humidity, and/or convective available potential energy. Corresponding expressions for the bulk fractional entrainment rate ε are also obtained. These expressions indicate rather complicated entrainment behavior of ascending updrafts, with local enhancement of ε up to a factor of ~2 associated with the aforementioned environmental air intrusions, consistent with recent large eddy simulation (LES) studies. These locally large entrainment rates contribute significantly to overall updraft dilution in thermal chain-like updrafts, while other regions within the updraft can remain relatively undilute. Part 2 of this study compares results from the theoretical expressions to idealized numerical simulations and LES.

54 ENVIRONMENTAL SCIENCES↗

Designing a Particle Flow Control Apparatus

Flow control within a particle-based Concentrated Solar Power (CSP) system is essential in determining the heat transfer coefficient, and therefore, the power generation capability of these systems. There are three areas where particle flow control is significant: the receivers, storage tanks, and particle-sCO2 heat exchangers. The focus of this work is on designing a new mechanism to control the flow in the particle-sCO2 heat exchangers due to the simplicity and potential cost savings when compared to the other areas of interest. The goal is for this new design to have quicker response times in terms of particle flowrate than a slide gate or flow control valve, which are designs currently used. The design resembles that of a chuck mechanism within a drill where a rotation of the sleeve elicits movement of the jaws both vertically and horizontally to close the outlet area of the nozzle. Additionally, this design will utilize the current actuator that is already used within these heat exchangers to reduce the complexity of implementation. The jaws are designed to be closed at an angle of 76˚ which is just slightly steeper than the hopper leading to the mechanism. Furthermore, this design can be tuned to limit particle bridging and other particle flow phenomena that result in blockages. The prototypes were 3D printed out of polylactic acid (PLA) and scaled up to 100%, 200%, and 400% to be able to observe the velocity profiles of the mechanism more clearly. Experiments are performed with this prototype to compare the inlet and outlet mass flow rates at different configurations of the jaws. The particles used in these experiments are 0.3mm HSP 40/70 that are commonly used in particle-based CSP systems.

14 SOLAR ENERGY↗

DESIGNING A PARTICLE FLOW CONTROL APPARATUS

Flow control within a particle-based Concentrated Solar Power (CSP) system is essential in determining the heat transfer coefficient, and therefore, the power generation capability of these systems. There are three areas where particle flow control is significant: the receivers, storage tanks, and particle-sCO2 heat exchangers. The focus of this paper is on designing a new mechanism to control the flow in the particle-sCO2 heat exchangers due to the simplicity and potential cost savings when compared to the other areas of interest. The goal is for this new design to have quicker response times in terms of particle flowrate than a slide gate or flow control valve, which are designs currently used. The design resembles that of a chuck mechanism within a drill where a rotation of the sleeve elicits movement of the jaws both vertically and horizontally to close the outlet area of the nozzle. Additionally, this design will utilize the current actuator that is already used within these heat exchangers to reduce the complexity of implementation. The jaws are designed to be closed at an angle of 76 ̊ which is just slightly steeper than the hopper leading to the mechanism. Furthermore, this design can be tuned to limit particle bridging and other particle flow phenomena that result in blockages. The prototypes were 3D printed out of polylactic acid (PLA) and scaled up to 100%, 200%, and 400% to be able to observe the velocity profiles of the mechanism more clearly. Experiments are performed with this prototype to compare the inlet and outlet mass flow rates at different configurations of the jaws. The particles used in these experiments are 0.3mm HSP 40/70 that are commonly used in particle-based CSP systems.

14 SOLAR ENERGY↗

Analysis of the Unsteady Wakes of Heavy Trucks in Platoon Formation and Their Potential Influence on Energy Savings

The authors present transient wind velocity measurements from two successive, well-documented truck platooning track-test campaigns to assess the wake-shedding behavior experienced by trucks in various platoon formations. Utilizing advanced analytics of data from fast-response (100-200-Hz) multi-hole pressure probes, this analysis examines aerodynamic flow features and their relationship to energy savings during close-following platoon formations. Applying Spectral analysis to the wind velocity signals, we identify the frequency content and vortex-shedding behavior from a forward truck trailer, which dominates the flow field encountered by the downstream trucks. The changes in dominant wake-shedding frequencies correlate with changes to the lead and follower truck fuel savings at short separation distances. All followers experience a reduction in wake-shedding frequency at following distances less than 20 m—the same zone where lead trucks experience fuel savings/drag reductions and where followers experience a reduction in platooning benefits. The dominant vortex-shedding frequency decreased by upwards of 35% in this region. Additionally, the third truck in a platoon is shown to experience a lower wake-shedding frequency than the second vehicle in either a two- or three-truck platoon. The correlation of changes in the vortex-shedding frequency to trends in drag reduction, and the associated energy savings, suggest that these phenomena may be related and may provide an opportunity for improving the efficiency potential of truck platooning.

47 OTHER INSTRUMENTATION↗

An immersed interface method for the 2D vorticity-velocity Navier-Stokes equations with multiple bodies

We present an immersed interface method for the vorticity-velocity form of the 2D Navier Stokes equations that directly addresses challenges posed by nonconvex immersed bodies, multiply connected domains, and the calculation of force distributions on immersed surfaces. The immersed interface method is re-interpreted as a polynomial extrapolation of flow quantities and boundary conditions into the immersed solid bodies, reducing computational cost and enabling simulations with nonconvex bodies that could not be discretized with previous immersed interface methods. In the flow, the vorticity transport equation is discretized using a conservative finite difference scheme and explicit Runge-Kutta time integration. The velocity reconstruction problem is transformed to a scalar Poisson equation that is discretized with conservative finite differences, and solved using an FFT-accelerated iterative algorithm. The use of conservative differencing throughout leads to exact enforcement of a discrete Kelvin's theorem, allowing for simulations with multiply connected domains and outflow boundaries that have challenged other immersed interface vortex methods. We also explore novel methods for recovering time-dependent pressure distributions on immersed bodies within a vorticity-based method and present a novel control volume formulation for recovering aerodynamic moments from only the vorticity and velocity fields. The method achieves second order spatial accuracy and third order temporal accuracy, and is validated on a variety of 2D flows in internal and free-space domains.

97 MATHEMATICS AND COMPUTING↗

Measurements of the Flow in the Vicinity of an Additively Manufactured Turbine Leading-Edge Using X-Ray Particle Tracking Velocimetry

X-ray particle tracking velocimetry (PTV) is performed, for the first time, to measure the velocity field inside a leading-edge of a turbine blade made by laser-additive-manufacturing (LAM) process. The traditional showerhead holes were replaced by a porous matrix in the leading-edge. The flow through such a leading-edge piece cannot be faithfully recreated by traditional prototype testing methods due to the surface roughness and imperfections caused by LAM process. Hence, direct measurement is the only option. However, it is difficult to measure flow inside such pieces with traditional velocimetry measurements due to the existence of metallic walls. Moreover, small internal size and high flow speeds call for a measurement technique with high spatial and temporal resolutions. To address these issues, we performed time-resolved X-ray PTV using the Advanced Photon Source (APS) synchrotron facility at the Argonne National Laboratory (ANL). A hydraulic system was constructed to run water, mixed with seeding particles, through the leading-edge piece. A high-speed camera captured the images of the seeding particles, which were later processed to create particle tracks. The time-averaged velocity field showed distinct pairs of vortices located in front of the porous outlet inside the leading-edge piece. The inlet channel showed reversed flow due to partial obstruction by the porous inlet of the test piece. Such knowledge of the flow field inside a leading-edge of a turbine blade will help us to design better cooling paths leading to higher cooling efficiency and increased life-span of a turbine blade.

42 ENGINEERING↗

A fourth order sharp immersed method for the incompressible Navier-Stokes equations with stationary and moving boundaries and interfaces

We propose a fourth order Navier-Stokes solver based on the immersed interface method (IIM), for flow problems with stationary and one-way coupled moving boundaries and interfaces. Our algorithm employs a Runge-Kutta-based projection method that maintains high-order temporal accuracy in both velocity and pressure for steady and unsteady velocity boundary conditions. Fourth order spatial accuracy is achieved through a novel fifth order IIM discretization scheme for the advection term, as well as existing high-order interface-corrected finite difference schemes for the other differential operators. Using a set of manufactured flow problems with stationary and moving boundaries, we demonstrate fourth order convergence of velocity and pressure in the infinity norm, both inside the domain and on the immersed boundaries. The solver’s performance is further validated through a range of practical flow simulations, highlighting its efficiency over a second order scheme. Finally, we showcase the ability of our immersed discretization scheme to handle interface-coupled multiphysics problems by solving a conjugate heat transfer problem with multiple immersed solids. Overall, the proposed approach robustly combines the efficiency of high order discretization schemes with the flexibility of immersed discretizations for flow problems with complex, moving boundaries and interfaces.

42 ENGINEERING↗

The Impact of Alfvénic Shear Flow on Magnetic Reconnection and Turbulence

Magnetic reconnection is a fundamental and omnipresent energy conversion process in plasma physics. Novel observations of fields and particles from Parker Solar Probe (PSP) have shown the absence of reconnection in a large number of current sheets in the near-Sun solar wind. Using near-Sun observations from PSP encounters 4–11 (2020 January–2022 March), we investigate whether reconnection onset might be suppressed by velocity shear. We compare estimates of the tearing mode growth rate in the presence of shear flow for time periods identified as containing reconnecting current sheets versus nonreconnecting times, finding systematically larger growth rates for reconnection periods. Upon examination of the parameters associated with reconnection onset, we find that 85% of the reconnection events are embedded in slow, non-Alfvénic wind streams. We compare with fast, slow non-Alfvénic, and slow Alfvénic streams, finding that the growth rate is suppressed in highly Alfvénic fast and slow wind, and reconnection is not seen in these wind types, as would be expected from our theoretical expressions. These wind streams have strong Alfvénic flow shear, consistent with the idea of reconnection suppression by such flows. This could help explain the frequent absence of reconnection events in the highly Alfvénic, near-Sun solar wind observed by PSP. Finally, we find a steepening of both the trace and magnitude magnetic field spectra within reconnection periods in comparison to ambient wind. We tie this to the dynamics of relatively balanced turbulence within these reconnection periods and the potential generation of compressible fluctuations.

slow solar wind↗

Velocity structure functions in multiphase turbulence: interpreting kinematics of Hα filaments in cool-core clusters

ABSTRACT The central regions of cool-core galaxy clusters harbour multiphase gas, with gas temperatures ranging from $10$ to $10^7\, \mathrm{K}$. Feedback from active galactic nuclei jets prevents the gas from undergoing a catastrophic cooling flow. However, the exact mechanism of this feedback energy input is unknown, mainly due to the lack of velocity measurements of the hot-phase gas. However, recent observations have measured the velocity structure functions (VSFs) of the cooler molecular (${\sim} 10\, \mathrm{K}$) and Hα filaments (${\sim} 10^4\, \mathrm{K}$) and used them to indirectly estimate the motions of the hot phase. In the first part of this study, we conduct high-resolution (3843–15363 resolution elements) simulations of homogeneous isotropic subsonic turbulence, without radiative cooling. We analyse the second-order velocity structure functions (VSF2) in these simulations and study the effects of varying spatial resolution, the introduction of magnetic fields, and the effect of projection along the line of sight (LOS) on it. In the second part of the study, we analyse high-resolution (7683 resolution elements) idealized simulations of multiphase turbulence in the intracluster medium from the companion study Mohapatra et al. We compare the VSF2 for both the hot ($T\sim 10^7\, \mathrm{K}$) and cold ($T\sim 10^4\, \mathrm{K}$) phases and find that their amplitude depends on the density contrast between the phases. They have similar scaling with separation, but introducing magnetic fields steepens the VSF2 of only the cold phase. We also find that projection along the LOS steepens the VSF2 for the hot phase and mostly flattens it for the cold phase.

Mohapatra, Rajsekhar (ORCID:0000000216007552)↗

Asymptotic coefficients of the attached-eddy model derived from an adiabatic atmosphere

The attached-eddy model (AEM) predicts that the mean streamwise velocity and streamwise velocity variance profiles follow a logarithmic shape, while the vertical velocity variance remains invariant with height in the overlap region of high Reynolds number wall-bounded turbulent flows. Moreover, the AEM coefficients are presumed to attain asymptotically constant values at very high Reynolds numbers. Here, the AEM predictions are examined using sonic anemometer measurements in the near-neutral atmospheric surface layer, with a focus on the logarithmic behaviour of the streamwise velocity variance. Utilizing an extensive 210-day dataset collected from a 62 m meteorological tower located in the Eastern Snake River Plain, Idaho, USA, the inertial sublayer is first identified by analysing the measured momentum flux and mean velocity profiles. The logarithmic behaviour of the streamwise velocity variance and the associated ‘−1’ scaling of the streamwise velocity energy spectra are then investigated. The findings indicate that the Townsend–Perry coefficient (A 1 ) is influenced by mild non-stationarity that manifests itself as a Reynolds number dependence. After excluding non-stationary runs, and requiring the bulk Reynolds number defined using the atmospheric boundary layer height to be larger than 4 × 10 7 , the inferred A 1 converges to values ranging between 1 and 1.25, consistent with laboratory experiments. Furthermore, nine benchmark cases selected through a restrictive quality control reveal a close relation between the ‘−1’ scaling in the streamwise velocity energy spectrum and the logarithmic behaviour of streamwise velocity variance. However, additional data are required to determine whether the plateau value of the pre-multiplied streamwise velocity energy spectrum is identical to A 1 .

atmospheric flows↗