Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence 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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sensitivity of fatigue reliability in wind turbines: effects of design turbulence and the Wöhler exponent

Fatigue assessment of wind turbines involves three main sources of uncertainty: material resistance, load, and the damage accumulation model. Many studies focus on increasing the accuracy of fatigue load assessment to improve the fatigue reliability. Probabilistic modeling of the wind's turbulence standard deviation is an example of an approach used for this purpose. Editions 3 and 4 of the IEC standard for the design of wind energy generation systems (IEC 61400-1) suggest different probability distributions as alternatives for the representative turbulence in the normal turbulence model (NTM) of edition 1. There are debates on whether the suggested distributions provide conservative reliability levels, as the established design safety factors are calibrated based on the representative turbulence approach. The current study addresses the debate by comparing annual reliability based on different scenarios of NTM using a probabilistic approach. More importantly, it elaborates on the relative importance of load assessment accuracy in defining the fatigue reliability. Using the DTU 10 MW reference wind turbine and the first-order reliability method (FORM), we study the changes in the annual reliability level and its sensitivity to the three main random inputs. We perform the study considering the blade root flapwise and the tower base fore–aft moments, assuming different fatigue exponents in each load channel. The results show that integration over distributions of turbulence in each mean wind speed results in less conservative annual reliability levels than representative turbulence. The difference in the reliability levels varies according to turbulence distribution and the fatigue exponent. In the case of the tower base, the difference in the annual reliability index after 20 years can be up to 50 %. However, the model and material uncertainty have much higher effects on the reliability levels compared to load uncertainty. Knowledge about such differences in the reliability levels due to the choice of turbulence distribution is especially important, as it impacts the extent of lifetime extension through reliability reassessments.

17 WIND ENERGY↗
DOE OSTI2025

Directional Variations in Tidal Flow Multifractality and Intermittency

We explored distinct directional variations of multifractal and intermittent characteristics of ebb and flood flow velocities at nodule point, WA, tidal energy site, and complementary inspection on the East River, NY using scaling exponents of the structure function, distribution flatness, detrending moving average (DMA) analysis, multifractal detrended fluctuation analysis (MF-DFA), and high-order spectral moments. Our findings reveal that tidal flow presents higher long-range dependence (LRD) and intermittent levels for the ebb flow sections, whereas the flood flow sections exhibit a higher degree of multifractality and greater sensitivity to larger magnitude of turbulent fluctuations. We demonstrate that long-range dependence predominantly contributes to multifractal behavior in both ebb and flood flows, as evidenced by the significantly reduced multifractal spectrum width for temporally randomly permuted time series. Moreover, spectral kurtosis analysis uncovers a higher intermittent level across all frequency scales for ebb flow sections and reveals a distinct pattern of tidal flow intermittency differing from the monotonically increasing intermittent level observed in wall-bounded and grid turbulence. Finally, we demonstrate that DMA, MF-DFA, and high-order spectral moments provide more comprehensive insights than structure function scaling exponents and PDF flatness methods.

Cheng, Shyuan [Univ. of Illinois at Urbana-Champai↗
DOE OSTI2024

Fluid-Dynamic Mechanisms Underlying Wind Turbine Wake Control with Strouhal-Timed Actuation

A reduction in wake effects in large wind farms through wake-aware control has considerable potential to improve farm efficiency. This work examines the success of several emerging, empirically derived control methods that modify wind turbine wakes (i.e., the pulse method, helix method, and related methods) based on Strouhal numbers on the O(0.3). Drawing on previous work in the literature for jet and bluff-body flows, the analyses leverage the normal-mode representation of wake instabilities to characterize the large-scale wake meandering observed in actuated wakes. Idealized large-eddy simulations (LES) using an actuator-line representation of the turbine blades indicate that the n=0 and ±1 modes, which correspond to the pulse and helix forcing strategies, respectively, have faster initial growth rates than higher-order modes, suggesting these lower-order modes are more appropriate for wake control. Exciting these lower-order modes with periodic pitching of the blades produces increased modal growth, higher entrainment into the wake, and faster wake recovery. Modal energy gain and the entrainment rate both increase with streamwise distance from the rotor until the intermediate wake. This suggests that the wake meandering dynamics, which share close ties with the relatively well-characterized meandering dynamics in jet and bluff-body flows, are an essential component of the success of wind turbine wake control methods. A spatial linear stability analysis is also performed on the wake flows and yields insights on the modal evolution. In the context of the normal-mode representation of wake instabilities, these findings represent the first literature examining the characteristics of the wake meandering stemming from intentional Strouhal-timed wake actuation, and they help guide the ongoing work to understand the fluid-dynamic origins of the success of the pulse, helix, and related methods.

17 WIND ENERGY↗
DOE OSTI2025

Assessing the Heat Transfer Modeling Capabilities of CFD Software for Involute-Shaped Plate Research Reactors

The ongoing efforts to convert High-Performance Research Reactors (HPRRs) using Highly Enriched Uranium (HEU) to Low-Enriched Uranium (LEU) fuel require reliable thermal–hydraulic assessments of modified core designs. The involute-shaped fuel plates used in several major HPRRs present unique modeling challenges due to their compact core geometries and high heat flux conditions. This study evaluates the capability of three commercial CFD tools, STAR-CCM+, COMSOL, and ANSYS CFX, to predict cladding-to-coolant heat transfer using Reynolds-Averaged Navier–Stokes (RANS) methods within the thermal–hydraulic regimes of involute-shaped plate reactors. Broad sensitivity analysis was conducted across a range of reactor-relevant parameters using two turbulence models (k−ϵ and k−ω SST) and different near-wall treatment strategies. The results were benchmarked against the Sieder–Tate correlation and experimental data from historic studies. The codes produced consistent results, showing good agreement with the empirical correlation of Sieder–Tate and the experimental measurements. The findings support the use of these commercial CFD codes as effective tools for assessing the thermal–hydraulic performance of involute-shaped plate HPRRs and guide future LEU core development.

CFD↗
DOE OSTI2021

Numerical study to predict the effect of surface roughness on the thermal and hydraulic performance of additively manufactured heat exchangers

Heat exchangers are used in many industries to transfer thermal energy between two or more objects with an associated thermal and hydraulic performance indicated by thermal effectiveness and pressure drop respectively. Many investigations have aimed to enhance the thermal effectiveness of heat exchangers by introducing turbulent flow inducers that increase the heat transfer rate at the expense of increasing pressure drop. Flow inducers like surface roughness can be controlled with advanced manufacturing methods like additive manufacturing for heat exchanger production. Selecting and manufacturing an optimal surface roughness for heat exchangers requires a more comprehensive approach that accounts for the surface roughness effect on the thermal and hydraulic performance of a sample. Since fabricating samples with additive manufacturing to a desired surface roughness can be an extensive experimental process, this study will fundamentally investigate the effect of surface roughness on the thermal and velocity boundary layers for various surface roughness profiles with physics-based computational modeling. To accurately model the surface roughness effect on thermohydraulic performance, correlations like roughness ratio, friction factor, and Nusselt number are investigated. This study expects to indicate correlations for different flow regimes which can successfully predict the effect of surface roughness on thermohydraulic performance. Results generated from modeling surface roughness effect with unique correlations are then compared to experimental data obtained from literature for validation. The results of this study will provide insight to understand the thermal fluid physics associated with surface roughness and design guidelines for heat exchangers produced with additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coupling of Laminar-Turbulent Transition with RANS Computational Fluid Dynamics

This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.

42 ENGINEERING↗
DOE OSTI2024

a priori uncertainty quantification of reacting turbulence closure models using Bayesian neural networks

While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulations (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such models requires reliable uncertainty estimates both in the data-informed and out-of-distribution regimes. Here, in this work, we employ Bayesian neural networks (BNNs) to capture both epistemic and aleatoric uncertainties in a reacting flow model. In particular, we model the filtered progress variable scalar dissipation rate which plays a key role in the dynamics of turbulent premixed flames. We demonstrate that BNN models can provide unique insights about the structure of uncertainty of the data-driven closure models. We also propose a method for the incorporation of out-of-distribution information in a BNN, which can be used for out-of-distribution query detection. The efficacy of the model is demonstrated by a priori evaluation on a dataset consisting of a variety of flame conditions and fuels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Filaments in and between galaxy clusters at low and mid-frequency with the SKA

Context. Understanding the magnetised Universe is a major challenge in modern astrophysics, and cosmic magnetism has been acknowledged as one of the key scientific drivers of the most ambitious radio instrument ever planned, the Square Kilometre Array. Aims. With this work, we aim to investigate the potential of the Square Kilometre Array and its precursors and pathfinders in the study of magnetic fields in galaxy clusters and filaments through diffuse synchrotron radio emission. Galaxy clusters and filaments of the cosmic web are indeed unique laboratories in which to investigate turbulent fluid motions and large-scale magnetic fields in action, and much of what is known about magnetic fields in galaxy clusters comes from sensitive radio observations. Methods. Based on cosmological magneto-hydrodynamic simulations, we predict radio properties (total intensity and polarisation) of a pair of galaxy clusters connected by a cosmic-web filament. Results. We use our theoretical expectations to explore the potential of polarimetric observations to study large-scale structure magnetic fields in the frequency ranges 50–350 MHz and 950–1760 MHz. We also present predictions for galaxy cluster polarimetric observations with the Square Kilometre Array precursors and pathfinders, such as the LOw frequency ARray 2.0 and the MeerKAT+ telescope. Conclusions. Our findings point out that polarisation observations are particularly powerful for the study of large-scale magnetic fields, since they are not significantly affected by confusion noise. The unprecedented sensitivity and spatial resolution of the intermediatefrequency radio telescopes make them the favourite instruments for the study of these sources through polarimetric data, potentially allowing us to understand if the energy density of relativistic electrons is in equipartition with the magnetic field or rather coupled with the thermal gas density. Our results show that low-frequency instruments also represent a precious tool to study diffuse synchrotron emission in total intensity and polarisation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Machine learning changes the rules for flux limiters

Learning to integrate non-linear equations from highly resolved direct numerical simulations has seen recent interest for reducing the computational load for fluid simulations. Here, we focus on determining a flux-limiter for shock capturing methods. Focusing on flux limiters provides a specific plug-and-play component for existing numerical methods. Since their introduction, an array of flux limiters has been designed. Using the coarse-grained Burgers' equation, we show that flux-limiters may be rank-ordered in terms of their log-error relative to high-resolution data. We then develop a theory to find an optimal flux-limiter and present flux-limiters that outperform others tested for integrating Burgers' equation on lattices with [Formula: see text], and 2x, 3x, 4x, and 8x coarse-grainings. We train a continuous piecewise linear limiter by minimizing the mean-squared misfit to six-grid point segments of high-resolution data, averaged over all segments. While flux limiters are generally designed to have an output of φ(r) = 1 at a flux ratio of r = 1, our limiters are not bound by this rule and yet produce a smaller error than standard limiters. Here we find that our machine learned limiters have distinctive features that may provide new rules-of-thumb for the development of improved limiters. Additionally, we use our theory to learn flux-limiters that outperform standard limiters across a range of values (as opposed to at a specific fixed value) of coarse-graining, number of discretized bins, and diffusion parameter. This demonstrates the ability to produce flux limiters that should be more broadly useful than standard limiters for general applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

GENERATING CORCOS-COHERENT SIGNAL SERIES USING PHASE PERTURBATION​

Phase perturbations are used to simulate the time-domain behavior of pressure loads from turbulent boundary layers while faithfully representing the coherence decay behaviors of turbulence-generated pressure fields. Turbulence "evolves" as it proceeds downstream, and coherence decays as the decaying exponential of the product of a flow-related constant with frequency and distance. Time domain methods for a line, a plane, a closed cyclic surface, and most of a cone were given.

Swift, Stephen Hales↗
DOE OSTI2026

Aircraft-derived particle fluxes distinguish entrainment zone and decoupled layer nucleation in marine boundary layers

The vertical distribution of freshly nucleated aerosol particles in the marine boundary layer remains poorly constrained, limiting our ability to represent new particle formation in climate models. Here we characterize 3–10 nm particle events, termed small particle events (SPEs), by deriving their vertical turbulent fluxes from aircraft measurements during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign. To overcome stationarity limitations of traditional eddy covariance methods, we applied continuous wavelet transform analysis to data collected during June–July 2017 and January–February 2018 flights over the Azores. Our flux-based analysis revealed two distinct SPE scenarios with different vertical structures and spatial extents. The first featured nucleation in the entrainment zone, where free tropospheric air entrains into the boundary layer. The second showed nucleation in the decoupled layer, a stratified region between the well-mixed surface layer and cloud-topped upper boundary layer. In both cases, convergence of air masses from different layers diluted preexisting aerosol surface area to very low levels, creating conditions favorable for nucleation and generating strong downward particle fluxes. SPEs occurred in 15 % of flights, challenging prevailing theoretical expectations that new particle formation should rarely occur in marine boundary layers due to high condensation and coagulation sink capacity of sea spray aerosols. Aircraft-derived particle fluxes provide first observational constraints on the vertical location and source strength of likely nucleation regions in the remote marine boundary layer, improving aerosol source representations in climate models and reducing uncertainties in aerosol-cloud interactions.

Rasheeda Satheesh, Ajmal [North Carolina State Uni↗
DOE OSTI2021

An on-demand, drop-on-drop method for studying enzyme catalysis by serial crystallography

Serial femtosecond crystallography has opened up many new opportunities in structural biology. In recent years, several approaches employing light-inducible systems have emerged to enable time-resolved experiments that reveal protein dynamics at high atomic and temporal resolutions. However, very few enzymes are light-dependent, whereas macromolecules requiring ligand diffusion into an active site are ubiquitous. In this work we present a drop-on-drop sample delivery system that enables the study of enzyme-catalyzed reactions in microcrystal slurries. The system delivers ligand solutions in bursts of multiple picoliter-sized drops on top of a larger crystal-containing drop inducing turbulent mixing and transports the mixture to the X-ray interaction region with temporal resolution. We demonstrate mixing using fluorescent dyes, numerical simulations and time-resolved serial femtosecond crystallography, which show rapid ligand diffusion through microdroplets. The drop-on-drop method has the potential to be widely applicable to serial crystallography studies, particularly of enzyme reactions with small molecule substrates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗
DOE OSTI2025

Nonlinear simulation of under-resolved flows with shocks

Here, we consider the numerical simulation of advection-dominated flows whose wide range of physical length scales exceed the memory capacity of finite computers. Simulating flows with shocks and turbulence presented challenges for the earliest computers that were quickly overcome by the development of new numerical methodology. Principal among those new ideas were artificial viscosity and finite volume methods, concepts that remain in common use today. We begin by describing the history of those methods, the innovators and their motivations. We then describe the development of finite scale theory, a reformulation of Navier–Stokes theory that exposes the physical principles on which artificial viscosity is based. We discuss the essential properties of the finite scale equations, the observer, unresolved kinetic energy and inviscid energy dissipation. We briefly consider the implementation of the finite scale equations on the computer from the point of view of Gisin’s conjectures about finite information.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Review of Physics-Informed Machine Learning in Fluid Mechanics

Physics-informed machine-learning (PIML) enables the integration of domain knowledge with machine learning (ML) algorithms, which results in higher data efficiency and more stable predictions. This provides opportunities for augmenting—and even replacing—high-fidelity numerical simulations of complex turbulent flows, which are often expensive due to the requirement of high temporal and spatial resolution. In this review, we (i) provide an introduction and historical perspective of ML methods, in particular neural networks (NN), (ii) examine existing PIML applications to fluid mechanics problems, especially in complex high Reynolds number flows, (iii) demonstrate the utility of PIML techniques through a case study, and (iv) discuss the challenges and opportunities of developing PIML for fluid mechanics.

42 ENGINEERING↗
DOE OSTI2022

Investigation of Gas Turbine Internal Cooling Using Supercritical CO 2 —Effect of Surface Roughness and Channel Aspect Ratio

In this paper, an experimental and numerical investigation of internal cooling channels with rib turbulators is presented with sCO 2 as the working fluid at process conditions (pressure-20.7 MPa and temperature up to 150 °C). The effect of channel aspect ratio up to 2:1 on thermal-hydraulic performance is explored in additively manufactured rectangular channels and square channels, both with and without 60 deg ribs on the top and bottom sides. The Wilson-plot method is employed to experimentally measure channel-averaged Nusselt number over a Reynolds number range up to 370,000. The friction factor is calculated from pressure drop and mass flow rate and additionally, the overall thermal performance factor (TPF) is reported. A companion computational fluid dynamics (CFD) simulation is performed for the rib turbulated cooling configurations reported in the experiments using the Reynolds average Navier–Stokes-based turbulence model. The objective of the numerical study is to gain insight into the local heat transfer augmentation in the ribbed channels as a result of varying the aspect ratio, channel configuration (square versus rectangular), operating conditions (Reynolds number) and the surface roughness, an inherent outcome of the additive manufacturing process. Surface roughness is simulated using sand grain roughness height (K S ) calculated from the experimental data, and a comparison is presented with the corresponding channel configuration with varying surface roughness heights starting from smooth surfaces (K S = 0). Experimental results indicate that the heat transfer augmentation is negligible in the rectangular channels with ribs on the long side compared to the square channel. However, it is enhanced by 60% in comparison to placing ribs on the shorter side. The TPF remains constant at around 1 for the entire range of Reynolds numbers consistent with prior work at the National Energy Technology Laboratory (NETL). Here, the simulation results highlight that increased surface roughness can have a favorable considerable influence on Nusselt number and overall thermal performance enhancement.

42 ENGINEERING↗
DOE OSTI2021

Evaluation of residual gas fraction estimation methods for cycle-to-cycle combustion variability analysis and modeling

Cycle-to-cycle combustion variability in spark-ignition engines during normal operation is mainly caused by random perturbations of the in-cylinder conditions such as the flow velocity field, homogeneity of the air-fuel distribution, spark energy discharge, and turbulence intensity of the flame front. Such perturbations translate into the variability of the energy released observed at the end of the combustion process. During normal operating conditions, the cycle-to-cycle variability (CCV) of the energy release behaves as random uncorrelated noise. However, during diluted combustion, in either the form of exhaust gas recirculation (EGR) or excess air (lean operation), the CCV tends to increase as dilution increases. Moreover, when the ignition limit is reached at high dilution levels, the combustion CCV is exacerbated by sporadic occurrences of incomplete combustion events, and the uncorrelation assumption no longer holds. The low or null energy released by partial burns and misfires has an impact on the following combustion event due to the residual gas that carries burned and unburned gases, which contributes to the deterministic coupling between engine cycles. Many residual gas fraction estimation methods, however, only address the nominal case where complete combustion occurs and combustion events are uncorrelated. Here we evaluate the efficacy of such methods on capturing the effects of partial burns and misfires on the residual gas estimate for high-EGR operation. The advantages and disadvantages of each method are discussed based on their ability to generate cycle-to-cycle estimates. Finally, a comparison between the different estimation techniques is presented based on their usefulness for control-oriented modeling.

42 ENGINEERING↗
DOE OSTI2024

Challenges in Rotor Aerodynamic Modeling for Non-Uniform Inflow Conditions

Within an international collaboration framework, the accuracy of rotor aerodynamic models used for design load calculations of wind turbines is being assessed. Where the use of high-fidelity computation fluid dynamics (CFD) and mid-fidelity free-vortex wake (FVW) models has become commonplace within the wind energy community, these still fail to meet the requirements in terms of execution time and computational cost needed for design load calculations. The fast but engineering fidelity blade-element/momentum (BEM) method can therefore still be considered the industry workhorse for design load simulations. At the same time, upscaling of wind turbine rotors makes inflow non-uniformities (e.g. shear, veer, turbulence) more important. The objective of this work is to assess model accuracy in non-uniform inflow conditions, which violate several BEM assumptions. Thereto a comparison in turbulent inflow has been executed including a wide variety of codes, focusing on the DanAero field measurements, where a 2.3-MW turbine was equipped with, among other sensors, a pressure measurement apparatus. The results indicate that, although average load patterns are in good agreement, this does not hold for the unsteady loads that drive fatigue damage and aero-elastic stability. A simplified comparison round in vertical shear was initiated to investigate the observed differences in a more controlled manner. A consistent offset in load amplitude was observed between CFD and free-vortex codes on the one hand and BEM-type codes on the other hand. To shed more light on the observations, dedicated efforts are ongoing to pinpoint the cause for these differences, in the end leading to guidelines for an improved BEM implementation.

17 WIND ENERGY↗