Search NASA⌕ Search

SEARCH · Search NASA

Results for “Very viscous flows”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Extracting the speed of sound in quark–gluon plasma with ultrarelativistic lead–lead collisions at the LHC

Abstract Ultrarelativistic nuclear collisions create a strongly interacting state of hot and dense quark–gluon matter that exhibits a remarkable collective flow behavior with minimal viscous dissipation. To gain deeper insights into its intrinsic nature and fundamental degrees of freedom, we determine the speed of sound in an extended volume of quark–gluon plasma using lead–lead (PbPb) collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data were recorded by the CMS experiment at the CERN LHC and correspond to an integrated luminosity of 0.607 nb −1 . The measurement is performed by studying the multiplicity dependence of the average transverse momentum of charged particles emitted in head-on PbPb collisions. Our findings reveal that the speed of sound in this matter is nearly half the speed of light, with a squared value of 0.241 ± 0.002 ( stat ) ± 0.016 ( syst ) in natural units. The effective medium temperature, estimated using the mean transverse momentum, is 219 ± 8 ( syst ) MeV . The measured squared speed of sound at this temperature aligns precisely with predictions from lattice quantum chromodynamic (QCD) calculations. This result provides a stringent constraint on the equation of state of the created medium and direct evidence for a deconfined QCD phase being attained in relativistic nuclear collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Imaging the breaking of electrostatic dams in graphene for ballistic and viscous fluids

The charge carriers in a material can, under special circumstances, behave as a viscous fluid. Here, in this work, we investigated such behavior by using scanning tunneling potentiometry to probe the nanometer-scale flow of electron fluids in graphene as they pass through channels defined by smooth and tunable in-plane p-n junction barriers. We observed that as the sample temperature and channel widths are increased, the electron fluid flow undergoes a Knudsen-to-Gurzhi transition from the ballistic to the viscous regime characterized by a channel conductance that exceeds the ballistic limit, as well as suppressed charge accumulation against the barriers. Our results are well modeled by finite element simulations of two-dimensional viscous current flow, and they illustrate how Fermi liquid flow evolves with carrier density, channel width, and temperature.

Science & Technology - Other Topics↗
DOE OSTI2023

An Efficient Three-Dimensional CFD-Based Numerical Wave Tank for a Wave Energy Converter in Extreme Irregular Waves

A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high-fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not predicted by potential flow engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.

CFD↗
DOE OSTI2025

Development of fiber orientation in injection molding: Comparison of glass fiber, carbon fiber and their hybrid composites

Fiber orientation distribution (FOD) in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for glass fiber (GF), carbon fiber (CF), and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varied with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design and Production for Maximum Structural Efficiency With Respect to Fiber Orientation With Increased Understanding of Hybrid Fiber Flow Behavior

Discontinuous fiber-reinforced thermoplastic composites have gained considerable attention in automotive, aerospace, and other industries, due to their high-rate of production combined with their ability to attain complex and intricate shapes. Among other high-rate thermoplastic manufacturing processes, injection-molding is one of the most common manufacturing methods due to fast production and high surface finishing of complex geometries. Fiber orientation in discontinuous fiber composites plays a pivotal role in determining the mechanical, electrical, and thermomechanical properties, underscoring the necessity to comprehend fiber orientation in injection molded parts. Among different fiber types, glass and carbon fibers are most common in the composite industries. The recent trend of hybrid composites comprising both glass fiber (GF) and carbon fiber (CF) is also gaining importance in the automotive industry. Hybrid fiber options allow designers to optimize the balance between glass and carbon fibers by leveraging the high durability and low cost of GF while the strength and lightweight properties of CF. Consequently, comparing the fiber oriented distribution (FOD) of injection molded composites containing GF, CF, and a hybrid of GF/CF is critical to investigating the local mechanical properties of intricate structures for high-end applications. In Phase I of this project, FOD in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for GF, CF, and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varies with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Temperature Effects on Droplet Oscillation Decay with Application to Fuel Property Measurement

Observation of oscillation decay in droplets has been shown to be an effective approach in determining physical properties such as viscosity and surface tension for emerging biofuels using μl quantities. Herein this work extends the approach to higher temperatures relevant to fuel injection conditions for internal combustion engines. Experiments are conducted that use high-resolution strobed imaging of moving, heated, μm-sized, fuel droplets to capture shape oscillation decay through image processing and analysis. Two fuels are investigated, iso-butanol and a primary reference fuel (PRF 84), which is a mixture of iso-octane and n-heptane. A piezoelectric droplet generator is used to generate a continuous train of single droplets, which are given an initial perturbation and subsequently undergo damped oscillations, captured using strobed imaging. Surface tension and viscosity calculated at four different temperatures ranging from 30°C-55°C using the frequency and decay time associated with the fundamental mode are found to be within ~ 10% of reference values obtained from the literature. Complementary numerical simulations are performed that utilize a volume-of-fluid approach to track transient droplet oscillation phenomena along with heat and mass transfer in a 2D axisymmetric domain. Simulations, where droplet size and temperature can be independently varied, capture an expected decrease in oscillation frequency and increase in decay time, with increase in fuel temperature. The simulations are further used to study the relative contribution to deviations in surface tension and viscosity predictions due to heat and mass transfer from the droplet, as well as viscous effects violating the inviscid flow assumption in the droplet oscillation theory. Mass loss effects are found to be negligible. Temperature change due to heat transfer is found to have the next highest sensitivity, particularly for the more volatile PRF 84, which has a lower viscosity and associated Ohnesorge number. For isobutanol, which has a higher viscosity and Ohnesorge number, viscous effects contribute the most to deviation in fuel property predictions.

09 BIOMASS FUELS↗
DOE OSTI2021

Viscous Behavior of Clay‐Rich Rocks and Its Role in Focused Fluid Flow

Abstract Focused fluid flow is common in sedimentary basins worldwide, where flow structures often penetrate through sandy reservoir rocks, and clay‐rich caprocks. To better understand the mechanisms forming such structures, the impacts of the viscoelastic deformation and strongly nonlinear porosity‐dependent permeability of clay‐rich materials are assessed from an experimental and numerical modeling perspective. The experimental methods to measure the poroviscoelastic and transport properties of intact and remolded shale have been developed, and the experimental data is used to constrain the numerical simulations. It is demonstrated that viscoelastic deformation combined with nonlinear porosity‐dependent permeability triggers the development of localized flow channels, often imaged as seismic chimneys. The permeability inside a channel increases by several orders of magnitude compared to the background values. In addition, the propagation time scale and the channel size strongly depend on the material properties of the fluid and the rock. The time‐dependent behavior of the clay‐rich rock may play a key role in the long‐term integrity of the subsurface formations.

58 GEOSCIENCES↗
DOE OSTI2020

Effects of Atwood and Reynolds numbers on the evolution of buoyancy-driven homogeneous variable-density turbulence

The evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at Atwood numbers up to 0.75 and large Reynolds numbers is studied by using high-resolution direct numerical simulations. To help understand the highly non-equilibrium nature of buoyancy-driven HVDT, the flow evolution is divided into four different regimes based on the behaviour of turbulent kinetic energy derivatives. The results show that each regime has a unique type of dependence on both Atwood and Reynolds numbers. It is found that the local statistics of the flow based on the flow composition are more sensitive to Atwood and Reynolds numbers compared to those based on the entire flow. It is also observed that, at higher Atwood numbers, different flow features reach their asymptotic Reynolds-number behaviour at different times. The energy spectrum defined based on the Favre fluctuations momentum has less large-scale contamination from viscous effects for variable-density flows with constant properties, compared to other forms used previously. The evolution of the energy spectrum highlights distinct dynamical features of the four flow regimes. Thus, the slope of the energy spectrum at intermediate to large scales evolves from -7/3 to -1 , as a function of the production-to-dissipation ratio. The classical Kolmogorov spectrum emerges at intermediate to high scales at the highest Reynolds numbers examined, after the turbulence starts to decay. Finally, the similarities and differences between buoyancy-driven HVDT and the more conventional stationary turbulence are discussed and new strategies and tools for analysis are proposed.

42 ENGINEERING↗
DOE OSTI2026

Fragile-to-strong transition in liquid As 2 S 3 under pressure: The effect of melt metallization

The well-known classification of glass-forming melts into fragile and strong liquids has several notable exceptions, including water, silica, and certain phase-change materials (PCMs). These exceptional fluid systems exhibit a fragile-to-strong transition (FST) upon cooling: a transformation from a high-temperature liquid with fast atomic dynamics, low viscosity, and low flow activation energy, to a viscous supercooled melt with high energy barriers near the glass transition temperature T g . This behavior is critically important for non-volatile memories, photonic tensor cores, reconfigurable metamaterials, and other devices, that use PCMs, enabling nanosecond-scale crystallization in the fragile regime and long data retention in the strong regime near or below T g . A significant structural transformation is expected between these two viscosity regimes, along with a semiconductor-metal (SC-M) transition upon heating, driven by high internal pressure and associated density increase. By applying high external pressure to the canonical low-conducting chalcogenide melt As 2 S 3 , we observed both the FST and the SC-M transition, occurring simultaneously within the same domain of the P, T−phase space. These findings suggest that the FST is not limited to a few exceptional liquids but is a common phenomenon, at least in systems that exhibit melt metallization within specific regions of their P, T−phase diagrams.

first-principles molecular dynamics↗
DOE OSTI2024

Emergence of Microphysical Bulk Viscosity in Binary Neutron Star Postmerger Dynamics

Abstract In nuclear matter in isolated neutron stars, the flavor content (e.g., proton fraction) is subject to weak interactions, establishing flavor ( β -)equilibrium. However, there can be deviations from this equilibrium during the merger of two neutron stars. We study the resulting out-of-equilibrium dynamics during the collision by incorporating direct and modified Urca processes (in the neutrino-transparent regime) into general-relativistic hydrodynamics simulations with a simplified neutrino transport scheme. We demonstrate how weak-interaction-driven bulk viscosity in postmerger simulations can emerge and assess the bulk viscous dynamics of the resulting flow. We further place limits on the impact of the postmerger gravitational-wave strain. Our results show that weak-interaction-driven bulk viscosity can potentially lead to a phase shift of the postmerger gravitational-wave spectrum, although the effect is currently on the same level as the numerical errors of our simulation.

Astronomy & Astrophysics↗
DOE OSTI2023

Probing turbulence with lasers and coherent light sources

Fluid and plasma turbulence consists of multi-scale interactions where the energy of the large scales is dissipated through a cascade to the smallest scales set by the microphysics. Turbulence and mixing can be detrimental to system performance, such as in inertial confinement fusion (ICF). Our current understanding of turbulence is mostly derived from studies of incompressible flows, where the intermediate-scaled eddies in the “inertial range” bridge the transfer of energy from the large eddies to the dissipative range, and that there is no communication between eddies of disparate scales. Whether a similar cascade process exists for compressible turbulence, however, remains to be seen; experiments have been limited by diagnostics and supersonic flow facilities. This project has focused on the development of techniques to measure the evolution of flows and microphysical transport properties using particle tracking velocimetry. This project has investigated three main areas: (1) particle-tracking experiments with x-ray radiography, (2) shock-particle interactions, (3) quantifying viscous effects on particles in shocked flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Influence of Wetting on Viscous Fingering Via 2D Lattice Boltzmann Simulations

We present simulations of two-phase flow using the Rothman and Keller colour gradient Lattice Boltzmann method to study viscous fingering when a “red fluid” invades a porous model initially filled with a “blue” fluid with different viscosity. We conducted eleven suites of 81 numerical experiments totalling 891 simulations, where each suite had a different random realization of the porous model and spanned viscosity ratios in the range $$M\in [0.01,100]$$ M ϵ [ 0.01 , 100 ] and wetting angles in the range $$\theta _w\in [180^\circ ,0^\circ ]$$ θ w ϵ [ 180 ° , 0 ° ] to allow us to study the effect of these parameters on the fluid-displacement morphology and saturation at breakthrough (sweep). Although sweep often increased with wettability, this was not always so and the sweep phase space landscape, defined as the difference in saturation at a given wetting angle relative to saturation for the non-wetting case, had hills, ridges and valleys. At low viscosity ratios, flow at breakthrough is localized through narrow fingers that span the model. After breakthrough, the flow field continues to evolve and the saturation continues to increase albeit at a reduced rate, and eventually exceeds 90% for both non-wetting and wetting cases. The existence of a complicated sweep phase space at breakthrough, and continued post-breakthrough evolution suggests the hydrodynamics and sweep is a complicated function of wetting angle, viscosity ratio and time, which has major potential implications to Enhanced Oil Recovery by water flooding, and hence, on estimates of global oil reserves. Validation of these results via experiments is required to ensure they translate to field studies.

Engineering↗
DOE OSTI2024

Λ polarization from vortex rings as the medium response for jet thermalization

We performed a systematic study on the formation of vorticity rings as the process for jet thermalization in the medium created in high-energy nuclear collisions. In this work, we expanded our previous analysis to a more realistic framework by considering noncentral events and fluctuations in the initial condition. We simulate the formation and evolution of the flow vortex structure in a relativistic viscous hydrodynamic model and study the sensitivity of the proposed “ring observable” (ℛ$^{𝑡}_{Λ}$) that can be measured experimentally through the polarization of Λ hyperons. We show that this observable is robust with respect to fluctuating initial conditions to capture the jet-induced vortex flow signal and further study its dependence on different model parameters, such as the jet's velocity, position, the fluid's shear viscosity, and the collision centrality. The proposed observable is associated with the formation of vorticity in a quark-gluon plasma, showing that the measurement of particle polarization can be a powerful tool to probe different properties of jet-medium interactions and to understand better the polarization induced by the transverse and longitudinal expansions of the medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Experiment and model for a Stokes layer in a strongly coupled dusty plasma

A Stokes layer, which is a flow pattern that arises in a viscous fluid adjacent to an oscillatory boundary, was observed in an experiment using a two-dimensional strongly coupled dusty plasma. Liquid conditions were maintained using laser heating, while a separate laser manipulation applied an oscillatory shear that was localized and sinusoidal. The evolution of the resulting flow was analyzed using space-time diagrams. These figures provide an intuitive visualization of a Stokes layer, including features such as the depth of penetration and wavelength. Another feature, the characteristic speed for the penetration of the oscillatory flow, also appears prominently in space-time diagrams. To model the experiment, the Maxwell-fluid model of a Stokes layer was generalized to describe a two-phase liquid. In our experiment, the phases were gas and dust, where the dust cloud was viscoelastic due to strong Coulomb coupling. Furthermore, the model is found to agree with the experiment, in the appearance of the space-time diagrams, and in the values of the characteristic speed, depth of penetration, and wavelength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

History Matching and Performance Prediction of a Polymer Flood Pilot in Heavy Oil Reservoir on Alaska North Slope

The first-ever polymer flood pilot to enhance heavy oil recovery on Alaska North Slope (ANS) is ongoing. After more than 2.5 years of polymer injection, significant benefit has been observed from the decrease in water cut from 65% to less than 15% in the project producers. The primary objective of this study is to develop a robust history-matched reservoir simulation model capable of predicting future polymer flood performance. In this work, the reservoir simulation model has been developed based on the geological model and available reservoir and fluid data. In particular, four high transmissibility strips were introduced to connect the injector-producer well pairs, simulating short-circuiting flow behavior that can be explained by viscous fingering and reproducing the water cut history. The strip transmissibilities were manually tuned to improve the history matching results during the waterflooding and polymer flooding periods, respectively. It has been found that higher strip transmissibilities match the sharp water cut increase very well in the waterflooding period. Then the strip transmissibilities need to be reduced with time to match the significant water cut reduction. The viscous fingering effect in the reservoir during waterflooding and the restoration of injection conformance during polymer flooding have been effectively represented. Based on the validated simulation model, numerical simulation tests have been conducted to investigate the oil recovery performance under different development strategies, with consideration for sensitivity to polymer parameter uncertainties. The oil recovery factor with polymer flooding can reach about 39% in 30 years, twice as much as forecasted with continued waterflooding. Besides, the updated reservoir model has been successfully employed to forecast polymer utilization, a valuable parameter to evaluate the pilot test’s economic efficiency. All the investigated development strategies indicate polymer utilization lower than 3.5 lbs/bbl in 30 years, which is economically attractive.

Wang, Xindan↗
DOE OSTI2024

An extension of Thwaites’ method for turbulent boundary layers

Thwaites developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations. The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

17 WIND ENERGY↗