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 883 records · Page 49

Euler and Navier-Stokes solutions for the leeside flow over delta wings at supersonic speeds

Distinctly different types of leeside flowfields over highly swept sharp leading edge delta wings in supersonic flow were numerically simulated using Euler and Navier-Stokes solvers. The Euler code was seen to be adequate only in predicting primary flow structures (leading edge vortex and cross flow shock) whereas the Navier-Stokes code was capable of predicting secondary flow structures (i.e., secondary vortex). A comparison of laminar and turbulent Navier-Stokes solutions indicated that the turbulent boundary layer model is more accurate in predicting the effect of the boundary layer model on the flowfield. Also, the Navier-Stokes code indicated detailed flow structures not observed in the qualitative experimental data available (i.e., vapor screen photographs) indicating a need for quantitative flow field data.

Mcmillin, S. N.↗

Nonlinear lift control at high speed and high angle of attack using vortex flow technology

Nonlinear lift control at subsonic, transonic and low supersonic speeds owes its origin to the separated but organized vortical flows interacting with the wing upper surface. Since most of this flow originates near the wing or control-surface leading edge, a variety of devices has been studied experimentally which interact with and/or control this flow in order to gain a beneficial effect. The benefits (effects) originally studied were associated only with lift enhancement. Whereas, now the studied benefits encompass performance increase, attention to changes in trimmed conditions and longitudinal stability, improvements in lateral stability, and the attendant variation with changing Mach number. For those devices that can be theoretically modeled, state-of-the-art computer codes have been used for device design and/or analysis. Comparisons at design and off-design conditions are presented for validation purposes.

Lamar, John E.↗

Introduction to the analysis of delamination related fracture processes in composites

This research concerns the analysis and prediction of delamination damage that occurs in composite structures on the sublaminate scale - that is, the scale of individual plies or groups of plies. The objective was to develop analytical models for fixed-mode delamination in composites. These include: (1) the influence of residual thermal and moisture strains; (2) local or transverse crack tip delamination originating at the tip of transverse matrix cracks; and (3) delamination in tapered composite under tensile loading. Computer codes based on the analytical models were developed and comparisons of predictions with available experimental and analytical results in the literature were performed.

Armanios, Erian A.↗

An efficient method for computing unsteady transonic aerodynamics of swept wings with control surfaces

A transonic equivalent strip (TES) method was further developed for unsteady flow computations of arbitrary wing planforms. The TES method consists of two consecutive correction steps to a given nonlinear code such as LTRAN2; namely, the chordwise mean flow correction and the spanwise phase correction. The computation procedure requires direct pressure input from other computed or measured data. Otherwise, it does not require airfoil shape or grid generation for given planforms. To validate the computed results, four swept wings of various aspect ratios, including those with control surfaces, are selected as computational examples. Overall trends in unsteady pressures are established with those obtained by XTRAN3S codes, Isogai's full potential code and measured data by NLR and RAE. In comparison with these methods, the TES has achieved considerable saving in computer time and reasonable accuracy which suggests immediate industrial applications.

Liu, D. D.↗

Analysis of generic scramjet external nozzle flowfields employing simulant gases

Experimental and computational results are presented for a generic scramjet nozzle/afterbody model. Mixtures of argon and Freon were employed as simulant scramjet exhaust gases because they correctly model the inviscid simulation parameters of actual scramjet combustion products, allowing test and analysis at significantly reduced temperatures. Air was also employed as an exhaust gas to provide experimental and computational comparisons with simulant gas effects. Several computational aerodynamic analysis codes, employing different levels of theoretical modeling, were applied to compute solutions about the scramjet geometry. The solutions are compared with measured surface static and flowfield pitot pressure data taken in a Mach 6 freestream.

Tatum, Kenneth E.↗

Navier-Stokes analysis of a liquid rocket engine disk cavity

This paper presents a Navier-Stokes analysis of hydrodynamic phenomena occurring in the aft disk cavity of a liquid rocket engine turbine. The cavity analyzed in the Space Shuttle Main Engine Alternate Turbopump currently being developed by NASA and Pratt and Whitney. Comparison of results obtained from the Navier-Stokes code for two rotating disk datasets available in the literature are presented as benchmark validations. The benchmark results obtained using the code show good agreement relative to experimental data, and the turbine disk cavity was analyzed with comparable grid resolution, dissipation levels, and turbulence models. Predicted temperatures in the cavity show that little mixing of hot and cold fluid occurs in the cavity and the flow is dominated by swirl and pumping up the rotating disk.

Benjamin, Theodore G.↗

Evaluation of MHOST analysis capabilities for a plate element

Results of the evaluation of the static, buckling, and free vibration analyses capabilities of MHOST for the plate elements are presented. Two large scale, general purpose finite element codes (MARC and MSC/NASTRAN) are used to validate MHOST. Comparisons of MHOST results with those from MARC and MSC/NASTRAN show good agreement and indicate that MHOST can be used with confidence to perform the aforementioned analyses using the plate element.

Lee, Ho-Jun↗

Study of unsteady flow field over a forward-looking endoatmospheric hit-to-kill interceptor

Forward-looking recessed aperture interceptor has significant aero-optical and aero-thermal advantages. Previous experimental studies have shown that the flow field in front of a forward-looking cavity is unsteady and the bow shock oscillates at the cavity fundamental resonant frequency. In this study, an advanced CFD code is applied to study the above unsteady phenomena. The code is first validated against the experiments and good comparisons are found. The numerical parametric study shows that the existence of oscillatory bow shock is very sensitive to the cavity geometry. At a FOV of 60 deg, the initial transient quickly dampens out to a steady state. With a decrease of FOV, an unsteady oscillatory flow field is sustained after initial transient and the amplitude of oscillation is a function of FOV. For FOV of 20 deg, the amplitude of pressure oscillation is 25 percent of the mean value in the cavity. For a FOV of 10 deg, it can be as high as 50 percent.

Yang, H. Q.↗

Computing Plasma Interactions Of A Spacecraft

NASCAP/LEO (NASA Spacecraft Charging Analyzer Program for Low Earth Orbit) implements collection of mathematical models and algorithms designed to study electrostatic interaction between cold, dense plasma and spacecraft surfaces. Computes variety of electrostatic, plasma, and flow effects. Appropriate for conditions in which temperature of plasma small in comparison with spacecraft-generated potentials and Debye screening length short in comparison with dimensions of spacecraft. Related NASCAP/GEO (LEW-12973) code (NASA Charging Analyzer Program for Geosynchronous Orbit, denoted as NASCAP) appropriate for conditions which spacecraft differential potentials result from interactions with hot plasma and Debye screening length larger than dimensions of spacecraft. Object-definition portion of NASCAP/GEO code provided as part of package. NASCAP/LEO written in FORTRAN 77 and C language.

Mandell, M. J.↗

Stress and Damage in Polymer Matrix Composite Materials Due to Material Degradation at High Temperatures

This report describes analytical methods for calculating stresses and damage caused by degradation of the matrix constituent in polymer matrix composite materials. Laminate geometry, material properties, and matrix degradation states are specified as functions of position and time. Matrix shrinkage and property changes are modeled as functions of the degradation states. The model is incorporated into an existing composite mechanics computer code. Stresses, strains, and deformations at the laminate, ply, and micro levels are calculated, and from these calculations it is determined if there is failure of any kind. The rationale for the model (based on published experimental work) is presented, its integration into the laminate analysis code is outlined, and example results are given, with comparisons to existing material and structural data. The mechanisms behind the changes in properties and in surface cracking during long-term aging of polyimide matrix composites are clarified. High-temperature-material test methods are also evaluated.

McManus, Hugh L.↗

Comparison of Columnar Water Vapor Measurements During The Fall 1997 ARM Intensive Observation Period: Solar Transmittance Methods

In the fall of 1997, during an Intensive Observation Period (IOP), the Atmospheric Radiation Measurement (ARM) program conducted a study of water vapor abundance measurement at its Southern Great Plains (SGP) site. Among a large number of instruments, four sun-tracking radiometers were present to measure the columnar water vapor (CWV). All four solar radiometers retrieve CWV by measuring total solar transmittance in the 0.94-gm water vapor absorption band and subtracting contributions due to Rayleigh, ozone and aerosol transmittances. The aerosol optical depth comparisons among the same four radiometers has been presented elsewhere (Geophys. Res. Lett., 26, 17, 2725-2728, 1999). We have used three different methods to retrieve CWV. In a first round of comparison no attempt was made to standardize on the same radiative transfer model and its underlying water vapor spectroscopy. In the second round of comparison we used the same line-by-line code (which includes recently corrected H2O spectroscopy) to retrieve CAN from all four suntracking radiometers. This decreased the mean CWV by 8% or 13%. The spread of 8% in the solar radiometer results found when using the same model is an indication of the other-than-model uncertainties involved in determining CWV from solar transmittance measurements with current instrumentation.

Schmid, B.↗

Improving Metallic Thermal Protection System Hypervelocity Impact Resistance Through Design of Experiments Approach

A design of experiments approach has been implemented using computational hypervelocity impact simulations to determine the most effective place to add mass to an existing metallic Thermal Protection System (TPS) to improve hypervelocity impact protection. Simulations were performed using axisymmetric models in CTH, a shock-physics code developed by Sandia National Laboratories, and validated by comparison with existing test data. The axisymmetric models were then used in a statistical sensitivity analysis to determine the influence of five design parameters on degree of hypervelocity particle dispersion. Several damage metrics were identified and evaluated. Damage metrics related to the extent of substructure damage were seen to produce misleading results, however damage metrics related to the degree of dispersion of the hypervelocity particle produced results that corresponded to physical intuition. Based on analysis of variance results it was concluded that the most effective way to increase hypervelocity impact resistance is to increase the thickness of the outer foil layer. Increasing the spacing between the outer surface and the substructure is also very effective at increasing dispersion.

Poteet, Carl C.↗

In Search of the Physics: NASA's Approach to Airframe Noise

An extensive numerical and experimental study of airframe noise mechanisms associated with a subsonic high-lift system has been performed at NASA Langley Research Center (LaRC). Investigations involving both steady and unsteady computations and experiments on small-scale models with part-span flaps and full-span flaps are presented. Both surface (steady and unsteady pressure measurements, hot films, oil flows, pressure sensitive paint) and off-surface (5 holeprobe, particle-imaged velocimetry, laser velocimetry, laser light sheet measurements) were taken in the LaRC Quiet Flow Facility (QFF) and several hard-wall tunnels. Experiments in the Low Turbulence Pressure Tunnel (LTPT) included Reynolds number variations up to flight conditions. Successful microphone array measurements were also taken providing both acoustic source maps on the model, and quantitative spectra. Critical directivity measurements were obtained in the QFF. NASA Langley unstructured and structured Reynolds-Averaged Navier-Stokes codes modeled the steady aspects of the flows. Excellent comparisons with surface and off-surface experimental data were obtained. Subsequently, these meanflow calculations were utilized in both linear stability and direct numerical simulations of the flow fields to calculate unsteady surface pressures and farfield acoustic spectra. Accurate calculations were critical in obtaining not only noise source characteristics, but shear layer correction data as well. Techniques utilized in these investigations as well as brief overviews of the results are given.

Macaraeg, Michele G.↗

Computations of Viking Lander Capsule Hypersonic Aerodynamics with Comparisons to Ground and Flight Data

Comparisons are made between the LAURA Navier-Stokes code and Viking Lander Capsule hypersonic aerodynamics data from ground and flight measurements. Wind tunnel data are available for a 3.48 percent scale model at Mach 6 and a 2.75 percent scale model at Mach 10.35, both under perfect gas air conditions. Viking Lander 1 aerodynamics flight data also exist from on-board instrumentation for velocities between 2900 and 4400 m/sec (Mach 14 to 23.3). LAURA flowfield solutions are obtained for the geometry as tested or flown, including sting effects at tunnel conditions and finite-rate chemistry effects in flight. Using the flight vehicle center-of-gravity location (trim angle approx. equals -11.1 deg), the computed trim angle at tunnel conditions is within 0.31 degrees of the angle derived from Mach 6 data and 0.13 degrees from the Mach 10.35 trim angle. LAURA Mach 6 trim lift and drag force coefficients are within 2 percent of measured data, and computed trim lift-to-drag ratio is within 4 percent of the data. Computed trim lift and drag force coefficients at Mach 10.35 are within 5 percent and 3 percent, respectively, of wind tunnel data. Computed trim lift-to-drag ratio is within 2 percent of the Mach 10.35 data. Using the nominal density profile and center-of-gravity location, LAURA trim angle at flight conditions is within 0.5 degrees of the total angle measured from on-board instrumentation. LAURA trim lift and drag force coefficients at flight conditions are within 7 and 5 percent, respectively, of the flight data. Computed trim lift-to-drag ratio is within 4 percent of the data. Computed aerodynamics sensitivities to center-of-gravity location, atmospheric density, and grid refinement are generally small. The results will enable a better estimate of aerodynamics uncertainties for future Mars entry vehicles where non-zero angle-of-attack is required.

Edquist, Karl T.↗

Axisymmetric Calculations of a Low-Boom Inlet in a Supersonic Wind Tunnel

This paper describes axisymmetric CFD predictions made of a supersonic low-boom inlet with a facility diffuser, cold pipe, and mass flow plug within wind tunnel walls, and compares the CFD calculations with the experimental data. The inlet was designed for use on a small supersonic aircraft that would cruise at Mach 1.6, with a Mach number over the wing of 1.7. The inlet was tested in the 8-ft by 6-ft Supersonic Wind Tunnel at NASA Glenn Research Center in the fall of 2010 to demonstrate the performance and stability of a practical flight design that included a novel bypass duct. The inlet design is discussed here briefly. Prior to the test, CFD calculations were made to predict the performance of the inlet and its associated wind tunnel hardware, and to estimate flow areas needed to throttle the inlet. The calculations were done with the Wind-US CFD code and are described in detail. After the test, comparisons were made between computed and measured shock patterns, total pressure recoveries, and centerline pressures. The results showed that the dual-stream inlet had excellent performance, with capture ratios near one, a peak core total pressure recovery of 96 percent, and a large stable operating range. Predicted core recovery agreed well with the experiment but predicted bypass recovery and maximum capture ratio were high. Calculations of offdesign performance of the inlet along a flight profile agreed well with measurements and previous calculations.

Chima, Rodrick V.↗

Solar Flare Ribbon Fronts. I. Constraining Flare Energy Deposition with IRIS Spectroscopy

Spectral lines formed at lower atmospheric layers show peculiar profiles at the "leading edge" of ribbons during solar flares. In particular, increased absorption of the BBSO/GST He I λ10830 line, as well as broad and centrally reversed profiles in the spectra of the Mg II and C II lines observed by the IRIS satellite, has been reported. In this work, we aim to understand the physical origin of such peculiar IRIS profiles, which seem to be common of many, if not all, flares. To achieve this, we quantify the spectral properties of the IRIS Mg II profiles at the ribbon leading edge during four large flares and perform a detailed comparison with a grid of radiative hydrodynamic models using the RADYN+FP code. We also studied their transition region (TR) counterparts, finding that these ribbon front locations are regions where TR emission and chromospheric evaporation are considerably weaker compared to other parts of the ribbons. Based on our comparison between the IRIS observations and modeling, our interpretation is that there are different heating regimes at play in the leading edge and the main bright part of the ribbons. More specifically, we suggest that bombardment of the chromosphere by more gradual and modest nonthermal electron energy fluxes can qualitatively explain the IRIS observations at the ribbon leading front, while stronger and more impulsive energy fluxes are required to drive chromospheric evaporation and more intense TR emission in the bright ribbon. Our results provide a possible physical origin for the peculiar behavior of the IRIS chromospheric lines in the ribbon leading edge and new constraints for the flare models.

Solar flares↗

Challenges in Rotor Aerodynamic Modeling for Non-Uniform Inflow Conditions

Within an international collaboration framework, the accuracy of rotor aerodynamic models used for design load calculations of wind turbines is being assessed. Where the use of high-fidelity computation fluid dynamics (CFD) and mid-fidelity free-vortex wake (FVW) models has become commonplace within the wind energy community, these still fail to meet the requirements in terms of execution time and computational cost needed for design load calculations. The fast but engineering fidelity blade-element/momentum (BEM) method can therefore still be considered the industry workhorse for design load simulations. At the same time, upscaling of wind turbine rotors makes inflow non-uniformities (e.g. shear, veer, turbulence) more important. The objective of this work is to assess model accuracy in non-uniform inflow conditions, which violate several BEM assumptions. Thereto a comparison in turbulent inflow has been executed including a wide variety of codes, focusing on the DanAero field measurements, where a 2.3-MW turbine was equipped with, among other sensors, a pressure measurement apparatus. The results indicate that, although average load patterns are in good agreement, this does not hold for the unsteady loads that drive fatigue damage and aero-elastic stability. A simplified comparison round in vertical shear was initiated to investigate the observed differences in a more controlled manner. A consistent offset in load amplitude was observed between CFD and free-vortex codes on the one hand and BEM-type codes on the other hand. To shed more light on the observations, dedicated efforts are ongoing to pinpoint the cause for these differences, in the end leading to guidelines for an improved BEM implementation.

17 WIND ENERGY↗

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