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 775 records · Page 43

Numerical solution to the glancing sidewall oblique shock wave/turbulent boundary layer interaction in three-dimension

A supersonic three-dimensional viscous forward-marching computer design code called PEPSIS is used to obtain a numerical solution of the three-dimensional problem of the interaction of a glancing sidewall oblique shock wave and a turbulent boundary layer. Very good results are obtained for a test case that was run to investigate the use of the wall-function boundary-condition approximation for a highly complex three-dimensional shock-boundary layer interaction. Two additional test cases (coarse mesh and medium mesh) are run to examine the question of near-wall resolution when no-slip boundary conditions are applied. A comparison with experimental data shows that the PEPSIS code gives excellent results in general and is practical for three-dimensional supersonic inlet calculations.

Anderson, B. H.↗

DSMC Simulations of Hypersonic Flows and Comparison With Experiments

This paper presents computational results obtained with the direct simulation Monte Carlo (DSMC) method for several biconic test cases in which shock interactions and flow separation-reattachment are key features of the flow. Recent ground-based experiments have been performed for several biconic configurations, and surface heating rate and pressure measurements have been proposed for code validation studies. The present focus is to expand on the current validating activities for a relatively new DSMC code called DS2V that Bird (second author) has developed. Comparisons with experiments and other computations help clarify the agreement currently being achieved between computations and experiments and to identify the range of measurement variability of the proposed validation data when benchmarked with respect to the current computations. For the test cases with significant vibrational nonequilibrium, the effect of the vibrational energy surface accommodation on heating and other quantities is demonstrated.

Moss, James N.↗

Wall functions for the kappa-epsilon turbulence model in generalized nonorthogonal curvilinear coordinates

A k-epsilon turbulence model suitable for compressible flow, including the new wall function formulation, has been incorporated into an existing compressible Reynolds-averaged Navier-Stokes code, F3D. The low Reynolds number k-epsilon model of Chien (1982) was added for comparison with the present method. A number of features were added to the F3D code including improved far-field boundary conditions and viscous terms in the streamwise direction. A series of computations of increasing complexity was run to test the effectiveness of the new formulation. Flow over a flat plate was computed by using both orthogonal and nonorthogonal grids, and the friction coefficients and velocity profiles compared with a semi-empirical equation. Flow over a body of revolution at zero angle of attack was then computed to test the method's ability to handle flow over a curved surface. Friction coefficients and velocity profiles were compared to test data. All models gave good results on a relatively fine grid, but only the wall function formulation was effective with coarser grids. Finally, in order to demonstrate the method's ability to handle complex flow fields, separated flow over a prolate spheroid at angle of attack was computed, and results were compared to test data. The results were also compared to a k-epsilon model by Kim and Patel (1991), in which one equation model patched in at the wall was employed. Both models gave reasonable solutions, but improvement is required for accurate prediction of friction coefficients in the separated regions.

Sondak, D. L.↗

Calibration of the NOAA 11 solar backscatter ultraviolet (SBUV/2) ozone data set from 1989 to 1993 using in-flight calibration data and SSBUV

Total ozone and ozone profiles are currently being measured by solar backscatter ultraviolet (SBUV/2) instruments onboard NOAA polar orbiting spacecraft using the backscattered ultraviolet technique. The NOAA 11 SBUV/2 operational data set, reprocessed from January 1989 to May 1993 (version 6), includes an updated algorithm and revised prelaunch and postlaunch calibrations of the geometrical albedo observations used to derive ozone values. The postlaunch calibration revisions remove time dependent errors in the ozone amounts due to instrument drift, while the revised prelaunch calibration corrects the absolute value of retrieved ozone. The prelaunch corrections are a result of calibration checks from in-orbit comparisons of ultraviolet geometric albedos measured by shuttle SBUV (SSBUV) and the NOAA 11 SBUV/2. Geometric albedo comparison data are further corrected using a radiative transfer code to account for the small difference in observing conditions between the two spacecraft. The postlaunch corrections rely on in-flight calibration and solar irradiance data to account for time dependent changes in instrument gain, thermal response, and instrument diffuser degradation over time. Comparison of data from three SSBUV flights, which occurred about one year apart, with concurrent SBUV/2 data provided an independent check of the time dependent change derived from the in-flight calibration data. Time independent corrections result in an increase of about 1% for total ozone, 5% for ozone at 1 mbar, and near 0% at 15 mbar. The time dependent corrections amount to an increase of 2% for total ozone, 10% for ozone near 1 mbar, and 3% at 15 mbar at the end of the current record in May 1993. Recent laboratory studies indicate that the absolute radiance calibrations may still be in error by a few percent which results in ozone profile values that are too low. The SBUV/2 total and ozone profile data are compared to the Nimbus SBUV data during the period when the data overlapped. Total ozone values agree to about 1%, while ozone profile differences range from -4% to +6%, depending on latitude and altitude, relative to SBUV. These differences are not statistically significant given the uncertainties of the two data sets.

Hilsenrath, E.↗

Low bit rate coding of Earth science images

In this paper, the authors discuss compression based on some new ideas in vector quantization and their incorporation in a sub-band coding framework. Several variations are considered, which collectively address many of the individual compression needs within the earth science community. The approach taken in this work is based on some recent advances in the area of variable rate residual vector quantization (RVQ). This new RVQ method is considered separately and in conjunction with sub-band image decomposition. Very good results are achieved in coding a variety of earth science images. The last section of the paper provides some comparisons that illustrate the improvement in performance attributable to this approach relative the the JPEG coding standard.

Kossentini, Faouzi↗

Two-dimensional computational analysis of a transport high-lift system and a comparison with flight-test results

Two currently available coupled inviscid/viscous multielement computational codes, including a relatively simple panel method and an Euler method, are used to analyze a high-lift system. The results are compared with two-dimensional wind-tunnel test results and then with the three-dimensional flight-test results obtained from the NASA Langley Transport Systems Research Vehicle five-element high-lift wing section. Comparisons were also made between the panel method, the Euler method, and flight data for two high-lift configurations, one representing a take-off configuration and the other an approach configuration. For the take-off configuration, both codes agreed reasonably well with experimental data, but both codes were found to overpredict the flap upper-surface pressures for the approach configuration.

Hardin, Jay D.↗

Experimental and computational ice shapes and resulting drag increase for a NACA 0012 airfoil

Tests were conducted in the Icing Research Tunnel (IRT) at LeRC to document the repeatability of the ice shape over the range of temperatures varying from -15 to 28 F. Measurements of drag increase due to the ice accretion were also made. The ice shape and drag coefficient data, with varying total temperatures at two different airspeeds, were compared with the computational predictions. The calculations were made with the 2D LEWICE/IBL code which is a combined code of LEWICE and the interactive boundary layer method developed for iced airfoils. Comparisons show good agreement with the experimental data in ice shapes. The calculations show the ability of the code to predict drag increases as the ice shape changes from a rime shape to a glaze shape.

Shin, Jaiwon↗

Experimental and computational ice shapes and resulting drag increase for a NACA 0012 airfoil

Tests were conducted in the Icing Research Tunnel (IRT) at the NASA Lewis Research Center to document the repeatability of the ice shape over the range of temperatures varying from -15 F to 28 F. Measurements of drag increase due to the ice accretion were also made. The ice shape and drag coefficient data, with varying total temperatures at two different airspeeds, were compared with the computational predictions. The calculations were made with the 2D LEWICE/IBL code which is a combined code of LEWICE and the interactive boundary layer method developed for iced airfoils. Comparisons show good agreement with the experimental data in ice shapes. The calculations show the ability of the code to predict drag increases as the ice shape changes from a rime shape to a glaze shape.

Shin, Jaiwon↗

Transonic flow about a thick circular-arc airfoil

An experimental and theoretical study of transonic flow over a thick airfoil, prompted by a need for adequately documented experiments that could provide rigorous verification of viscous flow simulation computer codes, is reported. Special attention is given to the shock-induced separation phenomenon in the turbulent regime. Measurements presented include surface pressures, streamline and flow separation patterns, and shadowgraphs. For a limited range of free-stream Mach numbers the airfoil flow field is found to be unsteady. Dynamic pressure measurements and high-speed shadowgraph movies were taken to investigate this phenomenon. Comparisons of experimentally determined and numerically simulated steady flows using a new viscous-turbulent code are also included. The comparisons show the importance of including an accurate turbulence model. When the shock-boundary layer interaction is weak the turbulence model employed appears adequate, but when the interaction is strong, and extensive regions of separation are present, the model is inadequate and needs further development.

Mcdevitt, J. B.↗

Comparison of Irradiated TRISO Fuel Radioactivity from Multiple Advanced Reactor Designs

The SCALE code system was used to model, deplete, and compare several different TRISO-fueled reactor designs: a helium-cooled prismatic reactor, a helium-cooled pebble-bed reactor (PBR), and a Fluoride-Lithium-Beryllium (FLIBE) molten-salt-cooled PBR. The purpose of this comparison was to understand how differences in the reactor designs affect the radioactivity of the fuel after discharge and whether those differences are significant. First, the various reactor designs were build and depleted in the TRITON module for each design and fuel enrichment. Then, the TRITON outputs were used to create burn-up dependent reactor libraries. These libraries were then used by ORIGEN to determine the activities of discharged fuel for each reactor, which were compared to generic Westinghouse 17x17 fuel. Overall, the results showed that short term decays are dominated by reactors with higher operating powers, and the reactor type, initial fuel enrichment, and maximum burn-up are of only secondary importance. Although this analysis only focuses on activities in becquerels, these dependencies are consistent with the expected behavior of decay heat. However, analysis of long- term time periods post-irradiation shows that the reactor type and maximum burn-up have strong impacts on the activities; initial fuel enrichment has a secondary impact while operating power is inconsequential. These results would be useful for analyses, such as dose assessment and modeling in post- release scenarios; normal fuel handling operations; and spent fuel transport, storage and disposal. Of particular interest, the results in this report show that analyses that focus on spent nuclear fuel of advanced reactors need to consider each parameter carefully. Unsurprisingly, if the correct operating power is not used in short term analyses, the results will not be correct. Perhaps unexpectedly, however, if the correct reactor type is not used, then the long term results will also be incorrect, especially for areas such as permanent disposal. Even though this report focuses on the total activity of nuclear fuel, it provides initial results on the effects of various input parameters and also provides a framework to extend the work into other analyses of spent fuel from advanced reactors, especially those employing TRISO fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Computational and numerical analysis of hypersonic nozzle flows with comparisons to wind tunnel calibration data

Navier-Stokes (NS), parabolized NS codes, and classical method of characteristics prediction techniques are employed to compute hypersonic nozzle flowfields for the 15 in Mach 6 High Temperature Tunnel and the 16 in Mach 17 Nitrogen Tunnel, which are part of the Langley Research Center. The study focuses on defining the accuracy of these techniques by comparing the computational results to wind tunnel pitot pressure measurements conducted in the newly designed nozzle of the Mach 6 tunnel and in the original nozzle of the Mach 17 tunnel. The comparisons between predicted results and test section pitot surveys for the Mach 6 nozzle are in good agreement and show highly uniform flow over a range of reservoir pressures and temperatures.

Hackett, Charles M.↗

Verification and validation of the open-source plasma fluid code: Zapdos

Zapdos is an open-source finite element plasma fluid solver based on the MOOSE multiphysics framework. This paper outlines Zapdos verification, benchmarking, and validation efforts for 1D and 2D RF capacitively coupled plasma discharge models for mid-range pressures (0.1 - 1 Torr). The verification process involved using the method of manufactured solutions to assess Zapdos spatial and temporal error convergence. L2 errors ranged from 10 -2 to 10 -4 , while the convergence’s slope were in agreement with the predicted slopes for the tested variable and time integration orders. The benchmarking process involved comparisons to previously results from the validated finite element code, LSODI. These works included 1D and 2D simulations for a range of plasma parameters (densities, temperatures, voltage, etc.). For the 1D cases, Zapdos and LSODI results were in very good agreement. In the 2D cases, variable behaviors matched, with slight discrepancies in peak values. The validation process involved comparisons to experimental works including electron density measurements by microwave interferometry and metastable density measurements by planar laser-induced fluorescence imaging. Results shown reasonable agreement at higher pressure, with results starting to diverge at low pressures. Furthermore, probable causes for this diverges are the limitation of the fluid assumption for plasmas at low pressure, or the need for more robust boundary conditions. Overall, Zapdos shown reasonable results for the verification, benchmarking, and validation efforts, and Zapdos can be downloaded at https://github.com/shannon-lab/zapdos.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Protostellar formation in rotating interstellar clouds. IV Nonisothermal collapse

Radiative transfer in the Eddington approximation is included in a multidimensional, self-gravitational, hydrodynamical computer code. Details of the numerical solution and thermodynamic relations are given. Comparison calculations with previous spherically symmetrical models of protostellar collapse are used to validate the basic approach and the artifices which allow the explicit hydrodynamics code to follow the accretion of gas onto a quasi-equilibrium core. A series of axisymmetric models is used to investigate the importance of rotation in collapsing clouds, as the initial amount of angular momentum is lowered, with an emphasis on the possible formation of rings. Rings readily form even in the nonisothermal regime except for very low initial angular momenta; even these clouds may experience ring formation prior to reaching stellar densities. The models imply that other effects (such as gravitational torques or turbulent viscosity) may be necesary to avoid binary formation and thus result in a presolar nebula consistent with the assumptions of either Safronov or Cameron.

Boss, A. P.↗

A Computational Study of an Oscillating VR-12 Airfoil with a Gurney Flap

Computations of the flow over an oscillating airfoil with a Gurney-flap are performed using a Reynolds Averaged Navier-Stokes code and compared with recent experimental data. The experimental results have been generated for different sizes of the Gurney flaps. The computations are focused mainly on a configuration. The baseline airfoil without a Gurney flap is computed and compared with the experiments in both steady and unsteady cases for the purpose of initial testing of the code performance. The are carried out with different turbulence models. Effects of the grid refinement are also examined and unsteady cases, in addition to the assessment of solver effects. The results of the comparisons of steady lift and drag computations indicate that the code is reasonably accurate in the attached flow of the steady condition but largely overpredicts the lift and underpredicts the drag in the higher angle steady flow.

Rhee, Myung↗

Closed coronal structures. V - Gasdynamic models of flaring loops and comparison with SMM observations

A time-dependent one-dimensional code incorporating energy, momentum and mass conservation equations, and taking the entire solar atmospheric structure into account, is used to investigate the hydrodynamic response of confined magnetic structures to strong heating perturbations. Model calculation results are compared with flare observations which include the light curves of spectral lines formed over a wide range of coronal flare temperatures, as well as determinations of Doppler shifts for the high temperature plasma. It is shown that the numerical simulation predictions are in good overall agreement with the observed flare coronal plasma evolution, correctly reproducing the temporal profile of X-ray spectral lines and their relative intensities. The predicted upflow velocities support the interpretation of the blueshifts as due to evaporation of chromospheric material.

Peres, G.↗

Comparisons of Flutter Analyses for an Experimental Fan

Two propulsion aeroelasticity codes were used to model the aeroelastic characteristics of an experimental forward-swept fan that encountered flutter during wind tunnel testing. Both of these three-dimensional codes model the unsteady flowfield due to blade vibrations using the Navier-Stokes equations. In the first approach, the unsteady flow equations are solved using an implicit time-marching approach. In the second approach, the unsteady flow equations are converted to a harmonic balance form and solved using a pseudo-time marching method. This paper describes the flutter calculations and compares the results to experimental measurements.

Bakhle, Milind A.↗

Dynamic Ground Effects Simulation Using OVERFLOW-D

This presentation is broken into 5 logical sections. The Background Information section describes the technical issues being address by this study. The Approach section describes the organization of the contract effort which was laid out as the most effective means of quantifying, with validated methods, the magnitude of dynamic ground effects for the TCA (Technology Concept Aircraft) configuration. The Validation Case section describes the analysis of the XB-70 configuration in both static and dynamic ground effect, with comparisons to wind tunnel and flight test data. The TCA Analysis section then describes the application of the same codes and methodologies to the TCA in both static and dynamic ground effect. Comparisons are made between the static and dynamic, as well as to early static data from a recent wind tunnel test on the TCA configuration. Finally, the work to date is summarized and the future direction of this study is outlined.

Dwyer, Bill↗

Performance analysis of three dimensional integral equation computations on a massively parallel computer

The purpose of this study is to investigate the performance of the integral equation computations using numerical source field-panel method in a massively parallel processing (MPP) environment. A comparative study of computational performance of the MPP CM-5 computer and conventional Cray-YMP supercomputer for a three-dimensional flow problem is made. A serial FORTRAN code is converted into a parallel CM-FORTRAN code. Some performance results are obtained on CM-5 with 32, 62, 128 nodes along with those on Cray-YMP with a single processor. The comparison of the performance indicates that the parallel CM-FORTRAN code near or out-performs the equivalent serial FORTRAN code for some cases.

Logan, Terry G.↗