Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Static Code Analysis”

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 163 records · Page 9

Effect of Microscopic Damage Events on Static and Ballistic Impact Strength of Triaxial Braid Composites

The reliability of impact simulations for aircraft components made with triaxial-braided carbon-fiber composites is currently limited by inadequate material property data and lack of validated material models for analysis. Methods to characterize the material properties used in the analytical models from a systematically obtained set of test data are also lacking. A macroscopic finite element based analytical model to analyze the impact response of these materials has been developed. The stiffness and strength properties utilized in the material model are obtained from a set of quasi-static in-plane tension, compression and shear coupon level tests. Full-field optical strain measurement techniques are applied in the testing, and the results are used to help in characterizing the model. The unit cell of the braided composite is modeled as a series of shell elements, where each element is modeled as a laminated composite. The braided architecture can thus be approximated within the analytical model. The transient dynamic finite element code LS-DYNA is utilized to conduct the finite element simulations, and an internal LS-DYNA constitutive model is utilized in the analysis. Methods to obtain the stiffness and strength properties required by the constitutive model from the available test data are developed. Simulations of quasi-static coupon tests and impact tests of a represented braided composite are conducted. Overall, the developed method shows promise, but improvements that are needed in test and analysis methods for better predictive capability are examined.

Littell, Justin D.↗

Precise and Efficient Static Array Bound Checking for Large Embedded C Programs

In this paper we describe the design and implementation of a static array-bound checker for a family of embedded programs: the flight control software of recent Mars missions. These codes are large (up to 250 KLOC), pointer intensive, heavily multithreaded and written in an object-oriented style, which makes their analysis very challenging. We designed a tool called C Global Surveyor (CGS) that can analyze the largest code in a couple of hours with a precision of 80%. The scalability and precision of the analyzer are achieved by using an incremental framework in which a pointer analysis and a numerical analysis of array indices mutually refine each other. CGS has been designed so that it can distribute the analysis over several processors in a cluster of machines. To the best of our knowledge this is the first distributed implementation of static analysis algorithms. Throughout the paper we will discuss the scalability setbacks that we encountered during the construction of the tool and their impact on the initial design decisions.

Venet, Arnaud↗

Stability analysis of a reinforced carbon carbon shell

This paper presents the development of a stability analysis for the nose cap of the NASA Space Shuttle Orbiter. Stability is evaluated by the differential stiffness analysis of the NASTRAN finite-element computer code, addressing those nonstandard characteristics in the nose cap such as nonuniform curvature, asymmetrical and nonuniform loads, support fixity, and various combinations of membrane and bending stresses. A full-sized nose cap, thinner than production, was statically tested and stability analyzed. The failing load level correlated to within 30%. The region and mode of buckling that occurred during test was accurately predicted by analysis. The criterion for predicting instability is based on the behavior of the nonlinear deflections. The deflections are nonlinear elastic in that the stresses are well within the elastic range of the material, but the geometry-load relationship produces nonlinear deflections. The load-deflection relationship is well defined by differential stiffness analysis up to the zero-slope portion of the curve, the point of neutral stability or where the shell 'snaps through' just prior to general instability.

Agan, W. E.↗

Performance analysis of fully explicit and fully implicit solvers within a spectral element shallow-water atmosphere model

In this study, explicit Runge–Kutta methods and implicit multistep methods utilizing a Newton–Krylov nonlinear solver are evaluated for a range of configurations of the shallow-water dynamical core of the spectral element community atmosphere model to evaluate their computational performance. These configurations are designed to explore the attributes of each method under different but relevant model usage scenarios including varied spectral order within an element, static regional refinement, and scaling to the largest problem sizes. This analysis is performed within the shallow-water dynamical core option of a full climate model code base to enable a wealth of simulations for study, with the aim of informing solver development within the more complete hydrostatic dynamical core used for climate research. The limitations and benefits to using explicit versus implicit methods, with different parameters and settings, are discussed in light of the trade-offs with Message Passing Interface (MPI) communication and memory and their inherent efficiency bottlenecks. Given the performance behavior across the configurations analyzed here, the recommendation for future work using the implicit solvers is conditional based on scale separation and the stiffness of the problem. For the regionally refined configurations, the implicit method has about the same efficiency as the explicit method, without considering efficiency gains from a preconditioner. The potential for improvement using a preconditioner is greatest for higher spectral order configurations, where more work is shifted to the linear solver. Finally, initial simulations with OpenACC directives to utilize a Graphics Processing Unit (GPU) when performing function evaluations show improvements locally, and that overall gains are possible with adjustments to data exchanges.

58 GEOSCIENCES↗

FASOR - A second generation shell of revolution code

An integrated computer program entitled Field Analysis of Shells of Revolution (FASOR) currently under development for NASA is described. When completed, this code will treat prebuckling, buckling, initial postbuckling and vibrations under axisymmetric static loads as well as linear response and bifurcation under asymmetric static loads. Although these modes of response are treated by existing programs, FASOR extends the class of problems treated to include general anisotropy and transverse shear deformations of stiffened laminated shells. At the same time, a primary goal is to develop a program which is free of the usual problems of modeling, numerical convergence and ill-conditioning, laborious problem setup, limitations on problem size and interpretation of output. The field method is briefly described, the shell differential equations are cast in a suitable form for solution by this method and essential aspects of the input format are presented. Numerical results are given for both unstiffened and stiffened anisotropic cylindrical shells and compared with previously published analytical solutions.

Cohen, G. A.↗

Reducing False Positives in Runtime Analysis of Deadlocks

This paper presents an improvement of a standard algorithm for detecting dead-lock potentials in multi-threaded programs, in that it reduces the number of false positives. The standard algorithm works as follows. The multi-threaded program under observation is executed, while lock and unlock events are observed. A graph of locks is built, with edges between locks symbolizing locking orders. Any cycle in the graph signifies a potential for a deadlock. The typical standard example is the group of dining philosophers sharing forks. The algorithm is interesting because it can catch deadlock potentials even though no deadlocks occur in the examined trace, and at the same time it scales very well in contrast t o more formal approaches to deadlock detection. The algorithm, however, can yield false positives (as well as false negatives). The extension of the algorithm described in this paper reduces the amount of false positives for three particular cases: when a gate lock protects a cycle, when a single thread introduces a cycle, and when the code segments in different threads that cause the cycle can actually not execute in parallel. The paper formalizes a theory for dynamic deadlock detection and compares it to model checking and static analysis techniques. It furthermore describes an implementation for analyzing Java programs and its application to two case studies: a planetary rover and a space craft altitude control system.

Bensalem, Saddek↗

Analysis of supersonic plug nozzle flowfield and heat transfer

A number of problems pertaining to the flowfield in a plug nozzle, designed as a supersonic thruster nozzle, with provision for cooling the plug with a coolant stream admitted parallel to the plug wall surface, were studied. First, an analysis was performed of the inviscid, nonturbulent, gas dynamic interaction between the primary hot stream and the secondary coolant stream. A numerical prediction code for establishing the resulting flowfield with a dividing surface between the two streams, for various combinations of stagnation and static properties of the two streams, was utilized for illustrating the nature of interactions. Secondly, skin friction coefficient, heat transfer coefficient and heat flux to the plug wall were analyzed under smooth flow conditions (without shocks or separation) for various coolant flow conditions. A numerical code was suitably modified and utilized for the determination of heat transfer parameters in a number of cases for which data are available. Thirdly, an analysis was initiated for modeling turbulence processes in transonic shock-boundary layer interaction without the appearance of flow separation.

Murthy, S. N. B.↗

Applications of a transonic wing design method

A method for designing wings and airfoils at transonic speeds using a predictor/corrector approach was developed. The procedure iterates between an aerodynamic code, which predicts the flow about a given geometry, and the design module, which compares the calculated and target pressure distributions and modifies the geometry using an algorithm that relates differences in pressure to a change in surface curvature. The modular nature of the design method makes it relatively simple to couple it to any analysis method. The iterative approach allows the design process and aerodynamic analysis to converge in parallel, significantly reducing the time required to reach a final design. Viscous and static aeroelastic effects can also be accounted for during the design or as a post-design correction. Results from several pilot design codes indicated that the method accurately reproduced pressure distributions as well as the coordinates of a given airfoil or wing by modifying an initial contour. The codes were applied to supercritical as well as conventional airfoils, forward- and aft-swept transport wings, and moderate-to-highly swept fighter wings. The design method was found to be robust and efficient, even for cases having fairly strong shocks.

Campbell, Richard L.↗

Evaluation of MHOST analysis capabilities for a plate element

Results of the evaluation of the static, buckling, and free vibration analyses capabilities of MHOST for the plate elements are presented. Two large scale, general purpose finite element codes (MARC and MSC/NASTRAN) are used to validate MHOST. Comparisons of MHOST results with those from MARC and MSC/NASTRAN show good agreement and indicate that MHOST can be used with confidence to perform the aforementioned analyses using the plate element.

Lee, Ho-Jun↗

Reliability Analysis of Brittle Material Structures - Including MEMS(?) - With the CARES/Life Program

Brittle materials are being used, or considered, for a wide variety of high tech applications that operate in harsh environments, including static and rotating turbine parts. thermal protection systems, dental prosthetics, fuel cells, oxygen transport membranes, radomes, and MEMS. Designing components to sustain repeated load without fracturing while using the minimum amount of material requires the use of a probabilistic design methodology. The CARES/Life code provides a general-purpose analysis tool that predicts the probability of failure of a ceramic component as a function of its time in service. For this presentation an interview of the CARES/Life program will be provided. Emphasis will be placed on describing the latest enhancements to the code for reliability analysis with time varying loads and temperatures (fully transient reliability analysis). Also, early efforts in investigating the validity of using Weibull statistics, the basis of the CARES/Life program, to characterize the strength of MEMS structures will be described as as well as the version of CARES/Life for MEMS (CARES/MEMS) being prepared which incorporates single crystal and edge flaw reliability analysis capability. It is hoped this talk will open a dialog for potential collaboration in the area of MEMS testing and life prediction.

Nemeth, Noel N.↗

Application of the ENS3DAE Euler/Navier-Stokes Aeroelastic method

This paper discusses recent applications of the ENS3DAE computational aeroelasticity method. In particular, it describes aeroelastic and unsteady aerodynamics calculations performed on wings with trailing edge control surfaces. These simulations include the investigation of control reversal for a structurally flexible wing with a deflected control surface, and a static and dynamic analysis of a rigid wing with an oscillating control surface. The two sets of calculations were performed independently on different wings using different grid topologies. The control reversal simulation represents an inviscid Euler static aeroelastic analysis of a thin wing with a rectangular planform. The geometry of this wing makes it suitable for computations using more approximate, inviscid aerodynamics methods. Thus, the results of the present Euler computations are compared with numerical data generated by a validated computational aeroelasticity code which uses a simpler aerodynamic formulation. The second illustrated case involves the simulation of a significantly more complex flowfield and the static and dynamic analyses of this geometry were performed using the viscous Navier-Stokes equation option in ENS3DAE. Results of both the steady and unsteady calculations on this wing are compared with existing experimental data.

Schuster, David M.↗

Uncertainty and Sensitivity Analysis of Reactivity Insertion Accidents using BISON Fuel Performance Code

The development of Accident Tolerant Fuels (ATF) has expanded into a considerable research program in the wake of the events at the Japanese Fukushima Daiichi nuclear power plant in 2011. The ATF goal is to insert a test assembly into a commercial reactor by 2022. To support this, testing of the ATF design under off-normal conditions is needed to assess the 'enhanced accident tolerance' capability of the fuel rod design over standard UO{sub 2}-Zircaloy concepts. Transient testing capabilities are necessary for off-normal testing of events such as loss of coolant accidents (LOCA) and reactivity insertion accidents (RIA). A transient testing program has been developed to support the testing of ATF designs. The Transient Reactor Test (TREAT) facility has been chosen as the site to perform the transient testing of the ATF designs. TREAT is a versatile test facility capable of subjecting fuel specimens to various off-normal conditions. Historic TREAT experiments were performed in a variety of experimental test vehicles including static environment capsules, once-through flowing systems, and also recirculating flowing systems. The planned TREAT experiments will be operated under PWR conditions. Therefore, a static capsule vehicle is under development to handle required temperatures and pressures. The static capsule, termed the Static Environment Rodlet Transient Test Apparatus (SERTTA), will be used for irradiation tests in TREAT to subject ATF fuel rodlets to a range of power pulse conditions and environments. To aid in the design and analysis of future tests in TREAT, the fuel performance code BISON is being used. Fuel performance codes capable of accident modeling have been developed for some time and they have shown an ability to reproduce experimental results with some degree of adequacy, but they still require extensive verification and validation. Future experiments in TREAT will add to a limited database of past transient tests performed in transient reactors to support LWR systems. These tests can be used for future verification and validation of fuel performance codes. Understanding the impact of the uncertainty of influential parameters is critical for not only analyzing the results of the experiment but also for modeling and validation of codes. For example, uncertainties associated with the initial manufactured state of the sample, instrument measurements, material property models, and planned experiment conditions can provide important insight into model performance, experimental data, and play an important role in experimental design. The objective of this work is to perform an uncertainty and sensitivity analysis on a simplified RIA case of a standard UO{sub 2}-Zircaloy sample using expected nominal conditions in the SERTTA capsule in TREAT. The BISON fuel performance code is used in conjunction with DAKOTA (Design Analysis Kit for Optimization and Tera-scale Applications) to perform the uncertainty and sensitivity analysis. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Mercury Model of the Molly-G Fast Burst Reactor (FBR)

In order to design neutron irradiation experiments at a given nuclear reactor, one requires a representative neutron spectrum for the configuration of the reactor that will exist during the experiment. This usually implies modeling the operation of the reactor using a neutron transport code. Since the modeling of pulsed neutron reactors is quite complex, requiring coupled modeling of neutronics and thermo-structural response, the typical methodology employs analysis of the neutron spectrum which is obtained from a static keff calculation. While this time-independent Eigenfunction is not truly representative of that during a reactor pulse, we postulate that it is a reasonable representation for our purposes. The Nuclear Survivability (NS) Program at LLNL is considering fielding such experiments at the Molly-G (‘Molybdenum Godiva-II’) Fast Burst Reactor (FBR) at the White Sands Missile Range (WSMR). The goal of these experiments is to (a) test / proof / calibrate diagnostics, and (b) test the efficacy of neutron shield configurations, before fielding at other nuclear reactors, such as the Annular Core Research Reactor (ACRR) at SNL/NM. After attempting to secure either drawings or a copy of an existing MCNP model of the Molly-G FBR from various institutions, it was decided to develop a new model for use with the Mercury Monte Carlo transport code. The model described herein will be used to design these experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Balancing Dynamic Strength of Spur Gears Operated at Extended Center Distance

This paper presents an analytical study on using hob offset to balance the dynamic tooth strength of spur gears operated at a center distance greater than the standard value. This study is an extension of a static study by Mabie and others. The study was limited to the offset values that assure the pinion and gear teeth will neither be undercut nor become pointed. The analysis presented in this paper was performed using DANST-PC, a new version of the NASA gear dynamics code. The operating speed of the transmission influences the amount of hob offset required to equalize the dynamic stresses in the pinion and gear. The optimum hob offset for the pinion was found to vary within a small range as the speed changes. The optimum value is generally greater than the optimum value found by static procedures. For gears that must operate over a wide range of speeds, an average offset value may be used.

Lin, Hsiang Hsi↗

Calculating Reuse Distance from Source Code

The efficient use of a system is of paramount importance in high-performance computing. Applications need to be engineered for future systems even before the architecture of such a system is clearly known. Static performance analysis that generates performance bounds is one way to approach the task of understanding application behavior. Performance bounds provide an upper limit on the performance of an application on a given architecture. Predicting cache hierarchy behavior and accesses to main memory is a requirement for accurate performance bounds. This work presents our static reuse distance algorithm to generate reuse distance histograms. We then use these histograms to predict cache miss rates. Experimental results for kernels studied show that the approach is accurate.

Cache memory↗

An Expert Assistant for Computer Aided Parallelization

The prototype implementation of an expert system was developed to assist the user in the computer aided parallelization process. The system interfaces to tools for automatic parallelization and performance analysis. By fusing static program structure information and dynamic performance analysis data the expert system can help the user to filter, correlate, and interpret the data gathered by the existing tools. Sections of the code that show poor performance and require further attention are rapidly identified and suggestions for improvements are presented to the user. In this paper we describe the components of the expert system and discuss its interface to the existing tools. We present a case study to demonstrate the successful use in full scale scientific applications.

Jost, Gabriele↗

An Aeroelastic Analysis of a Thin Flexible Membrane

Studies have shown that significant vehicle mass and cost savings are possible with the use of ballutes for aero-capture. Through NASA's In-Space Propulsion program, a preliminary examination of ballute sensitivity to geometry and Reynolds number was conducted, and a single-pass coupling between an aero code and a finite element solver was used to assess the static aeroelastic effects. There remain, however, a variety of open questions regarding the dynamic aeroelastic stability of membrane structures for aero-capture, with the primary challenge being the prediction of the membrane flutter onset. The purpose of this paper is to describe and begin addressing these issues. The paper includes a review of the literature associated with the structural analysis of membranes and membrane utter. Flow/structure analysis coupling and hypersonic flow solver options are also discussed. An approach is proposed for tackling this problem that starts with a relatively simple geometry and develops and evaluates analysis methods and procedures. This preliminary study considers a computationally manageable 2-dimensional problem. The membrane structural models used in the paper include a nonlinear finite-difference model for static and dynamic analysis and a NASTRAN finite element membrane model for nonlinear static and linear normal modes analysis. Both structural models are coupled with a structured compressible flow solver for static aeroelastic analysis. For dynamic aeroelastic analyses, the NASTRAN normal modes are used in the structured compressible flow solver and 3rd order piston theories were used with the finite difference membrane model to simulate utter onset. Results from the various static and dynamic aeroelastic analyses are compared.

Scott, Robert C.↗

Transient convective delayed neutron precursor of {sup 235}U for the molten salt reactor experiment

This work investigates the delayed neutron precursors from the fission of {sup 235}U in molten salt reactors. The six delayed neutron groups predict the spatial distribution of the advective transport for radioactive mass transfer in the molten-salt reactor experiment (MSRE). The Mole code, using the framework of the Multiphysics Object-Oriented Simulation Environment (MOOSE), was used to compute the delayed neutron concentration of the whole system. This approach can be used for analysis of reactivity in circulating conditions. The effect of delayed neutron concentrations on the static or dynamic behavior of the system in MSRE was analyzed. Using this approach, the pattern and validation of delayed neutron for a range of parameters of interest are reported and discussed herein. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗