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,333 records · Page 74

LLNL NESHAPs 2023 Annual Report

Lawrence Livermore National Security, LLC operates facilities at Lawrence Livermore National Laboratory (LLNL) in which radionuclides are handled and stored. These facilities are subject to the U.S. Environmental Protection Agency (EPA) National Emission Standards for Hazardous Air Pollutants (NESHAPs) in Code of Federal Regulations (CFR) Title 40, Part 61, Subpart H, which regulates radionuclide emissions to air from Department of Energy (DOE) facilities. Specifically, NESHAPs limits the emission of radionuclides to the ambient air to levels resulting in an annual effective dose equivalent of 10 mrem (100 µSv) to any member of the public. Using measured and calculated emissions, and building-specific and common parameters, LLNL personnel applied the EPA-approved computer code, CAP88-PC, Version 4.1.1.0, to calculate the dose to the maximally exposed individual member of the public for the Livermore Site and Site 300.

54 ENVIRONMENTAL SCIENCES↗

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.↗

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.↗

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↗

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.↗

SAS4A/SASSYS-1 Validation with the EBR-II Tests SHRT-17 and SHRT-45R

The EBR-II sodium-cooled fast reactor supported a wide range of fast reactor research activities. The final fifteen years of operations at EBR-II were used for experiments and tests to demonstrate the importance of passive safety in liquid metal reactors and the capability of such systems to provide strong passive responses during off-normal conditions. The Shutdown Heat Removal Test program at EBR-II provided test data supporting the validation of computer codes for design, licensing, and operation of LMRs, among other objectives. The program included nearly sixty tests, most of which were protected and unprotected loss of flow and loss of heat sink tests run at various initial powers and flow rates.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Kinetics calculation of fast periodic pulsed reactors using MCNP6

Fast periodic pulsed reactor is a type of reactor in which the fission bursts are formed entirely with external reactivity modulation with a specified time periodicity. This type of reactors could generate much larger intensity of neutron beams for experimental use, compared with the steady state reactors. In the design of fast periodic pulsed reactors, the time dependent simulation of the power pulse is majorly based on a point kinetic model, which is known to have limitations. A more accurate calculation method is desired for the design analyses of fast periodic pulsed reactors. Monte Carlo computer code MCNP6 is used for this task due to its three dimensional transport capability with a continuous energy library. Some new routines were added to simulate the rotation of the movable reflector parts in the time dependent calculation. Fast periodic pulsed reactor IBR-2M was utilized to validate the new routines. This reactor is periodically in prompt supercritical state, which lasts for similar to 400 mu s, during the equilibrium state. This generates long neutron fission chains, which requires tremendously large amount of computation time during Monte Carlo simulations. Russian Roulette was applied for these very long neutron chains in MCNP6 calculation, combined with other approaches to improve the efficiency of the simulations. In the power pulse of the IBR-2M at equilibrium state, there is some discrepancy between the experimental measurements and the calculated results using the point kinetics model. MCNP6 results matches better the experimental measurements, which shows the merit of using MCNP6 calculation relative to the point kinetics model.

Zhong, Z.↗

Computation of nonequilibrium three-dimensional inviscid flow over blunt-nosed bodies flying at supersonic speeds

A computer code based on the method of characteristics is described and applied to the study of two- and three-dimensional chemical nonequilibrium flow over sharp and blunt-nosed bodies. Nonequilibrium flow over a wedge is used to show the approach to equilibrium flow, and to demonstrate the nature of the reaction zone behind the bow shock wave. The structure and development of a blunt-body entropy layer in nonequilibrium flow is examined for a blunt cone at zero incidence. Three-dimensional computations for the space shuttle body at 30 deg angle of attack are presented. A nondimensional scaling parameter, the Damkoehler number, which is the ratio of flow time to chemical reaction time, is calculated and its significance discussed.

Rakich, J. V.↗

Non-dipolar magnetic field models and patterns of radio emission: Uranus and Neptune compared

The magnetic field geometries of Uranus and Neptune are superficially similar, and are similarly unlike those of other planets: the field strengths are similar, and they contain extraordinarily large non-dipolar components. As a corollary, the best dipolar field models of each of the two planets comprises a dipole that is considerably offset from the planetary center and tilted away from the rotational axis. However, in other respects the best field models of the two planets are quite different. Uranus has a quadrupole model in which all the terms are well determined and in which none of the higher order terms is determined. To represent the magnetometer data acquired during Voyager's Neptune encounter requires a model of order 8 (instead of Uranus' order 2), yet many of the coefficients are poorly determined. A second model, an octupole model comprising the terms up to order three of the order 8 model, has been suggested by the magnetometer team as being useful; its use, however, is limited only to the region outside of about 2R(exp N), whereas planetary radio emissions have their sources well inside this surface. Computer code has been written that permits an analysis of the detailed motion of low energy charged particles moving in general planetary magnetic fields. At Uranus, this code reveals the existence of an isolated region of the inner magnetosphere above the day side in which particles may be trapped, separate from the more general magnetospheric trapping. An examination of the so-call ordinary mode uranian radio emissions leads us to believe that these emissions are in fact extraordinary mode emissions coming from particles trapped in this isolated region. A similar attempt to discover trapping regions at Neptune has proved, unfortunately, to be impossible. This arises from three factors: (1) the computation needed to track particles in an eighth order field is more than an order of magnitude greater than that needed to perform a similar calculation in a quadrupole field, and is beyond the capacity of workstation-class computers; (2) the octupole field model is known to be in error by too large an ammount for it, or any similarly truncated version of the eighth order model, to produce trustworthy results; (3) the eighth order model can, in effect, be infinitely varied without affecting the field strength along the spacecraft trajectory.

Evans, D. R.↗