Search NASA⌕ Search

SEARCH · Search NASA

Results for “CHANNEL 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 253 records · Page 14

Computations of internal turbulent flow with large separated flow regions

An implicit two-equation turbulence solver in generalized coordinates has been developed and is used in conjunction with the three-dimensional incompressible Navier-Stokes solver, INS3D, to calculate the internal flow in one channel and in an additional channel with a sudden 2:3 expansion. A new and consistent boundary procedure for a low Reynolds number form of the kappa-epsilon turbulence model is chosen to integrate the equations up to the wall. The high Reynolds number form of the equations are integrated using wall functions. The latter approach yields a faster convergence to the steady state solution than the former. For the case of channel flow, both the wall function and wall boundary condition approaches yield results in good agreement with the experimental data. the back-step (sudden expansion) flow is calculated using wall function approach. The predictions are in reasonable agreement with the experimental data.

Kaul, U. K.↗

Computational simulation of transition to turbulence through inverse modeling

The present investigation has focused on a computational methodology for the fundamental case of transition in channel flow, in which recently published experimental data are utilized both as a stimulus and as a measure of merit of the method. The research has proceeded along three avenues in parallel. The first task has consisted of the development and verification of a computer code which calculates the mean evolution of flow in a channel similar to the one employed experimentally by Blair and Anderson. An analytical test case was created for the dual purposes of code verification and of highlighting the interactions between the Reynolds stress and the mean velocity profile. This test case generated a Reynolds stress by the residue in the momentum equation which is produced by a typical analytical velocity profile. By a substitution of this Reynolds stress into the appropriate code module, the correctness of the code may be verified, along with the accuracy of the computational method. The second task pursued has involved the development of a triple layer model for the Reynolds stress profile, which was suggested and derived from experimental velocity profiles. It is demonstrated that the innermost length scale is based on the local friction velocity, the intermediate layer corresponds to the usual logarithmic law of the wall region in which the normalized Reynolds stress is approximately unity, and the outermost layer is represented by a closed mathematical form depending explicitly on the velocity profile in the wake region. The third task was comprised of scrutiny of the excellent databases developed by Blair and others, and the planning of its incorporation into the transition analysis. These extensive measurements indicate that turbulent statistics in the transition regime may be considered to alternate between laminar and fully turbulent types, the proportions of which are quantified by a measured intermittency function.

Sepri, Paavo↗

A dynamic subgrid-scale model for compressible turbulence and scalar transport

The dynamic subgrid-scale (SGS) model of Germano et al. (1991) is generalized for the large eddy simulation (LES) of compressible flows and transport of a scalar. The model was applied to the LES of decaying isotropic turbulence, and the results are in excellent agreement with experimental data and direct numerical simulations. The expression for the SGS turbulent Prandtl number was evaluated using direct numerical simulation (DNS) data in isotropic turbulence, homogeneous shear flow, and turbulent channel flow. The qualitative behavior of the model for turbulent Prandtl number and its dependence on molecular Prandtl number, direction of scalar gradient, and distance from the wall are in accordance with the total turbulent Prandtl number from the DNS data.

Moin, P.↗

Evaporative Cooling and Dehumidification Garment for Portable Life Support Systems

This paper describes the design and development of an innovative thermal and humidity control system for future space suits. The system comprises an evaporation cooling and dehumidification garment (ECDG) and a lithium chloride absorber radiator (LCAR). The ECDG absorbs heat and water vapor from inside the suit pressure garment, while the LCAR rejects heat to space without venting water vapor. The ECDG is built from thin, flexible patches with coversheets made of non-porous, water-permeable membranes that -enclose arrays of vapor flow passages. Water vapor from inside the spacesuit diffuses across the water permeable membranes, enters the vapor flow channels, and then flows to the LCAR, thus dehumidifying the internal volume of the space suit pressure garment. Additional water evaporation inside the ECDG provides cooling for sensible heat loads. -The heat released from condensation and absorption in the LCAR is rejected to the environment by thermal radiation. We have assembled lightweight and flexible ECDG pouches from prototypical materials and measured their performance in a series of separate effects tests under well-controlled, prototypical conditions. Sweating hot plate tests at typical space suit pressures show that ECDG pouches can absorb over 60 W/ft of latent heat and 20 W/ft of sensible heat from the pressure garment environment. These results are in good agreement with the predictions of our analysis models.

Izenson, Michael↗

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗

Continuous wavelet analysis of coherent structures

We perform an analysis of planar cuts through three dimensional turbulent fields (planar channel flow and mixing layer) using the 2D continuous wavelet transform. We propose two new diagnostics: (1) a measure of intermittency I(r, vector x), which is the ratio of local energy and average energy at a given scale r; and (2) a local Reynolds number, defined on the local velocity contribution at a given scale, computed from the wavelet transform of the three velocity components, the scale of the transform, and molecular viscosity; this gives a representation of the local non-linearity of the flow viewed in both space and scale. We find, for the analyzed flows, strong small-scale intermittency located in the ejection regions for the channel flow and in the vortex core of the mixing layer.

Farge, M.↗

Turbulent flow in a channel with transverse rib heat transfer augmentation

Turbulent flow in a 2D channel with repeated rectangular rib roughness was numerically simulated using a low Reynolds number form of the k-epsilon turbulence model. Friction factors and average Stanton numbers were calculated for various pitch to rib height ratios and bulk Reynolds numbers. Comparisons with experiment were generally adequate, with the predictions of friction superior to those for heat transfer. The effect of variable properties for channel flow was investigated, and the results showed a greater effect for friction than for heat transfer. Comparison with experiment yielded no clear conclusions. The turbulence model was also validated for a related problem, that of flow downstream of an abrupt pipe expansion.

Chang, B. H.↗

On the measurement of lateral velocity derivatives in turbulent flows

Direct numerical simulation data for the lateral velocity derivative delta(u)/delta(y) at the centerline of a fully developed turbulent channel flow provide reasonable support for Wyngaard's analysis of the error involved in measuring this quantity using parallel hot wires. Numerical data in the wall region of the channel flow also provide a useful indication of how to select the separation between the wires. Justification for this choice is obtained by comparing several measured statistics of delta(u)/delta(y) with the corresponding numerical data.

Antonia, R. A.↗

Turbulent flow in a 180 deg bend: Modeling and computations

A low Reynolds number k-epsilon turbulence model was presented which yields accurate predictions of the kinetic energy near the wall. The model is validated with the experimental channel flow data of Kreplin and Eckelmann. The predictions are also compared with earlier results from direct simulation of turbulent channel flow. The model is especially useful for internal flows where the inflow boundary condition of epsilon is not easily prescribed. The model partly derives from some observations based on earlier direct simulation results of near-wall turbulence. The low Reynolds number turbulence model together with an existing curvature correction appropriate to spinning cylinder flows was used to simulate the flow in a U-bend with the same radius of curvature as the Space Shuttle Main Engine (SSME) Turn-Around Duct (TAD). The present computations indicate a space varying curvature correction parameter as opposed to a constant parameter as used in the spinning cylinder flows. Comparison with limited available experimental data is made. The comparison is favorable, but detailed experimental data is needed to further improve the curvature model.

Kaul, Upender K.↗

Second-order closure models for rotating turbulent flows

The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier-Stokes equations for the problem of fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler k-epsilon and k-l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.

Speziale, Charles G.↗

Second-order closure models for rotating turbulent flows

The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier-Stokes equations for the problem of a fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler K-epsilon and K-l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.

Speziale, C. G.↗

Predicting NonInertial Effects with Algebraic Stress Models which Account for Dissipation Rate Anisotropies

Three types of turbulence models which account for rotational effects in noninertial frames of reference are evaluated for the case of incompressible, fully developed rotating turbulent channel flow. The different types of models are a Coriolis-modified eddy-viscosity model, a realizable algebraic stress model, and an algebraic stress model which accounts for dissipation rate anisotropies. A direct numerical simulation of a rotating channel flow is used for the turbulent model validation. This simulation differs from previous studies in that significantly higher rotation numbers are investigated. Flows at these higher rotation numbers are characterized by a relaminarization on the cyclonic or suction side of the channel, and a linear velocity profile on the anticyclonic or pressure side of the channel. The predictive performance of the three types of models are examined in detail, and formulation deficiencies are identified which cause poor predictive performance for some of the models. Criteria are identified which allow for accurate prediction of such flows by algebraic stress models and their corresponding Reynolds stress formulations.

Jongen, T.↗

Advanced Manifolds for Improved Solid Oxide Electrolyzer Performance

An investigation was conducted to see if additive manufacturing could be used to fabricate more efficient manifold designs for improved flow, reduced stresses, and decreased number of joints to be sealed for a solid oxide electrolyzer used to convert carbon dioxide to oxygen. Computational flow and mechanical modeling were conducted on a NASA Glenn Research Center patented cell and stack design with the potential to achieve a 3-4 times mass reduction. Various manifold designs were modeled, and two were downselected to be fabricated and tested. Some benefit was seen in a baffled manifold design, which directed incoming flow more effectively into the flow channels, compared to the original design, where the flow spent more time within the manifold itself. Flow measurements indicated some non-uniformity of flow across the channels at higher flow rates, which were not predicted by the model. Some possible explanations for the differences are discussed.

Linne, Diane L.↗

An experimental investigation of shock wave-turbulent boundary layer interactions with and without boundary layer suction: A data summary report

Tabulated data from a series of experimental studies of the interaction of a shock wave with a turbulent boundary layer in axisymmetric flow configurations is presented. The studies were conducted at the walls of circular wind tunnels and on the cylindrical centerbody of an annular flow channel. Detailed pitot pressure profiles and wall static pressure profiles upstream of, within and downstream of the interaction region are given. Results are presented for flows at nominal freestream Mach Numbers of 2, 3 and 4. For studies at the tunnel sidewalls, the shock waves were produced by conical shock generators mounted on the centerline of the wind tunnel at zero angle of attack. The annular ring generator was used to produce the shock wave at the centerbody of the annular flow channel. The effects of boundary layer bleed were examined in the investigation. Both bleed rate and bleed location were studied. Most of the bleed studies were conducted with bleed holes drilled normal to the wall surface but the effects of slot suction were also examined. A summary of the principal results and conclusions is given.

Sun, C. C.↗

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

Scott, James R.↗

Boiling incipience and convective boiling of neon and nitrogen

Forced convection and subcooled boiling heat transfer data for liquid nitrogen and liquid neon were obtained in support of a design study for a 30 tesla cryomagnet cooled by forced convection of liquid neon. The cryogen data obtained over a range of system pressures, fluid flow rates, and applied heat fluxes were used to develop correlations for predicting boiling incipience and convective boiling heat transfer coefficients in uniformly heated flow channels. The accuracy of the correlating equations was then evaluated. A technique was also developed to calculate the position of boiling incipience in a uniformly heated flow channel. Comparisons made with the experimental data showed a prediction accuracy of + or - 15 percent.

Papell, S. S.↗

Reynolds stress closure modeling in wall-bounded flows

This report describes two projects. Firstly, a Reynolds stress closure for near-wall turbulence is described. It was motivated by the simpler k-epsilon-(v-bar(exp 2)) model described in last year's annual research brief. Direct Numerical Simulation of three-dimensional channel flow shows a curious decrease of the turbulent kinetic energy. The second topic of this report is a model which reproduces this effect. That model is described and used to discuss the relevance of the three dimensional channel flow simulation to swept wing boundary layers.

Durbin, Paul A.↗