Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow velocity”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model

The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Cylinders and jets in crossflow: Wake formations as a result of varying geometric proximities

The combined flow physics of several canonical flow configurations is experimentally studied. Here, we analyze an array of jets issuing into a crossflow, then immediately navigating past an array of cylinders. This is achieved with a 2 × 3 triangular pattern of jets and symmetric cylinders at three jets to crossflow velocity ratios, enabling near-complete optical access of each jet, with velocities measured by time-resolved particle image velocimetry. Jet trajectories reveal that each configuration adheres to a power-law trend and that greater penetration is achieved by the downstream and confined jets compared to the more conventional upstream one. Recirculation regions of the upstream and downstream jets are nearly independent, with the confined jet encompassing regions of overlap with both. Turbulent statistics reveal the influence of geometric placement and velocity ratio on the time-averaged velocity, anisotropy, and Reynolds stresses incurred by each jet. Galilean decomposition utilizes a supplemental crossflow-only velocity field to delineate the influence of each jet's low- and high-pressure regions on the otherwise uniform stream. Proper orthogonal decomposition suggests that increased jet penetration decreases the number of modes required for truncation in the investigated spanwise plane. Vortex identification algorithms are applied to the reconstructed flow fields, reaffirming that with increasing velocity ratio, the jets generate vortices of their own in similar statistical formations as the cylinders. This investigation provides a foundation to aid future modeling efforts toward characterizing flow physics of importance in designing and passively controlling transverse jets with varying blockage proximities in a crossflow.

Kristo, Paul J. (ORCID:0000000225605175)↗
DOE OSTI2022

Substructure at High Speed. II. The Local Escape Velocity and Milky Way Mass with Gaia eDR3

Abstract Measuring the escape velocity of the Milky Way is critical in obtaining the mass of the Milky Way, understanding the dark matter velocity distribution, and building the dark matter density profile. In Necib & Lin, we introduced a strategy to robustly measure the escape velocity. Our approach takes into account the presence of kinematic substructures by modeling the tail of the stellar distribution with multiple components, including the stellar halo and the debris flow called the Gaia Sausage (Enceladus). In doing so, we can test the robustness of the escape velocity measurement for different definitions of the “tail” of the velocity distribution and the consistency of the data with different underlying models. In this paper, we apply this method to the Gaia eDR3 data release and find that a model with two components is preferred, although results from a single-component fit are also consistent. Based on a fit to retrograde data with two bound components to account for the relaxed halo and the Gaia Sausage, we find the escape velocity of the Milky Way at the solar position to be v esc = 445 − 8 + 25 km s −1 . A fit with a single component to the same data gives v esc = 472 − 12 + 17 km s −1 . Assuming a Navarro−Frenck−White dark matter profile, we find a Milky Way concentration of c 200 = 19 − 7 + 11 and a mass of M 200 = 4.6 − 0.8 + 1.5 × 10 11 M ⊙ , which is considerably lighter than previous measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A remarkable CHF of 345W/cm 2 is achieved in a wicked-microchannel using HFE-7100

Boiling heat transfer of dielectric fluids is a promising cooling technique for thermal management of microelectronic systems. However, the critical heat flux achievable is generally low, because of the poor thermophysical properties of these fluids. Here, to address this dilemma, we propose a new cooling concept to substantially enhance global liquid supply during phase-change process using enhanced capillary-driven force. Additionally, dedicated vapor pathways are designed among a bank of micro-pillars to facilitate vapor removal. In this work, new wicks comprised of silicon micro-pinfin arrays are explored to significantly enhance the flow boiling heat transfer performance. To examine the functionalities of this wick, experiments on HFE-7100 were carried out with mass velocities varying from 247 to 3,465 kg/m 2 s. To explore the enhancement mechanisms and to analyze the capillary-assisted flow boiling process, visualization studies were conducted. The results indicate that sustainable evaporation induced by wick microstructures and efficient liquid supply are the enhancement mechanisms compared to parallel microchannels with solid walls. It is found that the overall heat transfer coefficient is substantially increased up to 75%. Remarkably, a high critical heat flux (CHF) of approximately 345 W/cm 2 is recorded at G = 3,465 kg/m 2 s at coolant inlet temperature of ~20 °C. Equally importantly, this noticeable enhancement of CHF value is associated with drastically decreased pressure drops compared to microchannels decorated with μ-pinfin fences.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Laser Measurements and Modeling of Shock Tunnel Freestream Velocity and Multispecies Thermal Nonequilibrium

Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N 2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N 2 and O 2 , which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. Furthermore, this extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N 2 exhibits the most nonequilibrium, followed by O 2 and NO. The CFD model captures this trend, although it consistently overpredicts N 2 and O 2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.

Aerodynamics↗
DOE OSTI2022

Observations on the deformation of metal microspheres in shock-driven polymer flows

We report that solid particles can be fragmented by a fast-moving fluid if their velocity difference is great enough, such as during the atmospheric entry of meteoroids or the shock compression of engineered particulate composites. The extent of particle deformation and breakup in such systems is poorly understood because the necessary extreme conditions make observation difficult and data scarce. To meet this need, experiments combining ultrafast synchrotron-based radiography with plate impact loading were performed at the dynamic compression sector at the advanced photon source. Metal microspheres of several densities and strengths (Au, Ta, and W) were placed inside a polymer matrix. A planar shock wave was then produced in the polymer by the impact of a gun-launched flyer plate. X-ray images of the resulting flow were collected at ~150ns intervals. These images document the progression of particle deformation across a range of flow conditions and particle materials. They show that the extent of deformation is sensitive to the ratio of drag stress to particle strength. The deforming particle's shape is determined by the initial shock–particle interaction, fluid stagnation pressure, and vorticity, each acting on its own timescale. A set of scaling relationships is presented to capture these observations and enable comparison with prior hydrodynamic data. The result is a framework for predicting the conditions under which strong particles are severely deformed by a shock-driven flow.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

A toroidal asymptotic matching model of the response of a tokamak plasma to a static resonant magnetic perturbation (RMP) is used to simulate the n = 3 RMP-induced edge-localized-mode-suppression windows in q 95 that are evident when the plasma current is slowly ramped in DIII-D discharge #145380. All quantities employed in the simulation are derived from experimental measurements, apart from the neutral particle data. Three cases are considered. In the first case, the natural frequencies of tearing modes resonant in the plasma are determined by the ion flows at the corresponding resonant surfaces, which is the prediction of nonlinear tearing mode theory. In the second case, the natural frequencies are determined by the local E x B velocities at the resonant surfaces. In the third case, the natural frequencies are determined by the electron flows at the resonant surfaces, which is the prediction of linear tearing mode theory. The second case gives the best agreement between the simulations and the experimental observations. The first and third cases only lead to partial agreement between the simulations and the observations. In the first case, the lack of complete agreement may be a consequence of using an inaccurate assumption for the neutral particle distribution in the pedestal. In the third case, the lack of complete agreement is probably due to the fact that the response of a tokamak plasma to an RMP is not accurately described by linear tearing mode theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory of the tertiary instability and the Dimits shift within a scalar model

The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry--Horton system has been proposed by St-Onge [J. Plasma Phys. {\bf 83}, 905830504 (2017)] as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$, where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U''$. Then, depending on $U''$, the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $U'$. Our analytic predictions are in agreement with direct numerical simulations of the modified Terry--Horton system.

fusion plasmas↗
DOE OSTI2021

The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

42 ENGINEERING↗
DOE OSTI2022

Multi-fidelity information fusion with concatenated neural networks

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

42 ENGINEERING↗
DOE OSTI2025

Thermal Performance of Triply Periodic Minimal Surface Lattice Structures in Single-Phase Dielectric Fluid Cooling of Power Electronics

Additive manufacturing has transformed thermal management by enabling the production of complex, optimized geometries that conventional manufacturing methods cannot achieve. This study investigates the single-phase convective heat transfer performance of gyroid triply periodic minimal surface (TPMS) lattice structures with functional porosity. TPMS structures provide high surface area to volume ratios and are amenable to 3D printing. A gyroid numerical model was created and validated against an existing experimental study with a similar feature size to the investigated geometries. The TPMS structure has a periodic width of 1.6 mm, a length of 10 mm, and a height of 4 mm, with a functional porosity ranging from 0.5 to 0.8, decreasing with distance from the heated surface. Three different flow configurations were examined for an inlet fluid temperature of 70 °C. The inlet velocities range from 0.01 to 1.2 m/s, corresponding to a Reynolds number range of 10–900 with a heat flux of 50 W/cm 2 applied at the base. AmpCool ® AC-110 dielectric fluid (Prandtl number 59.5) was used as the coolant. Thermal performance and friction characteristics were studied for the three flow orientations. The parallel flow configuration was identified as the most efficient for heat removal. A detailed analysis of the numerical results highlights the underlying physics behind the thermal performance differences among the flow configurations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Evaluation of external tubular configurations for a high-temperature chloride molten salt solar receiver operating above 700°C

Next-generation concentrating solar power (CSP) tower technologies target operating temperatures exceeding 700°C to increase the thermal-to-electric conversion efficiency. Molten chloride salts are one possible alternative to current commercial molten nitrate salts to enable the higher operating temperature. This paper analyzes the predicted optical, thermal-fluids, and structural performance of traditional external tubular solar receiver configurations applied with a chloride salt heat transfer fluid (HTF) and inlet/outlet temperatures of 500°C/735°C, and considers sensitivity analysis and optimization relative to receiver sizing, tube sizing, number of panels, flow circuit configurations, and solar flux concentration under constraints on internal velocity, pressure drop, wall thickness, and required creep-fatigue lifetime. The high temperature conditions increase the significance of inelastic deformation mechanisms such as creep relative to that expected in commercial 565°C nitrate salt designs. High-temperature creep and creep-fatigue damage in the metal alloy tubes are the key factors that limit allowable solar flux concentration and achievable receiver thermal efficiency at the near-800°C wall temperature conditions. For a traditional external cylindrical receiver configuration, the design parameters and conditions capable of satisfying all constraints produced, at best, a design point receiver efficiency of 78.2%, or 80.5% when excluding receiver intercept efficiency. Variation in the optimal receiver performance relative to uncertainty in the binding maximum velocity, minimum wall thickness, and minimum lifetime constraints is presented.

14 SOLAR ENERGY↗
DOE OSTI2022

Water‐Induced Diamond Formation at Earth's Core‐Mantle Boundary

Abstract The carbon and water cycles in the Earth's interior are linked to key planetary processes, such as mantle melting, degassing, chemical differentiation, and advection. However, the role of water in the carbon exchange between the mantle and core is not well known. Here, we show experimental results of a reaction between Fe 3 C and H 2 O at pressures and temperatures of the deep mantle and core‐mantle boundary (CMB). The reaction produces diamond, FeO, and FeH x , suggesting that water can liberate carbon from the core in the form of diamond (“core carbon extraction”) while the core gains hydrogen, if subducted water reaches to the CMB. Therefore, Earth's deep water and carbon cycles can be linked. The extracted core carbon can explain a significant amount of the present‐day mantle carbon. Also, if diamond can be collected by mantle flow in the region, it can result in unusually high seismic‐velocity structures.

Ko, Byeongkwan↗
DOE OSTI2025

Solidification cracking of laser melted commercial-purity tungsten

The high melting temperature of tungsten (W) makes it an attractive candidate for energy generation applications; however, its use is limited by its poor ductility at low temperatures. This limitation affects even melt-based additive manufacturing (AM) processes such as laser powder bed fusion (PBF-LB), with as-fabricated pure W exhibiting both longitudinal and transverse cracks. Metallurgical and processing factors that affect these cracks are still being explored. This work utilizes powderless single-tracks on pure W plates made over a range of power and velocity combinations in three unique PBF-LB setups: one with flowing argon, one with continuous vacuum, and an in situ high-speed synchrotron X-ray radiography setup. Each had different oxygen activities in their build environments. Longitudinal cracks with oxides appearing to exude from the crack were found for tracks deposited in a build environment with argon shield gas, while no such cracks were observed for builds conducted in vacuum. A calculation showed that direct oxygen ingress into the melt pool from the build environment is negligible when argon shield gas is present. Thus, incorporation of spatter is a likely mechanism for oxygen ingress into the melt pool. Radiography showed extensive keyholing, spatter generation, and crack formation at keyhole porosity. A heat transfer calculation showed the crack formation time was consistent with cooling below its ductile-to-brittle transition (DBT) temperature. This work identifies solidification cracking as a feasible mechanism in pure W beyond the well-known DBT related cracking.

Materials science↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Study of Wake Characteristics of Marine Turbine Arrays: Preprint

As the marine energy industry continues to develop and installations grow from single turbines to small arrays and larger farms, understanding and predicting wake behavior using modeling tools is necessary for array design and optimization. The wake characteristics of upstream turbines, including velocity deficit, wake swirling, tip vortices, wake shape and direction, and flow recovery, influence the inflow conditions experienced by downstream turbines in an array and, subsequently, affect the individual turbine loading and power generated by the array. These wake characteristics are influenced by many factors, including site geography, inflow shear profile, turbulence, waves, the strength of the current resource, gravitational forces, storm surges, and seasonal effects like temperature and wind conditions. Efficient array modeling tools that capture these wake effects and resulting turbine-to-turbine interactions need to be developed to optimize farm layouts and design individual turbines for operation within arrays. Unlike the wind industry, the marine energy research community is largely limited to the use of either numerical tools or tank/flume testing to understand marine turbine wake interactions. This paper will investigate recent advances in numerical and experimental marine turbine wake modeling, with a focus on recent publications and on both axial-flow and cross-flow marine turbines. A summary of recent research on these topics is provided. Additionally, a demonstration of turbine wake interactions using computational fluid dynamics is presented.

marine hydrokinetic turbine↗
DOE OSTI2024

TEAMER - Extreme Events Modeling for the MARMOK-OWC Wave Energy Converter

Through the TEAMER program, Sandia National Laboratories (SNL) collaborated with IDOM Incorporated to study their MARMOK-Oscillating Water Column (MARMOK-OWC) wave energy conversion device. The study yielded a quantitative understanding of hydrodynamic pressures on the oscillating water column (OWC) device surfaces, the mooring tensions, and the dynamic performance of the device under extreme ocean wave conditions. This project utilized a comprehensive multi-phase Navier-Stokes flow solver with an overset body-fit mesh to predict fluid velocities and hydrodynamic forces on the MARMOK-OWC device. Computational Fluid Dynamics (CFD) analysis were conducted using OpenFOAM. This data includes the OpenFOAM cases (setup and data) to run the extreme events developed during the project. This project is part of the TEAMER RFTS 4 (request for technical support) program.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Hydrodynamics of a multicomponent vesicle under strong confinement

Herein we numerically investigate the hydrodynamics and membrane dynamics of a multicomponent vesicle in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement we find that lipid coarsening is arrested under strong confinement, and resumes at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement.

97 MATHEMATICS AND COMPUTING↗