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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Plans and Example Results for the 2nd AIAA Aeroelastic Prediction Workshop

This paper summarizes the plans for the second AIAA Aeroelastic Prediction Workshop. The workshop is designed to assess the state-of-the-art of computational methods for predicting unsteady flow fields and aeroelastic response. The goals are to provide an impartial forum to evaluate the effectiveness of existing computer codes and modeling techniques, and to identify computational and experimental areas needing additional research and development. This paper provides guidelines and instructions for participants including the computational aerodynamic model, the structural dynamic properties, the experimental comparison data and the expected output data from simulations. The Benchmark Supercritical Wing (BSCW) has been chosen as the configuration for this workshop. The analyses to be performed will include aeroelastic flutter solutions of the wing mounted on a pitch-and-plunge apparatus.

Heeg, Jennifer↗
NASA NTRS2016

Simulation of Guided Wave Interaction with In-Plane Fiber Waviness

Reducing the timeline for certification of composite materials and enabling the expanded use of advanced composite materials for aerospace applications are two primary goals of NASA's Advanced Composites Project (ACP). A key a technical challenge area for accomplishing these goals is the development of rapid composite inspection methods with improved defect characterization capabilities. Ongoing work at NASA Langley is focused on expanding ultrasonic simulation capabilities for composite materials. Simulation tools can be used to guide the development of optimal inspection methods. Custom code based on elastodynamic finite integration technique is currently being developed and implemented to study ultrasonic wave interaction with manufacturing defects, such as in-plane fiber waviness (marcelling). This paper describes details of validation comparisons performed to enable simulation of guided wave propagation in composites containing fiber waviness. Simulation results for guided wave interaction with in-plane fiber waviness are also discussed. The results show that the wavefield is affected by the presence of waviness on both the surface containing fiber waviness, as well as the opposite surface to the location of waviness.

Leckey, Cara A. C.↗
NASA NTRS1988

Configuration trade and code validation study on a conical hypersonic vehicle

A study has been conducted on a generic wing-cone transatmospheric vehicle at Mach numbers form 2.5 to 4.5. The objectives of the study were to experimentally define the aerodynamic characteristics of the vehicle and evaluate several computational aerodynamic prediction methods through comparison with the experimental results. The baseline wing-cone configuration fuselage consisted of a 5 deg half-angle cone forebody, cylindrical midbody, and 9 deg truncated cone afterbody. The 4-percent-thick diamond airfoil wing had an aspect ratio of 1. Several configuration variables were investigated to provide trade information on canard, wing-position and incidence, vertical tail, and nose bluntness effects. Results of the study show that wing-position and wing-incidence effects on the longitudinal aerodynamic characteristics can be significantly influenced by wing-body interference. The use of positive wing incidence to provide favorable forebody orientation for possible inlet performance improvement is accompanied by trim drag and lift-drag ratio penalties. The lateral-directional stability characteristics were strongly influenced by the location of the vertical tails. The higher-order full-potential method provided better estimates of the aerodynamic characteristics than either the linearized supersonic potential method or the tangent-cone/tangent-wedge/shock-expansion on method.

Covell, Peter F.↗
NASA NTRS1983

Analytical and Experimental Evaluation of the Heat Transfer Distribution over the Surfaces of Turbine Vanes

Three airfoil data sets were selected for use in evaluating currently available analytical models for predicting airfoil surface heat transfer distributions in a 2-D flow field. Two additional airfoils, representative of highly loaded, low solidity airfoils currently being designed, were selected for cascade testing at simulated engine conditions. Some 2-D analytical methods were examined and a version of the STAN5 boundary layer code was chosen for modification. The final form of the method utilized a time dependent, transonic inviscid cascade code coupled to a modified version of the STAN5 boundary layer code featuring zero order turbulence modeling. The boundary layer code is structured to accommodate a full spectrum of empirical correlations addressing the coupled influences of pressure gradient, airfoil curvature, and free-stream turbulence on airfoil surface heat transfer distribution and boundary layer transitional behavior. Comparison of pedictions made with the model to the data base indicates a significant improvement in predictive capability.

Hylton, L. D.↗
NASA NTRS1984

Inverse transonic airfoil design methods including boundary layer and viscous interaction effects

This report covers the period 1 September 1983 to 31 January 1984. The primary task during this reporting period was the continued development of the massive separation model and computer code (SKANSEP). In particular, detailed investigations were conducted with the boundary layer displacement surface correction technique discovered near the end of the last reporting period. This report will present detailed results using this technique and show comparisons with experimental data.

Carlson, L. A.↗
NASA NTRS1993

A computational investigation of fuel mixing in a hypersonic scramjet

A parabolized, Navier-Stokes code, SHIP3D, is used to numerically investigate the mixing between air injection and hydrogen injection from a swept ramp injector configuration into either a mainstream low-enthalpy flow or a hypervelocity test flow. The mixing comparisons between air and hydrogen injection reveal the importance of matching injectant-to-mainstream mass flow ratios. In flows with the same injectant-to-mainstream dynamic pressure ratio, the mixing definition was altered for the air injection cases. Comparisons of the computed results indicate that the air injection cases overestimate the mixing performance associated with hydrogen injection simulation. A lifting length parameter, to account for the time a fluid particle transverses through the mixing region, is defined and used to establish a connection of injectant mixing in hypervelocity flows, based on nonreactive, low-enthalpy flows.

Fathauer, Brett W.↗
NASA NTRS2001

Two-Equation Turbulence Models for Prediction of Heat Transfer on a Transonic Turbine Blade

Two versions of the two-equation k-omega model and a shear stress transport (SST) model are used in a three-dimensional, multi-block, Navier-Stokes code to compare the detailed heat transfer measurements on a transonic turbine blade. It is found that the SST model resolves the passage vortex better on the suction side of the blade, thus yielding a better comparison with the experimental data than either of the k-w models. However, the comparison is still deficient on the suction side of the blade. Use of the SST model does require the computation of distance from a wall, which for a multiblock grid, such as in the present case, can be complicated. However, a relatively easy fix for this problem was devised. Also addressed are issues such as (1) computation of the production term in the turbulence equations for aerodynamic applications, and (2) the relation between the computational and experimental values for the turbulence length scale, and its influence on the passage vortex on the suction side of the turbine blade.

Garg, Vijay K.↗
NASA NTRS2001

MHD Instability and Turbulence in the Tachocline

In this quarter we have begun simulations on the Cray T3E at PSC and we are debugging our code on the TSC. The PSC simulations are examining stratified shear turbulence with a flow-aligned magnetic field and passive tracer particles. We have conducted analysis of neutral simulations to establish a firm basis of comparison. Second-order structure functions have been computed, fit, and compared to theoretical expressions relating the dissipation fields and the structure-function-fit parameters. Agreement with high-Reynolds number observations is excellent, giving us confidence that the lower-Re simulations are relevant to higher-Re flows. We have also evaluated the neutral layer anisotropy.

Werne, Joseph↗
NASA NTRS1995

Galileo Heat Shield Response Experiment Data Reduction and Analyses

The Galileo heat shield response experiment consists of 10 Analog Resistance Ablation Detectors (ARADs) embedded in the carbon phenolic heat shield. As the vehicle descends into the atmosphere of Jupiter, these ARADs will in effect measure the recession as a function of time at various locations on the heat shield. This recession data will be used to reconstruct the time-dependent shape of the vehicle. The shape data is of critical importance to the Atmospheric Reconstruction Experiment which must know mass and drag coefficient of the vehicle in order to determine the structure of the Jovian atmosphere. The data is also intrinsically useful as a database for evaluating the accuracy of our numerical codes and methodology for predicting highly radiating reentry flowfields and coupled heat shield material response. This paper will document the reduction of the data from the ARAD sensors, reconstruction of the heat shield shape history, and comparison of the actual heat shield response with previous calculations and experiments. If possible, we will also include results from new calculations which couple the shock layer flow with the transient material response.

Milos, Frank S.↗
NASA NTRS2004

Resolving the Issues with Flywheel Position Sensors

For the past few years, the Advanced Electrical Systems Branch here at NASA Glenn has been pursuing research in the area of flywheels. The purpose of these pursuits has been t o explore the potential for flywheels to replace current battery-powered systems in space. So far it has been learned that flywheels offer large momentum storage capacity, comparatively small volume, high durability, and near- complete discharge capabilities, all of which are advancements over the existing nickel hydrogen and nickel cadmium batteries. Another significant advantage of flywheels is the potential they offer for combining the function of attitude control with energy storage. During the summer of 2004, I worked with Dr. Barbara Kenny in the Advanced Electrical Systems Branch, supporting the work she is doing by analyzing and testing some new components for the new Generation-2 flywheel. To monitor the speed and angular position of the flywheel rotor, a once-around (OAR) signal along with a sensorless algorithm is used. The OAR signal is used for the magnetic bearings that keep the flywheel suspended for frictionless operation. The sensorless algorithm is used for the flywheel motor/generator control. The OAR is generated from position sensors that monitor a circular plate. The plate has a cut down the middle such that one half of the circle is on a slightly lower level than the other. Every half-turn, or 180, the sensors detect the "cut" on the plate, and trigger the OAR, telling the computer that the rotor has made half a revolution. This, however, doesn't provide needed detailed information about the angular position of the rotor, since it only provides a signal alert every half- revolution. This is enough information for the magnetic bearing control but is insufficient for the motor/generator control. A new resolver was designed such that it would give continuous angle information rather than the 180 degree information of the OAR. The new resolver has two separate observable pieces: a flat middle section to monitor vertical motion, and an angled section around the circumference, which, when observed from above, produces a sine-wave displacement through the entire 360" revolution. My first job when I arrived this summer was to calibrate the sensors that would be mounted on the inside of the flywheel casing to monitor the position (angular and vertical) of the shaft. After calibration, I used the sensors to evaluate voltage outputs created by position differences between two pairs of sensors on the angled portion of the resolver for eight different angular positions, moving the resolver vertically and laterally through its entire potential range of motion. The results of these tests will be used to determine the rotor angular (and axial) position from the sensor readings once the new flywheel unit is assembled. The sensorless algorighm mentioned above consists of two operations: the signal injection method and the back electro-motive force (EMF). The signal injection is meant to work at low speeds, while the back EMF algorithm is meant to work at higher speeds. Both work together to determine the correct estimate of rotor position and speed based on the measured motor/generator current. It was determined that we wanted to know exactly how accurate our estimation methods were, and so a resolver (a commercially available mechanical sensor mounted to the motor/generator shaft to measure rotor position and speed) and a "Resolver to Digital" (R2D) circuit board was purchased to make the comparison to the existing estimation. My work related to the R2D board has included the following: creating two connector cables (one to power the circuit and one to get readable output off the board), writing Simulink code to process the board's output, and building a dSpace panel to control and monitor the circuit. The next step in the process will be to perform tests to compare the estimated rotor position and speed from the sensorless algorithm to the actual rotor and speed from the resolver signal.

Fehrmann, Elizabeth A.↗
NASA NTRS2005

Effect of Impact Location on the Response of Shuttle Wing Leading Edge Panel 9

The objective of this paper is to compare the results of several simulations performed to determine the worst-case location for a foam impact on the Space Shuttle wing leading edge. The simulations were performed using the commercial non-linear transient dynamic finite element code, LS-DYNA. These simulations represent the first in a series of parametric studies performed to support the selection of the worst-case impact scenario. Panel 9 was selected for this study to enable comparisons with previous simulations performed during the Columbia Accident Investigation. The projectile for this study is a 5.5-in cube of typical external tank foam weighing 0.23 lb. Seven locations spanning the panel surface were impacted with the foam cube. For each of these cases, the foam was traveling at 1000 ft/s directly aft, along the orbiter X-axis. Results compared from the parametric studies included strains, contact forces, and material energies for various simulations. The results show that the worst case impact location was on the top surface, near the apex.

Lyle, Karen H.↗
NASA NTRS2011

Plans for Aeroelastic Prediction Workshop

This paper summarizes the plans for the first Aeroelastic Prediction Workshop. The workshop is designed to assess the state of the art of computational methods for predicting unsteady flow fields and aeroelastic response. The goals are to provide an impartial forum to evaluate the effectiveness of existing computer codes and modeling techniques, and to identify computational and experimental areas needing additional research and development. Three subject configurations have been chosen from existing wind tunnel data sets where there is pertinent experimental data available for comparison. For each case chosen, the wind tunnel testing was conducted using forced oscillation of the model at specified frequencies

Heeg, Jennifer↗
NASA NTRS2014

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Center's Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗
NASA NTRS2014

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Centers Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗
NASA NTRS2019

Weighted Least-Squares Cell-Average Gradient Construction Methods for the VULCAN-CFD Second-Order Accurate Unstructured-Grid Cell-Centered Finite-Volume Solver

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. The construction of cell-average gradients using a weighted linear least-squares method and the use of these gradients in the construction of the inviscid fluxes is the focus of this paper. A comparison of least-squares stencil construction methodologies is presented and approaches to augment the number of cells participating in the stencil while preserving accuracy are explored. Due to our interest in hypersonic flow, a robust multidimensional cell-average gradient limiter procedure that is consistent with the stencil used to construct the cell-average gradients is described and investigated. Canonical problems are computed to illustrate the challenges and investigate the accuracy, robustness and convergence behavior of the cell-average gradient methods on unstructured cell-centered finite-volume grids. Finally, thermally perfect, chemically frozen, Mach 8 turbulent flow of air around a blunt wedge is computed to demonstrate the robustness and convergence behavior of the new method for constructing stencils of use in a weighted linear least-squares gradient method for a hypersonic flow.

White, Jeffery A.↗
NASA NTRS2019

Efficient and Robust Weighted Least-Squares Cell-Average Gradient Construction Methods for the Simulation of Scramjet Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. The construction of cell-average gradients using a weighted linear least-squares method and the use of these gradients in the construction of the inviscid fluxes is the focus of this paper. A comparison of least-squares stencil construction methodologies is presented and approaches designed to minimize the number of cells used to augment/stabilize the least-squares stencil while preserving accuracy are explored. Due to our interest in hypersonic flow, a robust multidimensional cell-average gradient limiter procedure that is consistent with the stencil used to construct the cellaverage gradients is described. Canonical problems are computed to illustrate the challenges and investigate the accuracy, robustness and convergence behavior of the cell-average gradient methods on unstructured cell-centered finite-volume grids. Finally, thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet engine flowpath is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the preferred gradient method for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗
NASA NTRS1999

A Comparison of Single-Cycle Versus Multiple-Cycle Proof Testing Strategies

Single-cycle and multiple-cycle proof testing (SCPT and MCPT) strategies for reusable aerospace propulsion system components are critically evaluated and compared from a rigorous elastic-plastic fracture mechanics perspective. Earlier MCPT studies are briefly reviewed. New J-integral estimation methods for semielliptical surface cracks and cracks at notches are derived and validated. Engineering methods are developed to characterize crack growth rates during elastic-plastic fatigue crack growth (FCG) and the tear-fatigue interaction near instability. Surface crack growth experiments are conducted with Inconel 718 to characterize tearing resistance, FCG under small-scale yielding and elastic-plastic conditions, and crack growth during simulated MCPT. Fractography and acoustic emission studies provide additional insight. The relative merits of SCPT and MCPT are directly compared using a probabilistic analysis linked with an elastic-plastic crack growth computer code. The conditional probability of failure in service is computed for a population of components that have survived a previous proof test, based on an assumed distribution of initial crack depths. Parameter studies investigate the influence of proof factor, tearing resistance, crack shape, initial crack depth distribution, and notches on the MCPT versus SCPT comparison. The parameter studies provide a rational basis to formulate conclusions about the relative advantages and disadvantages of SCPT and MCPT. Practical engineering guidelines are proposed to help select the optimum proof test protocol in a given application.

McClung, R. C.↗
NASA NTRS1996

A Comparison of Single-Cycle Versus Multiple-Cycle Proof Testing Strategies

Single-cycle and multiple-cycle proof testing (SCPT and MCPT) strategies for reusable aerospace propulsion system components are critically evaluated and compared from a rigorous elastic-plastic fracture mechanics perspective. Earlier MCPT studies are briefly reviewed. New J-integral estimation methods for semi-elliptical surface cracks and cracks at notches are derived and validated. Engineering methods are developed to characterize crack growth rates during elastic-plastic fatigue crack growth (FCG) and the tear-fatigue interaction near instability. Surface crack growth experiments are conducted with Inconel 718 to characterize tearing resistance, FCG under small-scale yielding and elastic-plastic conditions, and crack growth during simulated MCPT. Fractography and acoustic emission studies provide additional insight. The relative merits of SCPT and MCPT are directly compared using a probabilistic analysis linked with an elastic-plastic crack growth computer code. The conditional probability of failure in service is computed for a population of components that have survived a previous proof test, based on an assumed distribution of initial crack depths. Parameter studies investigate the influence of proof factor, tearing resistance, crack shape, initial crack depth distribution, and notches on the MCPT vs. SCPT comparison. The parameter studies provide a rational basis to formulate conclusions about the relative advantages and disadvantages of SCPT and MCPT. Practical engineering guidelines are proposed to help select the optimum proof test protocol in a given application.

McClung, R. C.↗