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,135 records · Page 63

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations that include transverse-shear-deformation effects is presented. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Using several simplifying assumptions, linearized, stability equations are derived that describe the response of the plate just after bifurcation buckling occurs. These equations are then modified to allow the plate reference surface to be located a distance z(c), from the centroid surface which is convenient for modeling stiffened-plate assemblies. The implementation of the new theory into the VICONOPT buckling and vibration analysis and optimum design program code is described. Either classical plate theory (CPT) or first-order shear-deformation plate theory (SDPT) may be selected in VICONOPT. Comparisons of numerical results for several example problems with different loading states are made. Results from the new curved-plate analysis compare well with closed-form solution results and with results from known example problems in the literature. Finally, a design-optimization study of two different cylindrical shells subject to uniform axial compression is presented.

McGowan, David M.↗

A High Altitude-Low Reynolds Number Aerodynamic Flight Experiment

A sailplane is currently being developed at NASA's Dryden Flight Research Center to support a high altitude flight experiment. The purpose of the experiment is to measure the performance characteristics of an airfoil at altitudes between 100,000 and 70,000 feet at Mach numbers between 0.65 and 0.5. The airfoil lift and drag are measured from pilot and static pressures. The location of the separation bubble and vortex shedding are measured from a hot film strip. The details of the flight experiment are presented. A comparison of several estimates of the airfoil performance is also presented. The airfoil, APEX-16, was designed by Drela (MIT) with his MSES code. A two dimensional Navier-Stokes analysis has been performed by Tatineni and Zhong (UCLA) and another at the Dryden Flight Research Center. The role these analysis served to define the experiment is discussed.

Greer, Don↗

Relativistic Radiative and Auger Rates for Fe XXIV

As part of a project to compute improved atomic data for the spectral modeling of iron K lines, we report extensive calculations and comparisons of radiative and Auger rates for transitions involving the K-vacancy states in Fe XXIV. By making use of several computational codes, a detailed study is carried out of orbital representation, configuration interaction, relativistic corrections, cancellation effects, and fine tuning. It is shown that a formal treatment of the Breit interaction is essential to render the important magnetic correlations that take part in the decay pathways of this ion. As a result, the accuracy of the present A-values is firmly ranked at better than 10% while that of the Auger rates at only 15%.

Bautista, M. A.↗

Verification Assessment of Flow Boundary Conditions for CFD Analysis of Supersonic Inlet Flows

Boundary conditions for subsonic inflow, bleed, and subsonic outflow as implemented into the WIND CFD code are assessed with respect to verification for steady and unsteady flows associated with supersonic inlets. Verification procedures include grid convergence studies and comparisons to analytical data. The objective is to examine errors, limitations, capabilities, and behavior of the boundary conditions. Computational studies were performed on configurations derived from a "parameterized" supersonic inlet. These include steady supersonic flows with normal and oblique shocks, steady subsonic flow in a diffuser, and unsteady flow with the propagation and reflection of an acoustic disturbance.

Slater, John W.↗

Three-Dimensional High-Order Spectral Finite Volume Method for Unstructured Grids

Many areas require a very high-order accurate numerical solution of conservation laws for complex shapes. This paper deals with the extension to three dimensions of the Spectral Finite Volume (SV) method for unstructured grids, which was developed to solve such problems. We first summarize the limitations of traditional methods such as finite-difference, and finite-volume for both structured and unstructured grids. We then describe the basic formulation of the spectral finite volume method. What distinguishes the SV method from conventional high-order finite-volume methods for unstructured triangular or tetrahedral grids is the data reconstruction. Instead of using a large stencil of neighboring cells to perform a high-order reconstruction, the stencil is constructed by partitioning each grid cell, called a spectral volume (SV), into 'structured' sub-cells, called control volumes (CVs). One can show that if all the SV cells are partitioned into polygonal or polyhedral CV sub-cells in a geometrically similar manner, the reconstructions for all the SVs become universal, irrespective of their shapes, sizes, orientations, or locations. It follows that the reconstruction is reduced to a weighted sum of unknowns involving just a few simple adds and multiplies, and those weights are universal and can be pre-determined once for all. The method is thus very efficient, accurate, and yet geometrically flexible. The most critical part of the SV method is the partitioning of the SV into CVs. In this paper we present the partitioning of a tetrahedral SV into polyhedral CVs with one free parameter for polynomial reconstructions up to degree of precision five. (Note that the order of accuracy of the method is one order higher than the reconstruction degree of precision.) The free parameter will be determined by minimizing the Lebesgue constant of the reconstruction matrix or similar criteria to obtain optimized partitions. The details of an efficient, parallelizable code to solve three-dimensional problems for any order of accuracy are then presented. Important aspects of the data structure are discussed. Comparisons with the Discontinuous Galerkin (DG) method are made. Numerical examples for wave propagation problems are presented.

Liu, Yen↗

Micromechanics-Based Progressive Failure Analysis of Composite Laminates Using Different Constituent Failure Theories

Predicting failure in a composite can be done with ply level mechanisms and/or micro level mechanisms. This paper uses the Generalized Method of Cells and High-Fidelity Generalized Method of Cells micromechanics theories, coupled with classical lamination theory, as implemented within NASA's Micromechanics Analysis Code with Generalized Method of Cells. The code is able to implement different failure theories on the level of both the fiber and the matrix constituents within a laminate. A comparison is made among maximum stress, maximum strain, Tsai-Hill, and Tsai-Wu failure theories. To verify the failure theories the Worldwide Failure Exercise (WWFE) experiments have been used. The WWFE is a comprehensive study that covers a wide range of polymer matrix composite laminates. The numerical results indicate good correlation with the experimental results for most of the composite layups, but also point to the need for more accurate resin damage progression models.

Moncada, Albert M.↗

Longwave Band-by-band Cloud Radiative Effect and its Application in GCM Evaluation

The cloud radiative effect (CRE) of each longwave (LW) absorption band of a GCM fs radiation code is uniquely valuable for GCM evaluation because (1) comparing band-by-band CRE avoids the compensating biases in the broadband CRE comparison and (2) the fractional contribution of each band to the LW broadband CRE (f(sub CRE)) is sensitive to cloud top height but largely insensitive to cloud fraction, presenting thus a diagnostic metric to separate the two macroscopic properties of clouds. Recent studies led by the first author have established methods to derive such band ]by ]band quantities from collocated AIRS and CERES observations. We present here a study that compares the observed band-by-band CRE over the tropical oceans with those simulated by three different atmospheric GCMs (GFDL AM2, NASA GEOS-5, and CCCma CanAM4) forced by observed SST. The models agree with observation on the annual ]mean LW broadband CRE over the tropical oceans within +/-1W/sq m. However, the differences among these three GCMs in some bands can be as large as or even larger than +/-1W/sq m. Observed seasonal cycles of f(sub CRE) in major bands are shown to be consistent with the seasonal cycle of cloud top pressure for both the amplitude and the phase. However, while the three simulated seasonal cycles of f(sub CRE) agree with observations on the phase, the amplitudes are underestimated. Simulated interannual anomalies from GFDL AM2 and CCCma CanAM4 are in phase with observed anomalies. The spatial distribution of f(sub CRE) highlights the discrepancies between models and observation over the low-cloud regions and the compensating biases from different bands.

Huang, Xianglei↗

Supersonic Retropropulsion CFD Validation with Ames Unitary Plan Wind Tunnel Test Data

A validation study of Computational Fluid Dynamics (CFD) for Supersonic Retropropulsion (SRP) was conducted using three Navier-Stokes flow solvers (DPLR, FUN3D, and OVERFLOW). The study compared results from the CFD codes to each other and also to wind tunnel test data obtained in the NASA Ames Research Center 90 70 Unitary PlanWind Tunnel. Comparisons include surface pressure coefficient as well as unsteady plume effects, and cover a range of Mach numbers, levels of thrust, and angles of orientation. The comparisons show promising capability of CFD to simulate SRP, and best agreement with the tunnel data exists for the steadier cases of the 1-nozzle and high thrust 3-nozzle configurations.

Schauerhamer, Daniel G.↗

A Benchmark Example for Delamination Propagation Predictions Based on the Single Leg Bending Specimen Under Quasi-Static and Fatigue Loading

Benchmark examples based on Single Leg Bending (SLB) specimens with equal and unequal bending arm thicknesses were used to assess the performance of delamination prediction capabilities in finite element codes. First, the development of the quasi-static benchmark cases using the Virtual Crack Closure Technique (VCCT) is discussed in detail. Second, based on the quasi-static benchmark results, additional benchmark cases to assess delamination propagation under fatigue loading are created. Third, the application is demonstrated for the commercial finite element code Abaqus Standard 2018. The benchmark cases are compared to results obtained from VCCT-based, automated quasi-static propagation analysis. A comparison with results from automated fatigue propagation analysis was not performed at this point since the current version of Abaqus does not include this capability under variable mixed-mode conditions. In general, good agreement between the results obtained from the quasi-static propagation analysis and the benchmark results were achieved. Overall, the benchmarking procedure proved valuable for analysis verification.

Krueger, Ronald↗

A Benchmark Example for Delamination Propagation Predictions Based on the Single Leg Bending Specimen under Quasi-static and Fatigue Loading

Benchmark examples based on Single Leg Bending (SLB) specimens with equal and unequal bending arm thicknesses were used to assess the performance of delamination prediction capabilities in finite element codes. First, the development of the quasi-static benchmark cases using the Virtual Crack Closure Technique (VCCT) is discussed in detail. Second, based on the quasi-static benchmark results, additional benchmark cases to assess delamination propagation under fatigue loading are created. Third, the application is demonstrated for the commercial finite element code Abaqus Standard 2018. The benchmark cases are compared to results obtained from VCCT-based, automated quasi-static propagation analysis. A comparison with results from automated fatigue propagation analysis was not performed at this point since the current version of Abaqus does not include this capability under variable mixed-mode conditions. In general, good agreement between the results obtained from the quasi-static propagation analysis and the benchmark results were achieved. Overall, the benchmarking procedure proved valuable for analysis verification.

Ronald Krueger↗

A comparison of LOWTRAN-7 corrected Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data with ground spectral measurements

Atmospheric correction of imaging spectroscopy data is required for quantitative analysis. Different models were proposed for atmospheric correction of these data. LOWTRAN-7 is a low-resolution model and computer code for predicting atmospheric transmittance and background radiance from 0 to 50,00 cm(sup -1) which was developed by the Air Force Geophysics Laboratory. The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data used are radiometrically calibrated and include the 28 Sep. 1989 Providence Fan flight line segment 07, California. It includes a dark gravel surface defined as a calibration site by the Geologic Remote Sensing Field Experiment (GRSFE). Several ground measurements of portable spectrometer DAEDALUS AA440 Spectrafax were taken during the GRSFE, July 1989 field campaign. Comparisons of the LOWTRAN-7 corrected AVIRIS data with the ground spectrometer measurement were made.

Xu, Peng-Yang↗

Increased heat transfer to elliptical leading edges due to spanwise variations in the freestream momentum: Numerical and experimental results

A study of the effect of spanwise variation in momentum on leading edge heat transfer is discussed. Numerical and experimental results are presented for both a circular leading edge and a 3:1 elliptical leading edge. Reynolds numbers in the range of 10,000 to 240,000 based on leading edge diameter are investigated. The surface of the body is held at a constant uniform temperature. Numerical and experimental results with and without spanwise variations are presented. Direct comparison of the two-dimensional results, that is, with no spanwise variations, to the analytical results of Frossling is very good. The numerical calculation, which uses the PARC3D code, solves the three-dimensional Navier-Stokes equations, assuming steady laminar flow on the leading edge region. Experimentally, increases in the spanwise-averaged heat transfer coefficient as high as 50 percent above the two-dimensional value were observed. Numerically, the heat transfer coefficient was seen to increase by as much as 25 percent. In general, under the same flow conditions, the circular leading edge produced a higher heat transfer rate than the elliptical leading edge. As a percentage of the respective two-dimensional values, the circular and elliptical leading edges showed similar sensitivity to span wise variations in momentum. By equating the root mean square of the amplitude of the spanwise variation in momentum to the turbulence intensity, a qualitative comparison between the present work and turbulent results was possible. It is shown that increases in leading edge heat transfer due to spanwise variations in freestream momentum are comparable to those due to freestream turbulence.

Rigby, D. L.↗

Fast algorithms for combustion kinetics calculations: A comparison

To identify the fastest algorithm currently available for the numerical integration of chemical kinetic rate equations, several algorithms were examined. Findings to date are summarized. The algorithms examined include two general-purpose codes EPISODE and LSODE and three special-purpose (for chemical kinetic calculations) codes CHEMEQ, CRK1D, and GCKP84. In addition, an explicit Runge-Kutta-Merson differential equation solver (IMSL Routine DASCRU) is used to illustrate the problems associated with integrating chemical kinetic rate equations by a classical method. Algorithms were applied to two test problems drawn from combustion kinetics. These problems included all three combustion regimes: induction, heat release and equilibration. Variations of the temperature and species mole fraction are given with time for test problems 1 and 2, respectively. Both test problems were integrated over a time interval of 1 ms in order to obtain near-equilibration of all species and temperature. Of the codes examined in this study, only CREK1D and GCDP84 were written explicitly for integrating exothermic, non-isothermal combustion rate equations. These therefore have built-in procedures for calculating the temperature.

Radhakrishnan, K.↗

Lossless compression of AVIRIS data: Comparison of methods and instrument constraints

A family of lossless compression methods, allowing exact image reconstruction, are evaluated for compressing Airborne Visible/Infrared Imaging Spectrometers (AVIRIS) image data. The methods are used on Differential Pulse Code Modulation (DPCM). The compressed data have an entropy of order 6 bits/pixel. A theoretical model indicates that significantly better lossless compression is unlikely to be achieved because of limits caused by the noise in the AVIRIS channels. AVIRIS data differ from data produced by other visible/near-infrared sensors, such as LANDSAT-TM or SPOT, in several ways. Firstly, the data are recorded at a greater resolution (12 bits, though packed into 16-bit words). Secondly, the spectral channels are relatively narrow and provide continuous coverage of the spectrum so that the data in adjacent channels are generally highly correlated. Thirdly, the noise characteristics of the AVIRIS are defined by the channels' Noise Equivalent Radiances (NER's), and these NER's show that, at some wavelengths, the least significant 5 or 6 bits of data are essentially noise.

Roger, R. E.↗

Comparison of measured and predicted pure tone propagation levels from JAPE-1: An evaluation of the performance of ASOPRAT

Joint Acoustic Propagation Experiment Phase One (JAPE-1 ) short range propagation data has been used to evaluate the performance of the Advanced Sound Propagation in the Atmosphere (ASOPRAT) prediction code. The pure tone short range data was Fourier analyzed giving the propagated pressure levels as a function of frequency. Meteorological profiles measured at the experimental site were used as input for the acoustic prediction routine ASOPRAT. Predicted and measured propagation levels are compared in decibels (dB) relative to one of the measurement positions for receivers on the line passing between the two thirty meter towers. Agreement between predicted and measured levels is very good. Source strength data was not available, hence the comparisons show good agreement as to the shape of the propagation loss curve not necessarily the propagation levels.

Frederickson, Carl K.↗

Characterization of Aerodynamic Interactions with the Mars Science Laboratory Reaction Control System Using Computation and Experiment

On August 5 , 2012, the Mars Science Laboratory (MSL) entry capsule successfully entered Mars' atmosphere and landed the Curiosity rover in Gale Crater. The capsule used a reaction control system (RCS) consisting of four pairs of hydrazine thrusters to fly a guided entry. The RCS provided bank control to fly along a flight path commanded by an onboard computer and also damped unwanted rates due to atmospheric disturbances and any dynamic instabilities of the capsule. A preliminary assessment of the MSL's flight data from entry showed that the capsule flew much as predicted. This paper will describe how the MSL aerodynamics team used engineering analyses, computational codes and wind tunnel testing in concert to develop the RCS system and certify it for flight. Over the course of MSL's development, the RCS configuration underwent a number of design iterations to accommodate mechanical constraints, aeroheating concerns and excessive aero/RCS interactions. A brief overview of the MSL RCS configuration design evolution is provided. Then, a brief description is presented of how the computational predictions of RCS jet interactions were validated. The primary work to certify that the RCS interactions were acceptable for flight was centered on validating computational predictions at hypersonic speeds. A comparison of computational fluid dynamics (CFD) predictions to wind tunnel force and moment data gathered in the NASA Langley 31-Inch Mach 10 Tunnel was the lynch pin to validating the CFD codes used to predict aero/RCS interactions. Using the CFD predictions and experimental data, an interaction model was developed for Monte Carlo analyses using 6-degree-of-freedom trajectory simulation. The interaction model used in the flight simulation is presented.

Schoenenberger, Mark↗

Application of transport techniques to the analysis of NERVA shadow shields

A radiation shield internal to the NERVA nuclear rocket reactor required to limit the neutron and photon radiation levels at critical components located external to the reactor was evaluated. Two significantly different shield mockups were analyzed: BATH, a composite mixture of boron carbide, aluminum and titanium hydride, and a borated steel-liquid hydrogen system. Based on the comparisons between experimental and calculated neutron and photon radiation levels, the following conclusions were noted: (1) The ability of two-dimensional discrete ordinates code to predict the radiation levels internal to and at the surface of the shield mockups was clearly demonstrated. (2) Internal to the BATH shield mockups, the one-dimensional technique predicted the axial variation of neutron fluxes and photon dose rates; however, the magnitude of the neutron fluxes was about a factor of 1.8 lower than the two-dimensional analysis and the photon dose rate was a factor of 1.3 lower.

Capo, M. A.↗

Component aerodynamic characteristics of banked cruciform missiles with arbitrary control deflections

A comprehensive computer code has been developed for the prediction of the static aerodynamic characteristics of cruciform wing-body combinations with one or two sets of fins. The code may be used for subsonic, transonic, and supersonic speeds. Arbitrary bank angles and control deflections are allowed. The present version of the code may be used for angles of attack up to 20 deg. The code computes the forces and moments acting on each component of the missile including the effects of vortices shed by the body and fins. Comparisons between experiment and theory indicate good agreement for longitudinal characteristics and fair to good agreement for lateral characteristics.

Hemsch, M. J.↗