Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code Validation”

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

The Project for Intercomparison of Land-surface Parameterization Schemes

The Project for Intercomparison of Land-surface Parameterization Schemes (PILPS) is described and the first stage science plan outlined. PILPS is a project designed to improve the parameterization of the continental surface, especially the hydrological, energy, momentum, and carbon exchanges with the atmosphere. The PILPS Science Plan incorporates enhanced documentation, comparison, and validation of continental surface parameterization schemes by community participation. Potential participants include code developers, code users, and those who can provide datasets for validation and who have expertise of value in this exercise. PILPS is an important activity because existing intercomparisons, although piecemeal, demonstrate that there are significant differences in the formulation of individual processes in the available land surface schemes. These differences are comparable to other recognized differences among current global climate models such as cloud and convection parameterizations. It is also clear that too few sensitivity studies have been undertaken with the result that there is not yet enough information to indicate which simplifications or omissions are important for the near-surface continental climate, hydrology, and biogeochemistry. PILPS emphasizes sensitivity studies with and intercomparisons of existing land surface codes and the development of areally extensive datasets for their testing and validation.

Henderson-Sellers, A.↗

Fostering Better Collaboration in Software Development Cycles Between Scientists and Programmers to Ensure the Integrity of and Promote the Development of New Scientific Data Products.

Misaligned incentives lead to reduced interaction between scientists and programmers on modern NASA science data-product development teams. Typically, situations arise where the scientist is not incentivized to learn modern coding practices and the programmer does not understand the science algorithms in the code. A programmer is responsible for the deliverable code thus setting a tradeoff between the desire for code improvement versus fear of compromising the integrity of data-product while the scientist continues to rely on their legacy codebases owing to the complexity of using the delivered code outside the processing environment and lack of validation modules. The NASA/CERES-TISA project has adopted a collaborative approach, with scientists and programmers both utilizing the same software repository with multiple branches, some optimized for delivery to a processing datacenter and others for scientific product development and validation. A team of scientists and programmers jointly review any new science code updates for integration into the codebase and strive to improve practices through promoting algorithm understanding, better institutional knowledge exchange and documentation, modularization, and developing data processing flow-dictated validation and debugging methods. This leads to a reduction in the personnel single point failures and reduced development time for creation of new science data-products.

CERES↗

The Crucial Role of Error Correlation for Uncertainty Modeling of CFD-Based Aerodynamics Increments

The Ares I ascent aerodynamics database for Design Cycle 3 (DAC-3) was built from wind-tunnel test results and CFD solutions. The wind tunnel results were used to build the baseline response surfaces for wind-tunnel Reynolds numbers at power-off conditions. The CFD solutions were used to build increments to account for Reynolds number effects. We calculate the validation errors for the primary CFD code results at wind tunnel Reynolds number power-off conditions and would like to be able to use those errors to predict the validation errors for the CFD increments. However, the validation errors are large compared to the increments. We suggest a way forward that is consistent with common practice in wind tunnel testing which is to assume that systematic errors in the measurement process and/or the environment will subtract out when increments are calculated, thus making increments more reliable with smaller uncertainty than absolute values of the aerodynamic coefficients. A similar practice has arisen for the use of CFD to generate aerodynamic database increments. The basis of this practice is the assumption of strong correlation of the systematic errors inherent in each of the results used to generate an increment. The assumption of strong correlation is the inferential link between the observed validation uncertainties at wind-tunnel Reynolds numbers and the uncertainties to be predicted for flight. In this paper, we suggest a way to estimate the correlation coefficient and demonstrate the approach using code-to-code differences that were obtained for quality control purposes during the Ares I CFD campaign. Finally, since we can expect the increments to be relatively small compared to the baseline response surface and to be typically of the order of the baseline uncertainty, we find that it is necessary to be able to show that the correlation coefficients are close to unity to avoid overinflating the overall database uncertainty with the addition of the increments.

Hemsch, Michael J.↗

Benchmark Active Controls Technology (BACT) Wing CFD Results

The Benchmark Active Controls Technology (BACT) wing test (see chapter 8E) provides data for the validation of aerodynamic, aeroelastic, and active aeroelastic control simulation codes. These data provide a rich database for development and validation of computational aeroelastic and aeroservoelastic methods. In this vein, high-level viscous CFD analyses of the BACT wing have been performed for a subset of the test conditions available in the dataset. The computations presented in this section investigate the aerodynamic characteristics of the rigid clean wing configuration as well as simulations of the wing with a static and oscillating aileron and spoiler deflection. Two computational aeroelasticity codes extensively used at NASA Langley Research Center are implemented in this simulation. They are the ENS3DAE and CFL3DAE computational aeroelasticity programs. Both of these methods solve the three-dimensional compressible Navier-Stokes equations for both rigid and flexible vehicles, but they use significantly different approaches to the solution 6f the aerodynamic equations of motion. Detailed descriptions of both methods are presented in the following section.

Schuster, David M.↗

Navier-Stokes analysis of a liquid rocket engine disk cavity

This paper presents a Navier-Stokes analysis of hydrodynamic phenomena occurring in the aft disk cavity of a liquid rocket engine turbine. The cavity analyzed in the Space Shuttle Main Engine Alternate Turbopump currently being developed by NASA and Pratt and Whitney. Comparison of results obtained from the Navier-Stokes code for two rotating disk datasets available in the literature are presented as benchmark validations. The benchmark results obtained using the code show good agreement relative to experimental data, and the turbine disk cavity was analyzed with comparable grid resolution, dissipation levels, and turbulence models. Predicted temperatures in the cavity show that little mixing of hot and cold fluid occurs in the cavity and the flow is dominated by swirl and pumping up the rotating disk.

Benjamin, Theodore G.↗

Development of an Aeroelastic Code Based on an Euler/Navier-Stokes Aerodynamic Solver

This paper describes the development of an aeroelastic code (TURBO-AE) based on an Euler/Navier-Stokes unsteady aerodynamic analysis. A brief review of the relevant research in the area of propulsion aeroelasticity is presented. The paper briefly describes the original Euler/Navier-Stokes code (TURBO) and then details the development of the aeroelastic extensions. The aeroelastic formulation is described. The modeling of the dynamics of the blade using a modal approach is detailed, along with the grid deformation approach used to model the elastic deformation of the blade. The work-per-cycle approach used to evaluate aeroelastic stability is described. Representative results used to verify the code are presented. The paper concludes with an evaluation of the development thus far, and some plans for further development and validation of the TURBO-AE code.

Bakhle, Milind A.↗

Investigation of mixing in a turbofan exhaust duct. II Computer code application and verification

A three-dimensional analysis of turbofan forced mixer nozzle aerodynamics demonstrates that the complex flow structure is dominated by geometrically induced secondary flow rather than by turbulence. The test apparatus consisted of a fixed upstream model section and a rotating shroud. The Mach number of the fan and core streams at the mixing plane (lobe exit) was 0.45, the bypass ratio was about 4, and the Reynolds number based on the shroud radius was 1,100,000. The three velocity components near the exit plane of the lobes were measured using flow angularity probes to provide information about the mixer inflow conditions for turbulent computations. The validity of a previous computer code was demonstrated in a comparison of the nozzle exit temperature data with the computed temperature distributions. The mechanism most responsible for the generation of secondary flow within the lobes is due to the turning of the fan and core streams in opposite radial directions.

Povinelli, L. A.↗

Recent Updates to the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1

The Object Reentry Survival Analysis Tool (ORSAT) code is maintained and used by the NASA Orbital Debris Program Office (ODPO) and has been under continuous development and improvement since the mid-1990s. ORSAT is an object-oriented reentry simulation tool; it models a satellite as a collection of discrete components that follow independent trajectories upon the breakup of the parent object. Version 7.1 of the tool incorporates five years of new thermal and aerodynamic model development, multi-processor parametric study capability, codebase upgrades, and numerous bug-fixes. The thermal demise model was completely rewritten using a forward-time/central-space numerical stencil and incorporating a new pyrolysis model for fiber-reinforced plastic (FRP) materials. New aerodynamic and aeroheating models for hollow cylinders and hollow square prisms were developed using a combination of flow simulations in the direct simulation Monte Carlo (DSMC) Analysis Code (DAC) and Data Parallel Line Relaxation (DPLR) code and free-flight tests in the University of Texas at San Antonio’s Hypersonic Wind Tunnel. The latest version also incorporates a mechanical, strength-based demise model for FRP materials. Minor improvements include an update to the Fortran 2018 codebase; improved integration and speed with the Python-based, multi-core, parametric study tool, AutoORSAT; and fixes for many minor bugs. This new version of ORSAT will enable more accurate reentry risk assessments for modern satellites. This paper presents an overview of these changes and a summary of the verification and validation performed on the final code.

Benton R. Greene↗

Recent Updates to the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1

The Object Reentry Survival Analysis Tool (ORSAT) code is maintained and used by the NASA Orbital Debris Program Office (ODPO) and has been under continuous development and improvement since the mid-1990s. ORSAT is an object-oriented reentry simulation tool; it models a satellite as a collection of discrete components that follow independent trajectories upon the breakup of the parent object. Version 7.1 of the tool incorporates five years of new thermal and aerodynamic model development, multi-processor parametric study capability, codebase upgrades, and numerous bug-fixes. The thermal demise model was completely rewritten using a forward-time/central-space numerical stencil and incorporating a new pyrolysis model for fiber-reinforced plastic (FRP) materials. New aerodynamic and aeroheating models for hollow cylinders and hollow square prisms were developed using a combination of flow simulations in the direct simulation Monte Carlo (DSMC) Analysis Code (DAC) and Data Parallel Line Relaxation (DPLR) code and free-flight tests in the University of Texas at San Antonio’s Hypersonic Wind Tunnel. The latest version also incorporates a mechanical, strength-based demise model for FRP materials. Minor improvements include an update to the Fortran 2018 codebase; improved integration and speed with the Python-based, multi-core, parametric study tool, AutoORSAT; and fixes for many minor bugs. This new version of ORSAT will enable more accurate reentry risk assessments for modern satellites. This paper presents an overview of these changes and a summary of the verification and validation performed on the final code.

Benton R. Greene↗

High Speed Research (HSR) Multi-Year Summary Report for Calendar Years 1995-1999

The Aeroelasticity Task is intended to provide demonstrated technology readiness to predict and improve flutter characteristics of an HSCT configuration. This requires aerodynamic codes that are applicable to the wide range of flight regimes in which the HSCT will operate, and are suitable to provide the higher fidelity required for evaluation of aeroservoelastic coupling effects. Prediction of these characteristics will result in reduced airplane weight and risk associated with a highly flexible, low-aspect ratio supersonic airplane with narrow fuselage, relatively thin wings, and heavy engines. This Task is subdivided into three subtasks. The first subtask includes the design, fabrication, and testing of wind-tunnel models suitable to provide an experimental database relevant to HSCT configurations. The second subtask includes validation of candidate unsteady aerodynamic codes, applicable in the Mach and frequency ranges of interest for the HSCT, through analysis test correlation with the test data. The third subtask includes efforts to develop and enhance these codes for application to HSCT configurations. The wind tunnel models designed and constructed during this program furnished data which were useful for the analysis test correlation work but there were shortcomings. There was initial uncertainty in the proper tunnel configuration for testing, there was a need for higher quality measured model geometry, and there was a need for better measured model displacements in the test data. One of the models exhibited changes in its dynamic characteristics during testing. Model design efforts were hampered by a need for more and earlier analysis support and better knowledge of material properties. Success of the analysis test correlation work was somewhat muted by the uncertainties in the wind tunnel model data. The planned extent of the test data was not achieved, partly due to the delays in the model design and fabrication which could not be extended due to termination of the HSR program.

Baker, Myles↗

Potential flow calculations and preliminary wing design in support of an NLF variable sweep transition flight experiment

NASA Langley and NASA Ames-Dryden have defined a variable-sweep transition-flight experiment utilizing the F-14 aircraft to enhance understanding of the interaction of crossflow and Tollmien-Schlichting instabilities on a laminar-boundary-layer transition. The F-14 wing outer panel will be modified to generate favorable pressure gradients on the upper wing surface over a wide range of flight conditions. Extensive computations have been performed using two-dimensional and three-dimensional transonic analysis codes. Flight-test and computational data are compared and shown to validate the applicability of the three-dimensional codes (WBPPW and TAWFIVE). In addition, results from two preliminary glove designs derived from two different approaches to the design problem are presented. Advantages and disadvantages of each approach are identified, and it is concluded that coupling an analysis code with an automated design procedure yields a powerful code with distinct advantages over a 'cut-and-dry' approach.

Waggoner, E. G.↗

Computation of Sound Propagation by Boundary Element Method

This report documents the development of a Boundary Element Method (BEM) code for the computation of sound propagation in uniform mean flows. The basic formulation and implementation follow the standard BEM methodology; the convective wave equation and the boundary conditions on the surfaces of the bodies in the flow are formulated into an integral equation and the method of collocation is used to discretize this equation into a matrix equation to be solved numerically. New features discussed here include the formulation of the additional terms due to the effects of the mean flow and the treatment of the numerical singularities in the implementation by the method of collocation. The effects of mean flows introduce terms in the integral equation that contain the gradients of the unknown, which is undesirable if the gradients are treated as additional unknowns, greatly increasing the sizes of the matrix equation, or if numerical differentiation is used to approximate the gradients, introducing numerical error in the computation. It is shown that these terms can be reformulated in terms of the unknown itself, making the integral equation very similar to the case without mean flows and simple for numerical implementation. To avoid asymptotic analysis in the treatment of numerical singularities in the method of collocation, as is conventionally done, we perform the surface integrations in the integral equation by using sub-triangles so that the field point never coincide with the evaluation points on the surfaces. This simplifies the formulation and greatly facilitates the implementation. To validate the method and the code, three canonic problems are studied. They are respectively the sound scattering by a sphere, the sound reflection by a plate in uniform mean flows and the sound propagation over a hump of irregular shape in uniform flows. The first two have analytical solutions and the third is solved by the method of Computational Aeroacoustics (CAA), all of which are used to compare the BEM solutions. The comparisons show very good agreements and validate the accuracy of the BEM approach implemented here.

Guo, Yueping↗

Foam on Tile Impact Modeling for the STS-107 Investigation

Following the breakup of the Space Shuttle Columbia during reentry a NASA/Contractor investigation team was formed to examine the probable damage inflicted on Orbiter Thermal Protection System elements by impact of External Tank insulating foam projectiles. The authors formed a working subgroup within the larger team to apply the Smooth Particle Hydrodynamics code SPHC to the damage estimation problem. Numerical models of the Orbiter's tiles and of the Tank's foam were constructed and used as inputs into the code. Material properties needed to properly model the tiles and foam were obtained from other working subgroups who performed tests on these items for this purpose. Two- and three-dimensional models of the tiles were constructed, including the glass outer layer, the main body of LI-900 insulation, the densified lower layer of LI-900, the Nomex felt mounting layer, and the Aluminum 2024 vehicle skin. A model for the BX-250 foam including porous compression, elastic rebound, and surface erosion was developed. Code results for the tile damage and foam behavior were extensively validated through comparison with Southwest Research Institute foam-on-tile impact experiments carried out in 1999. These tests involved small projectiles striking individual tiles and small tile arrays. Following code and model validation we simulated impacts of larger foam projectiles on the examples of tile systems used on the Orbiter. Results for impacts on the main landing gear door are presented in this paper, including effects of impacts at several angles, and of rapidly rotating projectiles. General results suggest that foam impacts on tiles at about 500 mph could cause appreciable damage if the impact angle is greater than about 20 degrees. Some variations of the foam properties, such as increased brittleness or increased density could increase damage in some cases. Rotation up to 17 rps failed to increase the damage for the two cases considered. This does not rule out other cases in which the rotational energy might lead to an increase in tile damage, but suggests that in most cases rotation will not be an important factor.

Stellingwerf, R. F.↗

Concentration of off-axis radiation by solar concentrators for space power

Off-axis radiation is radiation from any direction not parallel to the axis of the solar concentrator. It will be reflected to regions other than the focus of the parabolic concentrator and possibly concentrated there. Four types of off-axis radiation are discussed. These are: (1) small off-axis angles during walk-off; (2) large off-axis angles; (3) an extended off-axis source such as Earth albedo; and (4) miscellaneous off-axis sources including radio frequency sources and local point sources. A previous analytical study used a computer code named PIXEL to predict concentration of off-axis radiation and a previous experimental study used an 11-m diameter multifaceted dish concentrator to validate the PIXEL analysis. The PIXEL code was limited in that it represented concentration by an ideal parabolic reflector of light from a point source. Another code named OFFSET has been developed to represent the solar concentrator being developed for Space Station Freedom. It is a detailed, ray tracing model which represents 50 ray originating points on the Sun and reflections from 10 points on each of the 456 concentrator facets. Results of this code are generally similar to the PIXEL results although there are small differences due to the more detailed representations of the Sun and concentrator that were used in the OFFSET code.

Jefferies, Kent S.↗

Comparison and Validation of FLUKA and HZETRN as Tools for Investigating the Secondary Neutron Production in Large Space Vehicles

NASA's exploration goals are focused on deep space travel and Mars surface operations. To accomplish these goals, large structures will be necessary to transport crew and logistics in the initial stages, and NASA will need to keep the crew and the vehicle safe during transport and any surface activities. One of the major challenges of deep space travel is the space radiation environment and its impacts on the crew, the electronics, and the vehicle materials. The primary radiation from the sun (solar particle events) and from outside the solar system (galactic cosmic rays) interact with materials of the vehicle. These interactions lead to some of the primary radiation being absorbed, being modified, or producing secondary radiation (primarily neutrons). With all vehicles, the high energy primary radiation is of most concern. However, with larger vehicles that have large shielding masses, there is more opportunity for secondary radiation production, and this secondary radiation can be significant enough to cause concern. When considering surface operations, there is also a secondary radiation source from the surface of the planet, known as albedo, with neutrons being one of the most significant species. Given new vehicle designs for deep space and Mars missions, the secondary radiation environment and the implications of that environment is currently not well understood. Thus, several studies are necessary to fill the knowledge gaps of this secondary radiation environment. In this paper, we put forth the initial steps to increasing our understanding of neutron production from large vehicles by comparing the neutron production resulting from our radiation transport codes and providing a preliminary validation of our results against flight data. This paper will review the details of these results and discuss the finer points of the analysis.

Rojdev, Kristina↗

Three Dimensional Viscous Flow Field in an Axial Flow Turbine Nozzle Passage

The objective of this investigation is experimental and computational study of three dimensional viscous flow field in the nozzle passage of an axial flow turbine stage. The nozzle passage flow field has been measured using a two sensor hot-wire probe at various axial and radial stations. In addition, two component LDV measurements at one axial station (x/c(sum m) = 0.56) were performed to measure the velocity field. Static pressure measurements and flow visualization, using a fluorescent oil technique, were also performed to obtain the location of transition and the endwall limiting streamlines. A three dimensional boundary layer code, with a simple intermittency transition model, was used to predict the viscous layers along the blade and endwall surfaces. The boundary layers on the blade surface were found to be very thin and mostly laminar, except on the suction surface downstream of 70% axial chord. Strong radial pressure gradient, especially close to the suction surface, induces strong cross flow components in the trailing edge regions of the blade. On the end-walls the boundary layers were much thicker, especially near the suction corner of the casing surface, caused by secondary flow. The secondary flow region near the suction-casing surface corner indicates the presence of the passage vortex detached from the blade surface. The corner vortex is found to be very weak. The presence of a closely spaced rotor downstream (20% of the nozzle vane chord) introduces unsteadiness in the blade passage. The measured instantaneous velocity signal was filtered using FFT square window to remove the periodic unsteadiness introduced by the downstream rotor and fans. The filtering decreased the free stream turbulence level from 2.1% to 0.9% but had no influence on the computed turbulence length scale. The computation of the three dimensional boundary layers is found to be accurate on the nozzle passage blade surfaces, away from the end-walls and the secondary flow region. On the nozzle passage endwall surfaces the presence of strong pressure gradients and secondary flow limit the validity of the boundary layer code.

Ristic, D.↗

High spatial resolution measurements in a single stage ram accelerator

High spatial resolution experimental tube wall pressure measurements of ram accelerator gas dynamic phenomena are presented in this paper. The ram accelerator is a ramjet-in-tube device which operates in a manner similar to that of a conventional ramjet. The projectile resembles the centerbody of a ramjet and travels supersonically through a tube filled with a combustible gaseous mixture, with the tube acting as the outer cowling. Pressure data are recorded as the projectile passes by sensors mounted in the tube wall at various locations along the tube. Utilization of special highly instrumented sections of tube has allowed the recording of gas dynamic phenomena with high resolution. High spatial resolution tube wall pressure data from the three regimes of propulsion studied to date (subdetonative, transdetonative, and superdetonative) in a single stage gas mixture are presented and reveal the three-dimensional character of the flow field induced by projectile fins and the canting of the fins and the canting of the projectile body relative to the tube wall. Also presented for comparison to the experimental data are calculations made with an inviscid, three-dimensional CFD code. The knowledge gained from these experiments and simulations is useful in understanding the underlying nature of ram accelerator propulsive regimes, as well as assisting in the validation of three-dimensional CFD coded which model unsteady, chemically reactive flows.

Hinkey, J. B.↗

The Proteus Navier-Stokes code

An effort is currently underway at NASA Lewis to develop two- and three-dimensional Navier-Stokes codes, called Proteus, for aerospace propulsion applications. The emphasis in the development of Proteus is not algorithm development or research on numerical methods, but rather the development of the code itself. The objective is to develop codes that are user-oriented, easily-modified, and well-documented. Well-proven, state-of-the-art solution algorithms are being used. Code readability, documentation (both internal and external), and validation are being emphasized. This paper is a status report on the Proteus development effort. The analysis and solution procedure are described briefly, and the various features in the code are summarized. The results from some of the validation cases that have been run are presented for both the two- and three-dimensional codes.

Towne, Charles E.↗