Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary layer flow”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Coupling Between Alfvén Wave and Kelvin–Helmholtz Waves in the Low Latitude Boundary Layer

The Kelvin–Helmholtz (KH) instability of magnetohydrodynamic surface waves at the low latitude boundary layer is examined using both an eigenfrequency analysis and a time-dependent wave simulation. The analysis includes the effects of sheared flow and Alfvén velocity gradient. When the magnetosheath flows are perpendicular to the ambient magnetic field direction, unstable KH waves that propagate obliquely to the sheared flow direction occur at the sheared flow surface when the Alfvén Mach number is higher than an instability threshold. Including a shear transition layer between the magnetosphere and magnetosheath leads to secondary KH waves (driven by the sheared flow) that are coupled to the resonant surface Alfvén wave. There are remarkable differences between the primary and the secondary KH waves, including wave frequency, the growth rate, and the ratio between the transverse and compressional components. The secondary KH wave energy is concentrated near the shear Alfvén wave frequency at the magnetosheath with a lower frequency than the primary KH waves. Although the growth rate of the secondary KH waves is lower than the primary KH waves, the threshold condition is lower, so it is expected that these types of waves will dominate at a lower Mach number. Because the transverse component of the secondary KH waves is stronger than that of the primary KH waves, more efficient wave energy transfer from the boundary layer to the inner magnetosphere is also predicted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Numerical investigation of inertial sphere impaction on a cylinder in subcritical flow

We investigated the accuracy of a Lagrangian scheme for simulating both the motion and wall impaction of neutrally buoyant inertial spheres (8.85 ≤ St ≤ 18.04) moving in the subcritical regime (Re = 10,972 and 22,366) of flows past a circular cylinder. The accuracy of the eddy-resolving simulations of the flow field was verified based on available data at nearly comparable flow conditions. The accuracy of the eddy-resolving simulations of the flow field was verified based on available data at nearly comparable flow conditions. The accuracy of the particle tracking scheme was verified with respect to corresponding experimental data in a laboratory flume. Particle dispersion en route to the circular cylinder is mostly driven by the flow disturbance from the particle injector tube. The trends in collision rates are largely explained by the rate of particle entertainment into the boundary layer around the cylinder and the inertial condition of the spheres in this proximity region—both of which can be better assessed by post-processing of the numerical outcomes. Here, the primary incentive to undertake this work was the need for validating a particle-based method for quantifying the hydraulic stressors leading to bodily damage of fish passing flows influenced by hydraulic infrastructure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Adaptive mesh based combustion simulations of direct fuel injection effects in a supersonic cavity flame-holder

Here, we present high-fidelity reacting simulations of a supersonic cavity flame-holder configuration. The focus of this work is on flame stabilization brought about by varying the location of fuel injection in a cavity stabilized supersonic flow of air. Central to our approach is a compressible multi-species reacting flow solver that uses adaptive-mesh-refinement (AMR), enabling the resolution of flame, shock-waves, boundary-layers, and small-scale structures in the computational domain. Our analysis indicates that fuel injection closer to the ramp at the aft end of the cavity allows for greater mixing and lower peak temperatures compared to fuel injection upstream that is closer to the backward facing step of the cavity. This difference is mainly due to greater turbulent fluctuations generated from the shear-layer towards the cavity ramp, thereby enhancing the mixing of fuel and air. A low frequency oscillatory behaviour in heat-release and pressure was also observed for the upstream injection case while a much higher-frequency phenomena was observed in the near-ramp injection case. By identifying the important physical determinants of the combustion processes, this study illustrates a promising pathway to design and optimize direct fuel injection strategies in supersonic cavity flame-holders that can improve flame stability, combustion efficiency, and reduce emissions.

42 ENGINEERING↗
DOE OSTI2024

Measurement-driven large-eddy simulations of a diurnal cycle during a wake-steering field campaign

Abstract. High-fidelity flow modeling with data assimilation enables accurate representation of the wind farm operating environment under realistic, nonstationary atmospheric conditions. Capturing the temporal evolution of the turbulent atmospheric boundary layer is critical to understanding the behavior of wind turbines under operating conditions with simultaneously varying inflow and control inputs. This paper has three parts: the identification of a case study during a field evaluation of wake steering; the development of a tailored mesoscale-to-microscale coupling strategy that resolved local flow conditions within a large-eddy simulation (LES), using observations that did not completely capture the wind and temperature fields throughout the simulation domain; and the application of this coupling strategy to validate high-fidelity aeroelastic predictions of turbine performance and wake interactions with and without wake steering. The case study spans 4.5 h after midnight local time, during which wake steering was toggled on and off five times, achieving yaw offset angles ranging from 0 to 17°. To resolve nonstationary nighttime conditions that exhibited shear instabilities, the turbulence field was evolved starting from the diurnal cycle of the previous day. These background conditions were then used to drive wind farm simulations with two different models: an LES with actuator disk turbines and a steady-state engineering wake model. Subsequent analysis identified two representative periods during which the up- and downstream turbines were most nearly aligned with the mean wind direction and had observed yaw offsets of 0 and 15°. Both periods corresponded to partial waking on the downstream turbine, which had errors in the LES-predicted power of 4 % and 6 %, with and without wake steering. The LES was also able to capture conditions during which an upstream turbine wake induced a speedup at a downstream turbine and increased power production by up to 13 %.

17 WIND ENERGY↗
DOE OSTI2025

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗
DOE OSTI2020

Effects of Tethered Polymers on Dynamics of Nanoparticles in Unentangled Polymer Melts

Polymer-tethered nanoparticles (NPs) are commonly added to a polymer matrix to improve the material properties. Critical to the fabrication and processing of such composites is the mobility of the tethered NPs. Here, we study the motion of tethered NPs in unentangled polymer melts using molecular dynamics simulations, which offer a precise control of the grafted chain length N g and the number z of grafted chains per particle. As Ng increases, there is a crossover from particle-dominated to tethered-chain-dominated terminal diffusion of NPs with the same z. The mean squared displacement of loosely tethered NPs in the case of tethered-chain-dominated terminal diffusion exhibits two subdiffusive regimes at intermediate time scales for small z. The first one at shorter time scales arises from the dynamical coupling of the particle and matrix chains, while the one at longer time scales is due to the participation of the particle in the dynamics of the tethered chains. The friction of loosely grafted chains in unentangled melts scales linearly with the total number of monomers in the chains, as the friction of individual monomers is additive in the absence of hydrodynamic coupling. As more chains are grafted to a particle, hydrodynamic interactions between grafted chains emerge. As a result, there is a nondraining layer of hydrodynamically coupled chain segments surrounding the bare particle. Outside the nondraining layer is a free-draining layer of grafted chain segments with no hydrodynamic coupling. The boundary of the two layers is the stick surface where the shear stress due to the relative melt flow is balanced by the friction between the grafted and melt chains in the interpenetration layer. The stick surface is located further away from the bare surface of the particle with higher grafting density.

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 OSTI2021

Bridging Land-Surface Fluxes and Aerosol Concentrations to Triggering Convective Rainfall

During the period August 24-28, 2016, a delayed component the GO-Amazon Project, we launched nine altitude-controlled free balloons (CMET). Smaller than typical rawinsondes, CMET are equipped with altitude control, global communication via Iridium satellite, and aspirated sensors. The aims of our effort were to examine the interactions among convective boundary layer and dual river breeze circulations near the confluence of the Tapajos and Amazon Rivers in the eastern Basin. The week-long field campaign was timed to examine the reestablishment of the breeze circulations shortly after the passage of a strong instability line on August 22. Nine CMET were launched at the Curua-Una hydroelectric dam (2.8S; 54.3W), timed to encounter the Tapajos river breeze front by late afternoon. Soundings were made to establish the thickness of interface between the easterly trade and westerly Tapajos breeze circulation. Careful use of sounding strategies allowed these free balloons to track along the northerly channeled flow in the lowest 300 m above the River. Following the river encounter, balloons tracked to the west, sounding to describe the diurnal course of boundary layer in the forest west of the Tapajos River. The longest flight traveled more than 770 km over three days and twice rested overnight in the rain forest canopy. Ancillary data from surface climate and flux stations as well as the Santarem radiosonde, satellite images will be used to illustrate how the breeze circulations are seen near the surface and how they were disrupted by larger-scale events. Comparisons with HYSPLIT trajectories illustrate how sensitive real trajectories are to the refraction that the encounter with the breeze effects. A second effort was intensive analyses of data from the eddy flux tower at the Tapajós National Forest near Santarém Pará Brazil, with special emphasis on the role of clouds and canopy presence on light absorbance and photosynthetic activity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Closed Bay-Breeze Circulation and Its Lifecycle From TRACER With a New Orienteering Tape Recorder Diagram

A closed bay-breeze circulation (BBC) followed by a gulf-breeze front (GBF) was observed on 10 September 2022 during the Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) TRACER campaign. Using high-resolution X-band Scanning Cloud Radar (XSACR) and a newly developed orienteering tape recorder diagram, the study analyzed radar reflectivity and Doppler velocity to identify anomalies and track the evolution of these circulations. The BBC, a mesoscale system approximately 30-km long, 30-km wide, and 2-km deep, formed from enhanced horizontal convective rolls along the Galveston Bay coast, progressing northwestward under a 6 m s−1 onshore flow anomaly and reaching 1.5 km depth with return flow aloft. The inland penetration speed was 2 m s−1, driven by an observed 6-9 K land-water temperature contrast. The GBF, coupled to the BBC, intensified with additional southerly flow, penetrating further inland after the BBC collapsed. Passing over the TRACER field site, both fronts significantly impacted boundary layer thermodynamics, dynamics, and aerosol concentration. The BBC event exhibited four lifecycle stages-formation, development, maturation, and dissipation-driven by solar heating, wind field rotation, and interactions with convective eddies and the GBF. This study provides insights into the inland evolution of coastal breeze circulations and their interactions with environmental processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Observed volatilization behavior of silicon carbide in flowing hydrogen above 2000 K

In this work, the intrinsic compatibility of silicon carbide (SiC) and hydrogen (H 2 ) at high temperatures (2000-2473 K) and pressure near one atmosphere was evaluated through a combination of thermodynamic calculations and hot hydrogen exposure testing. Thermodynamic calculations predict the decomposition of SiC in a hydrogen environment to form free silicon (Si) and free carbon (C). Free Si is predicted to vaporize from the surface as a volatile species, while free C may interact with H 2 to form the hydrocarbons CH 4 (T < 2100 K) or C 2 H 2 (T > 2100 K). Coupons of high purity chemical vapor deposition (CVD) β-SiC were exposed to slowly flowing hydrogen at temperatures ranging between 2000 and 2473 K. SiC experienced active attack as the result of H 2 exposure, exhibiting linear weight loss kinetics and an Arrhenius dependence of weight loss on exposure temperature. The linear volatilization constant was experimentally evaluated to correspond with an activation energy of 370 ± 18 kJ/mol. Due to the dependence of observed corrosion rates on gas velocity, corrosion of SiC in flowing H 2 was determined to be governed by external mass transfer of volatile Si species through the boundary layer. Experimentally derived mass losses were in good agreement with mass losses predicted by a boundary layer limited gas diffusion model.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Strong Impact of Weak Horizontal Convergence on Continental Shallow Convection

Idealized large-eddy simulation (LES) is a basic tool for studying three-dimensional turbulence in the planetary boundary layer. LES is capable of providing benchmark solutions for parameterization development efforts. However, real small-scale atmospheric flows develop in heterogeneous and transient environments with locally varying vertical motions inherent to open multiscale interactive dynamical systems. These variations are often too subtle to detect them by state-of-the-art remote and in situ measurements, and are typically excluded from idealized simulations. The present study addresses the impact of weak [i.e., O(10-6) s-1] short-lived low-level large-scale convergence/divergence perturbations on continental shallow convection. The results show a strong response of shallow nonprecipitating convection to the applied weak large-scale dynamical forcing. Evolutions of CAPE, mean liquid water path, and cloud-top heights are significantly affected by the imposed convergence/divergence. In contrast, evolving cloud-base properties, such as the area coverage and mass flux, are only weakly affected. To contrast those impacts with microphysical sensitivity, the baseline simulations are perturbed assuming different observationally based cloud droplet number concentrations and thus different rainfall. For the tested range of microphysical perturbations, the imposed convergence/divergence provides significantly larger impact than changes in the cloud microphysics. Simulation results presented here provide a stringent test for convection parameterizations, especially important for large-scale models progressing toward resolving some nonhydrostatic effects.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Simulating Wind Around Isolated Buildings with the System for Atmospheric Modeling

A method to represent the flow around obstacles such as buildings is implemented in the System for Atmospheric Modeling (SAM). The method, which we call a Quasi-Solid Box Method (QSBM), introduces a simple modification to the anelastic equations that forces the flow to stagnate within the obstacle's boundaries. The performance of the modified SAM is evaluated using CEDVAL (Compilation of Experimental Data for Validation of Microscale Dispersion Models) wind tunnel measurements of the wind and tracer dispersion around a single rectangular building. All major observed features of the flow are quite well reproduced. The dispersion of the tracer released at the building base is also simulated. The results of simulation of a flow around a cubic building rotated by 45 degrees relative to the flow, and around a building in the form of a cylinder of aspect ratio of one, are also presented.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2025

Photoionization of seeded combustion products as a method of enhancing the efficiency of magnetohydrodynamic power generators

Here, in this study, we performed an experimental and computational investigation into the feasibility of utilizing photoionization to enhance the electrical conductivity of seeded oxy-fuel combustion products and improve the performance of magnetohydrodynamic (MHD) power generators. We applied a variety of optical and microwave diagnostics to study the ionization and recombination processes of potassium excited by an excimer laser in a high-velocity oxy-fuel free jet. Computational fluid dynamic (CFD) simulations were performed to model the thermophysical properties and species densities of the free jet. The CFD results were validated with position-dependent potassium concentration measurements. Electron recombination exponential lifetimes were measured through time-resolved microwave transmission. The experimental electron lifetimes were compared with lifetimes calculated from CFD-predicted species densities and literature recombination rates. It was determined that K + or O 2 are the most likely recombination partners for photoionized electrons. Time-resolved fluorescence measurements provided evidence of an ionization pathway involving a two-photon ionization of KOH . Finally, a zero-dimensional chemical kinetic model was developed to assess the fundamental viability of inducing a non-equilibrium electron population to provide a net energy return in combustion-driven MHD power generators. We determined that a high energy return is feasible for targeting electrode boundary layers with ultraviolet photoionization. We also found that photoionization could potentially lower the required temperature of the bulk gas flow.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Mixing Driven by Critical Reflection of Near-Inertial Waves over the Texas–Louisiana Shelf

Studies of internal wave-driven mixing in the coastal ocean have been focused on internal tides, while wind-driven near-inertial waves (NIWs) have received less attention in this regard. Here, this study demonstrates a scenario of NIW-driven mixing over the Texas–Louisiana shelf. Supported by a high-resolution simulation over the shelf, the NIWs driven by land–sea breeze radiate downward at a sharp front and enhance the mixing in the bottom boundary layer where the NIWs are focused because of slantwise critical reflection. The criterion for slantwise critical reflection of NIWs is ω=$\sqrt{f{^2}+s^{2}_{bot}N{^2}(1-s_{ρ}/s_{bot})}$ (where ω is the wave frequency, s bot is the bottom slope, and s ρ is the isopycnal slope), under the assumption that the mean flow is in a thermal wind balance and only varies in the slope-normal direction. The mechanism driving the enhanced mixing is explored in an idealized simulation. During slantwise critical reflection, NIWs are amplified with enhanced shear and periodically destratify a bottom boundary layer via differential buoyancy advection, leading to periodically enhanced mixing. Turbulent transport of tracers is also enhanced during slantwise critical reflection of NIWs, which has implications for bottom hypoxia over the Texas–Louisiana shelf.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A scale-wise analysis of intermittent momentum transport in dense canopy flows

We investigate the intermittent dynamics of momentum transport and its underlying time scales in the near-wall region of the neutrally stratified atmospheric boundary layer in the presence of a vegetation canopy. This is achieved through an empirical analysis of the persistence time scales (periods between successive zero-crossings) of momentum flux events, and their connection to the ejection–sweep cycle. Using high-frequency measurements from the GoAmazon campaign, spanning multiple heights within and above a dense canopy, the analysis suggests that, when the persistence time scales ( $t_p$ ) of momentum flux events from four different quadrants are separately normalized by $\varGamma _{w}$ (integral time scale of the vertical velocity), their distributions $P(t_p/\varGamma _{w})$ remain height-invariant. This result points to a persistent memory imposed by canopy-induced coherent structures, and to their role as an efficient momentum-transporting mechanism between the canopy airspace and the region immediately above. Moreover, $P(t_p/\varGamma _{w})$ exhibits a power-law scaling at times $t_{p}<\varGamma _{w}$ , with an exponential tail appearing for $t_{p} \geq \varGamma _{w}$ . By separating the flux events based on $t_p$ , we discover that around 80 % of the momentum is transported through the long-lived events ( $t_{p} \geq \varGamma _{w}$ ) at heights immediately above the canopy, while the short-lived ones ( $t_{p} < \varGamma _{w}$ ) only contribute marginally ( $\approx 20\,\%$ ). To explain the role of instantaneous flux amplitudes in momentum transport, we compare the measurements with newly developed surrogate data and establish that the range of time scales involved with amplitude variations in the fluxes tends to increase as one transitions from within to above the canopy.

Mechanics↗