Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boundary layer thickness”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Editors’ Choice—Natural Convection Boundary Layer Thickness at Elevated Chloride Concentrations and Temperatures and the Effects on a Galvanic Couple

The natural convection boundary layer ( δ n c ) and its influence on cathodic current in a galvanic couple under varying electrolytes as a function of concentration (1 − 5.3 M NaCl) and temperature (25 °C−45 °C) were understood. Polarization scans were obtained under quiescent conditions and at defined boundary layer thicknesses using a rotating disk electrode on platinum and stainless steel 304L (SS304L); these were combined to determine δ n c . With increasing chloride concentration and temperature, δ n c decreased. Increased mass transport (Sherwood number) results in a decrease in δ n c , providing a means to predict this important boundary. Using Finite Element Modeling, the cathodic current was calculated for an aluminum alloy/SS304L galvanic couple as a function of water layer ( WL ) thickness and cathode length. Electrolyte domains were delineated, describing (i) dominance of ohmic resistance over mass transport under thin WL , (ii) the transition from thin film to bulk conditions at δ n c , and (iii) dominance of mass transport under thick WL . With increasing chloride concentration, cathodic current decreased due to decreases in mass transport. With increasing temperature, increased cathodic current was related to increases in mass transport and solution conductivity. This study has implications for sample sizing and corrosion prediction under changing environments.

25 ENERGY STORAGE↗
DOE OSTI2022

Aero-thermal numerical characterization of blunt fin-induced shock wave–boundary layer interaction and its control through leading-edge cooling injection

This study represents a novel evaluation of active flow control to alleviate the aerothermal penalties created by the blunt fin-induced shock wave–boundary layer interaction. The manuscript analyzes the effect of flow injection on a blunt fin-induced shock wave–boundary layer interaction via computational fluid dynamics simulations with various degrees of resolution. The impact on the mean flow topology and wall variables was investigated utilizing Reynolds-averaged Navier–Stokes simulations. Detached-eddy simulations revealed the low-frequency shock motion, shock wave–boundary layer, and horseshoe vortex interaction. The test article was exposed to two different incoming boundary layer thicknesses; the thicker boundary layer led to the appearance of larger turbulent scales. The Detached-eddy simulations revealed the time history of the shock wave–boundary layer interaction, focusing on the inception and development of the recirculated flow regions. Ultimately, spectral proper orthogonal decomposition was employed to identify the structures associated with the low-frequency shock motion caused by the shock wave–boundary layer interaction.

Lozano, Francisco (ORCID:0000000277871125)↗
DOE OSTI2021

Temporally and spatially resolved X-ray densitometry in a shock tube

Gas densities were measured by x-ray absorption spectroscopy to examine spatial inhomogeneity in density behind incident and reflected shock waves in a miniature (12.7 mm bore) high repetition rate shock tube. Shock waves were generated in argon. The pressure and temperature behind the incident shock were P 2 ~ 1.72 bar and T 2 ~ 800 K, while those behind the reflected shock were P 5 ~ 6.5 bar and T 5 ~ 1480 K. Time-resolved line-of-sight transmission measurements using x-rays at 9 keV photon energy were made along various axial and transverse locations covering the entire shock tube cross section. The x-ray transmission at each location was converted to pathlength-integrated densities, which in some conditions could be further converted to pathlength-averaged densities. The measured gas densities were compared with those calculated by the normal shock wave relations, and good agreement was found in the preshock and immediately behind the incident shock regions. Similar agreement was found for reflected shock conditions very near the endwall. However, increasingly large differences between the measured and calculated gas densities were found as distance from the endwall increased. Furthermore, reconstruction by Abel inversion allowed time and radially resolved densities to be obtained from the transverse measurements. While the profiles measured in this experiment are of insufficient signal-to noise level and resolution to definitively assign boundary layer thicknesses, future improvements may be able to do this. The radially resolved densities reveal growing sidewall thermal boundary layers in gases behind the reflected shock. These reconstructions were compared with a model of the thermal boundary layer and are consistent with values of the thermal boundary layer thickness being between 0.25 mm and 1 mm at 0.7 ms after shock reflection.

42 ENGINEERING↗
DOE OSTI2022

Direct numerical simulation of hypersonic turbulent boundary layers: effect of spatial evolution and Reynolds number

Direct numerical simulations (DNS) are performed to investigate the spatial evolution of flat-plate zero-pressure-gradient turbulent boundary layers over long streamwise domains ( >300δ i , with δ i the inflow boundary-layer thickness) at three different Mach numbers, 2.5, 4.9 and 10.9, with the surface temperatures ranging from quasiadiabatic to highly cooled conditions. The settlement of turbulence statistics into a fully developed equilibrium state of the turbulent boundary layer has been carefully monitored, either based on the satisfaction of the von Kármán integral equation or by comparing runs with different inflow turbulence generation techniques. The generated DNS database is used to characterize the streamwise evolution of multiple important variables in the high-Mach-number, cold-wall regime, including the skin friction, the Reynolds analogy factor, the shape factor, the Reynolds stresses, and the fluctuating wall quantities. The data confirm the validity of many classic and newer compressibility transformations at moderately high Reynolds numbers (up to friction Reynolds number Re τ ≈1200 ) and show that, with proper scaling, the sizes of the near-wall streaks and superstructures are insensitive to the Mach number and wall cooling conditions. The strong wall cooling in the hypersonic cold-wall case is found to cause a significant increase in the size of the near-wall turbulence eddies (relative to the boundary-layer thickness), which leads to a reduced-scale separation between the large and small turbulence scales, and in turn to a lack of an outer peak in the spanwise spectra of the streamwise velocity in the logarithmic region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

On the enhancement of boundary layer skin friction by turbulence: an angular momentum approach

Turbulence enhances the wall shear stress in boundary layers, significantly increasing the drag on streamlined bodies. Other flow features such as free stream pressure gradients and streamwise boundary layer growth also strongly influence the local skin friction. In this paper, an angular momentum integral (AMI) equation is introduced to quantify these effects by representing them as torques that alter the shape of the mean velocity profile. This approach uniquely isolates the skin friction of a Blasius boundary layer in a single term that depends only on the Reynolds number most relevant to the flow's engineering context, so that other torques are interpreted as augmentations relative to the laminar case having the same Reynolds number. The AMI equation for external flows shares this key property with the so-called FIK relation for internal flows (Fukagata et al. , Phys. Fluids , vol. 14, 2002, pp. L73–L76). Without a geometrically imposed boundary layer thickness, the length scale in the Reynolds number for the AMI equation may be chosen freely. After a brief demonstration using Falkner–Skan boundary layers, the AMI equation is applied as a diagnostic tool on four transitional and turbulent boundary layer direct numerical simulation datasets. Regions of negative wall-normal velocity are shown to play a key role in limiting the peak skin friction during the late stages of transition, and the relative strengths of terms in the AMI equation become independent of the transition mechanism a very short distance into the fully turbulent regime. The AMI equation establishes an intuitive, extensible framework for interpreting the impact of turbulence and flow control strategies on boundary layer skin friction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Design and Characterization of Highly Diffusive Turbine Vanes Suitable for Transonic Rotating Detonation Combustors

In rotating detonation engines the turbine inlet conditions may be transonic with unprecedented unsteady fluctuations. To ensure an acceptable engine performance, the turbine passages must be suited to these conditions. This article focuses on designing and characterizing highly diffusive turbine vanes to operate at any inlet Mach number up to Mach 1. First, the effect of pressure loss on the starting limit is presented. Afterward, a multi-objective optimization with steady RANS simulations, including the endwall and 3D vane design is performed. Compared to previous research, significant reductions in pressure loss and stator-induced rotor forcing are obtained, with an extended operating range and preserving high flow turning. Finally, the influence of the inlet boundary layer thickness on the vane performance is evaluated, inducing remarkable increases in pressure loss and downstream pressure distortion. Employing an optimization with a thicker inlet boundary layer, specific endwall design recommendations are found, providing a notable improvement in both objective functions.

Grasa, Sergio↗
DOE OSTI2022

On the Response of a Herschel–Bulkley Fluid Due to a Moving Plate

In this paper, we study the boundary-layer flow of a Herschel–Bulkley fluid due to a moving plate; this problem has been experimentally investigated by others, where the fluid was assumed to be Carbopol, which has similar properties to cement. The computational fluid dynamics finite volume method from the open-source toolbox/library OpenFOAM is used on structured quad grids to solve the mass and the linear momentum conservation equations using the solver “overInterDyMFoam” customized with non-Newtonian viscosity libraries. The governing equations are solved numerically by using regularization methods in the context of the overset meshing technique. The results indicate that there is a good comparison between the experimental data and the simulations. The boundary layer thicknesses are predicted within the uncertainties of the measurements. The simulations indicate strong sensitivities to the rheological properties of the fluid.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Alcove formation in dissolving cliffs driven by density inversion instability

We demonstrate conditions that give rise to cave-like features commonly found in dissolving cliffsides with a minimal two-phase physical model. Alcoves that are wider at the top and tapered at the bottom, with sharp-edged ceilings and sloping floors, are shown to develop on vertical solid surfaces dissolving in aqueous solvents. As evident from descending plumes, sufficiently large indentations evolve into alcoves as a result of the faster dissolution of the ceiling due to a solutal Rayleigh–Bénard density inversion instability. In contrast, defects of size below the boundary layer thickness set by the critical Rayleigh number smooth out, leading to stable planar interfaces. Furthermore, the ceiling recession rate and the alcove opening area evolution are shown to be given to first-order by the critical Rayleigh number. By tracking passive tracers in the fluid phase, we show that the alcoves are shaped by the detachment of the boundary layer flow and the appearance of a pinned vortex at the leading edge of the indentations. The attached boundary layer past the developing alcove is then found to lead to rounding of the other sides and the gradual sloping of the floor.

58 GEOSCIENCES↗
DOE OSTI2026

Mechanistic mass transfer in hollow fiber membrane solvent extraction for bio-based isobutanol

Membrane solvent extraction (MSE) has emerged as a promising method for selectively recovering bioproducts from complex aqueous streams. Bio-isobutanol, a next-generation feedstock for biofuel, remains challenging to recover because of its low concentration and the presence of inhibitory substances. This study explores the potential of hollow fiber (HF) MSE for bio-isobutanol recovery and systematically examines the coupled effects of fiber packing, shell-side flow dynamics, and aqueous chemistry on performance. A resistance-in-series model is applied to understand mass transfer in the HF MSE modules, quantify local resistances, and validate overall performance. The results show that increasing the fiber packing provides a larger interfacial area but induces poor flow distribution and channeling, hindering effective isobutanol transport. Meanwhile, increasing the shell-side velocity improves isobutanol recovery due to reductions in the boundary layer thickness. The presence of salts, added to mimic fermentation broth, increases the partition coefficient through salting-out effects, further improving isobutanol flux. A modified correlation for the shell-side mass transfer coefficient (k s,ϕ+v ), integrating geometric and hydrodynamic effects, was developed and validated. The proposed model achieves highly predictive accuracy (r 2 = 0.9808) across a wide range of conditions, outperforming previous models. The findings provide mechanistic insight into the interaction of geometric packing, hydrodynamics, and chemistry in governing mass transfer in HF MSE. Overall, this work demonstrates the potential of HF MSE for efficient bio-isobutanol recovery and also provides practical guidelines on critical factors (packing fraction, partition coefficient, and shell-side velocity), aiding in the design and scaling of MSE systems for resource recovery.

Aqueous chemistry↗
DOE OSTI2022

Heat transport and temperature boundary-layer profiles in closed turbulent Rayleigh–Bénard convection with slippery conducting surfaces

Here, we report direct numerical simulations (DNS) of the Nusselt number Nu , the vertical profiles of mean temperature Θ(z) and temperature variance Ω(z) across the thermal boundary layer (BL) in closed turbulent Rayleigh–Bénard convection (RBC) with slippery conducting surfaces ( z is the vertical distance from the bottom surface). The DNS study was conducted in three RBC samples: a three-dimensional cuboid with length L=H and width W=H/4 ( H is the sample height), and two-dimensional rectangles with aspect ratios Γ≡L/H=1 and 10. The slip length b for top and bottom plates varied from 0 to ∞ . The Rayleigh numbers Ra were in the range 10 6 $⩽$Ra$⩽$10 10 and the Prandtl number Pr was fixed at 4.3 . As b increases, the normalised Nu/Nu 0 ( Nu 0 is the global heat transport for b =0 ) from the three samples for different Ra and Γ can be well described by the same function Nu/Nu 0 = N 0 tanh( b /λ 0 )+1 , with N 0 =0.8±0.03 . Here λ 0 ≡ L /(2 Nu 0 ) is the thermal boundary layer thickness for b =0 . Considering the BL fluctuations for Pr >1, one can derive solutions of temperature profiles Θ(z) and Ω(z) near the thermal BL for b $⩾$ 0 . When b =0 , the solutions are equivalent to those reported by Shishkina et al. ( Phys. Rev. Lett. , vol. 114, 2015, 114302) and Wang et al. ( Phys. Rev. Fluids , vol. 1, 2016, 082301(R)), respectively, for no-slip plates. For b >0, the derived solutions are in excellent agreement with our DNS data for slippery plates.

42 ENGINEERING↗
DOE OSTI2025

Tracking Local pH Dynamics during Water Electrolysis via In-Line Continuous Flow Raman Spectroscopy

The performance of electrochemical devices, which play a critical role in decarbonization efforts, is often governed by proton-coupled electron transfer reactions at the electrode–electrolyte interface. These reactions are highly sensitive to the complex and dynamic microenvironment present at the electrode surface. However, characterizing this environment─particularly monitoring interfacial pH and its evolution under reaction conditions─remains challenging, necessitating the development of advanced analytical tools. Here, in this study, we introduce in-line continuous flow Raman spectroscopy (CFRS) as a spectroelectrochemical platform for quantifying interfacial pH swings generated during water-splitting. By monitoring phosphate ion speciation and controlling the hydrodynamics with a flow cell, we measure pH swings as a function of current density, flow rate, and distance from the electrode. Comparison with theoretical models reveals the impact of bulk pH, boundary layer thickness, and bubble dynamics at high current densities. Collectively, these findings establish CFRS as a platform for quantitatively investigating pH dynamics, offering critical insights for advancing electrochemical energy conversion technologies.

Marquez, Raul A. [Univ. of Texas, Austin, TX (Unit↗
DOE OSTI2026

Spatially Mapping the CO 2 Alkaline Sorbent Diffuse Microenvironment Using Operando Raman Spectroscopy

When designing a chemical process, the local balance of transport and kinetics, collectively referred to as the diffuse microenvironment, plays a critical role in performance but is difficult to directly observe. This work demonstrates a method of two-dimensional spatial chemical mapping of the diffuse microenvironment in the context of alkali metal hydroxide direct air capture of carbon dioxide using a custom operando gas-absorption flow cell along with confocal Raman spectroscopy. Notably, we observe the concentration boundary layer near the gas–liquid interface and elucidate the interplay of carbonate and bicarbonate ions within it while inferring local hydroxide depletion through continuum modeling. These first of their kind observations provide a technique to compare the performance of direct air capture solvents based on diffuse microenvironment dynamics while also providing metrics important for air contactor design such as boundary layer thickness. Overall, this work showcases a new experimental platform to study interfacial diffuse microenvironments in and outside of the field of direct air capture of carbon dioxide.

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

The Cascade Effectiveness of 3-Terminal Tandem Photocathode Architectures as Applied to CO 2 Reduction

Cascade catalysis for photoelectrochemical CO 2 reduction (CO 2 R) decouples the overall reaction into sequential steps occurring on separately optimized catalysts (for example, Ag and Cu) between which an intermediate species such as CO is transferred. A 3-terminal tandem (3TT) photovoltaic architecture advantageously holds two different catalytic regions at different potentials under a single illumination source, but its overall efficiency is low. Using a stochastic reaction-diffusion model, we have examined 3TT photocathode design principles, focusing on the coupling of surface chemistry and transport of CO both inside the boundary layer and outward, into the bulk electrolyte. We find that ensuring that the lateral diffusion distance within the boundary layer is short compared to the boundary layer thickness and controlling bulk flow are key, with interdigitated designs showing an overall conversion efficiency improvement by 2 orders of magnitude compared to the side-by-side case without flow. These findings can be adapted for other cascaded architectures.

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

Low-temperature Leidenfrost-like jumping of sessile droplets on microstructured surfaces

The Leidenfrost effect—the levitation and hovering of liquid droplets on hot solid surfaces—generally requires a sufficiently high substrate temperature to activate liquid vaporization. Here, in this paper, we report the modulation of Leidenfrost-like jumping of sessile water microdroplets on micropillared surfaces at a relatively low temperature. Compared to traditional Leidenfrost effect occurring above 230 °C, the fin-array-like micropillars enable water microdroplets to levitate and jump off the surface within milliseconds at a temperature of 130 °C by triggering the inertia-controlled growth of individual vapour bubbles at the droplet base. We demonstrate that droplet jumping, resulting from momentum interactions between the expanding vapour bubble and the droplet, can be modulated by tailoring of the thermal boundary layer thickness through pillar height. This enables regulation of the bubble expansion between the inertia-controlled mode and the heat-transfer-limited mode. The two bubble-growth modes give rise to distinct droplet jumping behaviours characterized by constant velocity and constant energy regimes, respectively. This heating strategy allows the straightforward purging of wetting liquid droplets on rough or structured surfaces in a controlled manner, with potential applications including the rapid removal of fouling media, even when located in surface cavities.

42 ENGINEERING↗
DOE OSTI2024

Validation of predictive performance models for supersonic gas-jet nozzles at the Laboratory for Laser Energetics

We present results characterizing the neutral-density distributions produced by the supersonic nozzles used in experiments on the OMEGA-60 and OMEGA-EP laser systems at the University of Rochester’s Laboratory for Laser Energetics (LLE). Axisymmetric Fluent® simulations using LLE nozzle specifications capture the viscous effects, gas expansion, and shock waves that complicate flow predictions for offsets above the nozzle exit. Here, these simulations show good agreement with neutral-density measurements obtained using a four-wave shearing interferometer. An analytical form is given for the plateau length. Fits to simulation data for boundary layer thickness, mean plateau density, and density ramps are given as functions of nozzle offset and nozzle backing pressure for a number of nozzles and gases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Computational Analysis of Water Film Thickness During Rain Events for Assessing Hydroplaning Risk Part 2: Rough Road Surfaces

Hydroplaning occurs when a water film forms on a road during a rain event and vehicles are traveling at a speed that does not provide sufficient time for the tires to push the water film out of the tire path. Under these conditions, the tire loses contact with the road and the driver may lose control over the vehicle. The present study is a part of a multiyear research effort aimed at gaining a better understanding of the factors that contribute to hydroplaning risk in order to minimize the occurrence of hydroplaning accidents. Three-dimensional computational fluid dynamics (CFD) analysis was used to investigate the water film thickness on multilane roads, one of the key parameters in evaluating hydroplaning risk. In the first part of the study, water films forming on wide roadways, i.e. multilane highways with 2, 4, and 6 lanes per side, were analyzed, with varying cross slope, longitudinal slope, and rainfall rate. Roadways with and without curbs and drainage were included in the analysis. The analysis in the first part was limited to nearly smooth asphalt or concrete surfaces because the maximum roughness height that can be specified in a CFD model using wall functions to determine shear stress at the road surface boundary is one half the thickness of computational cell layer adjacent to the boundary and the need to resolve the vertical velocity distribution in a thin water film requires cell layers that are a fraction of a millimeter thick. In this part of the study, the analysis was extended to consider the influence of pavement roughness over a range extending up 3.6 mm on the development of water film thickness. Various methods of modeling roughness were investigated and are presented: the roughness height model, a meshed-out geometry of the macrotexture of the pavement, and a porous region model. The computational results were analyzed and compared with experimental measurements performed by Gallaway et al.

42 ENGINEERING↗
DOE OSTI2026

Constraining the Dynamo Layers in Jupiter and Saturn with Observations and Scaling Laws

The dipole-dominated magnetic fields of Jupiter and Saturn provide evidence for active dynamos operating within their deep interiors, yet the depth of the convecting dynamo layers remains poorly constrained. While magnetic field observations, gravity data, and interior models each provide partial insight, they have not been combined into a single, self-consistent picture of the internal structure. Here, we develop a framework that links observed magnetic field strength with intrinsic heat flux and gravity-constrained interior structure using energy-based dynamo scaling laws. By relating the axial magnetic field strength to the convective power, we infer the radial thickness of the dynamo-generating region for both Jupiter and Saturn. The constants of proportionality in the scaling relations are derived using independent constraints from Earth observations, Jupiter observations, and numerical dynamo simulations. Applied to Jupiter, this framework shows how the inferred dynamo layer thickness is coupled to the outer boundary of the dynamo region. Thinner dynamo layers are predicted when the outer boundary shifts to shallower depths, and no solutions are possible when the outer boundary is less than 73% of Jupiter’s radius. These results constrain plausible geometries for future numerical dynamo simulations. Extending the analysis to Saturn, we find a thick, deep-seated dynamo layer with an outer radius at 42% of the radius to be most plausible. An alternative solution with an inner radius of the dynamo region at 60% of the planetary radius, as suggested by ring seismology models, requires a very thin dynamo layer, occupying only 2%–3% of the total radius.

Geosciences↗
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↗