Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface codes”

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 145 records · Page 8

Rapid State Space Modeling Tool for Rectangular Wing Aeroservoelastic Studies

This paper introduces a modeling and simulation tool for aeroservoelastic analysis of rectangular wings with trailing edge control surfaces. The inputs to the code are planform design parameters such as wing span, aspect ratio and number of control surfaces. A doublet lattice approach is taken to compute generalized forces. A rational function approximation is computed. The output, computed in a few seconds, is a state space aeroservoelastic model which can be used for analysis and control design. The tool is fully parameterized with default information so there is little required interaction with the model developer. Although, all parameters can be easily modified if desired.The focus of this paper is on tool presentation, verification and validation. This process is carried out in stages throughout the paper. The rational function approximation is verified against computed generalized forces for a plate model. A model composed of finite element plates is compared to a modal analysis from commercial software and an independently conducted experimental ground vibration test analysis. Aeroservoelastic analysis is the ultimate goal of this tool. Therefore the flutter speed and frequency for a clamped plate are computed using V-g and V-f analysis. The computational results are compared to a previously published computational analysis and wind tunnel results for the same structure. Finally a case study of a generic wing model with a single control surface is presented. Verification of the state space model is presented in comparison to V-g and V-f analysis. This also includes the analysis of the model in response to a 1-cos gust.

Modeling and Simulation↗

Rapid State Space Modeling Tool for Rectangular Wing Aeroservoelastic Studies

This report introduces a modeling and simulation tool for aeroservoelastic analysis of rectangular wings with trailing-edge control surfaces. The inputs to the code are planform design parameters such as wing span, aspect ratio, and number of control surfaces. Using this information, the generalized forces are computed using the doublet-lattice method. Using Roger's approximation, a rational function approximation is computed. The output, computed in a few seconds, is a state space aeroservoelastic model which can be used for analysis and control design. The tool is fully parameterized with default information so there is little required interaction with the model developer. All parameters can be easily modified if desired. The focus of this report is on tool presentation, verification, and validation. These processes are carried out in stages throughout the report. The rational function approximation is verified against computed generalized forces for a plate model. A model composed of finite element plates is compared to a modal analysis from commercial software and an independently conducted experimental ground vibration test analysis. Aeroservoelastic analysis is the ultimate goal of this tool, therefore, the flutter speed and frequency for a clamped plate are computed using damping-versus-velocity and frequency-versus-velocity analysis. The computational results are compared to a previously published computational analysis and wind-tunnel results for the same structure. A case study of a generic wing model with a single control surface is presented. Verification of the state space model is presented in comparison to damping-versus-velocity and frequency-versus-velocity analysis, including the analysis of the model in response to a 1-cos gust.

Verification and Validation↗

Potentials on large spacecraft in LEO

A system-oriented computer code is used to predict surface charging due to voltages generated within a satellite operating in the typical dense plasma environment of LEO. The use of this code is demonstrated by predicting the expansion of electric fields onto a kapton surface from a pinhole over a biased conductor in a LEO environment. The results are compared to a more-exact solution and experimental data.

Mandell, M. J.↗

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.↗

Foam on Tile Impact Modeling for the Space Shuttle Program

Following the breakup of the Space Shuttle Columbia during reentry a NASA-wide investigation team was formed to examine the probable damage inflicted on Orbiter Thermal Protection System (TPS) elements by impact of External Tank insulating foam projectiles. Our team was to apply rigorous, physics-based analysis techniques to help determine parameters of interest for an experimental test program, utilize validated codes to investigate the full range of impact scenarios, and use analysis derived models to predict aero-thermal-structural responses to entry conditions. We were to operate on a non-interference basis with the j Team, and were to supply significant findings to that team and to the Orbiter Vehicle Engineering Working Group, being responsive to any solicitations for support from these entities. The authors formed a working sub-group within the larger team to apply the Smooth Particle Hydrodynamics code SPHC to the damage estimation problem. Numerical models of the LI-900 TPS tiles and of the BX-250 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 sub-groups 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 densified lower layer of LI-900 insulation, the Nomex felt Strain Isolation Pad (SIP) mounting layer, and the underlying 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 the Southwest Research Institute (SwRI) 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 ET foam projectiles on the TPS tile systems used on the wings of the orbiter. Tiles used on the Wing Acreage, the Main Landing Gear Door, and the Carrier Panels near the front edge of the wing were modeled. Foam impacts shot for the CAB investigation were modeled, as well as impacts at larger angles, including rapid rotation of the projectile, and with varying foam properties. 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. Rapid (17 rps) rotation 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. Similar models will be applied for other impacting materials, other velocities, and other geometries as part of the Return to Flight process.

Stellingwerf, R. F.↗

The NASA Low-Pressure Turbine Flow Physics Program: A Review

An overview of the NASA Glenn Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The flow in the LPT is unique for the gas turbine. It is characterized by low Reynolds number and high freestream turbulence intensity and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The flow of most interest is on the suction surface, where large losses are generated due to separation. The LPT is a large, multistage, heavy, jet engine component that suffers efficiency degradation between takeoff and cruise conditions due to decrease in Reynolds number with altitude. The performance penalty is around 2 points for large commercial bypass engines and as much as 7 points for small, high cruise altitude, military engines. The gas-turbine industry is very interested in improving the performance of the LPT and in reducing its weight, part count and cost. Many improvements can be accomplished by improved airfoil design, mainly by increasing the airfoil loading that can yield reduction of airfoils and improved performance. In addition, there is a strong interest in reducing the design cycle time and cost. Key enablers of the needed improvements are computational tools that can accurately predict LPT flows. Current CFD tools in use cannot yet satisfactorily predict the unsteady, transitional and separated flow in the LPT. The main reasons are inadequate transition & turbulence models and incomplete understanding of the LPT flow physics. NASA Glenn has established its LPT program to answer these needs. The main goal of the program is to develop and assess models for unsteady CFD of LPT flows. An approach that consists of complementing and augmenting experimental and computational work elements has been adopted. The work is performed in-house and by several academic institutions, in cooperation and interaction with industry. The program was reviewed at the Minnowbrook II meeting in 1997. This review will summarize the progress that was made since and will introduce newly started projects. The LPT program is focused on three areas: acquisition of experimental and numerical databases and on modeling and computation. Priority was initially given to experiments. There are three classes of experiments: simulated LPT passages, linear cascade, both with and without wakes, and low-speed rotating rig. They are being conducted as follows: At NASA GRC on a flat surface with blade pressure distribution, at the US Naval Academy on a curved surface. The addition of wakes is studied at the University of Minnesota in a curved passage with a retractable wake generator, and at Texas A&M University in a linear cascade with continuously running wake generator. The pressure distribution of the Pratt & Whitney blade 'Pak B' is used in all these experiments. Experiments have been performed also in the GEAE Low-Speed Rotating Turbine (LSRT) rig with GE-designed airfoils. Work on numerically generated database is in progress at the University of Kentucky, using the DNS/LES code LESTool developed there. Turbulence/transition model assessment and development is performed also at the University of Kentucky, where a new intermittency transport model was developed and many experimental test cases have been numerically computed. Assessments of models using simulations of multistage LPT experiments were performed at Virginia Commonwealth University using the Corsair code. Work on suction surface separation delay, using passive and active flow-control, has also been initiated. Following the overview, Principal Investigators attending the workshop will present in detail several of the projects supported by NASA.

Ashpis, David E.↗

Analysis of a theoretically optimized transonic airfoil

Numerical optimization was used in conjunction with an inviscid, full potential equation, transonic flow analysis computer code to design an upper surface contour for a conventional airfoil to improve its supercritical performance. The modified airfoil was tested in a compressible flow wind tunnel. The modified airfoil's performance was evaluated by comparison with test data for the baseline airfoil and for an airfoil developed by optimization of leading edge of the baseline airfoil. While the leading edge modification performed as expected, the upper surface re-design did not produce all of the expected performance improvements. Theoretical solutions computed using a full potential, transonic airfoil code corrected for viscosity were compared to experimental data for the baseline airfoil and the upper surface modification. These correlations showed that the theory predicted the aerodynamics of the baseline airfoil fairly well, but failed to accurately compute drag characteristics for the upper surface modification.

Lores, M. E.↗

Research on Streamlines and Aerodynamic Heating for Unstructured Grids on High-Speed Vehicles

Engineering codes are needed which can calculate convective heating rates accurately and expeditiously on the surfaces of high-speed vehicles. One code which has proven to meet these needs is the Langley Approximate Three-Dimensional Convective Heating (LATCH) code. It uses the axisymmetric analogue in an integral boundary-layer method to calculate laminar and turbulent heating rates along inviscid surface streamlines. It requires the solution of the inviscid flow field to provide the surface properties needed to calculate the streamlines and streamline metrics. The LATCH code has been used with inviscid codes which calculated the flow field on structured grids, Several more recent inviscid codes calculate flow field properties on unstructured grids. The present research develops a method to calculate inviscid surface streamlines, the streamline metrics, and heating rates using the properties calculated from inviscid flow fields on unstructured grids. Mr. Chris Riley, prior to his departure from NASA LaRC, developed a preliminary code in the C language, called "UNLATCH", to accomplish these goals. No publication was made on his research. The present research extends and improves on the code developed by Riley. Particular attention is devoted to the stagnation region, and the method is intended for programming in the FORTRAN 90 language.

DeJarnette, Fred R.↗

A finite-element-based perturbation model for the rotordynamic analysis of shrouded pump impellers: Part 2: User's guide

This report describes the computational steps involved in executing a finite-element-based perturbation model for computing the rotor dynamic coefficients of a shrouded pump impeller or a simple seal. These arise from the fluid/rotor interaction in the clearance gap. In addition to the sample cases, the computational procedure also applies to a separate category of problems referred to as the 'seal-like' category. The problem, in this case, concerns a shrouded impeller, with the exception that the secondary, or leakage, passage is totally isolated from the primary-flow passage. The difference between this and the pump problem is that the former is analytically of the simple 'seal-like' configuration, with two (inlet and exit) flow-permeable stations, while the latter constitutes a double-entry / double-discharge flow problem. In all cases, the problem is that of a rotor clearance gap. The problem here is that of a rotor excitation in the form of a cylindrical whirl around the housing centerline for a smooth annular seal. In its centered operation mode, the rotor is assumed to give rise to an axisymmetric flow field in the clearance gap. As a result, problems involving longitudinal or helical grooves, in the rotor or housing surfaces, go beyond the code capabilities. Discarding, for the moment, the pre- and post-processing phases, the bulk of the computational procedure consists of two main steps. The first is aimed at producing the axisymmetric 'zeroth-order' flow solution in the given flow domain. Detailed description of this problem, including the flow-governing equations, turbulence closure, boundary conditions, and the finite-element formulation, was covered by Baskharone and Hensel. The second main step is where the perturbation model is implemented, with the input being the centered-rotor 'zeroth-order' flow solution and a prescribed whirl frequency ratio (whirl frequency divided by the impeller speed). The computational domain, in the latter case, is treated as three dimensional, with the number of computational planes in the circumferential direction being specified a priori. The reader is reminded that the deformations in the finite elements are all infinitesimally small because the rotor eccentricity itself is a virtual displacement. This explains why we have generically termed the perturbation model the 'virtually' deformable finite-element category. The primary outcome of implementing the perturbation model is the tangential and radial components, F(sub theta)(sup *) and F(sub r)(sup *) of the fluid-exerted force on the rotor surface due to the whirling motion. Repetitive execution of the perturbation model subprogram over a sufficient range of whirl frequency ratios, and subsequent interpolation of these fluid forces, using the least-square method, finally enable the user to compute the impeller rotor dynamic coefficients of the fluid/rotor interaction. These are the direct and cross-coupled stiffness, damping, and inertia effects of the fluid/rotor interaction.

Baskharone, Erian A.↗

Numerical computation of space shuttle heating and surface streamlines

Exact inviscid flow-field codes are used together with a quasi-three-dimensional boundary-layer analysis to provide estimates of the windward surface heating and streamline patterns of the shuttle orbiter vehicle. The accuracy and limitations of the methods are established by comparison with available wind-tunnel experiments and with more exact numerical solutions for simple flows. Flight predictions are presented showing the effects of finite-rate (nonequilibrium) chemical reactions, and the effects of varying boundary-layer edge conditions due to the growth of the boundary-layer into the inviscid flow (entropy layer swallowing). Differences between flow-field predictions at wind-tunnel and nominal flight conditions are discussed.

Rakich, J. V.↗

Coastal Zone Color Scanner (CZCS): Imagery of near-surface phytoplankton pigment concentrations from the first coastal ocean dynamics experiment (CODE-1), March - July 1981

As part of the first Coastal Ocean Dynamics Experiment, images of ocean color were collected from late March until late July, 1981, by the Coastal Zone Color Scanner aboard Nimbus-7. Images that had sufficient cloud-free area to be of interest were processed to yield near-surface phytoplankton pigment concentrations. These images were then remapped to a fixed equal-area grid. This report contains photographs of the digital images and a brief description of the processing methods.

Abbott, M. R.↗

Forced Response Aeromechanics Analysis in MATLAB®-Based Environment Code With Application to Distortion-Tolerant Fan R24 Blade Geometry

Researchers at the NASA Glenn Research Center have developed a research-level code for forced response analysis for turbomachinery that provides a streamlined framework for aeromechanics analysis, as well as for generating Goodman diagrams. The Forced Response Aeromechanics Analysis in a MATLAB-Based Environment (FRAAME, Version 1) code is designed to accept blade surface unsteady pressure time histories generated by three-dimensional (3D) Unsteady Reynolds-Averaged Navier-Stokes (URANS) code TURBO and modal displacements, modal stresses, and static stresses generated via expanded ANSYS® (Ansys, Inc.) cyclic symmetry analysis. The code’s looping structure allows for largescale analyses, including many blade surface pressure files, modal displacements, and modal stress files for full annulus forced response analysis, including a modal summation method for multiple engine orders and modes. This code is applied to the R24 fan geometry for the Boundary Layer Ingesting Inlet/Distortion-Tolerant Fan (BLI2DTF) project, a propulsion system being developed to increase the fuel efficiency of future aircraft, to serve as validation of experimental data as well as external simulation results. This tool’s functional, modular form is intended to allow users to make modifications easily. Those modifications can include requests for information for any step in the analysis process, as well as adding various functions to compute additional information of interest. It is also written, given sufficient input information, as a general blade geometry forced response solver that is not necessarily specific to the R24 geometry but is used here for verification of code functionality and accuracy. Although forced response analysis tools are also available in commercially available code packages, the post-processing FRAAME code offers reliable, fast, and scalable forced response analysis and Goodman diagram generation for large cases utilizing, specifically, TURBO and ANSYS® results information in an effort to streamline forced response and high cycle fatigue analyses.

Aeromechanics↗

Spacecraft Solar Particle Event (SPE) Shielding: Shielding Effectiveness as a Function of SPE model as Determined with the FLUKA Radiation Transport Code

Analysis of both satellite and surface neutron monitor data demonstrate that the widely utilized Exponential model of solar particle event (SPE) proton kinetic energy spectra can seriously underestimate SPE proton flux, especially at the highest kinetic energies. The more recently developed Band model produces better agreement with neutron monitor data ground level events (GLEs) and is believed to be considerably more accurate at high kinetic energies. Here, we report the results of modeling and simulation studies in which the radiation transport code FLUKA (FLUktuierende KAskade) is used to determine the changes in total ionizing dose (TID) and single-event environments (SEE) behind aluminum, polyethylene, carbon, and titanium shielding masses when the assumed form (i. e., Band or Exponential) of the solar particle event (SPE) kinetic energy spectra is changed. FLUKA simulations have fully three dimensions with an isotropic particle flux incident on a concentric spherical shell shielding mass and detector structure. The effects are reported for both energetic primary protons penetrating the shield mass and secondary particle showers caused by energetic primary protons colliding with shielding mass nuclei. Our results, in agreement with previous studies, show that use of the Exponential form of the event

Koontz, Steve↗

Users manual for the NASA Lewis three-dimensional ice accretion code (LEWICE 3D)

A description of the methodology, the algorithms, and the input and output data along with an example case for the NASA Lewis 3D ice accretion code (LEWICE3D) has been produced. The manual has been designed to help the user understand the capabilities, the methodologies, and the use of the code. The LEWICE3D code is a conglomeration of several codes for the purpose of calculating ice shapes on three-dimensional external surfaces. A three-dimensional external flow panel code is incorporated which has the capability of calculating flow about arbitrary 3D lifting and nonlifting bodies with external flow. A fourth order Runge-Kutta integration scheme is used to calculate arbitrary streamlines. An Adams type predictor-corrector trajectory integration scheme has been included to calculate arbitrary trajectories. Schemes for calculating tangent trajectories, collection efficiencies, and concentration factors for arbitrary regions of interest for single droplets or droplet distributions have been incorporated. A LEWICE 2D based heat transfer algorithm can be used to calculate ice accretions along surface streamlines. A geometry modification scheme is incorporated which calculates the new geometry based on the ice accretions generated at each section of interest. The three-dimensional ice accretion calculation is based on the LEWICE 2D calculation. Both codes calculate the flow, pressure distribution, and collection efficiency distribution along surface streamlines. For both codes the heat transfer calculation is divided into two regions, one above the stagnation point and one below the stagnation point, and solved for each region assuming a flat plate with pressure distribution. Water is assumed to follow the surface streamlines, hence starting at the stagnation zone any water that is not frozen out at a control volume is assumed to run back into the next control volume. After the amount of frozen water at each control volume has been calculated the geometry is modified by adding the ice at each control volume in the surface normal direction.

Bidwell, Colin S.↗

Surface relief structures for multiple beam LO generation

Linear and binary holograms for use in heterodyne detection with 10.6 micron imaging arrays are described. The devices match the amplitude and phase of the local oscillator to the received signal and thus maximize the system signal to noise ratio and resolution and minimize heat generation on the focal plane. In both the linear and binary approaches, the holographic surface-relief pattern is coded to generate a set of local oscillator beams when the relief pattern is illuminated by a single planewave. Each beam of this set has the same amplitude shape distribution as, and is collinear with, each single element wavefront illuminating array.

Veldkamp, W. B.↗

Interactive Graphics Analysis for Aircraft Design

Program uses higher-order far field drag minimization. Computer program WDES WDEM preliminary aerodynamic design tool for one or two interacting, subsonic lifting surfaces. Subcritical wing design code employs higher-order far-field drag minimization technique. Linearized aerodynamic theory used. Program written in FORTRAN IV.

Townsend, J. C.↗

Two-phase flows and heat transfer within systems with ambient pressure above the thermodynamic critical pressure

In systems where the design inlet and outlet pressures P sub amb are maintained above the thermodynamic critical pressure P sub c, it is often assumed that heat and mass transfer are governed by single-phase relations and that two-phase flows cannot occur. This simple rule of thumb is adequate in many low-power designs but is inadequate for high-performance turbomachines, boilers, and other systems where two-phase regions can exist even though P sub amb P sub c. Heat and mass transfer and rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two-phase zone can differ significantly from those for a single-phase zone. By using a laminar, variable-property bearing code and a rotating boiler code, pressure and temperature surfaces were determined that illustrate nesting of a two-phase region within a supercritical pressure region. The method of corresponding states is applied to bearings with reasonable rapport.

Hendricks, R. C.↗

Two-phase flows and heat transfer within systems with ambient pressure above the thermodynamic critical pressure

In systems where the design inlet and outlet pressure P sub amb are maintained above the thermodynamic critical pressure P sub c, it is often assumed that heat and mass transfer are governed by single-phase relations and that two-phase flows cannot occur. This simple rule of thumb is adequate in many low-power designs but is inadequate for high-performance turbomachines, boilers, and other systems where two-phase regions can exist even though P sub amb P sub c. Heat and mass transfer and rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two-phase zone can differ significantly from those for a single-phase zone. By using a laminar, variable-property bearing code and a rotating boiler code, pressure and temperature surfaces were determined that illustrate nesting of a two-phase region within a supercritical pressure region. The method of corresponding states is applied to bearings with reasonable rapport.

Hendricks, R. C.↗