Search NASA⌕ Search

SEARCH · Search NASA

Results for “NON-NEWTONIAN 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 109 records · Page 6

Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling

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 × .

Close-clearance impellers↗

Scientific machine learning for modeling and simulating complex fluids

The formulation of rheological constitutive equations—models that relate internal stresses and deformations in complex fluids—is a critical step in the engineering of systems involving soft materials. While data-driven models provide accessible alternatives to expensive first-principles models and less accurate empirical models in many engineering disciplines, the development of similar models for complex fluids has lagged. The diversity of techniques for characterizing non-Newtonian fluid dynamics creates a challenge for classical machine learning approaches, which require uniformly structured training data. Consequently, early machine-learning based constitutive equations have not been portable between different deformation protocols or mechanical observables. Here, we present a data-driven framework that resolves such issues, allowing rheologists to construct learnable models that incorporate essential physical information, while remaining agnostic to details regarding particular experimental protocols or flow kinematics. These scientific machine learning models incorporate a universal approximator within a materially objective tensorial constitutive framework. By construction, these models respect physical constraints, such as frame-invariance and tensor symmetry, required by continuum mechanics. We demonstrate that this framework facilitates the rapid discovery of accurate constitutive equations from limited data and that the learned models may be used to describe more kinematically complex flows. This inherent flexibility admits the application of these “digital fluid twins” to a range of material systems and engineering problems. We illustrate this flexibility by deploying a trained model within a multidimensional computational fluid dynamics simulation—a task that is not achievable using any previously developed data-driven rheological equation of state.

Science & Technology - Other Topics↗

On the Response of a Herschel–Bulkley Fluid Due to a Moving Plate

In this paper, we study the boundary-layer flow of a Herschel–Bulkley fluid due to a moving plate; this problem has been experimentally investigated by others, where the fluid was assumed to be Carbopol, which has similar properties to cement. The computational fluid dynamics finite volume method from the open-source toolbox/library OpenFOAM is used on structured quad grids to solve the mass and the linear momentum conservation equations using the solver “overInterDyMFoam” customized with non-Newtonian viscosity libraries. The governing equations are solved numerically by using regularization methods in the context of the overset meshing technique. The results indicate that there is a good comparison between the experimental data and the simulations. The boundary layer thicknesses are predicted within the uncertainties of the measurements. The simulations indicate strong sensitivities to the rheological properties of the fluid.

36 MATERIALS SCIENCE↗

The role of a pressure-dependent rheology in the dynamics of mantle circulation

A thermomechanical model for upper mantle convection was constructed such that the thickness and the structure of the lithosphere are determined self consistently by the heat transported by convection. In this study of the interaction between the lithosphere and upper mantle, strongly temperature and pressure dependent rheologies for both Newtonian and non-Newtonian creep mechanisms are employed. For a strictly temperature dependent rheology an insignificant amount of heat, less than 12.5 mW/sq m, can be transported convectively for an interior viscosity, 0(10 sup 21 Pas), compatible with post glacial rebound. On the other hand, for similar values of the interior viscosity, steady heat fluxes between 20 and 40 mW/sq m are produced by introducing pressure dependence into the rheology. For the temperature and pressure dependent flow law the horizontally averaged interior temperature displays very little variation with the amount of heat evacuated, once all of the rheological parameters are fixed. This finding may have important ramifications for parameterized convection.

Yuen, D. A.↗

Shear-Sensitive Liquid Crystal Coating Method: Surface-Inclination Effects on Shear Vector Measurements

The shear-sensitive liquid crystal coating (SSLCC) method is an image-based technique for both visualizing dynamic surface-flow phenomena, such as transition and separation, and for measuring the continuous shear-stress vector distribution acting on an aerodynamic surface. Under proper lighting and viewing conditions (discussed below), the coating changes color in response to an applied aerodynamic shear. This color-change response is continuous and reversible, with a response time of milliseconds, and is a function of both the shear magnitude and the shear vector orientation relative to the observer. The liquid crystal phase of matter is a weakly-ordered, viscous, non-Newtonian fluid state that exists between the nonuniform liquid phase and the ordered solid phase of certain organic compounds. Cholesteric liquid crystal compounds possess a helical molecular arrangement that selectively scatters white light, incident along the helical axis, as a three-dimensional spectrum. This property is linked to the helical pitch length, which is within the range of wavelengths in the visible spectrum. The pitch length, and hence the wavelength of the scattered light, is influenced by shear stress normal to the helical axis. This unique optical property produces a measurable color change in response to an applied shearing force. The full-surface shear stress vector measurement method, developed at NASA-Ames, is schematically illustrated. As with the visualization method, the coated test surface is illuminated from the normal direction with white light and the camera is positioned at an above-plane view angle of approximately 30 deg. Experiments have been initiated at NASA Ames to begin the process of quantifying surface-inclination (surface-curvature) effects on shear vector measurement accuracy. In preliminary experiments, surface-inclination angles theta(sub x), theta(sub y) of 0, +/-5, +/-10, and +/-15 deg were employed. In this arrangement, white-light illumination was positioned normal to the untilted test surface, and the camera above-plane view angle was set at 30 deg relative to the untilted test surface. As can be seen, vector-aligned lambda(sub d) values showed no dependence on theta(sub x) or theta(sub y) for absolute values of these tilt angles is less than or equal to 15 deg. Acquisition and analyses of full-surface color images are presently underway to definitively document the insensitivity limits of the shear vector measurement methodology to surface-slope variations.

Reda, Daniel C.↗

Onset of Convection in Ice I with Composite Newtonian and Non-Newtonian Rheology: Application to the Icy Galilean Satellites

Ice I exhibits a complex rheology at temperature and pressure conditions appropriate for the interiors of the ice I shells of Europa, Ganymede, and Callisto. We use numerical methods and existing parameterizations of the critical Rayleigh number to determine the conditions required to trigger convection in an ice I shell with the stress-, temperature- and grain size- dependent rheology measured in laboratory experiments by Goldsby and Kohlstedt [2001]. The critical Rayleigh number depends on the ice grain size and the amplitude and wavelength of temperature perturbation issued to an initially conductive ice I shell. If the shells have an assumed uniform grain size less than 0.4 mm, deformation during initial plume growth is accommodated by Newtonian volume diffusion. If the ice grain size is between 0.4 mm and 3 cm, deformation during plume growth is accommodated by weakly non-Newtonian grain boundary sliding, where the critical ice shell thickness for convection depends on the amplitude of temperature perturbation to the _0.5 power. If the ice grain size exceeds 2 cm, convection can not occur in the ice I shells of the Galilean satellites regardless of the amplitude or wavelength of temperature perturbation. If the grain size in a convecting ice I shell evolves to effective values greater than 2 cm, convection will cease. If the ice shell has a grain size large enough to permit flow by dislocation creep, the ice is too stiff to permit convection, even in the thickest possible ice I shell. Consideration of the composite rheology implies that estimates of the grain size in the satellites and knowledge of their initial thermal states are required when judging the convective instability of their ice I shells.

Barr, A. C.↗

Numerical Simulations of Non-Newtonian Convection in Ice: Application to Europa

Numerical simulations of solid state convection in Europa's ice shell have so far been limited to consideration of Newtonian flow laws, where the viscosity of ice is strongly dependent upon temperature, predicting that a stagnant lid should form at the top (10-40%) of a convecting ice shell. Such large thicknesses seem to contradict estimates of the effective elastic thickness of Europa s ice shell during its geologically active period. Recent laboratory experiments characterize the rheology of ice as the sum of contributions from several temperature and strain rate-dependent creep mechanisms. We present the results of numerical simulations of convection within Europa s ice shell using the finite-element model Citcom, applying the non-Newtonian rheology of grain boundry sliding. Our calculations suggest a shallower brittle/ductile transition and larger interior convective velocities compared to Newtonian rheology. The flow field is time-dependent, with small, localized upwellings and downwellings at the thermal boundary layers that have minimal topographic expression at the surface.

Barr, A. C.↗

Asperity lubrication and cavitation in face seals

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.

Walowitt, J. A.↗

Mantle plumes - A boundary layer approach for Newtonian and non-Newtonian temperature-dependent rheologies

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.

Yuen, D. A.↗

What can we learn about impact mechanics from large craters on Venus?

More than 50 unequivocal peak-ring craters and multiringed impact basins have been identified on Venus from Earth-based Arecibo, Venera 15/16, and Magellan radar images. These ringed craters are relatively pristine, and so serve as an important new dataset that will further understanding of the structural and rheological properties of the venusian surface and of impact mechanics in general. They are also the most direct analogues for craters formed on the Earth in Phanerozoic time. Finite-element simulations of basin collapse and ring formation were undertaken in collaboration with V. J. Hillgren (University of Arizona). These calculations used an axisymmetric version of the viscoelastic finite element code TECTON, modeled structures on the scale of Klenova or Meitner, and demonstrated two major points. First, viscous flow and ring formation are possible on the timescale of crater collapse for the sizes of multiringed basins seen on Venus and heat flows appropriate to the plant. Second, an elastic lithosphere overlying a Newtonian viscous asthenosphere results mainly in uplift beneath the crater. Inward asthenospheric flow mainly occurs at deeper levels. Lithospheric response is dominantly vertical and flexural. Tensional stress maxima occur and ring formation by normal faulting is predicted in some cases, but these predicted rings occur too far out to explain observed ring spacings on Venus (or on the Moon). Overall, these estimates and models suggest that multiringed basin formation is indeed possible at the scales observed on Venus. Furthermore, due to the strong inverse dependence of solid-state viscosity on stress, the absence of Cordilleran-style ring faulting in craters smaller than Meitner or Klenova makes sense. The apparent increase in viscosity of shock-fluidized rock with crater diameter, greater interior temperatures accessed by larger, deeper craters, and decreased non-Newtonian viscosity associated with larger craters may conspire to make the transition with diameter from peak-ring crater to Orientale-type multiringed basin rather abrupt.

Mckinnon, William B.↗

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗

Scaling of Shear Rheology of Concentrated Charged Colloidal Suspensions across Glass Transition

Electrostatic interparticle interactions are a key component in controlling and designing rheological characteristics of concentrated charged colloidal suspensions. Herein, we investigate electroviscous effects on shear rheology using highly charged silica particles. By fixing the volume fraction but varying the salinity, the system undergoes a glass transition as evidenced by the evolution of the yield stress and zero-shear viscosity. We show that the steady shear viscosities obey a critical scaling relation that scales the flow curves into a super- and a sub-critical branch with glass transition salinity serving as the bifurcation point; we also demonstrate an isoviscosity scaling that collapses all isoviscosity lines into a single master curve that exhibits no singularity. Based on each scaling relation, in conjunction with common modeling equations, the quantitative relationships between the shear viscosity, stress, and salinity are established. This study demonstrates a new framework to model the steady shear rheology of concentrated charged colloids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights from high-fidelity modeling of industrial rotary bell atomization

The global automotive industry sprayed over 2.6 billion liters of paint in 2018, much of which through electrostatic rotary bell atomization, a highly complex process involving the fluid mechanics of rapidly rotating thin films tearing apart into micrometer-thin filaments and droplets. Coating operations account for 65% of the energy usage in a typical automotive assembly plant, representing 10,000s of gigawatt-hours each year in the United States alone. Optimization of these processes would allow for improved robustness, reduced material waste, increased throughput, and significantly reduced energy usage. Here, we introduce a high-fidelity mathematical and algorithmic framework to analyze rotary bell atomization dynamics at industrially relevant conditions. Our approach couples laboratory experiment with the development of robust non-Newtonian fluid models; devises high-order accurate numerical methods to compute the coupled bell, paint, and gas dynamics; and efficiently exploits high-performance supercomputing architectures. These advances have yielded insight into key dynamics, including i) parametric trends in film, sheeting, and filament characteristics as a function of fluid rheology, delivery rates, and bell speed; ii) the impact of nonuniform film thicknesses on atomization performance; and iii) an understanding of spray composition via primary and secondary atomization. These findings result in coating design principles that are poised to improve energy- and cost-efficiency in a wide array of industrial and manufacturing settings.

97 MATHEMATICS AND COMPUTING↗

Flow-Driven Stretch Fluctuations Govern the Nonlinear Viscoelasticity of Elongating Associative Polymer Networks

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.

biomimetic & bio-inspired materials↗

Filament extension atomizers: A study of integration potential into spray dryers for high-viscosity and non-Newtonian fluids

Filament extension atomizers (FEAs) are an emerging class of spray nozzles designed to atomize high-viscosity and non-Newtonian fluids that are challenging for conventional pressure or two-fluid nozzles. In this study, we investigate the influence of roller geometry, surface velocity, roller material, and fluid rheology on the atomization performance of FEA systems using concentrated whey protein suspensions (50–70 wt. %). Extensional and shear rheology experiments, along with high-speed imaging and particle image velocimetry, reveal that filament breakup dynamics are governed by competition between inertial, capillary, and viscoelastic stresses. High roller rotational velocity leads to narrower spray cones, contradicting rheology experiments and suggests a significant inertial contribution to filament breakup. Smaller rollers operating at the same rotational velocity led to broader spray cones consistent with expectations. An FEA nozzle was integrated into a conventional dryer producing particles in the 100 μm range. Results suggest that FEA technology enables atomization of highly viscous fluids at industrially relevant spray cone angles similar to those generated by pressure nozzles, offering a pathway to improve energy efficiency in spray drying by enabling higher solids loading feedstocks. Furthermore, these insights provide critical guidance for optimizing FEA nozzle designs and process parameters across a range of applications.

Energy efficiency↗

Topographically driven crustal flow and its implication to the development of pinned oroclines

Pinned oroclines, a type of curved orogen which results from lateral pinning of a growing fold-thrust belt, tend to resemble parabolic Newtonian curvature modified by different degrees of flattening at the flow front. It is proposed that such curves can be generated by Newtonian crustal flow driven by topographic variations. In this model, regional topographic differences create a regional flow which produces a parabolic flow front on interaction with lateral bounding obstacles. Local topographic variations modify the parabolic curves and yield more flat-crested, non-Newtonian-type curvatures. A finite-difference thin-skin tectonic simulation demonstrates that both Newtonian and non-Newtonian curved orogens can be produced within a Newtonian crust.

Hsui, Albert T.↗

Hard turbulent thermal convection and thermal evolution of the mantle

Hard turbulent convection is investigated using laboratory experiments and numerical simulations. In Newtonian mantle convection, the appearance of disconnected plumes marks the transition from soft to hard turbulence. For non-Newtonian rheology, the transition to hard turbulence takes place at much lower Nusselt numbers than it does for Newtonian rheology. This has important ramifications for the mantle. Large curvatures are developed in the trajectories of non-Newtonian plumes in the hard turbulent regime, in contrast to the trajectories of Newtonian plumes. When phase transitions are considered, mantle convection tends to become more layered with increasing Rayleigh numbers. The manner of mantle convection might have changed with time from a layered to a more whole mantle type of flow. Superplume events could have been caused by catastrophic overturns associated with strong gravitational instabilities in the transition zone.

Yuen, D. A.↗

The Extensional Rheology of Non-Newtonian Materials

The evolution of the transient extensional stresses in dilute and semi-dilute viscoelastic polymer solutions are measured with a filament stretching rheometer of a design similar to that first introduced by Sridhar, et al. The solutions are polystyrene-based (PS) Boger fluids that are stretched at constant strain rates ranging from 0.6 less than or equal to epsilon(0) less than or equal to 4s(exp -1) and to Hencky strains of epsilon greater than 4. The test fluids all strain harden and Trouton ratios exceeding 1000 are obtained at high strains. The experimental data strain hardens at lower strain levels than predicted by bead-spring FENE models. In addition to measuring the transient tensile stress growth, we also monitor the decay of the tensile viscoelastic stress difference in the fluid column following cessation of uniaxial elongation as a function of the total imposed Hencky strain and the strain rate. The extensional stresses initially decay very rapidly upon cessation of uniaxial elongation followed by a slower viscoelastic relaxation, and deviate significantly from FENE relaxation predictions. The relaxation at long times t is greater than or equal to 5 s, is compromised by gravitational draining leading to non-uniform filament profiles. For the most elastic fluids, partial decohension of the fluid filament from the endplates of the rheometer is observed in tests conducted at high strain rates. This elastic instability is initiated near the rigid endplate fixtures of the device and it results in the progressive breakup of the fluid column into individual threads or 'fibrils' with a regular azimuthal spacing. These fibrils elongate and bifurcate as the fluid sample is elongated further. Flow visualization experiments using a modified stretching device show that the instability develops as a consequence of an axisymmetry-breaking meniscus instability in the nonhomogeneous region of highly deformed fluid near the rigid endplate.

Spiegelberg, Stephen H.↗