SEARCH · Search NASA

Results for “analysis codes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS2004

LASTRAC.3d: Transition Prediction in 3D Boundary Layers

Langley Stability and Transition Analysis Code (LASTRAC) is a general-purpose, physics-based transition prediction code released by NASA for laminar flow control studies and transition research. This paper describes the LASTRAC extension to general three-dimensional (3D) boundary layers such as finite swept wings, cones, or bodies at an angle of attack. The stability problem is formulated by using a body-fitted nonorthogonal curvilinear coordinate system constructed on the body surface. The nonorthogonal coordinate system offers a variety of marching paths and spanwise waveforms. In the extreme case of an infinite swept wing boundary layer, marching with a nonorthogonal coordinate produces identical solutions to those obtained with an orthogonal coordinate system using the earlier release of LASTRAC. Several methods to formulate the 3D parabolized stability equations (PSE) are discussed. A surface-marching procedure akin to that for 3D boundary layer equations may be used to solve the 3D parabolized disturbance equations. On the other hand, the local line-marching PSE method, formulated as an easy extension from its 2D counterpart and capable of handling the spanwise mean flow and disturbance variation, offers an alternative. A linear stability theory or parabolized stability equations based N-factor analysis carried out along the streamline direction with a fixed wavelength and downstream-varying spanwise direction constitutes an efficient engineering approach to study instability wave evolution in a 3D boundary layer. The surface-marching PSE method enables a consistent treatment of the disturbance evolution along both streamwise and spanwise directions but requires more stringent initial conditions. Both PSE methods and the traditional LST approach are implemented in the LASTRAC.3d code. Several test cases for tapered or finite swept wings and cones at an angle of attack are discussed.

Chang, Chau-Lyan↗

NASA NTRS1996

Gear crack propagation investigations

Analytical and experimental studies were performed to investigate the effect of gear rim thickness on crack propagation life. The FRANC (FRacture ANalysis Code) computer program was used to simulate crack propagation. The FRANC program used principles of linear elastic fracture mechanics, finite element modeling, and a unique re-meshing scheme to determine crack tip stress distributions, estimate stress intensity factors, and model crack propagation. Various fatigue crack growth models were used to estimate crack propagation life based on the calculated stress intensity factors. Experimental tests were performed in a gear fatigue rig to validate predicted crack propagation results. Test gears were installed with special crack propagation gages in the tooth fillet region to measure bending fatigue crack growth. Good correlation between predicted and measured crack growth was achieved when the fatigue crack closure concept was introduced into the analysis. As the gear rim thickness decreased, the compressive cyclic stress in the gear tooth fillet region increased. This retarded crack growth and increased the number of crack propagation cycles to failure.

Lewicki, David G.↗

NASA NTRS2019

Bumper: A Tool for Analyzing Spacecraft Micrometeoroid and Orbital Debris Risk

“Bumper” is NASA’s computer program for analyzing spacecraft micrometeoroid and orbital debris (MMOD) risk. Bumper was developed in the late-1980s and has been continuously used and maintained since. The user base has grown from a few government entities to now include numerous commercial entities as well. The NASA Johnson Space Center (JSC) Hypervelocity Impact Technology (HVIT) Team is responsible for all aspects of the Bumper software. Bumper has been used to characterize MMOD risk on hundreds of spacecraft. All of the International Space Station (ISS) modules, visiting vehicles and numerous external components and systems have been analyzed. Bumper was used to analyze each of the Space Shuttle missions since STS-50. The Orion Multi-Purpose Crew Vehicle (MPCV) MMOD shielding is being developed using Bumper as well. Bumper has also been used on numerous telescopes (Hubble, James Webb, and Fermi Gamma-ray Space Telescopes), scientific probes (Stardust, New Horizons, Parker Solar Probe), and Earth observation satellites (Landsat, Joint Polar Satellite System). Bumper is also being used to analyze the micrometeoroid risk and support design of the Deep Space Gateway (DSG) and Mars Sample Return (MSR) missions. The HVIT Bumper Configuration Control Board (CCB) ensures that all changes to the code are approved, reviewed, and documented. Most of the changes are made to add new MMOD damage “ballistic limit equations” (BLEs). BLEs are typically added in response to completion of a hypervelocity impact (HVI) test series and development of an associated BLE. Other less frequent changes include updates of the debris or meteoroid environment models, feature enhancements, and feature retirement. Some BLEs are commercially sensitive and/or proprietary, so the CCB also manages code user-version control and software distribution. The current version – “Bumper 3” – is a FORTRAN executable that utilizes a 64-bit architecture. Bumper 3 has numerous features that make it a powerful tool for analyzing spacecraft MMOD risk. Bumper uses the latest orbital debris and micrometeoroid environment models. Bumper also easily processes large spacecraft geometry models, recognizes hidden surfaces, permits BLE assignment by name or number, and conducts quality checks of the spacecraft geometry model. Bumper 3 can also be used to estimate the effects of particle penetration through thin, high-standoff distance hardware components such as solar arrays and radiators. This is done using a special HVIT-developed technique know as the “3-Part Analysis.” The paper introduces the Bumper 3 MMOD risk analysis code and provides an example MMOD risk assessment showing Bumper’s role in the overall MMOD protection design process.

Lear, Dana M.↗

NASA NTRS2020

Tuning Neural Network Models for Improved Prediction of Boundary Layer Transition

Boundary layer transition can strongly impact flight vehicle performance as it influences surface skin friction and aerodynamic heating, making accurate transition prediction a key to designing next generation aircraft. Artificial neural networks (ANNs) have shown promise toward predicting laminar-turbulent transition based on linear stability correlations. The computational efficiency of ANNs and the substantially reduced user involvement in relation to direct computations based on the linear stability theory (LST) makes them an attractive methodology for integrating the LST based correlations in computational fluid dynamics codes. Tollmien-Schlichting (TS) waves correspond to the dominant transition mechanism in 2D or weakly 3D subsonic boundary layers, such as those encountered in general aviation applications. Improvements to neural network model accuracy in predicting the amplification rates of TS instability waves have been investigated by leveraging recent machine learning developments in conjunction with surrogate optimization techniques and via suitable augmentation of the data used to train the networks. The optimized models trained on the modified dataset reduced the average transition location errors on different airfoils at several flow conditions by 51% of the original manually-tuned network’s errors on the same flow cases. The actual transition locations were derived from the Langley Stability and Transition Analysis Code (LASTRAC).

Machine Learning↗

DOE OSTI2017

Stabilizing effect of resistivity towards ELM-free H-mode discharge in lithium-conditioned NSTX

Linear stability analysis of the national spherical torus experiment (NSTX) Li-conditioned ELM-free H-mode equilibria is carried out in the context of the extended magneto-hydrodynamic (MHD) model in NIMROD. The purpose is to investigate the physical cause behind edge localized mode (ELM) suppression in experiment after the Li-coating of the divertor and the first wall of the NSTX tokamak. Besides ideal MHD modeling, including finite-Larmor radius effect and two-fluid Hall and electron diamagnetic drift contributions, a non-ideal resistivity model is employed, taking into account the increase of Z eff after Li-conditioning in ELM-free H-mode. Unlike an earlier conclusion from an eigenvalue code analysis of these equilibria, NIMROD results find that after reduced recycling from divertor plates, profile modification is necessary but insufficient to explain the mechanism behind complete ELMs suppression in ideal two-fluid MHD. After considering the higher plasma resistivity due to higher Z eff , the complete stabilization could be explained. A thorough analysis of both pre-lithium ELMy and with-lithium ELM-free cases using ideal and non-ideal MHD models is presented, after accurately including a vacuum-like cold halo region in NIMROD to investigate ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Drift-flux Analysis in Horizontal Two-phase Flow

While the most detailed approach to model two-phase flows is the two-fluid model, the drift-flux model may also provide effective methods for estimating one-dimensional two-phase flow parameters. For example, in the current nuclear system analysis codes, the interfacial drag force is a function of the void fraction, drift velocity, and area-averaged relative velocity, which are estimated based on drift-flux approach. While correlations for these parameters exist for vertical two-phase flows, correlations for horizontal two-phase flows have not been developed. (authors)

42 ENGINEERING↗

DOE OSTI2017

Stabilizing effect of resistivity towards ELM-free H-mode discharge in lithium-conditioned NSTX

Linear stability analysis of the national spherical torus experiment (NSTX) Li-conditioned ELM-free H-mode equilibria is carried out in the context of the extended magneto-hydrodynamic (MHD) model in NIMROD. Our purpose is to investigate the physical cause behind edge localized mode (ELM) suppression in experiment after the Li-coating of the divertor and the first wall of the NSTX tokamak. Besides ideal MHD modeling, including finite-Larmor radius effect and two-fluid Hall and electron diamagnetic drift contributions, a non-ideal resistivity model is employed, taking into account the increase of Z eff after Li-conditioning in ELM-free H-mode. And unlike an earlier conclusion from an eigenvalue code analysis of these equilibria, NIMROD results find that after reduced recycling from divertor plates, profile modification is necessary but insufficient to explain the mechanism behind complete ELMs suppression in ideal two-fluid MHD. After considering the higher plasma resistivity due to higher Z eff, the complete stabilization could be explained. Furthermore, a thorough analysis of both pre-lithium ELMy and with-lithium ELM-free cases using ideal and non-ideal MHD models is presented, after accurately including a vacuum-like cold halo region in NIMROD to investigate ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2023

CFD Integrated Transition Modeling for High-Speed Flows via Coupled OVERFLOW-LASTRAC Analysis

This work details ongoing efforts at the NASA Langley Research Center to develop and validate a general-use CFD tool that includes built-in predictions of boundary-layer transition in high-speed flows. Existing tools for the coupling of the NASA OVERFLOW structured overset RANS solver and LASTRAC stability analysis code have been extended to capture boundary-layer transition in high-speed flows driven by either Mack’s first-mode or secondmode instabilities and/or crossflow instabilities. The efficacy of this coupled approach has been demonstrated by examining a variety of supersonic and hypersonic test cases for which experimental validation data is available. The intermittency prescription parameters inherited from prior low-speed applications are found to be suitable for some, but not all high-speed flow scenarios. The method is found to be robust for straight cone configurations and improvements necessary to accurately capture transition on geometries with spatially complex amplification factor envelopes with intermediate regions of slow N-factor variation are examined. The automated, coupled analysis is demonstrated for a 3d supersonic test case and found to perform well within the limits of the linear stability analysis on which it relies. Finally, a preliminary investigation of the method’s robustness to non-ideal CFD meshes is conducted.

High-speed flow↗

DOE OSTI2017

DYNAMIC PRA OF A MULTI-UNIT PLANT

Dynamic Probabilistic Risk Analysis (PRA) methods couple stochastic methods (e.g., RAVEN, ADAPT, ADS, MCDET) with safety analysis codes (RELAP5-3D, MELCOR, MAAP) to determine risk associated to complex systems such as nuclear plants. Compared to classical PRA methods, which are based on static logic structures (e.g., Event-Trees, Fault-Trees), they can evaluate with higher resolution the safety impact of timing and sequencing of events on the accident progression. Recently, special attention has been given to nuclear plants which consist of multiple units and, in particular, on the safety impact of system dependencies, shared systems and common resources on core damage frequencies. In the literature, while classical PRA methods have been employed to model multi-unit plants, Dynamic PRA methods have never been applied to analyze a full multi-unit model. This paper presents a PRA analysis of a multi-unit plant using Dynamic PRA methods. We employ RAVEN as stochastic method coupled with RELAP5-3D. The plant under consideration consists of the three units and their associated spent fuel pools (SFPs). The studied initiating event is a seismic induced station blackout event. We will describe in detail how the multi-unit plant has been constructed and, in particular, how unit dependencies and shared resources are modeled.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1998

Minimum Climb to Cruise Noise Trajectories Modeled for the High Speed Civil Transport

The proposed U.S. High Speed Civil Transport (HSCT) will revolutionize commercial air travel by providing economical supersonic passenger service to destinations worldwide. Unlike the high-bypass turbofan engines that propel today's subsonic airliners, HSCT engines will have much higher jet exhaust speeds. Jet noise, caused by the turbulent mixing of high-speed exhaust with the surrounding air, poses a significant challenge for HSCT engine designers. To resolve this challenge, engineers have designed advanced mixer rejector nozzles that reduce HSCT jet noise to airport noise certification levels by entraining and mixing large quantities of ambient air with the engines' jet streams. Although this works well during the first several minutes of flight, far away from the airport, as the HSCT gains speed and climbs, poor ejector inlet recovery and ejector ram drag contribute to poor thrust, making it advantageous to turn off the ejector. Doing so prematurely, however, can cause unacceptable noise levels to propagate to the ground, even when the aircraft is many miles from the airport. This situation lends itself ideally to optimization, where the aircraft trajectory, throttle setting, and ejector setting can be varied (subject to practical aircraft constraints) to minimize the noise propagated to the ground. A method was developed at the NASA Lewis Research Center that employs a variation of the classic energy state approximation: a trajectory analysis technique historically used to minimize climb time or fuel burned in many aircraft problems. To minimize the noise on the ground at any given throttle setting, high aircraft altitudes are desirable; but the HSCT may either climb quickly to high altitudes using a high, noisy throttle setting or climb more slowly at a lower, quieter throttle setting. An optimizer has been programmed into NASA's existing aircraft and noise analysis codes to balance these options by dynamically choosing the best altitude-velocity path and throttle setting history. The noise level standard, or metric, used in the optimizer should be one that accurately reflects the subjective annoyance levels of ground-based observers under the flight path. A variety of noise metrics are available, many of which are practical for airport-vicinity noise certification. Unlike airport noise, however, the HSCT's climb noise will be characterized by relatively low noise levels, long durations, and low-frequency spectra. The noise metrics used in these calculations are based on the recommendations of researchers at the NASA Langley Research Center, who have correlated the flyover noise annoyance levels of actual laboratory subjects with a variety of measurements. Analysis of data from this optimizer has shown that significant reductions in noise may be obtained with trajectory optimization. And since throttling operations are performed in the subsonic portion of the climb path (where thrust is plentiful), only small penalties in HSCT range or fuel performance occur.

Berton, Jeffrey J.↗

NASA NTRS1989

An experimental investigation of multi-element airfoil ice accretion and resulting performance degradation

An investigation of the ice accretion pattern and performance characteristics of a multi-element airfoil was undertaken in the NASA Lewis 6- by 9-Foot Icing Research Tunnel. Several configurations of main airfoil, slat, and flaps were employed to examine the effects of ice accretion and provide further experimental information for code validation purposes. The text matrix consisted of glaze, rime, and mixed icing conditions. Airflow and icing cloud conditions were set to correspond to those typical of the operating environment anticipated tor a commercial transport vehicle. Results obtained included ice profile tracings, photographs of the ice accretions, and force balance measurements obtained both during the accretion process and in a post-accretion evaluation over a range of angles of attack. The tracings and photographs indicated significant accretions on the slat leading edge, in gaps between slat or flaps and the main wing, on the flap leading-edge surfaces, and on flap lower surfaces. Force measurments indicate the possibility of severe performance degradation, especially near C sub Lmax, for both light and heavy ice accretion and performance analysis codes presently in use. The LEWICE code was used to evaluate the ice accretion shape developed during one of the rime ice tests. The actual ice shape was then evaluated, using a Navier-Strokes code, for changes in performance characteristics. These predicted results were compared to the measured results and indicate very good agreement.

Potapczuk, Mark G.↗

NASA NTRS1989

An experimental investigation of multi-element airfoil ice accretion and resulting performance degradation

An investigation of the ice accretion pattern and performance characteristics of a multi-element airfoil was undertaken in the NASA Lewis 6- by 9-Foot Icing Research Tunnel. Several configurations of main airfoil, slat, and flaps were employed to examine the effects of ice accretion and provide further experimental information for code validation purposes. The text matrix consisted of glaze, rime, and mixed icing conditions. Airflow and icing cloud conditions were set to correspond to those typical of the operating environment anticipated for a commercial transport vehicle. Results obtained included ice profile tracings, photographs of the ice accretions, and force balance measurements obtained both during the accretion process and in a past-accretion evaluation over a range of angles of attack. The tracings and photographs indicated significant accretions on the slat leading edge, in gaps between slat or flaps and the main wing, on the flap leading-edge surfaces, and on flap lower surfaces. Force measurements indicate the possibility of severe performance degradation, especially near C sub Lmax, for both light and heavy ice accretion and performance analysis codes presently in use. The LEWICE code was used to evaluate the ice accretion shape developed during one of the rime ice tests. The actual ice shape was then evaluated, using a Navier-Strokes code, for changes in performance characteristics. These predicted results were compared to the measured results and indicate very good agreement.

Potapczuk, Mark G.↗

NASA NTRSDate not supplied

Modeling Boundary-Layer Transition in Subsonic Flow over a Swept Wing

Predicting the onset of boundary-layer transition is often more accurate using physics-based models that directly compute disturbance growth rather than phenomenological models often implemented into industrial CFD codes. The aim of this ongoing study is to calibrate linear, physics-based computations of transition in subsonic flows over swept wings against a large set of experimental data. Advancing the calibration of linear models of transition contributes to the CFD-Vision-2030 goal of automated boundary-layer transition prediction. This progress report uses the dual N-factor method to model transition over the swept NACA 64-2-015A wing. The flow conditions match selected test conditions from an extensive experimental dataset acquired from the NASA Ames 12-ft Pressure Tunnel. The OVERFLOW 2.4b flow solver is used to obtain laminar basic states based on an infinite-span assumption. Stability analyses are performed on 365 distinct configurations with linear stability theory (LST) and parabolized stability equations (PSE) from the Langley Stability and Transition Analysis Codes (LASTRAC), modeling the growth of Tollmien-Schlichting (TS) and stationary crossflow (SCF) disturbances. From a total of 67 data points for unswept, i.e., TS-dominant configurations, the critical N-factor based on PSE is found to be N_TS = 9. The SCF critical N-factor is found to be near 8 for the highly swept, SCF-dominant configurations. Dual N-factor curves for both LST and PSE computations demonstrate a high level of interaction between TS and SCF. It may be worthwhile to investigate an alternate metric to visualize maximal SCF amplification upstream of the transition location to account for the growth of SCF modes near the leading edge, which is not considered in the conventional applications of the dual N-factor criterion.

boundary-layer transition↗

NASA NTRSDate not supplied

Modeling Boundary-Layer Transition in Subsonic Flow over a Swept Wing

computational modeling↗

NASA NTRS1995

Wave Journal Bearing. Part 2: Experimental Pressure Measurements and Fractional Frequency Whirl Threshold for Wave and Plain Journal Bearings

A new hydrodynamic bearing concept, the wave journal bearing, is being developed because it has better stability characteristics than plain journal bearings while maintaining similar load capacity. An analysis code to predict the steady state and dynamic performance of the wave journal bearing is also part of the development. To verify numerical predictions and contrast the wave journal bearing's stability characteristics to a plain journal bearing, tests were conducted at NASA Lewis Research Center using an air bearing test rig. Bearing film pressures were measured at 16 ports located around the bearing circumference at the middle of the bearing length. The pressure measurements for both a plain journal bearing and a wave journal bearing compared favorably with numerical predictions. Both bearings were tested with no radial load to determine the speed threshold for self-excited fractional frequency whirl. The plain journal bearing started to whirl immediately upon shaft start-up. The wave journal did not incur self-excited whirl until 800 to 900 rpm as predicted by the analysis. Furthermore, the wave bearing's geometry limited the whirl orbit to less than the bearing's clearance. In contrast, the plain journal bearing did not limit the whirl orbit, causing it to rub.

Walker, James F.↗

DOE OSTI2019

Measurement Error Variance Estimation in Gamma-Spectroscopy Data using FRAM

The FRAM gamma-ray isotopic analysis code has recently been applied to four Pu data sets and four U data sets, each with relatively well-characterized standard items, in order to estimate the variances of random errors and of item-specific systematic errors. Assay results are available for Pu238, Pu239, Pu240, Pu241, Pu242 (operator entered, not independently measured), Am241, Specific Power, and Pu240eff. For the four U data sets, results are available for the measurands U234, U235, U236 (operator entered, not independently measured), and U238. For example, in Pu data set 1, there are 33 standard reference Pu items, which have well-known nominal/true values (except for a few items). FRAM was used to make approximately 15 repeat measurements for each of the 33 items, so empirical standard deviations can be calculated. In addition, FRAM internally estimates its random error variance using propagation of variance. This paper applies analysis of variance (ANOVA) to separately estimate the repeatability standard deviation and the standard deviation of the item-specific biases δ s in each of the 8 data sets.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2000

Modeling Woven Polymer Matrix Composites with MAC/GMC

NASA's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) is used to predict the elastic properties of plain weave polymer matrix composites (PMCs). The traditional one step three-dimensional homogertization procedure that has been used in conjunction with MAC/GMC for modeling woven composites in the past is inaccurate due to the lack of shear coupling inherent to the model. However, by performing a two step homogenization procedure in which the woven composite repeating unit cell is homogenized independently in the through-thickness direction prior to homogenization in the plane of the weave, MAC/GMC can now accurately model woven PMCs. This two step procedure is outlined and implemented, and predictions are compared with results from the traditional one step approach and other models and experiments from the literature. Full coupling of this two step technique with MAC/ GMC will result in a widely applicable, efficient, and accurate tool for the design and analysis of woven composite materials and structures.

Bednarcyk, Brett A.↗

DOE OSTI2022

HYSPLIT/MACCS Atmospheric Dispersion Model Technical Documentation and Benchmark Analysis

The nuclear accident consequence analysis code MACCS has traditionally modeled dispersion during downwind transport using a Gaussian plume segment model. MACCS is designed to estimate consequence measures such as air concentrations and ground depositions, radiological doses, and health and economic impacts on a statistical basis over the course of a year to produce annualaveraged output measures. The objective of this work is to supplement the Gaussian atmospheric transport and diffusion (ATD) model currently in MACCS with a new option using the HYSPLIT model. HYSPLIT/MACCS coupling has been implemented, with HYSPLIT as an alternative ATD option. The subsequent calculations in MACCS use the HYSPLIT-generated air concentration, and ground deposition values to calculate the same range of output quantities (dose, health effects, risks, etc.) that can be generated when using the MACCS Gaussian ATD model. Based on the results from the verification test cases, the implementation of the HYSPLIT/MACCS coupling is confirmed. This report contains technical details of the HYSPLIT/MACCS coupling and presents a benchmark analysis using the HYSPLIT/MACCS coupling system. The benchmark analysis, which involves running specific scenarios and sensitivity studies designed to examine how the results generated by the traditional MACCS Gaussian plume segment model compare to the new, higher fidelity HYSPLIT/MACCS modeling option, demonstrates the modeling results that can be obtained by using this new option. The comparisons provided herein can also help decision-makers evaluate the potential benefit of using results based on higher fidelity modeling with the additional computational burden needed to perform the calculations. Three sensitivity studies to investigate the potential impact of alternative modeling options, regarding 1) input meteorological data set, 2) method to estimate stability class, and 3) plume dispersion model for larger distances, on consequence results were also performed. The results of these analyses are provided and discussed in this report.

97 MATHEMATICS AND COMPUTING↗