Search NASA⌕ Search

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.

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

At least 1,387 records · Page 77

Innovative design of composite structures: The use of curvilinear fiber format in structural design of composites

The influence is investigated of a curvilinear fiber format on load carrying capacity of a layered fiber reinforced plate with a centrally located hole. A curvilinear fiber format is descriptive of layers in a laminate having fibers which are aligned with the principal stress directions in those layers. Laminates of five curvilinear fiber format designs and four straightline fiber format designs are considered. A quasi-isotropic laminate having a straightline fiber format is used to define a baseline design for comparison with the other laminate designs. Four different plate geometries are considered and differentiated by two values of hole diameter/plate width equal to 1/6 and 1/3, and two values of plate length/plate width equal to 2 and 1. With the plates under uniaxial tensile loading on two opposing edges, alignment of fibers in the curvilinear layers with the principal stress directions is determined analytically by an iteration procedure. In-plane tensile load capacity is computed for all of the laminate designs using a finite element analysis method. A maximum strain failure criterion and the Tsai-Wu failure criterion are applied to determine failure loads and failure modes. Resistance to buckling of the laminate designs to uniaxial compressive loading is analyzed using the commercial code Engineering Analysis Language. Results indicate that the curvilinear fiber format laminates have higher in-plane tensile load capacity and comparable buckling resistance relative to the straightline fiber format laminates.

Charette, R. F.↗

Quantitative Analysis of Color Differences within High Contrast, Low Power Reversible Electrophoretic Displays

Here, quantitative color analysis is performed on videos of high contrast, low power reversible electrophoretic deposition (EPD)-based displays operated under different applied voltages. This analysis is coded in an open-source software, relies on a color differentiation metric, ΔE * 00 , derived from digital video, and provides an intuitive relationship between the operating conditions of the devices and their performance. Time-dependent ΔE * 00 color analysis reveals color relaxation behavior, recoverability for different voltage sequences, and operating conditions that can lead to optimal performance.

36 MATERIALS SCIENCE↗

Conjugate Heat Transfer Analysis in Pressurized Helium Gas Cooling Channels

Los Alamos National Laboratory (LANL) is working with NorthStar Medical Radioisotopes (NMR) on their efforts to produce Mo-99 from the irradiation of Mo-100 targets. The NMR target consists of an Inconel window that allows the electron beam to penetrate and irradiate a stack of Mo-100 disks. The irradiation process generates large amounts of heat and is cooled using pressurized helium gas flowing at 400 g/s. LANL provides both numerical analysis support and experimental support. Computational fluid dynamics (CFD) and conjugate heat transfer (CHT) analysis are performed on the NMR target to analyze the helium flow and volumetric heating through the NMR system. In this iteration of numerical analysis, we update the Mo-100 disk thickness to 0.74 mm and the gap between the Inconel window and the first disk to 0.69 mm. Estimates of the heat flux from the electron beams from Monte Carlo N-Particle code (MCNP) analysis were combined with CFD to perform the CHT analysis. The pressure drop through the targets, force applied on the first disk, beam window and target disk temperatures are some of the quantities estimated using the simulations. The pressure drop through the target was estimated as 0.198 MPa, while the force due to venturi effects across the first disk was estimated to be 100 N. The results show maximum window temperatures over 600 °C and peak disk temperatures over 1600 °C. Past work has shown that beyond 650-700 °C, the Inconel window becomes susceptible to yield and rupture. The increased temperature in these results reduces the margin of error possible on beam spot size. The disk temperatures showed a sinusoidal distribution with the 6 th – 8 th disk from each window exhibiting highest disk temperatures. The peak temperature of the stainless steel housing was close to 316 °C and that of the laminations and spacers was approximately 582 °C.

43 PARTICLE ACCELERATORS↗

Risk-Informed Mechanistic Source Term Calculation (Final CRADA Report)

As part of an FY18 U.S. Department of Energy (DOE) Gateway for Accelerated Innovation in Nuclear (GAIN) award, a CRADA agreement was established between Argonne National Laboratory (Argonne) and Oklo. The CRADA focused on the following three tasks to assist Oklo with the regulatory licensing of their innovative nuclear reactor design: • Task 1 – Coupling and Demonstration Analysis using The Nek5000, Diablo, and OpenMC Codes: Establishing a coupling between the neutronic, structural, and thermohydraulic analysis codes and utilizing the code package for a high-fidelity demonstration analysis of radial core expansion due to temperature changes and its impact on core reactivity. • Task 2 – Mechanistic Source Term Analysis: Performing a realistic evaluation of the potential release of radionuclides to the environment during postulated transient event sequences. • Task 3 – Data in Support of Metal Fuel Qualification: An assessment of past postirradiation examination data of metal fuel utilized in the Argonne EBR-II reactor to assist in the potential qualification of the fuel typer for regulatory licensing. The details regarding each analysis and the associated conclusions were provided to Oklo in three separate project reports. Due to the utilization of proprietary and export controlled information, the current summary report was created, which provides an uncontrolled overview of each of the three tasks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Progressive Damage Analysis of Laminated Composite (PDALC)-A Computational Model Implemented in the NASA COMET Finite Element Code

A method for analysis of progressive failure in the Computational Structural Mechanics Testbed is presented in this report. The relationship employed in this analysis describes the matrix crack damage and fiber fracture via kinematics-based volume-averaged variables. Damage accumulation during monotonic and cyclic loads is predicted by damage evolution laws for tensile load conditions. The implementation of this damage model required the development of two testbed processors. While this report concentrates on the theory and usage of these processors, a complete list of all testbed processors and inputs that are required for this analysis are included. Sample calculations for laminates subjected to monotonic and cyclic loads were performed to illustrate the damage accumulation, stress redistribution, and changes to the global response that occur during the load history. Residual strength predictions made with this information compared favorably with experimental measurements.

Lo, David C.↗

Progressive Damage Analysis of Laminated Composite (PDALC) (A Computational Model Implemented in the NASA COMET Finite Element Code)

A method for analysis of progressive failure in the Computational Structural Mechanics Testbed is presented in this report. The relationship employed in this analysis describes the matrix crack damage and fiber fracture via kinematics-based volume-averaged damage variables. Damage accumulation during monotonic and cyclic loads is predicted by damage evolution laws for tensile load conditions. The implementation of this damage model required the development of two testbed processors. While this report concentrates on the theory and usage of these processors, a complete listing of all testbed processors and inputs that are required for this analysis are included. Sample calculations for laminates subjected to monotonic and cyclic loads were performed to illustrate the damage accumulation, stress redistribution, and changes to the global response that occurs during the loading history. Residual strength predictions made with this information compared favorably with experimental measurements.

Coats, Timothy W.↗

Mixing with applications to inertial-confinement-fusion implosions

Approximate one-dimensional (1D) as well as 2D and 3D simulations are playing an important supporting role in the design and analysis of future experiments at National Ignition Facility. This work is mainly concerned with 1D simulations, used extensively in design and optimization. We couple a 1D buoyancy-drag mix model for the mixing zone edges with a 1D inertial confinement fusion simulation code. This analysis predicts that National Ignition Campaign (NIC) designs are located close to a performance cliff, so modeling errors, design features (fill tube and tent) and additional, unmodeled instabilities could lead to significant levels of mix. The performance cliff we identify is associated with multimode plastic ablator (CH) mix into the hot-spot deuterium and tritium (DT). The buoyancy-drag mix model is mode number independent and selects implicitly a range of maximum growth modes. Our main conclusion is that single effect instabilities are predicted not to lead to hot-spot mix, while combined mode mixing effects are predicted to affect hot-spot thermodynamics and possibly hot-spot mix. Combined with the stagnation Rayleigh-Taylor instability, we find the potential for mix effects in combination with the ice-to-gas DT boundary, numerical effects of Eulerian species CH concentration diffusion, and ablation-driven instabilities. With the help of a convenient package of plasma transport parameters developed here, we give an approximate determination of these quantities in the regime relevant to the NIC experiments, while ruling out a variety of mix possibilities. Plasma transport parameters affect the 1D buoyancy-drag mix model primarily through its phenomenological drag coefficient as well as the 1D hydro model to which the buoyancy-drag equation is coupled.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Report on Automated Semantic Analysis of Scientific and Engineering Codes

The loss of the Mars Climate Orbiter due to a software error reveals what insiders know: software development is difficult and risky because, in part, current practices do not readily handle the complex details of software. Yet, for scientific software development the MCO mishap represents the tip of the iceberg; few errors are so public, and many errors are avoided with a combination of expertise, care, and testing during development and modification. Further, this effort consumes valuable time and resources even when hardware costs and execution time continually decrease. Software development could use better tools! This lack of tools has motivated the semantic analysis work explained in this report. However, this work has a distinguishing emphasis; the tool focuses on automated recognition of the fundamental mathematical and physical meaning of scientific code. Further, its comprehension is measured by quantitatively evaluating overall recognition with practical codes. This emphasis is necessary if software errors-like the MCO error-are to be quickly and inexpensively avoided in the future. This report evaluates the progress made with this problem. It presents recommendations, describes the approach, the tool's status, the challenges, related research, and a development strategy.

Stewart. Maark E. M.↗

Cooperation on Improved Isotopic Identification and Analysis Software for Portable, Electrically Cooled High-Resolution Gamma Spectrometry Systems Final Report

Under a 2006 agreement between the Department of Energy (DOE) of the United States of America and the Institut de Radioprotection et de Sûreté Nucléaire (IRSN) of France, the National Nuclear Security Administration (NNSA) within DOE and IRSN initiated a collaboration to improve isotopic identification and analysis of nuclear material [i.e., plutonium (Pu) and uranium (U)]. The specific aim of the collaborative project was to develop new versions of two types of isotopic identification and analysis software: (1) the fixed-energy response-function analysis for multiple energies (FRAM) codes and (2) multi-group analysis (MGA) codes. The project is entitled Action Sheet 4 – Cooperation on Improved Isotopic Identification and Analysis Software for Portable, Electrically Cooled, High-Resolution Gamma Spectrometry Systems (Action Sheet 4). FRAM and MGA/U235HI are software codes used to analyze isotopic ratios of U and Pu. FRAM is an application that uses parameter sets for the analysis of U or Pu. MGA and U235HI are two separate applications that analyze Pu or U, respectively. They have traditionally been used by safeguards practitioners to analyze gamma spectra acquired with high-resolution gamma spectrometry (HRGS) systems that are cooled by liquid nitrogen. However, it was discovered that these analysis programs were not as accurate when used on spectra acquired with a newer generation of more portable, electrically cooled HRGS (ECHRGS) systems. In response to this need, DOE/NNSA and IRSN collaborated to update the FRAM and U235HI codes to improve their performance with newer ECHRGS systems. Lawrence Livermore National Laboratory (LLNL) and Los Alamos National Laboratory (LANL) performed this work for DOE/NNSA.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Laser Anemometer Measurements of the Flow Field in a 4:1 Pressure Ratio Centrifugal Impeller

A laser-doppler anemometer was used to obtain flow-field velocity measurements in a 4:1 pressure ratio, 4.54 kg/s (10 lbm/s), centrifugal impeller, with splitter blades and backsweep, which was configured with a vaneless diffuser. Measured through-flow velocities are reported for ten quasi-orthogonal survey planes at locations ranging from 1% to 99% of main blade chord. Measured through-flow velocities are compared to those predicted by a 3-D viscous steady flow analysis (Dawes) code. The measurements show the development and progression through the impeller and vaneless diffuser of a through-flow velocity deficit which results from the tip clearance flow and accumulation of low momentum fluid centrifuged from the blade and hub surfaces. Flow traces from the CFD analysis show the origin of this deficit which begins to grow in the inlet region of the impeller where it is first detected near the suction surface side of the passage. It then moves toward the pressure side of the channel, due to the movement of tip clearance flow across the impeller passage, where it is cut by the splitter blade leading edge. As blade loading increases toward the rear of the channel the deficit region is driven back toward the suction surface by the cross-passage pressure gradient. There is no evidence of a large wake region that might result from flow separation and the impeller efficiency is relatively high. The flow field in this impeller is quite similar to that documented previously by NASA Lewis in a large low-speed backswept impeller.

Skoch, G. J.↗

Efficient Multidisciplinary Analysis Procedure Using Multi-Level Parallelization Approach

Multidisciplinary applications are suitable for parallel computing environment by adopting the domain decomposition method. Immediately, a multidisciplinary application can be parallelized by solving each discipline separately. In order to perform coupled multidisciplinary analysis, coupling of each discipline can be accomplished by exchanging boundary data at the interfaces. This is regarded as discipline-level parallelization. Next level could be a "coarse-grain" parallelization of each discipline, which mainly depends on the physical geometry and nature of each discipline. For example, it is almost impossible for structured-grid based computational fluid dynamics codes to do flow analysis of an aircraft by using a single grid because of the complexity of its configuration. Thus, multi-block grid is commonly used to describe the details of complex geometry. Similarly, in structural analysis, the structure is frequently subdivided into substructures. Thus, the computation of each subdomain can be easily parallelized since each subdomain is solved separately independent of other domains. The parallelization is accomplished by solving each subdomain separately on a separate processor and exchanging the boundary conditions at domain interfaces periodically. However, the physical decomposition of the domain introduces explicit boundary conditions at the domain interfaces. This is not desirable for critical areas such as those containing shock waves or flow separations. Thus, a "fine-grain" parallelization is introduced to overcome this problem. The "fine-grain" parallelization is one that solves exactly the same system of equations of a subdomain by using more than one processors without introducing any explicit boundary conditions. An efficient multidisciplinary analysis procedure can be accomplished by successfully combining the above multi-level parallelism. A multidisciplinary analysis code, ENSAERO developed at NASA Ames Research Center is used in this study to implement the proposed approach. The communication data structure required for the proposed approach will be studied in detail. This work will demonstrate the feasibility of using multi-level parallelization approach in multidisciplinary analysis applications.

Byun, Chansup↗

Kingdom-Wide Analysis of Fungal Protein-Coding and tRNA Genes Reveals Conserved Patterns of Adaptive Evolution

Protein-coding genes evolved codon usage bias due to the combined but uneven effects of adaptive and nonadaptive influences. Studies in model fungi agree on codon usage bias as an adaptation for fine-tuning gene expression levels; however, such knowledge is lacking for most other fungi. Our comparative genomics analysis of over 450 species supports codon usage and transfer RNAs (tRNAs) as coadapted for translation speed and this is most likely a realization of convergent evolution. Rather than drift, phylogenetic reconstruction inferred adaptive radiation as the best explanation for the variation of interspecific codon usage bias. Although the phylogenetic signals for individual codon and tRNAs frequencies are lower than expected by genetic drift, we found remarkable conservation of highly expressed genes being codon optimized for translation by the most abundant tRNAs, especially by inosine-modified tRNAs. As an application, we present a sequence-to-expression neural network that uses codons to reliably predict highly expressed transcripts. The kingdom Fungi, with over a million species, includes many key players in various ecosystems and good targets for biotechnology. Collectively, our results have implications for better understanding the evolutionary success of fungi, as well as informing the biosynthetic manipulation of fungal genes.

59 BASIC BIOLOGICAL SCIENCES↗

Quantitative radiography for determining density fluctuations in HED experiments

We have developed a method to extract density fluctuation measurements from x-ray radiographs of high-energy density (HED) instability growth and turbulence experiments. We use this information to calculate density fluctuation statistics for constraining the performance of turbulent mix models in HED systems. The density calculation combines image filtering, removal of systemic effects such as backlighter variation, calculation of transmission across multiple materials, and use of tracer materials to generate an approximate single-material density field. From the density map, we calculate both average density and a variance-like moment b (density-specific-volume covariance), which we compare to our models. We infer both quantities from a single image, which is significantly more information than the historic single scalar mix width measurements. We also develop a method of analyzing simulation outputs that incorporate both the density fluctuation metric from a turbulence model and the bulk material maps from the hydrodynamic code. This analysis helps address the question of how to initialize the simulations for best comparison to data from systems with large separations of scale in the mixing perturbation initial condition. We find that our data analysis method yields 1D average density and b curves with similar morphology and amplitudes as those from preliminary simulation comparisons.

47 OTHER INSTRUMENTATION↗

Rotordynamics on the PC: Further Capabilities of ARDS

Rotordynamics codes for personal computers are now becoming available. One of the most capable codes is Analysis of RotorDynamic Systems (ARDS) which uses the component mode synthesis method to analyze a system of up to 5 rotating shafts. ARDS was originally written for a mainframe computer but has been successfully ported to a PC; its basic capabilities for steady-state and transient analysis were reported in an earlier paper. Additional functions have now been added to the PC version of ARDS. These functions include: 1) Estimation of the peak response following blade loss without resorting to a full transient analysis; 2) Calculation of response sensitivity to input parameters; 3) Formulation of optimum rotor and damper designs to place critical speeds in desirable ranges or minimize bearing loads; 4) Production of Poincard plots so the presence of chaotic motion can be ascertained. ARDS produces printed and plotted output. The executable code uses the full array sizes of the mainframe version and fits on a high density floppy disc. Examples of all program capabilities are presented and discussed.

Fleming, David P.↗

FPA Depot - Web Application

Lines of code (LOC) analysis is one of the methods used to measure programmer productivity and estimate schedules of programming projects. The Launch Control System (LCS) had previously used this method to estimate the amount of work and to plan development efforts. The disadvantage of using LOC as a measure of effort is that one can only measure 30% to 35% of the total effort of software projects involves coding [8]. In the application, instead of using the LOC we are using function point for a better estimation of hours in each software to develop. Because of these disadvantages, Jamie Szafran of the System Software Branch of Control And Data Systems (NE-C3) at Kennedy Space Canter developed a web application called Function Point Analysis (FPA) Depot. The objective of this web application is that the LCS software architecture team can use the data to more accurately estimate the effort required to implement customer requirements. This paper describes the evolution of the domain model used for function point analysis as project managers continually strive to generate more accurate estimates.

Avila, Edwin M. Martinez↗

Solution of the Burnett equations for hypersonic flows near the continuum limit

The INCA code, a three-dimensional Navier-Stokes code for analysis of hypersonic flowfields, was modified to analyze the lower reaches of the continuum transition regime, where the Navier-Stokes equations become inaccurate and Monte Carlo methods become too computationally expensive. The two-dimensional Burnett equations and the three-dimensional rotational energy transport equation were added to the code and one- and two-dimensional calculations were performed. For the structure of normal shock waves, the Burnett equations give consistently better results than Navier-Stokes equations and compare reasonably well with Monte Carlo methods. For two-dimensional flow of Nitrogen past a circular cylinder the Burnett equations predict the total drag reasonably well. Care must be taken, however, not to exceed the range of validity of the Burnett equations.

Imlay, Scott T.↗

Metallurgical Analysis of the High Flux Isotope Reactor (HFIR) Carrier Lifting Bails (Rev.1)

The dissolution rates of the aluminum alloys in the High Flux Isotope Reactor (HFIR) element carriers and the Material Test Reactor (MTR) L-bundles in the H-Canyon facility have been identified as the possible cause of extended dissolutions that result in significant time and financial expenditures. A study, carried out by Savannah River National Laboratory (SRNL) to determine relationships between the dissolution rates and the metallurgical properties of the aluminum alloy materials of construction of the HFIR carriers and the L-bundles, considered the dissolution rates of aluminum alloy (AA) series 1100, 6061, and 6063. The study determined that the aluminum alloy compositions played a principal role in the dissolution rate of the carrier/bundle components. Higher dissolution rates were correlated with lower concentrations of the minor element additions in the alloys and with specific element concentrations. Aluminum alloys 1100 and 6063 were found to have similar dissolution rates that were approximately two orders of magnitude (100X) greater than those of AA6061. Based on the results of the dissolution behavior study, a Technical Assistance Request (TAR) was first issued to determine if the replacement of AA6061-T6 with AA6063-T6 is feasible for the HFIR carrier lifting bails. A Technical Task Request was then issued to consider AA6063-T5 as well as other alloys to improve possible supply chain issues. The metallurgical properties of the L-bundle (specifically the end caps) were not evaluated in this report because L-Bundle drawings already allow for the use of AA6063-T6 in all structural components. The HFIR carriers are composed of thin-walled components with significant surface areas that allow for relatively quick overall dissolution times. Conversely, the carrier lifting bails and the supporting constituents are composed of solid bars and thick plate regions with relatively small surface areas that experience longer overall dissolution times. While the MTR L-bundle design includes allowances for the materials of construction to be either AA6061-T6 or AA6063-T6, the HFIR carriers are specified to be constructed fully with AA6061-T6 alloy. This report analyzes the recommendations of the dissolution behavior study to replace the materials of construction of the HFIR carrier lifting bails. The analysis considers the operational requirements of the lifting bail and its supporting structures. To decrease dissolution times, the analysis considers direct replacement of the material as well as reductions in the thicknesses of the components to decrease the mass of the elements. Material reductions are considered on options for using either AA6061 and/or AA6063. The calculations are based on specifications from the American Society of Mechanical Engineer (ASME) and The Aluminum Association, Inc. design codes. The analysis finds that direct replacement of the lifting bail material of construction with AA6063-T6, and AA6063-T5 as well as reductions in the dimensions of the lifting bail components are acceptable. Note that this study considers the structural suitability of the alloys. It does not consider their dissolution rates in the dissolvers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

1025:1 Area Ratio Nozzle Evaluated at High Combustion Chamber Pressure

A recently completed experimental test program obtained performance data on an optimally contoured nozzle with an exit-to-throat area ratio of 1025:1 and on a truncated version of this nozzle with an area ratio of 440:1. The nozzles were tested with gaseous hydrogen and liquid oxygen propellants at combustion chamber pressures of 12.4 to 16.5 mPa (1800 to 2400 psia). Testing was conducted in the altitude test capsule at the NASA Lewis Research Center's Rocket Engine Test Facility (RETF), and results were compared with analytical performance predictions. This testing builds on previous work with this nozzle at Lewis, where testing was completed at nominal chamber pressure of 350 psia. High-area-ratio nozzles have long been sought as a means to increase the performance of spacebased rocket engines. However, as the area ratio increases, the physical size and weight of the nozzle also increase. As a result, engine and vehicle designers must make tradeoffs between nozzle size and performance enhancement. Until this test program, very little experimental data existed on the performance of the high-area-ratio nozzles used in rocket engine designs. The computer codes being used by rocket engine designers rely on data extrapolated from tests of low-area-ratio nozzles, and these extrapolations do not always provide the accuracy needed for a reliable design assessment. Therefore, we conducted this high-area-ratio nozzle testing program to provide performance data for use in rocket engine design and analysis computer codes. The nozzle had a nominal 2.54-cm- (1-in.-) diameter throat, an exit diameter of 81.3-cm (32.0-in.) at an exit-to-throat area ratio of 1025, and a length of 128.6 cm (50.6 in.). Testing was conducted in an altitude test capsule to simulate the static pressure at altitude by vacuum pumping. Data such as propellant mass flow, oxidizer-to-fuel mixture, and thrust were measured. These measurements were then used to calculate performance factors such as the thrust coefficient, the characteristic exhaust velocity efficiency, and the vacuum specific impulse. In addition, the nozzle temperature was measured to calculate the amount of heat transferred from the combustion gases to the nozzle.

Source record↗