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 523 records · Page 29

Three dimensional full potential method for the aeroelastic modeling of propfans

A three-dimensional time dependent full potential aerodynamic analysis of single rotation propellers has been developed. The primary purpose of the code is to provide a capability of doing propfan aeroelastic analysis in the nonlinear transonic regime. A secondary purpose is to provide a validation of the unsteady lifting surface panel method that has been have developed. Results will be shown for steady state aerodynamic loading, unsteady aerodynamic response to forced aeroelastic deformations, and free aeroelastic response. Comparisons are made to experimental data and corresponding panel code results.

Williams, Marc H.↗

Analysis of STS-2 experimental heating rates and transition data

Experimental laminar and 'fully' turbulent windward-ray heating rates obtained from the second Space Shuttle flight are compared with predicted rates based on equilibrium-air approximate and detailed analyses. A comparison of the results of the prediction techniques yield discrepancies of approximately 10 percent. The experimental laminar heating rates for altitudes greater than 67 km are as much as 35 percent lower than the results of the approximate code over the first 40 percent of the shuttle length. The approximate equilibrium predictions are in good agreement with these laminar data beyond the 40-percent station and over the entire shuttle at altitudes lower than 67 km. However, recent results of a detailed viscous-shock-layer nonequilibrium code indicate significant departures, especially at high altitudes, from the equilibrium state. Since the comparison of these predicted values and the data result in some currently unresolved questions, an engineering correlation of the predictions for use in the approximate code is not feasible at this time. Factors which may affect the comparisons of data and the nonequilibrium predictions are considered. The turbulent heating comparisons are good and boundary-layer transition data are compared with a current boundary-layer transition criterion.

Zoby, E. V.↗

A Non-Cut Cell Immersed Boundary Method for Use in Icing Simulations

This paper describes a computational fluid dynamic method used for modelling changes in aircraft geometry due to icing. While an aircraft undergoes icing, the accumulated ice results in a geometric alteration of the aerodynamic surfaces. In computational simulations for icing, it is necessary that the corresponding geometric change is taken into consideration. The method used, herein, for the representation of the geometric change due to icing is a non-cut cell Immersed Boundary Method (IBM). Computational cells that are in a body fitted grid of a clean aerodynamic geometry that are inside a predicted ice formation are identified. An IBM is then used to change these cells from being active computational cells to having properties of viscous solid bodies. This method has been implemented in the NASA developed node centered, finite volume computational fluid dynamics code, FUN3D. The presented capability is tested for two-dimensional airfoils including a clean airfoil, an iced airfoil, and an airfoil in harmonic pitching motion about its quarter chord. For these simulations velocity contours, pressure distributions, coefficients of lift, coefficients of drag, and coefficients of pitching moment about the airfoil's quarter chord are computed and used for comparison against experimental results, a higher order panel method code with viscous effects, XFOIL, and the results from FUN3D's original solution process. The results of the IBM simulations show that the accuracy of the IBM compares satisfactorily with the experimental results, XFOIL results, and the results from FUN3D's original solution process.

Sarofeen, Christian M.↗

Method for the prediction of the installation aerodynamics of a propfan at subsonic speeds: User manual

This document is the user's guide for the method developed earlier for predicting the slipstream wing interaction at subsonic speeds. The analysis involves a subsonic panel code (HESS code) modified to handle the propeller onset flow. The propfan slipstream effects are superimposed on the normal flow boundary condition and are applied over the surface washed by the slipstream. The effects of the propeller slipstream are to increase the axial induced velocity, tangential velocity, and a total pressure rise in the wake of the propeller. Principles based on blade performance theory, momentum theory, and vortex theory were used to evaluate the slipstream effects. The code can be applied to any arbitrary three dimensional geometry, expressed in the form of HESS input format. The code can handle a propeller alone configuration or a propeller/nacelle/airframe configuration, operating up to high subcritical Mach numbers over a range of angles of attack. Inclusion of a viscous modelling is briefly outlined. Wind tunnel results/theory comparisons are included as examples for the application of the code to a generic supercritical wing/overwing Nacelle with a powered propfan. A sample input/output listing is provided.

Chandrasekaran, B.↗

Improved calculations for the C III 1907,1909 and Si III 1883,1892 electron density sensitive emission-line ratios, and a comparison with IUE observations

Atomic data are used in conjunction with the statistical equilibrium code of Dufton (1977) to calculate relative C III and Si III level populations, and hence emission-line strengths for a range of electron temperatures and densities. It is assumed that photoexcitation and deexcitaton rates are negligible in comparison with the corresponding collisional rates, that ionization to and recombination from other ionic levels are slow compared with bound-bound rates, and that all transitions are optically thin. The observed values of R1 and R2 for several planetary nebulae and a symbiotic star, measured from high-resolution spectra obtained with the IUE satellite, lead to electron densities that are compatible, and are also in good agreement with those deduced from line ratios in other species.

Keenan, F. P.↗

Computation of Thermally Perfect Oblique Shock Wave Properties

A set of compressible flow relations describing flow properties across oblique shock waves, derived for a thermally perfect, calorically imperfect gas, is applied within the existing thermally perfect gas (TPG) computer code. The relations are based upon the specific heat expressed as a polynomial function of temperature. The updated code produces tables of compressible flow properties of oblique shock waves, as well as the original properties of normal shock waves and basic isentropic flow, in a format similar to the tables for normal shock waves found in NACA Rep. 1135. The code results are validated in both the calorically perfect and the calorically imperfect, thermally perfect temperature regimes through comparisons with the theoretical methods of NACA Rep. 1135. The advantages of the TPG code for oblique shock wave calculations, as well as for the properties of isentropic flow and normal shock waves, are its ease of use and its applicability to any type of gas (monatomic, diatomic, triatomic, polyatomic, or any specified mixture thereof).

Tatum, Kenneth E.↗

Coupled Solid Rocket Motor Ballistics and Trajectory Modeling for Higher Fidelity Launch Vehicle Design

Multi-stage launch vehicles with solid rocket motors (SRMs) face design optimization challenges, especially when the mission scope changes frequently. Significant performance benefits can be realized if the solid rocket motors are optimized to the changing requirements. While SRMs represent a fixed performance at launch, rapid design iterations enable flexibility at design time, yielding significant performance gains. The streamlining and integration of SRM design and analysis can be achieved with improved analysis tools. While powerful and versatile, the Solid Performance Program (SPP) is not conducive to rapid design iteration. Performing a design iteration with SPP and a trajectory solver is a labor intensive process. To enable a better workflow, SPP, the Program to Optimize Simulated Trajectories (POST), and the interfaces between them have been improved and automated, and a graphical user interface (GUI) has been developed. The GUI enables real-time visual feedback of grain and nozzle design inputs, enforces parameter dependencies, removes redundancies, and simplifies manipulation of SPP and POST's numerous options. Automating the analysis also simplifies batch analyses and trade studies. Finally, the GUI provides post-processing, visualization, and comparison of results. Wrapping legacy high-fidelity analysis codes with modern software provides the improved interface necessary to enable rapid coupled SRM ballistics and vehicle trajectory analysis. Low cost trade studies demonstrate the sensitivities of flight performance metrics to propulsion characteristics. Incorporating high fidelity analysis from SPP into vehicle design reduces performance margins and improves reliability. By flying an SRM designed with the same assumptions as the rest of the vehicle, accurate comparisons can be made between competing architectures. In summary, this flexible workflow is a critical component to designing a versatile launch vehicle model that can accommodate a volatile mission scope.

Ables, Brett↗

Rapid Aeroelastic Analysis of Blade Flutter in Turbomachines

The LINFLUX-AE computer code predicts flutter and forced responses of blades and vanes in turbomachines under subsonic, transonic, and supersonic flow conditions. The code solves the Euler equations of unsteady flow in a blade passage under the assumption that the blades vibrate harmonically at small amplitudes. The steady-state nonlinear Euler equations are solved by a separate program, then equations for unsteady flow components are obtained through linearization around the steady-state solution. A structural-dynamics analysis (see figure) is performed to determine the frequencies and mode shapes of blade vibrations, a preprocessor interpolates mode shapes from the structural-dynamics mesh onto the LINFLUX computational-fluid-dynamics mesh, and an interface code is used to convert the steady-state flow solution to a form required by LINFLUX. Then LINFLUX solves the linearized equations in the frequency domain to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. A post-processor uses the unsteady pressures to calculate generalized aerodynamic forces, response amplitudes, and eigenvalues (which determine the flutter frequency and damping). In comparison with the TURBO-AE aeroelastic-analysis code, which solves the equations in the time domain, LINFLUX-AE is 6 to 7 times faster.

Trudell, J. J.↗

Improved Compression of Wavelet-Transformed Images

A recently developed data-compression method is an adaptive technique for coding quantized wavelet-transformed data, nominally as part of a complete image-data compressor. Unlike some other approaches, this method admits a simple implementation and does not rely on the use of large code tables. A common data compression approach, particularly for images, is to perform a wavelet transform on the input data, and then losslessly compress a quantized version of the wavelet-transformed data. Under this compression approach, it is common for the quantized data to include long sequences, or runs, of zeros. The new coding method uses prefixfree codes for the nonnegative integers as part of an adaptive algorithm for compressing the quantized wavelet-transformed data by run-length coding. In the form of run-length coding used here, the data sequence to be encoded is parsed into strings consisting of some number (possibly 0) of zeros, followed by a nonzero value. The nonzero value and the length of the run of zeros are encoded. For a data stream that contains a sufficiently high frequency of zeros, this method is known to be more effective than using a single variable length code to encode each symbol. The specific prefix-free codes used are from two classes of variable-length codes: a class known as Golomb codes, and a class known as exponential-Golomb codes. The codes within each class are indexed by a single integer parameter. The present method uses exponential-Golomb codes for the lengths of the runs of zeros, and Golomb codes for the nonzero values. The code parameters within each code class are determined adaptively on the fly as compression proceeds, on the basis of statistics from previously encoded values. In particular, a simple adaptive method has been devised to select the parameter identifying the particular exponential-Golomb code to use. The method tracks the average number of bits used to encode recent runlengths, and takes the difference between this average length and the code parameter. When this difference falls outside a fixed range, the code parameter is updated (increased or decreased). The Golomb code parameter is selected based on the average magnitude of recently encoded nonzero samples. The coding method requires no floating- point operations, and more readily adapts to local statistics than other methods. The method can also accommodate arbitrarily large input values and arbitrarily long runs of zeros. In practice, this means that changes in the dynamic range or size of the input data set would not require a change to the compressor. The algorithm has been tested in computational experiments on test images. A comparison with a previously developed algorithm that uses large code tables (generated via Huffman coding on training data) suggests that the data-compression effectiveness of the present algorithm is comparable to the best performance achievable by the previously developed algorithm.

Kiely, Aaron↗

Some experience with arc-heater simulation of outer planet entry radiation

An electric arc heater was operated at 800 amperes and 100,000 pa (1 atm) with hydrogen, helium, and two mixtures of hydrogen and helium. A VUV-scanning monochromator was used to record the spectra from an end view while a second spectrometer was used to determine the plasma temperature using hydrogen continuum radiation at 562 nm. Except for pure helium, the plasma temperature was found to be too low to produce significant helium radiation, and the measured spectra were primarily the hydrogen spectra with the highest intensity in the pure hydrogen case. A radiation computer code was used to compute the spectra for comparison to the measurements and to extend the study to simulation of outer planet entry radiation. Conductive cooling prevented ablation of phenolic carbon material samples mounted inside the arc heater during a cursory attempt to produce radiation absorption by ablation gases.

Wells, W. L.↗

A benchmark for galactic cosmic ray transport codes

A nontrivial analytic benchmark solution for galactic cosmic ray transport is presented for use in transport code validation. Computational accuracy for a previously-developed cosmic ray transport code is established to within one percent by comparison with this exact benchmark. Hence, solution accuracy for the transport problem is mainly limited by inaccuracies in the input spectra, input interaction databases, and the use of a straight ahead/velocity-conserving approximation.

Wilson, John W.↗

Recent advances in transonic computational aeroelasticity

A transonic unsteady aerodynamic and aeroelasticity code called CAP-TSD was developed for application to realistic aircraft configurations. The code permits the calculation of steady and unsteady flows about complete aircraft configurations for aeroelastic analysis in the flutter critical transonic speed range. The CAP-TSD code uses a time accurate approximate factorization algorithm for solution of the unsteady transonic small disturbance potential equation. An overview is given of the CAP-TSD code development effort and results are presented which demonstrate various capabilities of the code. Calculations are presented for several configurations including the General Dynamics 1/9 scale F-16 aircraft model and the ONERA M6 wing. Calculations are also presented from a flutter analysis of a 45 deg sweptback wing which agrees well with the experimental data. Descriptions are presented of the CAP-TSD code and algorithm details along with results and comparisons which demonstrate these recent developments in transonic computational aeroelasticity.

Batina, John T.↗

A benchmark for galactic cosmic-ray transport codes

A nontrivial analytic benchmark solution for galactic cosmic ray transport is presented for use in transport code validation. Computational accuracy for a previously-developed cosmic ray transport code is established to within one percent by comparison with this exact benchmark. Hence, solution accuracy for the transport problem is mainly limited by inaccuracies in the input spectra, input interaction databases, and the use of a straight ahead/velocity-conserving approximation.

Wilson, John W.↗

Steady and unsteady transonic small disturbance analysis of realistic aircraft configurations

A transonic unsteady aerodynamic and aeroelastic code called CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) was developed for application to realistic aircraft configurations. It permits the calculation of steady and unsteady flows about complete aircraft configurations for aeroelastic analysis of the flutter critical transonic speed range. The CAP-TSD code uses a time accurate approximate factorization algorithm for solution of the unsteady transonic small disturbance potential equation. An overview is given of the CAP-TSD code development effort along with recent algorithm modifications which are listed and discussed. Calculations are presented for several configurations including the General Dynamics 1/9th scale F-16C aircraft model to evaluate the algorithm and hence the reliability of the CAP-TSD code in general. Calculations are also presented for a flutter analysis of a 45 deg sweptback wing which agree well with the experimental data. Descriptions are presented of the CAP-TSD code and algorithm details along with results and comparisons which demonstrate the stability, accuracy, efficiency, and utility of CAP-TSD.

Batina, John T.↗

Recent advances in transonic computational aeroelasticity

A transonic unsteady aerodynamic and aeroelasticity code called CAP-TSD was developed for application to realistic aircraft configurations. The code permits the calculation of steady and unsteady flows about complete aircraft configurations for aeroelastic analysis in the flutter critical transonic speed range. The CAP-TSD code uses a time accurate approximate factorization algorithm for solution of the unsteady transonic small disturbance potential equation. An overview is given of the CAP-TSD code development effort and results are presented which demonstrate various capabilities of the code. Calculations are presented for several configurations including the General Dynamics 1/9 scale F-16 aircraft model and the ONERA M6 wing. Calculations are also presented from a flutter analysis of a 45 deg sweptback wing which agrees well with the experimental data. Descriptions are presented of the CAP-TSD code and algorithm details along with results and comparisons which demonstrate these recent developments in transonic computational aeroelasticity.

Batina, John T.↗

Adaptive grid embedding for the two-dimensional flux-split Euler equations

A numerical algorithm is presented for solving the 2-D flux-split Euler equations using a multigrid method with adaptive grid embedding. The method uses an unstructured data set along with a system of pointers for communication on the irregularly shaped grid topologies. An explicit two-stage time advancement scheme is implemented. A multigrid algorithm is used to provide grid level communication and to accelerate the convergence of the solution to steady state. Results are presented for a subcritical airfoil and a transonic airfoil with 3 levels of adaptation. Comparisons are made with a structured upwind Euler code which uses the same flux integration techniques of the present algorithm. Good agreement is obtained with converged surface pressure coefficients. The lift coefficients of the adaptive code are within 2 1/2 percent of the structured code for the sub-critical case and within 4 1/2 percent of the structured code for the transonic case using approximately one-third the number of grid points.

Warren, Gary Patrick↗

Simulation of real-gas effects on pressure distributions for a proposed Aeroassist Flight Experiment vehicle and comparison to prediction

Pressure distributions measured on a 60-deg elliptic cone, raked off at a 73-deg angle and having an ellipsoid nose (ellipticity equal to 2.0), are presented for a range of angle of attack from -10 to 15 deg. The high normal shock density ratio aspect of a real gas was simulated by testing in Mach-6 air (normal shock density ratio equal to 5.25) and Mach 6 CF4 (normal shock density ratio equal to 12.0). The effects of Reynolds number, angle of attack, and normal shock density ratio on these measurements are examined, and comparisons are made to an inviscid flowfield computer code known as HALIS. A significant effect of density ratio on pressure distributions on the cone section of the configuration was observed; the magnitude of this effect decreased with increased angle of attack. The effect of Reynolds number on pressures was negligible for forebody pressure distributions, but a measurable effect was noted on base pressures. In general, the HALIS code accurately predicted the measured pressure distributions in air and CF4.

Micol, John R.↗

The NASA aircraft noise prediction program improved propeller analysis system

The improvements and the modifications of the NASA Aircraft Noise Prediction Program (ANOPP) and the Propeller Analysis System (PAS) are described. Comparisons of the predictions and the test data are included in the case studies for the flat plate model in the Boundary Layer Module, for the effects of applying compressibility corrections to the lift and pressure coefficients, for the use of different weight factors in the Propeller Performance Module, for the use of the improved retarded time equation solution, and for the effect of the number grids in the Transonic Propeller Noise Module. The DNW tunnel test data of a propeller at different angles of attack and the Dowty Rotol data are compared with ANOPP predictions. The effect of the number of grids on the Transonic Propeller Noise Module predictions and the comparison of ANOPP TPN and DFP-ATP codes are studied. In addition to the above impact studies, the transonic propeller noise predictions for the SR-7, the UDF front rotor, and the support of the enroute noise test program are included.

Nguyen, L. Cathy↗