SEARCH · Engineering Papers

Results for “evolving boundary layer”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2021

High-Fidelity Simulation of Turbulent Flow Past a Gaussian Bump

A spanwise-periodic computation of a turbulent flow past a Gaussian bump is performed in the form of a hybrid direct numerical simulation and wall-resolved large-eddy simulation. A fourth-order spatially-accurate flow solver is employed to perform the simulation, using 10.2 billion grid points for a Reynolds number of 170000 based on the bump height. The key findings from the simulation are reported in the acceleration and deceleration flow regions associated with the bump shape. Significant anisotropy in the normal Reynolds stresses, along both the wall-normal and streamwise directions, is observed within the acceleration region. The ratio between the Reynolds shear stress and turbulent kinetic energy in that region also experiences significant deviations from the norms of a zero pressure gradient turbulent boundary layer. The chosen Reynolds number generates strong flow separation in the adverse pressure gradient region, which is in contrast with a previous simulation at half the Reynolds number that only indicated incipient separation. An internal layer generated in the acceleration region evolves into a free shear layer that develops in the deceleration region and separates. Proper modeling of this inner layer appears crucial to predict the flow separation. Surface curvature effects on the attached flow development are also discussed.

Turbulent Boundary Layer↗

NASA NTRS2021

Tracking Smoke from a Prescribed Fire and its Impacts on Local Air Quality using Temporally Resolved GOES-16 ABI Aerosol Optical Depth (AOD)

Aerosol optical depth (AOD) retrieved from the GOES-16 Advanced Baseline Imager (ABI) was used to track a smoke plume from a prescribed fire in northeastern Virginia on 8 March 2020. Weather and atmospheric conditions created a favorable environment to transport the plume through the Washington, D.C., and Baltimore, Maryland, metro areas in the afternoon and concentrate smoke near the surface, degrading air quality for several hours. ABI AOD with 5-min temporal resolution and 2-km spatial resolution definitively identified the timing and geographic extent of the plume during daylight hours. Comparison to AERONET AOD indicates that ABI AOD captured the relative change in AOD due to passage of the smoke, with a mean absolute error of 0.047. Ground-based measurements of fine particulate matter (PM2.5) confirm deteriorations in air quality coincident with the progression of the smoke. Ceilometer aerosol backscatter profiles verify plume transport timing and indicate that smoke aerosols were well mixed in a shallow boundary layer. This event illustrates the advantages of using multiple datasets to analyze the impacts of aerosols on ambient air quality. Given the quickly evolving nature of the event over several hours, ABI AOD provided information for the public and decision-makers that was not available from any other source, including polar-orbiting satellite sensors. This study suggests that PM2.5 concentrations estimated from ABI AOD can be used to fill in the gaps in nationwide regulatory PM2.5 monitor networks and may be a valuable addition to EPA’s PM2.5 NowCast of current air quality conditions.

satellite observations↗

DOE OSTI2024

The Association Between Cloud Droplet Number over the Summer Southern Ocean and Air Mass History

The cloud properties and governing processes in Southern Ocean marine boundary layer clouds have emerged as a central issue in understanding the Earth's climate sensitivity. While our understanding of Southern Ocean cloud feedbacks have evolved in the most recent climate model intercomparison, the background properties of simulated summertime clouds in the Southern Ocean are not consistent with measurements due to known biases in simulating cloud condensation nuclei concentrations. This paper presents several case studies collected during the Capricorn 2 and Marcus campaigns held aboard Australian research vessels in the Austral Summer of 2018. Combining the surface–observed cases with MODIS data along forward and backward air mass trajectories, we demonstrate the evolution of cloud properties with time. These cases are consistent with multi–year statistics showing that long trajectories of air masses over the Antarctic ice sheet are critical to creating high droplet number clouds in the high latitude summer Southern Ocean. We speculate that secondary aerosol production via the oxidation of biogenically derived aerosol precursor gasses over the high actinic flux region of the high latitude ice sheets is fundamental to maintaining relatively high droplet numbers in Southern Ocean clouds during Summer.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Towards a Marine Stratus Climatology on Drizzle Occurrence from CALIPSO

Marine stratus are a predominant feature of our planet with the annual mean coverage exceeding 20%. They strongly reflect sunlight, yet exert only a modest effect on outgoing infrared radiation, providing a significant net cooling to the Earth’s radiative balance. Their formation is coupled to boundary layer circulations that are driven, in part, by cloud top radiative cooling and evaporative cooling from precipitation in downdrafts. Understanding how these cloud systems evolve as the climate changes is a key question that requires additional information on their lifecycle and microphysical properties to accurately represent their behavior in global circulation models. From a large-scale perspective, insight into the microphysical properties of marine stratus at cloud top can be realized through estimates of the effective radius (Re) of the droplet size distributions derived from MODIS observations. Estimates on the occurrence of rain/drizzle are available from CloudSat. Together these observations indicate that precipitation frequently occurs in clouds with higher cloud top Re. This relationship is consistent with the well documented shift in cloud top droplet size distributions towards fewer, yet larger droplets prior the onset of precipitation. Here we report on a new and complementary set observations from the CALIPSO mission. The approach derives an extinction-to-backscatter ratio (Sc, also known as the cloud lidar ratio) using an established relationship that depends on observations of the lidar attenuated backscatter and volume depolarization ratio within the cloud. Because Sc is strongly and inversely related to Re, a change in the derived Sc from higher to lower values corresponds to a change in the droplet size distribution as seen by MODIS. This change in the lidar signals at cloud top clearly identifies clouds that are capable of precipitation. The presentation provides a brief overview of the approach for deriving Sc and compares CALIOP-derived Sc with observations from other techniques. CALIOP classifications of drizzling clouds, based on the retrieved values Sc, are compared to independent, collocated assessments of drizzle occurrence reported in the standard CloudSat data products. Regional and seasonal comparisons highlight the strengths and weaknesses of the two sensors. A machine learning approach that combines information from both CALIOP and CloudSat showcases possible improvements in the global identification of scenes likely to contain rain-bearing clouds.

CALIPSO↗

NASA NTRS1990

Turbulent supersonic convection in three dimensions

Previous numerical calculations of two-dimensional, compressible convection are extended to three dimensions, using a higher order Godunov scheme. The results show that the flow readily becomes supersonic in the upper boundary layer, where shock structures form intermittently in the vicinity of the strong downflow lanes. The convection as a whole is strongly time-dependent and evolves on a time scale comparable to the sound crossing time. The motions in the upper layers are characterized by the rapid expansion of the upward-moving fluid elements. In the interior, most of the heat is carried by a small fraction of the fluid residing in strong, highly coherent downflows. The remaining fluid is dominated by small-scale, disorganized turbulent motions.

Malagoli, Andrea↗

NASA NTRS1991

Twenty-five years of aerodynamic research with IR imaging: A survey

Infrared imaging used in aerodynamic research evolved during the last 25 years into a rewarding experimental technique for investigation of body-flow viscous interactions, such as heat flux determination and boundary layer transition. The technique of infrared imaging matched well its capability to produce useful results, with the expansion of testing conditions in the entire spectrum of wind tunnels, from hypersonic high-enthalpy facilities to cryogenic transonic wind tunnels. With unique achievements credited to its past, the current trend suggests a change in attitude towards this technique: from the perception as an exotic, project-oriented tool, to the status of a routine experimental procedure.

Gartenberg, Ehud↗

NASA NTRS2007

On the Evolution of Localized Disturbances and their Spanwise Interactions Leading to Breakdown

Localized disturbances in a laminar boundary layer represent a more realistic model of transition than the extensively studies, two or quasi three-dimensional perturbations regardless of the fact if they evolve in a linear manner or not. Localized disturbances can originate by surface imperfections, insects or dust. The disturbances can be harmonic (i.e. containing a single frequency and a complete set of spanwise wave numbers) or Pulsed (i.e. containing a band of streamwise and spanwise wave numbers). At sufficiently low amplitudes localized disturbances behave according to a linear stability model. It is highly probably that in a natural transition process such localized disturbances will overslap and interact. These interactions could either delay transition because of a partial wave cancellation resulting in an attenuation of the disturbance, or adversely enhance it by promoting nonlinear interactions. The nonlinearity could be simply amplitude dependent or cause a triad resonance. Nonlinear processes in a wave packet lead to breakdown and to the formation of turbulent spots. When the amplitude of the harmonic disturbance saturates, nonlinear processes widen the band of the lower amplified frequencies adjacent to the frequency of excitation. Experimental results describing the spanwise interactions of harmonic and pulsed localized disturbances leading to breakdown will be presented and discussed. A comparison to the evolution and breakdown of a single localized disturbance will be provided.

Seifert, Avi↗

DOE OSTI2026

Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon

The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.

Lyu, Yanjun [Department of Materials Science and E↗

NASA NTRSDate not supplied

Mesh Deformation Boundary Conditions for Three-Dimensional Ablation Solvers

A method for determining boundary node displacements for three-dimensional ablation solvers is developed. The approach is applicable to a wide range of mesh motion algorithms used to determine internal node displacements. The method leverages radial basis functions to analytically define the geometry to which the boundary nodes should be constrained. Constrained optimization, elasticity, and spring analogy solvers are utilized to redistribute the nodes on the evolved geometry surface while maintaining mesh quality. The method is implemented in the CHarring Ablator Response (CHAR) code, and a realistic demonstration case is presented on the Boundary Layer Transition (BOLT) flight experiment configuration. The method is shown to be robust for highly-curved geometries with large domain deformations.

Adam J Amar↗

NASA NTRS1985

Evolution of a curved vortex filament into a vortex ring

The deformation of a hairpin-shaped vortex filament under self-induction and in the presence of shear is studied numerically using the Biot-Savart law. It is shown that the tip region of an elongated hairpin vortex evolves into a vortex ring and that the presence of mean shear impedes the process. Evolution of a finite-thickness vortex sheet under self-induction is also investigated using the Navier-Stokes equations. The layer evolves into a hairpin vortex which in turn produces a vortex ring of high Reynolds stress content. These results indicate a mechanism for the generation of ring vortices in turbulent shear flows, and a link between the experimental and numerical observation of hairpin vortices and the observation of ring vortices in the outer regions of turbulent boundary layers.

Moin, P.↗

NASA NTRS1986

Evolution of a curved vortex filament into a vortex ring

Moin, P.↗

NASA NTRS1994

Auroral activity associated with Kelvin-Helmholtz instability at the inner edge of the low-latitude boundary layer

Auroral activity occurred in the late afternoon sector (approx. 16 MLT) in the northern hemisphere during the passage at Earth of an interplanetary magnetic cloud on January 14, 1988. The auroral activity consisted of a very dynamic display which was preceded and followed by quiet auroral displays. During the quiet displays, discrete rayed arcs aligned along the geomagnetic L shells were observed. In the active stage, rapidly evolving spiral forms centered on magnetic zenith were evident. The activity persisted for many minutes and was characterized by the absence of directed motion. They were strongly suggestive of intense filaments of upward field-aligned currents embedded in the large-scale region 1 current system. Distortions of the flux ropes as they connect from the equatorial magnetosphere to the ionosphere were witnessed. We assess as possible generating mechanisms three nonlocal sources known to be associated with field-aligned currents. Of these, partial compressions of the magnetosphere due to variations of solar wind dynamic pressure seem an unlikely source. The possibility that the auroral forms are due to reconnection is investigated but is excluded because the active aurora were observed on the closed field line region just equatorward of the convection reversal boundary. To support this conclusion further, we apply recent results on the mapping of ionospheric regions to the equatorial plane based on the Tsyganenko 1989 model (Kaufmann et al., 1993). We find that for comparable magnetic activity the aurora map to the equatorial plane at X(sub GSM) = approx. 3 R(sub E) and approx. 2 R(sub E) inward of the magnetopause, that is, the inner edge of the boundary layer close to dusk. Since the auroral forms are manifestly associated with magnetic field shear, a vortical motion at the equatorial end of the flux rope is indicated, making the Kelvin-Helmholtz instability acting at the inner edge of the low-latitude boundary layer the most probable generating source.

Farrugia, C. J.↗

NASA NTRS2015

A-Train Based Observational Metrics for Model Evaluation in Extratropical Cyclones

Extratropical cyclones contribute most of the precipitation in the midlatitudes, i.e. up to 70 during winter in the northern hemisphere, and can generate flooding, extreme winds, blizzards and have large socio-economic impacts. As such, it is important that general circulation models (GCMs) accurately represent these systems so their evolution in a warming climate can be understood. However, there are still uncertainties on whether warming will increase their frequency of occurrence, their intensity and how much rain or snow they bring. Part of the issue is that models have trouble representing their strength, but models also have biases in the amount of clouds and precipitation they produce. This is caused by potential issues in various aspects of the models: convection, boundary layer, and cloud scheme to only mention a few. In order to pinpoint which aspects of the models need improvement for a better representation of extratropical cyclone precipitation and cloudiness, we will present A-train based observational metrics: cyclone-centered, warm and cold frontal composites of cloud amount and type, precipitation rate and frequency of occurrence. Using the same method to extract similar fields from the model, we will present an evaluation of the GISS-ModelE2 and the IPSL-LMDZ-5B models, based on their AR5 and more recent versions. The AR5 version of the GISS model underestimates cloud cover in extratropical cyclones while the IPSL AR5 version overestimates it. In addition, we will show how the observed CloudSat-CALIPSO cloud vertical distribution across cold fronts changes with moisture amount and cyclone strength, and test if the two models successfully represent these changes. We will also show how CloudSat-CALIPSO derived cloud type (i.e. convective vs. stratiform) evolves across warm fronts as cyclones age, and again how this is represented in the models. Our third process-based analysis concerns cumulus clouds in the post-cold frontal region and how their amount relates to the stability of the boundary layer. This test uses Aqua cloud and vertical atmospheric profiles and when applied to the model output can help assess the accuracy of the convection, boundary layer and cloud scheme.

Convection clouds↗

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

Irradiation-Induced Structural Disorder and Its Influence on the Mechanical Response of Polycrystalline MoS2

Molybdenum disulfide (MoS2) thin films are widely used as dry-film lubricants and protective coatings in aerospace and other radiation-exposed environments. Conventional synthesis routes produce polycrystalline films whose grain boundaries and other native defects cause their mechanical and tribological behavior to differ substantially from that of ideal single crystals. Under irradiation, these films progressively evolve from polycrystalline structures, composed of layered MoS2 grains, into highly disordered and eventually amorphous structures, altering both their tribological performance and mechanical integrity. Here, we employ reactive atomistic simulations to investigate irradiation-driven structural evolution in bulk polycrystalline MoS2. Using controlled primary knock-on atom (PKA) events, we characterize the progressive transition from a polycrystalline microstructure to an amorphous network by tracking defect accumulation and structural disorder. We then establish how this transition modifies the dominant deformation mechanisms and the temperature-dependent tensile response. Specifically, irradiation suppresses interlayer sliding and delamination, mechanisms which facilitate the deformation of the pristine polycrystal, resulting in defect-induced hardening. Broadly, our results establish direct process–structure–property relationships linking irradiation-induced defect accumulation, microstructural evolution, deformation mechanisms, and mechanical behavior, providing an atomistic framework for understanding the structural integrity and long-term reliability of irradiated MoS2 coatings.

Moore, Daniel [Sandia National Laboratories (SNL)]↗

NASA NTRS1994

Direct simulation of a self-similar plane wake

Direct simulations of two time-developing turbulent wakes have been performed. Initial conditions for the simulations were obtained from two realizations of a direct simulation of a turbulent boundary layer at momentum thickness Reynolds number 670. In addition, extra two dimensional disturbances were added in one of the cases to mimic two dimensional forcing. The unforced wake is allowed to evolve long enough to attain self similarity. The mass-flux Reynolds number (equivalent to the momentum thickness Reynolds number in spatially developing wakes) is 2000, which is high enough for a short k(exp -5/3) range to be evident in the streamwise one dimensional velocity spectrum. Several turbulence statistics have been computed by averaging in space and over the self-similar period in time. The growth rate in the unforced flow is low compared to experiments, but when this growth-rate difference is accounted for, the statistics of the unforced case are in reasonable agreement with experiments. However, the forced case is significantly different. The growth rate, turbulence Reynolds number, and turbulence intensities are as much as ten times larger in the forced case. In addition, the forced flow exhibits large-scale structures similar to those observed in transitional wakes, while the unforced flow does not.

Moser, Robert D.↗

NASA NTRS1993

Applying and validating the RANS-3D flow-solver for evaluating a subsonic serpentine diffuser geometry

Subsonic inlet ducts for advanced, high-performance aircraft are evolving towards complex three-dimensional shapes for reasons of overall integration and weight. These factors lead to diffuser geometries that may sacrifice inlet performance, unless careful attention to design details and boundary layer management techniques are employed. The ability of viscous computational fluid dynamic (CFD) analysis of such geometries to aid the aircraft configurator in this complex design problem is herein examined. The RANS-3D Reynolds-Averaged Navier-Stokes solver is applied to model the complex flowfield occurring in a representative diffuser geometry and the solutions are compared to experimental results from a static test of the inlet duct. The computational results are shown to compare very favorably with experimental results over a range of mass flow rates, including those involving large amounts of separation in the diffuser. In addition, a novel grid topology is presented, and two turbulence models are evaluated in this study as part of the RANS-3D code.

Fletcher, Michael J.↗

DOE OSTI2024

Asynchronicity in opposed-piston RCMs: Does it matter?

Rapid Compression Machines (RCMs) are widely utilized to study combustion phenomena at engine-relevant conditions, and significant efforts are typically made to create a quiescent environment, particularly for investigations of autoignition chemistry. Opposed-piston configurations can be advantageous due to shorter compression times and reduced surface area to volume ratios. Each side must be actuated simultaneously, but this can be challenging in practice. These devices, like most RCMs, utilize hydraulics for actuation, speed control and arrestation of the piston at the end of the stroke; there is no mechanical control or linkage of the two piston trajectories. To quantify the magnitudes and effects of piston asynchronous behavior, this work employs both detailed experimental measurements and, for the first time, high-fidelity, Direct Numerical Simulation (DNS). The boundary conditions are carefully considered applying insight from high-resolution linear variable differential transformer (LVDT) measurements of the piston trajectory and a zero-dimensional kinematics model of the piston-shaft assembly. Sufficient resolution in the piston crevice region is used. The complicated fluid dynamical behavior that can evolve during piston compression and the ensuing delay processes due to offset timings from t offset = 0-10 ms is elucidated. It is found that near t offset = 6 ms and beyond, the boundary layer on the face of the first-seating piston can be sufficiently perturbed, due initially to reemergence of gas from the crevice of the firstseating piston, so that the adiabatic core can become degraded at long ignition delay times. Substantial mixing of colder gas into the interior of the reaction chamber can alter the measurements, similar to effects previously observed for improper piston crevice configuration. In conclusion, experimental techniques to mitigate asynchronous behavior are discussed and demonstrated.

33 ADVANCED PROPULSION SYSTEMS↗