Engineering Papers⌕ Search

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.

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

At least 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗
DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗
DOE OSTI2022

Analysis of contact conditions and microstructure evolution in shear assisted processing and extrusion using smoothed particle hydrodynamics method

Shear assisted processing and extrusion (ShAPE) is a solid-phase processing technique that adds an additional shear force as compared with a conventional extrusion approach. Recently, ShAPE has demonstrated the capability of extruding high-performance aluminum alloy 7075 (AA7075) tubes at speeds up to 12.2 m/min without surface tearing. However, the relationship among the ShAPE processing parameters, thermomechanical conditions, contact conditions, heat generation, and microstructure evolution remains primarily empirical because an insightful understanding of the associated physics is still lacking. To help elucidate these relationships, this work proposes a thermomechanical meshfree model for the first time for ShAPE processing of AA7075 using the smoothed particle hydrodynamics (SPH) method. The meshfree model is first validated thoroughly by experimental data in terms of material flow, die face temperature, and extrusion force with various processing parameters. The validated model is then used to analyze the steady-state contact conditions and heat generation rates during ShAPE processing. Distributions of the average grain size of AA7075 being extruded are calculated using the SPH model output. The meshfree model results reveal that extrusions conducted at lower temperatures and higher strain rates yield more refined grains and possibly higher material strength, which is also consistent with the experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Transport of Compressed Woody Biomass: Correlating Rheology and Microcompounder Measurements

Transport of milled solid biomass can be problematic in biomass conversion processes due to the material’s high cohesiveness, low bulk density, and poor flowability. Here, the goal is to investigate the flow behavior of compressed biomass under conditions similar to compression-screw feeders. To mimic the flow of woody biomass inside a compression-screw feeder, two lignocellulosic biomass feedstocks were fed to a co-rotating twin-screw microcompounder. Force and screw speed from the microcompounder were correlated to viscosity in a traditional rheometer, where low-density polyethylene served as a control. The deformation of compressed biomass in a lab-scale compounder is strongly dependent on screw speed and moisture content. Shear stress in woody biomass decreased with the shear rate, indicating negative plastic viscosity. Moreover, viscosity decreased one order of magnitude as screw speed increased from 10 to 90 rpm. Increasing moisture from 10 to 50% reduced the viscosity of the compressed biomass by 60%.

09 BIOMASS FUELS↗
DOE OSTI2021

Influences on Discharge Partitioning on a Large River Delta: Case Study of the Mississippi-Atchafalaya Diversion, 1926-1950

The modern Mississippi River (M.R.) Delta is plumbed by the Mississippi and Atchafalaya rivers, setting water and sediment dispersal pathways for Earth’s fifth-largest river system. The Atchafalaya River’s (A.R.) partial annexation of discharge from the M.R., particularly between 1926 and 1950, prompted warnings of a rapid river avulsion and the construction of the Old River Control Structure to regulate flow. Natural and anthropogenic causes of this discharge annexation are difficult to disentangle. Here, we develop and validate a hydrodynamic model of flow partitioning through the historic channel network. We then isolate how several key changes to the system affected discharge partitioning and stage at the A.R.-M.R. diversion. Simulations show that erosion of the upper A.R. can account for 73% of the water discharge increase. Dredging in the lower A.R. between 1932 and 1950 can account for 35% of the water discharge increase, and was also an important control on shear stress distribution. The lower M.R. was slightly erosional during this period, and therefore hindered the A.R. discharge increase. Significant lacustrine delta deposition in A.R. had little effect on partitioning. These findings highlight the importance of A.R. enlargement processes on avulsion dynamics at this site. Given the essential nature of this river junction to the society, transportation, and commerce of the United States, improved attribution of discharge increases may lead to future management strategies that are broadly impactful.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

PARTICLE RESOLVED DNS STUDY OF TURBULENCE EFFECTS ON HYPORHEIC MIXING IN RANDOMLY PACKED SEDIMENT BEDS

Pore-resolved direct numerical simulations (DNS) are used to investigate the interactions between stream-water flow turbulence and groundwater flow through a porous sediment bed in the hyporheic zone. Two permeability Reynolds numbers (2.56 and 5.17), representative of aquatic systems and representing ratio of permeability to viscous length scales, were simulated to understand its influence on the momentum exchange at the sediment-water interface (SWI). A doubleaveraging methodology is used to compute the Reynolds stresses, form-induced stresses, and pressure fluctuations. It is observed that both shear layer and turbulent shear stress penetration increases with ReK. Reynolds and form-induced bed-normal stresses increase with ReK. The peak values of the form-induced stresses for the lower (2.56) and higher (5.17) ReK happen within the top layer of the sediment bed. The sum of turbulent and form-induced pressure fluctuations, analyzed at their respective zero-displacement planes, are statistically similar and can be well approximated by a t location-scale distribution fit providing with a model that could potentially be used to impose boundary conditions at the SWI in reach scale simulations.

Karra, Shashank↗
DOE OSTI2023

Investigation of the yielding transition in concentrated colloidal systems via rheo-XPCS

We probe the microstructural yielding dynamics of a concentrated colloidal system by performing creep/recovery tests with simultaneous collection of coherent scattering data via X-ray Photon Correlation Spectroscopy (XPCS). This combination of rheology and scattering allows for time-resolved observations of the microstructural dynamics as yielding occurs, which can be linked back to the applied rheological deformation to form structure–property relations. Under sufficiently small applied creep stresses, examination of the correlation in the flow direction reveals that the scattering response recorrelates with its predeformed state, indicating nearly complete microstructural recovery, and the dynamics of the system under these conditions slows considerably. Conversely, larger creep stresses increase the speed of the dynamics under both applied creep and recovery. The data show a strong connection between the microstructural dynamics and the acquisition of unrecoverable strain. By comparing this relationship to that predicted from homogeneous, affine shearing, we find that the yielding transition in concentrated colloidal systems is highly heterogeneous on the microstructural level.

42 ENGINEERING↗
DOE OSTI2021

Learning unknown physics of non-Newtonian fluids

Here, we use physics-informed neural networks (PINNs) to learn viscosity models of two non-Newtonian systems (polymer melts and suspensions of particles) using only velocity measurements. For synthetic velocity data generated with the power-law viscosity model, the PINN-inferred viscosity model agrees with the analytical model for shear rates with large absolute values but deviates for shear rates near zero where the analytical model has an unphysical singularity. Once the viscosity model is learned the PINN method can solve the momentum conservation equation using only the boundary conditions.

42 ENGINEERING↗
DOE OSTI2024

Direct numerical simulation of a high-pressure hydrogen micromix combustor: Flame structure and stabilisation mechanism

A high-pressure hydrogen micromix combustor has been investigated using direct numerical simulation with detailed chemistry to examine the flame structure and stabilisation mechanism. The configuration of the combustor was based on the design by Schefer et al., using numerical periodicity to mimic a large square array. A precursor simulation of an opposed jet-in-crossflow was first conducted to generate appropriate partially-premixed inflow boundary conditions for the subsequent reacting simulation. The resulting flame can be described as an predominantly-lean inhomogeneously-premixed lifted jet flame. Five main zones were identified: a jet mixing region, a core flame, a peripheral flame, a recirculation zone, and combustion products. The core flame, situated over the jet mixing region, was found to burn as a thin reaction front, responsible for over 85% of the total fuel consumption. The peripheral flame shrouded the core flame, had low mean flow with high turbulence, and burned at very lean conditions (in the distributed burning regime). It was shown that turbulent premixed flame propagation was an order-of-magnitude too slow to stabilise the flame at these conditions. Stabilisation was identified to be due to ignition events resulting from turbulent mixing of fuel from the jet into mean recirculation of very lean hot products. Ignition events were found to correlate with shear-driven Kelvin-Helmholtz vortices, and increased in likelihood with streamwise distance. At the flame base, isolated events were observed, which developed into rapidly burning flame kernels that were blown downstream. Further downstream, near-simultaneous spatially-distributed ignition events were observed, which appeared more like ignition sheets. The paper concludes with a broader discussion that considers generalising from the conditions considered here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Investigations of the linear and non-linear flow harmonics using the a multi-phase transport model

The multi-phase transport model (AMPT) is used to study the effects of the parton-scattering cross-sections (σ pp ) and hadronic re-scattering on the linear contributions to the flow harmonic v 4 , the non-linear response coefficients, and the correlations between different order flow symmetry planes in Au + Au collisions at 200 GeV. The model results, which agree with current experimental measurements, indicate that the higher-order flow harmonics are sensitive to the σ pp variations. However, the non-linear response coefficients and the correlations between different order flow symmetry planes are σ pp independent. Furthermore, these results suggest that further detailed experimental measurements which span a broad range of collision systems and beam energies could serve as an additional constraint for the theoretical models' calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Modified equilibrium distributions for Cooper-Frye particlization

We introduce a positive-definite single-particle distribution that is suitable for describing the transition from a macroscopic hydrodynamic to a microscopic kinetic description during the late stages of heavy-ion collisions in the presence of moderately large viscous corrections. The modified equilibrium distribution function can be constructed with hydrodynamic input from either relativistic viscous fluid dynamics or anisotropic fluid dynamics. We test the modified equilibrium distribution's hydrodynamic output for a stationary hadron resonance gas subject to either shear stress, bulk pressure, or baryon diffusion current at a given freeze-out temperature and baryon chemical potential. While it does not reproduce all components of the net baryon current and energy-momentum tensor exactly, it significantly improves upon the customary linearized approximations for the nonequilibrium correction δf n which typically lead to unphysical negative distribution functions at large particle momenta. A comparison of particle spectra and p T –differential elliptic flow coefficients from the Cooper–Frye formula computed with the modified equilibrium distribution and with linearized δf n corrections is presented, for two different (2 + 1)–dimensional hypersurfaces corresponding to central and noncentral Pb + Pb collisions at the Large Hadron Collider (LHC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Burst Pressure Solutions of Thin and Thick-Walled Cylindrical Vessels

Pressure vessels (PVs) are widely used in the energy industry. Accurate burst pressure is critical to structural design and safe operation for both thin and thick-walled PVs. The traditional strength theories utilized a single-parameter material property, such as the yield stress or the ultimate tensile stress (UTS) to develop failure models for determining the yield or ultimate pressure carrying capacity in the PV design. The UTS-based Barlow formula is a typical burst pressure model developed from the Tresca strength theory that provides the basis for developing regulation rules and failure models for different industry design codes, such as ASME BPVC, ASME B31.3, and ASME B31G, among others. In order to reduce the conservatism of the Tresca strength model, ASME BPVC recently adapted failure models developed from the von Mises strength theory for the PV design and analysis. It has been commonly accepted that the burst pressure of pipelines depends on the UTS and strain hardening exponent, n, of the pipeline steel. An average shear stress yield theory was thus developed, and the Zhu-Leis solution of burst pressure was obtained as a function of UTS and n for thin-walled line pipes. Experiments showed that the Zhu-Leis solution provides an accurate, reliable prediction of burst pressure for defect-free thin-walled pipes. In order to extend the Zhu–Leis solution to thick-walled cylindrical PVs, this paper defined three new flow stresses, modified the traditional strength theories, and obtained three new burst pressure solutions that are valid for both thin and thick-walled cylindrical vessels. The proposed flow stresses are able to describe the tensile strength and the plastic flow response of PVs for a strain hardening steel. The associated strength theories were then developed in terms of the Tresca, von Mises, and Zhu-Leis yield criteria. From these new strength theories, three burst pressure solutions were obtained for thick-walled cylinders, where the von Mises solution is an upper bound prediction, the Tresca solution is a lower bound prediction, and the Zhu-Leis solution is an intermediate prediction of burst pressure for thick-walled cylinders. Lastly, the proposed burst pressure solutions were evaluated and validated by two large datasets of full-scale burst tests for thick-walled tubes and for thin-walled pipes.

42 ENGINEERING↗
DOE OSTI2023

Universal behavior in fragmenting brittle, isotropic solids across material properties

A bonded particle model is used to explore how variations in the material properties of brittle, isotropic solids affect critical behavior in fragmentation. To control material properties, a model is proposed which includes breakable two- and three-body particle interactions to calibrate elastic moduli and mode I and mode II fracture toughnesses. In the quasistatic limit, fragmentation leads to a power-law distribution of grain sizes which is truncated at a maximum grain mass that grows as a nontrivial power of system size. In the high-rate limit, truncation occurs at a mass that decreases as a power of increasing rate. A scaling description is used to characterize this behavior by collapsing the mean-square grain mass across rates and system sizes. Consistent scaling persists across all material properties studied, although there are differences in the evolution of grain size distributions with strain as the initial number of grains at fracture and their subsequent rate of production depend on Poisson's ratio. Importantly, this evolving granular structure is found to induce a unique rheology where the ratio of the shear stress to pressure, an internal friction coefficient, decays approximately as the logarithm of increasing strain rate. The stress ratio also decreases at all rates with increasing strain as fragmentation progresses and depends on elastic properties of the solid.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solid–liquid separation of lignocellulosic sugars from biomass by rotating ceramic disc filtration

In many biomass conversion processes, the separation of cellulosic sugars from residual, lignin-rich solids is a critical step, and achieving high recovery yields of sugars by conventional tangential crossflow and vacuum filtration is challenged by the presence of biomass solids, which rapidly foul filters, resulting in decreased throughput. Considering the performance limitations of traditional filtration methods, dynamic filtration, which generates high shear at the membrane surface to decrease fouling, is emerging as a viable alternative for demanding solid–liquid separations. For high solids separations, there is little available information regarding the performance, limitations, and energy consumption of dynamic filtration. To that end, here we characterized the performance of a dynamic filtration module, specifically a rotating ceramic disc (RCD) filter, for the aseptic recovery of cellulosic sugars from biomass solids following pretreatment and enzymatic hydrolysis. We show how RCD rotational velocity and percent biomass solids impact the filter throughput. Additionally, we used computational fluid dynamics (CFD) simulations to estimate the shear rate at the membrane surface and to visualize hydrodynamic profiles within the module. With the combined CFD simulations and experimental results, we estimated the energy demand and operating expenses for a viable dynamic filtration system operating with a lignocellulosic feed slurry. Our results indicate that an RCD filter can achieve ≥95% recovery of sugars and produce a retentate slurry containing 12 wt% insoluble solids with low energy consumption (a 2.2-fold improvement over cross-flow filtration) and low operating costs ($\$$0.06 per kg sugars). These results show a viable path towards operationally reliable, energy efficient, and cost-effective separations of sterilized cellulosic sugars from biomass solids and highlight the potential of dynamic filtration systems for challenging solid–liquid separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Computational Fluid Dynamics analysis of a jet-induced attrition

The flow dynamics, collisional behavior and attrition rates of a grid jet attrition unit are computationally investigated using a kinetic theory-based approach. Profiles of the flow statistical properties are presented and discussed to describe the jet dynamics. Attrition behavior and rates are analyzed from the local reconstruction of the collision energy spectra. The commonly used jet-induced attrition mechanisms which involve the entrainment of the particles into the jet followed by their acceleration and collision at the jet tip is reproduced by the CFD calculations in this work. Here, supported by the local attrition rates, it is shown that the shear-based collisions around the jets, due to counter-current particle axial motions, are an additional jet-induced mechanism generating more fines than at the jet tips.

42 ENGINEERING↗
DOE OSTI2022

Summertime Marine Boundary Layer Cloud, Thermodynamic, and Drizzle Morphology over the Eastern North Atlantic: A Four-Year Study

Abstract Summertime remote sensor and in situ data from 2016 to 2019 collected at the ARM Eastern North Atlantic (ENA) Observatory are combined with aircraft measurements from the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign to quantify marine boundary layer (MBL) cloud, thermodynamic, and drizzle morphology in the region. A radar reflectivity–rainfall rate relationship ( Z – R ) is developed from aircraft data and 6-h cloud morphological regimes are identified from ENA data using a k -means algorithm driven by three independent inputs quantifying cloud thickness, drizzle intensity, and cloud field geometric complexity. Four separate MBL structural regimes representing non- or weakly drizzling single-layer stratocumulus, drizzling stratocumulus and cumulus-coupled stratocumulus, deep convection, and broken clouds embedded in northerly flow are identified. Single-layer stratocumulus is indicated when weak subtropical anticyclones are significantly west of the ENA site, and the MBL is cooler and drier than when drizzling and cumulus-coupled stratocumulus and broken clouds are observed. Drizzling and cumulus-coupled stratocumulus clouds are observed on the eastern flank of strong subtropical anticyclones in deep warm moist air masses with wind speeds exceeding 7 m s −1 and strong near-surface wind shear. Broken clouds exhibit strong wind shear near the inversion, while single-layer stratocumulus clouds have lower wind speeds and minimal shear. Net latent heat fluxes in the subcloud layer resulting from a combination of the ocean surface heat flux and evaporating drizzle average near zero over long periods in drizzling and cumulus-coupled stratocumulus. The ECMWF reanalysis version 5 (ERA5) is found to accurately represent single-layer stratocumulus properties, while producing significant discrepancies when drizzling stratocumulus and cumulus-coupled stratocumulus are observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The EGS Collab Project: Status and Accomplishments

The EGS Collab project, supported by the US Department of Energy, is addressing challenges in implementing enhanced geothermal systems (EGS). This includes improving understanding of the stimulation of crystalline rock to create appropriate flow pathways, and the ability to effectively simulate both the stimulation and the flow and transport processes in the resulting fracture network. The project is performing intensively monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory. Data and observations from the field test are compared to simulations to understand processes and to build confidence in numerical modeling of the processes. In Experiment 1, we examined hydraulic fracturing an underground test bed at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, at a depth of approximately 1.5 km. We drilled eight sub-horizontal boreholes in a well-characterized phyllite. Six of the boreholes were instrumented with many sensor types to allow careful monitoring of stimulation events and flow tests, and the other two boreholes were used for water injection and production. We performed a number of stimulations and flow tests in the testbed. Our monitoring systems allowed detailed observations and collection of numerous data sets of processes occurring during stimulation and during dynamic flow tests. Long-term ambient temperature and chilled water flow tests were performed in addition to many tracer tests to examine system behavior. Data were rapidly analyzed, allowing adaptive control of the tests. Numerical simulation was used to answer key experimental design questions, to forecast fracture propagation trajectories and extents, and to analyze and evaluate results. Many simulations were performed in near-real-time in conjunction with the field experiments, with more detailed process study simulations performed on a longer timeframe. Experiment 2 will examine hydraulic shearing in a test bed being built at the SURF at a depth of about 1.25 km in amphibolite under a different set of stress and fracture conditions than Experiment 1. Five sets of fracture orientations were considered in design, and three orientations seem to be consistently observed.

EGS Collab↗