Search NASA⌕ Search

Engineering topics

Hankey, W. L.

Publications and source records attributed to Hankey, W. L..

Computation of turbulent flow in a thin liquid layer of fluid involving a hydraulic jump

Numerically computed flow fields and free surface height distributions are presented for the flow of a thin layer of liquid adjacent to a solid horizontal surface that encounters a hydraulic jump. Two kinds of flow configurations are considered: two-dimensional plane flow and axisymmetric radial flow. The computations used a boundary-fitted moving grid method with a k-epsilon model for the closure of turbulence. The free surface height was determined by an optimization procedure which minimized the error in the pressure distribution on the free surface. It was also checked against an approximate procedure involving integration of the governing equations and use of the MacCormack predictor-corrector method. The computed film height also compared reasonably well with previous experiments. A region of recirculating flow was found to be present adjacent to the solid boundary near the location of the jump, which was caused by a rapid deceleration of the flow.

Rahman, M. M.↗

Analysis of the fluid flow and heat transfer in a thin liquid film in the presence and absence of gravity

The hydrodynamic and thermal behavior of a thin liquid film flowing over a solid horizontal surface is analyzed for both plane and radially spreading flows. The situations where the gravitational force is completely absent and where it is significant are analyzed separately and their practical relevance to a micro-gravity environment is discussed. In the presence of gravity, in addition to Reynolds number, the Froude number of the film is found to be an important parameter that determines the supercritical and subcritical flow regimes and any associated hydraulic jump. A closed-form solution is possible under some flow situations, whereas others require numerical integration of ordinary differential equations. The approximate analytical results are found to compare well with the available two-dimensional numerical solutions.

Rahman, M. M.↗

Prediction of heat transfer to a thin liquid film in plane and radially spreading flows

The energy equation is incorporated in the solution algorithm of Rahman et al. (1990) to compute the heat transfer to a thin film in the presence or absence of gravity. For a plane flow under zero gravity, it is found that, for both isothermal and uniformly heated walls, the heat transfer coefficient gradually decreases downstream, with Nu-asterisk (the Nusselt number in terms of film height) remaining approximately constant except for regions very close to the entrance. In the case of radial flow under zero gravity, Nu-asterisk is found to decrease monotonically when the plane is uniformly heated. Two different flow regimes are identified in the presence of gravity for both plane and radial flows. The results of the study may be applicable to the design of space-based cooling systems.

Rahman, M. M.↗

The flow of a thin liquid film on a stationary and rotating disk. II - Theoretical prediction

The existing theoretical models are improved and a systematic procedure to compute the free surface flow of a thin liquid film is suggested. The solutions for axisymmetric radial flow on a stationary horizontal disk and for the disk rotating around its axis are presented. The theoretical predictions are compared with the experimental data presented in Part I of this report. The analysis shows results for both supercritical and subcritical flows and the flow structure in the vicinity of a hydraulic jump which isolates these two flow types. The detailed flow structure in a hydraulic jump was computed and shown to contain regions of separation including a 'surface roller'. The effects of surface tension are found to be important near the outer edge of the disk where the fluid experiences a free fall. At other locations, the surface tension is negligible. For a rotating disk, the frictional resistance in the angular direction is found to be as important as that in the radial direction.

Rahman, M. M.↗

A numerical investigation of the influence of surface roughness on heat transfer in ice accretion

The flowfield and resulting heat transfer rate over a series of ice accretion shapes is obtained by solving the Navier-Stokes equations. The influence of surface roughness on surface heat transfer is examined by including blockage, form drag, and stagnation heating effects as source terms in the governing equations. The results indicate increases of a factor of three in cooling rates due to distributed roughness compared to smotth surfaces. In addition, droplet impingement efficiencies are studied for the same series of ice accretion shapes using a time-dependent solution procedure.

Scott, J. N.↗

Heat transfer to a thin liquid film with a free surface

The numerically computed flow field and heat transfer coefficient are presented for the free surface flow of a thin liquid film in the presence or absence of a gravitational body force. The results of an approximate analysis using the Pohlhausen integral method is also discussed. The flow systems studied here include both plane and radial film flows in the presence or absence of a gravitational body force. The heating conditions include isothermal and uniformly heated surfaces. The transport conditions considered at the free surface are an adiabatic condition and an evaporative free surface maintained at its saturation temperature. The height of the free surface, flow field and heat transfer coefficient, were found to be strongly affected by the gravitational body force. They were also found to depend on the Reynolds number and Froude number of the incoming fluid. In the presence of gravity, a hydraulic jump was found to occur under some flow conditions.

Rahman, M. M.↗

Navier-Stokes solutions of flowfield characteristics produced by ice accretion

The flowfield and resultant heat transfer rates over a series of ice accretion shapes have been obtained through numerical solutions of the Navier-Stokes equations. The influence of roughness is modeled by including blockage, form drag and stagnation heating effects as source terms in the governing equations. Using the flowfield information obtained from the Navier-Stokes equations the droplet impingement efficiencies are computed using a PNS-type solving scheme. Good agreement is achieved between the numerical results and experimental data.

Scott, J. N.↗

Numerical simulation of excited jet mixing layers

A numerical simulation of unsteady flow in jet mixing layers, both with and without external excitation, has been performed by solving the time-dependent compressible Navier-Stokes equations. Computations were performed on a CRAY X-MP computer using MacCormick's explicit finite difference algorithm. Different excitation methods were investigated and were shown to be very effective in controlling the well organized periodic production, shedding and pairing of large scale vortex structures. It is found that pressure excitation was generally more effective than temperature excitation, and that grid refinement results in substantial improvement in the resolution of unsteady features. The location and orientation, in addition to the frequency, of the excitation source are shown to have a significant influence on the production and interaction of large scale vortex structures in the jet mixing layer.

Scott, J. N.↗

Navier Stokes solution of the flowfield over ice accretion shapes

The numerical simulation of flow about ice accretion shapes has been accomplished by solving the Navier-Stokes equations using MacCormack's explicit finite difference scheme. The computations were performed on a CRAY-XMP computer. The influence of turbulence is taken into account by means of an algebraic eddy-viscosity model. In order to optimize the grid spacing and to achieve near orthogonality at the surface of the complex ice shapes, a hyperbolic grid generation scheme is utilized. Particular attention is given to the heat transfer process for which good agreement between the numerical and experimental results is achieved. In addition, liquid water droplet trajectories are coupled within the flowfield along with the resulting collection efficiencies using a parabolized Navier-Stokes formulation.

Scott, J. N.↗

Numerical simulation of unsteady flow in an axisymmetric shear layer

Numerical simulation of unsteady flow in an axisymmetric subsonic shear layer is accomplished by solving the time-dependent compressible Navier-Stokes equations. The objective of the effort is to investigate by numerical means the influence of various flow parameters on the shear layer behavior. The parameters investigated include the velocity ratio of two streams, total temperature, and nozzle lip thickness. The computations are performed on a CRAY-IS computer using McCormack's explicit finite difference scheme. The computed results generally show qualitative agreement with experimental data.

Scott, J. N.↗

Flow oscillations of spike-tipped bodies

The oscillating flow field around a spike-tipped body has been simulated by the numerical solution of the 3-D compressible Navier-Stokes equations at a Mach number of 3.0 and a nominal Reynolds number of 7.87 x 10 to the 6th/m. Computations were performed using a vectorized 3-D Navier-Stokes program on the STAR 100 computer. Numerical solutions confirmed the experimental result that the self-sustained oscillation occurs within a limited range of the protruded spike length to shoulder height ratio. The numerical result predicted correctly the discrete frequency range as well as the rms pressure intensity. The detailed flow structure is also presented and discussed.

Shang, J. S.↗