Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Agreements”

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 109 records · Page 6

A comparison between implicit and hybrid methods for the calculation of steady and unsteady inlet flows

Numerical simulation of steady and unsteady transonic diffuser flows using two different computer codes are discussed and compared with experimental data. The codes solve the Reynolds-averaged, compressible, Navier-Stokes equations using various turbulence models. One of the codes has been applied extensively to diffuser flows and uses the hybrid method of MacCormack. This code is relatively inefficient numerically. The second code, which was developed more recently, is fully implicit and is relatively efficient numerically. Simulations of steady flows using the implicit code are shown to be in good agreement with simulations using the hybrid code. Both simulations are in good agreement with experimental results. Simulations of unsteady flows using the two codes are in good qualitative agreement with each other, although the quantitative agreement is not as good as in the steady flow cases. The implicit code is shown to be eight times faster than the hybrid code for unsteady flow calculations and up to 32 times faster for steady flow calculations. Results of calculations using alternative turbulence models are also discussed.

Coakley, T. J.↗

Evaluation of the ERBE scene identification algo-rithm

The sensitivity of the radiation budget parameters and the scene selection process to different factors is evaluated. The use of ERB-7 CLE models provides instantaneous albedo and longwave flux values which are in essential agreement in the 70 deg cutoff case with those from the SAB method. The increase in albedo with increasing satellite zenith angle cutoff is not apparent at the target area level for different surface types. GOES models seem to show an increase in instantaneous albedo at the global level with satellite zenith angle of the same nature as the ERB-7 models investigated in a previous report and are, therefore, probably not an improvement. The use of recently derived NCLE models did not make a noticeable change in the budget parameters; but the cloud classification did show hemispherical differences and caused a day-night readjustment of cloudiness amounts. Preliminary results for an additional month (December 1979) indicate good agreement between the SAB and MLE methods. Additional work is required to establish the agreement for all seasons. Comparison of derived cloud amounts with other data sets such as THIR/TOMS indicates good zonal agreement. Sampling adequacy investigations at different temporal averaging intervals with the SAB method, indicate large uncertainties in small time average cases. For a 6-day averaging period, sampling is very poor with only 23% of the globe contributing toward the global mean albedo. Effect of modifications to the MLE procedure seem to have little effect on the derived budget quantities in an averaged sense. Significant differences in the flux values due to differences in scene selection are apparent in individual target areas studies. A modification of the MLE procedure to mimic the perpendicular bisector algorithm indicated no effect on the gross radiation budget quantities.

Vemury, S. K.↗

Comparison of nonequilibrium viscous-shock-layer solutions with windward surface Shuttle heating data

A three-dimensional viscous-shock-layer code has been modified and used to compute the nonequilibrium flowfield over the windward surface of the Space Shuttle Orbiter for reentry conditions between 75 and 60 km. Effects of the modifications are demonstrated by comparison of present heat-transfer predictions with earlier results. The present predictions show a roughly 30 percent decrease in heat transfer along the windward centerline in comparison with previous results and are in agreement with two-dimensional viscous-shock-layer results. The latter agreement indicates that three-dimensional effects are not as significant as previously reported. Windward symmetry plane and off-centerline heating predictions with the modified code are compared with flight data from STS-2 and STS-3. Windward-centerline heating predictions obtained with a recent expression for oxygen recombination at the Shuttle surface were found to be in generally good agreement with the flight data. Comparisons of heating predictions on the off-centerline windward surface were also in good agreement with the data, and the calculated results followed trends in the crossflow heating distributions.

Thompson, R. A.↗

Euler equation computations for the flow over a hovering helicopter rotor

A numerical solution technique is developed for computing the flow field around an isolated helicopter rotor in hover. The flow is governed by the compressible Euler equations which are integrated using a finite volume approach. The Euler equations are coupled to a free wake model of the rotary wing vortical wake. This wake model is incorporated into the finite volume solver using a prescribed flow, or perturbation, technique which eliminates the numerical diffusion of vorticity due to the artificial viscosity of the scheme. The work is divided into three major parts: (1) comparisons of Euler solutions to experimental data for the flow around isolated wings show good agreement with the surface pressures, but poor agreement with the vortical wake structure; (2) the perturbation method is developed and used to compute the interaction of a streamwise vortex with a semispan wing. The rapid diffusion of the vortex when only the basic Euler solver is used is illustrated, and excellent agreement with experimental section lift coefficients is demonstrated when using the perturbation approach; and (3) the free wake solution technique is described and the coupling of the wake to the Euler solver for an isolated rotor is presented. Comparisons with experimental blade load data for several cases show good agreement, with discrepancies largely attributable to the neglect of viscous effects. The computed wake geometries agree less well with experiment, the primary difference being that too rapid a wake contraction is predicted for all the cases.

Roberts, Thomas Wesley↗

Comparison of experimental and calculated attachment rate constants for CFCl3 and CCl4 in the temperature range 294-500 K

Electron-attachment cross sections and rate constants have been measured and calculated for the dissociative attachment processes e + CFCl3 - Cl(-) + CFCl2 and e + CCl4 - Cl(-) + CCl3. Good agreement over the electron-energy range 1-200 meV is found in energy dependence between present calculated cross sections and experimental cross sections at 300 K. The same calculation, with suitable adjustment of thermal populations, was used to calculate electron-attachment rate constants in the range 50-600 K. Experimental rate constants for CFCl3 and CCl4 were measured at temperatures of 294, 404, and 496 K (CFCl3) and 294, 400, and 500 K (CCl4) using the Cavalleri electron-density sampling method. Good agreement is found between present measurements and calculations, poor agreement with flowing-afterglow Langmuir-probe (FALP) data in CFCl3 at the higher temperatures, and reasonable agreement with FALP data for CCl4.

Orient, O. J.↗

A Detailed Investigation of Staged Normal Injection into a Mach 2 Flow

A study of the staged injection of two jets of air behind a rearward facing step into a Mach 2 flow was performed using the SPARK 3-D Navier-Stokes code. Calculated mole fraction distributions were compared with an extensive set of planar mole fraction measurements made with a laser induced iodine fluorescence technique. A statistical measure, the standard deviation, was used to help assess agreement between calculation and experiment. Overall, good agreement was found between calculated and measured values. Generally, agreement was better in the far field of the injectors. The effect of grid resolution was investigated by calculating solutions on grids of 60,000, 200,000, and 450,000 points. Differences in the solutions on the two finer grids were small. However, the mole fraction distributions were distinguishable. The effect of turbulence modeling was investigated by employing three different algebraic models for the jet turbulence: the Baldwin-Lomax model, the Prandtl mixing length model, and the Eggers mixing length model. Overall, the Eggers mixing length model was found to be superior for this case. Finally, the effect of the jet exit conditions was examined. A recently proposed Mach number distribution at the jet exit was found to slightly improve agreement between measurement and calculation.

Eklund, Dean R.↗

Assessment of a prototype sunphotometer for network applications, phase 1

Data taken in Tucson, AZ, for comparing the CE sunphotometer's results against the Univ. of Arizona's no. 3 solar radiometer are in very good agreement. The spectral voltage intercepts obtained on 15 different days taken with the CE instrument during a seven week period are stable and follow the same daily trends as the voltage intercepts for similar wavelengths contained in the Univ. of Arizona's no. 3 radiometer. The spectral optical depths obtained from both instruments are also in very good agreement and follow the same trends. The only problem that exists in the CE instrument does not have a heater to stabilize the detector temperature which is important around 1000 nm. Zenith sky radiance scans are in good agreement, but both instruments are biased by direct solar rays entering the instrument's respective entrance apertures. All diffuse light data shown were, therefore, taken with the instruments shielded from the direct sun. The CE sunphotometer results pertaining to the water vapor wavelength's are also in good agreement with the Univ. of Arizona's no. 3 and should prove useful for water vapor retrieval.

Reagan, J. A.↗

Mergers, cooling flows, and evaporation

Mergers (the capture of cold gas, especially) can have a profound influence on the hot coronal gas of early-type galaxies and clusters, potentially inducing symptoms hitherto attributed to a cooling flow, if thermal conduction is operative in the coronal plasma. Heat can be conducted from the hot phase into the cold phase, simultaneously ionizing the cold gas to make optical filaments, while locally cooling the coronal gas to mimic a cooling-flow. If there is heat conduction, though, there is no standard cooling-flow since radiative losses are balanced by conduction and not mass deposition. Amongst the strongest observational support for the existence of cooling-flows is the presence of intermediate temperature gas with x-ray emission-line strengths in agreement with cooling-flow models. Here, x-ray line strengths are calculated for this alternative model, in which mergers are responsible for the observed optical and x-ray properties. Since gas around 10(exp 4) K is thermally stable, the cold cloud need not necessarily evaporate and hydrostatic solutions exist. Good agreement with the x-ray data is obtained. The relative strengths of intermediate temperature x-ray emission lines are in significantly better agreement with a simple conduction model than with published cooling-flow models. The good agreement of the conduction model with optical, infrared and x-ray data indicates that significantly more theoretical effort into this type of solution would be profitable.

Sparks, W. B.↗

Dynamic stall of an oscillating wing. Part 1: Evaluation of turbulence models

Unsteady flowfields of a two-dimensional oscillating wing are calculated using an implicit, finite-difference, Navier-Stokes numerical scheme using five widely used turbulence models. The objective of this study is to identify an appropriate turbulence model for accurate simulation of three-dimensional dynamic stall. Three unsteady flow conditions corresponding to attached flow, light-stall, and deep-stall of an oscillating wing experiment were chosen as test cases for computations. Results of unsteady airload hysteresis curves, harmonics of unsteady pressures, and instantaneous flow pictures are presented. Comparison of unsteady airloads with experiment show that all models are deficient in some sense and not a single model predicts all airloads consistently and in agreement with experiment for all flow conditions. For the attached flow condition, the Renormalization Group Theory (RNG), the Johnoson-King (J-K), and the Spalart-Allmaras (S-A) models have better performance. The Baldwin-Lomax (B-L) and the Baldwin-Barth (B-B) models fair poorly. At the light-stall condition, the results for the RNG, the J-K, and S-A models are in agreement with experiment for the upstroke but they all over predict the separation shown by the experiment and therefore have bigger hysteresis loops than experimental results. The B-B model results are also in good agreement for upstroke but have poor lift hysteresis for downstroke. It has superior drag and pitching-moment predictions. For deep-stall conditions, the airloads for the RNG, the B -B, and the S-A models have fair agreement with experiment, but the B-B model performed better at the extreme deep-stall condition. Overall, the RNG model provides significant improvement over the B-L model in all flow regimes with no additional computational cost. The Baldwin-Barth model is the most expensive of the models considered here, costing about 2.5 times that of the Baldwin-Lomax model. Finally, a brief discussion of the effects of grid density, time-step size, and numerical dissipation on the unsteady solutions are also presented.

Srinivasan, G. R.↗

Planetary Boundary Layer Simulation Using TASS

Boundary conditions to an existing large-eddy simulation model have been changed in order to simulate turbulence in the atmospheric boundary layer. Several options are now available, including the use of a surface energy balance. In addition, we compare convective boundary layer simulations with the Wangara and Minnesota field experiments as well as with other model results. We find excellent agreement of modelled mean profiles of wind and temperature with observations and good agreement for velocity variances. Neutral boundary simulation results are compared with theory and with previously used models. Agreement with theory is reasonable, while agreement with previous models is excellent.

Schowalter, David G.↗

Molecular Processes in Comets

The research focused on molecular hydrogen and its response to ultraviolet radiation, photoelectron impact excitation and X-ray radiation and on the interpretation of the ultraviolet spectra of the Jupiter dayglow and auroras. A systematic effort was made to obtain reliable rate coefficients for rovibrational energy transfer of H2, particularly in collisions with hydrogen atoms. We carried out elaborate quantum-mechanical calculations of the scattering of H and H2 on what had been characterized as a reliable potential energy surface of the H3 molecular system. These calculations took into account reactive channels and rate coefficients for ortho-para transitions were obtained. Quantal calculations are too complex to be applied to all the possible rovibrational transitions and we turned to a semi-classical method. With it we calculated rate coefficients for transitions, reactive and non-reactive, for all the rovibrational levels. We carried out the calculations for three of the available H3 potential energy surfaces. We discovered an unexpected sensitivity of the rate coefficients for the non-reactive channels to the potential energy surface. This discovery stimulated more thorough investigations of the potential energy surface elsewhere and to the construction of a new surface. We have used it in further semi-classical calculations in work that will shortly be 2 completed and which, together with new quantum-mechanical calculations should comprise a set of reliable rate coefficients that can be used in discussions of H2 on the Jovian planets. We carried out a detailed study of the Jovian ultraviolet dayglow. There has been a long-running argument about the dayglow on Jupiter. There are two sources of excitation: fluorescence and photoelectron impact excitation. It had been argued that a third source "the electron glow" was needed to bring theory and observation into agreement. We believe we have shown conclusively that the third source is unnecessary. We have achieved a close quantitative agreement between the predicted spectrum arising from fluorescence and photoelectrons and the measured spectrum. We also demonstrated a method by which the presence of HD could be established observationary. Similar calculations were carried out of the ultraviolet spectra of the Jovian auroras. We again found close agreement. Indeed, the agreement was so detailed we were able to derive the temperature of the atmosphere. We found, contrary to the standard model of the time, a high temperature between 400 K and 600 K and established the presence of a significant temperature gradient. A large temperature gradient was indeed found by direct sampling with the Jupiter probe.

Dalgarno, A.↗

Simulation of Long Lived Tracers Using an Improved Empirically Based Two-Dimensional Model Transport Algorithm

We have developed a new empirically-based transport algorithm for use in our GSFC two-dimensional transport and chemistry model. The new algorithm contains planetary wave statistics, and parameterizations to account for the effects due to gravity waves and equatorial Kelvin waves. As such, this scheme utilizes significantly more information compared to our previous algorithm which was based only on zonal mean temperatures and heating rates. The new model transport captures much of the qualitative structure and seasonal variability observed in long lived tracers, such as: isolation of the tropics and the southern hemisphere winter polar vortex; the well mixed surf-zone region of the winter sub-tropics and mid-latitudes; the latitudinal and seasonal variations of total ozone; and the seasonal variations of mesospheric H2O. The model also indicates a double peaked structure in methane associated with the semiannual oscillation in the tropical upper stratosphere. This feature is similar in phase but is significantly weaker in amplitude compared to the observations. The model simulations of carbon-14 and strontium-90 are in good agreement with observations, both in simulating the peak in mixing ratio at 20-25 km, and the decrease with altitude in mixing ratio above 25 km. We also find mostly good agreement between modeled and observed age of air determined from SF6 outside of the northern hemisphere polar vortex. However, observations inside the vortex reveal significantly older air compared to the model. This is consistent with the model deficiencies in simulating CH4 in the northern hemisphere winter high latitudes and illustrates the limitations of the current climatological zonal mean model formulation. The propagation of seasonal signals in water vapor and CO2 in the lower stratosphere showed general agreement in phase, and the model qualitatively captured the observed amplitude decrease in CO2 from the tropics to midlatitudes. However, the simulated seasonal amplitudes were attenuated too rapidly with altitude in the tropics. Overall, the simulations with the new transport formulation are in substantially better agreement with observations compared with our previous model transport.

Fleming, E. L.↗

Application of an Unstructured Grid Navier-Stokes Solver to a Generic Helicopter Boby: Comparison of Unstructured Grid Results with Structured Grid Results and Experimental Results

An unstructured-grid Navier-Stokes solver was used to predict the surface pressure distribution, the off-body flow field, the surface flow pattern, and integrated lift and drag coefficients on the ROBIN configuration (a generic helicopter) without a rotor at four angles of attack. The results are compared to those predicted by two structured- grid Navier-Stokes solvers and to experimental surface pressure distributions. The surface pressure distributions from the unstructured-grid Navier-Stokes solver are in good agreement with the results from the structured-grid Navier-Stokes solvers. Agreement with the experimental pressure coefficients is good over the forward portion of the body. However, agreement is poor on the lower portion of the mid-section of the body. Comparison of the predicted surface flow patterns showed similar regions of separated flow. Predicted lift and drag coefficients were in fair agreement with each other.

Mineck, Raymond E.↗

Vertical Profiles of Latent Heat Release and Their Retrieval for TOGA COARE Convective Systems Using a Cloud Resolving Model, SSM/I, and Ship-borne Radar Data

Latent heating profiles associated with three (TOGA COARE) Tropical Ocean and Global Atmosphere Coupled Ocean Atmosphere Response Experiment active convective episodes (December 10-17 1992; December 19-27 1992; and February 9-13 1993) are examined using the Goddard Cumulus Ensemble (GCE) Model and retrieved by using the Goddard Convective and Stratiform Heating (CSH) algorithm . The following sources of rainfall information are input into the CSH algorithm: Special Sensor Microwave Imager (SSM/1), Radar and the GCE model. Diagnostically determined latent heating profiles calculated using 6 hourly soundings are used for validation. The GCE model simulated rainfall and latent heating profiles are in excellent agreement with those estimated by soundings. In addition, the typical convective and stratiform heating structures (or shapes) are well captured by the GCE model. Radar measured rainfall is smaller than that both estimated by the GCE model and SSM/I in all three different COARE IFA periods. SSM/I derived rainfall is more than the GCE model simulated for the December 19-27 and February 9-13 periods, but is in excellent agreement with the GCE model for the December 10-17 period. The GCE model estimated stratiform amount is about 50% for December 19-27, 42% for December 11-17 and 56% for the February 9-13 case. These results are consistent with large-scale analyses. The accurate estimates of stratiform amount is needed for good latent heating retrieval. A higher (lower) percentage of stratiform rain can imply a maximum heating rate at a higher (lower) altitude. The GCE model always simulates more stratiform rain (10 to 20%) than the radar for all three convective episodes. SSM/I derived stratiform amount is about 37% for December 19-27, 48% for December 11-17 and 41% for the February 9-13 case. Temporal variability of CSH algorithm retrieved latent heating profiles using either GCE model simulated or radar estimated rainfall and stratiform amount is in good agreement with that diagnostically determined for all three periods. However, less rainfall and a smaller stratiform percentage estimated by radar resulted in a weaker (underestimated) latent heating profile and a lower maximum latent heating level compared to those determined diagnostically. Rainfall information from SSM/I can not retrieve individual convective events due to poor temporal sampling. Nevertheless, this study suggests that a good 4r, rainfall retrieval from SSM/I for a convective event always leads to a good latent heating retrieval. Sensitivity testing has been performed and the results indicate that the SSM/I derived time averaged stratiform amount may be underestimated for December 19-27. Time averaged heating profiles derived from SSM/I, however, are not in bad agreement with those derived by soundings for the December 10-17 convective period. The heating retrievals may be more accurate for longer time scales provided there is no bias in the sampling.

Tao, Wei-Kuo↗

A Numerical Evaluation of Icing Effects on a Natural Laminar Flow Airfoil

As a part of CFD code validation efforts within the Icing Branch of NASA Glenn Research Center, computations were performed for natural laminar flow (NLF) airfoil, NLF-0414. with 6 and 22.5 minute ice accretions. Both 3-D ice castings and 2-D machine-generated ice shapes were used in wind tunnel tests to study the effects of natural ice is well as simulated ice. They were mounted in the test section of the Low Turbulence Pressure Tunnel (LTPT) at NASA Langley that the 2-dimensionality of the flow can be maintained. Aerodynamic properties predicted by computations were compared to data obtained through the experiment by the authors at the LTPT. Computations were performed only in 2-D and in the case of 3-D ice, the digitized ice shape obtained at one spanwise location was used. The comparisons were mainly concentrated on the lift characteristics over Reynolds numbers ranging from 3 to 10 million and Mach numbers ranging from 0.12 to 0.29. WIND code computations indicated that the predicted stall angles were in agreement with experiment within one or two degrees. The maximum lift values obtained by computations were in good agreement with those of the experiment for the 6 minute ice shapes and the minute 3-D ice, but were somewhat lower in the case of the 22.5 minute 2-D ice. In general, the Reynolds number variation did not cause much change in the lift values while the variation of Mach number showed more change in the lift. The Spalart-Allmaras (S-A) turbulence model was the best performing model for the airfoil with the 22.5 minute ice and the Shear Stress Turbulence (SST) turbulence model was the best for the airfoil with the 6 minute ice and also for the clean airfoil. The pressure distribution on the surface of the iced airfoil showed good agreement for the 6 minute ice. However, relatively poor agreement of the pressure distribution on the upper surface aft of the leading edge horn for the 22.5 minute ice suggests that improvements are needed in the grid or turbulence models.

Chung, James J.↗

Airborne Sunphotometer Measurements of Aerosol Optical Depth and Water Vapor in ACE-Asia and Their Comparisons to Correlative Measurements

In the Spring 2001 phase of the Asian Pacific Regional Aerosol Characterization Experiment (ACE-Asia), the 6-channel NASA Ames Airborne Tracking Sunphotometer (AATS-6) operated on 15 of the 19 research flights of the NCAR C-130, while its 14-channel counterpart (AATS-14) flew successfully on all 19 research flights of the CIRPAS Twin Otter. ACE-Asia studied aerosol outflow from the Asian continent to the Pacific basin. It was designed to integrate suborbital and satellite measurements and models to reduce the uncertainty in calculations of the climate forcing due to aerosols. AATS-6 and AATS-14 measured solar beam transmission at six and 14 wavelengths (380-1021 and 354-1558 nm, respectively), yielding aerosol optical depth (AOD) spectra and columnar water vapor (CWV). Vertical differentiation in profiles yielded aerosol extinction spectra and water vapor concentration. In this paper, we plan to present examples of the following, preliminary findings that are based in part on our airborne sunphotometer measurements: (1) The wavelength dependence of sunphotometer-derived AOD and extinction indicates that supermicron dust was often a major component of the aerosol, frequently extending to high altitudes. The percentage of full-column AOD (525 nm) that Jay above 3 km was typically 34+/-13%. In contrast, the analogous percentage of columnar water vapor was only 10+/-4%; (2) Initial comparison studies between AOD data obtained by AATS-6 and AATS-14 during coordinated low-level flight legs show agreement well within the instruments' error bars; (3) Aerosol extinction has been derived from airborne in situ measurements of scattering (nephelometers) and absorption (particle soot/ absorption photometer, PSAP) or calculated from particle size distribution measurements (mobility analyzers and aerodynamic particle sizers). Comparison with corresponding extinction values derived from the Ames airborne sunphotometer measurements shows good agreement for the vertical distribution of aerosol layers. However, the level of agreement in absolute magnitude of the derived aerosol extinction/optical depth varied among the aerosol layers sampled; (4) Initial comparisons of sunphotometer and satellite-derived AOD using SeaWiFS, MISR and AVHRR show promising results. We also plan to include comparisons with MODIS and TOMS: (5) Initial comparisons of sunphotometer-derived AOD and aerosol extinction profiles with lidars in Tokyo and on a ship show reasonable agreement.

Schmid, B.↗

Experimental Study of the Richtmyer-Meshkov Instability of Incompressible Fluids

The Richtmyer-Meshkov instability of a low Atwood number, miscible, two-liquid system is investigated experimentally. The initially stratified fluids are contained within a rectangular tank mounted to a sled that rides on a vertical set of rails. The instability is generated by dropping the sled onto a coil spring, producing a nearly impulsive upward acceleration. The subsequent freefall that occurs as the container travels upward and then downward on the rails allows the instability to evolve in the absence of gravity. The interface separating the two liquids initially has a well-defined, sinusoidal perturbation that quickly inverts and then grows in amplitude after undergoing the impulsive acceleration. Disturbance amplitudes are measured and compared to theoretical predictions. Linear stability theory gives excellent agreement with the measured initial growth rate, a(sub 0), for single-mode perturbations with the predicted amplitudes differing by less than 10% from experimental measurements up to a nondimensional time ka(sub 0)t = 0.7, where k is the wavenumber. Linear stability theory also provides excellent agreement for the individual mode amplitudes of multi-mode initial perturbations up until the interface becomes multi-valued. Comparison with previously published weakly nonlinear single-mode models shows good agreement up to ka(sub 0)t = 3, while published nonlinear single-mode models provide good agreement up to ka(sub 0)t = 30. The effects of Reynolds number on the vortex core evolution and overall growth rate of the interface are also investigated. Measurements of the overall amplitude are found to be unaffected by the Reynolds number for the range of values studied here. However, experiments carried out at lower values of Reynolds numbers were found to have decreased vortex core rotation rates. In addition, an instability in the vortex cores is observed.

Niederhaus, Charles↗

Computations of Internal and External Axisymmetric Nozzle Aerodynamics at Transonic Speeds

Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles have been completed in support of NASA's Next Generation Launch Technology Program to investigate the effects of high-speed nozzle geometries on the nozzle internal flow and the surrounding boattail regions. These computations span the very difficult transonic flight regime, with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined in Wind, including an Explicit Algebraic Stress model (EASM). Computations on two nozzle geometries have been completed at freestream Mach numbers ranging from 0.6 to 0.9, driven by nozzle pressure ratios (NPR) ranging from 2.9 to 5. Results obtained on converging-only geometry indicate reasonable agreement to experimental data, with the EASM and Shear Stress Transport (SST) turbulence models providing the best agreement. Calculations completed on a converging-diverging geometry involving large-scale internal flow separation did not converge to a true steady-state solution when run with variable timestepping (steady-state). Calculations obtained using constant timestepping (time-accurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was found to be the result of difficulties in using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show fairly good agreement with experimental data. The SST turbulence model demonstrates the best over-all agreement with experimental data.

Dalbello, Teryn↗