Results for “code comparisons” — Page 36

NASA NTRS1993

The annual pressure cycle on Mars: Results from the LMD Martian atmospheric general circulation model

We have been developing a General Circulation Model (GCM) of the martian atmosphere since 1989. The model has been described rather extensively elsewhere and only the main characteristics are given here. The dynamical part of the model, adapted from the LMD terrestrial climate model, is based on a finite-difference formulation of the classical 'primitive equations of meteorology.' The radiative transfer code includes absorption and emission by CO2 (carefully validated by comparison to line-by-line calculations) and dust in the thermal range and absorption and scattering by dust in the visible range. Other physical parameterizations are included: modeling of vertical turbulent mixing, dry convective adjustment (in order to prevent vertical unstable temperature profiles), and a multilayer model of the thermal conduction in the soil. Finally, the condensation-sublimation of CO2 is introduced through specification of a pressure-dependent condensation temperature. The atmospheric and surface temperatures are prevented from falling below this critical temperature by condensation and direct precipitation onto the surface of atmospheric CO2. The only prespecified spatial fields are the surface thermal inertia, albedo, and topography.

Hourdin, Frederic↗

NASA NTRS1998

Summary of Inlet Characteristics of the F/A-18A High Alpha Research Vehicle

Effects of high-angle-of-attack flight on aircraft inlet aerodynamic characteristics were investigated at NASA Dryden Flight Research Center as part of NASA's High Alpha Technology Program. The highly instrumented F/A-18A High Alpha Research Vehicle was used for this research. A newly designed inlet total-pressure rake was installed in front of the right-hand F404-GE-400 engine to measure inlet recovery and distortion characteristics. Objectives included: (1) determining the inlet total-pressure characteristics at steady high-angle-of-attack conditions; (2) assessing if inlet distortion is significantly different between rapid angle-of-attack maneuvers and corresponding steady aerodynamic conditions; (3) assessing inlet characteristics during aircraft departures; (4) providing data for developing and verifying computational fluid dynamic codes; and (5) calculating engine airflow using four methods for comparison with a reference method. This paper describes the results obtained from this investigation. These data and the associated database were rigorously validated to establish the foundation for understanding inlet characteristics at high angle of attack.

Walsh, Kevin↗

NASA NTRS1984

The Prediction of Transonic Flows on an Advancing Rotor

The object of this investigation was to find and demonstrate a means of performing efficient finite-difference computations of rotor loading for a trimmed rotor in high-speed, forward flight. The essence of the scheme that was developed is a loose-coupled iteration procedure between a finite difference and a comprehensive integral rotor code. The coupling involves a transfer of appropriate load and inflow data on the advancing side between the two codes such that consistency maintained. Sample computations, including a limited comparison with model rotor data, are presented. The scheme converges rapidly. However, even one iteration with this scheme can provide sufficient accuracy for many purposes.

Tung, Chee↗

NASA NTRS1999

GHRS Observations of Cool, Low-Gravity Stars: The Outer Atmosphere and Wind of the Nearby K Supergiant Lambda Velorum - 5

UV spectra of lambda Velorum taken with the Goddard High Resolution Spectrograph (GHRS) on the Hubble Space Telescope are used to probe the structure of the outer atmospheric layers and wind and to estimate the mass-loss rate from this K5 lb-II supergiant. VLA radio observations at lambda = 3.6 cm are used to obtain an independent check on the wind velocity and mass-loss rate inferred from the UV observations, Parameters of the chromospheric structure are estimated from measurements of UV line widths, positions, and fluxes and from the UV continuum flux distribution. The ratios of optically thin C II] emission lines indicate a mean chromospheric electron density of log N(sub e) approximately equal 8.9 +/- 0.2 /cc. The profiles of these lines indicate a chromospheric turbulence (v(sub 0) approximately equal 25-36 km/s), which greatly exceeds that seen in either the photosphere or wind. The centroids of optically thin emission lines of Fe II and of the emission wings of self-reversed Fe II lines indicate that they are formed in plasma approximately at rest with respect to the photosphere of the star. This suggests that the acceleration of the wind occurs above the chromospheric regions in which these emission line photons are created. The UV continuum detected by the GHRS clearly traces the mean flux-formation temperature as it increases with height in the chromosphere from a well-defined temperature minimum of 3200 K up to about 4600 K. Emission seen in lines of C III] and Si III] provides evidence of material at higher than chromospheric temperatures in the outer atmosphere of this noncoronal star. The photon-scattering wind produces self-reversals in the strong chromospheric emission lines, which allow us to probe the velocity field of the wind. The velocities to which these self-absorptions extend increase with intrinsic line strength, and thus height in the wind, and therefore directly map the wind acceleration. The width and shape of these self-absorptions reflect a wind turbulence of approximately equal 9-21 km/s. We further characterize the wind by comparing the observations with synthetic profiles generated with the Lamers et al. Sobolev with Exact Integration (SEI) radiative transfer code, assuming simple models of the outer atmospheric structure. These comparisons indicate that the wind in 1994 can be described by a model with a wind acceleration parameter beta approximately 0.9, a terminal velocity of 29-33 km/s, and a mass-loss rate approximately 3 x 10(exp -9) solar M/yr. Modeling of the 3.6 cm radio flux observed in 1997 suggests a more slowly accelerating wind (higher beta) and/or a higher mass-loss rate than inferred from the UV line profiles. These differences may be due to temporal variations in the wind or from limitations in one or both of the models. The discrepancy is currently under investigation.

Carpenter, Kenneth G.↗

NASA NTRS1994

Investigation of Methods for the Structural Weight Analysis of a Mach 2.4 Axisymmetric Inlet

Structural design and analysis tools appropriate for estimating the structural weight of an axisymmetric inlet designed for Mach 2.4 cruise were evaluated. Little information regarding the inlet mechanical design is available in the preliminary design phase, so it is necessary to first develop a reasonable structural design before estimating the inlet weight. The Internally Pressurized Structure Synthesis and Optimization (IPSSO) program, employing an analytical approach, was chosen for evaluation due to its combined design and analysis capabilities. The inlet design produced by IPSSO was then analyzed using the NASTRAN finite element program. The finite element analysis was performed to help identify the limitations of the analytically based code as well as to evaluate NASTRAN for this application. Comparison between the IPSSO inlet weight and that of a similar inlet developed by the Boeing Commercial Airplane Group was also made. Program evaluation concluded that the combined use of IPSSO to create an initial design and NASTRAN to perform a numerical analysis would provide the capability to evaluate a limited number of inlet design The development of a new tool for the minimum weight design and analysis of inlet structures would be required for greater flexibility in evaluating inlet conceptual designs.

Nadell, Shari-Beth↗

NASA NTRS2002

Electrical Performance of the International Space Station U.S. Photovoltaic Array During Bifacial Illumination

With the first United States (U.S.) photovoltaic array (PVA) activated on International Space Station (ISS) in December 2000, on-orbit data can now be compared to analytical predictions. Due to ISS operational constraints, it is not always possible to point the front side of the arrays at the Sun. Thus, in many cases, sunlight directly illuminates the backside of the PVA as well as albedo illumination on either the front or the back. During this time, appreciable power is produced since the solar cells are mounted on a thin, solar transparent substrate. It is important to present accurate predictions for both front and backside power generation for mission planning, certification of flight readiness for a given mission, and on-orbit mission support. To provide a more detailed assessment of the ISS power production capability, the authors developed a PVA electrical performance model applicable to generalized bifacial illumination conditions. On-orbit PVA performance data were also collected and analyzed. This paper describes the ISS PVA performance model, and the methods used to reduce orbital performance data. Analyses were performed using SPACE. a NASA-GRC developed computer code for the ISS program office. Results showed a excellent comparison of on-orbit performance data and analytical results.

Delleur, Ann M.↗

NASA NTRS2004

Numerical Modeling of Active Flow Control in a Boundary Layer Ingesting Offset Inlet

This investigation evaluates the numerical prediction of flow distortion and pressure recovery for a boundary layer ingesting offset inlet with active flow control devices. The numerical simulations are computed using a Reynolds averaged Navier-Stokes code developed at NASA. The numerical results are validated by comparison to experimental wind tunnel tests conducted at NASA Langley Research Center at both low and high Mach numbers. Baseline comparisons showed good agreement between numerical and experimental results. Numerical simulations for the inlet with passive and active flow control also showed good agreement at low Mach numbers where experimental data has already been acquired. Numerical simulations of the inlet at high Mach numbers with flow control jets showed an improvement of the flow distortion. Studies on the location of the jet actuators, for the high Mach number case, were conducted to provide guidance for the design of a future experimental wind tunnel test.

Allan, Brian G.↗

NASA NTRS2005

Measurement and Prediction of the Thermomechanical Response of Shape Memory Alloy Hybrid Composite Beams

Previous work at NASA Langley Research Center (LaRC) involved fabrication and testing of composite beams with embedded, pre-strained shape memory alloy (SMA) ribbons within the beam structures. That study also provided comparison of experimental results with numerical predictions from a research code making use of a new thermoelastic model for shape memory alloy hybrid composite (SMAHC) structures. The previous work showed qualitative validation of the numerical model. However, deficiencies in the experimental-numerical correlation were noted and hypotheses for the discrepancies were given for further investigation. The goal of this work is to refine the experimental measurement and numerical modeling approaches in order to better understand the discrepancies, improve the correlation between prediction and measurement, and provide rigorous quantitative validation of the numerical analysis/design tool. The experimental investigation is refined by a more thorough test procedure and incorporation of higher fidelity measurements such as infrared thermography and projection moire interferometry. The numerical results are produced by a recently commercialized version of the constitutive model as implemented in ABAQUS and are refined by incorporation of additional measured parameters such as geometric imperfection. Thermal buckling, post-buckling, and random responses to thermal and inertial (base acceleration) loads are studied. The results demonstrate the effectiveness of SMAHC structures in controlling static and dynamic responses by adaptive stiffening. Excellent agreement is achieved between the predicted and measured results of the static and dynamic thermomechanical response, thereby providing quantitative validation of the numerical tool.

Davis, Brian↗

NASA NTRS2005

Investigation of the Thermomechanical Response of Shape Memory Alloy Hybrid Composite Beams

Previous work at NASA Langley Research Center (LaRC) involved fabrication and testing of composite beams with embedded, pre-strained shape memory alloy (SMA) ribbons. That study also provided comparison of experimental results with numerical predictions from a research code making use of a new thermoelastic model for shape memory alloy hybrid composite (SMAHC) structures. The previous work showed qualitative validation of the numerical model. However, deficiencies in the experimental-numerical correlation were noted and hypotheses for the discrepancies were given for further investigation. The goal of this work is to refine the experimental measurement and numerical modeling approaches in order to better understand the discrepancies, improve the correlation between prediction and measurement, and provide rigorous quantitative validation of the numerical model. Thermal buckling, post-buckling, and random responses to thermal and inertial (base acceleration) loads are studied. Excellent agreement is achieved between the predicted and measured results, thereby quantitatively validating the numerical tool.

Davis, Brian A.↗

NASA NTRS2006

A Quick Method for Estimating Vehicle Characteristics Appropriate for Continuous Thrust Round Trip Missions Within the Solar System

A simple method of estimating vehicle parameters appropriate for interplanetary travel can provide a useful tool for evaluating the suitability of particular propulsion systems to various space missions. Although detailed mission analyses for interplanetary travel can be quite complex, it is possible to derive hirly simple correlations which will provide reasonable trip time estimates to the planets. In the present work, it is assumed that a constant thrust propulsion system propels a spacecraft on a round trip mission having equidistant outbound and inbound legs in which the spacecraft accelerates during the first portion of each leg of the journey and decelerates during the last portion of each leg of the journey. Comparisons are made with numerical calculations from low thrust trajectory codes to estimate the range of applicability of the simplified correlations.

Emrich, Bill↗

NASA NTRS2008

Projectile and Lab Frame Differential Cross Sections for Electromagnetic Dissociation

Differential cross sections for electromagnetic dissociation in nuclear collisions are calculated for the first time. In order to be useful for three - dimensional transport codes, these cross sections have been calculated in both the projectile and lab frames. The formulas for these cross sections are such that they can be immediately used in space radiation transport codes. Only a limited amount of data exists, but the comparison between theory and experiment is good.

Norbury, John W.↗

NASA NTRS2013

Launch Vehicle Debris Models and Crew Vehicle Ascent Abort Risk

For manned space launch systems, a reliable abort system is required to reduce the risks associated with a launch vehicle failure during ascent. Understanding the risks associated with failure environments can be achieved through the use of physics-based models of these environments. Debris fields due to destruction of the launch vehicle is one such environment. To better analyze the risk posed by debris, a physics-based model for generating launch vehicle debris catalogs has been developed. The model predicts the mass distribution of the debris field based on formulae developed from analysis of explosions. Imparted velocity distributions are computed using a shock-physics code to model the explosions within the launch vehicle. A comparison of the debris catalog with an existing catalog for the Shuttle external tank show good comparison in the debris characteristics and the predicted debris strike probability. The model is used to analyze the effects of number of debris pieces and velocity distributions on the strike probability and risk.

Gee, Ken↗

NASA NTRSDate not supplied

End-To-End Uncertainty Quantification with Analytical Derivatives for Design Under Uncertainty

Uncertainty quantification (UQ) is a rapidly growing and evolving discipline, especially within the aerospace community. Performing analysis with UQ can provide decision makers with a wealth of information about a candidate design. However, the value of UQ is fully realized when the information gained during UQ analysis is leveraged in a feedback loop of a design optimization process, often referred to as design under uncertainty. Although design under uncertainty can be a powerful risk mitigation technique, there are a number of roadblocks that prevent its implementation. Two primary factors are computational costs and added complexity of the analysis. High fidelity simulations on the order tens of uncertain variables quickly become computationally infeasible. Also, implementing UQ into an existing multidisciplinary design and optimization (MDO) process often requires extensive knowledge of the UQ methods and careful treatment of the problem formulation. The objective of this work is to address these two primary roadblocks and enable practitioners to efficiently perform design under uncertainty with limited knowledge of the UQ discipline. Methods outlined in this paper demonstrate MDO incorporating UQ into the design process, leveraging an analytic derivative tool chain through the entire optimization. The proposed approach leverages machine learning techniques to generate a differentiable confidence interval output from polynomial chaos models. This technique, coupled with the incorporation of analytical derivatives through the Polynomial Chaos Expansion (PCE) process, eliminates the need to estimate derivatives which are usually obtained from finite difference, complex step, or similar methods. Developing a differentiable confidence interval allows mixed uncertainty problems (both epistemic and aleatory) to be modeled. Without such modeling, these problems cannot accurately predict objective functions containing statistical quantities such as mean and variance. The addition of analytic derivatives to a polynomial chaos-based UQ method decreases the computational costs of performing design under uncertainty by orders of magnitude in comparison with methods such as complex step. The method and codes developed are modular in nature and are a drop-in solution for design under uncertainty within existing MDO problems. A low-fidelity analytical multidisciplinary optimization under uncertainty for a wing design in OpenMDAO is detailed in this paper. This demonstration case will include both objective functions and constraints which are influenced by uncertain parameters.

Ben D Phillips↗

NASA NTRSDate not supplied

End-To-End Uncertainty Quantification with Analytical Derivatives for Design Under Uncertainty

Uncertainty quantification (UQ) is a rapidly growing and evolving discipline, especially within the aerospace community. Performing analysis with UQ can provide decision makers with a wealth of information about a candidate design. However, the value of UQ is fully realized when the information gained during UQ analysis is leveraged in a feedback loop of a design optimization process, often referred to as design under uncertainty. Although design under uncertainty can be a powerful risk mitigation technique, there are a number of roadblocks that prevent its implementation. Two primary factors are computational costs and added complexity of the analysis. High fidelity simulations on the order tens of uncertain variables quickly become computationally infeasible. Also, implementing UQ into an existing multidisciplinary design and optimization (MDO) process often requires extensive knowledge of the UQ methods and careful treatment of the problem formulation. The objective of this work is to address these two primary roadblocks and enable practitioners to efficiently perform design under uncertainty with limited knowledge of the UQ discipline. Methods outlined in this paper demonstrate MDO incorporating UQ into the design process, leveraging an analytic derivative tool chain through the entire optimization. The proposed approach leverages machine learning techniques to generate a differentiable confidence interval output from polynomial chaos models. This technique, coupled with the incorporation of analytical derivatives through the Polynomial Chaos Expansion (PCE) process, eliminates the need to estimate derivatives which are usually obtained from finite difference, complex step, or similar methods. Developing a differentiable confidence interval allows mixed uncertainty problems (both epistemic and aleatory) to be modeled. Without such modeling, these problems cannot accurately predict objective functions containing statistical quantities such as mean and variance. The addition of analytic derivatives to a polynomial chaos-based UQ method decreases the computational costs of performing design under uncertainty by orders of magnitude in comparison with methods such as complex step. The method and codes developed are modular in nature and are a drop-in solution for design under uncertainty within existing MDO problems. A low-fidelity analytical multidisciplinary optimization under uncertainty for a wing design in OpenMDAO is detailed in this paper. This demonstration case will include both objective functions and constraints which are influenced by uncertain parameters.

Ben Phillips↗

NASA NTRS1991

Comparison of heating rate calculations with experimental data on a modified Shuttle orbiter at Mach 6

Heating-rate calculations from an 'engineering' code and a 'benchmark' Navier-Stokes code are compared with experimental wind tunnel data obtained on a 'modified' 0.0075-scale Shuttle orbiter at a Mach number of 6 and a freestream Reynolds number of 1.9 x 10 to the 6th/m. Comparisons are presented along the windward symmetry plane, in a circumferential direction around the body, and along the wing leading edge at angles of attack of 30 and 40 deg. Both codes provide accurate predictions of the heating rate over most of the body, but because the run time for the engineering code is relatively short, it is an ideal tool for both parametric and design studies.

Hamilton, H. H., II↗

NASA NTRSDate not supplied

Validation of the Hypersolve CFD Solver for Entry Descent and Landing Applications

The functional equivalence of the HyperSolve unstructured edge-based, finite-volume computational fluid dynamics code to the Langley Aerothermodynamic Upwind Relaxation Algorithm multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing community. A suite of cases using a range of thermochemical gas models on relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flowfield quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. The functional equivalence of the HyperSolve unstructured edge-based finite-volume computational fluid dynamics (CFD) code to the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing (EDL) community. A suite of cases using a range of thermochemical gas models on EDL-relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flow field quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. In general, HyperSolve predictions for surface pressure and surface heat flux are in close agreement with those predicted by LAURA.

hypersolve↗

NASA NTRSDate not supplied

Validation of the HyperSolve CFD Solver for Entry Descent and Landing Applications

The functional equivalence of the HyperSolve unstructured edge-based, finite-volume computational fluid dynamics code to the Langley Aerothermodynamic Upwind Relaxation Algorithm multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing community. A suite of cases using a range of thermochemical gas models on relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flowfield quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. The functional equivalence of the HyperSolve unstructured edge-based finite-volume computational fluid dynamics (CFD) code to the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing (EDL) community. A suite of cases using a range of thermochemical gas models on EDL-relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flow field quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. In general, HyperSolve predictions for surface pressure and surface heat flux are in close agreement with those predicted by LAURA.

hypersolve↗

NASA NTRS1992

Fast methods to numerically integrate the Reynolds equation for gas fluid films

The alternating direction implicit (ADI) method is adopted, modified, and applied to the Reynolds equation for thin, gas fluid films. An efficient code is developed to predict both the steady-state and dynamic performance of an aerodynamic journal bearing. An alternative approach is shown for hybrid journal gas bearings by using Liebmann's iterative solution (LIS) for elliptic partial differential equations. The results are compared with known design criteria from experimental data. The developed methods show good accuracy and very short computer running time in comparison with methods based on an inverting of a matrix. The computer codes need a small amount of memory and can be run on either personal computers or on mainframe systems.

Dimofte, Florin↗