High-Performance Computational Electromagnetic Modeling Using Low-Cost Parallel Computers
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Predictive computational models have been developed to support the technical basis for extended dry storage of aluminum-clad spent nuclear fuel (ASNF) in helium-backfilled cannisters. To date, these models have been optimized on a variety of irradiation experiments designed to elucidate the radiation-induced formation of molecular hydrogen gas, a radiolysis product that is potentially problematic for the safe storage of ASNF. However, the yield of nitric acid (HNO3) has not been determined, despite conservative predictions of its formation (300–4000 ppm) in 1% residual air environments irradiated in contact with ASNF. HNO3, another problematic radiolysis product, can lead to enhanced corrosion and potentially compromise storage canister integrity. Thus, to support the validation of predictive computer models, we report the measurement and quantification of HNO3 from the gamma irradiation of aluminum alloy coupons in humid air.
Abstract Motivated by a computer model calibration problem from the oil and gas industry, involving the design of a honeycomb seal, we develop a new Bayesian methodology to cope with limitations in the canonical apparatus stemming from several factors. We propose a new strategy of on‐site design and surrogate modeling for a computer simulator acting on a high‐dimensional input space that, although relatively speedy, is prone to numerical instabilities, missing data, and nonstationary dynamics. Our aim is to strike a balance between data‐faithful modeling and computational tractability in a calibration framework—tailoring the computer model to a limited field experiment. Situating our on‐site surrogates within the canonical calibration apparatus requires updates to that framework. We describe a novel yet intuitive Bayesian setup that carefully decomposes otherwise prohibitively large matrices by exploiting the sparse blockwise structure. Empirical illustrations demonstrate that this approach performs well on toy data and our motivating honeycomb example.
Ionospheric computer model, describing electron and ion densities in middle latitudes as function of latitude, longitude, altitude, season, time and solar activity
The computer model for erythropoietic control was adapted to the mouse system by altering system parameters originally given for the human to those which more realistically represent the mouse. Parameter values were obtained from a variety of literature sources. Using the mouse model, the mouse was studied as a potential experimental model for spaceflight. Simulation studies of dehydration and hypoxia were performed. A comparison of system parameters for the mouse and human models is presented. Aside from the obvious differences expected in fluid volumes, blood flows and metabolic rates, larger differences were observed in the following: erythrocyte life span, erythropoietin half-life, and normal arterial pO2.
A computer model was developed to assess the noise impact of an airport on the community which it serves. Assessments are made using the Fractional Impact Method by which a single number describes the community aircraft noise environment in terms of exposed population and multiple event noise level. The model is comprised of three elements: a conventional noise footprint model, a site specific population distribution model, and a dose response transfer function. The footprint model provides the noise distribution for a given aircraft operating scenario. This is combined with the site specific population distribution obtained from a national census data base to yield the number of residents exposed to a given level of noise. The dose response relationship relates noise exposure levels to the percentage of individuals highly annoyed by those levels.
Numerical computations of stationary solar coronal loop atmospheres are used to extend earlier analytical work. Two classes of loops are examined, namely symmetric loops with a temperature maximum at the top but now having a length greater than the pressure scale height and loops which have a local temperature minimum at the top. For the first class, new scaling laws are found which relate the base pressure and loop length to the base heating, the heating deposition scale height, and the pressure scale height. It is found that loops for which the length is greater than about two to three times the pressure scale height do not have stable solutions unless they have a temperature minimum at the top. Computed models with a temperature inversion at the top are permitted in a wider range of heating deposition scale height values than are loops with a temperature maximum at the top. These results are discussed in relation to observations showing a dependence of prominence formation and stability on the state of evolution of magnetic structures, and a general scenario is suggested for the understanding of loop evolution from emergence in active regions through the large-scale structure phase to opening in coronal holes.
A computer model is described which uses mixed-integer linear programming to provide optimal DSN spacecraft schedules given a mission set and specified scheduling requirements. A solution technique is proposed which uses Bender's Method and a heuristic starting algorithm.
A computer model has been constructed to simulate the compliance and load sharing in a spur gear mesh. The model adds the effect of rim deflections to previously developed state-of-the-art gear tooth deflection models. The effects of deflections on mesh compliance and load sharing are examined. The model can treat gear meshes composed to two external gears or an external gear driving an internal gear. The model includes deflection contributions from the bending and shear in the teeth, the Hertzian contact deformations, and primary and secondary rotations of the gear rims. The model shows that rimmed gears increase mesh compliance and, in some cases, improve load sharing.
A computer model has been constructed to simulate the compliance and load sharing in a spur gear mesh. The model adds the effect of rim deflections to previously developed state-of-the-art gear tooth deflection models. The effects of deflections on mesh compliance and load sharing are examined. The model can treat gear meshes composed of two external gears or an external gear driving an internal gear. The model includes deflection contributions from the bending and shear in the teeth, the Hertzian contact deformations, and primary and secondary rotations of the gear rims. The model shows that rimmed gears increase mesh compliance and, in some cases, improve load sharing.
Computing cosmic ray cutoff rigidity for several geomagnetic field models
Computer program for calculating of power-density spectrum (PDS) from data base generated by Advanced Continuous Simulation Language (ACSL) uses algorithm that employs fast Fourier transform (FFT) to calculate PDS of variable. Accomplished by first estimating autocovariance function of variable and then taking FFT of smoothed autocovariance function to obtain PDS. Fast-Fourier-transform technique conserves computer resources.
IBR-2M is a fast research reactor located in Dubna (Russia), which operates in supercritical condition for ~800 μs every 200 ms. Two reflector parts are rotating in opposite directions generating 1.8 GWth peak power when the two parts align with the fuel zone. The reactor core uses plutonium fuel and sodium coolant. The control rods are located in stationary reflector zones. This report describes the computational modeling of the IBR-2M reactor. Two independent computational models were generated using MCNP and SERPENT computer programs. The MCNP model uses combinatorial geometry based on macro bodies, e.g. boxes, cylinders, hexagons, spheres, cones, etc., and surfaces. The SERPENT model is based on the STereo Lithographic (STL) geometry representation that is used by 3D printers. The STL geometry was constructed using the CUBIT computer program. The water cooling system of the stationary and the rotating reflectors has a very complicated geometry, involving serpentines and irregular surfaces profiles. The materials in these areas have been homogenized using the GIMP image editor according to the material distributions obtained from the colors of the design drawings. The MCNP and the SERPENT inputs consist of ~1500 and ~1750 lines, respectively. The calculated effective multiplication factors from these two Monte Carlo computer programs with different models agree within 50 pcm.
We present a computational model of a closed-loop, pilot-aircraft-visual scene-taxiway system created to shed light on possible sources of taxi error. Creating the cognitive aspects of the model using ACT-R required us to conduct studies with subject matter experts to identify experiential adaptations pilots bring to taxiing. Five decision strategies were found, ranging from cognitively-intensive but precise, to fast, frugal but robust. We provide evidence for the model by comparing its behavior to a NASA Ames Research Center simulation of Chicago O'Hare surface operations. Decision horizons were highly variable; the model selected the most accurate strategy given time available. We found a signature in the simulation data of the use of globally robust heuristics to cope with short decision horizons as revealed by errors occurring most frequently at atypical taxiway geometries or clearance routes. These data provided empirical support for the model.
Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.
Queuing network models of a computer system operating with a single workload type are presented. Program algorithms are adapted for use on the Texas Instruments SR-52 programmable calculator. By slightly altering the algorithm to process the G and H matrices row by row instead of column by column, six devices and an unlimited job/terminal population could be handled on the SR-52. Techniques are also introduced for handling a simple load dependent server and for studying interactive systems with fixed multiprogramming limits.
The velocity field at locations above the plane of the rotor disk of a generic helicopter configuration has been measured experimentally and modeled using two computational models. The first computational model represents only the portion of the velocity field due to the fuselage, using only the potential due to freestream onset flow. The second model represents all elements of the problem in an incompressible field: the rotor, the rotor wake, and the non-lifting fuselage. The first model is used to assess the fidelity requirement for fuselage paneling. Comparisons are made between a much simplified body shape, an accurate body shape, a body shape with a crude hub/shaft model and wind-tunnel measurements. The second model is used to assess the relative magnitude of the velocity perturbations at the rotor disk due to the effects of the fuselage. A fully interactive time-stepping wake method is used to predict the velocities at the rotor disk. These velocities are computed with and without the fuselage in the model and compared with measured wind-tunnel inflow velocities. The difference between the velocities computed with and without the fuselage is then compared with the perturbation velocity field computed by the first method.
This document contains presentations given at Workshop on Computational Turbulence Modeling held 15-16 Sep. 1993. The purpose of the meeting was to discuss the current status and future development of turbulence modeling in computational fluid dynamics for aerospace propulsion systems. Papers cover the following topics: turbulence modeling activities at the Center for Modeling of Turbulence and Transition (CMOTT); heat transfer and turbomachinery flow physics; aerothermochemistry and computational methods for space systems; computational fluid dynamics and the k-epsilon turbulence model; propulsion systems; and inlet, duct, and nozzle flow.