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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Experimental and numerical study of swept ramp injection into a supersonic flowfield

Time-averaged measurements of pressure, temperature, velocity, and injectant mole fraction are presented using the planar laser-induced iodine fluorescence technique in the complex three-dimensional compressible flowfield around a swept ramp fuel injector. Within the range of thermodynamic conditions present in the test case studied, the technique's accuracy is estimated to be 4% for pressure, temperature, and velocity and 3% for injectant mole fraction. Comparisons with numerical simulations using the SPARK three-dimensional Navier-Stokes computer code with an algebraic turbulence model are made at the centerplane of the flowfield as well as on three crossflow planes downstream of the injector. Calculations and measurements are in good agreement throughout the flowfield, with deviations on the order of 5%; however, in specific regions, such as in the base of the ramp, deviations are larger. A weak asymmetry in the incoming flowfield appears to be amplified by boundary-layer separation occurring when the ramp-generated shock reflects off the tunnel walls. Ramp-generated vortices are weaker in the calculated results due to the effects of numerical viscosity in the vortex cores. This leads to less turning and mixing of the jet plume than observed in the experiments. The rate of decay of the maximum injectant mole fraction with streamwise distance is greater for the present ramp injection scheme than for previously measured transverse injection schemes. In this recirculation region at the base of the injector, laminar calculations show better agreement with the measurements than turbulent calculations.

Donohue, James M.↗
NASA NTRS1999

Realistic NLTE Radiative Transfer for Modeling Stellar Winds

This NASA grant supported the development of codes to solve the non-LTE multi-level spherical radiative transfer problem in the presence of velocity fields. Much of this work was done in collaboration with Graham Harper (CASA, University of Colorado). These codes were developed for application to the cool, low-velocity winds of evolved late-type stars. Particular emphasis was placed on modeling the wind of lambda Velorum (K4 lb), the brightest K supergiant in the sky, based on extensive observations of the ultraviolet spectrum with the HST/GHRS from GO program 5307. Several solution techniques were examined, including the Eddington factor Approach described in detail by Bennett & Harper (1997). An Eddington factor variant of Harper's S-MULTI code (Harper 1994) for stationary atmospheres was developed and implemented, although full convergence was not realized. The ratio of wind terminal velocity to turbulent velocity is large (approx. 0.3-0.5) in these cool star winds so this assumption of stationarity provides reasonable starting models. Final models, incorporating specified wind laws, were converged using the comoving CRD S-MULTI code. Details of the solution procedure were published by Bennett & Harper (1997). Our analysis of the wind of lambda Vel, based on wind absorption superimposed on chromospheric emission lines in the ultraviolet, can be found in Carpenter et al. (1999). In this paper, we compare observed wind absorption features to an exact CRD calculation in the comoving frame, and also to a much quicker, but approximate, method using the SEI (Sobolev with Exact Integration) code of Lamers, Cerruti-Sola, & Perinotto (1987). Carpenter et al. (1999) provide detailed comparisons of the exact CRD and approximate SEI results and discuss when SEI is adequate to use for computing wind line profiles. Unfortunately, the observational material is insufficient to unambiguously determine the wind acceleration law for lambda Vel. Relatively few unblended Fe II lines of optical depth sensitive to the wind acceleration region are present in the existing HST/GHRS data set. Most of the Fe II lines are either too optically thick (resulting in a board, black wind absorption profile) or too optically thin (resulting in no wind absorption feature present). Also, most of the ultraviolet spectra obtained from HST GO-5307 was at medium resolution (R approx. 40,000, corresponding to a velocity resolution of 7.5 km/s). This is simply inadequate to resolve the turbulence in the outer wind; a key parameter in theoretical wind models. We can now say that an unambiguous determination of the wind velocity law in lambda Vel will require complete coverage of the ultraviolet spectrum at high dispersion (R approx. 10(exp 5), or 3 km/s). This is now feasible usin, the STIS echelle spectrograph on-board HST.

Bennett, Philip D.↗
NASA NTRS2003

Improvements to Wall Corrections at the NASA Langley 14 x 22-Ft Subsonic Tunnel

The new wall pressure measurement system and the TWICS wall correction system for the 14x22-Ft subsonic tunnel are described. Results from a recent semispan test and a full-span test are presented. Comparison with existing classical methods of correction is shown. A modification of the TWICS code to treat the effect due to a deflected wake from a high-lift wing is also discussed. The current implementation of TWICS for the 14x22-Ft tunnel is shown to be an improvement over existing methods.

Iyer, Venkit↗
NASA NTRS2004

Discrete Spring Model for Predicting Delamination Growth in Z-Fiber Reinforced DCB Specimens

Beam theory analysis was applied to predict delamination growth in DCB specimens reinforced in the thickness direction with pultruded pins, known as Z-fibers. The specimen arms were modeled as cantilever beams supported by discrete springs, which were included to represent the pins. A bi-linear, irreversible damage law was used to represent Z-fiber damage, the parameters of which were obtained from previous experiments. Closed-form solutions were developed for specimen compliance and displacements corresponding to Z-fiber row locations. A solution strategy was formulated to predict delamination growth, in which the parent laminate mode I fracture toughness was used as the criterion for delamination growth. The solution procedure was coded into FORTRAN 90, giving a dedicated software tool for performing the delamination prediction. Comparison of analysis results with previous analysis and experiment showed good agreement, yielding an initial verification for the analytical procedure.

Ratcliffe, James G.↗
NASA NTRS1992

Dose rate, dose-equivalent rate, and quality factor in SLS-1

A tissue-equivalent proportional counter (TEPC) sensitive to the lineal energy range of 0.26-300 keV micrometer-1 was flown on STS-40 (39 degrees x 278 km x 296 km) inside the Spacelab. This instrument was previously flown on STS-31 but was modified to provide a finer resolution at lower lineal energies to better map the South Atlantic Anomaly (SAA) protons. The instrument was turned on 6 June 1991, and operated for 7470 min (124.5 h). The flight duration was characterized by a very large number of X-ray solar flares and enhanced magnetic field fluctuations; however, no significant dose from the solar particles was measured at the location of this instrument. The flight data can be separated into trapped and galactic cosmic radiation parts. The dose rate, dose-equivalent rate and quality factor for trapped radiation were 4.21 +/- 0.03 mrad day-1, 7.72 +/- 0.05 mrem day-1, and 1.83 +/- 0.1, respectively. The dose rate, dose-equivalent rate, and quality factor for galactic cosmic radiation were 5.34 +/- 0.03 mrad day-1, 14.63 +/- 0.06 mrem day-1, and 2.74 +/- 0.1, respectively. The overall quality factor for the flight was 2.38. The dose from the GCR is higher than from SAA protons because of the high inclination and low altitude of this flight. The AP8MAX model of the trapped radiation gives a dose rate of 2.43 mrad day-1 and a quality factor of 1.77. The CREME solar maximum model of galactic cosmic radiation gives a dose rate of 2.54 mrad day-1 and a quality factor of 2.91. Thus the AP8MAX model underestimates the dose by a factor of 1.8 whereas the CREME model leads to an underestimation of the dose by a factor of 2. A comparison of the LET spectra using the AP8MAX model and galactic cosmic radiation transport codes shows only a qualitative agreement.

NASA Discipline Number 00-00↗
NASA NTRS2004

Discrete Spring Model for Predicting Delamination Growth in Z-Fiber Reinforced DCB Specimens

Beam theory analysis was applied to predict delamination growth in Double Cantilever Beam (DCB) specimens reinforced in the thickness direction with pultruded pins, known as Z-fibers. The specimen arms were modeled as cantilever beams supported by discrete springs, which were included to represent the pins. A bi-linear, irreversible damage law was used to represent Z-fiber damage, the parameters of which were obtained from previous experiments. Closed-form solutions were developed for specimen compliance and displacements corresponding to Z-fiber row locations. A solution strategy was formulated to predict delamination growth, in which the parent laminate mode I critical strain energy release rate was used as the criterion for delamination growth. The solution procedure was coded into FORTRAN 90, giving a dedicated software tool for performing the delamination prediction. Comparison of analysis results with previous analysis and experiment showed good agreement, yielding an initial verification for the analytical procedure.

Ratcliffe, James G.↗
NASA NTRS2010

Predictions of Supersonic Jet Mixing and Shock-Associated Noise Compared With Measured Far-Field Data

Codes for predicting supersonic jet mixing and broadband shock-associated noise were assessed using a database containing noise measurements of a jet issuing from a convergent nozzle. Two types of codes were used to make predictions. Fast running codes containing empirical models were used to compute both the mixing noise component and the shock-associated noise component of the jet noise spectrum. One Reynolds-averaged, Navier-Stokes-based code was used to compute only the shock-associated noise. To enable the comparisons of the predicted component spectra with data, the measured total jet noise spectra were separated into mixing noise and shock-associated noise components. Comparisons were made for 1/3-octave spectra and some power spectral densities using data from jets operating at 24 conditions covering essentially 6 fully expanded Mach numbers with 4 total temperature ratios.

Dahl, Milo D.↗
NASA NTRS2012

Unsteady Aero Computation of a 1 1/2 Stage Large Scale Rotating Turbine

This report is the documentation of the work performed for the Subsonic Rotary Wing Project under the NASA s Fundamental Aeronautics Program. It was funded through Task Number NNC10E420T under GESS-2 Contract NNC06BA07B in the period of 10/1/2010 to 8/31/2011. The objective of the task is to provide support for the development of variable speed power turbine technology through application of computational fluid dynamics analyses. This includes work elements in mesh generation, multistage URANS simulations, and post-processing of the simulation results for comparison with the experimental data. The unsteady CFD calculations were performed with the TURBO code running in multistage single passage (phase lag) mode. Meshes for the blade rows were generated with the NASA developed TCGRID code. The CFD performance is assessed and improvements are recommended for future research in this area. For that, the United Technologies Research Center's 1 1/2 stage Large Scale Rotating Turbine was selected to be the candidate engine configuration for this computational effort because of the completeness and availability of the data.

To, Wai-Ming↗
NASA NTRS2013

High-Contrast Near-Infrared Imaging Polarimetry of the Protoplanetary Disk around RY Tau

We present near-infrared coronagraphic imaging polarimetry of RY Tau. The scattered light in the circumstellar environment was imaged at H-band at a high resolution (approx. 0.05) for the first time, using Subaru-HiCIAO. The observed polarized intensity (PI) distribution shows a butterfly-like distribution of bright emission with an angular scale similar to the disk observed at millimeter wavelengths. This distribution is offset toward the blueshifted jet, indicating the presence of a geometrically thick disk or a remnant envelope, and therefore the earliest stage of the Class II evolutionary phase. We perform comparisons between the observed PI distribution and disk models with: (1) full radiative transfer code, using the spectral energy distribution (SED) to constrain the disk parameters; and (2) monochromatic simulations of scattered light which explore a wide range of parameters space to constrain the disk and dust parameters. We show that these models cannot consistently explain the observed PI distribution, SED, and the viewing angle inferred by millimeter interferometry. We suggest that the scattered light in the near-infrared is associated with an optically thin and geometrically thick layer above the disk surface, with the surface responsible for the infrared SED. Half of the scattered light and thermal radiation in this layer illuminates the disk surface, and this process may significantly affect the thermal structure of the disk.

protoplanetary disks↗
NASA NTRS2015

Impact Testing and Simulation of a Sinusoid Foam Sandwich Energy Absorber

A sinusoidal-shaped foam sandwich energy absorber was developed and evaluated at NASA Langley Research Center through multi-level testing and simulation performed under the Transport Rotorcraft Airframe Crash Testbed (TRACT) research project. The energy absorber, designated the "sinusoid," consisted of hybrid carbon- Kevlar® plain weave fabric face sheets, two layers for each face sheet oriented at +/-45deg with respect to the vertical or crush direction, and a closed-cell ELFOAM(TradeMark) P200 polyisocyanurate (2.0-lb/ft3) foam core. The design goal for the energy absorber was to achieve an average floor-level acceleration of between 25- and 40-g during the full-scale crash test of a retrofitted CH-46E helicopter airframe, designated TRACT 2. Variations in the design were assessed through quasi-static and dynamic crush testing of component specimens. Once the design was finalized, a 5-ft-long subfloor beam was fabricated and retrofitted into a barrel section of a CH-46E helicopter. A vertical drop test of the barrel section was conducted onto concrete to evaluate the performance of the energy absorber prior to retrofit into TRACT 2. Finite element models were developed of all test articles and simulations were performed using LSDYNA ®, a commercial nonlinear explicit transient dynamic finite element code. Test analysis results are presented for the sinusoid foam sandwich energy absorber as comparisons of load-displacement and acceleration-time-history responses, as well as predicted and experimental structural deformations and progressive damage for each evaluation level (component testing through barrel section drop testing).

Jackson, Karen E.↗
NASA NTRS1990

Parallel processors and nonlinear structural dynamics algorithms and software

Techniques are discussed for the implementation and improvement of vectorization and concurrency in nonlinear explicit structural finite element codes. In explicit integration methods, the computation of the element internal force vector consumes the bulk of the computer time. The program can be efficiently vectorized by subdividing the elements into blocks and executing all computations in vector mode. The structuring of elements into blocks also provides a convenient way to implement concurrency by creating tasks which can be assigned to available processors for evaluation. The techniques were implemented in a 3-D nonlinear program with one-point quadrature shell elements. Concurrency and vectorization were first implemented in a single time step version of the program. Techniques were developed to minimize processor idle time and to select the optimal vector length. A comparison of run times between the program executed in scalar, serial mode and the fully vectorized code executed concurrently using eight processors shows speed-ups of over 25. Conjugate gradient methods for solving nonlinear algebraic equations are also readily adapted to a parallel environment. A new technique for improving convergence properties of conjugate gradients in nonlinear problems is developed in conjunction with other techniques such as diagonal scaling. A significant reduction in the number of iterations required for convergence is shown for a statically loaded rigid bar suspended by three equally spaced springs.

Belytschko, Ted↗
NASA NTRS1994

Modeling of stress/strain behavior of fiber-reinforced ceramic matrix composites including stress redistribution

A computational simulation procedure is presented for nonlinear analyses which incorporates microstress redistribution due to progressive fracture in ceramic matrix composites. This procedure facilitates an accurate simulation of the stress-strain behavior of ceramic matrix composites up to failure. The nonlinearity in the material behavior is accounted for at the constituent (fiber/matrix/interphase) level. This computational procedure is a part of recent upgrades to CEMCAN (Ceramic Matrix Composite Analyzer) computer code. The fiber substructuring technique in CEMCAN is used to monitor the damage initiation and progression as the load increases. The room-temperature tensile stress-strain curves for SiC fiber reinforced reaction-bonded silicon nitride (RBSN) matrix unidirectional and angle-ply laminates are simulated and compared with experimentally observed stress-strain behavior. Comparison between the predicted stress/strain behavior and experimental stress/strain curves is good. Collectively the results demonstrate that CEMCAN computer code provides the user with an effective computational tool to simulate the behavior of ceramic matrix composites.

Mital, Subodh K.↗
NASA NTRS1994

An assessment of the adaptive unstructured tetrahedral grid, Euler Flow Solver Code FELISA

A three-dimensional solution-adaptive Euler flow solver for unstructured tetrahedral meshes is assessed, and the accuracy and efficiency of the method for predicting sonic boom pressure signatures about simple generic models are demonstrated. Comparison of computational and wind tunnel data and enhancement of numerical solutions by means of grid adaptivity are discussed. The mesh generation is based on the advancing front technique. The FELISA code consists of two solvers, the Taylor-Galerkin and the Runge-Kutta-Galerkin schemes, both of which are spacially discretized by the usual Galerkin weighted residual finite-element methods but with different explicit time-marching schemes to steady state. The solution-adaptive grid procedure is based on either remeshing or mesh refinement techniques. An alternative geometry adaptive procedure is also incorporated.

Djomehri, M. Jahed↗
NASA NTRS1999

Spaceflight Validation of Hzetrn Code

HZETRN is being developed as a fast deterministic radiation transport code applicable to neutrons, protons, and multiply charged ions in the space environment. It was recently applied to 50 hours of IMP8 data measured during the August 4, 1972 solar event to map the hourly exposures within the human body under several shield configurations. This calculation required only 18 hours on a VAX 4000 machine. A similar calculation using the Monte Carlo method would have required two years of dedicated computer time. The code has been benchmarked against well documented and tested Monte Carlo proton transport codes with good success. The code will allow important trade studies to be made with relative ease due to the computational speed and will be useful in assessing design alternatives in an integrated system software environment. Since there are no well tested Monte Carlo codes for HZE particles, we have been engaged in flight validation of the HZETRN results. To date we have made comparison with TEPC, CR-39, charge particle telescopes, and Bonner spheres. This broad range of detectors allows us to test a number of functions related to differing physical processes which add to the complicated radiation fields within a spacecraft or the human body, which functions can be calculated by the HZETRN code system. In the present report we will review these results.

Wilson, J. W.↗
NASA NTRS2010

Air-Loads Prediction of a UH-60A Rotor inside the 40- by 80-Foot Wind Tunnel

The presented research extends the capability of a loose coupling computational fluid dynamics (CFD) and computational structure dynamics (CSD) code to calculate the flow-field around a rotor and test stand mounted inside a wind tunnel. Comparison of predicted air-load results for a full-scale UH-60A rotor recently tested inside the National Full-Scale Aerodynamics Complex (NFAC) 40- by 80-Foot Wind Tunnel at Ames Research Center and in free-air flight are made for three challenging flight data points from the earlier conducted UH-60A Air-loads Program. Overall results show that the extension of the coupled CFD/CSD code to the wind-tunnel environment is generally successful.

Chang, I-Chung↗
NASA NTRS2005

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The paper discusses the following: 1. No-cost Extension. The no-cost extension is required to complete the work on the unified model codes (both hydrodynamic and kinetic Monte Carlo) as described in the initial proposal and previous annual reports. 2. Scientific Accomplishments during the Report Period. We completed a comprehensive model of Alfvtn wave reflection that spans the full distance from the photosphere to the distant heliosphere. 3. Comparison of Accomplishments with Proposed Goals. The proposal contained two specific objectives for Year 3: (1) to complete the unified model code, and (2) to apply it to various kinds of coronal holes (and polar plumes within coronal holes). Although the anticipated route toward these two final goals has changed (see accomplishments 2a and 2b above), they remain the major milestones for the extended period of performance. Accomplishments la and IC were necessary prerequisites for the derivation of "physically relevant transport and mode-coupling terms" for the unified model codes (as stated in the proposal Year 3 goals). We have fulfilled the proposed "core work" to study 4 general types of physical processes; in previous years we studied turbulence, mode coupling (Le., non-WKB reflection), and kinetic wave damping, and accomplishment lb provides the fourth topic: nonlinear steepening.

Wagner, William↗
NASA NTRS2003

Some Advanced Concepts in Discrete Aerodynamic Sensitivity Analysis

An efficient incremental iterative approach for differentiating advanced flow codes is successfully demonstrated on a two-dimensional inviscid model problem. The method employs the reverse-mode capability of the automatic differentiation software tool ADIFOR 3.0 and is proven to yield accurate first-order aerodynamic sensitivity derivatives. A substantial reduction in CPU time and computer memory is demonstrated in comparison with results from a straightforward, black-box reverse-mode applicaiton of ADIFOR 3.0 to the same flow code. An ADIFOR-assisted procedure for accurate second-rder aerodynamic sensitivity derivatives is successfully verified on an inviscid transonic lifting airfoil example problem. The method requires that first-order derivatives are calculated first using both the forward (direct) and reverse (adjoinct) procedures; then, a very efficient noniterative calculation of all second-order derivatives can be accomplished. Accurate second derivatives (i.e., the complete Hesian matrices) of lift, wave drag, and pitching-moment coefficients are calculated with respect to geometric shape, angle of attack, and freestream Mach number.

Taylor, Arthur C., III↗
NASA NTRS1988

Comparison of Space Shuttle Hot Gas Manifold analysis to air flow data

This paper summarizes several recent analyses of the Space Shuttle Main Engine Hot Gas Manifold and compares predicted flow environments to air flow data. Codes used in these analyses include INS3D, PAGE, PHOENICS, and VAST. Both laminar (Re = 250, M = 0.30) and turbulent (Re = 1.9 million, M = 0.30) results are discussed, with the latter being compared to data for system losses, outer wall static pressures, and manifold exit Mach number profiles. Comparison of predicted results for the turbulent case to air flow data shows that the analysis using INS3D predicted system losses within 1 percent error, while the PHOENICS, PAGE, and VAST codes erred by 31, 35, and 47 percent, respectively. The INS3D, PHOENICS, and PAGE codes did a reasonable job of predicting outer wall static pressure, while the PHOENICS code predicted exit Mach number profiles with acceptable accuracy. INS3D was approximately an order of magnitude more efficient than the other codes in terms of code speed and memory requirements. In general, it is seen that complex internal flows in manifold-like geometries can be predicted with a limited degree of confidence, and further development is necessary to improve both efficiency and accuracy of codes if they are to be used as design tools for complex three-dimensional geometries.

Mcconnaughey, P. K.↗