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 577 records · Page 32

Miniature Laser Doppler Velocimeter for Measuring Wall Shear

A miniature optoelectronic instrument has been invented as a nonintrusive means of measuring a velocity gradient proportional to a shear stress in a flow near a wall. The instrument, which can be mounted flush with the wall, is a variant of a basic laser Doppler velocimeter. The laser Doppler probe volume can be located close enough to the wall (as little as 100 micron from the surface) to lie within the viscosity-dominated sublayer of a turbulent boundary layer. The instrument includes a diode laser, the output of which is shaped by a diffractive optical element (DOE) into two beams that have elliptical cross sections with very high aspect ratios.

Gharib, Morteza↗

Atmospheric flow over two-dimensional bluff surface obstructions

The phenomenon of atmospheric flow over a two-dimensional surface obstruction, such as a building (modeled as a rectangular block, a fence or a forward-facing step), is analyzed by three methods: (1) an inviscid free streamline approach, (2) a turbulent boundary layer approach using an eddy viscosity turbulence model and a horizontal pressure gradient determined by the inviscid model, and (3) an approach using the full Navier-Stokes equations with three turbulence models; i.e., an eddy viscosity model, a turbulence kinetic-energy model and a two-equation model with an additional transport equation for the turbulence length scale. A comparison of the performance of the different turbulence models is given, indicating that only the two-equation model adequately accounts for the convective character of turbulence. Turbulence flow property predictions obtained from the turbulence kinetic-energy model with prescribed length scale are only insignificantly better than those obtained from the eddy viscosity model. A parametric study includes the effects of the variation of the characteristics parameters of the assumed logarithmic approach velocity profile. For the case of the forward-facing step, it is shown that in the downstream flow region an increase of the surface roughness gives rise to higher turbulence levels in the shear layer originating from the step corner.

Bitte, J.↗

Nonlinear evolution of subsonic and supersonic disturbances on a compressible free shear layer

The effects of a nonlinear-nonequilibrium-viscous critical layer on the spatial evolution of subsonic and supersonic instability modes on a compressible free shear layer is considered. It is shown that the instability wave amplitude is governed by an integrodifferential equation with cubic-type nonlinearity. Numerical and asymptotic solutions to this equation show that the amplitude either ends in a singularity at a finite downstream distance or reaches an equilibrium value, depending on the Prandtl number, viscosity law, viscous parameter and a real parameter which is determined by the linear inviscid stability theory. A necessary condition for the existence of the equilibrium solution is derived, and whether or not this condition is met is determined numerically for a wide range of physical parameters including both subsonic and supersonic disturbances. it is found that no equilibrium solution exists for the subsonic modes unless the temperature ratio of the low-to-high-speed streams exceeds a critical value, while equilibrium solutions for the most rapidly growing supersonic mode exist over most of the parameter range examined.

Leib, S. J.↗

Rheology guiding the design and printability of aqueous colloidal composites for additive manufacturing

Vat photopolymerization (VP) and direct ink write (DIW) additive manufacturing (AM) provide complex geometries with precise spatial control employing a vast array of photo-reactive polymeric systems. Although VP is recognized for superior resolution and surface finish, DIW provides versatility for higher viscosity systems. However, each AM platform presents specific rheological requirements that are essential for successful 3D printing. First, viscosity requirements constrain VP polymeric materials to viscosities below 10 Pa s. Thus, this requirement presents a challenging paradox that must be overcome to attain the physical performance of high molecular weight polymers while maintaining suitable viscosities for VP polymeric materials. Second, the necessary rheological complexity that is required for DIW pastes requires additional rheological measurements to ensure desirable thixotropic behavior. Further, this manuscript describes the importance of rheological measurements when designing polymeric latexes for AM. Latexes effectively decouple the dependency of viscosity on molecular weight, thus enabling high molecular weight polymers with low viscosities. Photo-crosslinking of water-soluble monomers and telechelic oligomeric diacrylates in the presence of the latex enables the fabrication of a scaffold, which is restricted to the continuous aqueous phase and effectively surrounds the latex nanoparticles enabling the printing of otherwise inaccessible high molecular weight polymers. Rheological testing, including both steady and oscillatory shear experiments, provides insights into system properties and provides predictability for successful printing. This perspective article aims to provide an understanding of both chemical functionality (photo- and thermal-reactivity) and rheological response and their importance for the successful design and evaluation of VP and DIW processable latex formulations.

36 MATERIALS SCIENCE↗

Measurement and prediction of flow through a replica segment of a mildly atherosclerotic coronary artery of man

Pressure distributions were measured along a hollow vascular axisymmetric replica of a segment of the left circumflex coronary artery of man with mildly atherosclerotic diffuse disease. A large range of physiological Reynolds numbers from about 60 to 500, including hyperemic response, was spanned in the flows investigation using a fluid simulating blood kinematic viscosity. Predicted pressure distributions from the numerical solution of the Navier-Stokes equations were similar in trend and magnitude to the measurements. Large variations in the predicted velocity profiles occurred along the lumen. The influence of the smaller scale multiple flow obstacles along the wall (lesion variations) led to sharp spikes in the predicted wall shear stresses. Reynolds number similarity was discussed, and estimates of what time averaged in vivo pressure drop and shear stress might be were given for a vessel segment.

Back, L. H.↗

Large eddy simulation of turbulent channel flow: ILLIAC 4 calculation

The three-dimensional time dependent equations of motion were numerically integrated for fully-developed turbulent channel flow. A large scale flow field was obtained directly from the solution of these equations, and small scale field motions were simulated through an eddy viscosity model. The calculations were carried out on the ILLIAC 4 computer. The computed flow patterns show that the wall layer consists of coherent structures of low speed and high speed streaks alternating in the spanwise direction. These structures were absent in the regions away from the wall. Hot spots, small localized regions of very large turbulent shear stress, were frequently observed. The profiles of the pressure velocity-gradient correlations show a significant transfer of energy from the normal to the spanwise component of turbulent kinetic energy in the immediate neighborhood of the wall ('the splatting effect').

Kim, J.↗

Synthesis and Properties of Phenylethynyl-Terminated, Star-Branched, Phenylquinoxaline Oligomers

The primary objective of this work was to prepare readily melt and solution processable phenylquinoxaline (PQ) oligomers that could be thermally crosslinked to solvent-resistant resins. Thus, a mixture of 2-(4-hydroxyphenyl)-3-phenyl-6-fluoroquinoxaline and 3-(4-hydroxyphenyl)-2-phenyl-6-fluoroquinoxaline (HPFQ) was used to prepare star-branched PQ oligomers end-capped with 4-fluoro-4-phenylethynylbenzophenone (FPEB). 1,1,1-Tris(4-hydroxyphenyl)ethane (THPE) was used as the branching unit. The oligomer number-average molecular weights (M (bar) (sub n) S) as determined by size exclusion chromatography (SEC) were close to the calculated values of 2922, 4698, 6474, and 13,578 g/mol, and their intrinsic viscosities ranged from 0.16 to 0.57 dl/g (m-cresol at 30 C). The oligomers, which were quite soluble in common organic solvents, had glass transition temperatures (T (sub g) S) that ranged from 181 to 233 C (DSC, DELTA T = 20 C/min). They also underwent an exothermic cure with maxima from 377 to 443 C. The T (sub g) S of the cured oligomers ranged from 259 to 284 C depending on the oligomer M (bar) (sub n) and the curing conditions. The oligomers had low melt viscosities, e.g. an oligomer (SPQ-46) with an M (bar) (sub n) of 4816 g/mol (SEC) had a melt viscosity of 150 Pa s at 348 C. A cured thin film of SPQ-46, which was insoluble in common organic solvents, had a room temperature (RT) tensile strength of 100 MPa, a RT modulus of 2358 MPa, and a RT elongation of 5.9%. A cured sample of SPQ-46 displayed a RT titanium-titanium lap shear tensile strength of 35.2 MPa. SPQ-46/carbon fiber(IM-7) composites, were prepared that displayed a RT flexural strength of 1902 MPa, a RT modulus of 1.38 GPa and a RT open hole compressive strength of 433 MPa.

Ooi, I. H.↗

Turbulence models and Reynolds analogy for two-dimensional supersonic compression ramp flow

Results of the application of turbulence models and the Reynolds analogy to the Navier-Stokes computations of Mach 2.9 two-dimensional compression ramp flows are presented. The Baldwin-Lomax eddy viscosity model and the kappa-epsilon turbulence transport equations for the turbulent momentum flux modeling in the Navier-Stokes equations are studied. The Reynolds analogy for the turbulent heat flux modeling in the energy equation was also studied. The Navier-Stokes equations and the energy equation were numerically solved for the flow properties. The Reynolds shear stress, the skin friction factor, and the surface heat transfer rate were calculated and compared with their measurements. It was concluded that with a hybrid kappa-epsilon turbulence model for turbulence modeling, the present computations predicted the skin friction factors of the 8 deg and 16 deg compression ramp flows and with the turbulent Prandtl number Pr(sub t) = 0.93 and the ratio of the turbulent thermal and momentum transport coefficients mu(sub q)/mu(sub t) = 2/Prt, the present computations also predicted the surface heat transfer rates beneath the boundary layer flow of the 16 compression ramp.

Wang, Chi R.↗

Computation of aircraft component flow fields at transonic Mach numbers using a three-dimensional Navier-Stokes algorithm

A computer analysis was developed for calculating steady (or unsteady) three-dimensional aircraft component flow fields. This algorithm, called ENS3D, can compute the flow field for the following configurations: diffuser duct/thrust nozzle, isolated wing, isolated fuselage, wing/fuselage with or without integrated inlet and exhaust, nacelle/inlet, nacelle (fuselage) afterbody/exhaust jet, complete transport engine installation, and multicomponent configurations using zonal grid generation technique. Solutions can be obtained for subsonic, transonic, or hypersonic freestream speeds. The algorithm can solve either the Euler equations for inviscid flow, the thin shear layer Navier-Stokes equations for viscous flow, or the full Navier-Stokes equations for viscous flow. The flow field solution is determined on a body-fitted computational grid. A fully-implicit alternating direction implicit method is employed for the solution of the finite difference equations. For viscous computations, either a two layer eddy-viscosity turbulence model or the k-epsilon two equation transport model can be used to achieve mathematical closure.

Shrewsbury, George D.↗

Turbulent Diffusion and Concentration of "Chondrules" in the Protoplanetary Nebula

Turbulence is known to possess structure on many scales. The largest or integral scale L is usually taken to be the largest relevant dimension of the system; for the nebula as a whole, it would be on the order of a vertical gas scale height H approximately 0.1R where R is the distance from the sun. Turbulent kinetic energy cascades towards smaller scales; for fully developed, homogeneous turbulence, the turbulent kinetic energy peaks at the integral scale and is distributed within its "inertial range" according to the Kolmogorov energy spectrum dE(k) = E(sub 0)(kL)(sup -5/3)dk where k is the wavenumber or inverse wavelength. The smallest scale in a turbulent regime is the Kolmogorov scale (n)(sometimes referred to as the inner scale). On this scale, molecular viscosity can dissipate turbulent kinetic energy . Turbulent flows are characterized by their Reynolds number Re = UL/v, where U is a characteristic velocity scale and v is molecular viscosity. More energetic (higher Re) flows drive turbulence to smaller Kolmogorov scales: n = LRe(sup -3/4). Two different turbulent regimes are likely to be of importance; early stage convective turbulence in a hot nebula probably extends over the entire nebula scale height, and probably has typical eddy frequency comparable to the orbit frequency. At a later stage, midplane shear layer turbulence can be driven by the presence of a differentially rotating, densely settled particle layer, with different length and time scales.

Cuzzi, Jeffrey N.↗

The coupling of interfacial instabilities and the stabilization of two-layer annular flows

In this paper the stability of annular pressure-driven parallel flows of two liquids sandwiching a free cylindrical interface is considered. For small to moderate Reynolds numbers, the interface is susceptible to capillary and interfacial wave instabilities, the latter instability caused by a jump in viscosity at the interface. It is shown that favorable velocity profiles in both liquids may stabilize capillary breakup of the interface and suppress the axisymmetric interfacial wave instability. A long-wave analysis leads to the physical mechanism responsible for stabilization of capillary breakup. This physical mechanism is a generalization of that by which capillary breakup is stabilized by interfacial shear in an annular film of a single liquid. Stabilization of intermediate wavelengths is studied with a mechanical energy analysis, which leads to a description of the energetic processes at work.

Dijkstra, Henk A.↗

Properties of Desert Sand and CMAS Glass

As-received desert sand from a Middle East country has been characterized for its phase composition and thermal stability. X-ray diffraction analysis showed the presence of quartz (SiO2), calcite (CaCO3), gypsum (CaSO4.2H2O), and NaAlSi3O8 phases in as-received desert sand and showed weight loss of approx. 35 percent due to decomposition of CaCO3 and CaSO4.2H2O when heated to 1400 C. A batch of as-received desert sand was melted into calcium magnesium aluminosilicate (CMAS) glass at approx. 1500 C. From inductively coupled plasma-atomic emission spectrometry, chemical composition of the CMAS glass was analyzed to be 27.8CaO-4MgO-5Al2O3-61.6SiO2-0.6Fe2O3-1K2O (mole percent). Various physical, thermal and mechanical properties of the glass have been evaluated. Bulk density of CMAS glass was 2.69 g/cc, Young's modulus 92 GPa, Shear modulus 36 GPa, Poisson's ratio 0.28, dilatometric glass transition temperature (T (sub g)) 706 C, softening point (T (sub d)) 764 C, Vickers microhardness 6.3 +/- 0.4 GPa, indentation fracture toughness 0.75 +/- 0.15 MPa.m (sup 1/2), and coefficient of thermal expansion (CTE) 9.8 x 10 (exp -6)/degC in the temperature range 25 to 700 C. Temperature dependence of viscosity has also been estimated from various reference points of the CMAS glass using the Vogel-Fulcher-Tamman (VFT) equation. The glass remained amorphous after heat treating at 850 C for 10 hr but crystallized into CaSiO3 and Ca-Mg-Al silicate phases at 900 C or higher temperatures. Crystallization kinetics of the CMAS glass has also been investigated by differential thermal analysis (DTA). Activation energies for the crystallization of two different phases in the glass were calculated to be 403 and 483 kJ/mol, respectively.

Desert sand↗

Filler/ Polycarbosilane Systems as CMC Matrix Precursors

Pyrolytic conversion of polymeric precursors to ceramics is accompanied by loss of volatiles and large volume changes. Infiltration of a low viscosity polymer into a fiber preform will fill small spaces within fiber tows by capillary forces, but create large matrix cracks within large, intertow areas. One approach to minimizing shrinkage and reducing the number of required infiltration cycles is to use particulate fillers. In this study, Starfire allylhydridopolycarbosilane (AHPCS) was blended with a silicon carbide powder, with and without dispersant, using shear mixing. The polymer and polymer/particle interactions were characterized using nuclear magnetic resonance, differential scanning calorimetry, thermogravimetric analysis and rheometry. Polymer/particulate slurries and suspensions were used to infiltrate a figidized preform of an eight ply five harness satin CG Nicalon fiber having a dual layer BN/SiC interface coating, and the resulting composites characterized by optical and scanning electron microscopy.

Hurwitz, Frances I.↗

Pterodactyl: Aerodynamic and Aerothermal Modeling for a Symmetric Deployable Earth Entry Vehicle with Flaps

NASA’s Pterodactyl project has investigated a deployable atmospheric entry vehicle integrated with a flap control system that provides precision targeting during reentry. The control system consists of eight flaps mounted at the edge of a heatshield that can deflect in and out of the flow. An aerodynamic and aerothermodynamic analysis process was developed for preliminary vehicle design and aerodynamic database generation using engineering and CFD tools with varying levels of fidelity. The objective of this analysis was to: 1) understand the inherent aerodynamics, 2) provide an aerodynamics database for stability and control analysis, and 3) provide qualitative and quantitative aeroheating analysis for Thermal Protection System modeling of the flaps. A high fidelity Euler code, Cart3D, was used to resolve complex flow features such as secondary shocks and shock impingement. A wide range of supersonic (Mach 2) and hypersonic Mach numbers (up to Mach 40) were tested with an Earth atmosphere model. A process was developed to utilize the adaptive volume mesh generator utility of the perfect gas model of Cart3D to create meshes with high cell efficiency and numerical stability prior to running a 2nd order accurate solution with a real gas model. High-fidelity aerothermal CFD simulations were performed using US3D to further improve the aerothermal analysis on hypersonic flow around a complex entry vehicle shape including viscosity, chemical reactions of air species, vibrational energy, and catalytic surface reactions. In this process, it was found that the flap control system provided multi-axis control that can be utilized for entry precision targeting. Additionally, the Pterodactyl vehicle can achieve up to a trim L/D of 0.2. Finally, the increased fidelity of the aerothermal heating environments revealed that the shear stress contributes to increasing the heating on the flaps.

Deployable Entry Vehicle↗

Pterodactyl: Aerodynamic and Aerothermal Modeling for a Symmetric Deployable Earth Entry Vehicle with Flaps

NASA’s Pterodactyl project has investigated a deployable atmospheric entry vehicle integrated with a flap control system that provides precision targeting during reentry. The control system consists of eight flaps mounted at the edge of a heatshield that can deflect in and out of the flow. An aerodynamic and aerothermodynamic analysis process was developed for preliminary vehicle design and aerodynamic database generation using engineering and CFD tools with varying levels of fidelity. The objective of this analysis was to: 1) understand the inherent aerodynamics, 2) provide an aerodynamics database for stability and control analysis, and 3) provide qualitative and quantitative aeroheating analysis for Thermal Protection System modeling of the flaps. A high fidelity Euler code, Cart3D, was used to resolve complex flow features such as secondary shocks and shock impingement. A wide range of supersonic (Mach 2) and hypersonic Mach numbers (up to Mach 40) were tested with an Earth atmosphere model. A process was developed to utilize the adaptive volume mesh generator utility of the perfect gas model of Cart3D to create meshes with high cell efficiency and numerical stability prior to running a 2nd order accurate solution with a real gas model. High-fidelity aerothermal CFD simulations were performed using US3D to further improve the aerothermal analysis on hypersonic flow around a complex entry vehicle shape including viscosity, chemical reactions of air species, vibrational energy, and catalytic surface reactions. In this process, it was found that the flap control system provided multi-axis control that can be utilized for entry precision targeting. Additionally, the Pterodactyl vehicle can achieve up to a trim L/D of 0.2. Finally, the increased fidelity of the aerothermal heating environments revealed that the shear stress contributes to increasing the heating on the flaps.

Deployable Entry Vehicle↗

Turbulence Model Validation through Joint Experimental / Computational Studies of Separated Flow Over a Three-Dimensional Tapered Bump: Part II - Computational Investigation

This report serves as documentation of activities related to Computational Fluid Dynamics (CFD) analysis applied to the Smooth Body Separation Experiment designed to produce validation-quality experimental data from wind-tunnel tests performed at the University of Notre Dame. All CFD analysis utilized the Boeing CFD (BCFD) flow solver. BCFD was used to analyze flow over the Boeing Speed Bump with second-order accuracy in space and time. Inviscid fluxes were computed using the HLLE scheme for steady-state simulations, and the Low-Dissipation Roe scheme for time-accurate simulations. The working fluid was modeled as calorically perfect air using Sutherland's law for viscosity. Thermal diffusivity was determined by assuming a constant molecular Prandtl number of 0.72, and a turbulent Prandtl number of 0.9. Reynolds Averaged Navier-Stokes (RANS) simulations were conducted using the one-equation Spalart-Allmaras turbulence model with rotation correction3 and the quadratic constitutive relation (SA-RC-QCR), as well as Menter's two-equation, Shear-Stress Transport (SST) model. Time-accurate computations were conducted using Delayed Detached eddy Simulation (DDES) with the one-equation RANS model, and Large Eddy Stress balance (LESb) with the two-equation RANS model. Implicit dual-time stepping was used to achieve second-order temporal accuracy. The experimental investigation conducted during the course of this study is described in detail in an accompanying report.

Flow separation↗

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The rigid shell component for superrotation in planetary atmospheres: Angular momentum budget, mechanical analog and simulation of the spin up process

An analysis of superrotation in the atmosphere of planets, with rotation axis perpendicular to the orbital plane is presented. As the atmosphere expands, Hadley cells develop producing a redistribution of mass and angular momentum. A three dimensional thermally driven zonally symmetric spectral model and Laplace transformation simulate the time evolution of a fluid leading from corotation under globally uniform heating to superrotation under globally nonuniform heating. For high viscosities the rigid shell component of atmospheric superrotation can be understood in analogy with a pirouette. During spin up angular momentum is transferred to the planet. For low iscosities, the process is reversed. A tendency toward geostrophy, combined with increase of surface pressure toward the poles (due to meridional mass transport), induces the atmosphere to subrotate temporarily at lower altitudes. Resultant viscous shear near the surface permits angular momentum to flow from the planet into the atmosphere propagating upwards to produce high altitude superrotation rates.

Mayr, H. G.↗