Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer 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.

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

At least 1,549 records · Page 86

The Uses and Abuses of the Acoustic Analogy in Helicopter Rotor Noise Prediction

This paper is theoretical in nature and addresses applications of the acoustic analogy in helicopter rotor noise prediction. It is argued that in many instances the acoustic analogy has not been used with care in rotor noise studies. By this it is meant that approximate or inappropriate formulations have been used. By considering various mechanisms of noise generation, such abuses are identified and the remedy is suggested. The mechanisms discussed are thickness, loading, quadrupole, and blade-vortex interaction noise. The quadrupole term of the Ffowcs Williams-Hawkings equation is written in a new form which separates the contributions of regions of high gradients such as shock surfaces. It is shown by order of magnitude studies that such regions are capable of producing noise with the same directivity as the thickness noise. The inclusion of this part of quadrupole sources in current acoustic codes is quite practical. Some of the difficulties with the use of loading noise formulations of the first author in predictions of blade-vortex interaction noise are discussed. It appears that there is a need for development of new theoretical results based on the acoustic analogy in this area. Because of the impulsive character of the blade surface pressure, a time scale of integration different from that used in loading and thickness computations must he used in a computer code for prediction of blade-vortex interaction noise.

Farassat, F.↗

Experimental Data Sets for CFD Calculations of Flow and Heat Transfer in Deformed Fuel Assemblies

Liquid-metal fast-breeder reactors (LMFBRs) have traditionally used fuel assemblies consisting of a hexagonal array of pins enclosed in a hexagonally shaped duct. Each fuel pin is surrounded by a wire that is wrapped in a helical configuration, which is used to separate the pins sufficiently to allow sodium coolant to flow through the assembly. The helical wire wrap also promotes mixing of the coolant. Experience with these assemblies has shown that they begin to deform with time during normal operation leading to significant changes in core thermal-hydraulics. For example, deformations of assemblies in the Experimental Breeder Reactor-2 and the Fast Flux Test Facility led to pin overheating resulting from pin-to-duct and pin-to-pin contact. This problem is well known and was investigated extensively in the 1970's and early 1980's. This was a period, however, before the development of modern Computational Fluid Dynamics (CFD) computer codes. Therefore, the experiments performed at the time did not record their results in sufficient detail to be useful for compiling validation datasets for CFD codes. Most of the experiments focused only on bulk measurements, such as the pressure drop across the bundle. When the effects of deformation of the fuel assembly were investigated, the studies performed to quantify these effects considered only bare pins without the wire wrap. With the resurgence of interest in LMFBRs in recent years, particularly in designs that leverage long-life performance of the fuel, the need has arisen for data that are suitable for the validation and benchmarking of CFD codes to predict the thermal-hydraulic behavior of flow in wire-wrapped fuel bundles. Such data must include high-resolution, transient, full-field measurements of the fluid velocity and temperature and an accurate record of the test conditions, including a precise description of the geometry. This paper describes a project to produce a validation-quality set of data that can be used to benchmark CFD codes for modeling flow and thermal performance in wire-wrapped fuel assemblies of the kind used by LMFBR designs. In addition, this effort investigates the changes to the flow and other phenomena that occur in the deformed geometries that can form later in a fuel assembly's life. This work is a joint effort of four organizations. The experimental work is being performed by AREVA and Texas A and M University. The computational work, both to support the experiments and to perform blind benchmark calculations, is being done by TerraPower, LLC, and Argonne National Laboratory (ANL). The experiments use a geometry that consists of a complete assembly of 61 wire-wrapped pins, which is surrounded by a hexagonal duct. This configuration is considered the minimum size necessary to adequately capture the important results of interest, such as the thermal gradients across the assembly and the velocity field effects in the assembly resulting from the surrounding duct. These types of effects cannot be captured by a geometry consisting of only a small group of fuel pins. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A method for smoothing multiple yield functions

A method is presented that produces a single C 2 differentiable and convex yield function from a plastic model that contains multiple yield surfaces that are individually C 2 differentiable and convex. The method is suitable for implementation in numerical computer codes. Here, the error incurred through the smoothing procedure is quantified. The method is illustrated through a number of examples. Two simple finite-element models are solved to compare this method with existing approaches.

42 ENGINEERING↗

Simulation of charge exchange plasma propagation near an ion thruster propelled spacecraft

A model describing the charge exchange plasma and its propagation is discussed, along with a computer code based on the model. The geometry of an idealized spacecraft having an ion thruster is outlined, with attention given to the assumptions used in modeling the ion beam. Also presented is the distribution function describing charge exchange production. The barometric equation is used in relating the variation in plasma potential to the variation in plasma density. The numerical methods and approximations employed in the calculations are discussed, and comparisons are made between the computer simulation and experimental data. An analytical solution of a simple configuration is also used in verifying the model.

Robinson, R. S.↗

Design for cyclic loading endurance of composites

The present paper describes the application of the computer code IPACS (Integrated Probabilistic Assessment of Composite Structures) to air craft wing type structures. The code performs a complete probabilistic structural analysis for composites taking into account the uncertainties in geometry, boundary conditions, material properties, laminate lay-ups and loads. Results of the analysis are presented in terms of cumulative distribution functions (CDF) and probability density function (PDF) of life of a wing type composite structure under different hygrothermal environments subjected to random pressure. The sensitivity of fatigue life to a number of critical structural/material variables is also computed from the analysis.

Shiao, Michael C.↗

Influence of wing geometry on leading-edge vortices and vortex-induced aerodynamics at supersonic speeds

An assessment of the influence of wing geometry on wing leading-edge vortex flows at supersonic speeds is discussed as well as the applicability of various aerodynamic codes for predicting these results. A series of delta-wing wind-tunnel models were tested in the NASA Langley Research Center Unitary Plan Wind Tunnel over a Mach number range from 1.6 to 4.6. The data show that wing airfoil has a significant impact on the localized loading on the wing. The experimental data for the flat wings were compared with results from full-potential, Euler, and Parabolized Navier-Stokes (PNS) computer codes. The theoretical evaluation showed that the full-potential analysis predicted accurate results for the attached-flow (alpha = 0 deg) conditions and that the Euler and PNS analyses made reasonable predictions for both attached and separated flow conditions.

Wood, Richard M.↗

Development of a Smart Signal Detection Method for Cyclic Voltammetry via Artificial Neural Intelligence

The electrorefiner (ER) is the heart of pyroprocessing technology which is a high-temperature method for separating uranium from Used Nuclear Fuel (UNF). It is important to improve this technology with respect to nuclear materials detection and accountability. Artificial Neural Intelligence (ANI) is a novel data analysis and simulation method that can be applied to electrochemical data sets. This computational code, which has been performed using the commercial software MATLAB, can be trained to generalize adequate electrical current and potential simulated data sets for the unseen data with a high accuracy of prediction. For this purpose, a massive collection of cyclic voltammetry (CV) data sets by Hoover (2014), for 0.5, 1, 2.5, and 5 wt% of zirconium chloride in LiCl-KCl molten salt with different scan rates at 773 K has been considered. The computer is trained via ANI to predict the unseen data after providing suitable hidden layers and validation numbers. In addition, this work can trace the CV plot for a blind condition by interpolating between two simulated data set. The different hidden layers with various neurons (from 5 to 30) at several validation numbers (from 5 to 30) has been studied and the average percent error between experimental and theoretical data for 0.5 wt% with 200 and 450 mV/s has been calculated. Preliminary results demonstrate that if the number of hidden layers increases from one to three, the average error falls down from 44% to around 8%. The best condition which gives a minimum average percent error has been discovered and the simulated CV graph has been compared with the experimental data. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design for cyclic loading endurance of composites

The application of the computer code IPACS (Integrated Probabilistic Assessment of Composite Structures) to aircraft wing type structures is described. The code performs a complete probabilistic analysis for composites taking into account the uncertainties in geometry, boundary conditions, material properties, laminate lay-ups, and loads. Results of the analysis are presented in terms of cumulative distribution functions (CDF) and probability density function (PDF) of the fatigue life of a wing type composite structure under different hygrothermal environments subjected to the random pressure. The sensitivity of the fatigue life to a number of critical structural/material variables is also computed from the analysis.

Shiao, Michael C.↗

Simulation of Conformal Spiral Slot Antennas on Composite Platforms

During the course of the grant, we wrote and distributed about 12 reports and an equal number of journal papers supported fully or in part by this grant. The list of reports (title & abstract) and papers are given in Appendices A and B. This grant has indeed been instrumental in developing a robust hybrid finite element method for the analysis of complex broadband antennas on doubly curved platforms. Previous to the grant, our capability was limited to simple printed patch antennas on mostly planar platforms. More specifically: (1) mixed element formulations were developed and new edge-based prisms were introduced; (2) these elements were important in permitting flexibility in geometry gridding for most antennas of interest; (3) new perfectly matched absorbers were introduced for mesh truncations associated with highly curved surfaces; (4) fast integral algorithms were introduced for boundary integral truncations reducing CPU time from O(N-2) down to O(N-1.5) or less; (5) frequency extrapolation schemes were developed for efficient broadband performance evaluations. This activity has been successfully continued by NASA researchers; (6) computer codes were developed and extensively tested for several broadband configurations. These include FEMA-CYL, FEMA-PRISM and FEMA-TETRA written by L. Kempel, T. Ozdemir and J. Gong, respectively; (7) a new infinite balun feed was designed nearly constant impedance over the 800-3000 MHz operational band; (8) a complete slot spiral antenna was developed, fabricated and tested at NASA Langley. This new design is a culmination of the projects goals and integrates the computational and experimental efforts. this antenna design resulted in a U.S. patent and was revised three times to achieve the desired bandwidth and gain requirements from 800-3000 MHz.

Volakis, J. L.↗

Particle-In-Cell Computer Simulations of Beam Spill in the Scorpius Injector

Computer simulations of beam transport in the Scorpius injector were performed using the particle-in-cell computer code LSP. The simulations were used to predict spot size growth due to electron stimulated emission of ions from beam spill on the anode stalk and the entrance transport region.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Three-dimensional transonic potential flow about complex 3-dimensional configurations

An analysis has been developed and a computer code written to predict three-dimensional subsonic or transonic potential flow fields about lifting or nonlifting configurations. Possible condfigurations include inlets, nacelles, nacelles with ground planes, S-ducts, turboprop nacelles, wings, and wing-pylon-nacelle combinations. The solution of the full partial differential equation for compressible potential flow written in terms of a velocity potential is obtained using finite differences, line relaxation, and multigrid. The analysis uses either a cylindrical or Cartesian coordinate system. The computational mesh is not body fitted. The analysis has been programmed in FORTRAN for both the CDC CYBER 203 and the CRAY-1 computers. Comparisons of computed results with experimental measurement are presented. Descriptions of the program input and output formats are included.

Reyhner, T. A.↗

Tool for Generation of MAC/GMC Representative Unit Cell for CMC/PMC Analysis

This document describes a recently developed analysis tool that enhances the resident capabilities of the Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC) 4.0. This tool is especially useful in analyzing ceramic matrix composites (CMCs), where higher fidelity with improved accuracy of local response is needed. The tool, however, can be used for analyzing polymer matrix composites (PMCs) as well. MAC/GMC 4.0 is a composite material and laminate analysis software developed at NASA Glenn Research Center. The software package has been built around the concept of the generalized method of cells (GMC). The computer code is developed with a user friendly framework, along with a library of local inelastic, damage, and failure models. Further, application of simulated thermomechanical loading, generation of output results, and selection of architectures to represent the composite material have been automated to increase the user friendliness, as well as to make it more robust in terms of input preparation and code execution. Finally, classical lamination theory has been implemented within the software, wherein GMC is used to model the composite material response of each ply. Thus, the full range of GMC composite material capabilities is available for analysis of arbitrary laminate configurations as well. The primary focus of the current effort is to provide a graphical user interface (GUI) capability that generates a number of different user-defined repeating unit cells (RUCs). In addition, the code has provisions for generation of a MAC/GMC-compatible input text file that can be merged with any MAC/GMC input file tailored to analyze composite materials. Although the primary intention was to address the three different constituents and phases that are usually present in CMCs-namely, fibers, matrix, and interphase-it can be easily modified to address two-phase polymer matrix composite (PMC) materials where an interphase is absent. Currently, the tool capability includes generation of RUCs for square packing, hexagonal packing, and random fiber packing as well as RUCs based on actual composite micrographs. All these options have the fibers modeled as having a circular cross-sectional area. In addition, a simplified version of RUC is provided where the fibers are treated as having a square cross section and are distributed randomly. This RUC facilitates a speedy analysis using the higher fidelity version of GMC known as HFGMC. The first four mentioned options above support uniform subcell discretization. The last one has variable subcell sizes due to the primary intention of keeping the RUC size to a minimum to gain the speed ups using the higher fidelity version of MAC. The code is implemented within the MATLAB (The Mathworks, Inc., Natick, MA) developmental framework; however, a standalone application that does not need a priori MATLAB installation is also created with the aid of the MATLAB compiler.

Materials Engineering↗

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Post-Remediation Radiological Dose Assessment, Painesville Site, Painesville, Ohio

A post-remediation radiological dose assessment was conducted for the Formerly Utilized Sites Remedial Action Program (FUSRAP) Painesville Site by using the measured residual concentrations of the radionuclides of concern following the completion of remedial action. The site’s FUSRAP-related constituents of concern (COCs) are radionuclides associated with past activities conducted by the Manhattan Engineer District (MED) and Atomic Energy Commission (AEC) in support of the Nation’s early atomic energy and weapons program, and include radium-226 (Ra-226), thorium-230 (Th-230), thorium-232 (Th-232), and total uranium (Total U). Remedial actions were conducted at the site in accordance with the Final Record of Decision for the Painesville Site, Painesville, Ohio (ROD) (USACE 2006). In the ROD, the U.S. Army Corps of Engineers (USACE) determined that the U.S. Nuclear Regulatory Commission (NRC) standards for decommissioning of licensed facilities found in Title 10, Part 20.1402 of the Code of Federal Regulations (10 CFR 20.1402) and in the Ohio Administrative Code (OAC) 3701:1-38-22(B) are applicable or relevant and appropriate requirements (ARARs) for the site. On the basis of the ARAR analysis, a total effective dose equivalent (TEDE) goal of 25 millirem per year (mrem/yr) was assumed for the site with a construction worker considered as the average member of the critical group. Derived concentration guideline levels (DCGLs) for residual radioactivity of the COCs in soil were developed for the Painesville Site using the RESRAD (Version 6.22) computer code (Yu et al. 2001), which was designed to estimate radiation doses and risks from RESidual RADioactive materials. The DCGLs represent activity levels above site average background activity corresponding to the 25-mrem/yr criterion for a construction worker scenario. The DCGLs are above the background activity concentrations and represent average concentration guidelines. The DCGLs and average background activity concentrations are provided in Table 1. The Final Feasibility Study Addendum (USACE 2005) provides a complete description of the derivation of the DCGLs for the Painesville Site. A post-remediation radiological dose assessment was performed to determine the potential radiation doses to a construction worker, consistent with the approach used to develop the DCGLs for the Painesville Site given in the Final Feasibility Study Addendum (USACE 2005) and ROD (USACE 2006). A second more conservative dose assessment was performed to determine the potential radiation doses to an urban resident who is assumed to live and garden on the site. To further evaluate the long-term protectiveness of the remedial action, a third dose assessment was conducted to determine the potential radiation doses to a hypothetical subsistence farmer. The primary difference between exposure to an urban resident versus a conservative subsistence farmer would be an increase in the amount of foodstuffs (homegrown produce, milk, and beef) produced/cultivated and consumed on-site.

54 ENVIRONMENTAL SCIENCES↗

Computational Modeling, Simulation and Analysis of the RPI Walthousen Reactor Critical Facility (RCF)

The Walthousen Reactor Critical Facility (RCF), located in Schenectady, New York, is currently owned and operated by Rensselaer Polytechnic Institute (RPI). The facility was originally constructed by the American Locomotive Company (ALCO) for the research and development of small, portable reactors. Originally designed for highly-enriched uranium fuel, the prototype reactor first came online in August 26, 1956, and was transferred to RPI ownership in 1963. Several changes to the reactor's design and operation have been made, the most notable of which was the conversion from high-enriched uranium to low-enriched uranium fuel. The RCF is an open-pool critical assembly that operates at a low power with a license limit of less than 100 watts. This low operating power permits a much greater flexibility in configuration and operation as compared to other research reactors. Because of this, experiments may be designed that would be unfeasible for research reactors that operate at much higher powers. The RCF's functional capabilities provide opportunities to perform several reactor experiments and measurements, such as: core criticality experiments, differential and integral rod worth measurements, moderator temperature coefficient of reactivity calculations, critical rod height and critical water height measurements, etc.. These experiments are regularly performed in conjunction with RPI course-work education and Senior Reactor Operator (SRO) training. The unique capabilities and experiments that may be performed in the RCF are particularly useful for validating certain aspects of computational analysis codes. A computational model of the RCF reactor has been previously generated and validated using the Monte Carlo N-Particle (MCNP) code. For the purpose of future experiments, this model is also being adapted for use with the SHARP Toolset, which includes the PROTEUS neutronics code, the Nek5000 fluid dynamics code, and the Diablo structural mechanics code. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Effect of a localized minimum in equatorial field strength on resistive tearing instability in the geomagnetotail

A two-dimensional, resistive-MHD computer code is used to investigate the spontaneous reconnection of magnetotaillike configurations. The initial conditions adopted in the simulations are of two types: (1) in which the equatorial normal magnetic field component B(ze) declines monotonically down the tail, and (2) in which B(ze) exhibits a deep minimum in the near-earth plasma sheet. To represent the case where the earthward convection stops before the X line forms, zero-flow boundary conditions are imposed at the edges of the computational box. The initial configurations are in equilibrium and table within ideal MHD. The dynamic evolution of the system starts after the resistivity is turned on. The main results of these simulations basically support the neutral-line model of substorms and confirm Birn's (1980) computer studies. Spontaneous formation of an X-type neutral point and a single O-type plasmoid with strong tailward flow on the tailward side of the X point is found. The time interval from the turning on of the resistivity to the formation of a plasmoid is much shorter in the case where there is an initial deep minimum. A simple analytic calculation is also carried out to demonstrate why the configuration with a deep minimum is more susceptible to the development of the neutral point.

Hau, L.-N.↗

Pre-Flight Dark Forward Electrical Testing of the Mir Cooperative Solar Array

The Mir Cooperative Solar Array (MCSA) was developed jointly by the United States (US) and Russia to provide approximately 6 kW of photovoltaic power to the Russian space station Mir. After final assembly in Russia, the MCSA was shipped to the NASA Kennedy Space Center (KSC) in the summer of 1995 and launched to Mir in November 1995. Program managers were concerned of the potential for MCSA damage during the transatlantic shipment and the associated handling operations. To address this concern, NASA Lewis Research Center (LERC) developed an innovative dark-forward electrical test program to assess the gross electrical condition of each generator following shipment from Russia. The use of dark test techniques, which allow the array to remain in the stowed configuration, greatly simplifies the checkout of large area solar arrays. MCSA dark electrical testing was successfully performed at KSC in July 1995 following transatlantic shipment. Data from this testing enabled engineers to quantify the effects of potential MCSA physical damage that would degrade on-orbit electrical performance. In this paper, an overview of the principles and heritage of photovoltaic array dark testing is given. The specific MCSA dark test program is also described including the hardware, software, testing procedures and test results. The current-voltage (4) response of both solar cell circuitry and by-pass diode circuitry was obtained. To guide the development of dark test hardware, software and procedures, a dedicated FORTRAN computer code was developed to predict the dark 4 responses of generators with a variety of feasible damage modes. By comparing the actual test data with the predictions, the physical condition of the generator could be inferred. Based on this data analysis, no electrical short-circuits or open-circuits were detected. This suggested the MCSA did not sustain physical damage that affected electrical performance during handling and shipment from Russia to the US. Good agreement between the test data and computational predictions indicated MCSA electrical performance was amenable to accurate analysis and was well understood.

Kerslake, Thomas W.↗

NSRD-06. Computational Capability to Substantiate DOE-HDBK-3010 Data

Safety basis analysts throughout the U.S. Department of Energy (DOE) complex rely heavily on the information provided in the DOE Hand book, DOE-HDBK-3010, Airborne Release Fractions/Rates and Resp irable Fractions for Nonreactor Nuclear Facilities , to determine source terms. In calcula ting source terms, analysts tend to use the DOE Handbook's bounding values on airbor ne release fractions (ARFs) and respirable fractions (RFs) for various cat egories of insults (representing potential accident release categories). This is typica lly due to both time constraints and the avoidance of regulatory critique. Unfort unately, these bounding ARFs/RFs represent extremely conservative values. Moreover, th ey were derived from very limited small- scale table-top and bench/labo ratory experiments and/or fr om engineered judgment. Thus the basis for the data may not be re presentative to the actual unique accident conditions and configura tions being evaluated. The goal of this res earch is to develop a more ac curate method to identify bounding values for the DOE Handbook using the st ate-of-art multi-physics-based high performance computer codes. This enable s us to better understand the fundamental physics and phenomena associated with the ty pes of accidents for the data described in it. This research has examined two of the DOE Handbook's liquid fire experiments to substantiate the airborne release frac tion data. We found th at additional physical phenomena (i.e., resuspension) need to be included to derive bounding values. For the specific cases of solid powder under pre ssurized condition and mechanical insult conditions the codes demonstrated that we can simulate the phenomena. This work thus provides a low-cost method to establis h physics-justified sa fety bounds by taking into account specific geometri es and conditions that may not have been previously measured and/or are too costly to do so.

97 MATHEMATICS AND COMPUTING↗