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 847 records · Page 47

Effects of rotation on coolant passage heat transfer. Volume 1: Coolant passages with smooth walls

An experimental program was conducted to investigate heat transfer and pressure loss characteristics of rotating multipass passages, for configurations and dimensions typical of modern turbine blades. The immediate objective was the generation of a data base of heat transfer and pressure loss data required to develop heat transfer correlations and to assess computational fluid dynamic techniques for rotating coolant passages. Experiments were conducted in a smooth wall large scale heat transfer model.

Hajek, T. J.↗

Evaluation of a three-dimensional Navier-Stokes solver for transonic flow over a low-aspect-ratio wing mounted in a wind tunnel

A three dimensional Navier-Stokes solver is evaluated for transonic flow over a thin, swept, low-aspect ratio wing. The computational study was undertaken in support of a wind tunnel experimental program. The computational results are compared to experimental surface pressure data obtained in a cryogenic wind tunnel with an adaptive wall test section. The results show favorable agreement over a wide range of conditions, further the numerical results provide additional data of the complex three-dimensional flow field. Differences in the predictions and experiment suggest a need to conduct further experiments to evaluate the adaptive wall testing technique, and to model the tunnel sidewall boundary layer in the computations.

Milholen, William E., II↗

Low Reynolds number multiple-time-scale turbulence model and calculations of steady and pulsating shear layers

A low Reynolds number multiple-time-scale turbulence model (LMS) and its application to fully developed turbulent channel flows and pulsating pipe flows are presented. The LMS can describe the inequilibrium turbulence phenomena down to the viscous sublayer. The calculated fluid flow and turbulence fields for the channel flows are in better agreement with the direct numerical simulation (DNS) results than those obtained using a Reynolds stress turbulence model, and the calculated near-wall dissipation rates are in qualitatively correct agreement with the DNS results. The LMS also successfully predicts the rapidly varying phase-lead of the wall shearing stress that occurs in a narrow range of the dimensionless frequency omega (+) = (omega x nu)/(upsilon x tau(exp2)) for the pulsating pipe flows while various other turbulence models fail to predict this phenomenon, and the LMS yields significantly improved numerical results for a wide range of the dimensionless frequency compared with those obtained using a rapid distortion theory (RDT).

Kim, Sang-Wook↗

Unsteady flows

The development status of tactical missile dynamic behavior prediction is evaluated, with a view to the formulation of guidelines for high angle-of-attack tests where highly nonlinear aerodynamics are encountered and strong cross-coupling arises between lateral and longitudinal degrees of freedom. Attention is given to the Ericsson (1979) inviscid computer code, which is entirely general within the geometric constraints of axisymmetry and low aspect ratio lifting surfaces. A serious problem in the prediction of full-scale high-alpha dynamics is the almost total dependence on subscale experimental data, whose use is severely restricted by Reynolds number scaling as well as by model support and wind tunnel wall interference.

Ericsson, L. E.↗

Solution strategies and heat transfer calculations for three-dimensional configurations at hypersonic speeds

A procedure which reduces the memory requirements for computing the viscous flow over a modified Orbiter geometry at a hypersonic flight condition is presented. The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) code which incorporates a thermochemical nonequilibrium chemistry model, a finite rate catalytic wall boundary condition and wall temperature distribution based on radiation equilibrium is used in this study. In addition, the effect of choice of 'min mod' function, eigenvalue limiter and grid density on surface heating is investigated. The surface heating from a flowfield calculation at Mach number 22, altitude of 230,000 ft and 40 deg angle of attack is compared with flight data from three Orbiter flights.

Weilmuenster, K. J.↗

Dynamics and control of coherent structures in the turbulent wall layer: An overview

We expand the velocity field in the vicinity of the wall in empirical eigenfunctions obtained from experiment. Truncating our system and using Galerkin projection, we obtain a closed set of non-linear ordinary differential equations with ten degrees of freedom. We find a rich dynamical behavior, including in particular a heteroclinic attracting orbit giving rise to intermittency. The intermittent jump from one attracting point to the other resembles in many respects the bursts observed in experiments. Specifically, the time between jumps and the duration of the jumps, is approximately that observed in a burst; the jump begins with the formation of a narrowed and intensified updraft, like the ejection phase of a burst, and is followed by a gentle, diffuse downdraft like the sweep phase of a burst. The magnitude of the Reynolds stress spike produced during a burst is limited by our truncation. The behavior is quite robust, much of it being due to the symmetries present (Aubry's group has examined dimensions up to 128 with persistence of the global behavior). We have examined eigenvalues and coefficients obtained from experiment, and from exact simulation, which differ in magnitude. Similar behavior is obtained in both cases; in the latter case, the heteroclinic orbits connect limit cycles instead of fixed points, corresponding to cross-stream waving of the streamwise rolls. The bifurcation diagram remains structurally similar, but somewhat distorted. The role of the pressure term is made clear - it triggers the intermittent jumps, which otherwise would occur at longer and longer intervals, as the system trajectory is attracted closer and closer to the heteroclinic cycle. The pressure term results in the jumps occurring at essentially random times, and the magnitude of the signal determines the average timing. Stretching of the wall region shows that the model is consistent with observations of polymer drag reduction. Change of the third order coefficients, corresponding to acceleration or deceleration of the mean flow, changes the heteroclinic cycles from attracting to repelling, increasing or decreasing the stability, in agreement with observations. The existence of fixed points is an artifact introduced by the projection; however, a decoupled model still displays the rich dynamics. Numerous assumptions made in Aubry et al. (1988) can now be proved exactly. Feeding back eigenfuncitons with the proper phase can delay the bursting, (the heteroclinic jump to the other fixed point), decreasing the drag. It is also possible to speed up the bursting, increasing mixing to control separation. Our approach is optimal for short time tracking in control.

Berkooz, Gal↗

Turbulence modeling of free shear layers for high-performance aircraft

The High Performance Aircraft (HPA) Grand Challenge of the High Performance Computing and Communications (HPCC) program involves the computation of the flow over a high performance aircraft. A variety of free shear layers, including mixing layers over cavities, impinging jets, blown flaps, and exhaust plumes, may be encountered in such flowfields. Since these free shear layers are usually turbulent, appropriate turbulence models must be utilized in computations in order to accurately simulate these flow features. The HPCC program is relying heavily on parallel computers. A Navier-Stokes solver (POVERFLOW) utilizing the Baldwin-Lomax algebraic turbulence model was developed and tested on a 128-node Intel iPSC/860. Algebraic turbulence models run very fast, and give good results for many flowfields. For complex flowfields such as those mentioned above, however, they are often inadequate. It was therefore deemed that a two-equation turbulence model will be required for the HPA computations. The k-epsilon two-equation turbulence model was implemented on the Intel iPSC/860. Both the Chien low-Reynolds-number model and a generalized wall-function formulation were included.

Sondak, Douglas L.↗

Effects of rotation on coolant passage heat transfer. Volume 2: Coolant passages with trips normal and skewed to the flow

An experimental program was conducted to investigate heat transfer and pressure loss characteristics of rotating multipass passages, for configurations and dimensions typical of modem turbine blades. This experimental program is one part of the NASA Hot Section Technology (HOST) Initiative, which has as its overall objective the development and verification of improved analysis methods that will form the basis for a design system that will produce turbine components with improved durability. The objective of this program was the generation of a data base of heat transfer and pressure loss data required to develop heat transfer correlations and to assess computational fluid dynamic techniques for rotating coolant passages. The experimental work was broken down into two phases. Phase 1 consists of experiments conducted in a smooth wall large scale heat transfer model. A detailed discussion of these results was presented in volume 1 of a NASA Report. In Phase 2 the large scale model was modified to investigate the effects of skewed and normal passage turbulators. The results of Phase 2 along with comparison to Phase 1 is the subject of this Volume 2 NASA Report.

Johnson, B. V.↗

A representation for the turbulent mass flux contribution to Reynolds-stress and two-equation closures for compressible turbulence

The turbulent mass flux, or equivalently the fluctuating Favre velocity mean, appears in the first and second moment equations of compressible kappa-epsilon and Reynolds stress closures. Mathematically it is the difference between the unweighted and density-weighted averages of the velocity field and is therefore a measure of the effects of compressibility through variations in density. It appears to be fundamental to an inhomogeneous compressible turbulence, in which it characterizes the effects of the mean density gradients, in the same way the anisotropy tensor characterizes the effects of the mean velocity gradients. An evolution equation for the turbulent mass flux is derived. A truncation of this equation produces an algebraic expression for the mass flux. The mass flux is found to be proportional to the mean density gradients with a tensor eddy-viscosity that depends on both the mean deformation and the Reynolds stresses. The model is tested in a wall bounded DNS at Mach 4.5 with notable results.

Ristorcelli, J. R.↗

Investigation of Advanced Counterrotation Blade Configuration Concepts for High Speed Turboprop Systems. Task 8: Cooling Flow/heat Transfer Analysis

The focus of this task was to validate the ADPAC code for heat transfer calculations. To accomplish this goal, the ADPAC code was modified to allow for a Cartesian coordinate system capability and to add boundary conditions to handle spanwise periodicity and transpiration boundaries. The primary validation case was the film cooled C3X vane. The cooling hole modeling included both a porous region and grid in each discrete hold. Predictions for these models as well as smooth wall compared well with the experimental data.

Hall, Edward J.↗

Flutter analysis of composite box beams

The dynamic aeroelastic instability of flutter is an important factor in the design of modern high-speed, flexible aircraft. The current trend is toward the creative use of composites to delay flutter. To obtain an optimum design, we need an accurate as well as efficient model. As a first step towards this goal, flutter analysis is carried out for an unswept composite box beam using a linear structural model and Theodorsen's unsteady aerodynamic theory. Structurally, the wing was modeled as a thin-walled box-beam of rectangular cross section. Theodorsen's theory was used to get 2-D unsteady aerodynamic forces, which were integrated over the span. A free-vibration analysis is carried out. These fundamental modes are used to get the flutter solution using the V-g method. Future work is intended to build on this foundation.

Hodges, Dewey H.↗

TranAir: A full-potential, solution-adaptive, rectangular grid code for predicting subsonic, transonic, and supersonic flows about arbitrary configurations. User's manual

The TranAir computer program calculates transonic flow about arbitrary configurations at subsonic, transonic, and supersonic freestream Mach numbers. TranAir solves the nonlinear full potential equations subject to a variety of boundary conditions modeling wakes, inlets, exhausts, porous walls, and impermeable surfaces. Regions with different total temperature and pressure can be represented. The user's manual describes how to run the TranAir program and its graphical support programs.

Johnson, F. T.↗

A Concept for Directly Coupled Pulsed Electromagnetic Acceleration of Plasmas

Plasma jets with high momentum flux density are required for a variety of applications in propulsion research. Methods of producing these plasma jets are being investigated at NASA Marshall Space Flight Center. The experimental goal in the immediate future is to develop plasma accelerators which are capable of producing plasma jets with momentum flux density represented by velocities up to 200 km/s and ion density up to 10(exp 24) per cu m, with sufficient precision and reproducibility in their properties, and with sufficiently high efficiency. The jets must be sufficiently focused to allow them to be transported over several meters. A plasma accelerator concept is presented that might be able to meet these requirements. It is a self-switching, shaped coaxial pulsed plasma thruster, with focusing of the plasma flow by shaping muzzle current distribution as in plasma focus devices, and by mechanical tapering of the gun walls. Some 2-D MHD modeling in support of the conceptual design will be presented.

Thio, Y.C. Francis↗

Splashing Droplets

Current data on droplet breakup is scarce for the sizes and velocities typical of practical applications such as in spray combustion processes and coating processes. While much more representative of practical applications, the small spatial scales and rapid time-scales prevent detailed measurement of the internal fluid dynamics and liquid property gradients produced by impinging upon surfaces. Realized through the extended spatial and temporal scales afforded by a microgravity environment, an improved understanding of drop breakup dynamics is sought to understand and ultimately control the impingement dynamics of droplets upon surfaces in practical situations. The primary objective of this research will be to mark the onset of different 'splashing modes' and to determine their temperature, pressure and angle dependence for impinging droplets representative of practical fluids. In addition, we are modeling the evolution of droplets that do not initially splash but rather undergo a 'fingering' evolution observed on the spreading fluid front and the transformation of these fingers into splashed products. An example of our experimental data is presented below. These images are of Isopar V impacting a mirror-polished surface. They were acquired using a high-speed camera at 1000 frames per second. They show the spreading of a single droplet after impact and ensuing finger instabilities. Normal gravity experimental data such as this will guide low gravity measurements in the 2.2 second drop tower and KC-135 aircraft as available. Presently we are in the process of comparing the experimental data of droplet shape evolution to numerical models, which can also capture the internal fluid dynamics and liquid property gradients such as produced by impingement upon a heated surface. To-date isothermal numerical data has been modeled using direct numerical simulations of representative splashing droplets. The data obtained so far indicates that the present model describes well the droplet wall interactions to a point in time just before splash. Additional information is included in the original extended abstract.

VanderWal, Randall L.↗

Volatile-rich Crater Interior Deposits in the Polar Regions of Mars: Evidence for Ice Cap Advance and Retreat

Many craters on Mars are partially filled by distinctive material emplaced by post-impact processes. This crater fill material is an interior mound which is generally separated from the walls of the crater by a trough that may be continuous along the crater circumference (i.e. a ring-shaped trough), or which may only partially contact the crater walls (i.e. a crescent-shaped trough). The fill deposit is frequently offset from the crater center and may be asymmetric in plan view. Populations of such craters include those in the circum-south polar cap region, in Arabia Terra, associated with the Medusae Fossae Formation, and in the northern lowlands proximal to the north polar cap. We focus on those craters in circumpolar regions and assess their relationship to polar cap advance and retreat, especially the possibility that fill material represents remnants of a formerly larger contiguous cap. Volatile-rich deposits have the property of being modifiable by the local stability of the solid volatile, which is governed by local energy balance. Here we test the hypothesis that asymmetries in volatile fill shape, profile, and center-location within a crater result from asymmetries in local energy balance within the crater, due mainly to variation of solar insolation and radiative effects of the crater walls over the crater interior. Model profiles of crater fill are compared with MOLA topographic profiles to assess this hypothesis. If asymmetry in morphology and location of crater fill are consistent with radiative-dominated asymmetries in energy budget within the crater, then 1) the volatile-rich composition of the fill is supported (this process should not be effective at shaping volcanic or sedimentary deposits), and 2) the dominant factor determining the observed shape of volatile-rich crater fill is the local radiative energy budget (and erosive processes such as eolian deflation are secondary or unnecessary). We also use a geographic and energy model approach to specifically test the idea that material in partially filled craters around the south pole may once have been contiguous to the cap and may have been sustained and modified by radiative processes specific to the crater environment (as opposed to the surrounding plains) as the cap retreated.

Russell, Patrick S.↗

Detached-Eddy Simulations of Separated Flow Around Wings With Ice Accretions: Year One Report

A computational investigation was performed to assess the effectiveness of Detached-Eddy Simulation (DES) as a tool for predicting icing effects. The AVUS code, a public domain flow solver, was employed to compute solutions for an iced wing configuration using DES and steady Reynolds Averaged Navier-Stokes (RANS) equation methodologies. The configuration was an extruded GLC305/944-ice shape section with a rectangular planform. The model was mounted between two walls so no tip effects were considered. The numerical results were validated by comparison with experimental data for the same configuration. The time-averaged DES computations showed some improvement in lift and drag results near stall when compared to steady RANS results. However, comparisons of the flow field details did not show the level of agreement suggested by the integrated quantities. Based on our results, we believe that DES may prove useful in a limited sense to provide analysis of iced wing configurations when there is significant flow separation, e.g., near stall, where steady RANS computations are demonstrably ineffective. However, more validation is needed to determine what role DES can play as part of an overall icing effects prediction strategy. We conclude the report with an assessment of existing computational tools for application to the iced wing problem and a discussion of issues that merit further study.

Choo, Yung K.↗

Efficient Fluid Dynamic Design Optimization Using Cartesian Grids

This report is subdivided in three parts. The first one reviews a new approach to the computation of inviscid flows using Cartesian grid methods. The crux of the method is the curvature-corrected symmetry technique (CCST) developed by the present authors for body-fitted grids. The method introduces ghost cells near the boundaries whose values are developed from an assumed flow-field model in vicinity of the wall consisting of a vortex flow, which satisfies the normal momentum equation and the non-penetration condition. The CCST boundary condition was shown to be substantially more accurate than traditional boundary condition approaches. This improved boundary condition is adapted to a Cartesian mesh formulation, which we call the Ghost Body-Cell Method (GBCM). In this approach, all cell centers exterior to the body are computed with fluxes at the four surrounding cell edges. There is no need for special treatment corresponding to cut cells which complicate other Cartesian mesh methods.

Dadone, A.↗

RANS and DES Computations for a Wing with Ice Accretion

A computational investigation was performed to assess the effectiveness of Detached Eddy Simulation (DES) as a tool for predicting icing effects. The AVUS code, a public domain flow solver, was employed to compute solutions for an iced wing configuration using DES and steady Reynolds Averaged Navier-Stokes (RANS) equation methodologies. The wing section considered here was a business jet airfoil (GLC305) with a 22.5-minute glaze ice accretion (944-ice shape). The section was extruded to form a rectangular planform. The model was mounted between two walls so no tip effects were considered. The numerical results were validated by comparison with experimental data for the same configuration. The time-averaged DES computations showed some improvement in lift and drag results near stall when compared to steady RANS results. However, comparisons of the flow field details did not show the level of agreement suggested by the integrated quantities. More validation is needed to determine what role DES can play as part of an overall icing effects prediction strategy.

Mogili, P.↗