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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The plasma–sheath transition and Bohm criterion in a high recycling divertor

The high recycling regime of a divertor is characterized by high plasma particle fluxes and low temperature at the target, where a strong hydrogen recirculation loop exists. Atomic processes in the high recycling regime, such as ion–neutral friction and radiation, can affect the plasma momentum and energy transport in the sheath transition region. Here, in this work, the plasma–sheath transition near a high recycling wall is investigated. The Bohm speed, which constraints the ion exit flow speed, is evaluated from a transport model that accounts for the effect of the anisotropic transport and atomic collisions in the transition layer. A first principles kinetic code vector particle-in-cell with the atomic collision package is used to investigate a 1D self-consistent slab plasma with a high recycling boundary for the tungsten and carbon divertors. The results demonstrate the accuracy of the Bohm speed model in predicting the ion exit flow speed in the transition region, as well as the reduction of the Bohm speed due to the ion–neutral friction. The effect of different wall materials, tungsten, and carbon, on the Bohm speed and near-wall plasma profile is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Unstable turbulent channel flow response to spanwise-heterogeneous heat fluxes: Prandtl's secondary flow of the third kind

Turbulent secondary flows are defined as Prandtl's secondary flow of the first or second kind, the former produced by stretching and/or tilting of vorticity, the latter produced via spatial heterogeneity of Reynolds stresses. Both mechanisms are instantaneously active within inertia-dominated wall turbulence; Reynolds stress spatial heterogeneity is required for Reynolds-averaged secondary flows. Spanwise-variable surface roughness can induce turbulent stress spatial heterogeneity in the spanwise–wall-normal plane and provide sustenance for streamwise-aligned mean secondary flows. Herein, we demonstrate that turbulent secondary flows can also be sustained by spanwise variability in the surface heat flux in unstably stratified turbulent channels, defined hereafter as Prandtl's secondary flow of the third kind. Support for this mechanism is established with scaling arguments, while large-eddy simulation is used to model inertia-dominated channel turbulence responding to a lower boundary with uniform aerodynamic/hydrodynamic roughness but spanwise-variable surface heat flux. Transport equations for streamwise vorticity and turbulent kinetic energy, k, outline the conditions needed for third-kind production: shear and buoyancy production over the elevated heat flux regions necessitates lateral entrainment of low- k fluid, inducing mean counter-rotating secondary cells aligned such that upwelling and downwelling occur over the high and low heat flux regions, respectively. Buoyancy-driven production of k alters aggregate flow response and thus is a distinctly different mechanism responsible for sustenance of secondary flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Fluctuations and Self-Organization in Plasma Boundary Layers (Final Report)

Final technical report on an experimental study of the boundary of an ionized gas (plasma) where it interacts with solid surfaces. The plasma self-organizes into a sheath which experiences plasma flow and is naturally unstable to various fluctuations which can be observed with advanced laser-spectroscopy techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING↗
DOE OSTI2021

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2024

Effect of Wettability Gradient on the Scale Formation in Falling Film Flow

Mitigation of scale formation and performance degradation remains a vital challenge for falling film evaporators in various industries. In this work, an experimental study of falling film flow on a horizontal tube is conducted to investigate the effects of wettability gradients on thermal, hydraulic, and fouling behavior. It is revealed that certain hydrophobic coating patterns, such as strip, ring, and grid patterns, lead to unwetted heat transfer area, which results in decreased heat transfer compared to fully wetted plain tube. By adjusting the geometry and position of the wettability gradient, the hybrid coating demonstrates improved heat transfer performance. Based on the characteristics of horizontal tube falling film flow, impinging jet, thin film flow, and liquid retention at the tube bottom, a hybrid coating pattern is developed to improve surface wetting and mitigate the scaling coverage. It is revealed that scale deposition is regulated by wettability gradient. Crystals tend to be dense and compact in hydrophilic areas, while they appear scattered or even absent in hydrophobic regions, depending on the dimension of the hydrophobic area. While at the hydrophilic/hydrophobic boundary, a noticeable scale thickness step is observed, which raises the potential for self-cleaning. Furthermore, the balance of minimization of scaling layer coverage and maximization of wetting area requires an optimal design in coating dimensions, for which a systemic study of both flow dynamics and fouling characteristics on the falling film is necessary in the future.

42 ENGINEERING↗
DOE OSTI2022

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Generalized Integral Boundary Layer Equations in Streamline Curvilinear Coordinates

The MEIT equations are derived in orthogonal streamline curvilinear coordinates defined by the inviscid velocity at the wall and the surface-normal direction. The resulting formulations is (1) is not restricted to flat-plate or axisymmetric geometries, (2) accommodates vortical inviscid flows (e.g., induced by curved shocks) and wall mass blowing, and (3) allows for an arbitrary definition of the edge state. The streamline ODEs are also extended to a surface PDE framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Grain Boundary-Limited Thermal Transport in Suspended Thin Graphite across an Unexplored Thickness Regime

Here, we present systematic thermal conductivity (κ) measurements of suspended thin graphite ribbons, 234–527 nm thick, using a four-probe 3ω method. Unlike recent reports of phonon hydrodynamics and exceptionally high κ in micrometer-thick graphite ( Science, 2020), we observe significantly lower κ and no signatures of collective phonon flow in this intermediate thickness regime. Instead, our measured κ lies between few-layer graphene and bulk graphite. These results agree with a first-principles-informed Peierls–Boltzmann transport model with spatially resolved Monte Carlo sampling. Additionally, the temperature for the peak κ shifts lower with increasing thickness, due to the interplay of phonon-boundary and phonon-isotope scattering. Incorporating grain boundary scattering into simulations is necessary to replicate the experimental trends. These findings delineate the boundary between ballistic, hydrodynamic, and diffusive transport regimes in graphite and underscore the dominant role of disorder and geometry in phonon transport in quasi-two-dimensional materials, offering insights for nanoscale thermal management.

Boltzmann transport model↗
DOE OSTI2024

Contact solitons as spin pistons in one-dimensional ferromagnetic channels

Ferromagnetic channels subject to spin injection have been theoretically shown to sustain dissipative exchange flows (DEFs). In the strong injection regime, a soliton is stabilized at the injection site, which has been termed a contact soliton DEF or CS-DEF. Here, in this work, we investigate the modulation of CS-DEFs as a mechanism to inject magnons into a DEF. By varying the injection current, the parameters connecting the soliton with the DEF via a boundary layer are varied. These lead to a modification in the DEF which is interpreted as a spin piston. The injected magnons follow the expected dispersion relation of DEFs, akin to the Bogoliubov–de Gennes dispersion relation. As a result, this work demonstrates that changes in the injected current can pump magnons along a ferromagnetic channel via dissipative exchange flows.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗
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 OSTI2024

Bubble departure and sliding in high-pressure flow boiling of water

Bubble growth, departure and sliding in low-pressure flow boiling has received considerable attention in the past. However, most applications of boiling heat transfer rely on high-pressure flow boiling, for which very little is known, as experimental data are scarce and very difficult to obtain. In this work, we conduct an experiment using high-resolution optical techniques. By combining backlit shadowgraphy and phase-detection imaging, we track bubble shape and physical footprint with high spatial ($6\,\mathrm {\mu }{\rm m}$) and temporal ($33\,\mathrm {\mu }{\rm s}$) resolutions, as well as bubble size and position as bubbles nucleate and slide on top of the heated surface. We show that at pressures above 1 MPa bubbles retain a spherical shape throughout the growth and sliding process. We analytically derive non-dimensional numbers to correlate bubble velocity and liquid velocity throughout the turbulent boundary layer and predict the sliding of bubbles on the surface, solely from physical properties and the bubble growth rate. We also show that these non-dimensional solutions can be leveraged to formulate elementary criteria that predict the effect of pressure and flow rate on bubble departure diameter and growth time.

Mechanics↗
DOE OSTI2021

Tuning plasma parameters to control reactive species fluxes to substrates in the context of plasma catalysis

The key reactive plasma-produced gas phase species responsible for the enhanced conversion of chemicals in plasma catalysis compared to thermal catalysis have to date not been identified. This outstanding question is mainly due to the inherent large variety of plasma-produced species and the challenge of controlling and measuring the flux of each constituent of the cocktail of reactive species to a (catalytic) substrate. In this paper, we explore the possibility to control the dominant reactive species fluxes, relevant for plasma-catalysis, to a substrate in the effluent of an RF driven Ar-O 2 plasma jet. The absolute species densities of the major reactive species (O, O 2 (a 1 Δ g ), O 3 and ions) were quantified by molecular beam mass spectrometry (MBMS) to assess the possibility of using treatment distance, O 2 admixture concentration, plasma dissipated power, RF modulation frequency and duty cycle as well as the feed gas flow rate to alter the dominant species densities. Selected experimental results were also compared with a pseudo-1D plug flow model. The short-lived and long-lived species can be effectively separated by changing the treatment distance and the RF modulation frequency. Furthermore, adjusting the O 2 admixture concentration enables to change the ratio of the O 2 (a 1 Δ g ) and O 3 density. The changes in the trend of ion and O flux were found to be very similar for nearly all investigated parameters. Nonetheless the gas flow rate was able to significantly change the ratio of the O and ion density in the plasma jet effluent. Here, the impact of the surface-dependent loss probability and boundary layer reactions on the species flux to a substrate and how this qualitatively relates to the MBMS density measurements is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗