Search NASA⌕ Search

SEARCH · Search NASA

Results for “wall modeling”

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 217 records · Page 12

An Improved K-Epsilon Model for Near-Wall Turbulence and Comparison with Direct Numerical Simulation

An improved k-epsilon model for low Reynolds number turbulence near a wall is presented. The near-wall asymptotic behavior of the eddy viscosity and the pressure transport term in the turbulent kinetic energy equation is analyzed. Based on this analysis, a modified eddy viscosity model, having correct near-wall behavior, is suggested, and a model for the pressure transport term in the k-equation is proposed. In addition, a modeled dissipation rate equation is reformulated. Fully developed channel flows were used for model testing. The calculations using various k-epsilon models are compared with direct numerical simulations. The results show that the present k-epsilon model performs well in predicting the behavior of near-wall turbulence. Significant improvement over previous k-epsilon models is obtained.

Shih, T. H.↗

On modelling the pressure-strain correlations in wall bounded flows

Turbulence models for the pressure-strain term of the Reynolds-stress equations in the vicinity of a moving wall are evaluated for a high Reynolds number flow using decaying grid turbulence as a model problem. The data of Thomas and Hancock are used as a base for evaluating the different turbulence models. In particular, the Rotta model for return-to-isotropy is evaluated both in its inclusion into the Reynolds-stress equation model and in comparison to a nonlinear model advanced by Sarkar and Speziale. Further, models for the wall correction to the transfer term advanced by Launder et al., Shir, and Shih and Lumley are compared. Initial data using the decaying grid turbulence experiment as a base suggests that the coefficients proposed for these models are high perhaps by as much as an order of magnitude. The Shih and Lumley model which satisfies realizability constraints, in particular, seems to hold promise in adequately modeling the Reynolds stress components of this flow. Extensions of this work are to include testing the homogeneous transfer model by Shih and Lumley and the testing of the wall transfer models using their proposed coefficients and the coefficients chosen from this work in a flow with mean shear component.

Peltier, L. J.↗

A One-Dimensional Global-Scaling Erosive Burning Model Informed by Blowing Wall Turbulence

A derivation of turbulent flow parameters, combined with data from erosive burning test motors and blowing wall tests results in erosive burning model candidates useful in one-dimensional internal ballistics analysis capable of scaling across wide ranges of motor size. The real-time burn rate data comes from three test campaigns of subscale segmented solid rocket motors tested at two facilities. The flow theory admits the important effect of the blowing wall on the turbulent friction coefficient by using blowing wall data to determine the blowing wall friction coefficient. The erosive burning behavior of full-scale motors is now predicted more closely than with other recent models.

Kibbey, Timothy P.↗

Computation of wind tunnel wall effects for complex models using a low-order panel method

A technique for determining wind tunnel wall effects for complex models using the low-order, three dimensional panel method PMARC (Panel Method Ames Research Center) has been developed. Initial validation of the technique was performed using lift-coefficient data in the linear lift range from tests of a large-scale STOVL fighter model in the National Full-Scale Aerodynamics Complex (NFAC) facility. The data from these tests served as an ideal database for validating the technique because the same model was tested in two wind tunnel test sections with widely different dimensions. The lift-coefficient data obtained for the same model configuration in the two test sections were different, indicating a significant influence of the presence of the tunnel walls and mounting hardware on the lift coefficient in at least one of the two test sections. The wind tunnel wall effects were computed using PMARC and then subtracted from the measured data to yield corrected lift-coefficient versus angle-of-attack curves. The corrected lift-coefficient curves from the two wind tunnel test sections matched very well. Detailed pressure distributions computed by PMARC on the wing lower surface helped identify the source of large strut interference effects in one of the wind tunnel test sections. Extension of the technique to analysis of wind tunnel wall effects on the lift coefficient in the nonlinear lift range and on drag coefficient will require the addition of boundary-layer and separated-flow models to PMARC.

Ashby, Dale L.↗

Investigation of turbulent flow in highly curved ducts with application to turbomachinery components

Numerical predictions were performed using a semi-elliptic calculation procedure for the case of turbulent flow in passage through a 90 deg bend of square cross section. Two versions of the isotropic turbulent viscosity two equation k-epsilon model were used. The first, the wall function model (WFM), employs the logarithmic law-of-the-wall relation and the notion of equilibrium flow to set all the necessary boundary conditions at the first grid point adjacent to a solid wall. The second model, the Van Driest model (VDM), employs Prandtl's original mixing length formulation, in conjunction with Van Driest's semi-empirical relation for the mixing length, to calculate the turbulent viscosity in the near wall regions of the flow. In this case, boundary conditions for k and epsilon, required to calculate these quantities in the core of the flow, are obtained by matching the mixing length and Reynolds number model formulations in an overlapping region of the flow near the walls. In both cases the results obtained show an improvement over earlier calculations using an elliptic numerical procedure. This is attributed to the finer grids possible in the present work. Of the two models, the VDM formulation shows better overall conformity with the mean flow measurements. Neither model reproduces well the details of the stress distribution as a result of the implied isotropic turbulent viscosity.

Humphrey, J. A. C.↗

Electroless-plating technique for fabricating thin-wall convective heat-transfer models

A technique for fabricating uniform thin-wall metallic heat-transfer models and which simulates a Shuttle thermal protection system tile is described. Two 6- by 6- by 2.5-in. tiles were fabricated to obtain local heat transfer rates. The fabrication process is not limited to any particular geometry and results in a seamless thin-wall heat-transfer model which uses a one-wire thermocouple to obtain local cold-wall heat-transfer rates. The tile is relatively fragile because of the brittle nature of the material and the structural weakness of the flat-sided configuration; however, a method was developed and used for repairing a cracked tile.

Avery, D. E.↗

Determination of equivalent model geometry for tunnel wall interference assessment/correction

A formula for the determination of equivalent model geometry with two variables measured at the interface is derived, based on two dimensional subsonic flow. This predicted model profile is a reasonable initial estimate for transonic flow as long as the sonic region does not reach the interface. A general formula is given in two forms. One is in terms of complex variable functions and the other is an integral equation. The complex-function formula has the advantage of using analytic expressions. The integral equation form requires a numerical solution after assuming the model geometry as a polynomial function. Examples are given to illustrate the application of the formulas.

Lo, C. F.↗

A near-wall four-equation turbulence model for compressible boundary layers

A near-wall four-equation turbulence model is developed for the calculation of high-speed compressible turbulent boundary layers. The four equations used are the k-epsilon equations and the theta(exp 2)-epsilon(sub theta) equations. These equations are used to define the turbulent diffusivities for momentum and heat fluxes, thus allowing the assumption of dynamic similarity between momentum and heat transport to be relaxed. The Favre-averaged equations of motion are solved in conjunction with the four transport equations. Calculations are compared with measurements and with another model's predictions where the assumption of the constant turbulent Prandtl number is invoked. Compressible flat plate turbulent boundary layers with both adiabatic and constant temperature wall boundary conditions are considered. Results for the range of low Mach numbers and temperature ratios investigated are essentially the same as those obtained using an identical near-wall k-epsilon model. In general, the numerical predictions are in very good agreement with measurements and there are significant improvements in the predictions of mean flow properties at high Mach numbers.

Sommer, T. P.↗

A Wall-Distance Method for Turbulence Modeling

The distance from a grid point to the closest wall surface, wall distance, is a funda- mental quantity in turbulence modeling. Efficiency of wall-distance calculations has become more critical as the size of computational grids has significantly increased in recent years. This paper reports on an initial implementation of a new search-based wall-distance method that is suitable for general unstructured computational fluid dynamics (CFD) grids and tailored for requirements specific for turbulence modeling. The method represents a two-step approach to calculate the wall distance. In the first step, the wall distance is approximated for each grid point as the minimum distance from this point to a vertex of a triangular face at the wall. The point-to-vertex distance calculation is relatively inexpensive but may lead to a significant error in the wall-distance ap- proximation, especially for grid points near the wall. In the second step, for grid points located within a predefined distance ( threshold ) from the wall, the wall distance is computed as the minimum distance to wall faces. As a result, the wall distance is exact for all grid points within the threshold. This two-step approach reduces the computational cost yet achieves high and controllable accuracy in the evaluation of the wall distance. Algorithmic enhancements are presented to improve efficiency of wall-distance computations. Comprehensive assessment of the new method is reported for large-scale unstructured CFD grids generated for the Fifth AIAA CFD High-Lift Prediction Workshop. The performance of the new wall-distance method compares favorably with performance of two established methods implemented in high-performance CFD codes.

Wall Distance↗

Stagnation-point heat-transfer rate predictions at aeroassist flight conditions

The results are presented for the stagnation-point heat-transfer rates used in the design process of the Aeroassist Flight Experiment (AFE) vehicle over its entire aeropass trajectory. The prediction methods used in this investigation demonstrate the application of computational fluid dynamics (CFD) techniques to a wide range of flight conditions and their usefulness in a design process. The heating rates were computed by a viscous-shock-layer (VSL) code at the lower altitudes and by a Navier-Stokes (N-S) code for the higher altitude cases. For both methods, finite-rate chemically reacting gas was considered, and a temperature-dependent wall-catalysis model was used. The wall temperature for each case was assumed to be radiative equilibrium temperature, based on total heating. The radiative heating was estimated by using a correlation equation. Wall slip was included in the N-S calculation method, and this method implicitly accounts for shock slip. The N-S/VSL combination of projection methods was established by comparison with the published benchmark flow-field code LAURA results at lower altitudes, and the direct simulation Monte Carlo results at higher altitude cases. To obtain the design heating rate over the entire forward face of the vehicle, a boundary-layer method (BLIMP code) that employs reacting chemistry and surface catalysis was used. The ratio of the VSL or N-S method prediction to that obtained from the boundary-layer method code at the stagnation point is used to define an adjustment factor, which accounts for the errors involved in using the boundary-layer method.

Gupta, Roop N.↗

Kolmogorov Behavior of Near-Wall Turbulence and Its Application in Turbulence Modeling

The near-wall behavior of turbulence is re-examined in a way different from that proposed by Hanjalic and Launder and followers. It is shown that at a certain distance from the wall, all energetic large eddies will reduce to Kolmogorov eddies (the smallest eddies in turbulence). All the important wall parameters, such as friction velocity, viscous length scale, and mean strain rate at the wall, are characterized by Kolmogorov microscales. According to this Kolmogorov behavior of near-wall turbulence, the turbulence quantities, such as turbulent kinetic energy, dissipation rate, etc. at the location where the large eddies become Kolmogorov eddies, can be estimated by using both direct numerical simulation (DNS) data and asymptotic analysis of near-wall turbulence. This information will provide useful boundary conditions for the turbulent transport equations. As an example, the concept is incorporated in the standard k-epsilon model which is then applied to channel and boundary flows. Using appropriate boundary conditions (based on Kolmogorov behavior of near-wall turbulence), there is no need for any wall-modification to the k-epsilon equations (including model constants). Results compare very well with the DNS and experimental data.

Shih, Tsan-Hsing↗

Experimental and theoretical investigation of a slotted noise shield model for wind tunnel walls

Tests at Mach 6 and 8 have been conducted on a model of a wind tunnel sound radiation shield which consists of interchangeable circular and wedge-shaped longitudinal rods with gaps between the rods for boundary layer removal. In view of the large effective sweep angle of the circular rods with respect to the local flow, the boundary layer behavior was surprisingly well represented by swept cylinder theory. Laminar flow was maintained up to a maximum local unit Reynolds number of 8 million per foot which would presumably be independent of model length for a given rod diameter by analogy with results on swept cylinders. Reduction in wind-tunnel sound levels within the shielded region of the circular rod model were about 90% of the maximum possible reduction for an ideal shield.

Harvey, W. D.↗

A One-Dimensional Global-Scaling Erosive Burning Model Informed by Blowing Wall Turbulence

This paper makes no attempt to comprehensively review erosive burning models or the data collected in pursuit of them; the interested reader could begin with Landsbaum for a historical summary. However, a discussion and comparison to recent work by McDonald and Rettenmaier and Heister will be included, along with data generated by Strand, et. al. Suffice it to say that the search for a way to predict erosive burning in any size motor with formulas cleanly applicable to a typical 1D ballistics analysis has been long thwarted. Some models were based on testing that failed to adequately simulate the solid rocket motor environment. In most cases, no real‐time burn rate measurement was available. Two popular models, even when calibrated to recent motor‐like real‐time burn rate data obtained by Furfaro, were shown by McMillin to be inadequate at modeling erosive burning in the Space Shuttle Reusable Solid Rocket Motor (RSRM), the Space Launch Systems' Five‐Segment RSRM (RSRMV), and the five‐segment Engineering Test Motor (ETM)‐3. Subsequently to the data cited from Strand and Furfaro, additional motors of the same kind as Furfaro's were tested with RSRMV propellant, utilizing 7 segments per motor and 3 throat sizes. By measuring propellant web thickness with ultrasonic gages, the burn rate was determined at cross‐flow Mach numbers up to Mach 0.8. Furthermore, because of the different throat sizes in otherwise identical motors, this provides a unique look at the effect of pressure and base burn rate on the erosive response. Figure 1 shows example of the data pertaining to the high Mach motor, where the port area is initially less than the throat area. The burn rate data was processed using a smoothing method developed to reduce the noise without too severely introducing end effects that limit the range of useful data. Then, an empirical ballistics scheme was used to estimate the flow condition based on the burn rate measurements and pressure measured between each segment.

Kibbey, Timothy P.↗

Use of the Ames Check Standard Model for the Validation of Wall Interference Corrections

The new check standard model of the NASA Ames 11-ft Transonic Wind Tunnel was chosen for a future validation of the facility's wall interference correction system. The chosen validation approach takes advantage of the fact that test conditions experienced by a large model in the slotted part of the tunnel's test section will change significantly if a subset of the slots is temporarily sealed. Therefore, the model's aerodynamic coefficients have to be recorded, corrected, and compared for two different test section configurations in order to perform the validation. Test section configurations with highly accurate Mach number and dynamic pressure calibrations were selected for the validation. First, the model is tested with all test section slots in open configuration while keeping the model's center of rotation on the tunnel centerline. In the next step, slots on the test section floor are sealed and the model is moved to a new center of rotation that is 33 inches below the tunnel centerline. Then, the original angle of attack sweeps are repeated. Afterwards, wall interference corrections are applied to both test data sets and response surface models of the resulting aerodynamic coefficients in interference-free flow are generated. Finally, the response surface models are used to predict the aerodynamic coefficients for a family of angles of attack while keeping dynamic pressure, Mach number, and Reynolds number constant. The validation is considered successful if the corrected aerodynamic coefficients obtained from the related response surface model pair show good agreement. Residual differences between the corrected coefficient sets will be analyzed as well because they are an indicator of the overall accuracy of the facility's wall interference correction process.

wind tunnel wall interferenc↗

Calculations of Supersonic and Hypersonic Flows using Compressible Wall Functions

The present paper presents a numerical procedure to calculate supersonic and hypersonic flows using the compressible law of the wall. The turbulence models under consideration include the Launder-Reece-Rodi-Gibson Reynolds-stress model and the k-epsilon model. The models coupled with the proposed wall function technique have been tested in both separated and unseparated flows. The flows include (1) an insulated flat plate flow over a range of Mach numbers, (2) a Mach 5 flat plate flow with cold wall conditions, (3) a two dimensional supersonic compression corner flow, (4) a hypersonic flow over an axisymmetric flare, and (5) a hypersonic flow over a 2-D compression corner. Results indicate that the wall function technique gives improved predictions of skin friction and heat transfer in separated flows compared with models using wall dampers. Predictions of the extent of separation are not improved over the wall damper models except with the Reynolds-stress model for the supersonic compression corner flow case.

Huang, P. G.↗