Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code of Accounts”

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,099 records · Page 61

Critical speeds and forced response solutions for active magnetic bearing turbomachinery, part 2

The need for better performance of turbomachinery with active magnetic bearings has necessitated a study of such systems for accurate prediction of their vibrational characteristics. A modification of existing transfer matrix methods for rotor analysis is presented to predict the response of rotor systems with active magnetic bearings. The position of the magnetic bearing sensors is taken into account and the effect of changing sensor position on the vibrational characteristics of the rotor system is studied. The modified algorithm is validated using a simpler Jeffcott model described previously. The effect of changing from a rotating unbalance excitation to a constant excitation in a single plane is also studied. A typical eight stage centrifugal compressor rotor is analyzed using the modified transfer matrix code. The results for a two mass Jeffcott model were presented previously. The results obtained by running this model with the transfer matrix method were compared with the results of the Jeffcott analysis for the purposes of verification. Also included are plots of amplitude versus frequency for the eight stage centrifugal compressor rotor. These plots demonstrate the significant influence that sensor location has on the amplitude and critical frequencies of the rotor system.

Rawal, D.↗

Determination of cosmic bodies size-velocity distribution by observation of current impacts on Mars

Collisions of cosmic bodies with terrestrial planets involve many physical processes such as deceleration and ablation during their flight through an atmosphere, the impact at a surface accompanied by cratering, melting and evaporation of surface material, generation of shock waves, etc. If body velocity is high enough then a thermal radiation is very important. All these processes on Mars proceed differently than on the other planets because of the low density of its atmosphere. In particular, this leads to the fact that smaller bodies of sizes of the order of 0.1-10 m strike the planet surface without being decelerated and perform some effects which may be detected by equipment placed on a board of artificial satellites, by a network of stations at the surface of Mars and even from the Earth. These observations can be used to determine size-velocity distribution of such bodies in the Solar System. Numerical simulation of the impacts at the surface of Mars have been carried out using two-dimensional gas dynamic code with detailed consideration of the thermal radiative transfer. This work is an extension of our previous paper. We have expanded a range of projectile sizes up to r sub 0 = 100 m. For such a large-scale body, the initial stage of the impact, involving crate ring and ejection of surface material, is very important. Thus, these effects have been taken into account.

Nemchinov, I. V.↗

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.↗

Viscous-flow analysis of a subsonic transport aircraft high-lift system and correlation with flight data

High-lift system aerodynamics has been gaining attention in recent years. In an effort to improve aircraft performance, comprehensive studies of multi-element airfoil systems are being undertaken in wind-tunnel and flight experiments. Recent developments in Computational Fluid Dynamics (CFD) offer a relatively inexpensive alternative for studying complex viscous flows by numerically solving the Navier-Stokes (N-S) equations. Current limitations in computer resources restrict practical high-lift N-S computations to two dimensions, but CFD predictions can yield tremendous insight into flow structure, interactions between airfoil elements, and effects of changes in airfoil geometry or free-stream conditions. These codes are very accurate when compared to strictly 2D data provided by wind-tunnel testing, as will be shown here. Yet, additional challenges must be faced in the analysis of a production aircraft wing section, such as that of the NASA Langley Transport Systems Research Vehicle (TSRV). A primary issue is the sweep theory used to correlate 2D predictions with 3D flight results, accounting for sweep, taper, and finite wing effects. Other computational issues addressed here include the effects of surface roughness of the geometry, cove shape modeling, grid topology, and transition specification. The sensitivity of the flow to changing free-stream conditions is investigated. In addition, the effects of Gurney flaps on the aerodynamic characteristics of the airfoil system are predicted.

Potter, R. C.↗

An Improved Analytic Model for Microdosimeter Response

An analytic model used to predict energy deposition fluctuations in a microvolume by ions through direct events is improved to include indirect delta ray events. The new model can now account for the increase in flux at low lineal energy when the ions are of very high energy. Good agreement is obtained between the calculated results and available data for laboratory ion beams. Comparison of GCR (galactic cosmic ray) flux between Shuttle TEPC (tissue equivalent proportional counter) flight data and current calculations draws a different assessment of developmental work required for the GCR transport code (HZETRN) than previously concluded.

Shinn, Judy L.↗

Numerical Study Comparing RANS and LES Approaches on a Circulation Control Airfoil

A numerical study over a nominally two-dimensional circulation control airfoil is performed using a large-eddy simulation code and two Reynolds-averaged Navier-Stokes codes. Different Coanda jet blowing conditions are investigated. In addition to investigating the influence of grid density, a comparison is made between incompressible and compressible flow solvers. The incompressible equations are found to yield negligible differences from the compressible equations up to at least a jet exit Mach number of 0.64. The effects of different turbulence models are also studied. Models that do not account for streamline curvature effects tend to predict jet separation from the Coanda surface too late, and can produce non-physical solutions at high blowing rates. Three different turbulence models that account for streamline curvature are compared with each other and with large eddy simulation solutions. All three models are found to predict the Coanda jet separation location reasonably well, but one of the models predicts specific flow field details near the Coanda surface prior to separation much better than the other two. All Reynolds-averaged Navier-Stokes computations produce higher circulation than large eddy simulation computations, with different stagnation point location and greater flow acceleration around the nose onto the upper surface. The precise reasons for the higher circulation are not clear, although it is not solely a function of predicting the jet separation location correctly.

Rumsey, Christopher L.↗

A Visual Atlas on Strategic Trade

Strategic trade concerns global trade flows of items with strategic military value. Items include goods with exclusive military use, but also dual-use items that can have both military and civil uses. Many countries require licenses for trade in strategic items as a measure of non-proliferation. For example, items listed for trade control by multilateral export control arrangements, such as the Nuclear Suppliers Group (NSG), typically are considered strategic goods. Analysts and officials involved in strategic trade control often lack quantitative global understanding of relevant trade patterns that would help them focus control and enforcement efforts, select companies for outreach and audit, target transactions for verification and inspection, and assess the potential economic impact of trade control policies. Many are unaware of available open source data on international trade that can inform these issues. Many do not know how to make use of such data or lack the tools necessary to do so. A book of charts dubbed a ‘Strategic Trade Atlas’ was developed to illustrate global trade flows of strategic goods and raise awareness of the availability and utility of open source trade data. The Atlas was prepared for the World Customs Organization (WCO) and its Member States to support their planning and execution of Strategic Trade Control Enforcement (STCE) operations. It offers macroscopic graphical representations of global trade flows classified under Harmonized System (HS) codes associated by the WCO with commodities described in its STCE Implementation Guide. While these HS codes are understood to also encompass many non-strategic items, they nevertheless provide a good overview of potentially strategic trade flows. These include HS codes covering nuclear and nuclear-related commodities listed by the NSG, such as uranium, heavy water, graphite, carbon fiber, nickel powder, beryllium, aluminum tubes, vacuum pumps, valves, machine tools, cameras, and many others. The Atlas uses trade data from the Base pour l’Analyse du Commerce International (BACI), a statistical elaboration of public UN Comtrade data. BACI data reconciles trade flows, taking into account asymmetries in import and export trade valuation. The data visualization approach facilitates rapid sense making, pattern recognition, and anomaly detection. Country-based profiles summarize vast amounts of data into readily understandable graphical representations of each country’s primary strategic imports and exports and their main trading partners. Commodity-based profiles provide a global view of major exporters, importers, and trade patterns for each commodity. These insights can inform the risk management approach employed by customs to allocate enforcement and compliance efforts.

customs controls↗

Enabling BWR fuel rod analysis in the BISON fuel performance code

Nuclear fuel vendors around the world are pursuing approaches to sustain the existing nuclear reactor fleet consisting primarily of pressurized-water reactors (PWRs) and boiling-water reactors (BWRs). To support the industry's efforts, advanced modeling and simulation tools need to be capable of analyzing both legacy reactor concepts. BWR fuel rods are significantly different than those used in PWRs, which can affect fuel performance analysis. BWR fuel rods include: (1) an extensive use of gadolinia dopant as a burnable absorber, (2) an axial variation in fuel enrichment and gadolinia content, (3) the inclusion of a liner on the inner cladding surface to mitigate the impact of pellet-clad mechanical interaction (which impacts hydrogen and hydride distribution), (4) a lower initial fill gas pressure, (5) bottom-entry control rods, and (6) a lower coolant pressure that results in the two-phase flow boiling phenomenon. Although the primary focus of BISON has been in the area of PWR and advanced reactor fuel analyses, this paper presents the developments in BISON to support BWR fuel performance analysis. An overview of the models that account for the effects of gadolinia is highlighted. Internal mesh generation capabilities to include a liner is presented. Normal operating and transient (reactivity insertion accident) demonstration problems are presented to illustrate the impact of gadolinia, the hydrogen and hydride evolution due to the presence of the liner, and BISON's ability to simulate axially varying enrichments and dopant concentration. Bottom-entry control effects are captured by the axial power peaking factors present in the demonstration cases. Comparisons to integral experiments from the Halden IFA-681 experiments are discussed as initial validation. Reasonable comparisons are obtained for fuel centerline temperature and rod internal pressure as a function of time. In conclusion, simulations of additional experiments containing Gd 2 O 3 -bearing fuel are necessary to completely validate the code for BWR applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

IRIS Reimagined: Advancements in Intelligent Runtime System for Task-Based Programming

Task-based programming models are gaining traction in scientific computing. IRIS is a portable runtime system that exploits multiple heterogeneous programming systems and can discover available resources and manage multiple diverse programming systems (e.g., CUDA, Hexagon, HIP, Level Zero, OpenCL, and OpenMP) simultaneously. It accounts for the constraints of task dependencies and provides customizable scheduling policies to map those tasks to heterogeneous devices. In this paper, we present new capabilities added to IRIS to improve its portability for heterogeneous programming, build-friendliness, and performance efficiency. The new additions include vendor-specific kernel support, a runtime system with a foreign function interface to eliminate writing wrapper or boilerplate code for heterogeneous kernels, an easy-to-use and configurable CMake-based build environment, automatic and efficient data transfers and orchestration, and the Hunter and DAGGER toolchains to evaluate IRIS’s task scheduling algorithms.

Miniskar, Narasinga Rao↗

The impact of incorporating shell-corrections to energy loss in silicon

Modern silicon tracking detectors based on hybrid or fully integrated CMOS technology are continuing to push to thinner sensors. The ionization energy loss fluctuations in very thin silicon sensors significantly deviates from the Landau distribution. Therefore, we have developed a charge deposition setup that implements the Bichsel straggling function, which accounts for shell-effects. This enhanced simulation is important for comparing with testbeam or collision data with thin sensors as demonstrated by reproducing more realistically the degraded position resolution compared with naïve ionization models based on simple Landau-like fluctuation. Our implementation of the Bichsel model agrees well with the multipurpose photo absorption ionization (PAI) model in Geant4 and is significantly faster. The code is made publicly available as part of the Allpix software package in order to facilitate predictions for new detector designs and comparisons with testbeam data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Implementation of the D1S Methodology for Shutdown Dose Rate Calculations in the OpenMC Monte Carlo Particle Transport Code

We present an implementation of the direct one-step (D1S) methodology for shutdown dose rate (SDR) calculations in the OpenMC Monte Carlo particle transport code. In addition to being the first fully open-source D1S implementation, it is also the first to require no ad hoc source code or nuclear data library modifications. The code can seamlessly switch between production of prompt and decay photons based on a user input parameter, and the decay data needed for decay photon generation are made available through a depletion chain file, which is already used for OpenMC’s built-in depletion/activation solver. A set of Python functions significantly eases the burden of computing and applying time correction factors needed to properly account for the time dependence of radionuclide activity. To assess the accuracy of the D1S implementation, SDR calculations have been carried out for three problems: a prism of iron irradiated by 14-MeV neutrons, the ITER port plug computational benchmark, and the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For each of these problems, comparisons were made to calculations using the rigorous two-step (R2S) method. The results on the iron prism problem illustrate how the D1S method achieves superior spatial resolution compared to the R2S method without the need for spatial discretization of the activation regions. The D1S and R2S results for the ITER port plug benchmark agree well with previously reported results in the literature. While the D1S results are 10% to 15% lower than the R2S results, this may be due to stochastic uncertainty and/or spatial discretization in the R2S calculations. On the FNG dose rate benchmark problem, the D1S method produces dose rate estimates that are within 4% of the dose rates predicted using a cell-based R2S workflow. The D1S estimates of the SDR are also in reasonable agreement with the experimental measurements and show the same basic trends that have been observed in previous works. A qualitative analysis of the execution time and uncertainty for the R2S and D1S workflows suggests that the D1S method would attain a higher figure of merit.

D1S method↗

Map Applications to Target Exascale Architecture with Machine-Specific Performance Analysis, Including Challenges and Projections

This Exascale Computing Project (ECP) milestone report summarizes the status of all 30 ECP Applications Development (AD) subprojects at the end of FY20. In October and November of 2020, a comprehensive assessment of AD projects was conducted by the ECP leadership. Reviews occurred virtually between October 27, 2020 and November 12, 2020. The review committee—consisting of the AD lead, deputy, and L3—was tasked with evaluating each subproject’s progress in porting their codes to early exascale architectures considered precursors to the planned exascale machines. This includes characterizing which modules have been ported to multi-accelerator nodes, initial performance analyses, the status of software integration, and a current vision of successes, obstacles, and next steps. As such, this report contains not only an accurate snapshot of each subproject’s current status but also represents an unprecedentedly broad account of experiences in porting large scientific applications to next-generation high-performance computing architectures.

97 MATHEMATICS AND COMPUTING↗

Enhance BISON Metal Fuel Transient Models and Enable Assessment Based on Out-of-Pile Transient Experiments

In this report, the activities and achievements made by Argonne National Laboratory for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) BISON code metallic fuel validation and verification (V&V) project in FY2023, with a focus on the capability development and assessment case development for the out-of-pile transient experiments, are summarized. The liquid cladding penetration model was improved based on a comprehensive set of fuel behavior test apparatus (FBTA) test results. The new model was implemented into BISON along with a simple model to account for fuel melting effects. The metallic fuel meshing capability was also upgraded to allow local refinement using unstructured triangular mesh. Based on the previous success in developing the BISON-FIPD integration powered BISON assessment case, this database enabled BISON evaluation and V&V capability was expanded to the out-of-pile transient experiments by introducing BISON-OPTD integration. A series of initial out-of-pile transient assessment cases were implemented into BISON to provide powerful platforms to evaluate current BISON’s capabilities in simulating metallic fuel performance during out-of-pile transients. The established capabilities will also benefit the calibration and V&V of future advanced mechanistic models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Transonic flow visualization using holographic interferometry

An account is made of some of the applications of holographic interferometry to the visualization of transonic flows. In the case of the compressor shock visualization, the method is used regularly and has moved from being a research department invention to a design test tool. With the implementation of automatic processing and simple digitization systems, holographic vibrational analysis has also moved into routine nondestructive testing. The code verification interferograms were instructive, but the main turbomachinery interest is now in 3 dimensional flows. A major data interpretation effort will be required to compute tomographically the 3 dimensional flow around the leading or the trailing edges of a rotating blade row. The bolt on approach shows the potential application to current unsteady flows of interest. In particular that of the rotor passing and vortex interaction effects is experienced by the new generation of unducted fans. The turbocharger tests presents a new area for the application of holography.

Bryanston-Cross, Peter J.↗

Computations of spray, fuel-air mixing, and combustion in a lean-premixed-prevaporized combustor

A code was developed for computing the multidimensional flow, spray, combustion, and pollutant formation inside gas turbine combustors. The code developed is based on a Lagrangian-Eulerian formulation and utilizes an implicit finite-volume method. The focus of this paper is on the spray part of the code (both formulation and algorithm), and a number of issues related to the computation of sprays and fuel-air mixing in a lean-premixed-prevaporized combustor. The issues addressed include: (1) how grid spacings affect the diffusion of evaporated fuel, and (2) how spurious modes can arise through modelling of the spray in the Lagrangian computations. An upwind interpolation scheme is proposed to account for some effects of grid spacing on the artificial diffusion of the evaporated fuel. Also, some guidelines are presented to minimize errors associated with the spurious modes.

Dasgupta, A.↗

The Acoustic Influence of Cell Depth on the Rotordynamic Characteristics of Smooth-Rotor/Honeycomb-Stator Annular Gas Seals

A two-control volume is employed for honeycomb-stator/smooth-rotor seals, with a conventional control-volume used for the through flow and a 'capacitance accumulator' model for the honeycomb cells. The control volume for the honeycomb cells is shown to cause a dramatic reduction in the effective acoustic velocity of the main flow, dropping the lowest acoustic frequency into the frequency range of interest for rotordynamics. In these circumstances, the impedance functions for the seals can not be modeled with conventional (frequency-independent) stiffness, damping, and mass coefficients. More general transfer functions are required to account for the reaction forces, and calculated here as a lead-lag term for the direct force function and a lag term for the cross-coupled function. These first order functions are simple compared to transfer functions for magnetic bearings or foundations, For synchronous response to imbalance, they can be approximated by running-speed-dependent stiffness and damping coefficients in conventional rotordynamic codes. Correct predictions for stability and transient response will require more general algorithms, pressumably using a state-space format.

Childs, Dara W.↗

Studying the MLT by a Combined Analysis of SABER/TIMED and Lidar Measurements

The SABER instrument on board the TIMED Satellite is a limb scanning infrared radiometer designed to measure temperature and minor constituent vertical profiles and energetics parameters in the mesosphere and lower thermosphere (MLT). The measurements have been performed continuously since January 25, 2002 to provide excellent coverage for both hemispheres. The Leibniz-Institute of Atmospheric Physics (LAP) at Kuehlungsborn, Germany (54N, 12E) operates two lidar instruments, using three different temperature measurement methods, optimized for three altitude ranges. The total altitude range of the lidar installation lies from 1 to 105 km. Another instrument used for intercomparison is the ALOMAR RMR lidar, located at Andoya, Norway (69N, 16E). We have searched the SABER and lidar datasets for coincidental common volume measurements within plus or minus 1 degree in latitude, plus or minus 2 degrees in longitude and approx. 1 hour in time for the sake of (a) comparison of measured temperatures; (b) validation of the models used in SABER data analysis; and (c) extracting new information about MLT parameters. In this work we applied the non-LTE ALI-ARMS code designed to calculate the nonequilibrium radiance in different viewing geometries to the analysis of measurements which satisfied these search criteria. The results of this analysis (a) support the application of higher value of CO2-O quenching rate (6e-12 cubic centimeters per second) by the non-LTE temperature retrievals from the SABER 15 micrometer limb radiance data, and (b) demonstrate the importance of accounting for the vibrational-vibrational energy exchange among the CO2 isotopes for accurate temperature retrievals. Using temperature profiles obtained in lidar measurements as inputs for the retrieval algorithm we also retrieved the nighttime CO2 densities from the SABER 15 micrometer limb radiances and compared them with the model and climatology CO2 data used in the SABER nighttime temperature retrievals.

Feofilov, A. G.↗

3D Visualization of Monte-Carlo Simulation's of HZE Track Structure and Initial Chemical Species

Heavy ions biophysics is important for space radiation risk assessment [1] and hadron-therapy [2]. The characteristic of heavy ions tracks include a very high energy deposition region close to the track (<20 nm) denoted as the track core, and an outer penumbra region consisting of individual secondary electrons (6-rays). A still open question is the radiobiological effects of 6- rays relative to the track core. Of importance is the induction of double-strand breaks (DSB) [3] and oxidative damage to the biomolecules and the tissue matrix, considered the most important lesions for acute and long term effects of radiation. In this work, we have simulated a 56Fe26+ ion track of 1 GeV/amu with our Monte-Carlo code RITRACKS [4]. The simulation results have been used to calculate the energy depiction and initial chemical species in a "voxelized" space, which is then visualized in 3D. Several voxels with dose >1000 Gy are found in the penumbra, some located ~0.1 mm from the track core. In computational models, the DSB induction probability is calculated with radial dose [6], which may not take into account the higher RBE of electron track ends for DSB induction. Therefore, these simulations should help improve models of DSB induction and our understanding of heavy ions biophysics.

Plante, Ianik↗