Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “viscosity, shear”

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 451 records · Page 25

Zonal Two Equation Kappa-Omega Turbulence Models for Aerodynamic Flows

Two new versions of the kappa-omega two-equation turbulence model will be presented. The new Baseline (BSL) model is designed to give results similar to those of the original kappa-omega model of Wilcox, but without its strong dependency on arbitrary freestream values. The BSL model is identical to the Wilcox model in the inner 50% of the boundary-layer but changes gradually to the standard kappa-epsilon model (in a kappa- omega formulation) towards the boundary-layer edge. The free shear layers. The second version of the model is called Shear-Stress Transport (SST) model. It is a variation of the BSL model with the additional ability to account for the transport of the principal turbulent shear stress in adverse pressure gradient boundary-layers. The model is based on Bradshaw's assumption that the principal shear-stress is proportional to the turbulent kinetic energy, which is introduced into the definition of the eddy-viscosity. Both models are tested for a large number of different flowfields. The results of the BSL model are similar to those of the original kappa-omega model, but without the undesirable freestream dependency. The predictions of the SST model are also independent of the freestream values but show better agreement with experimental data for adverse pressure gradient boundary-layer flows.

Menter, Florian R.↗

Synthesis and Characterization of Modified Phenylethynyl Imides

As an ongoing effort to develop structural adhesives for high-performance aerospace applications, recent work has focused on phenylethynyl terminated imide (PETI) oligomers. The work reported herein involves the synthesis and characterization of a series of phenylethynyl containing oligomers designated LARC(TM) MPEI (modified phenylethynyl imide). These oligomers presumably contain mixtures of linear, branched and star-shaped molecules. The fully imidized polymers exhibited minimum melt viscosities as low as 600 poise at 335 C, significantly lower than equivalent molecular weight linear materials. Ti/Ti lap shear specimens processed at 288 C under 15 psi showed tensile shear strengths as high as approx. 6000 psi and 5200 psi at ambient temperature and 177 C respectively. The chemistry and properties of these new MPEIs are presented and compared with an optimized linear PETI, LARC(TM)PETI-5.

Jensen, B. J.↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D for Axisymmetric Jet Flow

Two full seven-equation turbulence models available in the FUN3D code are evaluated for their ability to improve the computation of challenging mixing flows encountered in aerospace propulsion. These models are the SSG/LRR and Wilcox full second-moment Reynolds stress models. They solve equations for the six components of the Reynolds stress and a seventh equation for the turbulent length scale. Two standard eddy viscosity models are also evaluated for comparison, the Spalart-Allmaras (SA) one-equation model and the Menter Shear Stress Transport (SST-V) two-equation turbulence model. Flow through an axisymmetric reference nozzle is examined at three flow conditions: subsonic unheated, subsonic heated, and near sonic unheated. Centerline profiles of velocity and turbulent kinetic energy and radial profiles of velocity, turbulent kinetic energy and turbulent stresses are examined. Results showed that the SA model did well at predicting the jet potential core length, but over-mixed the downstream flow, whereas the SST-V model over-predicted the potential core length. The Wilcox-model significantly over-predicted the potential core length and under-predicted the mixing and was not well-suited for the jet flows evaluated, however the SSG/LRR Reynolds stress model did well at predicting the mixing rate and mean velocity for all cases examined.

mixing layers↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D for Axisymmetric Jet Flow

Two full seven-equation turbulence models have been implemented into the FUN3D code to evaluate their ability to improve the computation of challenging flows encountered in aerospace propulsion, including mixing flows. These models are the SSG/LRR and Wilcox full second-moment Reynolds stress models. They solve equations for the six components of the Reynolds stress and a seventh equation for the mixing length. Two standard eddy viscosity models are also evaluated for comparison, the Spalart-Allmaras (SA) one-equation model and the Menter Shear Stress Transport (SST-V) two-equation turbulence model. Flow through an axisymmetric reference nozzle is examined at three flow conditions: subsonic unheated, subsonic heated, and near sonic unheated. Centerline profiles of velocity and turbulent kinetic energy and radial profiles of velocity, turbulent kinetic energy and turbulent stresses are examined. characteristics, no significant changes in the downstream flow behavior compared to the baseline case are observed. Furthermore, the total power consumed by the fans for different incoming flow conditions also remain marginally the same. It is hoped that the results, albeit obtained at very low speeds. would serve as a database for this technologically interesting flow field that has not been explored adequately before.

computational fluid dynamics↗

Magnetic field draping at the low-latitude magnetopause

Magnetohydrodynamic simulations are used to investigate the structure of the low-latitude magnetopause for interplanetary magnetic field conditions with a dominant southward component. The structure is self-consistently calculated as an initial-value problem in which the system is allowed to evolve into a quasi-steady state. All components of the 3D velocity and magnetic field as well as compressibility, resistivity, and viscosity are included in the 2D calculation. The simulation model shows that magnetic field draping can occur at the magnetopause boundary when magnetic merging takes place in the presence of a tangential shear flow. For 'normal' (positive Bx) draping, the higher-latitude portion of the field lines are curved toward the sun on the magnetospheric side of the magnetopause and away from the sun on the magnetosheath side. The thickness of the normal draping structure scales with the viscosity. The field-aligned current system that accompanies normal magnetic draping is consistent with the sense of the region 1 currents that flow into the dayside ionosphere.

Richard, R. L.↗

Nonlinear evolution of tidally distorted accretion disks: Two-dimensional simulations

According to a previously published linear analysis, the tidal distortion of accretion disks in binary star systems produces a local hydrodynamic instability to m = 1 internal waves, which may have arbitrarily small wavelengths in the absence of viscosity. The instability is three-dimensional and approximately incrompressible. To explore the nonlinear outcome of this instability, we develop a shearing-sheet approximation on scales comparable to the disk thickness. The large-scale azimuthal variation of the disk is represented by varying the local metric with time (local orbital phase). The hydrodynamic equations can then be posed two-dimensionally on local meridional planes. We solve these equations with a second-order gasdynamical code based on the Total-Variation-Diminishing scheme. Our simulations confirm the predicted linear growth rate. The modes saturate chaotically at velocities scaling as the product of the linear growth rate and the wavelength. If the wavelength is small compared with the disk thickness, the modes remain nearly incompressible even when nonlinear. The two-dimensional power spectrum of velocities after saturation is roughly isotropic and extends over a broad range of scales in an approximately power-law fashion. We measure the heating rate associated with the nonlinear dissipation of the modes. The dissipation implies a secular torque on the disk and a return of angular momentum to the secondary star via the tidal potential. The estimated torque is somewhat larger than the tidal torque produced by maximal disk viscosity (alpha approximately 1). At least in these two-dimensional simulations, however, there is no significant angular momentum flux within the disk.

Ryu, Dongsu↗

A unified theory of zonal flow shears and density corrugations in drift wave turbulence

A unified theory of zonal flow shears and density corrugations in drift wave turbulence is presented. Polarization and density advection beat excitation are studied in combination with modulational response. Noise is driven by two-time flux correlation. While the effective zonal flow eddy viscosity can go negative, the zonal diffusivity is positive definite. There is no inverse cascade of density corrugation. The connection between avalanches and corrugations is discussed. Here, the zonal cross-correlation is identified and calculated. Conditions for alignment of zonal shears and corrugation gradients are determined, and the implications for staircase structure are discussed. We show that the synergy of beat noise and modulational effects is stronger than either alone. Strong zonal flows can be excited well below the modulational instability threshold. In the context of L–H transition, zonal noise quenches turbulence overshoot by eliminating the threshold for zonal flow excitation. The power threshold for L–H transition is lowered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Connecting particle interactions to agglomerate morphology and rheology of boehmite nanocrystal suspensions

Rheology imposes significant challenges on processing of complex suspensions such as nuclear waste slurries at the Hanford and Savannah River sites. Understanding rheology connecting to microstructures and underlying particle interactions in complex slurries is therefore important for both fundamental knowledge and practical applications. Here, we use suspensions of aluminum oxyhydroxide minerals in the form of boehmite as an analog of the radioactive waste slurry to gain physical insights on the correlation between particle interactions, microstructures, and slurry rheology. Specifically, we use a combination of Couette rheometry and small-angle scattering techniques (independently and simultaneously) to understand how the slurry microstructure changes under flow and how these structural changes manifest themselves in the bulk rheology of the suspensions. Our experiments show that the boehmite slurries are thixotropic, with the rheology and structure of the suspensions changing with increasing exposure to flow. In the slurries, particle aggregates begin as loose, system-spanning clusters, but exposure to moderate shear rates causes the aggregates to irreversibly consolidate into denser clusters of finite size. The microstructural changes directly influence the rheological properties of the slurries such as viscosity and viscoelasticity. More importantly, our study shows that solution pH affects the amount of structural rearrangement and the kinetics of the rearrangement process, with an increase in pH leading to faster and more dramatic changes in the bulk rheology. Such dynamic microstructural changes and resultant rheology are understood via correlations between particle interactions and strength of particle network, coupled with surface chemistry and anisotropic nature of particle interactions. Nearly identical structural changes are also observed in Poiseuille flow geometries, implying that the observed changes are relevant in the pipe flow conditions present during waste processing.

Weston, Javen S.↗

Impact of Drying on Meso- and Nanoscale Structures of Citrus Fiber: A Study by SFG, ATR-IR, XRD, and DLS

Citrus fibers are a sidestream of the pectin extraction process from citrus peel. This sidestream can be converted into a functional ingredient through a shear-induced homogenization process. One technical challenge with this material is that dehydration and subsequent rehydration result in reduction of viscosity compared to the original product. In this study, various drying methods were compared with never-dried fibers to investigate the structural changes underlying the viscosity loss. Infrared and X-ray diffraction analyses confirmed no changes in chemical composition and crystalline structure of citrus fibers. Here, the dynamic light scattering and sum frequency generation analyses of citrus fiber suspension showed that the rehydration process could not fully disperse aggregated fibers, which appears to be the main cause for the viscosity loss.

09 BIOMASS FUELS↗

Learning how structures form in drift-wave turbulence

Drift-wave turbulence produces anomalous transport via cross-correlations between fluctuations. This transport has profound implications for confinement, structure formation, and virtually all aspects of the non-linear turbulent dynamics. Herein, we use a data-driven method based on deep learning in order to study turbulent transport in the 2D Hasegawa–Wakatani system and infer a reduced mean-field model from numerical solution. In addition to the usual turbulent diffusion, we find an effect which couples the particle flux to the local gradient of vorticity, which tends to modulate the density profile. The direct coupling to the shear is relatively weak. In addition, the deep learning method finds a model for spontaneous zonal flow generation by negative viscosity, stabilized by non-linear and hyperviscous terms. We compare these results to analytic calculations using quasilinear theory and wave kinetics, finding qualitative agreement, though the calculations miss certain higher-order effects. A simplified, 1-D model for the evolution of the profile, flow, and intensity based on the deep learning results is solved numerically and compared to previous models for staircasing based on bistability. We see that the physics uncovered by the deep learning method provided simple explanations for the formation of zonal structures in the density, flow, and turbulence fields. We highlight the important role of symmetry in the deep learning method and speculate on the portability of the method to other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surface temperatures and glassy state investigations in tribology, part 4

Measurements were made of the limiting shear stress for two naphthenic oils of differing molecular weight and three blends of the lower molecular weight oil and polyalkylmethacrylate polymers of differing molecular weight. The two base oils reached the same limiting shear stress for the same temperature and pressure. This was also true for all the polymer solutions although the polymer reduced the limiting shear stress by about 15 percent. It is shown that limiting stress is more a function of material type than viscosity or molecular weight. A new falling body viscometer was constructed to operate to 230 C and 0.6 GPa. Another viscometer was constructed to extend the pressure range to 1.1 GPa. A concentrated contact simulator was developed which allows recording of the traction force while the slide-roll ratio is continuously varied and the rolling speed is maintained essentially constant by a single drive motor. The configuration is that of a crowned roller against a disk. Measurement of lubricant minimum film thickness of elliptical EHD contacts of various aspect ratios were made by optical interferometry. The data collected were used to evaluate the Hamrock and Dowson minimum film thickness model over a range of contract ellipticity ratio where the major axis of the contact ellipse was aligned both parallel and perpendicular to the direction of motion. A statistical analysis of the measured film thickness data showed that on the average the experimental data were 30 percent greater than the film thickness predicted by the model. Preliminary development of the application of a scanning infrared radiation system to a tribo-system was completed.

Bair, S. S.↗

Experimental investigation of flow about a strut-endwall configuration

An experimental study was conducted to investigate incompressible flow about a strut-endwall configuration positioned within a constant area duct. The endwal boundary layer was tripped, but natural transition occurred on the strut surface. The results indicate that spanwise varying transition on this surface leads to the formation of a secondary vortex which coexists with the horseshoe vortex generated by endwall flow separation upstream of the strut leading edge. Both vortices are similar in strength downstream of the strut trailing edge, and both distort the primary flow and local turbulence structure in the wake-endwall region. The level of distortion is demonstrated by means of axial mean velocity contours, turbulence kinetic energy contours, and Reynolds shear stress contours measured in the cross plane at two streamwise locations. Analysis of the results shows that conventional eddy viscosity and k-epsilon transport equation models are not wholly adequate for predicting this flow situation.

Chang, P. S.↗

Interaction between a normal shock wave and a turbulent boundary layer at high transonic speeds. II - Wall shear stress

Asymptotic methods are used to calculate the shear stress at the wall for the interaction between a normal shock wave and a turbulent boundary layer on a flat plate. A mixing length model is used for the eddy viscosity. The shock wave is taken to be strong enough that the sonic line is deep in the boundary layer and the upstream influence is thus very small. It is shown that unlike the result found for laminar flow an asymptotic criterion for separation is not found; however, conditions for incipient separation are computed numerically using the derived solution for the shear stress at the wall. Results are compared with available experimental measurements.

Liou, M. S.↗

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↗

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↗

Chemical heterogeneity and imperfect mixing in the solar nebula

The 'standard' accretion-disk model of the solar nebula is examined together with its basic equations, which are solved for the steady-state case. It is shown that, for this model of the solar nebula, the efficiency of mixing the products of thermochemical processing from small to large disk radii depends not on the magnitude of the eddy diffusivity but only on the ratio k between the eddy diffusivity and the eddy viscosity. In most cases, this ratio k is probably of order of unity or less, implying that most of the solar nebula was not contaminated by the consequences of thermochemical equilibria that were established at 'small' disk radii. It is noted that, while this condition is almost certainly satisfied if the physical process responsible for viscosity is thermal convection or waves, it has uncertain validity during the early phases of disk evolution where accretion-induced shear instabilities may dominate.

Stevenson, David J.↗

Investigation of turbulent processes in magnetospheric boundary layers

A self-consistent non-evolving two dimensional slab model of a viscous low-latitude boundary layer (LLBL) coupled to the ionosphere was developed by Phan, et al., (1989). Numerical results from the model and possible use of observations to determine the model parameters are discussed. The dynamical model developed by Lotko, et al., (1987) was used by Lotko and Shen (1991) to examine dynamical processes relevant to the LLBL with particular application to post-noon auroral shear layers. Initial results from a magnetohydrodynamic study of flank-side mangetopause boundary configuration are described. Effects of compressibility, scalar viscosity, and electrical resistivity are included in the MHD equations.

Lotko, William↗

A mechanism for tectonic deformation on Venus

In the absence of identifiable physiographic features directly associated with plate tectonics, alternate mechanisms are sought for the intense tectonic deformation observed in radar images of Venus. One possible mechanism is direct coupling into an elastic lithosphere of the stresses associated with convective flow in the interior. Spectral Green's function solutions have been obtained for stresses in an elastic lithosphere overlying a Newtonian interior with an exponential depth dependence of viscosity, and a specified surface-density distribution driving the flow. At long wavelengths and for a rigid elastic/fluid boundary condition, horizontal normal stresses in the elastic lid are controlled by the vertical shear stress gradient and are directly proportional to the depth of the density disturbance in the underlying fluid. The depth and strength of density anomalies in the Venusian interior inferred by analyses of long wavelength gravity data suggest that stresses in excess of 100 MPa would be generated in a 10 km thick elastic lid unless a low viscosity channel occurring beneath the lid or a positive viscosity gradient uncouples the flow stresses. The great apparent depth of compensation of topographic features argues against this, however, thus supporting the importance of the coupling mechanism. If there is no elastic lid, stresses will also be very high near the surface, providing also that the viscosity gradient is negative.

Phillips, Roger J.↗