Apparent departures from newtonian behavior in liquids caused by viscous heating
Departures from Newtonian behavior in liquids caused by viscous heating - ideal Couette flow of incompressible Newtonian liquids
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.
Departures from Newtonian behavior in liquids caused by viscous heating - ideal Couette flow of incompressible Newtonian liquids
A higher order accurate numerical procedure has been developed for solving incompressible Navier-Stokes equations for 2D or 3D fluid flow problems. It is based on low-storage Runge-Kutta schemes for temporal discretization and fourth and sixth order compact finite-difference schemes for spatial discretization. The particular difficulty of satisfying the divergence-free velocity field required in incompressible fluid flow is resolved by solving a Poisson equation for pressure. It is demonstrated that for consistent global accuracy, it is necessary to employ the same order of accuracy in the discretization of the Poisson equation. Special care is also required to achieve the formal temporal accuracy of the Runge-Kutta schemes. The accuracy of the present procedure is demonstrated by application to several pertinent benchmark problems.
Over the past thirty years, numerical methods and simulation tools for incompressible flows have been advanced as a subset of the computational fluid dynamics (CFD) discipline. Although incompressible flows are encountered in many areas of engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to the rather stringent requirements for predicting aerodynamic performance characteristics of flight vehicles, while flow devices involving low-speed or incompressible flow could be reasonably well designed without resorting to accurate numerical simulations. As flow devices are required to be more sophisticated and highly efficient CFD took become increasingly important in fluid engineering for incompressible and low-speed flow. This paper reviews some of the successes made possible by advances in computational technologies during the same period, and discusses some of the current challenges faced in computing incompressible flows.
A through-flow (hub-to-shroud) truly inverse method is proposed and described. It uses an imposition of mean swirl, i.e., radius times mean tangential velocity, given throughout the meridional section of the turbomachine as an initial design specification. In the present implementation, it is assumed that the fluid is inviscid, incompressible, and irrotational at inlet and that the blades are supposed to have zero thickness. Only blade rows that impart to the fluid a constant work along the space are considered. An application of this procedure to design the rotor of a mixed-flow pump is described in detail. The strategy used to find a suitable mean swirl distribution and the other design inputs is also described. The final blade shape and pressure distributions on the blade surface are presented, showing that it is possible to obtain feasible designs using this technique. Another advantage of this technique is the fact that it does not require large amounts of CPU time.
During the past thirty years, numerical methods and simulation tools for incompressible flows have been advanced as a subset of CFD discipline. Even though incompressible flows are encountered in many areas of engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to rather stringent requirements for predicting aerodynamic performance characteristics of flight vehicles, while flow devices involving low speed or incompressible flow could be reasonably well designed without resorting to accurate numerical simulations. As flow devices are required to be more sophisticated and highly efficient, CFD tools become indispensable in fluid engineering for incompressible and low speed flow. This paper is intended to review some of the successes made possible by advances in computational technologies during the same period, and discuss some of the current challenges.
A constitutive relation for an equivalent, homogeneous fluid is developed for the anisotropic viscous flow of an oriented assembly of discontinuous fibers suspended in a viscous fluid. The anisotropic viscous compliance matrix can be expressed in terms of three constants by assuming the equivalent fluid to be incompressible and the microstructure to have axial symmetry (transversely isotropic). By means of a micromechanics analysis, the three terms of the constitutive relation are expressed in terms of the viscosity of the matrix fluid, the fiber aspect ratio, and the fiber volume fraction. A comparison of the viscosity terms reveals that the elongational viscosity in the fiber direction varies as the square of the fiber aspect ratio and a complex function of the fiber volume fraction. Furthermore, the ratio of the axial elongational viscosity to the transverse elongational viscosity and both axial and transverse shear viscosities was shown to be 10 exp 4 - 10 exp 6 for fiber aspect ratio of 100-1000, except at extreme values of the fiber volume fraction.
A NASTRAN analysis of the external tank (ET) substructure of the 1/8-scale space shuttle structural dynamics model is described. The NASTRAN hydroelastic procedures were used to form a model of the liquid oxygen portion of the ET. Large computer storage requirements and running times were required unless these procedures were modified. Several possibilities were demonstrated including the substitution of the real for the complex eigenvalue routine and the use of the OMIT capability to reduce the number of fluid coordinates. A NASTRAN model of the complete ET was formed and reduced to 252 degrees of freedom. A review of the eigenvector extracted, using the unsymmetrical inverse power method, indicated that the structural OMIT's resulted in unsatisfactory modal deflections. Modifications to the NASTRAN hydroelastic capability appeared necessary to generate a successful ET model. An approach is described which, by assuming incompressibility, reduces the fluid representation to a symmetric mass matrix which can be added to the structural mass. The problem is solved using faster and more efficient eigenvalue routines. Using this approach the ET NASTRAN model was analyzed for three separate weight conditions.
The refinement was investigated of a method for evaluating the pressure distribution on a body surface of arbitrary shape in incompressible flow. The solution was obtained in terms of the velocity potential, through numerical approximations which require the use of a high speed digital computer. The box method and the modal method are described in detail, and were applied to a very thin, rectangular wing in incompressible, steady flow. The box method is found to be more practical as it is applicable to more general geometries (the modal method requires a new set of functions for each geometry), and requires less computer time (fifty percent of that required by the modal method for the same problem).
Free laminar boundary layer flow between parallel streams of different magnetic fields and temperatures for incompressible, viscous, thermal, and electric conducting fluids
Vortex flow analysis of viscous, incompressible, electrically conducting fluid between two plates under applied axial magnetic field
Series solution for unbounded mixing of two incompressible homogeneous coaxial fluids with constant properties, using successive approximations method
Two second-order-closure turbulence models were devised that are suitable for predicting properties of complex turbulent flow fields in both incompressible and compressible fluids. One model is of the "two-equation" variety in which closure is accomplished by introducing an eddy viscosity which depends on both a turbulent mixing energy and a dissipation rate per unit energy, that is, a specific dissipation rate. The other model is a "Reynolds stress equation" (RSE) formulation in which all components of the Reynolds stress tensor and turbulent heat-flux vector are computed directly and are scaled by the specific dissipation rate. Computations based on these models are compared with measurements for the following flow fields: (a) low speed, high Reynolds number channel flows with plane strain or uniform shear; (b) equilibrium turbulent boundary layers with and without pressure gradients or effects of compressibility; and (c) flow over a convex surface with and without a pressure gradient.
Consistent forms for the interfacial force, the interfacial pressure, the Reynolds stresses and the particle stress have been derived for the inviscid, irrotational incompressible flow of fluid in a dilute suspension of spheres. The particles are assumed to have a velocity distribution, giving rise to an effective pressure and stress in the particle phase. The velocity fluctuations also contribute in the fluid Reynolds stress and in the (elastic) stress field inside the spheres. The relation of these constitutive equations to the force on an individual sphere is discussed.
The role of the internal gravity modes in mediating the growth of nonaxisymmetric instabilities is investigated by studying the instability of stratified incompressible differentially rotating fluid cylinders to global nonaxisymmetric modes. The results indicate that, in addition to a modified version of the well-known principal branch mediated by surface modes of the system (analogous to f-modes of stars), there exist unstable branches of the dispersion relation mediated by internal gravity modes of the system (similar to the g-modes of stars). These branches arise due to the interaction between the g-modes. It is shown that the maximum growth rate on one of the new branches can sometimes equal or exceed that on the principal branch, thus modifying the principal branch.
The use of preconditioning methods to accelerate the convergence to a steady state for both the incompressible and compressible fluid dynamic equations are considered. The relation between them for both the continuous problem and the finite difference approximation is also considered. The analysis relies on the inviscid equations. The preconditioning consists of a matrix multiplying the time derivatives. Hence, the steady state of the preconditioned system is the same as the steady state of the original system. For finite difference methods the preconditioning can change and improve the steady state solutions. An application to flow around an airfoil is presented.
Analysis of unsteady incompressible flow out of different tank configurations to determine flow rates and position of liquid levels - fluid mechanics
Current interest in microgravity materials processing has focused attention upon the finite fluid column. This configuration is used in the modeling of float zones. In this paper the incompressible inviscid finite length fluid column is subjected to an axial, time-dependent disturbance. The response of the fluid system and the resulting interface location are determined via Laplace transform methods.
The use of a magnetic field to control the motion of electrically conducting fluids is studied. The incompressible boundary-layer solutions are found for flow over a flat plate when the magnetic field is fixed relative to the plate or to the fluid. The equations are integrated numerically for the effect of the transverse magnetic field on the velocity and temperature profiles, and hence, the skin friction and rate of heat transfer. It is concluded that the skin friction and the heat-transfer rate are reduced when the transverse magnetic field is fixed relative to the plate and increased when fixed relative to the fluid. The total drag is increased in all of the areas.