Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer code”

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 613 records · Page 34

Computational design of an experimental laser-powered thruster

An extensive numerical experiment, using the developed computer code, was conducted to design an optimized laser-sustained hydrogen plasma thruster. The plasma was sustained using a 30 kW CO2 laser beam operated at 10.6 micrometers focused inside the thruster. The adopted physical model considers two-dimensional compressible Navier-Stokes equations coupled with the laser power absorption process, geometric ray tracing for the laser beam, and the thermodynamically equilibrium (LTE) assumption for the plasma thermophysical and optical properties. A pressure based Navier-Stokes solver using body-fitted coordinate was used to calculate the laser-supported rocket flow which consists of both recirculating and transonic flow regions. The computer code was used to study the behavior of laser-sustained plasmas within a pipe over a wide range of forced convection and optical arrangements before it was applied to the thruster design, and these theoretical calculations agree well with existing experimental results. Several different throat size thrusters operated at 150 and 300 kPa chamber pressure were evaluated in the numerical experiment. It is found that the thruster performance (vacuum specific impulse) is highly dependent on the operating conditions, and that an adequately designed laser-supported thruster can have a specific impulse around 1500 sec. The heat loading on the wall of the calculated thrusters were also estimated, and it is comparable to heat loading on the conventional chemical rocket. It was also found that the specific impulse of the calculated thrusters can be reduced by 200 secs due to the finite chemical reaction rate.

Jeng, San-Mou↗

Steady and unsteady transonic flow

An investigation of the transonic flow over a circular arc airfoil was conducted to obtain basic information for turbulence modeling of shock-induced separated flows and to verify numerical computer codes which are being developed to simulate such flows. The investigation included the employment of a laser velocimeter to obtain data concerning the mean velocity, the shear stress, and the turbulent kinetic energy profiles in the flowfield downstream of the airfoil midchord where the flow features are more complex. Depending on the free-stream Mach number, the flowfield developed was either steady with shock-wave-induced separation extending from the foot of the shock wave to beyond the trailing edge or unsteady with a periodic motion also undergoing shock-induced separation. The experimental data were compared with the results of numerical simulations in which a computer code was employed that solved the time-dependent Reynolds' averaged compressible Navier-Stokes equations.

Seegmiller, H. L.↗

Investigation of Vortex Flaps and Other Flow Control Devices on Generic High-Speed Civil Transport Planforms

A parametric cranked delta planform study has been conducted in the Langley 14- by 22-Foot Subsonic Tunnel with the following objectives: (1) to evaluate the vortex flap design methodology for cranked delta wings, (2) to determine the influence of leading-edge sweep and the outboard wing on vortex flap effectiveness, (3) to evaluate novel flow control concepts, and (4) to validate unstructured grid Euler computer code predictions with modeled vortex and trailing-edge flaps. Two families of cranked delta planforms were investigated. One family had constant aspect ratio, while the other had a constant nondimensional semispan location of the leading-edge break. The inboard leading-edge sweep of the planforms was varied between 68 deg., 71 deg., and 74 deg., while outboard leading-edge sweep was varied between 48 deg. and 61 deg. Vortex flaps for the different planforms were designed by an analytical vortex flap design method. The results indicate that the effectiveness of the vortex flaps was only slightly influenced by the variations in the parametric planforms. The unstructured grid Euler computer code was successfully used to model the configurations with vortex flaps. The vortex trap concept was successfully demonstrated.

Kjerstad, Kevin J.↗

Unit cell geometry of multiaxial preforms for structural composites

The objective of this study is to investigate the yarn geometry of multiaxial preforms. The importance of multiaxial preforms for structural composites is well recognized by the industry but, to exploit their full potential, engineering design rules must be established. This study is a step in that direction. In this work the preform geometry for knitted and braided preforms was studied by making a range of well designed samples and studying them by photo microscopy. The structural geometry of the preforms is related to the processing parameters. Based on solid modeling and B-spline methodology a software package is developed. This computer code enables real time structural representations of complex fiber architecture based on the rule of preform manufacturing. The code has the capability of zooming and section plotting. These capabilities provide a powerful means to study the effect of processing variables on the preform geometry. the code also can be extended to an auto mesh generator for downstream structural analysis using finite element method. This report is organized into six sections. In the first section the scope and background of this work is elaborated. In section two the unit cell geometries of braided and multi-axial warp knitted preforms is discussed. The theoretical frame work of yarn path modeling and solid modeling is presented in section three. The thin section microscopy carried out to observe the structural geometry of the preforms is the subject in section four. The structural geometry is related to the processing parameters in section five. Section six documents the implementation of the modeling techniques into the computer code MP-CAD. A user manual for the software is also presented here. The source codes and published papers are listed in the Appendices.

Ko, Frank↗

ICAN/DAMP-integrated composite analyzer with damping analysis capabilities: User's manual

This manual describes the use of the computer code ICAN/DAMP (Integrated Composite Analyzer with Damping Analysis Capabilities) for the prediction of damping in polymer-matrix composites. The code is written in FORTRAN 77 and is a version of the ICAN (Integrated Composite ANalyzer) computer program. The code incorporates a new module for synthesizing the material damping from micromechanics to laminate level. Explicit micromechanics equations based on hysteretic damping are programmed relating the on-axis damping capacities to the fiber and matrix properties and fiber volume ratio. The damping capacities of unidirectional composites subjected to off-axis loading are synthesized from on-axis damping values. The hygrothermal effect on the damping performance of unidirectional composites caused by temperature and moisture variation is modeled along with the damping contributions from interfacial friction between broken fibers and matrix. The temperature rise is continuously vibrating composite plies and composite laminates is also estimated. The ICAN/DAMP user's manual provides descriptions of the damping analysis module's functions, structure, input requirements, output interpretation, and execution requirements. It only addresses the changes required to conduct the damping analysis and is used in conjunction with the 'Second Generation Integrated Composite Analyzer (ICAN) Computer Code' user's manual (NASA TP-3290).

Saravanos, Dimitrious A.↗

Benchmark Specification for FFTF LOFWOS Test #13

The Fast Flux Test Facility (FFTF) at the Hanford site in Washington was designed by the Westinghouse Electric Corporation for the U.S. Department of Energy. FFTF was a 400 MW thermal, oxide-fueled, liquid sodium cooled test reactor, built to assist development and testing of advanced fuels and materials for fast breeder reactors. After reaching criticality in 1980, FFTF operated until 1992, providing the U.S. Department of Energy (DOE) with the means to test fuels, materials, and other components in a fast neutron flux environment. In July 1986, a series of unprotected transients (with the plant protection system intentionally disabled) were performed in FFTF as part of the passive safety demonstration program. Among these were thirteen loss of flow without scram (LOFWOS) tests. The goals of this program included confirming the liquid metal reactor safety margins, providing data for computer code validation, and demonstrating the inherent and passive safety benefits of specific design features. The test defined in this benchmark is LOFWOS Test #13, which was initiated at 50% power and 100% flow with the pump pony motors turned off. This benchmark specification is intended to support collaborative efforts within international partnerships on the validation of simulation tools and models in the area of Sodium-cooled Fast Reactor (SFR) safety. Validated tools and models are needed to evaluate SFR inherent safety characteristics and assess the impact of passive design features in response to accident initiators. Comparisons with experimental data and the results of safety analyses from other groups create unique opportunities to improve predictive capabilities of computational codes and methods for SFR modeling and simulation. The conditions of the LOFWOS test along with the feedback from FFTF’s limited free bow core restraint system and the novel passive reactivity control Gas-Expansion Modules (GEMs) pose a very challenging and uniquely valuable benchmark exercise.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Numerical Method of Calculating Propeller Noise Including Acoustic Nonlinear Effects

Using the transonic flow fields(s) generated by the NASPROP-E computer code for an eight blade SR3-series propeller, a theoretical method is investigated to calculate the total noise values and frequency content in the acoustic near and far field without using the Ffowcs Williams - Hawkings equation. The flow field is numerically generated using an implicit three dimensional Euler equation solver in weak conservation law form. Numerical damping is required by the differencing method for stability in three dimensions, and the influence of the damping on the calculated acoustic values is investigated. The acoustic near field is solved by integrating with respect to time the pressure oscillations induced at a stationary observer location. The acoustic far field is calculated from the near field primitive variables as generated by NASPROP-E computer code using a method involving a perturbation velocity potential as suggested by Hawkings in the calculation of the acoustic pressure time-history at a specified far field observed location. the methodologies described are valid for calculating total noise levels and are applicable to any propeller geometry for which a flow field solution is available.

Korkan, K. D.↗

Investigations of radiation-induced and carrier-enhanced conductivity

A steady-state carrier computer code, PECK (Parker Enhanced Carrier Kinetics), that predicts the radiation-induced conductivity (RIC) produced in a dielectric by an electron beam was developed. The model, which assumes instantly-trapped holes, was then applied to experimental measurements on thin Kapton samples penetrated by an electron beam. Measurements at high bias were matched in the model by an appropriate choice for the trap-modulated electron mobility. A fractional split between front and rear currents measured at zone bias is explained on the basis of beam-scattering. The effects of carrier-enhanced conductivity (CEC) on data obtained for thick, free-surface Kapton samples is described by using an analytical model that incorporates field injection of carriers from the RIC region. The computer code, LWPCHARGE, modified for carrier transport, is also used to predict partial penetration effects associated with CEC in the unirradiated region. Experimental currents and surface voltages, when incorporated in the appropriate models, provide a value for the trap modulated mobility that is in essential agreement with the RIC results.

Meulenberg, A., Jr.↗

Development of phenomena identification and ranking table for Westinghouse lead fast reactor’s safety

The Westinghouse Lead-cooled Fast Reactor (LFR) is a medium-size, passively safe, economic, Gen-IV nuclear reactor. An important effort within the Westinghouse LFR program is the development of the safety analysis methodology, which comprises computer code development, model development, and experimental testing. A key initial task associated with the development of the safety analysis methodology is the identification of processes and phenomena that affect the plant's capability to meet selected safety performance indicators. This is accomplished through the development of a Phenomena Identification and Ranking Table (PIRT) for selected accident scenarios which, for this specific PIRT effort, included selected postulated design basis accidents and hypothetical beyond design basis accidents in LFRs. This paper describes the role of PIRT in the development of the Westinghouse LFR safety analysis methodology and the process used in the PIRT development. Specifically, the Westinghouse LFR PIRT assessed importance of pertinent phenomena and identified gaps in their knowledge-base by evaluating current modeling capabilities and data available for validation. The ultimate goal was to provide guidance on computer code development and validation efforts and to prioritize testing to support LFR design and licensing. The key phenomena and processes that are deemed highly important for the safety performance indicators, but for which the state of knowledge is low, are presented. The testing program and analyses development are planned to address significant phenomena in the PIRT.

Lead fast reactor↗

Transferring ecosystem simulation codes to supercomputers

Many ecosystem simulation computer codes have been developed in the last twenty-five years. This development took place initially on main-frame computers, then mini-computers, and more recently, on micro-computers and workstations. Supercomputing platforms (both parallel and distributed systems) have been largely unused, however, because of the perceived difficulty in accessing and using the machines. Also, significant differences in the system architectures of sequential, scalar computers and parallel and/or vector supercomputers must be considered. We have transferred a grassland simulation model (developed on a VAX) to a Cray Y-MP/C90. We describe porting the model to the Cray and the changes we made to exploit the parallelism in the application and improve code execution. The Cray executed the model 30 times faster than the VAX and 10 times faster than a Unix workstation. We achieved an additional speedup of 30 percent by using the compiler's vectoring and 'in-line' capabilities. The code runs at only about 5 percent of the Cray's peak speed because it ineffectively uses the vector and parallel processing capabilities of the Cray. We expect that by restructuring the code, it could execute an additional six to ten times faster.

Skiles, J. W.↗

Oxidation Behavior and Property Degradation of Nuclear Graphites

During its multidecade operation in the core of nuclear reactors, graphite components are subjected to aggressive and continuous exposure to a high field of ionizing and neutron irradiation, high temperature, and various types of present and postulated chemical attacks. High density, high crystallinity polygranular synthetic graphite is unique among other materials for its extraordinary capacity of resisting and adapting to the aggression inflicted by high temperature, high energy neutron bombardment and ionizing gamma radiation. But, as a carbonaceous material, even though of very high purity, graphite is reactive towards common oxidizing agents: oxygen, carbon dioxide, water. Safe operation of HTGRs relies, among other aspects, on engineered safeguard systems for efficient and continuous protection of graphite components against oxidation. Graphite oxidation behavior was, and continues to be, an important direction of theoretical and experimental research, engineering analyses, models and simulations, and design and safety regulations. The avalanche of publications, reports, experimental data, computer codes, and regulatory documents related to oxidation behavior of nuclear graphite is now accelerating to new levels, prompted by the increased interest for nuclear energy as a clean, carbon-free energy source. Even though public’s perception of nuclear energy advantages may still be influenced by the memories of past accidents of nuclear reactors from generations II and III, the community of informed scientists and engineers, regulators and statemen knows that generation IV of nuclear reactors is designed at very high safety standards, doubled by great advances of scientific knowledge and technological progress. One of routes of these recent advances is directed at better understanding of graphite oxidation behavior, its relationship with graphite manufacturing and microstructural properties, along with the effects of various environmental factors and process variables. Together, the recent progress in manufacturing, properties characterization, and modeling of intricated physical and chemical processes that concur to the oxidation behavior led to development of powerful simulation codes able to analyze various scenarios of normal operation and hypothetical off-normal events, and thus to clearly specify the allowable parameters envelopes for the designers, constructors, and operators of current and future modular HTGRs. This review begins with an introduction on manufacturing methods, structure, and properties of nuclear graphite, including basic requirements that this specialty graphite type must satisfy for nuclear use. It continues with a chapter on environmental effects on nuclear graphite, where emphasis is placed less on irradiation and much more on oxidation phenomena, their safety implications, and the basic traits of chronic and acute oxidation by air (oxygen) and water (humidity, steam). Particular attention is placed on the three graphite grades of interest for this document (IG-110, NBG-18, PCEA). A chapter on properties degradation induced by oxidation follows, with focus on density, dimensional, and mechanical properties changes. The next chapter is intended as a brief review of various approaches used for modeling of graphite oxidation behavior. It summarizes the progress of oxidation models, from the early attempts to complex computational approaches interfaced with specialized computer codes designed for nuclear reactor simulations. Last, a list is presented of knowledge gaps where more research is needed. A short summary concludes the review.

36 MATERIALS SCIENCE↗

Experimental Investigation of Water Droplet Impingement on Airfoils, Finite Wings, and an S-duct Engine Inlet

Validation of trajectory computer codes, for icing analysis, requires experimental water droplet impingement data for a wide range of aircraft geometries as well as flow and icing conditions. This report presents improved experimental and data reduction methods for obtaining water droplet impingement data and provides a comprehensive water droplet impingement database for a range of test geometries including an MS(1)-0317 airfoil, a GLC-305 airfoil, an NACA 65(sub 2)-415 airfoil, a commercial transport tail section, a 36-inch chord natural laminar flow NLF(1)-0414 airfoil, a 48-inch NLF(1)-0414 section with a 25 percent chord simple flap, a state-of-the-art three-element high lift system, a NACA 64A008 finite span swept business jet tail, a full-scale business jet horizontal tail section, a 25 percent-scale business jet empennage, and an S-duct turboprop engine inlet. The experimental results were obtained at the NASA Glenn Icing Research Tunnel (IRT) for spray clouds with median volumetric diameter (MVD) of 11, 11.5, 21, 92, and 94 microns and for a range of angles of attack. The majority of the impingement experiments were conducted at an air speed of 175 mph corresponding to a Reynolds number of approximately 1.6 million per foot. The maximum difference of repeated tests from the average ranged from 0.24 to 12 percent for most of the experimental results presented. This represents a significant improvement in test repeatability compared to previous experimental studies. The increase in test repeatability was attributed to improvements made to the experimental and data reduction methods. Computations performed with the LEWICE-2D and LEWICE-3D computer codes for all test configurations are presented in this report. For the test cases involving median volumetric diameters of 11 and 21 microns, the correlation between the analytical and experimental impingement efficiency distributions was good. For the median volumetric diameters of 92 and 94-micron cases, however, the analysis produced higher impingement efficiencies and larger impingement limits than the experiment. It is speculated that this discrepancy is due to droplet splashing and breakup experienced by large droplets during impingement.

Papadakis, Michael↗

Shielding materials for highly penetrating space radiations

Interplanetary travel involves the transfer from an Earth orbit to a solar orbit. Once outside the Earth's magnetosphere, the major sources of particulate radiation are solar cosmic rays (SCR's) and galactic cosmic rays (GCR's). Intense fluxes of SCR's come from solar flares and consist primarily of protons with energies up to 1 GeV. The GCR consists of a low flux of nuclei with energies up to 10(exp 10) GeV. About 70 percent of the GCR are protons, but a small amount (0.6 percent) are nuclei with atomic numbers greater than 10. High energy charged particles (HZE) interact with matter by transferring energy to atomic electrons in a Coulomb process and by reacting with an atomic nucleus. Energy transferred in the first process increases with the square of the atomic number, so particles with high atomic numbers would be expected to lose large amounts of energy by this process. Nuclear reactions produced by (HZE) particles produce high-energy secondary particles which in turn lose energy to the material. The HZE nuclei are a major concern for radiation protection of humans during interplanetary missions because of the very high specific ionization of both primary and secondary particles. Computer codes have been developed to calculate the deposition of energy by very energetic charged particles in various materials. Calculations show that there is a significant buildup of secondary particles from nuclear fragmentation and Coulomb dissociation processes. A large portion of these particles are neutrons. Since neutrons carry no charge, they only lose energy by collision or reaction with a nucleus. Neutrons with high energies transfer large amounts of energy by inelastic collisions with nuclei. However, as the neutron energy decreases, elastic collisions become much more effective for energy loss. The lighter the nucleus, the greater the fraction of the neutron's kinetic energy that can be lost in an elastic collision. Thus, hydrogen-containing materials such as polymers are most effective in reducing the energy of neutrons. Once neutrons are reduced to very low energies, the probability for undergoing a reaction with a nucleus (the cross section) becomes very high. The product of such a reaction is often radioactive and can involve the release of a significant amount of energy. Thus, it is important to provide protection from low energy neutrons during a long duration space flight. Among the light elements, lithium and boron each have an isotope with a large thermal neutron capture cross section, Li-6 and B-10. However, B-10 is more abundant in the naturally-occurring element than Li-6, has a thermal neutron capture cross section four times that of Li-6, and produces the stable products, He-4 and Li-7 in the interaction while Li-6 produces radioactive tritium (H-3). Thus, boron is the best light-weight material for thermal neutron absorption in spacecraft. The work on this project was focused in two areas: computer design where existing computer codes were used, and in some cases modified, to calculate the propagation and interactions of high energy charged particles through various media, and materials development where boron was incorporated into high performance materials.

Kiefer, Richard L.↗

Coding the Computing Continuum: Fluid Function Execution in Heterogeneous Computing Environments

Advances in network technologies have greatly decreased barriers to accessing physically distributed computers. This newfound accessibility coincides with increasing hardware specialization, creating exciting new opportunities to dispatch workloads to the best resource for a specific purpose, rather than those that are closest or most easily accessible. We present Delta, a service designed to intelligently schedule function-based workloads across a distributed set of heterogeneous computing resources. Delta implements an extensible architecture in which different predictors and scheduling algorithms can be integrated to provide dynamically evolving estimates of function execution times on different resources-estimates that can be used to determine the most appropriate location for execution. We describe predictors for function runtime, data transfer time, and cold-start resource provisioning and configuration delay; dynamic learning methods that update predictor models over time; and scheduling strategies that take into account both function and endpoint information. We show that these methods can halve workload makespan when compared with a strategy that selects the fastest resource, and decrease makespan by a factor of five when compared to a round robin strategy, when deployed on a heterogeneous testbed with resources ranging from a Raspberry Pi to a GPU node in an academic cloud.

Computing continuum↗

Establishing capabilities for quantum computing and simulations for energy applications

Quantum information science (QIS) is creating potential transformative opportunities to exploit intricate quantum mechanical phenomena in new ways for obtaining and processing information to advance many areas of science and engineering. Since the National Quantum Initiative Act was signed into law in 2018, developing QIS capability and competency is one of the most urgent tasks of DOE to make sure the US win the quantum race. The QIS contains four pillars: quantum computing, quantum simulations, quantum sensing, and quantum networking. To apply QIS in energy related applications, the key is to develop the capability of quantum computing & simulation tools. In this project, we propose to develop the capability of quantum computing and simulations at NETL to target fossil energy related problems. We will install a simulator (e. g. IBM qiskit) on the NETL supercomputer to simulate the environments of quantum computer. Based on the current available quantum algorithms for quantum chemistry, we will develop quantum computing codes to perform simulation which will focus on fossil energy sector challenges, including CO2 capture & conversion, sensing, and fuel conversion. We then will seek opportunities to run our codes on real quantum computers (such as IBM-Q, google Sycamore, etc.). Through this project, the new capability of quantum computing and simulations will be established at NETL. In addition, the NETL workforce in this area will be trained ready to conduct more complicated tasks in line with NETL missions to enhance the nation’s energy foundation.

97 MATHEMATICS AND COMPUTING↗