Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bar 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 145 records · Page 8

PCI: A PATRAN-NASTRAN model translator

The amount of programming required to develop a PATRAN-NASTRAN translator was surprisingly small. The approach taken produced a highly flexible translator comparable with the PATNAS translator and superior to the PATCOS translator. The coding required varied from around ten lines for a shell element to around thirty for a bar element, and the time required to add a feature to the program is typically less than an hour. The use of a lookup table for element names makes the translator also applicable to other versions of NASTRAN. The saving in time as a result of using PDA's Gateway utilities was considerable. During the writing of the program it became apparent that, with a somewhat more complex structure, it would be possible to extend the element data file to contain all data required to define the translation from PATRAN to NASTRAN by mapping of data between formats. Similar data files on property, material and grid formats would produce a completely universal translator from PATRAN to any FEA program, or indeed any CAE system.

Sheerer, T. J.↗

ESTIMATION OF DISTANCES TO STARS WITH STELLAR PARAMETERS FROM LAMOST

We present a method to estimate distances to stars with spectroscopically derived stellar parameters. The technique is a Bayesian approach with likelihood estimated via comparison of measured parameters to a grid of stellar isochrones, and returns a posterior probability density function for each star’s absolute magnitude. This technique is tailored specifically to data from the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) survey. Because LAMOST obtains roughly 3000 stellar spectra simultaneously within each ∼5° diameter “plate” that is observed, we can use the stellar parameters of the observed stars to account for the stellar luminosity function and target selection effects. This removes biasing assumptions about the underlying populations, both due to predictions of the luminosity function from stellar evolution modeling, and from Galactic models of stellar populations along each line of sight. Using calibration data of stars with known distances and stellar parameters, we show that our method recovers distances for most stars within ∼20%, but with some systematic overestimation of distances to halo giants. We apply our code to the LAMOST database, and show that the current precision of LAMOST stellar parameters permits measurements of distances with ∼40% error bars. This precision should improve as the LAMOST data pipelines continue to be refined.

79 ASTRONOMY AND ASTROPHYSICS↗

Watts Bar I Ex-Core Analyses using VERA

The integration and direct coupling of the Oak Ridge National Laboratory (ORNL) Shift Monte Carlo code within the Consortium for Advanced Simulation of Light Water Reactors (CASL) Virtual Environment for Reactor Applications (VERA) offers unique capabilities that combine high-fidelity in-core and ex-core radiation transport. In this recent development activity, the deterministic neutronics code MPACT performs the incore radiation transport with temperature feedback using COBRA-TF (CTF), and provides the fission source to Shift for a follow-on fixed source radiation transport calculation to determine ex-core quantities of interest. The coupling of MPACT to Shift allows the neutron source for the ex-core transport to be calculated from the time-dependent, fully coupled MPACT solution. This paper presents the first use of VERA to perform coupon irradiation for validation of its ex-core capabilities.

Davidson, Eva↗

Watts Bar I ex-core analyses using VERA - 25430

The integration and direct coupling of the Oak Ridge National Laboratory (ORNL) Shift Monte Carlo code within the Consortium for Advanced Simulation of Light Water Reactors (CASL) Virtual Environment for Reactor Applications (VERA) offers unique capabilities that combine high-fidelity in-core and ex-core radiation transport. In this recent development activity, the deterministic neutronics code MPACT performs the in-core radiation transport with temperature feedback using COBRA-TF (CTF), and provides the fission source to Shift for a follow-on fixed source radiation transport calculation to determine ex-core quantities of interest. The coupling of MPACT to Shift allows the neutron source for the ex-core transport to be calculated from the time-dependent, fully coupled MPACT solution. This paper presents the first use of VERA to perform coupon irradiation for code-to-code comparison of its ex-core capabilities. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Opus radiative opacity code, a new implementation of the super-transition array method

This report describes the development of our new radiative opacity code Opus, submitted as part of a 2020 Level 2 milestone requirement. This code aims to incorporate the capabilities in existing but separate opacity codes previously developed at LLNL (Opal, Topaz, and Tycho). The first phase in the development, which was to implement the super-transition (STA) method of Bar-Shalom and co-workers is now completed. In this work, we describe the underlying physics, implementation, usage, and regression testing as well as providing numerical examples for the opacity tables of boron, chromium, iron, and nickel. The latter are benchmarked against the LLNL legacy STA code Tycho, with an average absolute relative difference of ~1.5 and 0.5% for Planck and Rosseland opacity means, respectively. The cases where differences are larger are well understood in terms of numerical instabilities near discontinuities of the average degree of ionization inherent to the assumed opacity/equation-of-state model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The dynamics of stellar discs in live dark-matter haloes

Recent developments in computer hardware and software enable researchers to simulate the self-gravitating evolution of galaxies at a resolution comparable to the actual number of stars. Here we present the results of a series of such simulations. In this study, we performed N-body simulations of disc galaxies with between 100 and 500 million particles over a wide range of initial conditions. Our calculations include a live bulge, disc, and dark-matter halo, each of which is represented by self-gravitating particles in the N-body code. The simulations are performed using the gravitational N-body tree-code BONSAI running on the Piz Daint supercomputer. We find that the time-scale over which the bar forms increases exponentially with decreasing disc-mass fraction and that the bar formation epoch exceeds a Hubble time when the disc-mass fraction is ~0.35. These results can be explained with the swing-amplification theory. The condition for the formation of m = 2 spirals is consistent with that for the formation of the bar, which is also an m = 2 phenomenon. We further argue that the non-barred grand-design spiral galaxies are transitional, and that they evolve to barred galaxies on a dynamical time-scale. We also confirm that the disc-mass fraction and shear rate are important parameters for the morphology of disc galaxies. The former affects the number of spiral arms and the bar formation epoch, and the latter determines the pitch angle of the spiral arms.

79 ASTRONOMY AND ASTROPHYSICS↗

ARCANE Reweighting Demo

ARCANE reweighting is a Monte Carlo technique for tackling the negative weights problem in simulations in high energy physics. The technique works by redistributing (via an additive reweighting) the contributions of different pathways within an event generator that lead to the same final event. The technique is exact and does not introduce any biases in the distributions of physical observables. This repository contains the code and data associated with a demonstration of ARCANE reweighting, for the generation of $e^+ e^- \longrightarrow q\bar{q} + 1\,\mathrm{jet}$ events at next-to-leading-order accuracy in perturbation theory. ARCANE reweighting was shown to almost completely eliminate the negative weights problem for this process.

Shyamsundar, Prasanth [Fermi National Accelerator ↗

N-BAR PROBLEMS AS APPROXIMATIONS TO THE BREE PROBLEM

The Bree solution to the problem of a ratcheting cylinder under constant pressure and cyclic thermal load is a fundamental result in nuclear engineering and widely used as the technical basis for the ASME Boiler and Pressure Vessel Code and other design methods. However, because the loading conditions in the Bree problem are difficult to achieve experimentally there have been relatively few works experimentally examining the problem and extending it to other relevant design situations, for example cladded components. In contrast, 2-bar problems are widely studied experimentally and are relatively easy to setup. These 2-bar problems are thought to be representative of Bree-type geometries, but a formal connection has not been demonstrated. This work formally establishes the connection between the Bree cylinder and an n-bar problem – a coupled bar experiment with, in general, more than two bars linked in parallel. The connection suggests that n-bar experiments using a fairly limited number of bars might be an experimentally-accessible setup that better represents ratcheting phenomenon in actual nuclear pressurized components. Such experiments could test surrogate cladded or multi-material components by using bars of different materials. Finally, this work suggests control schemes that yield optimally efficient n-bar experiments – experiments that best replicates a Bree cylinder with a limited number of bars.

Messner, Mark C.↗

{ital n}-bar problems as approximations to the Bree problem

The Bree solution to the problem of a ratcheting cylinder under constant pressure and cyclic thermal load is a fundamental result in nuclear engineering and widely used as the technical basis for the ASME Boiler and Pressure Vessel Code and other design methods. However, because the loading conditions in the Bree problem are difficult to achieve experimentally there have been relatively few works experimentally examining the problem and extending it to other relevant design situations, for example cladded components. In contrast, 2-bar problems are widely studied experimentally and are relatively easy to setup. These 2-bar problems are thought to be representative of Bree-type geometries, but a formal connection has not been demonstrated. This work formally establishes the connection between the Bree cylinder and an n-bar problem a coupled bar experiment with, in general, more than two bars linked in parallel. The connection suggests that n-bar experiments using a fairly limited number of bars might be an experimentally-accessible setup that better represents ratcheting phenomenon in actual nuclear pressurized components. Such experiments could test surrogate cladded or multi-material components by using bars of different materials. Finally, this work suggests control schemes that yield optimally efficient n-bar experiments experiments that best replicates a Bree cylinder with a limited number of bars.

Messner, M. C.↗

Ensemble Monte Carlo calculations with five novel moves

We introduce five novel types of Monte Carlo (MC) moves that brings the number of moves of ensemble MC calculations from three to eight. So far such calculations have relied on affine invariant stretch moves that were originally introduced by Christen (2007), walk moves by Goodman and Weare (2010) and quadratic moves by Militzer (2023). Ensemble MC methods have been very popular because they harness information about the fitness landscape from a population of walkers rather than relying on expert knowledge. Here we modified the affine method and employed a simplex of points to set the stretch direction. We adopt the simplex concept to quadratic moves. We also generalize quadratic moves to arbitrary order. Finally, we introduce directed moves that employ the values of the probability density while all other types of moves rely solely on the location of the walkers. We apply all algorithms to the Rosenbrock density in 2 and 20 dimensions and to the ring potential in 12 and 24 dimensions. We evaluate their efficiency by comparing error bars, autocorrelation time, travel time, and the level of cohesion that measures whether any walkers were left behind. Our code is open source.

97 MATHEMATICS AND COMPUTING↗

Identifying microturbulence regimes in a TCV discharge making use of physical constraints on particle and heat fluxes

Reducing the uncertainty on physical input parameters derived from experimental measurements is essential towards improving the reliability of gyrokinetic turbulence simulations. This can be achieved by introducing physical constraints. Amongst them, the zero particle flux condition is considered here. A first attempt is also made to match as well the experimental ion/electron heat flux ratio. This procedure is applied to the analysis of a particular Tokamak à Configuration Variable discharge. A detailed reconstruction of the zero particle flux hyper-surface in the multi-dimensional physical parameter space at fixed time of the discharge is presented, including the effect of carbon as the main impurity. Both collisionless and collisional regimes are considered. Hyper-surface points within the experimental error bars are found. In conclusion, the analysis is done performing gyrokinetic simulations with the local version of the GENE code, computing the fluxes with a Quasi-Linear (QL) model and validating the QL results with non-linear simulations in a subset of cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance of an Array of 6 X 6 X 20 mm3 Virtual Frisch-Grid CdZnTe Detectors With A Waveform Sampling Readout System

Arrays of CdZnTe (CZT) bars with high spatial and energy resolution are an enabling technology for future space missions, such as the image plane detector of the Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO). GECCO will explore the medium-energy (100 keV – 10 MeV) gamma-ray sky, which is among the least covered windows to the Universe with a massive potential for discovery. We have integrated an array of 6x6x20 mm3 virtual Frisch-grid (VFG) CZT detectors, measuring signals from the cathodes, anodes, and pads on each of the four sides of each bar. The VFG CZT detector is read out with an IDEAS GDS-100 waveform sampling readout system to optimize the signal processing and event reconstruction. We evaluated the energy and position resolution of the integrated VFG CZT array with a 137Cs source and successfully demonstrated better than 1% energy resolution and finer than 1 mm position resolution at 662 keV. Furthermore, we carried out a beam test at the High Intensity Gamma-ray Source (HIGS) at Duke University to investigate the performance of the system for higher energy gamma-rays. Together with the basic performance obtained from the 137Cs source, we present here those results from the beam test with various gamma-ray energies up to 8 MeV.

Makoto Sasaki↗

Design and Modeling of an 80 bar Oxy-Combustor for Direct Fired Supercritical CO2 Applications

As supercritical CO2 power cycles for fossil energy power generation continue to generate interest, there are significant issues and unanswered questions regarding injector design, flame stabilization, wall heat transfer, CO emissions, combustion dynamics and other combustion phenomenon. For natural gas, direct-fired cycles with carbon capture it is believed that Computational Fluid Dynamics (CFD) modeling will play an essential role in the combustor design process. To accurately model turbulent reacting flows at these unique conditions, experimental data is needed to validate CFD codes and sub-models and is currently lacking at conditions relevant for these cycles. This paper presents the conceptual design and CFD simulations of an experimental 80 bar oxy-combustion facility and test article currently under construction at NETL. The facility is targeted towards the testing of a single injector, direct-fired sCO2 combustor at the 100 kW thermal output level. While these conditions do not reflect the actual Allam cycle operating conditions (300 bar) they are viewed as a stepping stone in the model validation process at supercritical conditions. Reynolds Averaged Navier-Stokes as well as Large Eddy Simulation are used to model the turbulent combustion process and aid in the design of the combustor and injector. Process parameters including oxidizer preheat temperature and combustor flowrate are investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The coupling between flame surface dynamics and species mass conservation in premixed turbulent combustion

Current flamelot models based on a description of the flame surface dynamics require the closure of two inter-related equations: a transport equation for the mean reaction progress variable, (tilde)c, and a transport equation for the flame surface density, Sigma. The coupling between these two equations is investigated using direct numerical simulations (DNS) with emphasis on the correlation between the turbulent fluxes of (tilde)c, bar(pu''c''), and Sigma, (u'')(sub S)Sigma. Two different DNS databases are used in the present work: a database developed at CTR by A. Trouve and a database developed by C. J. Rutland using a different code. Both databases correspond to statistically one-dimensional premixed flames in isotropic turbulent flow. The run parameters, however, are significantly different, and the two databases correspond to different combustion regimes. It is found that in all simulated flames, the correlation between bar(pu''c'') and (u'')(sub S)Sigma is always strong. The sign, however, of the turbulent flux of (tilde)c or Sigma with respect to the mean gradients, delta(tilde)c/delta(x) or delta(Sigma)/delta(x), is case-dependent. The CTR database is found to exhibit gradient turbulent transport of (tilde)c and Sigma, whereas the Rutland DNS features counter-gradient diffusion. The two databases are analyzed and compared using various tools (a local analysis of the flow field near the flame, a classical analysis of the conservation equation for (tilde)(u''c''), and a thin flame theoretical analysis). A mechanism is then proposed to explain the discrepancies between the two databases and a preliminary simple criterion is derived to predict the occurrence of gradient/counter-gradient turbulent diffusion.

Trouve, A.↗

Boundary Layer Transition in the Leading Edge Region of a Swept Cylinder in High Speed Flow

Experiments were conducted on a 76 degree swept cylinder to establish the behavior of the attachment line transition process in a low-disturbance level, Mach number 1.6 flow. For a near adiabatic wall condition, the attachment-line boundary layer remained laminar up to the highest attainable Reynolds number. The attachment-line boundary layer transition under the influence of trip wires depended on wind tunnel disturbance level, and a transition onset condition for this flow is established. Internal heating raised the surface temperature of the attachment line to induce boundary layer instabilities. This was demonstrated experimentally for the first time and the frequencies of the most amplified disturbances were determined over a range of temperature settings. Results were in excellent agreement to those predicted by a linear stability code, and provide the first experimental verification of theory. Transition onset along the heated attachment line at an R-bar of 800 under quiet tunnel conditions was found to correlate with an N factor of 13.2. Increased tunnel disturbance levels caused the transition onset to occur at lower cylinder surface temperatures and was found to correlate with an approximate N factor of 1 1.9, so demonstrating that the attachment-line boundary layer is receptive to increases in the tunnel disturbance level.

Coleman, Colin P.↗

Report on an Assessment of the Application of EPP Results from the Strain Limit Evaluation Procedure to the Prediction of Cyclic Life Based on the SMT Methodology

The goal of the proposed integrated Elastic Perfectly-Plastic (EPP) and Simplified Model Test (SMT) methodology is to incorporate an SMT data based approach for creep-fatigue damage evaluation into the EPP methodology to avoid the separate evaluation of creep and fatigue damage and eliminate the requirement for stress classification in current methods; thus greatly simplifying evaluation of elevated temperature cyclic service. This methodology should minimize over-conservatism while properly accounting for localized defects and stress risers. To support the implementation of the proposed methodology and to verify the applicability of the code rules, analytical studies and evaluation of thermomechanical test results continued in FY17. This report presents the results of those studies. An EPP strain limits methodology assessment was based on recent two-bar thermal ratcheting test results on 316H stainless steel in the temperature range of 405 to 7050C. Strain range predictions from the EPP evaluation of the two-bar tests were also evaluated and compared with the experimental results. The role of sustained primary loading on cyclic life was assessed using the results of pressurized SMT data from tests on Alloy 617 at 9500C. A viscoplastic material model was used in an analytic simulation of two-bar tests to compare with EPP strain limits assessments using isochronous stress strain curves that are consistent with the viscoplastic material model. A finite element model of a prior 304H stainless steel Oak Ridge National Laboratory (ORNL) nozzle-to-sphere test was developed and used for an EPP strain limits and creep-fatigue code case damage evaluations. A theoretical treatment of a recurring issue with convergence criteria for plastic shakedown illustrated the role of computer machine precision in EPP calculations.

36 MATERIALS SCIENCE↗

Finite element analysis of notch behavior using a state variable constitutive equation

The state variable constitutive equation of Bodner and Partom was used to calculate the load-strain response of Inconel 718 at 649 C in the root of a notch. The constitutive equation was used with the Bodner-Partom evolution equation and with a second evolution equation that was derived from a potential function of the stress and state variable. Data used in determining constants for the constitutive models was from one-dimensional smooth bar tests. The response was calculated for a plane stress condition at the root of the notch with a finite element code using constant strain triangular elements. Results from both evolution equations compared favorably with the observed experimental response. The accuracy and efficiency of the finite element calculations also compared favorably to existing methods.

Dame, L. T.↗

Fifteenth NASTRAN (R) Users' Colloquium

Numerous applications of the NASA Structural Analysis (NASTRAN) computer program, a general purpose finite element code, are discussed. Additional features that can be added to NASTRAN, interactive plotting of NASTRAN data on microcomputers, mass modeling for bars, the design of wind tunnel models, the analysis of ship structures subjected to underwater explosions, and buckling analysis of radio antennas are among the topics discussed.

Source record↗