Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “G 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 73 records · Page 4

Fatigue Life Analysis of Tapered Hybrid Composite Flexbeams

Nonlinear-tapered flexbeam laminates from a full-size composite helicopter rotor hub flexbeam were tested under combined constant axial tension and cyclic bending loads. The two different graphite/glass hybrid configurations tested under cyclic loading failed by delamination in the tapered region. A 2-D finite element model was developed which closely approximated the flexbeam geometry, boundary conditions, and loading. The analysis results from two geometrically nonlinear finite element codes, ANSYS and ABAQUS, are presented and compared. Strain energy release rates (G) obtained from the above codes using the virtual crack closure technique (VCCT) at a resin crack location in the flexbeams are presented for both hybrid material types. These results compare well with each other and suggest that the initial delamination growth from the resin crack toward the thick region of the flexbeam is strongly mode II. The peak calculated G values were used with material characterization data to calculate fatigue life curves and compared with test data. A curve relating maximum surface strain to number of loading cycles at delamination onset compared reasonably well with the test results.

Murri, Gretchen B.↗

Computational modeling of jet induced mixing of cryogenic propellants in low-G

The SOLA-ECLIPSE Code is being developed to enable computational prediction of jet induced mixing in cryogenic propellant tanks in a low-gravity environment. Velocity fields, predicted for scale model tanks, are presented which compare favorably with the available experimental data. A full scale liquid hydrogen tank for a typical Orbit Transfer Vehicle is analyzed with the conclusion that coupling an axial mixing jet with a thermodynamic vent system appears to be a viable concept for the control of tank pressure.

Hochstein, J. I.↗

Computational modeling of jet induced mixing of cryogenic propellants in low-G

The SOLA-ECLIPSE Code is being developed to enable computational prediction of jet induced mixing in cryogenic propellant tanks in a low-gravity environment. Velocity fields, predicted for scale model tanks, are presented which compare favorably with the available experimental data. A full scale liquid hydrogen tank for a typical Orbit Transfer Vehicle is analyzed with the conclusion that coupling an axial mixing jet with a thermodynamic vent system appears tobe a viable concept for the control of tank pressure. Previously announced in STAR as N84-25000

Hochstein, J. I.↗

Plasma rotation and diamagnetic drift effects on the resistive wall modes in the negative triangularity tokamaks

Abstract It was found previously that the negative triangularity (NT) configuration is more MHD-unstable for low n modes than the positive triangularity (PT) case, although the situation is reversed for intermediate n modes and the NT configuration becomes more stable for intermediate n modes ( n = 3 − 10 ) (Zheng et al 2021 Nucl. Fusion 61 116014). Here, n is the toroidal mode number. In this work, we extend the studies to include the rotation effects, as well as the diamagnetic drift effects, to see how the resistive wall modes (RWMs) in the NT configuration are affected as compared with the PT configuration. This is particularly motivated by noting that the wall interface with the plasma is quite different between the NT and PT configurations. It affects the plasma rotation and diamagnetic drift effects on the low n RWM. We consider the DIII-D-NT-experiment equilibrium reconstructed by the EFIT code. Based on the equilibrium g-file, the extended equilibria are constructed with the VMEC code by varying the beta values while keeping the pressure and poloidal current flux profiles basically unchanged. The bootstrap current contribution to the equilibria is taken into account with the Sauter formula. The MHD stability is then computed using the AEGIS code with the rotation and diamagnetic drift effects taken into account. We found that, although the NT configuration is less stable for n = 1 MHD modes, the rotation and diamagnetic drift stabilization effects on RWMs are more effective in the NT configuration than in the PT one. Note that even in the PT case, the stabilization of RWMs by the rotation and kinetic effects is critical. Because the low-n RWMs in the regular NT case are more unstable, the rotation and diamagnetic drift stabilization effects found in this research are important for the NT tokamak concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NuScale Pressure and Temperature Limits Methodology Using Finite Element Analysis

Per 10 CFR 50 Appendix G, the pressure-temperature (P-T) limits curves and minimum temperature must be established to provide adequate margins for ferritic pressure-retaining components of the reactor coolant pressure boundary; this is to protect against brittle failure during any normal operating conditions, including anticipated operational occurrences and system hydrostatic tests, to which the pressure boundary may be subjected over its service lifetime. Specifically, ASME Code Section XI Appendix G procedures must be used for P-T limits calculation considering the pressure and temperature at various operating transient conditions. However, the elastic fracture mechanics solutions in Section XI are only suitable for cylindrical reactor pressure vessel (RPV) beltline without geometric discontinuities. Hence, these solutions are not suitable for postulated flaws near the core support blocks attached to the NuScale RPV inside surface, which is part of the beltline. As a result, NuScale has used finite element analysis to calculate thermal stress and stress intensity factor for the postulated flaws. The solutions using finite-element analysis have been validated using the formulations for straight cylinders for both axial and circumferential flaws. In addition, special considerations are given to the RPV beltline nil-ductility transition temperature (RTNDT) due to neutron irradiation occurring at lower temperatures than conventional plants. The paper summarizes the methodology and finite-element models used to develop P-T limits curves for NuScale RPV at the end of its 60-year design life.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hand Calculation Methods for Nuclear Criticality Safety

This primer provides an overview of the most common hand calculation methods used for criticality safety calculations. The most widely used tools available to a nuclear criticality safety (NCS) practitioner are probably the common Monte Carlo or deterministic criticality safety codes, which can be used to model very complex systems. However, use of these codes can obscure the parameters to which a particular fissile system may be sensitive, whereas the hand calculation methods can be used to delve into the ways each parameter may affect the reactivity of a fissile material system. Furthermore, practitioners must avoid using computer codes as devices that take inputs and simply provide outputs (i.e., a “black box”). Many years ago, pioneers such as Joe Thomas, David Smith, and Hugh Paxton, among others in the field of nuclear criticality safety, took the time before the advent of high-speed desktop computers to create simple hand methods for criticality safety analyses. Some of the methods can be used for single fissile units; others are applicable to fissile units arranged into simple array configurations. This primer discusses the applicability of the various methods, illustrates how they are used, and provides an interpretation of the various results. The NCS practitioner will need to spend time to master the methods that could be most useful; however, they can provide the practitioner with fast and accurate answers to criticality safety problems if they are used correctly and if critical data exist for the problem at hand. Hand calculation methods can be used as a starting point for more advanced calculations, and in many circumstances, they can provide sensitivity and perturbation information more quickly than using a criticality code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ultraviolet photometry of globular clusters

Filter photometry measurements of the integrated light from several globular clusters between 1500 and 4000 A are discussed. The techniques and problems associated with carrying out photometry of faint extended sources with the OAO-2 instrumentation are described first. The energy distribution for globular clusters determined from these data are corrected for interstellar extinction on the basis of B-V color excesses and a mean ultraviolet extinction curve. The results are compared with synthetic cluster calculations. The ultraviolet energy distribution is primarily determined by the number of stars at the blue end of the horizontal branch.

Welch, G. A.↗

Ultraviolet spectropolarimetry of the Be star PP Carinae with the Wisconsin Ultraviolet Photo-Polarimeter Experiment

We present the first ultraviolet spectropolarimetric observations of the Be star PP Car, obtained with the Wisconsin Ultraviolet Photo-Polarimeter Experiment (WUPPE) aboard the Astro 1 mission. Usable polarization data were obtained from 1400 to 2330 A, along with a good spectrum from 1400 to 3200 A. These data show a lower polarization shortward of the Balmer jump than had been predicted by standard models, and a broad UV polarization dip around 1900 A is seen. These results are in agreement with those found from the WUPPE observations of two other Be stars, Xi Tau and Pi Aqr, which were published earlier. All these observations are an important probe of the Be circumstellar envelopes and demonstrate the need for the inclusion of metal-line effects in circumstellar disk models of Be star UV polarization.

Bjorkman, K. S.↗

Paving the Way for Net Zero Energy Codes through Performance Based Approaches

The prescriptive path is the most widely used approach for commercial code compliance in the United States. However, to achieve near-term net zero energy performance that many stakeholders in the buildings industry desire, it is clear that energy codes will need to transition from prescriptive to performance-based approaches. However, that transition is not without potential pitfalls. Among the concerns expressed by stakeholders are added complexity and cost, lack of confidence in energy modeling results, gamesmanship on the part of applicants, lack or qualified reviewers, challenges in promoting de-carbonization, and the inequity of trading long lived envelope efficiency for short lived measures such as building controls. This paper will discuss multiple approaches being evaluated and implemented by Pacific Northwest National Laboratory to improve the usability of the performance-based approach for both code compliance and beyond code programs to meet aggressive energy savings targets. A 'system performance' approach provides a simpler approach compared to whole building performance, while keeping tradeoffs limited to specific building systems. A simplified whole-building performance-based approach for small commercial buildings can cost effectively achieve deeper savings for buildings that typically follow the prescriptive compliance approach. Improved reporting and verification processes applied to the traditional performance path can provide greater confidence in simulation results and facilitate adoption of performance-based approaches. Prescriptive packages provide prescriptive solutions while targeting specific performance thresholds. This paper discusses the progress being made with each of these approaches and provides examples of their implementation in energy codes and beyond code programs.

Performance Based Codes, Appendix G, PRM, TSPR, Sy↗

Simulation of Combustion Systems with Realistic g-jitter

In this project a transient, fully three-dimensional computer simulation code was developed to simulate the effects of realistic g-jitter on a number of combustion systems. The simulation code is capable of simulating flame spread on a solid and nonpremixed or premixed gaseous combustion in nonturbulent flow with simple combustion models. Simple combustion models were used to preserve computational efficiency since this is meant to be an engineering code. Also, the use of sophisticated turbulence models was not pursued (a simple Smagorinsky type model can be implemented if deemed appropriate) because if flow velocities are large enough for turbulence to develop in a reduced gravity combustion scenario it is unlikely that g-jitter disturbances (in NASA's reduced gravity facilities) will play an important role in the flame dynamics. Acceleration disturbances of realistic orientation, magnitude, and time dependence can be easily included in the simulation. The simulation algorithm was based on techniques used in an existing large eddy simulation code which has successfully simulated fire dynamics in complex domains. A series of simulations with measured and predicted acceleration disturbances on the International Space Station (ISS) are presented. The results of this series of simulations suggested a passive isolation system and appropriate scheduling of crew activity would provide a sufficiently "quiet" acceleration environment for spherical diffusion flames.

Mell, William E.↗

The LLNL nuclear data infrastructure for the GNDS data format

The next generation of nuclear data infrastructure tools at the Livermore National Laboratory (LLNL) consists of pipeline of codes that read and process nuclear data from evaluated files saved in the new GNDS (Generalised Nuclear Data Structure) nuclear data format. The processing code FUDGE (For Updating Data and Generating Evaluations) is at the front-end of this pipeline as it reads and process the evaluated data for use in downstream transport codes. FUDGE is Python based with C and C++ extensions for computationally intensive tasks. As is the case for the evaluated data, the processed output is also saved in the GNDS format and the GIDI+ API is provided as the interface between the processed data and the transport codes. GIDI+ is a C++ based suite of codes and it includes GIDI (General Interaction Data Interface), a library for reading and writing GNDS data, and MCGIDI which is the cross section lookup, and reaction and product distribution sampling interface between Monte Carlo transport codes and the GNDS data. GIDI provides methods for easy access to the multi-group processed GNDS data and this is demonstrated through its implementation in ARDRA, the LLNL deterministic transport code. The evaluation and sampling methods in MCGIDI are available as both CPU and GPU methods which facilitates the use of MCGIDI in both traditional CPU-based as well as the next generation mixed model computational architectures. This is demonstrated through the GIDI+ implementation in MERCURY, the LLNL Monte Carlo transport code. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Sidelobe reduction of Barker codes, part 7.1C

Barker codes can be realized as simple digital preprocessors, or with the use of a delay line, can be realized as an analogue device. Assuming the process observed has an autocorrelation time that is long compared to the code length, a gain in signal to noise, G = N where N is the number of elements in the code, may be realized. Typically Barker codes are implemented by shifting the phase of the transmitted signal by 0 or 180 degrees according to the code pattern. A typical detection scheme is shown where the output of the radar receiver is delayed by an analogue delay line or a shift register. The signal in each element of the delay line or shift register is continuously multiplied by the code and all elements are summed. A principal disadvantage was the sidelobe response. In the case of MST or ST radar echoes where there is a dynamic range as large as 60 dB, these unwanted sidelobes cause an increase in the apparent width of atmospheric layers, or some narrow layer to appear at several altitudes. The sidelobes of a Barker code may be reduced by tapering the response of the decoder. Optimum tapers, are shown. It is found that the sidelobe response is greatly reduced with very little loss in G.

Green, J. L.↗

Damage and Failure in a Statistical Crack Model

In this paper I will take a close look at a statistical crack model (SCM) as is used in engineering computer codes to simulate fracture at high strain rates. My general goal is to understand the macroscopic behavior effected by the microphysical processes incorporated into an SCM. More specifically, I will assess the importance of including local interactions between cracks into the growth laws of an SCM. My strategy will be to construct a numerical laboratory that represents a single computational cell containing a realization of a statistical distribution of cracks. The cracks will evolve by the microphysical models of the SCM, leading to quantifiable damage and failure of the computational cell. I will use the numerical data generated by randomly generated ensembles of the fracture process to establish scaling laws that will modify and simplify the implementation of the SCM into large scale engineering codes.

42 ENGINEERING↗

A connection between block and convolutional codes

Convolutional codes of any rate and any constraint length give rise to a sequence of quasi-cyclic codes. Conversely, any quasi-cyclic code may be convolutionally encoded. Among the quasi-cyclic codes are the quadratic residue codes, Reed-Solomon codes and optimal BCH codes. The constraint length K for the convolutional encoding of many of these codes (Golay, (48, 24) QR, etc.) turns out to be surprisingly small. Thus using the soft decoding techniques for convolutional decoding, a new maximum likelihood decoding algorithm for many block codes is established. Conversely an optimal quasi-cyclic code will yield a convolutional encoding with optimal local properties and therefore with good infinite convolutional coding properties.

Solomon, G.↗

A Computer Code for Swirling Turbulent Axisymmetric Recirculating Flows in Practical Isothermal Combustor Geometries

A primitive pressure-velocity variable finite difference computer code was developed to predict swirling recirculating inert turbulent flows in axisymmetric combustors in general, and for application to a specific idealized combustion chamber with sudden or gradual expansion. The technique involves a staggered grid system for axial and radial velocities, a line relaxation procedure for efficient solution of the equations, a two-equation k-epsilon turbulence model, a stairstep boundary representation of the expansion flow, and realistic accommodation of swirl effects. A user's manual, dealing with the computational problem, showing how the mathematical basis and computational scheme may be translated into a computer program is presented. A flow chart, FORTRAN IV listing, notes about various subroutines and a user's guide are supplied as an aid to prospective users of the code.

Lilley, D. G.↗

Portable EDITOR (PEDITOR): A portable image processing system

The PEDITOR image processing system was created to be readily transferable from one type of computer system to another. While nearly identical in function and operation to its predecessor, EDITOR, PEDITOR employs additional techniques which greatly enhance its portability. These cover system structure and processing. In order to confirm the portability of the software system, two different types of computer systems running greatly differing operating systems were used as target machines. A DEC-20 computer running the TOPS-20 operating system and using a Pascal Compiler was utilized for initial code development. The remaining programmers used a Motorola Corporation 68000-based Forward Technology FT-3000 supermicrocomputer running the UNIX-based XENIX operating system and using the Silicon Valley Software Pascal compiler and the XENIX C compiler for their initial code development.

Angelici, G.↗