Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code comparisons”

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 559 records · Page 31

An improved far field drag calculation method for nonlinear CFD codes

An improved method to calculate drag based on far field conditions is presented and demonstrated. The method is illustrated using two CFD codes, PARC and CFL3D, and examining their ability to preserve lift and drag coefficients along grid line defined integration paths. The flow fields were generated by solving the Euler equations for a NACA 0012 airfoil at a free stream Mach number of 0.8 and angle of attack of 1.25 degrees. In comparison to force coefficients obtained by surface pressure integration, neither code acceptably preserved both force coefficients throughout the near and far fields. An investigation into the relationship between numerical prediction error and calculated force coefficients revealed a direct connection between solution mass conservation error and force coefficients error. A method to correct the predicted force coefficients based on integrated mass conservation error is described and demonstrated. The corrected force coefficients for both codes are shown to be more accurate than the uncorrected values. The correction method is applicable to both two and three dimensions and is independent of the algorithm used to generate the flow field.

Sirbaugh, James R.↗

The effects of finite rate chemical processes on high enthalpy nozzle performance - A comparison between SPARK and SEAGULL

The generalized chemistry version of the computer code SPARK is extended to include two higher-order numerical schemes, yielding fourth-order spatial accuracy for the inviscid terms. The new and old formulations are used to study the influences of finite rate chemical processes on nozzle performance. A determination is made of the computationally optimum reaction scheme for use in high-enthalpy nozzles. Finite rate calculations are compared with the frozen and equilibrium limits to assess the validity of each formulation. In addition, the finite rate SPARK results are compared with the constant ratio of specific heats (gamma) SEAGULL code, to determine its accuracy in variable gamma flow situations. Finally, the higher-order SPARK code is used to calculate nozzle flows having species stratification. Flame quenching occurs at low nozzle pressures, while for high pressures, significant burning continues in the nozzle.

Carpenter, M. H.↗

Comparison of Stopping Power and Range Databases for Radiation Transport Study

The codes used to calculate stopping power and range for the space radiation shielding program at the Langley Research Center are based on the work of Ziegler but with modifications. As more experience is gained from experiments at heavy ion accelerators, prudence dictates a reevaluation of the current databases. Numerical values of stopping power and range calculated from four different codes currently in use are presented for selected ions and materials in the energy domain suitable for space radiation transport. This study of radiation transport has found that for most collision systems and for intermediate particle energies, agreement is less than 1 percent, in general, among all the codes. However, greater discrepancies are seen for heavy systems, especially at low particle energies.

Tai, H.↗

Numerical Calculations of 3-D High-Lift Flows and Comparison with Experiment

Solutions were obtained with the Navier-Stokes CFD code TLNS3D to predict the flow about the NASA Trapezoidal Wing, a high-lift wing composed of three elements: the main-wing element, a deployed leading-edge slat, and a deployed trailing-edge flap. Turbulence was modeled by the Spalart-Allmaras one-equation turbulence model. One case with massive separation was repeated using Menter's two-equation SST (Menter's Shear Stress Transport) k-omega turbulence model in an attempt to improve the agreement with experiment. The investigation was conducted at a free stream Mach number of 0.2, and at angles of attack ranging from 10.004 degrees to 34.858 degrees. The Reynolds number based on the mean aerodynamic chord of the wing was 4.3 x 10 (sup 6). Compared to experiment, the numerical procedure predicted the surface pressures very well at angles of attack in the linear range of the lift. However, computed maximum lift was 5% low. Drag was mainly under predicted. The procedure correctly predicted several well-known trends and features of high-lift flows, such as off-body separation. The two turbulence models yielded significantly different solutions for the repeated case.

Compton, William B, III↗

Analysis of an advanced ducted propeller subsonic inlet

It is shown that a time marching Navier-Stokes code called PARC can be utilized to provide a reasonable prediction of the flow field within an inlet for an advanced ducted propeller. The code validation was implemented for a nonseparated flow condition associated with the inlet functioning at angles-of-attack of zero and 25 deg. Comparison of the computational results with the test data shows that the PARC code with the propeller face fixed flow properties boundary conditions (BC) provided a better prediction of the inlet surface static pressures than the prediction when the mass flow BC was employed.

Iek, Chanthy↗

Summary of the GASP code application and evaluation effort for scramjet combustor flowfields

Numerical simulations of 3D turbulent mixing and reacting flows have been systematically evaluated by comparison with experimental data and other numerical solutions to provide confidence in the General Aerodynamic Simulation Program (GASP) code version 1.3. Specifically, the GASP flow solver has been used to model high speed flow through scram jet combustors. Unit injector problems evaluated include (1) UVA staged sonic normal injection of air into Mach 2 air stream, (2) the VPI Mach 1.7 low angled (15 deg) flush wall helium injection into an unconfined Mach 6 'cold air' stream, and (3) the HYPULSE angled (30 deg) flush wall hydrogen injection into a Mach 6, high enthalpy (flight Mach 17 simulation) nitrogen and air test gases. Results of the numerical simulation for the first two test cases have been compared with detailed in-stream measurements and with SPARK CFD solutions. The results for the last two cases have been compared with measured wall pressure and heat flux data and with SPARK CFD solutions. The present GASP solutions compare favorably with both the experimental data and the SPARK solutions.

Srinivasan, Shivakumar↗

Navier-Stokes predictions of multifunction nozzle flows

A two-dimensional, Navier-Stokes code developed by Imlay based on the implicit, finite-volume method of MacCormack has been applied to the prediction of the flow fields and performance of several nonaxisymmetric, convergent-divergent nozzles with and without thrust vectoring. Comparisons of predictions with experiment show that the Navier-Stokes code can accurately predict both the flow fields and performance for nonaxisymmetric nozzles where the flow is predominantly two-dimensional and at nozzle pressure ratios at or above the design values. Discrepancies between predictions and experiment are noted at lower nozzle pressure ratios where separation typically occurs in portions of the nozzle. The overall trends versus parameters such as nozzle pressure ratio, flap angle, and vector angle were generally predicted correctly.

Wilmoth, Richard G.↗

Code validation for the simulation of supersonic viscous flow about the F-16XL

The viewgraphs and discussion on code validation for the simulation of supersonic viscous flow about the F-16XL are provided. Because of the large potential gains related to laminar flow on the swept wings of supersonic aircraft, interest in the applications of laminar flow control (LFC) techniques in the supersonic regime has increased. A supersonic laminar flow control (SLFC) technology program is currently underway within NASA. The objective of this program is to develop the data base and design methods that are critical to the development of laminar flow control technology for application to supersonic transport aircraft design. Towards this end, the program integrates computational investigations underway at NASA Ames-Moffett and NASA Langley with flight-test investigations being conducted on the F-16XL at the NASA Ames-Dryden Research Facility in cooperation with Rockwell International. The computational goal at NASA Ames-Moffett is to integrate a thin-layer Reynolds averaged Navier-Stokes flow solver with a stability analysis code. The flow solver would provide boundary layer profiles to the stability analysis code which in turn would predict transition on the F-16XL wing. To utilize the stability analysis codes, reliable boundary layer data is necessary at off-design cases. Previously, much of the prediction of boundary layer transition has been accomplished through the coupling of boundary layer codes with stability theory. However, boundary layer codes may have difficulties at high Reynolds numbers, of the order of 100 million, and with the current complex geometry in question. Therefore, a reliable code which solves the thin-layer Reynolds averaged Navier-Stokes equations is needed. Two objectives are discussed, the first in greater depth. The first objective is method verification, via comparisons of computations with experiment, of the reliability and robustness of the code. To successfully implement LFC techniques to the F-16XL wing, the flow about the leading edge must be maintained as laminar flow. Therefore, the second objective is to focus on a series of numerical simulations with different values of angle of attack, alpha, and Reynolds numbers. The purpose of the simulations is to study their effects on the two main factors which precipitate transition to turbulence at leading edges of highly swept wings (e.g., 'spanwise contamination' and 'crossflow instability').

Flores, Jolen↗

Verification and Validation of the k-kL Turbulence Model in FUN3D and CFL3D Codes

The implementation of the k-kL turbulence model using multiple computational uid dy- namics (CFD) codes is reported herein. The k-kL model is a two-equation turbulence model based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's equation is well suited for term-by-term mod- eling and displays useful features compared to other two-equation models. An important di erence is that this formulation leads to the inclusion of higher-order velocity derivatives in the source terms of the scale equations. This can enhance the ability of the Reynolds- averaged Navier-Stokes (RANS) solvers to simulate unsteady ows. The present report documents the formulation of the model as implemented in the CFD codes Fun3D and CFL3D. Methodology, veri cation and validation examples are shown. Attached and sepa- rated ow cases are documented and compared with experimental data. The results show generally very good comparisons with canonical and experimental data, as well as matching results code-to-code. The results from this formulation are similar or better than results using the SST turbulence model.

Abdol-Hamid, Khaled S.↗

Analysis of thermo-chemical nonequilibrium models for carbon dioxide flows

The aerothermodynamics of thermochemical nonequilibrium carbon dioxide flows is studied. The chemical kinetics models of McKenzie and Park are implemented in separate three-dimensional computational fluid dynamics codes. The codes incorporate a five-species gas model characterized by a translational-rotational and a vibrational temperature. Solutions are obtained for flow over finite length elliptical and circular cylinders. The computed flowfields are then employed to calculate Mach-Zehnder interferograms for comparison with experimental data. The accuracy of the chemical kinetics models is determined through this comparison. Also, the methodology of the three-dimensional thermochemical nonequilibrium code is verified by the reproduction of the experiments.

Rock, Stacey G.↗

Grid-Adapted FUN3D Computations for the Second High Lift Prediction Workshop

Contributions of the unstructured Reynolds-averaged Navier-Stokes code FUN3D to the 2nd AIAA CFD High Lift Prediction Workshop are described, and detailed comparisons are made with experimental data. Using workshop-supplied grids, results for the clean wing configuration are compared with results from the structured code CFL3D Using the same turbulence model, both codes compare reasonably well in terms of total forces and moments, and the maximum lift is similarly over-predicted for both codes compared to experiment. By including more representative geometry features such as slat and flap brackets and slat pressure tube bundles, FUN3D captures the general effects of the Reynolds number variation, but under-predicts maximum lift on workshop-supplied grids in comparison with the experimental data, due to excessive separation. However, when output-based, off-body grid adaptation in FUN3D is employed, results improve considerably. In particular, when the geometry includes both brackets and the pressure tube bundles, grid adaptation results in a more accurate prediction of lift near stall in comparison with the wind-tunnel data. Furthermore, a rotation-corrected turbulence model shows improved pressure predictions on the outboard span when using adapted grids.

Lee-Rausch, E. M.↗

A review of current finite difference rotor flow methods

Rotary-wing computational fluid dynamics is reaching a point where many three-dimensional, unsteady, finite-difference codes are becoming available. This paper gives a brief review of five such codes, which treat the small disturbance, conservative and nonconservative full-potential, and Euler flow models. A discussion of the methods of applying these codes to the rotor environment (including wake and trim considerations) is followed by a comparison with various available data. These data include tests of advancing lifting and nonlifting, and hovering model rotors with significant supercritical flow regions. The codes are also compared for computational efficiency.

Caradonna, F. X.↗

Foam on Tile Impact Modeling for the STS-107 Investigation

Following the breakup of the Space Shuttle Columbia during reentry a NASA/Contractor investigation team was formed to examine the probable damage inflicted on Orbiter Thermal Protection System elements by impact of External Tank insulating foam projectiles. The authors formed a working subgroup within the larger team to apply the Smooth Particle Hydrodynamics code SPHC to the damage estimation problem. Numerical models of the Orbiter's tiles and of the Tank's foam were constructed and used as inputs into the code. Material properties needed to properly model the tiles and foam were obtained from other working subgroups who performed tests on these items for this purpose. Two- and three-dimensional models of the tiles were constructed, including the glass outer layer, the main body of LI-900 insulation, the densified lower layer of LI-900, the Nomex felt mounting layer, and the Aluminum 2024 vehicle skin. A model for the BX-250 foam including porous compression, elastic rebound, and surface erosion was developed. Code results for the tile damage and foam behavior were extensively validated through comparison with Southwest Research Institute foam-on-tile impact experiments carried out in 1999. These tests involved small projectiles striking individual tiles and small tile arrays. Following code and model validation we simulated impacts of larger foam projectiles on the examples of tile systems used on the Orbiter. Results for impacts on the main landing gear door are presented in this paper, including effects of impacts at several angles, and of rapidly rotating projectiles. General results suggest that foam impacts on tiles at about 500 mph could cause appreciable damage if the impact angle is greater than about 20 degrees. Some variations of the foam properties, such as increased brittleness or increased density could increase damage in some cases. Rotation up to 17 rps failed to increase the damage for the two cases considered. This does not rule out other cases in which the rotational energy might lead to an increase in tile damage, but suggests that in most cases rotation will not be an important factor.

Stellingwerf, R. F.↗

Euler flutter analysis of airfoils using unstructured dynamic meshes

Modifications to a two-dimensional unsteady Euler code for the aeroelastic analysis of airfoils are described. The modifications involve including the structural equations of motion and their simultaneous time-integration with the governing flow equations. A novel aspect of the capability is that the solutions are obtained using unstructured grids made up of triangles. Comparisons are made with parallel calculations performed using linear theory and a structured grid Euler code to assess the accuracy of the unstructured grid Euler results. Results are presented for a flat plate airfoil and the NACA 0012 airfoil to demonstrate applications of the Euler code for generalized force computations and aeroelastic analysis. In these comparisons, two different finite-volume discretizations of the Euler equations on unstructured meshes were employed. Sensitivity of the Euler results to changes in numerical parameters were also investigated.

Rausch, Russ D.↗

As-Flown Shuttle Micro-Meteoroid Orbital Debris (MMOD) Assessment History

The final flight of the Space Shuttle program was completed on July 21st of 2011. The reusable nature of the Orbiter vehicles has provided NASA with a unique opportunity to inspect and sample spacecraft surfaces exposed to the MMOD environment in low Earth orbit. This paper will provide details of MMOD damages found on Orbiter surfaces after recent flights, as well as putting these damages in context by providing historical comparisons of recent damage to previous flights. The Bumper threat assessment computer code is used by NASA to determine spacecraft MMOD penetration risk. The tool was used before each mission to estimate pre-flight MMOD risk for the Shuttle Program Office using planned vehicles attitudes and exposure times. This paper summarizes the efforts of the authors to assess the expected number MMOD impacts of selected particles sizes of 52 shuttle missions ranging from STS-50 (1992) through STS-133 (2011) using as-flown attitude data. The missions encompass the Shuttle/Mir precursors to the International Space Station (ISS) as well as the ISS assembly flights. The paper includes a comparison of observed MMOD impact damage to the predicted number of impacts from Bumper code assessments. This data can be used to validate engineering models of the orbital debris and micrometeoroid environment, such NASA s ORDEM debris model and MEM meteoroid model.

Hyde, James L.↗

Neutron ages computed from experimental activation data

Computer program reduces time manually required to compute neutron age and to provide definite plan of procedural choices. Program allows convenient comparison of several fitting and error analysis procedures. Program code provides for error estimation of various integration options.

Semler, T. T.↗

Demonstration of the range over which the Langley Research Center digital computer charring ablation program (CHAP) can be used with confidence: Comparisons of CHAP predictions and test data for three ablation materials

Comparisons of ablation calculations with the charring ablation computer code and ablation test data are presented over a wide range of environmental conditions in air for three materials: low-density nylon phenolic, Avcoat 5026-39HC/G, and a filled silicon elastomer. Heat fluxes considered range from over 500 Btu/sq ft-sec to less than 50 Btu/sq ft-sec. Pressures range from 0.5 atm to .004 atm. Enthalpies range from about 2000 Btu/lb to 18000 Btu/lb. Predictions of recession, pyrolysis penetration, and thermocouple responses are considered. Recession predictions for nylon phenolic are good as steady state is approached, but strongly transient cases are underpredicted. Pyrolysis penetrations and thermocouple responses are very well predicted. Recession amounts for Avcoat and silicone elastomer are less well predicted, although high heat flux cases near steady state are fairly satisfactory. Pyrolysis penetrations and thermocouple responses are very well predicted.

Moyer, C. B.↗