Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOUNDARY LAYER, LAMINAR - WALLS”

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.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Aeroacoustics Noise Model of OpenFAST

The report describes theory and application of a newly released module of OpenFAST to simulate the aeroacoustic noise generated by the rotor of an arbitrary wind turbine. OpenFAST is a fully open-source publicly available wind turbine analysis tool actively developed at the National Renewable Energy Laboratory. The aeroacoustic module, which is also fully open-source and publicly available, is based on work performed over the past three decades. Frequency-based models for turbulent inflow, turbulent boundary layer – trailing edge, laminar boundary layer – vortex shedding, tip vortex, and trailing edge bluntness – vortex shedding noise mechanisms are included. A simple directivity model is also included. The appendices of this report describe in detail the inputs to the module and how to find them in the input files of OpenFAST. The noise models are exercised simulating the aeroacoustic noise emissions of the IEA Wind Task 37 land-based reference wind turbine. A code-to-code comparison between the implementation presented here and the implementation available at the Wind Energy Institute of the Technical University of Munich, Germany, is also presented.

17 WIND ENERGY↗
DOE OSTI2024

An extension of Thwaites’ method for turbulent boundary layers

Thwaites developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations. The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

17 WIND ENERGY↗
DOE OSTI2022

Critical Shear Stress for Erosion Under Laminar Jet Flow

Particle erosion is an important process in both natural and manmade environments. For example, in natural environments particle or sediment erosion rates determine the course of rivers, depth of streams, and fate of deltas. The interplay between erosion rates and flow patterns can have dramatic impact on shorelines. Where hydrodynamic jet flows or rip currents change direction, sediments may accumulate or disperse altering the shape of coastlines and the ecosystems and economies that rely on them. In man-made systems, jet flows may be used to mobilize particles from surfaces or within industrial scale mixing tanks. For example, in the nuclear industry, impinging jet mixers may be used to mobilize and suspend sludges at the bottom of waste tanks. In each of these examples, the Shields diagram is useful to determine the conditions under which erosion occurs. This diagram plots the scaled critical shear stress for erosion as a function of the particle size as embedded in particle Reynolds numbers, Archimedes numbers, or their equivalents. Although the data, that carries a large spread, is largely in agreement with mathematical models, for the smallest particles, model predictions diverge from each other and from experimental data. Indeed, there is substantial uncertainty about the Shields diagrams at relatively small particle sizes and flow rates. Here we evaluate critical shear stresses for erosion under laminar jet flow. To date critical shear stresses for erosion, as embodied in Shields diagrams, have largely been developed for turbulent flow both experimental and mathematically. However, curves on the Shields diagrams diverge among various models at lower particle Reynolds numbers. Indeed, how critical shear stresses for erosion develop under fully laminar jet flow conditions (not simply laminar boundary layers under turbulent flow) remains unclear. Here we address this gap, which is particularly important for mobilization of the smallest particles. We find the Shields parameter to be inversely proportional to the particle Reynolds numbers under laminar jet flow conditions similar to turbulent conditions.

critical shear stresses, erosion rates, effective ↗
DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗
DOE OSTI2020

HTR solver: An open-source exascale-oriented task-based multi-GPU high-order code for hypersonic aerothermodynamics

In this study, the open-source Hypersonics Task-based Research (HTR) solver for hypersonic aerothermodynamics is described. The physical formulation of the code includes thermochemical effects induced by high temperatures (vibrational excitation and chemical dissociation). The HTR solver uses high-order TENO-based spatial discretization on structured grids and efficient time integrators for stiff systems, is highly scalable in GPU-based supercomputers as a result of its implementation in the Regent/Legion stack, and is designed for direct numerical simulations of canonical hypersonic flows at high Reynolds numbers. Additionally, the performance of the HTR solver is tested with benchmark cases including inviscid vortex advection, low- and high-speed laminar boundary layers, inviscid one-dimensional compressible flows in shock tubes, supersonic turbulent channel flows, and hypersonic transitional boundary layers of both calorically perfect gases and dissociating air.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Pressure-Gradient-Based RANS Model for Predicting Separation in Transitional and Turbulent Flows

Predicting flow separation poses a significant challenge for RANS models, particularly in transitional flows over airfoils. We propose a novel improvement to RANS models to predict incipient separation in both transitional and fully turbulent flows. Our approach modifies the eddy viscosity model in regions indicated by a pressure-gradient criterion that accounts for intermittency - determining whether the boundary layer is laminar or turbulent. This model demonstrates robust generalization across various airfoil shapes and Reynolds numbers. Applied to the NREL Phase VI wind turbine rotor, our model shows improved aerodynamic performance predictions compared to the baseline RANS model.

k-omega SST↗
DOE OSTI2025

Assessment and Improvement of the SST-Gamma Transition Model in Nalu-Wind

We conduct laminar–turbulent boundary-layer transition simulations using a local correlation-based transition model for two-dimensional incompressible flow and present enhancements to improve the accuracy of transition predictions. Menter’s Galilean-invariant 𝛾 transition model is implemented in the incompressible, unstructured-grid flow solver Nalu-Wind and is validated against experimental data and results from NASA’s flow solvers. The test cases of the AIAA Transition Prediction and Modeling Workshop are investigated, namely, the T3A/T3B flat plates and the NLF(1)-0416 and S809 airfoils. Based on the results, best practices for transition simulations, particularly for an unstructured-grid flow solver, are identified. Additional airfoil simulations are conducted for two wind turbine airfoils, S822 at Reynolds numbers of 𝑂⁡(10 5 ) and DU00-W-212 at Reynolds numbers of 𝑂⁡(10 7 ), to assess the model at low and high Reynolds numbers. Furthermore, through this work, we propose several approaches to enhance transition simulations, including 1) enforcing positivity of the implicit operator for the source terms of the transition model, 2) employing a constant turbulence intensity in stationary external flow simulations, and 3) recommending meshing for unstructured-grid flow solvers. Finally, we provide detailed documentation of the validation and data for the canonical cases to the transition modeling community.

17 WIND ENERGY↗
DOE OSTI2022

Multi-fidelity information fusion with concatenated neural networks

Recently, computational modeling has shifted towards the use of statistical inference, deep learning, and other data-driven modeling frameworks. Although this shift in modeling holds promise in many applications like design optimization and real-time control by lowering the computational burden, training deep learning models needs a huge amount of data. This big data is not always available for scientific problems and leads to poorly generalizable data-driven models. This gap can be furnished by leveraging information from physics-based models. Exploiting prior knowledge about the problem at hand, this study puts forth a physics-guided machine learning (PGML) approach to build more tailored, effective, and efficient surrogate models. For our analysis, without losing its generalizability and modularity, we focus on the development of predictive models for laminar and turbulent boundary layer flows. In particular, we combine the self-similarity solution and power-law velocity profile (low-fidelity models) with the noisy data obtained either from experiments or computational fluid dynamics simulations (high-fidelity models) through a concatenated neural network. We illustrate how the knowledge from these simplified models results in reducing uncertainties associated with deep learning models applied to boundary layer flow prediction problems. The proposed multi-fidelity information fusion framework produces physically consistent models that attempt to achieve better generalization than data-driven models obtained purely based on data. While we demonstrate our framework for a problem relevant to fluid mechanics, its workflow and principles can be adopted for many scientific problems where empirical, analytical, or simplified models are prevalent. In line with grand demands in novel PGML principles, this work builds a bridge between extensive physics-based theories and data-driven modeling paradigms and paves the way for using hybrid physics and machine learning modeling approaches for next-generation digital twin technologies.

42 ENGINEERING↗
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 OSTI2022

Medium-Scale Methanol Pool Fire Model Validation

In this work, medium scale (30 cm diameter) methanol pool fires were simulated using the latest fire modeling suite implemented in Sierra/Fuego, a low Mach number multiphysics reacting flow code. The sensitivity of model outputs to various model parameters was studied with the objective of providing model validation. This work also assesses model performance relative to other recently published large eddy simulations (LES) of the same validation case. Two pool surface boundary conditions were simulated. The first was a prescribed fuel mass flux and the second used an algorithm to predict mass flux based on a mass and energy balance at the fuel surface. Gray gas radiation model parameters (absorption coefficients and gas radiation sources) were varied to assess radiant heat losses to the surroundings and pool surface. The radiation model was calibrated by comparing the simulated radiant fraction of the plume to experimental data. The effects of mesh resolution were also quantified starting with a grid resolution representative of engineering type fire calculations and then uniformly refining that mesh in the plume region. Simulation data were compared to experimental data collected at the University of Waterloo and the National Institute of Standards and Technology (NIST). Validation data included plume temperature, radial and axial velocities, velocity temperature turbulent correlations, velocity velocity turbulent correlations, radiant and convective heat fluxes to the pool surface, and plume radiant fraction. Additional analyses were performed in the pool boundary layer to assess simulated flame anchoring and the effect on convective heat fluxes. This work assesses the capability of the latest Fuego physics and chemistry model suite and provides additional insight into pool fire modeling for nonluminous, non-sooting flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Evaluating Scaling Effects of Steam Condensation Physics-based Model for Small Modular Reactor

This study focused on assessing condensation heat transfer (CHT) models for enhancing the passive safety of small modular reactors (SMRs). It also encompassed a review of the theoretical and semi-theoretical CHT models and correlations. Condensation models such as the degradation factor, heat and mass transfer analogy (HMTA), and boundary layer models are used to study advanced nuclear reactor PCCS. The previous CHT models were developed for small tube geometries, which differ from the containment of SMR. This study examined the scaled relations of a physics-based CHT model for analyzing laminar and turbulent film condensation downward flow conditions inside scaled tube geometries.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Presentation: Evaluating scaling effects of steam condensation physics-based model for small modular reactors safety

This study focused on assessing condensation heat transfer (CHT) models for enhancing the passive safety of small modular reactors (SMRs). It also encompassed a review of the theoretical and semi-theoretical CHT models and correlations. Condensation models such as the degradation factor, heat and mass transfer analogy (HMTA), and boundary layer models are used to study advanced nuclear reactor PCCS. The previous CHT models were developed for small tube geometries, which differ from the containment of SMR. This study examined the scaled relations of a physics-based CHT model for analyzing laminar and turbulent film condensation downward flow conditions inside scaled tube geometries.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

On the interaction of a wind turbine wake with a conventionally neutral atmospheric boundary layer

In this work, we investigate the dynamics of wind turbine tip-vortex breakdown in a conventionally neutral atmospheric boundary layer (ABL). To this end, high-resolution data are collected from large-eddy simulations of a wind turbine operating within a neutral ABL and studied by means of proper orthogonal decomposition (POD) and Fourier analysis. The high resolution of the generated data in both space and time allows us to gain insight into the tip-vortex breakdown mechanisms by (i) capturing the energy modes of the coherent structures, (ii) studying their contribution to the tip-vortex breakdown through their power spectra functions and mean kinetic energy (MKE) flux, and (iii) analysing the growth rate of each contributing perturbation frequency along tip vortices. Our analysis shows that under a fully turbulent scenario, the growth rate of perturbations along the tip vortices is largest for low wave numbers, i.e. long-wave perturbations. Additionally, the MKE flux reaches its highest value at two diameters downstream of the rotor plane, a behaviour that can be attributed to the coexistence of multiple interacting POD modes, with the streamwise vortex roller mode being the primary contributor to the total MKE flux budget, contributing approximately . Finally, comparisons with a laminar, uniform flow scenario subject to a single-frequency perturbation highlight the differences between the two ambient flow conditions. In the non-turbulent, uniform flow scenario, the growth rate attains its maximum value at a wave number corresponding to the out-of-phase mutual-inductance mechanism, whereas the MKE flux exhibits local minima and maxima along the wake and at different downstream locations depending on the perturbation frequency. Our analyses suggest that the breakdown of the wind turbine tip vortices under a fully turbulent neutral ABL inflow is due to complex interactions across a range of excitation frequencies, in which the mutual-inductance instability may not be the dominant one.

17 WIND ENERGY↗
DOE OSTI2023

Expanded verification and validation studies of hypersonic aerodynamics with multiple physics-fidelity models

Hypersonic aerothermodynamics is an important domain of modern multiphysics simulation. The Multi-Fidelity Toolkit is a simulation tool being developed at Sandia National Laboratories to predict aerodynamic properties for compressible flows from a range of physics fidelities and computational speeds. These models include the Reynolds-averaged Navier–Stokes (RANS) equations, the Euler equations with momentum-energy integral technique (MEIT), and modified Newtonian aerodynamics with flat-plate boundary layer (MNA+FPBL) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities, as well as model validation with uncertainty considerations. Code verification activities on the MNA+FPBL model build on previous work by focusing on the viscous portion of the model. Viscous quantities of interest are compared against those from an analytical solution for flat-plate, inclined-plate, and cone geometries. The code verification methodology for the MEIT model is also presented. Test setup and results of code verification tests on the laminar and turbulent models within MEIT are shown. Solution-verification activities include grid-refinement studies on simulations that model the HIFiRE-1 wind tunnel experiments. These experiments are used for validation of all model fidelities. A thorough validation comparison with prediction error and uncertainty is also presented. Three additional HIFiRE-1 experimental runs are simulated in this study, and the solution verification and validation work examines the effects of the associated parameter changes on model performance. Finally, a study is presented that compares the computational costs and fidelities from each of the different models.

42 ENGINEERING↗
DOE OSTI2025

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗
DOE OSTI2022

Onset of instabilities in rotating flows by direct numerical simulation

A rotating disk is the canonical experiment for measuring surface reaction rates in geochemical and electrochemical systems. Using the similarity solution for laminar flow around an infinite disk, the mass transfer coefficient can be simply related to the intrinsic reaction rate at the surface. However, measurements of mass transfer rates use a finite-size disk within a larger container of solution; here the flow is no longer strictly laminar, but there must always be some recirculation. Our interest was initially in the assumption of a uniform radial concentration field, how this breaks down near the perimeter of the disk, and what effect that might have on the measured mass transfer rates. However, our numerical simulations suggest that the flow around a finite-size disk becomes time dependent at Reynolds number ($\textit{Re}$) below 1000, which is much smaller than the typical values in mass-transfer measurements ($\textit{Re}$~10 4 ). In this work, we observe the formation of coherent structures in the flow, which suggest a non-uniform mass transfer at the disk surface. The rotating-disk flow follows a similar sequence of instabilities to the Taylor–Couette flow: a centrifugal instability leading an axisymmetric, time-invariant flow, followed by a Hopf bifurcation to a time-periodic flow. To minimise the possibility that our results are a numerical artefact, we have also simulated the instability in the stationary boundary layer of a rotor–stator flow, comparing with self-similar solutions at low $\textit{Re}$ and with spectral methods near the critical Reynolds number.

42 ENGINEERING↗