Turbulent flow of an elasticoviscous non- Newtonian fluid.
Friction factors and velocity profiles in turbulent flow of viscoelastic non-Newtonian fluid, noting correlation between frictional characteristics and Reynolds number
SEARCH · Search NASA
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.
Friction factors and velocity profiles in turbulent flow of viscoelastic non-Newtonian fluid, noting correlation between frictional characteristics and Reynolds number
This study investigates the extrusion processes of deuterium and protium using ANSYS-Polyflow. The geometries and computational fluid dynamics (CFD) settings closely replicate the experimental setups and data acquired from the extruder experiments at Oak Ridge National Laboratory (ORNL) for validation purposes. We explore the impacts of (1) slip versus non-slip boundary conditions and (2) the use of constant, temperature-, and shear rate–dependent viscosities, concluding that the implementation of non-slip wall boundary conditions combined with shear rate–dependent viscosity produced more accurate predictions. The simulations achieved excellent agreement with the experimental data, with relative differences of only 5% for deuterium, and 3% to 6% for protium. This is the first time that experimental extrusion data at ORNL have been accurately predicted through high-fidelity CFD modeling. In conclusion, the advancements offer valuable insights and a foundational modeling tool for optimizing pellet injectors for ITER and other future reactor-scale devices.
Solvent purity effects on non-Newtonian power law viscosity
Turbulent flow of dilute viscoelastic non- Newtonian fluids in pipes, noting frictional characteristics, elastic fluid parameters, etc
An experimental technique which allows visual observations of the wall area in turbulent pipe flow is described in detail. It requires neither the introduction of any injection or measuring device into the flow nor the presence of a two-phase flow or of a non-Newtonian fluid. The technique involves suspending solid MgO particles of colloidal size in trichloroethylene and photographing their motions near the wall with a high speed movie camera moving with the flow. Trichloroethylene was chosen in order to eliminate the index of refraction problem in a curved wall. Evaluation of the technique including a discussion of limitations is included. Also the technique is compared with previous methods of visual observations of turbulent flow.
A study was made of the behavior of polyisobutylene under motion at a constant stretch history for both strip biaxial extensional flow and simple extensional flow. Steady-state non-Newtonian viscosities were observed at various constant stretch histories. Newtonian viscosities for both strip biaxial and simple extensional flow were found to be in agreement with the classical theory. The results of the study provide an essential part of the experimental background necessary for the development of a new general stress-strain-time relation for uncrosslinked and lightly crosslinked polymers.
Simultaneous effects of temperature and concentration on tensile viscosity are investigated using a semi-empirical scaling law established to relate the concentration, temperature and deformation rate by a reduced deformation rate variable. Results show the onset of high viscosity occurring in the vicinity of the onset of non-Newtonian viscosity observed in a simple shear flow experiment. High viscosity is also compared with steady and dynamic shear data. Finally, the asymptotic solution of Hassager's kinetic theory of a bead-rod model for macromolecular solutions is compared with the data, indicating that tensile viscosity depends strongly on the molecular weight, the higher the molecular weight, the more sensitive the tensile viscosity in terms of the deformation rate.
A kinetic rate theory previously presented for describing non-Newtonian phenomena has been further modified to predict the flow behavior of viscoelastic materials under constant stress conditions. The thixotropic shear stress or shear rate is predicted by the kinetic theory, and the experimental stress or shear rate is obtained by modifying the thixotropic value by a stress or shear rate retardation term. The retardation term stems from a Maxwellian approach for stress retardation. In order to test the validity of this approach, transient and steady-state data were obtained for two solutions of polymethylmethacrylate in diethylphthalate. Both constant stress measurements and constant shear rate data were taken over a broad range.
Analysis predicts decrease in fluid viscosity by comparing pressure profile of fluid/gas mixture with that of power-law fluid. Fluid is taken to be viscous, non-Newtonian, and incompressible; the gas to be ideal; the flow to be inertia-free, isothermal, and one dimensional. Analysis assists in design of flow systems for petroleum, coal, polymers, and other materials.
In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.
Tutorials for modeling of slot-die and slide-die coating flows with Goma 7, an open source finite element code, are presented. The tutorials cover the workflow to attaining steady state solutions for these flows, and continuation strategies for navigating the operating windows. Advanced topics of coating window prediction, automated multiparameter continuation, non-Newtonian rheology, dynamic contact line modeling, and some more solution strategies are also covered.
A finite element computer code for the analysis of mantle convection is described. The coupled equations for creeping viscous flow and heat transfer can be solved for either a transient analysis or steady-state analysis. For transient analyses, either a control volume or a control mass approach can be used. Non-Newtonian fluids with viscosities which have thermal and spacial dependencies can be easily incorporated. All material parameters may be written as function statements by the user or simply specified as constants. A wide range of boundary conditions, both for the thermal analysis and the viscous flow analysis can be specified. For steady-state analyses, elastic strain rates can be included. Although this manual was specifically written for users interested in mantle convection, the code is equally well suited for analysis in a number of other areas including metal forming, glacial flows, and creep of rock and soil.
A procedure is outlined for the numerical solution of the complete elastohydrodynamic lubrication of rectangular contacts incorporating a non-Newtonian fluid model. The approach uses a Newtonian model as long as the shear stress is less than a limiting shear stress. If the shear stress exceeds the limiting value, the shear stress is set equal to the limiting value. The numerical solution requires the coupled solution of the pressure, film shape, and fluid rheology equations from the inlet to the outlet. Isothermal and no-side-leakage assumptions were imposed in the analysis. The influence of dimensionless speed, load, materials, and sliding velocity and limiting-shear-strength proportionality constant on dimensionless minimum film thickness was investigated. Fourteen cases were used in obtaining the minimum-film-thickness equation for an elastohydrodynamically lubricated rectangular contact incorporating a non-Newtonian fluid model. Computer plots are also presented that indicate in detail pressure distribution, film shape, shear stress at the surfaces, and flow throughout the conjunction.
The hydrodynamic stability of plane Poiseuille flow of second-order fluids to finite amplitude disturbances is examined using the method of Stuart and Watson as extended by Reynolds and Potter. For slightly non-Newtonian fluids subcritical instabilities are predicted. No supercritical equilibrium states are expected if the entire spectrum of disturbance wavelengths is present. Possible implications with respect to the Toms phenomenon are discussed.
Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .
The basic aspect of parallel-surface lubrication that distinguishes it from other areas of lubrication technology is that classical lubrication theory does not predict the existence of a stable hydrodynamic film for steady-state, isothermal, incompressible flow between smooth, parallel surfaces. Hydrodynamic films between apparently parallel surfaces have been observed in practice and are often essential for the reliable performance of thrust bearings and seals. In order to account for this fortunate discrepancy between classical theory and experiment, load-support mechanisms relating to accidential features characteristic of seal performance, which relax one or more of the assumptions in the classical theory and permit the theoretical prediction of load support are proposed. Some of the features that have been analyzed are vibratory effects such as wobble and bounce, surface waviness, nonsymmetric rotation resulting from various types of misalignment, lubricant density change, non-Newtonian lubricant effects, and surface roughness.
Stress is placed on the temperature dependence of both a linear Newtonian rheology and a nonlinear olivine rheology in accounting for narrow mantle flow structures. The boundary-layer theory developed incorporates an arbitrary temperature-dependent power-law rheology for the medium, in order to facilitate the study of mantle plume dynamics under real conditions. Thermal, kinematic, and dynamic structures of mantle plumes are modelled by a two-dimensional natural-convection boundary layer rising in a fluid with a temperature-dependent power-law relationship between shear stress and strain rate. An analytic similarity solution is arrived at for upwelling adjacent to a vertical isothermal stress-free plane. Newtonian creep as a deformation mechanism, thermal anomalies resulting from chemical heterogeneity, the behavior of plumes in non-Newtonian (olivine) mantles, and differences in the dynamics of wet and dry olivine are discussed.
We use nonequilibrium molecular dynamics simulations to verify recent tube-model predictions that associative polymer networks exhibit broad stretch fluctuations during elongational flow. Simulations further show that these fluctuating dynamics give rise to the rate-dependent extensional viscosity 𝜂 𝐸 measured in filament stretching experiments on H-bonding networks. Simulations model bivalent associative networks with a reactive bead-spring model for varying association strength and extensional strain rate. We observe that stretch fluctuations are driven by a new form of chain tumbling, where chains continually collapse and elongate as their associations break and reform within the advecting network. This produces a broad, nearly uniform distribution of chain stretch over a wide range of strain rates, manifesting as a rate-independent plateau in the extensional stress. Our results show that the nonlinear viscoelasticity of associative networks is dominated by large fluctuations in molecular response, which cannot be captured by current mean-field models.