Search NASASearch

SEARCH · Search NASA

Results for “Gradient 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 73 records · Page 4

Transonic Navier-Stokes wing solutions using a zonal approach. Part 2: High angle-of-attack simulation

A computer code is under development whereby the thin-layer Reynolds-averaged Navier-Stokes equations are to be applied to realistic fighter-aircraft configurations. This transonic Navier-Stokes code (TNS) utilizes a zonal approach in order to treat complex geometries and satisfy in-core computer memory constraints. The zonal approach has been applied to isolated wing geometries in order to facilitate code development. Part 1 of this paper addresses the TNS finite-difference algorithm, zonal methodology, and code validation with experimental data. Part 2 of this paper addresses some numerical issues such as code robustness, efficiency, and accuracy at high angles of attack. Special free-stream-preserving metrics proved an effective way to treat H-mesh singularities over a large range of severe flow conditions, including strong leading-edge flow gradients, massive shock-induced separation, and stall. Furthermore, lift and drag coefficients have been computed for a wing up through CLmax. Numerical oil flow patterns and particle trajectories are presented both for subcritical and transonic flow. These flow simulations are rich with complex separated flow physics and demonstrate the efficiency and robustness of the zonal approach.

Chaderjian, N. M.

An experimental and numerical investigation of shock wave induced turbulent boundary layer separation at hypersonic speeds

A thoroughly documented experiment is reported that was specifically designed to test and guide computations of the interaction of an impinging shock wave with a turbulent boundary layer. Detailed mean flow field and surface data are presented for two shock strengths which resulted in attached and separated flows, respectively. Numerical computations are used to illustrate the dependence of the computations on the particulars of the turbulence models. Models appropriate for zero pressure gradient flows predicted the overall features of the flow fields, but were deficient in predicting many of the details of the interaction regions. Improvements to the turbulence model parameters were sought through a combination of detailed data analysis and computer simulations which tested the sensitivity of the solutions to model parameter changes. Computer simulations using these improvements are presented and discussed.

Marvin, J. G.

Effects of Periodic Unsteady Wake Flow and Pressure Gradient on Boundary Layer Transition Along the Concave Surface of a Curved Plate

Boundary layer transition and development on a turbomachinery blade is subjected to highly periodic unsteady turbulent flow, pressure gradient in longitudinal as well as lateral direction, and surface curvature. To study the effects of periodic unsteady wakes on the concave surface of a turbine blade, a curved plate was utilized. On the concave surface of this plate, detailed experimental investigations were carried out under zero and negative pressure gradient. The measurements were performed in an unsteady flow research facility using a rotating cascade of rods positioned upstream of the curved plate. Boundary layer measurements using a hot-wire probe were analyzed by the ensemble-averaging technique. The results presented in the temporal-spatial domain display the transition and further development of the boundary layer, specifically the ensemble-averaged velocity and turbulence intensity. As the results show, the turbulent patches generated by the wakes have different leading and trailing edge velocities and merge with the boundary layer resulting in a strong deformation and generation of a high turbulence intensity core. After the turbulent patch has totally penetrated into the boundary layer, pronounced becalmed regions were formed behind the turbulent patch and were extended far beyond the point they would occur in the corresponding undisturbed steady boundary layer.

Schobeiri, M. T.

Calculation of skin-friction coefficients for low Reynolds number turbulent boundary layer flows

An analysis is presented of the reliability of various generally accepted empirical expressions for the prediction of the skin-friction coefficient C/sub f/ of turbulent boundary layers at low Reynolds numbers in zero-pressure-gradient flows on a smooth flat plate. The skin-friction coefficients predicted from these expressions were compared to the skin-friction coefficients of experimental profiles that were determined from a graphical method formulated from the law of the wall. These expressions are found to predict values that are consistently different than those obtained from the graphical method over the range 600 Re/sub theta 2000. A curve-fitted empirical relationship was developed from the present data and yields a better estimated value of C/sub f/ in this range. The data, covering the range 200 Re/sub theta 7000, provide insight into the nature of transitional flows. They show that fully developed turbulent boundary layers occur at Reynolds numbers Re/sub theta/ down to 425. Below this level there appears to be a well-ordered evolutionary process from the laminar to the turbulent profiles. These profiles clearly display the development of the turbulent core region and the shrinking of the laminar sublayer with increasing values of Re/sub theta/.

Barr, P. K.

Unstructured and adaptive mesh generation for high Reynolds number viscous flows

A method for generating and adaptively refining a highly stretched unstructured mesh suitable for the computation of high-Reynolds-number viscous flows about arbitrary two-dimensional geometries was developed. The method is based on the Delaunay triangulation of a predetermined set of points and employs a local mapping in order to achieve the high stretching rates required in the boundary-layer and wake regions. The initial mesh-point distribution is determined in a geometry-adaptive manner which clusters points in regions of high curvature and sharp corners. Adaptive mesh refinement is achieved by adding new points in regions of large flow gradients, and locally retriangulating; thus, obviating the need for global mesh regeneration. Initial and adapted meshes about complex multi-element airfoil geometries are shown and compressible flow solutions are computed on these meshes.

Mavriplis, Dimitri J.

Unstructured and adaptive mesh generation for high Reynolds number viscous flows

A method for generating and adaptively refining a highly stretched unstructured mesh suitable for the computation of high-Reynolds-number viscous flows about arbitrary two-dimensional geometries was developed. The method is based on the Delaunay triangulation of a predetermined set of points and employs a local mapping in order to achieve the high stretching rates required in the boundary-layer and wake regions. The initial mesh-point distribution is determined in a geometry-adaptive manner which clusters points in regions of high curvature and sharp corners. Adaptive mesh refinement is achieved by adding new points in regions of large flow gradients, and locally retriangulating; thus, obviating the need for global mesh regeneration. Initial and adapted meshes about complex multi-element airfoil geometries are shown and compressible flow solutions are computed on these meshes.

Mavriplis, D. J.

A Grid-resolved Analysis of Base Flowfield for a Four-Engine Clustered Nozzle Configuration

The objective of this study is to propose a computational methodology that can effectively anchor the base flowfield of a four-engine clustered nozzle configuration. This computational methodology is based on a three-dimensional, viscous flow, pressure-based computational fluid dynamics (CFD) formulation. For efficient CFD calculation, a Prandtl-Meyer solution treatment is applied to the algebraic grid lines for initial plume expansion resolution. As the solution evolves, the computational grid is adapted to the pertinent flow gradients. The CFD model employs an upwind scheme in which second- and fourth-order central differencing schemes with artificial dissipation are used. The computed quantitative base flow properties such as the radial base pressure distributions, model centerline static pressure, Mach number and impact pressure variations, and base pressure characteristic curve agreed reasonably well with those of the measurement.

Wang, Ten-See

Applications of velocity potential function to acoustic duct propagation and radiation from inlets using finite element theory

A finite element velocity potential program was developed to study acoustic wave propagation in complex geometries. For irrotational flows, relatively low sound frequencies, and plane wave input, the finite element solutions showed significant effects of inlet curvature and flow gradients on the attenuation of a given acoustic liner in a realistic variable area turbofan inlet. The velocity potential approach can not be used to estimate the effects of rotational flow on acoustic propagation, since the potential acoustic disturbances propagate at the speed of the media in sheared flow. Approaches are discussed that are being considered for extending the finite element solution to include the far field, as well as the internal portion of the duct. A new matrix partitioning approach is presented that can be incorporated in previously developed programs to allow the finite element calculation to be marched into the far field. The partitioning approach provided a large reduction in computer storage and running times.

Baumeister, K. J.

SAGE: The Self-Adaptive Grid Code

The multi-dimensional self-adaptive grid code, SAGE, is an important tool in the field of computational fluid dynamics (CFD). It provides an efficient method to improve the accuracy of flow solutions while simultaneously reducing computer processing time. Briefly, SAGE enhances an initial computational grid by redistributing the mesh points into more appropriate locations. The movement of these points is driven by an equal-error-distribution algorithm that utilizes the relationship between high flow gradients and excessive solution errors. The method also provides a balance between clustering points in the high gradient regions and maintaining the smoothness and continuity of the adapted grid, The latest version, Version 3, includes the ability to change the boundaries of a given grid to more efficiently enclose flow structures and provides alternative redistribution algorithms.

Davies, Carol B.

Numerical simulation of aerothermal loads in hypersonic engine inlets due to shock impingement

The effect of shock impingement on an axial corner simulating the inlet of a hypersonic vehicle engine is modeled using a finite-difference procedure. A three-dimensional dynamic grid adaptation procedure is utilized to move the grids to regions with strong flow gradients. The adaptation procedure uses a grid relocation stencil that is valid at both the interior and boundary points of the finite-difference grid. A linear combination of spatial derivatives of specific flow variables, calculated with finite-element interpolation functions, are used as adaptation measures. This computational procedure is used to study laminar and turbulent Mach 6 flows in the axial corner. The description of flow physics and qualitative measures of heat transfer distributions on cowl and strut surfaces obtained from the analysis are compared with experimental observations. Conclusions are drawn regarding the capability of the numerical scheme for enhanced modeling of high-speed compressible flows.

Ramakrishnan, R.

Dynamics Explorer measurements of particles, fields, and plasma drifts over a horse-collar auroral pattern

As shown from ground-based measurements and satellite-borne imagers, one type of global auroral pattern characteristic of quiet (usually northward IMF) intervals is that of a contracted but thickened emission region of a pattern referred to as 'horse-collar' aurora (Hones et al., 1989). In this report we use the Dynamics Explorer data set to examine a case in which this horse-collar pattern was observed by the DE-1 auroral imager, while at the same time DE-2, at lower altitude, measured precipitating particles, electric and magnetic fields, and plasma drifts. Our analysis shows that, in general, there is close agreement between the optical signatures and the particle precipitation patterns. In many instances, over scales ranging from tens to a few hundred kilometers, electron precipitation features and upward field-aligned currents are observed at locations where the plasma flow gradients indicate negative V-average x E. The particle, plasma, and field measurements made along the satellite track and the 2D perspective of the imager provide a means of determining the configuration of convective flows in the high-latitude ionosphere during this interval of northward IMF. Recent mapping studies are used to relate the low-altitude observations to possible magnetospheric source regions.

Sharber, J. R.

Refraction of high frequency noise in an arbitrary jet flow

Refraction of high frequency noise by mean flow gradients in a jet is studied using the ray-tracing methods of geometrical acoustics. Both the two-dimensional (2D) and three-dimensional (3D) formulations are considered. In the former case, the mean flow is assumed parallel and the governing propagation equations are described by a system of four first order ordinary differential equations. The 3D formulation, on the other hand, accounts for the jet spreading as well as the axial flow development. In this case, a system of six first order differential equations are solved to trace a ray from its source location to an observer in the far field. For subsonic jets with a small spreading angle both methods lead to similar results outside the zone of silence. However, with increasing jet speed the two prediction models diverge to the point where the parallel flow assumption is no longer justified. The Doppler factor of supersonic jets as influenced by the refraction effects is discussed and compared with the conventional modified Doppler factor.

Khavaran, Abbas

Applications of velocity potential function to acoustic duct propagation and radiation from inlets using finite element theory

A finite element velocity potential program has been developed to study acoustic wave propagation in complex geometries. For irrotational flows, relatively low sound frequencies, and plane wave input, the finite element solutions show significant effects of inlet curvature and flow gradients on the attenuation of a given acoustic liner in a realistic variable area turbofan inlet. In addition, as shown in the paper, the velocity potential approach can not be used to estimate the effects of rotational flow on acoustic propagation since the potential acoustic disturbances propagate at the speed of the media in sheared flow. Approaches are discussed that are being considered for extending the finite element solution to include the far field as well as the internal portion of the duct. A new matrix partitioning approach is presented that can be incorporated in previously developed programs to allow the finite element calculation to be marched into the far field. The partitioning approach provides a large reduction in computer storage and running times.

Baumeister, K. J.

Calculation of convective heat transfer on highly blunt bodies at flow incidence

An implicit finite-difference code is used to study three-dimensional viscous heat-conducting flows over the forebody of hypersonic vehicles. In the method, adaptive grids are generated to the shock and body contour, and local flow gradients and total enthalpy are used to control numerical dissipation. Typical axisymmetric configurations of a sphere, ellipsoid, and flat-face disk are considered, along with the cases of a blunt 70-deg cone at 0 and 20 deg flow incidences and an asymmetric 60-deg cone raked off at a 73-deg angle.

Li, C. P.

Flow prediction over a transport multi-element high-lift system and comparison with flight measurements

Accurate prediction of surface-pressure distributions, merging boundary-layers, and separated-flow regions over multi-element high-lift airfoils is required to design advanced high-lift systems for efficient subsonic transport aircraft. The availability of detailed measurements of pressure distributions and both averaged and time-dependent boundary-layer flow parameters at flight Reynolds numbers is critical to evaluate computational methods and to model the turbulence structure for closure of the flow equations. Several detailed wind-tunnel measurements at subscale Reynolds numbers were conducted to obtain detailed flow information including the Reynolds-stress component. As part of a subsonic-transport high-lift research program, flight experiments are conducted using the NASA-Langley B737-100 research aircraft to obtain detailed flow characteristics for support of computational and wind-tunnel efforts. Planned flight measurements include pressure distributions at several spanwise locations, boundary-layer transition and separation locations, surface skin friction, as well as boundary-layer profiles and Reynolds stresses in adverse pressure-gradient flow.

Vijgen, P. M. H. W.

Recent Turbulence Model Advances Applied to Multielement Airfoil Computations

A one-equation linear turbulence model and a two-equation nonlinear explicit algebraic stress model (EASM) are applied to the flow over a multielement airfoil. The effect of the K-epsilon and K-omega forms of the two-equation model are explored, and the K-epsilon form is shown to be deficient in the wall-bounded regions of adverse pressure gradient flows. A new K-omega form of EASM is introduced. Nonlinear terms present in EASM are shown to improve predictions of turbulent shear stress behind the trailing edge of the main element and near midflap. Curvature corrections are applied to both the one- and two-equation turbulence models and yield only relatively small local differences in the flap region, where the flow field undergoes the greatest curvature. Predictions of maximum lift are essentially unaffected by the turbulence model variations studied.

Rumsey, Christopher L.

An experimental and numerical investigation of shock-wave induced turbulent boundary-layer separation at hypersonic speeds

An experiment designed to test and guide computations of the interaction of an impinging shock wave with a turbulent boundary layer is described. Detailed mean flow-field and surface data are presented for two shock strengths which resulted in attached and separated flows, respectively. Numerical computations, employing the complete time-averaged Navier-Stokes equations along with algebraic eddy-viscosity and turbulent Prandtl number models to describe shear stress and heat flux, are used to illustrate the dependence of the computations on the particulars of the turbulence models. Models appropriate for zero-pressure-gradient flows predicted the overall features of the flow fields, but were deficient in predicting many of the details of the interaction regions. Improvements to the turbulence model parameters were sought through a combination of detailed data analysis and computer simulations which tested the sensitivity of the solutions to model parameter changes. Computer simulations using these improvements are presented and discussed.

Marvin, J. G.

A prediction method for velocity and temperature profiles in a two-dimensional nominally steady turbulent boundary layer

This paper describes a recently developed boundary-layer prediction method for a variable property compressible flow, in which heat transfer takes place primarily by forced convection and for which the mainstream Mach number is small. The leading order terms, in asymptotic expansions for large Reynolds numbers, are obtained for the mean velocity and temperature distribution in both the inner and outer layer of the turbulent boundary layer. Closure in the inner layer is achieved using an analytical model for the mean profiles which is based on the observed coherent structure of the time-dependent inner layer flow. For the outer layer, simple eddy viscosity and conductivity models are developed without recourse to the Reynolds analogy. In the prediction method a numerical solution of the outer layer equations is matched to the analytical inner layer profiles as the computation procedes downstream. Calculations are presented for a range of adverse and favorable pressure gradient flows and the predicted results compare well with existing data.

Weigand, G. G.