Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coded computing”

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,297 records · Page 72

A general and computationally fast formulation for radiative transfer with scattering

A general formulation of monocromatic radiative transfer with scattering has been developed for plane-parallel geometry. The inhomogeneous and nonisothermal medium absorbs, emits, and anisotropically scatters radiation. Surfaces can emit and scatter radiation in any specified manner. The solution procedure uses the fact that phase incoherent scattering is linear in radiative sources. Certain basic scattering functions are then defined and calculated by an adding computer code using matrix algebra. These scattering functions are weighted by the temperature field and summed (superimposed) to obtain the solution for any specific problem. Numerical results for exiting intensities and one-sided heat fluxes from general media bound by one arbitrary surface are presented. These parametric studies demonstrate the effects of scattering particles and surfaces on radiative transfer from inhomogeneous and nonisothermal media. Application of the formulation to radiative equilibrium is also discussed. The conclusion is that all problems in plane-parallel radiative transfer with scattering can be solved by a common and computationally fast algorithm based on this formulation.

Cogley, A. C.↗

Numerical predictions of residence times behind a rearward facing step with transverse injection

A numerical investigation of the two-dimensional recirculation zone formed between a rear-ward facing step and a transverse hydrogen jet in a supersonic flow is presented. Calculations were made using an elliptic flow computer code with a time-relaxation finite difference integration algorithm. Cases were examined for jet static pressures between 2 and 5 times mainstream with jets located 3 to 7 step heights downstream of the step. Nominal mainstream flow conditions were Mach 2.4 air at 1 atm pressure and 1000K (1800R) static temperature. Results in the form of mean properties in the recirculation zone and residence time are presented. Calculations of the decay of hydrogen mass from the steady state are presented for selected cases. The combustion potential for each case is assessed by comparison of computed conditions in the recirculation zone with the results from a well-stirred reactor theory.

Rogers, R. C.↗

The Scientific Impact of the Exascale Computing Project

The recent arrival of the Frontier Supercomputer at Oak Ridge National Laboratory officially marked the dawn of the exascale computing era. Its successful deployment coincided with the culmination of the U.S. Department of Energy Exascale Computing Project (ECP), an ambitious, complex, and risky research and development effort that integrated contributions from a broad and diverse subset of the high-performance computing community. The success of ECP will ultimately be judged by the scientific and engineering advances that it enabled. In conclusion, this Special Issue is focused on showcasing early successes in the use of exascale resources to enable breakthroughs in key areas of science in engineering.

97 MATHEMATICS AND COMPUTING↗

Effectiveness evaluation of STOL transport operations

A short-takeoff and landing (STOL) systems simulation model has been developed and implemented in a computer code (known as STOL OPS) which permits evaluation of the operation of a STOL aircraft and its avionics in a commercial airline operating environment. STOL OPS concentrated on the avionics functions of navigation, guidance, control, communication, hazard aviodance, and systems management. External world factors influencing the operation of the STOL aircraft include each airport and its geometry, air traffic at each airport, air traffic control equipment and procedures, weather (including winds and visibility), and the flight path between each airport served by the route. The development of the STOL OPS program provides NASA a set of computer programs which can be used for detailed analysis of a STOL aircraft and its avionics and permit establishment of system requirements as a function of airline mission performance goals.

Hitt, E. F.↗

Three-dimensional computer model for the atmospheric general circulation experiment

An efficient, flexible, three-dimensional, hydrodynamic, computer code has been developed for a spherical cap geometry. The code will be used to simulate NASA's Atmospheric General Circulation Experiment (AGCE). The AGCE is a spherical, baroclinic experiment which will model the large-scale dynamics of our atmosphere; it has been proposed to NASA for future Spacelab flights. In the AGCE a radial dielectric body force will simulate gravity, with hot fluid tending to move outwards. In order that this force be dominant, the AGCE must be operated in a low gravity environment such as Spacelab. The full potential of the AGCE will only be realized by working in conjunction with an accurate computer model. Proposed experimental parameter settings will be checked first using model runs. Then actual experimental results will be compared with the model predictions. This interaction between experiment and theory will be very valuable in determining the nature of the AGCE flows and hence their relationship to analytical theories and actual atmospheric dynamics.

Roberts, G. O.↗

An incompressible Navier-Stokes flow solver in three-dimensional curvilinear coordinate systems using primitive variables

An implicit, finite-difference computer code has been developed to solve the incompressible Navier-Stokes equations in a three-dimensional, curvilinear coordinate system. The pressure-field solution is based on the pseudo compressibility approach in which the time derivative pressure term is introduced into the mass conservation equation to form a set of hyperbolic equations. The solution procedure employs an implicit, approximate factorization scheme. The Reynolds stresses, that are uncoupled from the implicit scheme, are lagged by one time-step to facilitate implementing various levels of the turbulence model. Test problems for external and internal flows are computed, and the results are compared with existing experimental data. The application of this technique for general three-dimensional problems is then demonstrated.

Kwak, D.↗

Numerical predictions of residence times behind a rearward facing step with transverse injection

A numerical investigation of the two-dimensional recirculation zone formed between a rearward facing step and a transverse hydrogen jet in a supersonic flow is presented. Calculations were made using an elliptic flow computer code with a time-relaxation finite difference integration algorithm. Cases were examined for jet static pressures between 2 and 5 times mainstream with jets located 3 to 7 step heights downstream of the step. Nominal mainstream flow conditions were Mach 2.4 air at 1 atm pressure and 1000 K (1800 R) static temperature. Results in the form of mean properties in the recirculation zone and residence time are presented. Calculations of the decay of hydrogen mass from the steady state are presented for selected cases. The combustion potential for each case is assessed by comparison of computed conditions in the recirculation zone with the results from a well-stirred reactor theory.

Rogers, R. C.↗

Development of unsteady aerodynamic analyses for turbomachinery aeroelastic and aeroacoustic applications

Theoretical analyses and computer codes are being developed for predicting compressible unsteady inviscid and viscous flows through blade rows. Such analyses are needed to determine the impact of unsteady flow phenomena on the structural durability and noise generation characteristics of turbomachinery blading. Emphasis is being placed on developing analyses based on asymptotic representations of unsteady flow phenomena. Thus, flow driven by small-amplitude unsteady excitations in which viscous effects are concentrated in thin layers are being considered. The resulting analyses should apply in many practical situations, lead to a better understanding of the relevent physics, and they will be efficient computationally, and therefore, appropriate for aeroelastic and aeroacoustic design applications. Under the present phase (Task 3), the effort was focused on providing inviscid and viscid prediction capabilities for subsonic unsteady cascade flows.

Verdon, Joseph M.↗

Assessment of Critical Experiment Benchmark Applicability to a Large-Capacity HALEU Transportation Package Concept

This work presents an assessment of the applicability of existing benchmark critical experiments to the criticality safety code validation for a large-capacity high-assay low-enriched uranium (HALEU) transportation package concept. Numerous next-generation nuclear reactor designs require HALEU fuel, which is characterized by an enrichment between 5 and 20 wt% 235 U. The U.S. Department of Energy (DOE) has proposed to recover and downblend highly enriched uranium from DOE-owned used nuclear fuel to accelerate the demonstration of commercially viable microreactor technologies. One element of the infrastructure needed to demonstrate HALEU-fueled reactors is the ability to safely transport enriched product to be used for fuel fabrication. There is uncertainty as to whether existing critical benchmark experiment data are sufficient to support criticality safety code validation for HALEU transportation applications. The anticipated chemical form of the HALEU in the proposed transportation concept is UO 2 with 20 wt% 235 U/U. The concept uses a combination of an existing transportation packaging design and a novel basket design, including borated aluminum flux traps. The basket provides space for 18 reusable, stainless steel canisters that contain the HALEU. In 10 CFR 71, normal conditions of transport (NCTs) and hypothetical accident conditions (HACs) are defined for fissile material transportation packages. NCT and HAC KENO-VI models of the transportation package were developed using the Standardized Computer Analyses for Licensing Evaluation (SCALE) 6.2.3 computer code package, and optimum moderation conditions were determined using the SCALE SAMPLER sequence. The SCALE Tools for Sensitivity and Uncertainty Analysis Methodology Implementation (TSUNAMI) sequences were then used to compare the neutronic characteristics of 1584 International Criticality Safety Benchmark Evaluation Project benchmark critical experiments with the NCT and HAC HALEU transportation models. The TSUNAMI integral correlation coefficient c k was the criterion used to rank neutronic similarity. Thirty-four experiments were identified as similar (c k ≥ 0.9) to the NCT model, and 55 experiments were identified as similar to the HAC model. Hundreds of experiments were also identified as at least marginally similar (c k ≥ 0.8) to both models. The results indicate that additional critical experiments are unlikely to be needed to support HALEU transportation criticality safety analyses for package concepts similar to the concept package analyzed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computer programs for smoothing and scaling airfoil coordinates

Detailed descriptions are given of the theoretical methods and associated computer codes of a program to smooth and a program to scale arbitrary airfoil coordinates. The smoothing program utilizes both least-squares polynomial and least-squares cubic spline techniques to smooth interatively the second derivatives of the y-axis airfoil coordinates with respect to a transformed x-axis system which unwraps the airfoil and stretches the nose and trailing-edge regions. The corresponding smooth airfoil coordinates are then determined by solving a tridiagonal matrix of simultaneous cubic-spline equations relating the y-axis coordinates and their corresponding second derivatives. A technique for computing the camber and thickness distribution of the smoothed airfoil is also discussed. The scaling program can then be used to scale the thickness distribution generated by the smoothing program to a specific maximum thickness which is then combined with the camber distribution to obtain the final scaled airfoil contour. Computer listings of the smoothing and scaling programs are included.

Morgan, H. L., Jr.↗

Finite-volume Euler and Navier-Stokes solvers for three-dimensional and conical vortex flows over delta wings

A unified central-difference finite-volume Euler and Navier-Stokes solver with four-stage Runge-Kutta time stepping is presented. The computer code developed for this purpose is capable of solving the standard set and nonstandard sets (zero-total-pressure loss) of Euler equations and the thin-layer and full Navier-Stokes equations. Applications are presented for conical supersonic flows with weak shocks using the standard and nonstandard sets of Euler equations, and the thin-layer and full Navier-Stokes equations for sharp and round-edged delta wings. Applications are also presented for three-dimensional transonic and subsonic flows using the standard set of Euler equations for sharp-edged delta wings. The computational results of the different sets of equations are compared with each other and with the experimental results and conclusions on the validity of these sets to these applications, are presented.

Kandil, Osama A.↗

Overview of the Tolerance Limit Calculations with Application to TSURFER

To establish confidence in the results of computerized physics models, a key regulatory requirement is to develop a scientifically defendable process. The methods employed for confidence, characterization, and consolidation, or C3, are statistically involved and are often accessible only to avid statisticians. This manuscript serves as a pedagogical presentation of the C3 process to all stakeholders—including researchers, industrial practitioners, and regulators—to impart an intuitive understanding of the key concepts and mathematical methods entailed by C3. The primary focus is on calculation of tolerance limits, which is the overall goal of the C3 process. Tolerance limits encode the confidence in the calculation results as communicated to the regulator. Understanding the C3 process is especially critical today, as the nuclear industry is considering more innovative ways to assess new technologies, including new reactor and fuel concepts, via an integrated approach that optimally combines modeling and simulation and minimal targeted validation experiments. This manuscript employs intuitive, analytical, numerical, and visual representations to explain how tolerance limits may be calculated for a wide range of configurations, and it also describes how their values may be interpreted. Various verification tests have been developed to test the calculated tolerance limits and to help delineate their values. The manuscript demonstrates the calculation of tolerance limits for TSURFER, a computer code developed by the Oak Ridge National Laboratory for criticality safety applications. The goal is to evaluate the tolerance limit for TSURFER-determined criticality biases to support the determination of upper, subcritical limits for regulatory purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Unsteady Aerodynamic Models for Turbomachinery Aeroelastic and Aeroacoustic Applications

Theoretical analyses and computer codes are being developed for predicting compressible unsteady inviscid and viscous flows through blade rows of axial-flow turbomachines. Such analyses are needed to determine the impact of unsteady flow phenomena on the structural durability and noise generation characteristics of the blading. The emphasis has been placed on developing analyses based on asymptotic representations of unsteady flow phenomena. Thus, high Reynolds number flows driven by small amplitude unsteady excitations have been considered. The resulting analyses should apply in many practical situations and lead to a better understanding of the relevant flow physics. In addition, they will be efficient computationally, and therefore, appropriate for use in aeroelastic and aeroacoustic design studies. Under the present effort, inviscid interaction and linearized inviscid unsteady flow models have been formulated, and inviscid and viscid prediction capabilities for subsonic steady and unsteady cascade flows have been developed. In this report, we describe the linearized inviscid unsteady analysis, LINFLO, the steady inviscid/viscid interaction analysis, SFLOW-IVI, and the unsteady viscous layer analysis, UNSVIS. These analyses are demonstrated via application to unsteady flows through compressor and turbine cascades that are excited by prescribed vortical and acoustic excitations and by prescribed blade vibrations. Recommendations are also given for the future research needed for extending and improving the foregoing asymptotic analyses, and to meet the goal of providing efficient inviscid/viscid interaction capabilities for subsonic and transonic unsteady cascade flows.

Verdon, Joseph M.↗

Integrated control and health management. Orbit transfer rocket engine technology program

To insure controllability of the baseline design for a 7500 pound thrust, 10:1 throttleable, dual expanded cycle, Hydrogen-Oxygen, orbit transfer rocket engine, an Integrated Controls and Health Monitoring concept was developed. This included: (1) Dynamic engine simulations using a TUTSIM derived computer code; (2) analysis of various control methods; (3) Failure Modes Analysis to identify critical sensors; (4) Survey of applicable sensors technology; and, (5) Study of Health Monitoring philosophies. The engine design was found to be controllable over the full throttling range by using 13 valves, including an oxygen turbine bypass valve to control mixture ratio, and a hydrogen turbine bypass valve, used in conjunction with the oxygen bypass to control thrust. Classic feedback control methods are proposed along with specific requirements for valves, sensors, and the controller. Expanding on the control system, a Health Monitoring system is proposed including suggested computing methods and the following recommended sensors: (1) Fiber optic and silicon bearing deflectometers; (2) Capacitive shaft displacement sensors; and (3) Hot spot thermocouple arrays. Further work is needed to refine and verify the dynamic simulations and control algorithms, to advance sensor capabilities, and to develop the Health Monitoring computational methods.

Holzmann, Wilfried A.↗

Effects of airfoil shape, thickness, camber, and angle of attack on calculated transonic unsteady airloads

The effects of airfoil shape, thickness, camber, and mean angle of attack on transonic unsteady airloads were investigated as calculated by the transonic small-disturbance computer code XTRAN2L. Shape effects were investigated by examining the pressure distributions, shock locations, and unsteady airloads for three 10 percent thick airfoils. NACA 0010, NACA 64A010, and parabolic arc. Thickness effects were determined by studying a single airfoil shape with three different thicknesses: NACA 0008, NACA 0010, and NACA 0012. Angle of attack and camber effects were studied by including mean angle of attack or by adding a simple parabolic camber distribution to the originally symmetric airfoils. Comparisons of unsteady airloads for different airfoil configurations show similar results caused by variations in airfoil shape, thickness, camber, or mean angle of attack. Computer costs can be reduced by limiting the number of transonic unsteady aerodynamic calculations for small changes in airfoil geometry or angle of attack.

Batina, J. T.↗

Closed Loop Geothermal Working Group: GeoCLUSTER App, Subsurface Simulation Results, and Publications

To better understand the heat production, electricity generation performance, and economic viability of closed-loop geothermal systems in hot-dry rock, the Closed-Loop Geothermal Working Group -- a consortium of several national labs and academic institutions has tabulated time-dependent numerical solutions and levelized cost results of two popular closed-loop heat exchanger designs (u-tube and co-axial). The heat exchanger designs were evaluated for two working fluids (water and supercritical CO2) while varying seven continuous independent parameters of interest (mass flow rate, vertical depth, horizontal extent, borehole diameter, formation gradient, formation conductivity, and injection temperature). The corresponding numerical solutions (approximately 1.2 million per heat exchanger design) are stored as multi-dimensional HDF5 datasets and can be queried at off-grid points using multi-dimensional linear interpolation. A Python script was developed to query this database and estimate time-dependent electricity generation using an organic Rankine cycle (for water) or direct turbine expansion cycle (for CO2) and perform a cost assessment. This document aims to give an overview of the HDF5 database file and highlights how to read, visualize, and query quantities of interest (e.g., levelized cost of electricity, levelized cost of heat) using the accompanying Python scripts. Details regarding the capital, operation, and maintenance and levelized cost calculation using the techno-economic analysis script are provided. This data submission will contain results from the Closed Loop Geothermal Working Group study that are within the public domain, including publications, simulation results, databases, and computer codes. GeoCLUSTER is a Python-based web application created using Dash, an open-source framework built on top of Flask that streamlines the building of data dashboards. GeoCLUSTER provides users with a collection of interactive methods for streamlining the exploration and visualization of an HDF5 dataset. The GeoCluster app and database are contained in the compressed file geocluster_vx.zip, where the "x" refers to the version number. For example, geocluster_v1.zip is Version 1 of the app. This zip file also contains installation instructions. **To use the GeoCLUSTER app in the cloud, click the link to "GeoCLUSTER on AWS" in the Resources section below. To use the GeoCLUSTER app locally, download the geocluster_vx.zip to your computer and uncompress this file. When uncompressed this file comprises two directories and the geocluster_installation.pdf file. The geo-data app contains the HDF5 database in condensed format, and the GeoCLUSTER directory contains the GeoCLUSTER app in the subdirectory dash_app, as app.py. The geocluster_installation.pdf file provides instructions on installing Python, the needed Python modules, and then executing the app.

15 GEOTHERMAL ENERGY↗

Verification of the PERSENT Software

Ongoing commercial design activities require a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system and substantial work has been done to verify its accuracy on several identified commercial needs. This manuscript details the verification work done on PERSENT which relies upon the DIF3D code for its forward and adjoint flux solution. Previous work identified the PERSENT features required to be verified to support commercial design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying PERSENT’s ability to correctly calculate four key quantities: perturbation worth distributions, kinetics parameters, sensitivity coefficients, and cross section uncertainty quantification. This manuscript provides the verification tasks and their results with respect to these quantities needed for commercial design activities. For the perturbation worth distributions, hand calculations are deployed to verify the PERSENT calculated results. Similarly, hand calculation of the PERSENT computed kinetics parameters is also used to verify the PERSENT results. In both of these, the input to PERSENT is manipulated to ensure the hand calculation exactly matches the equations PERSENT is calculating. The sensitivity coefficients involve calculating the derivatives of a parameter (such as reactivity worth), with respect to the cross section data. Direct finite difference calculations with DIF3D are used to verify the PERSENT calculated results. For the uncertainty quantification, manufactured input to PERSENT is used to allow an exact hand calculation to reproduce the PERSENT calculated results. The work detailed in this report verified that significant issues were identified for earlier versions of PERSENT for sensitivity coefficients which were corrected in this work and thus version 12.1.0 of PERSENT must be used to reproduce all of the verified work in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

OMV/VTE variable thrust engine analysis

The objective of the present work is to develop a predictive CFD based analytical tool for the variable thrust engine (VTE) in the orbital maneuvering vehicle (OMV). This objective is being accomplished within the framework of the Los Alamos KIVA computer code for chemically reactive flows with sprays. For the OMV application, the main structure of KIVA is to be retained while reformulating many of the phenomenological submodels, enhancing some of the numerics, and adding more features. The analytical model consists of the general conservation equations for two phase reactive flows and of submodels for turbulence, chemical reactions, and bipropellant sprays. Tailoring this model to the OMV engine brings about the added complexities of combustion and flow processes that occur in a liquid hypergolic propellant rocket chamber. This report exposes the foundation upon which the analytical tool is being constructed and developed. Results from a cursory computational exercise involving the simulation of the flow and combustion processes in a hypothetical N2H4/N2O4 rocket engine thrust chamber is presented and discussed.

Larosillere, Louis↗