SEARCH · Engineering Papers

Results for “shear 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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Field Sensitivity Analysis of Turbulence Model Parameters for Flow Over a Wing

Reynolds-averaged-Navier-Stokes (RANS) turbulence models are a critical tool in computational-fluid-dynamics simulations of aerodynamic systems, but simulation results can be highly sensitive to RANS-model parameter choices. Sensitivity analysis can be used to quantify these impacts, and the objective of this study is to demonstrate field sensitivity analysis with respect to ten parameters in the 2003 Menter shear-stress-transport (SST) turbulence model. Here, the analysis is demonstrated for an application relevant to wind energy, namely, flow over a NACA 0015 wing at 12 degree angle of attack and a Reynolds number of 1.5 x 10 6 . We quantify sensitivity using Sobol indices and the mean-squared gradient, which are estimated using polynomial chaos and active subspace models, respectively. Our results indicate that there are substantial spatial variations in parameter sensitivities, with different sets of most-sensitive parameters near the wing, as well as in the downstream wake, consistent with the physical interpretations of the turbulence model inputs. We show that, for this particular turbulence model and flow, simultaneous dimension reduction is possible across all quantities of interest, enabling efficient exploration of model outcomes. Ultimately, this analysis provides new insights into turbulence model parameter sensitivities in incompressible flows, and also demonstrates the implementation of field sensitivity analysis for applications relevant to aerodynamics simulations.

17 WIND ENERGY↗

DOE OSTI2025

Novel Preprocessing Techniques for Enhancing Flowability of Miscanthus

Purpose-grown energy crops possess significant potential as promising feedstocks for biofuel and biochemical production. However, similar to other biomass feedstocks, their utilization has been hindered by handling and feeding challenges such as clogging and segregation, which prevent biorefineries from operating at full capacity. These challenges arise primarily due to particle attributes including irregular shape, high flexibility, and high compressibility, leading to interlocking, high compaction, and significant friction under stress consolidation. To enhance operational reliability, it is crucial to focus on preprocessing techniques to address these issues, alongside accurately characterizing the mechanical and physical properties and flow performance of these feedstocks. This study investigates novel preprocessing techniques, including pelletization and torrefaction, to fundamentally alter these particle attributes and improve flowability. Miscanthus, an herbaceous energy crop, was the primary focus of this study. Samples were processed using a Forest Concepts Crumbler rotary shear system and an oscillating multi-stack screen bed to achieve two nominal (baseline) particle sizes: 6mm and 2mm. The study evaluated the efficiency and energy consumption of size reduction, revealing that while 2mm particles required more energy, they exhibited better flowability and higher yield strength compared to 6mm particles. Shear and compression tests assessed the material's shear strength, internal friction, and elastic modulus, while wedge-shaped hopper tests measured the critical arching distance and mass flow rate. The results provide valuable insights into the mechanical properties and processing efficiency of Miscanthus, contributing to the development of a preprocessing framework that enhances energy efficiency and material flowability in bioenergy production. However, further preprocessing and development are required to optimize this framework fully.

09 - BIOMASS FUELS↗

DOE OSTI2021

Vortex Flow and Cavitation in Liquid Injection: A Comparison between High-Fidelity CFD Simulations and Experimental Visualizations on Transparent Nozzle Replicas

Experimental instantaneous shadowgraph visualizations on transparent glass nozzle replicas of high-pressure fuel injectors have been used to validate a novel in-house high-fidelity LES-VOF multiphase solver, to study the evolution of vortex flow and fuel cavitation. Both experiments and simulations capture the formation of an unsteady vapor structure inside the nozzle volume, which is referred to as 'string-cavitation'; strings are found at the core of the recirculation zones. Furthermore, high-fidelity simulations provide a very detailed insight into the vortex generation in the injector nozzle; strings appear within the time scales that are relevant for fast injection events (on the order of 0.1 milliseconds) and, for the problem under consideration, their generation seems mostly related to the flow pattern in the sac. It is also shown that vortexes interact, merge till they disrupt and favor the temporary inception of shear cavitation.

3-phase LES-VOF solver↗

DOE OSTI2022

Computational studies of hydraulic stressors for biological performance assessment in a hydropower plant with Kaplan turbine

We report Hydropower is currently one of the preeminent sources of renewable energy in the United States and globally. Hydropower plants also have detrimental impacts on the environment and ecology, including direct impacts to anadromous fish populations. The computational fluid dynamics (CFD) – based Biological Performance Assessment (BioPA) toolset is used for biological evaluations of fish passage through hydropower plants. The hydraulic environment of a hydropower plant was evaluated using CFD coupled with discrete element method (DEM) simulations. The predicted flow field and particle collision rate were validated against the experimental data in a water flume that has an idealized hydroturbine distributor geometry. Flow simulations were conducted to evaluate the hydraulic stressors, such as nadir pressure, fluid shear, runner collision, in a physical scale in a hydropower plant with Kaplan turbine which are responsible for injury and mortality of fish in a downstream migration. The cumulative exposure probability for the nadir pressure and collision with turbine runner was found to decrease with increased discharge rate. The lowest discharge rate shows the higher value of cumulative shear exposure probability. The maximum value of collision velocity increases with increased discharge rate. We offer the conclusions that will help in understanding various hydraulic stressors for biological assessment for environmentally sustainable hydroturbine passage.

13 HYDRO ENERGY↗

DOE OSTI2020

Low Mach number fluctuating hydrodynamics model for ionic liquids

We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively-studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.

42 ENGINEERING↗

DOE OSTI2025

TEAMER – Enhanced Flow Measurement for Aquantis Tidal Turbine Test

The AQ10 is a floating, two-bladed, passive yawing tidal turbine developed by Aquantis that has a 10-meterrotor diameter, 160 kW rating, and employs reliable off-the-shelf powertrain and power conversion hardware. Aquantis is planning on-water turbine power performance and loads (blade loading and thrust)testing, where the turbine will be pushed through still water up to 4 knots and placed in a ‘station keeping’ tow in a tidal race up to its rated speed. On-water testing will be conducted using vessels and floating platforms on the sea surface to improve ease of testing and reduce disturbance to the environment. In this TEAMER project, Pacific Northwest National Laboratory (PNNL) will conduct water velocity and turbulence measurements in front of the turbine during on-water testing using acoustic Doppler instrumentation. By measuring both the tidal current flowing past the turbine and the resulting electrical power output, test results will provide a power curve (power vs flow speed) for the turbine up to rated power. Measurements of turbulence and velocity shear in front of the rotor will also provide information to assess the structural response of the rotor blades. With this analysis, Aquantis can use the performance and loads data to validate Tidal Bladed and OpenFAST simulations of the measured operating conditions. Measuring the power performance of a prototype turbine is a valuable step to improving device development and conducting a complete power performance assessment to IEC/TS 62600-200 standards in the future.

16 TIDAL AND WAVE POWER↗

DOE OSTI2022

The influence of temperature on the strain-hardening behavior of Fe-22/25/28Mn-3Al-3Si TRIP/TWIP steels

In this work, the influence of temperature and stacking fault energy (SFE) on the strain-hardening behavior and critical resolved shear stress for twinning was investigated for three Fe–22/25/28Mn–3Al–3Si wt.% transformation- and twinning-induced plasticity (TRIP/TWIP) steels. The SFEs were calculated by two different methods, density functional theory and statistical thermodynamic modeling. The dislocation structure, observed at low levels of plastic deformation, transitions from “planar” to “wavy” dislocation glide with an increase in temperature, Mn content, and/or SFE. The change in dislocation glide mechanisms from planar to wavy reduces the strain hardening rate, in part due to fewer planar obstacles and greater cross slip activity. Additionally, the alloys exhibit a large decrease in strength and ductility with increasing temperature from 25 to 200 °C, attributed to a substantial reduction in the thermally activated component of the flow stress, predominate suppression of TRIP and TWIP, and a significant increase in the critical resolved shear stress for mechanical twinning. Interestingly, the increase in SFE with temperature had a rather minor influence on the critical resolved shear stress for mechanical twinning, and other temperature dependent factors which likely play a more dominant role are discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Prediction of Complex Observed Shear Wave Splitting Patterns at Ryukyu Subduction Zone Using a Strong Intra-Slab Anisotropy Model

Complex shear wave splitting (SWS) patterns in subduction zones are often interpreted geodynamically as resulting from complex mantle flow; however, this may not always be necessary. We analyzed 7,093 high-quality SWS measurements from teleseismic S waves recorded by Hi-net stations across the Ryukyu arc in Japan. Our findings show a systematic rotation of the fast S polarization from trench-parallel to trench-perpendicular depending on the earthquake backazimuth. For the same earthquake, the measured splitting patterns also vary spatially across the southwest Japan. Using full-wave seismic modeling, we showed that a dipping slab with ~30% shear anisotropy of the tilted transverse isotropy (TTI) type, with a symmetry axis perpendicular to the slab interface, can predict the observed delay times and polarization rotation. Our results highlight the importance of considering dipping anisotropic slabs in interpreting SWS at subduction zones.

58 GEOSCIENCES↗

DOE OSTI2025

Effect of nozzle translation on a bounded vortex flow

A bounded vortex flow is generated by a nozzle that combines azimuthally tilted downward jets in a circular array and a central suction port. When the nozzle is directed toward a bottom surface (called the impingement surface), the flow develops a strong intake vortex spanning from the impingement surface to the suction port, which acts to generate strong shear stress on the impingement surface. This strong shear stress can aid in removing particles attached to the impingement surface or in enhancement of heat or mass transfer, where in typical applications, the nozzle would be translated along the impingement surface. This paper reports on a computational study of the effect of nozzle translation on the bounded vortex flow field. The nozzle translation induces a bending of the bounded vortex in the downstream flow direction, which becomes larger with increase in the nozzle translation speed. We identify an interesting oscillatory phenomenon by which the bounded vortex can resist the bending induced by the translational flow and rotate about a mean position downstream of the suction port. At sufficiently high translation speeds, the oscillatory process can no longer resist the downwash effect from nozzle translation and the vortex eventually detaches from the impingement surface, limiting the effectiveness of this flow device.

42 ENGINEERING↗

DOE OSTI2020

Threading–Unthreading Transition of Linear-Ring Polymer Blends in Extensional Flow

Adding small amounts of ring polymers to a matrix of their linear counterparts is known to increase the zero-shear-rate viscosity because of linear-ring threading. Uniaxial extensional rheology measurements show that, unlike its pure linear and ring constituents, the blend exhibits an overshoot in the stress growth coefficient. By combining these measurements with ex-situ small-angle neutron scattering and nonequilibrium molecular dynamics simulations, this overshoot is shown here to be driven by a transient threading–unthreading transition of rings embedded within the linear entanglement network. Prior to unthreading, embedded rings deform affinely with the linear entanglement network and produce a measurably stronger elongation of the linear chains in the blend compared to the pure linear melt. Thus, rings uniquely alter the mechanisms of transient elongation in linear polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Assessment of Machine Learning Wall Modeling Approaches for Large Eddy Simulation of Gas Turbine Film Cooling Flows: An a Priori Study

Here, in this work, a priori analysis of machine learning (ML) strategies is carried out with the goal of data-driven wall modeling for large eddy simulation (LES) of gas turbine film cooling flows. High-fidelity flow datasets are extracted from wall-resolved LES (WRLES) of flow over a flat plate interacting with the coolant flow supplied by a single row of 7-7-7 shaped cooling holes inclined at 30 degrees with the flat plate at different blowing ratios (BR). The WRLES are performed using the high-order Nek5000 spectral element computational fluid dynamics (CFD) solver. Light gradient boosting machine (LightGBM) is employed as the ML algorithm for the data-driven wall model. Parametric tests are conducted to systematically assess the influence of a wide range of input flow features (velocity components, velocity gradients, pressure gradients, and fluid properties) on the accuracy of ML wall model with respect to prediction of wall shear stress. In addition, the use of spatial stencil and time delay is also explored within the ML wall modeling framework. It is shown that features associated with gradients of the streamwise and spanwise velocity components have a major impact on the prediction fidelity of wall model, while the effect of gradients of wall-normal velocity component is found to be negligible. Moreover, adding flow feature information from an x-y-z spatial stencil significantly improves the ML model accuracy and generalizability compared to just using local flow features from the matching location. Overall, highest prediction accuracy is achieved when both spatial stencil and time delay features are incorporated within the data-driven wall modeling paradigm.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2023

Predicted viscosity of liquid HMX up to 40 GPa

Viscous flow serves as a significant heating mechanism during the formation of hot spots, but the shear viscosity which determines this response is poorly characterized for most high explosives. Most recently, a model was proposed for the shear viscosity of liquid HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) that was fit to pressures reaching 5 GPa, but this work uncovered uncertainties in the viscosity at 0 GPa and remains untested at higher pressures. We use molecular dynamics (MD) simulations and the Green-Kubo formalism to predict the temperature- and pressure-dependent shear viscosity of HMX over the pressure interval 0 GPa≤P≤40 GPa. Reassessment of the viscosity at 0 GPa rules out several potential explanations for discrepancies between earlier reports; we attribute these differences to details of MD trajectory integration impacting molecular flexibility. The shear viscosity of HMX exhibits an Arrhenius temperature dependence at each pressure considered, with exponential prefactor and activation energy terms that are also strong functions of pressure. An analytic form for the viscosity is developed based on an extension of the well-known Andrade equation that simultaneously captures the temperature and pressure dependencies in the MD data up to 40 GPa. Comparison against a recently developed model for the viscosity of liquid RDX (1,3,5-trinitro-1,3,5-triazinane) shows that both materials exhibit similar functional dependencies with the viscosity of HMX being higher by roughly an order of magnitude at a given temperature-pressure state.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Toward Optimizing Electrospun Nanofiber Fuel Cell Catalyst Layers: Polymer–Particle Interactions and Spinnability

We investigate the effect of the poly(acrylic acid) (PAA) carrier polymer concentration on the microstructure and rheological properties of catalyst inks for electrospun polymer–electrolyte membrane fuel-cell catalyst layers. Characterization of an ink microstructure using oscillatory shear rheology showed that the catalyst particles (platinum on carbon) are significantly agglomerated in the absence of PAA or an ionomer. Both the ionomer and PAA promoted the stability of the particles against agglomeration via electrosteric stabilization by adsorbing onto the particle surface. Increasing the PAA concentration increased the stability of the particles (or reduced the agglomerated structure) due to increasing PAA coverage onto the free surface area of the particles. However, beyond a certain increase in concentration, PAA was found to predominantly remain as an excess free polymer in the ink due to an insufficient free/available surface area on the particles for further PAA coverage. Extensional rheology measurements demonstrated that PAA enhances the extensional viscosities of the inks. Consequently, increasing the PAA concentration in the ink promoted the evolution of uniform nanofibers. However, beyond a certain concentration, a significant increase in the shear viscosities of the inks led to defective fiber morphologies because of the onset of flow instabilities. Electrochemical performance comparisons between catalyst layers with different PAA concentrations showed maximum performance at the PAA concentration that led to the least agglomerated structure of the catalyst, most uniform fiber morphologies, and low concentrations of free (non-adsorbing) PAA in the electrode. These results provide a rationale for optimization of electrospun catalyst nanofibers for both spinnability and electrochemical performance.

25 ENERGY STORAGE↗

DOE OSTI2021

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗

DOE OSTI2020

Drag viscosity of metals and its connection to Coulomb drag

Recent years have seen a surge of interest in studies of hydrodynamic transport in electronic systems. We investigate the electron viscosity of metals and find a component that is closely related to Coulomb drag. By using linear-response theory, viscosity, which is a transport coefficient for momentum, can be extracted from the retarded correlation function of the momentum flux, i.e., the stress tensor. There exists a previously overlooked contribution to the shear viscosity from the interacting part of the stress tensor which accounts for the momentum flow induced by interactions. This contribution, which we dub drag viscosity, is caused by the frictional drag force due to long-range interactions. It is therefore linked to Coulomb drag which also originates from the interaction-induced drag force. Starting from the Kubo formula and using the Keldysh technique, we compute the drag viscosity of two- and three-dimensional metals along with the drag resistivity of double-layer two-dimensional electronic systems. Both the drag resistivity and drag viscosity exhibit a crossover from quadratic-in-$\textit{T}$ behavior at low temperatures to a linear behavior at higher temperatures. Although the drag viscosity appears relatively small compared with the normal Drude component for the clean metals, it may dominate hydrodynamic transport in some systems, which are discussed in the conclusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Physical properties, internal structure, and the three‐dimensional petrography of CI chondrites

physical properties and the nature of their breccation, we investigated nine samples of the Ivuna and Orgueil CI chondrites ranging in size from 1 mm to 4 cm in approximate diameter. The combined mass of unique material investigated in this work is 113 g. For our investigations, we use ideal gas pycnometry, 3-D laser scanning, x-ray computed microtomography (μCT), and accompanying digital data extraction techniques. We found that the bulk density of the samples ranged from 1.61 to 2.10 g cm −3 . Larger samples tend to have a lower bulk density. Grain density (ranging from 2.44 to 2.55 g cm −3 ) is significantly less variable than the bulk density in our samples and the quantity of porosity (ranging from 14.6% to 33.8%) is the dominant factor in determining the bulk density of CI chondrite material. Our μCT results show that the visible porosity across all sizes of our CI chondrite samples is in the form of cracks, but these cracks can account for less than two-thirds of the porosity in the CI chondrites. Other porosity is not visible, even at μCT resolutions of 2.7 μm voxel edge −1 and we conclude that it is sub-micron in nature. It is not clear if the cracks seen in our samples are indigenous to the chondrites or are a result of terrestrial processes. We also find that the CI chondrites are excellent examples of the fractal-like nature of brecciation, where clasts can be observed at all scales we imaged. The breccias are composed of sub-equant-shaped and sub-rounded-textured clasts like melt-free impact breccias on other solar system bodies. From our μCT volume and digital data extraction, we determine that the Ivuna CI chondrite breccia is organized: the mostly sub-equant clasts within our ~2 cm chunk of Ivuna have a mean diameter of 1.33 mm and their aligned longest axes define a lineation structure. We speculate that the lineation was imparted after fragmentation of the clasts by slight shear on the parent asteroid which could be the result of seismic-related granular flow or mild non-axial impact-related compaction. These data will help to place returned asteroidal material from asteroids 162173 Ryugu and 101955 Bennu and the CI chondrites into a mutual geological context.

CI chondrite↗

DOE OSTI2021

2019 Groundwater Monitoring Report Project Shoal Area: Subsurface Corrective Action Unit 447

The Project Shoal Area is a site in Nevada where an underground nuclear test was conducted in 1963. It later came to be known as the Shoal, Nevada, Site. Surface contamination at the site has been remediated, but investigation of groundwater contamination resulting from the test is still in the corrective action process. Annual sampling and water-level monitoring are conducted as part of the subsurface corrective action strategy, which has focused on revising the site conceptual model and evaluating the adequacy of the monitoring well network. It has also included enhancements to the monitoring well network to address uncertainties in the groundwater flow direction and the cause of rising water levels in site wells west of the shear zone since the first hydrologic characterization (HC) wells were installed in 1996. Revisions to the site conceptual model and enhancements to the monitoring strategy were provided to Nevada Division of Environmental Protection (NDEP) in the Addendum to: Corrective Action Decision Document/Corrective Action Plan (CADD/CAP) for the Subsurface Corrective Action Unit 447 Shoal, Nevada, Site. NDEP approved the addendum to the CADD/CAP, which included an expanded contaminant boundary and compliance boundary for the site.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Grassmannian Shape Representations for Aerodynamic Applications: Preprint

Airfoil shape design is a classical problem in engineering, science, and manufacturing. Our motivation is to combine principled physics-based considerations for the shape design problem with modern computational techniques informed by a data-driven approach. Traditional analyses of airfoil shapes emphasize a flow-based sensitivity to deformations which can be represented generally by affine transformations (rotation, scaling, shearing, shifting). We present a novel representation of shapes which decouples affine-style deformations from a rich set of data-driven deformations over a submanifold of the Grassmannian. The Grassmannian representation, informed by a database of physically relevant airfoils, offers (i) a rich set of novel 2D airfoil deformations not previously captured in the data, (ii) improved low-dimensional parameter domain for inferential statistics, and (iii) consistent 3D blade representation and perturbation over a sequence of nominal shapes.

blade representation↗