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,603 records · Page 89

On the utility of threads for data parallel programming

Threads provide a useful programming model for asynchronous behavior because of their ability to encapsulate units of work that can then be scheduled for execution at runtime, based on the dynamic state of a system. Recently, the threaded model has been applied to the domain of data parallel scientific codes, and initial reports indicate that the threaded model can produce performance gains over non-threaded approaches, primarily through the use of overlapping useful computation with communication latency. However, overlapping computation with communication is possible without the benefit of threads if the communication system supports asynchronous primitives, and this comparison has not been made in previous papers. This paper provides a critical look at the utility of lightweight threads as applied to data parallel scientific programming.

Fahringer, Thomas↗

Installed Transonic 2D Nozzle Nacelle Boattail Drag Study

The Transonic Nozzle Boattail Drag Study was initiated in 1995 to develop an understanding of how external nozzle transonic aerodynamics effect airplane performance and how strongly those effects are dependent on nozzle configuration (2D vs. axisymmetric). MDC analyzed the axisymmetric nozzle. Boeing subcontracted Northrop-Grumman to analyze the 2D nozzle. AU participants analyzed the AGARD nozzle as a check-out and validation case. Once the codes were checked out and the gridding resolution necessary for modeling the separated flow in this region determined, the analysis moved to the installed wing/body/nacelle/diverter cases. The boat tail drag validation case was the AGARD B.4 rectangular nozzle. This test case offered both test data and previous CFD analyses for comparison. Results were obtained for test cases B.4.1 (M=0.6) and B.4.2 (M=0.938) and compared very well with the experimental data. Once the validation was complete a CFD grid was constructed for the full Ref. H configuration (wing/body/nacelle/diverter) using a combination of patched and overlapped (Chimera) grids. This was done to ensure that the grid topologies and density would be adequate for the full model. The use of overlapped grids allowed the same grids from the full configuration model to be used for the wing/body alone cases, thus eliminating the risk of grid differences affecting the determination of the installation effects. Once the full configuration model was run and deemed to be suitable the nacelle/diverter grids were removed and the wing/body analysis performed. Reference H wing/body results were completed for M=0.9 (a=0.0, 2.0, 4.0, 6.0 and 8.0), M=1.1 (a=4.0 and 6.0) and M=2.4 (a=0.0, 2.0, 4.4, 6.0 and 8.0). Comparisons of the M=0.9 and M=2.4 cases were made with available wind tunnel data and overall comparisons were good. The axi-inlet/2D nozzle nacelle was analyzed isolated. The isolated nacelle data coupled with the wing/body result enabled the interference effects of the installed nacelles to be determined. Isolated nacelle mm were made at M=0.9 and M=1.1 for both the supersonic and transonic nozzle settings. AU of the isolated nacelle cases were run at alpha=0. Full configuration runs were to be made at Mach numbers of 0.9, 1.1, and 2.4 (the same as the wing/body and isolated nacelles). Both the isolated nacelles and installed nacelles were run with inlet conditions designed to give zero spillage. This was to be done in order to isolate the boattail effects as much as possible. Full configuration runs with the supersonic nozzles were completed for M=0.9 and 1.1 at a=4.0 and 6.0 (4 runs total) and with the transonic nozzles at M=0.9 and 1.1 at a=2.0, 4.0 and 6.0 (6 runs total). Drag breakdowns were completed for the M=0.9 and M= 1.1 showing favorable interference drag for both cases.

Malone, Michael B.↗

Structural Modeling of a Five-Meter Thin Film Inflatable Antenna/Concentrator With Rigidized Support Struts

Dynamic characterization of a non-rigidized thin film inflatable antenna/solar concentrator structure with rigidized composite support struts is described in detail. A two-step finite element modeling approach in MSC/NASTRAN is utilized, consisting of: (1) a nonlinear static pressurization procedure used to obtain the updated stiffness matrix, and (2) a modal "restart" eigen solution that uses the modified stiffness matrix. Unique problems encountered in modeling of this large 5-m lightweight inflatable are identified, including considerable difficulty in obtaining convergence in the nonlinear pressurization solution. It was found that the extremely thin polyimide film material (.001 in or I mil) presents tremendous problems in obtaining a converged solution when internal pressure loading is applied. It was concluded that the ratios of film thickness to other geometric dimensions such as torus cross-sectional and ring diameter and lenticular diameter are the critical parameters for convergence of the pressurization procedure. Comparison of finite element predictions for frequency and mode shapes with experimental results indicated reasonable agreement considering the complexity of the structure, the film-to-air interaction, and the nonlinear material properties of the film. It was also concluded that analysis should be done using different finite element to codes to determine if a more robust and stable solution can be obtained.

Smalley, Kurt B.↗

A Green's function method for high charge and energy ion transport

A heavy-ion transport code using Green's function methods is developed. The low-order perturbation terms exhibiting the greatest energy variation are used as dominant energy-dependent terms, and the higher order collision terms are evaluated using nonperturbative methods. The recently revised NUCFRG database is used to evaluate the solution for comparison with experimental data for 625A MeV 20Ne and 517A MeV 40Ar ion beams. Improved agreements with the attenuation characteristics for neon ions are found, and reasonable agreement is obtained for the transport of argon ions in water.

Energy Transfer↗

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

Boundary layer transition↗

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

On the sensitivity of a residual circulation model to differences in input temperature data

The residual mean circulation (RMC) formulation of zonally averaged transport in the middle atmosphere produces a circulation which depends on the distributions of net diabatic heating and temperature. Such circulations are from two temperature data sets, using the same radiative transfer code (Rosenfield et al. 1987). These circulations are then used to transport N2O in a photochemical model. The circulations and the resulting N2O distributions are notably different during the Northern Hemisphere winter, with that based on the NMC temperatures producing too much upward transport in the tropical stratosphere, as judged by comparison with the stratospheric and mesoscale sounder data. The experiment demonstrates that model calculations, in general, and perturbation assessments, in particular, are likely to be quite sensitive to the choice of input temperature data (where this is not computed self-consistently). It also reveals what appears to be a seasonally dependent bias in NMC zonally averaged temperatures with respect to those obtained from the LIMS instrument during 1978/1979.

Guthrie, Paul D.↗

A Mach 6 external nozzle experiment with Argon-Freon exhaust simulation

A scramjet exhaust simulation technique for hypersonic wind tunnel testing has been developed. Mixtures of Argon and Freon correctly match the inviscid simulation parameters of Mach number, static-pressure ratio, and the ratio of specific heats at the combustor exit location; this simulation is accomplished at significantly reduced temperatures and without combustion. An investigation of nozzle parametrics in a Mach 6 freestream showed that the external nozzle ramp angle, the cowl trailing-edge angle, an external nozzle flow fence and the nozzle static-pressure ratio significantly affected the external nozzle thrust and pitching moment as measured by the integration of surface-pressure data. A comparison of Argon-Freon and air exhaust simulation showed that the external nozzle thrust and pitching moment were in error by roughly a factor of 2 using air due to the incorrect match of the ratio of specific heats. An assessment of two-dimensional Euler and Navier-Stokes codes for predicting external nozzle aerodynamic characteristics was made by comparing computed and experimental results.

Pittman, James L.↗

A comparison between theoretical prediction and experimental measurement of the dynamic behavior of spur gears

A comparison was made between computer model predictions of gear dynamics behavior and experimental results. The experimental data were derived from the NASA gear noise rig, which was used to record dynamic tooth loads and vibration. The experimental results were compared with predictions from the DSTO Aeronautical Research Laboratory's gear dynamics code for a matrix of 28 load speed points. At high torque the peak dynamic load predictions agree with the experimental results with an average error of 5 percent in the speed range 800 to 6000 rpm. Tooth separation (or bounce), which was observed in the experimental data for light torque, high speed conditions, was simulated by the computer model. The model was also successful in simulating the degree of load sharing between gear teeth in the multiple tooth contact region.

Rebbechi, Brian↗

Preliminary Results of NASA's First Autonomous Formation Flying Experiment: Earth Observing-1 (EO-1)

NASA's first autonomous formation flying mission is completing a primary goal of demonstrating an advanced technology called enhanced formation flying. To enable this technology, the Guidance, Navigation, and Control center at the Goddard Space Flight Center has implemented an autonomous universal three-axis formation flying algorithm in executive flight code onboard the New Millennium Program's (NMP) Earth Observing-1 (EO-1) spacecraft. This paper describes the mathematical background of the autonomous formation flying algorithm and the onboard design and presents the preliminary validation results of this unique system. Results from functionality assessment and autonomous maneuver control are presented as comparisons between the onboard EO-1 operational autonomous control system called AutoCon(tm), its ground-based predecessor, and a stand-alone algorithm.

Folta, David↗

Full Scale Rotor Aeroacoustic Predictions and the Link to Model Scale Rotor Data

The NASA Aeroacoustic Prediction System (NAPS) is used to establish a link between model-scale and full-scale rotor predictions and is partially validated against measured wind tunnel and flight aeroacoustic data. The prediction approach of NAPS couples a comprehensive rotorcraft analysis with acoustic source noise and propagation codes. The comprehensive analysis selected for this study is CAMRAD-II, which provides the performance/trim/wake solution for a given rotor or flight condition. The post-trim capabilities of CAMRAD-II are used to compute high-resolution sectional airloads for the acoustic tone noise analysis, WOPMOD. The tone noise is propagated to observers on the ground with the propagation code, RNM (Rotor Noise Model). Aeroacoustic predictions are made with NAPS for an isolated rotor and compared to results of the second Harmonic Aeroacoustic Rotor Test (HART-II) program, which tested a 40% dynamically and Mach-scaled BO-105 main rotor at the DNW. The NAPS is validated with comparisons for three rotor conditions: a baseline condition and two Higher Harmonic Control (HHC) conditions. To establish a link between model and full-scale rotor predictions, a full-scale BO-105 main rotor input deck for NAPS is created from the 40% scale rotor input deck. The full-scale isolated rotor predictions are then compared to the model predictions. The comparisons include aerodynamic loading, acoustic levels, and acoustic pressure time histories for each of the three conditions. With this link established, full-scale predictions are made for a range of descent flight conditions and compared with measured trends from the recent Rotorcraft Operational Noise Abatement Procedures (RONAP) flight test conducted by DLR and ONERA. Additionally, the effectiveness of two HHC conditions from the HART-II program is demonstrated for the full-scale rotor in flight.

Boyd, D. Douglas, Jr.↗

Arcjet Semi-Elliptic Nozzle Simulations and Validation in Support of TPS Testing

Much of the ground based testing of advanced thermal protection system (TPS) components for the X33 program is done in arc-heated wind tunnels such as those located in the Arc-Jet Complex at NASA Ames Research Center. These facilities are capable of simulating the high temperature, chemically reacting flow environment experienced by the vehicle during flight. This allows one to test critical design issues such as maximum reuse temperatures, seals, gaps, and increases in heating due to interfaces between different materials. Computational fluid dynamics (CFD) has evolved to the point where it now can be used in the vehicle design process for accurate and timely prediction of trajectory based aerothermal heating environments for re-entry vehicles. It can also be used for simulation of the flow environments in ground based facilities such as arcjets. By utilization of the same CFD code and solution methodology, the important differences between ground test and flight may be quantified. The goal of this paper is to utilize CFD to provide validated simulations of the flow environment in the NASA-Ames semi elliptic nozzle arcjet facilities. The validation of the ground simulations will come From comparison to existing calibration data. Specific tests in support of the X33 TPS test program will ilso be simulated. In this manner, the differences between the ground test simulation and the flight environment can be identified for a measure of ground test to flight traceability.

Loomis, Mark↗

Fatigue Life Methodology for Tapered Hybrid Composite Flexbeams

Nonlinear-tapered flexbeam specimens from a full-size composite helicopter rotor hub flexbeam were tested under combined constant axial tension and cyclic bending loads. Two different graphite/glass hybrid configurations tested under cyclic loading failed by delamination in the tapered region. A 2-D finite element model was developed which closely approximated the flexbeam geometry, boundary conditions, and loading. The analysis results from two geometrically nonlinear finite element codes, ANSYS and ABAQUS, are presented and compared. Strain energy release rates (G) associated with simulated delamination growth in the flexbeams are presented from both codes. These results compare well with each other and suggest that the initial delamination growth from the tip of the ply-drop toward the thick region of the flexbeam is strongly mode II. The peak calculated G values were used with material characterization data to calculate fatigue life curves for comparison with test data. A curve relating maximum surface strain to number of loading cycles at delamination onset compared well with the test results.

urri, Gretchen B.↗

Structural-Thermal-Optical-Performance (STOP) Model Development and Analysis of a Field-widened Michelson Interferometer

An integrated Structural-Thermal-Optical-Performance (STOP) model was developed for a field-widened Michelson interferometer which is being built and tested for the High Spectral Resolution Lidar (HSRL) project at NASA Langley Research Center (LaRC). The performance of the interferometer is highly sensitive to thermal expansion, changes in refractive index with temperature, temperature gradients, and deformation due to mounting stresses. Hand calculations can only predict system performance for uniform temperature changes, under the assumption that coefficient of thermal expansion (CTE) mismatch effects are negligible. An integrated STOP model was developed to investigate the effects of design modifications on the performance of the interferometer in detail, including CTE mismatch, and other three- dimensional effects. The model will be used to improve the design for a future spaceflight version of the interferometer. The STOP model was developed using the Comet SimApp'TM' Authoring Workspace which performs automated integration between Pro-Engineer®, Thermal Desktop®, MSC Nastran'TM', SigFit'TM', Code V'TM', and MATLAB®. This is the first flight project for which LaRC has utilized Comet, and it allows a larger trade space to be studied in a shorter time than would be possible in a traditional STOP analysis. This paper describes the development of the STOP model, presents a comparison of STOP results for simple cases with hand calculations, and presents results of the correlation effort to bench-top testing of the interferometer. A trade study conducted with the STOP model which demonstrates a few simple design changes that can improve the performance seen in the lab is also presented.

Scola, Salvatore J.↗

Comparisons of cross-section predictions for relativistic iron and argon beams with semiempirical fragmentation models

Cross-section predictions with semi-empirical nuclear fragmentation models from the Langley Research Center and the Naval Research Laboratory are compared with experimental data for the breakup of relativistic iron and argon projectile nuclei in various targets. Both these models are commonly used to provide fragmentation cross-section inputs into galactic cosmic ray transport codes for shielding and exposure analyses. Overall, the Langley model appears to yield better agreement with the experimental data.

Townsend, Lawrence W.↗

Dual-Code Solution Strategy for Chemically-Reacting Hypersonic Flows

A new procedure seeks to combine the thin-layer Navier-Stokes solver LAURA with the parabolized Navier-Stokes solver UPS for the aerothermodynamic solution of chemically-reacting air flow fields. The interface protocol is presented and the method is applied to two slender, blunted shapes. Both axisymmetric and three-dimensional solutions are included with surface pressure and heat transfer comparisons between the present method and previously published results. The case of Mach 25 flow over an axisymmetric six degree sphere-cone with a non-catalytic wall is considered to 100 nose radii. A stability bound on the marching step size was observed with this case and is attributed to chemistry effects resulting from the non-catalytic wall boundary condition. A second case with Mach 28 flow over a sphere-cone-cylinder-flare configuration is computed at both two and five degree angles of attack with a fully-catalytic wall. Surface pressures are seen to be within five percent with the present method compared to the baseline LAURA solution and heat transfers are within 10 percent. The effect of grid resolution is investigated in both the radial and streamwise directions. The procedure demonstrates significant, order of magnitude reductions in solution time and required memory for the three-dimensional case in comparison to an all thin-layer Navier-Stokes solution.

Wood, William A.↗

A Comparison of EAST Shock-Tube Radiation Measurements with a New Air Radiation Model

This paper presents a comparison between the recent EAST shock tube radiation measurements (Grinstead et al., AIAA 2008-1244) and the HARA radiation model. The equilibrium and nonequilibrium radiation measurements are studied for conditions relevant to lunar-return shock-layers; specifically shock velocities ranging from 9 to 11 kilometers per second at initial pressures of 0.1 and 0.3 Torr. The simulated shock-tube flow is assumed one-dimensional and is calculated using the LAURA code, while a detailed nonequilibrium radiation prediction is obtained in an uncoupled manner from the HARA code. The measured and predicted intensities are separated into several spectral ranges to isolate significant spectral features, mainly strong atomic line multiplets. The equations and physical data required for the prediction of these strong atomic lines are reviewed and their uncertainties identified. The 700-1020 nm wavelength range, which accounts for roughly 30% of the radiative flux to a peak-heating lunar return shock-layer, is studied in detail and the measurements and predictions are shown to agree within 15% in equilibrium. The plus or minus 1.5% uncertainty on the measured shock velocity is shown to cause up to a plus or minus 30% difference in the predicted radiation. This band of predictions contains the measured values in almost all cases. For the highly nonequilibrium 0.1 Torr cases, the nonequilibrium radiation peaks are under-predicted by about half. This under-prediction is considered acceptable when compared to the order-of-magnitude over-prediction obtained using a Boltzmann population of electronic states. The reasonable comparison in the nonequilibrium regions provides validation for both the non-Boltzmann modeling in HARA and the thermochemical nonequilibrium modeling in LAURA. The N2 (+)(1-) and N2(2+) molecular band systems are studied in the 290 480 nm wavelength range for both equilibrium and nonequilibrium regimes. The non-Boltzmann rate models for these systems, which have significant uncertainties, are tuned to improve the comparison with measurements.

Johnston, Christopher O.↗