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 1,729 records · Page 96

A Revised Validation Process for Ice Accretion Codes

A research project is underway at NASA Glenn to produce computer software that can accurately predict ice growth under any meteorological conditions for any aircraft surface. This report will present results from the latest LEWICE release, version 3.5. This program differs from previous releases in its ability to model mixed phase and ice crystal conditions such as those encountered inside an engine. It also has expanded capability to use structured grids and a new capability to use results from unstructured grid flow solvers. A quantitative comparison of the results against a database of ice shapes that have been generated in the NASA Glenn Icing Research Tunnel (IRT) has also been performed. This paper will extend the comparison of ice shapes between LEWICE 3.5 and experimental data from a previous paper. Comparisons of lift and drag are made between experimentally collected data from experimentally obtained ice shapes and simulated (CFD) data on simulated (LEWICE) ice shapes. Comparisons are also made between experimentally collected and simulated performance data on select experimental ice shapes to ensure the CFD solver, FUN3D, is valid within the flight regime. The results show that the predicted results are within the accuracy limits of the experimental data for the majority of cases.

turbulent heat flux↗

Minimum-mass designs of stiffened graphite/polyimide compression panels

This study presents results from an analytic investigation to determine minimum-mass designs of stiffened Gr/Pi compression panels for a wide range of uniaxial compression loads. Four panel configurations are considered: (1) hat stiffened-laminated skin, (2) hat stiffened-honeycomb core sandwich skin, (3) blade stiffened-laminated skin and (4) blade stiffened-honeycomb core sandwich skin panels. Designs are generated by an automatic optimization computer code entitled PASCO (Panel Analysis and Sizing COde). The analysis portion of PASCO contains a generalized buckling analysis which calculates all critical modes without assuming any of them a priori. Mass-strength curves are presented for the four configurations along with specific designs generated for low, moderate and high values of compression load. A mass-strength curve for a titanium hat stiffened panel is provided for comparison. A 30% to 50% mass savings over titanium panels is shown for moderately and heavily loaded Gr/Pi panels.

Weaver, G. G., II↗

Conservative full-potential calculations for axisymmetric, transonic flow

A conservative, finite-difference, full-potential relaxation code has been developed to solve transonic flow around axisymmetric or nonlifting, planar, two-dimensional bodies. The program utilizes the artificial compressibility method to provide an upwind bias in supersonic regions. Calculated examples include a wide variety of axisymmetric and planar two-dimensional shapes with various blunt, pointed and open ends in subsonic, transonic, and low supersonic free streams. Comparisons between conservative and nonconservative full-potential calculations show perfect agreement at convergence when the flow is entirely subsonic. Noticeable differences exist, however, between the conservative and nonconservative solutions for transonic flows with shocks.

Green, L. L.↗

Aerodynamic design for improved manueverability by use of three-dimensional transonic theory

Improvements in transonic maneuver performance by the use of three-dimensional transonic theory and a transonic design procedure were examined. The FLO-27 code of Jameson and Caughey was used to design a new wing for a fighter configuration with lower drag at transonic maneuver conditions. The wing airfoil sections were altered to reduce the upper-surface shock strength by means of a design procedure which is based on the iterative application of the FLO-27 code. The plan form of the fighter configuration was fixed and had a leading edge sweep of 45 deg and an aspect ratio of 3.28. Wind-tunnel tests were conducted on this configuration at Mach numbers from 0.60 to 0.95 and angles of attack from -2 deg to 17 deg. The transonic maneuver performance of this configuration was evaluated by comparison with a wing designed by empirical methods and a wing designed primarily by two-dimensional transonic theory. The configuration designed by the use of FLO-27 had the same or lower drag than the empirical wing and, for some conditions, lower drag than the two-dimensional design. From some maneuver conditions, the drag of the two-dimensional design was somewhat lower.

Mann, M. J.↗

Flight vehicle propellant surface motion by numerical method

The marker and cell technique is applied to three idealized hydrodynamic situations of interest to the Space Shuttle External Tank (ET). These applications illustrate the versatility of this technique as well as its usefulness in solving similar problems in the aerospace industry. The codes used are the simplified marker and cell (SMAC) code and the variable grid simplified marker and cell code named ERIE. The examples focus on the dropout height in the LH2 tank, the slosh behavior of propellant in the bottom dome of the ET during thrust termination, both of which are required to accurately determine the amount of usable propellants, and the effect on liquid surface due to LO2 tank wall motions occurring during Space Shuttle lift-off. Comparisons with experimental results and linear theory are also presented.

Mehta, G.↗

Modelling of E-region auroral winds

A numerical code, describing the interactions between chemistry and dynamics in the coupled thermosphere and ionosphere systems, has been used to investigate the morphology of wind systems of the auroral E-region and, in particular, their response to geomagnetic forcing. The structure and dynamics of both the thermosphere and the ionosphere show very strong interactions at all levels of geomagnetic disturbance. By comparison, other perturbations, particularly those resulting from tidal and gravity wave sources within the lower atmosphere, which are important and easily identified during quiet periods, are always overwhelmed by the magnitude of the effects of large geomagnetic disturbances above about 110 km altitude. The model simulations will be used to characterize the geomagnetic response of the E-region neutral wind, temperature, and composition in different parts of the polar cap and the auroral oval, and the major changes resulting from seasonal variations.

Rees, David↗

Numerical simulation of self-field MPD thrusters

A fully two dimensional magnetohydrodynamics code was developed to predict self-field, steady-state MPD thruster performance. The governing equations are outlined, and methods of solution are presented. Model predictions are compared with experimental data for two thruster geometries over a range of discharge currents and mass flow rates. Model limitations are evaluated, and issues concerning quasi-steady versus steady-state thruster comparisons are discussed.

Lapointe, Michael R.↗

Transonic wind tunnel wall interference analysis for the Space Shuttle Launch Vehicle

The AEDC perforated-wall boundary condition has been incorporated into the NASA Ames OVERFLOW code and has been used to calculate the tunnel flow field in the AEDC Aerodynamic Wind Tunnel (16T) over a three-percent model of the Space Shuttle Launch Vehicle. Tunnel flow-field calculations have been performed for two supersonic flow-field conditions, including one near sonic condition. Wall interference assessments and analyses have been performed by making detailed comparisons of the tunnel flow-field calculations with corresponding free-air flow-field calculations and with tunnel data from different scale models. A capability to assess wall interference at low supersonic conditions and to improve data correlations with CFD results by adding the tunnel wall boundary condition is demonstrated.

Martin, F. W., Jr.↗

Horizontal axis wind turbine post stall airfoil characteristics synthesization

Blade-element/momentum performance prediction codes are routinely used for wind turbine design and analysis. A weakness of these codes is their inability to consistently predict peak power upon which the machine structural design and cost are strongly dependent. The purpose of this study was to compare post-stall airfoil characteristics synthesization theory to a systematically acquired wind tunnel data set in which the effects of aspect ratio, airfoil thickness, and Reynolds number were investigated. The results of this comparison identified discrepancies between current theory and the wind tunnel data which could not be resolved. Other factors not previously investigated may account for these discrepancies and have a significant effect on peak power prediction.

Tangler, James L.↗

Application of a Navier-Stokes Solver to the Analysis of Multielement Airfoils and Wings Using Multizonal Grid Techniques

A computational study was performed to determine the predictive capability of a Reynolds averaged Navier-Stokes code (CFL3D) for two-dimensional and three-dimensional multielement high-lift systems. Three configurations were analyzed: a three-element airfoil, a wing with a full span flap and a wing with a partial span flap. In order to accurately model these complex geometries, two different multizonal structured grid techniques were employed. For the airfoil and full span wing configurations, a chimera or overset grid technique was used. The results of the airfoil analysis illustrated that although the absolute values of lift were somewhat in error, the code was able to predict reasonably well the variation with Reynolds number and flap position. The full span flap analysis demonstrated good agreement with experimental surface pressure data over the wing and flap. Multiblock patched grids were used to model the partial span flap wing. A modification to an existing patched- grid algorithm was required to analyze the configuration as modeled. Comparisons with experimental data were very good, indicating the applicability of the patched-grid technique to analyses of these complex geometries.

Jones, Kenneth M.↗

A Comparison of Ffowcs Williams-Hawkings Solvers for Airframe Noise Applications

This paper presents a comparison between two implementations of the Ffowcs Williams and Hawkings equation for airframe noise applications. Airframe systems are generally moving at constant speed and not rotating, so these conditions are used in the current investigation. Efficient and easily implemented forms of the equations applicable to subsonic, rectilinear motion of all acoustic sources are used. The assumptions allow the derivation of a simple form of the equations in the frequency-domain, and the time-domain method uses the restrictions on the motion to reduce the work required to find the emission time. The comparison between the frequency domain method and the retarded time formulation reveals some of the advantages of the different approaches. Both methods are still capable of predicting the far-field noise from nonlinear near-field flow quantities. Because of the large input data sets and potentially large numbers of observer positions of interest in three-dimensional problems, both codes utilize the message passing interface to divide the problem among different processors. Example problems are used to demonstrate the usefulness and efficiency of the two schemes.

Lockard, David P.↗

Addressing Control Research Issues Leading to Piloted Simulations in Support of the IFCS F-15

This report summarizes the research effort by a team of researchers at West Virginia University in support of the NASA Intelligent Flight Control System (IFCS) F-15 program. In particular, WVU researchers assisted NASA Dryden researchers in the following technical tasks leading to piloted simulation of the 'Gen_2' IFCS control laws. Task #1- Performance comparison of different neural network (NN) augmentation for the Dynamic Inversion (DI) -based VCAS 'Gen_2' control laws. Task #2- Development of safety monitor criteria for transition to research control laws with and without failure during flight test. Task #3- Fine-tuning of the 'Gen_2' control laws for cross-coupling reduction at post-failure conditions. Matlab/Simulink-based simulation codes were provided to the technical monitor on a regular basis throughout the duration of the project. Additional deliverables for the project were Power Point-based slides prepared for different project meetings. This document provides a description of the methodology and discusses the general conclusions from the simulation results.

Napolitano, Marcello↗

Computational Aerothermodynamic Simulation Issues on Unstructured Grids

The synthesis of physical models for gas chemistry and turbulence from the structured grid codes LAURA and VULCAN into the unstructured grid code FUN3D is described. A directionally Symmetric, Total Variation Diminishing (STVD) algorithm and an entropy fix (eigenvalue limiter) keyed to local cell Reynolds number are introduced to improve solution quality for hypersonic aeroheating applications. A simple grid-adaptation procedure is incorporated within the flow solver. Simulations of flow over an ellipsoid (perfect gas, inviscid), Shuttle Orbiter (viscous, chemical nonequilibrium) and comparisons to the structured grid solvers LAURA (cylinder, Shuttle Orbiter) and VULCAN (flat plate) are presented to show current capabilities. The quality of heating in 3D stagnation regions is very sensitive to algorithm options in general, high aspect ratio tetrahedral elements complicate the simulation of high Reynolds number, viscous flow as compared to locally structured meshes aligned with the flow.

Gnoffo, Peter A.↗

Predicting Failure Progression and Failure Loads in Composite Open-Hole Tension Coupons

Failure types and failure loads in carbon-epoxy [45n/90n/-45n/0n]ms laminate coupons with central circular holes subjected to tensile load are simulated using progressive failure analysis (PFA) methodology. The progressive failure methodology is implemented using VUMAT subroutine within the ABAQUS(TradeMark)/Explicit nonlinear finite element code. The degradation model adopted in the present PFA methodology uses an instantaneous complete stress reduction (COSTR) approach to simulate damage at a material point when failure occurs. In-plane modeling parameters such as element size and shape are held constant in the finite element models, irrespective of laminate thickness and hole size, to predict failure loads and failure progression. Comparison to published test data indicates that this methodology accurately simulates brittle, pull-out and delamination failure types. The sensitivity of the failure progression and the failure load to analytical loading rates and solvers precision is demonstrated.

Arunkumar, Satyanarayana↗

Simplified Vicarious Radiometric Calibration

A measurement-based radiance estimation approach for vicarious radiometric calibration of spaceborne multispectral remote sensing systems has been developed. This simplified process eliminates the use of radiative transfer codes and reduces the number of atmospheric assumptions required to perform sensor calibrations. Like prior approaches, the simplified method involves the collection of ground truth data coincident with the overpass of the remote sensing system being calibrated, but this approach differs from the prior techniques in both the nature of the data collected and the manner in which the data are processed. In traditional vicarious radiometric calibration, ground truth data are gathered using ground-viewing spectroradiometers and one or more sun photometer( s), among other instruments, located at a ground target area. The measured data from the ground-based instruments are used in radiative transfer models to estimate the top-of-atmosphere (TOA) target radiances at the time of satellite overpass. These TOA radiances are compared with the satellite sensor readings to radiometrically calibrate the sensor. Traditional vicarious radiometric calibration methods require that an atmospheric model be defined such that the ground-based observations of solar transmission and diffuse-to-global ratios are in close agreement with the radiative transfer code estimation of these parameters. This process is labor-intensive and complex, and can be prone to errors. The errors can be compounded because of approximations in the model and inaccurate assumptions about the radiative coupling between the atmosphere and the terrain. The errors can increase the uncertainty of the TOA radiance estimates used to perform the radiometric calibration. In comparison, the simplified approach does not use atmospheric radiative transfer models and involves fewer assumptions concerning the radiative transfer properties of the atmosphere. This new technique uses two neighboring uniform ground target areas having different reflectance values. The target areas can be natural or artificial and must be large enough to minimize adjacent-pixel contamination effects. The radiative coupling between the atmosphere and the terrain needs to be approximately the same for the two targets. This condition can be met for relatively uniform backgrounds when the distance between the targets is within a few hundred meters. For each target area, the radiance leaving the ground in the direction of the satellite is measured with a radiometrically calibrated spectroradiometer. Using the radiance measurements from the two targets, atmospheric adjacency and atmospheric scattering effects can be subtracted, thereby eliminating many assumptions about the atmosphere and the radiative interaction between the atmosphere and the terrain. In addition, the radiometrically calibrated spectroradiometer can be used with a known reflectance target to estimate atmospheric transmission and diffuse- to-global ratios without the need for ancillary sun photometers. Several comparisons between the simplified method and traditional techniques were found to agree within a few percent. Hence, the simplified method reduces the overall complexity of performing vicarious calibrations and can serve as a method for validating traditional radiative transfer models

Stanley, Thomas↗

Development and Documentation of kL-Based Linear, Nonlinear, and Full Reynolds Stress Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) turbulence models are reported herein. The kL is 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 kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important di erences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady ows in URANS mode. The present report documents the formulation of three model levels of turbulence models as implemented in the CFD code FUN3D. Methodology and calibration examples are shown in detail. The levels are the linear k-kL and the two-equation algebraic Reynolds stress model (ARSM) as well as the full Reynolds Stress Model (RSM). Attached, separated and corner ow cases are documented and compared with experimental, theoretical and other turbulence model data. The results show generally very good comparisons with canonical and experimental data. The results from this formulation are similar or better than results using the SST two- equation turbulence model. ARSM shows great promise with similar level of computational resources as general two equation turbulence models.

Abdol-Hamid, Khaled S.↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D for Axisymmetric Jet Flow

Two full seven-equation turbulence models available in the FUN3D code are evaluated for their ability to improve the computation of challenging mixing flows encountered in aerospace propulsion. These models are the SSG/LRR and Wilcox full second-moment Reynolds stress models. They solve equations for the six components of the Reynolds stress and a seventh equation for the turbulent length scale. Two standard eddy viscosity models are also evaluated for comparison, the Spalart-Allmaras (SA) one-equation model and the Menter Shear Stress Transport (SST-V) two-equation turbulence model. Flow through an axisymmetric reference nozzle is examined at three flow conditions: subsonic unheated, subsonic heated, and near sonic unheated. Centerline profiles of velocity and turbulent kinetic energy and radial profiles of velocity, turbulent kinetic energy and turbulent stresses are examined. Results showed that the SA model did well at predicting the jet potential core length, but over-mixed the downstream flow, whereas the SST-V model over-predicted the potential core length. The Wilcox-model significantly over-predicted the potential core length and under-predicted the mixing and was not well-suited for the jet flows evaluated, however the SSG/LRR Reynolds stress model did well at predicting the mixing rate and mean velocity for all cases examined.

mixing layers↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D for Axisymmetric Jet Flow

Two full seven-equation turbulence models have been implemented into the FUN3D code to evaluate their ability to improve the computation of challenging flows encountered in aerospace propulsion, including mixing flows. These models are the SSG/LRR and Wilcox full second-moment Reynolds stress models. They solve equations for the six components of the Reynolds stress and a seventh equation for the mixing length. Two standard eddy viscosity models are also evaluated for comparison, the Spalart-Allmaras (SA) one-equation model and the Menter Shear Stress Transport (SST-V) two-equation turbulence model. Flow through an axisymmetric reference nozzle is examined at three flow conditions: subsonic unheated, subsonic heated, and near sonic unheated. Centerline profiles of velocity and turbulent kinetic energy and radial profiles of velocity, turbulent kinetic energy and turbulent stresses are examined. characteristics, no significant changes in the downstream flow behavior compared to the baseline case are observed. Furthermore, the total power consumed by the fans for different incoming flow conditions also remain marginally the same. It is hoped that the results, albeit obtained at very low speeds. would serve as a database for this technologically interesting flow field that has not been explored adequately before.

computational fluid dynamics↗