Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer programming”

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 289 records · Page 16

Center for Extended Magnetohydrodynamic Modeling (CEMM), Final Technical Report for UW-Madison Engineering Physics

The Center for Extended Magnetohydrodynamic Modeling (CEMM) focused on the further development, verification, validation, and application of two scientific computation programs, the NIMROD and M3D-C1 codes, for modeling macroscopic dynamics in magnetized plasma. Both NIMROD and M3D-C1 solve the mathematical equations that model two-fluid 3D magnetohydrodynamic (MHD). They are designed and optimized to be efficient and accurate for toroidal confinement physics calculations. The effort improved the fidelity of the physics models, the realism of the treatment of the surrounding vessel, coils, and other structures, the accuracy and efficiency of the algorithms, and the ability to scale calculations to thousands of processors in order to take advantage of leading-edge computers. The primary application areas of Resonant Magnetic Perturbations (RMP), Edge Harmonic Oscillations (EHO), Edge Localized Modes (ELMs), Disruption Studies, Sawteeth and Stationary States, and Kinetic MHD and Neoclassical Tearing Modes (NTM) are of concern for the international ITER experiment, which is under construction in France. The effort by the UW-Madison Engineering Physics component of CEMM improved the numerical methods in the NIMROD code with respect to modeling MHD interchange and two-fluid effects. It also improved the algorithm for solving large algebraic systems to enhance NIMROD's computational performance when running the code in parallel on increasing numbers of processor cores. The modeling efforts included a study of the nonlinear evolution of ballooning, and the findings were supported by analytical theory and computation. Numerical MHD modeling also examined the influences of two-fluid effects on ELMs and RMPs, where "two-fluid" refers to separation of electron dynamics from ion dynamics. Another related aspect considered the effects of a localized layer of electrical current density near the edge of the plasma confinement region. Two-fluid effects on resonant reconnecting instabilities were verified and applied to computations of sawtooth and drift-tearing dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Center for Extended Magnetohydrodynamic Modeling (CEMM, UW-Madison Physics)

The Center for Extended Magnetohydrodynamic Modeling (CEMM) focused on the further development, verification, validation, and application of two scientific computation programs, the NIMROD and M3D-C1 codes, for modeling macroscopic dynamics in magnetized plasma. Both NIMROD and M3D-C1 solve the mathematical equations that model two-fluid 3D magnetohydrodynamic (MHD). They are designed and optimized to be efficient and accurate for toroidal confinement physics calculations. The effort improved the fidelity of the physics models, the realism of the treatment of the surrounding vessel, coils, and other structures, the accuracy and efficiency of the algorithms, and the ability to scale calculations to thousands of processors in order to take advantage of leading-edge computers. The primary application areas of Resonant Magnetic Perturbations (RMP), Edge Harmonic Oscillations (EHO), Edge Localized Modes (ELMs), Disruption Studies, Sawteeth and Stationary States, and Kinetic MHD and Neoclassical Tearing Modes (NTM) are of concern for the international ITER experiment, which is under construction in France. The component of CEMM that is described in this report applied recently developed aspects of NIMROD to sawtooth and EHO dynamics, which occur in tokamak plasma confinement systems. In particular, the effort researched energetic-particle, two-fluid, and plasma-flow effects. The work included an effort to validate energetic-particle modeling using laboratory measurements and previous analysis of giant sawtooth oscillations, where energetic particles temporarily stabilize the system. A side effort that grew out of this work was benchmarking of the slab-geometry ion temperature gradient (ITG) instability computed by NIMROD's two-fluid modeling against analytical results. These fluid-model results were also compared with results from kinetic modeling to help quantify the range of validity for the two-fluid model. The EHO simulations show that sheared plasma flow plays a crucial role in saturating the underlying MHD instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PEBBLES Operation and Theory Manual

The PEBBLES manual describes the PEBBLES code. The PEBBLES code is a computer program designed to simulation the motion, packing and vibration of spheres that undergo various mechanical forces including gravitation, Hooke’s law force and various friction forces. The frictional forces include true static friction that allows non-zero angles of repose. Each pebble is individually simulated using the distinct element method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Determination of time zero from a charged particle detector

A method, system and computer program is used to determine a linear track having a good fit to a most likely or expected path of charged particle passing through a charged particle detector having a plurality of drift cells. Hit signals from the charged particle detector are associated with a particular charged particle track. An initial estimate of time zero is made from these hit signals and linear tracks are then fit to drift radii for each particular time-zero estimate. The linear track having the best fit is then searched and selected and errors in fit and tracking parameters computed. The use of large and expensive fast detectors needed to time zero in the charged particle detectors can be avoided by adopting this method and system.

47 OTHER INSTRUMENTATION↗

Radiation shielding properties of barite coated fabric by computer programme

With the development of technology radiation started to be used in variety of different fields. As the radiation is hazardous for human health, it is important to keep radiation dose as low as possible. This is done mainly using shielding materials. Barite is one of the important materials in this purpose. As the barite is not used directly it can be used in some other materials such as fabric. For this purposes barite has been coated on fabric in order to improve radiation shielding properties of fabric. Determination of radiation shielding properties of coated fabric has been done by using computer program written C# language. With this program the images obtained from digital Rontgen films is used to determine radiation shielding properties in terms of image processing numerical values. Those values define radiation shielding and in this way the coated barite effect on radiation shielding properties of fabric has been obtained.

07 ISOTOPE AND RADIATION SOURCES↗

Clustering procedures for the optimal selection of data sets from multiple crystals in macromolecular crystallography

A systematic approach to the scaling and merging of data from multiple crystals in macromolecular crystallography is introduced and explained. The availability of intense microbeam macromolecular crystallography beamlines at third-generation synchrotron sources has enabled data collection and structure solution from microcrystals of <10 µm in size. The increased likelihood of severe radiation damage where microcrystals or particularly sensitive crystals are used forces crystallographers to acquire large numbers of data sets from many crystals of the same protein structure. The associated analysis and merging of multi-crystal data is currently a manual and time-consuming step. Here, a computer program, BLEND, that has been written to assist with and automate many of the steps in this process is described. It is demonstrated how BLEND has successfully been used in the solution of a novel membrane protein.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multigroup computation of the temperature-dependent Resonance Scattering Model (RSM) and its implementation

A multi-group formulation for the exact neutron elastic scattering kernel is developed. This formulation is intended for implementation into a lattice physics code. The correct accounting for the crystal lattice effects influences the estimated values for the probability of neutron absorption and scattering, which in turn affect the estimation of core reactivity and burnup characteristics. A computer program has been written to test the formulation for various nuclides. Results of the multi-group code have been verified against the correct analytic scattering kernel. In both cases neutrons were started at various energies and temperatures and the corresponding scattering kernels were tallied. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Surface Modification of ICF Target Capsules by Pulsed Laser Ablation

Topographical modifications of spherical surfaces are imprinted on National Ignition Facility (NIF) target capsules by extending the capabilities of a recently developed full surface (4π) laser ablation and mapping apparatus. The laser ablation method combines the precision, energy density and long reach of a focused laser beam to pre-impose sinusoidal modulations on the outside surface of High Density Carbon (HDC) capsules and the inside surface of Glow Discharge Polymer (GDP) capsules. Sinusoidal modulations described in this paper have sub-micron to 10’s of microns vertical scale and wavelengths as small as 30 μm and as large as 200 μm. The modulated patterns are created by rastering a focused laser fired at discrete capsule surface locations for a specified number of pulses. The computer program developed to create these raster patterns uses inputs such as laser beam intensity profile, the material removal function, the starting surface figure and the desired surface figure. The patterns are optimized to minimize surface roughness. Lastly, in this paper, simulated surfaces are compared with actual ablated surfaces measured using confocal microscopy.

36 MATERIALS SCIENCE↗

NREL Develops Diagnostic Test Cases to Improve Building Energy Simulation Programs (Fact Sheet)

This technical highlight describes NREL research to develop a set of diagnostic test cases for building energy simulations in order to achieve more accurate energy use and savings predictions. The National Renewable Energy Laboratory (NREL) Residential and Commercial Buildings research groups developed a set of diagnostic test cases for building energy simulations. Eight test cases were developed to test surface conduction heat transfer algorithms of building envelopes in building energy simulation programs. These algorithms are used to predict energy flow through external opaque surfaces such as walls, ceilings, and floors. The test cases consist of analytical and vetted numerical heat transfer solutions that have been available for decades, which increases confidence in test results. NREL researchers adapted these solutions for comparisons with building energy simulation results. Testing the new cases with EnergyPlus identified issues with the conduction finite difference (CondFD) heat transfer algorithm in versions 5 and 6. NREL researchers resolved these issues for EnergyPlus version 7. The new test cases will help users and developers of EnergyPlus and other building energy tools to identify and fix problems associated with solid conduction heat transfer algorithms of building envelopes and their boundary conditions. In the long term, improvements to software algorithms will result in more accurate energy use and savings predictions. NREL researchers plan to document the set of test cases and make them available for future consideration by validation standards such as ASHRAE Standard 140: Standard Method of Test for the Evaluation of Building Energy Analysis Computer Programs. EnergyPlus users will also have access to the improved CondFD model in version 7 after its next scheduled release.

14 SOLAR ENERGY↗

NLOX, a one-loop provider for Standard Model processes

NLOX is a computer program for calculations in high-energy particle physics. It provides fully renormalized scattering matrix elements in the Standard Model of particle physics, up to one-loop accuracy for all possible coupling-power combinations in the strong and electroweak couplings, and for processes with up to six external particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Beam Based RF Voltage Measurements and Longitudinal Beam Tomography at the Fermilab Booster

Increasing proton beam power on neutrino production targets is one of the major goals of the Fermilab long term accelerator programs. In this effort, the Fermilab 8 GeV Booster synchrotron plays a critical role for at least the next two decades. Therefore, understanding the Booster in great detail is important as we continue to improve its performance. For example, it is important to know accurately the available RF power in the Booster by carrying out beam-based measurements in order to specify the needed upgrades to the Booster RF system. Since the Booster magnetic field is changing continuously measuring/calibrating the RF voltage is not a trivial task. Here, we present a beam based method for the RF voltage measurements. Data analysis is carried out using computer programs developed in Python and MATLAB. The method presented here is applicable to any RCS which do not have flat-bottom and flat-top in the acceleration magnetic ramps. We have also carried out longitudinal beam tomography at injection and extraction energies with the data used for RF voltage measurements. Beam based RF voltage measurements and beam tomography were never done before for the Fermilab Booster. The results from these investigations will be very useful in future intensity upgrades.

43 PARTICLE ACCELERATORS↗

Method and apparatus for real-time measurement of fuel gas compositions and heating values

An exemplary embodiment can be an apparatus for real-time, in situ measurement of gas compositions and heating values. The apparatus includes a near infrared sensor for measuring concentrations of hydrocarbons and carbon dioxide, a mid infrared sensor for measuring concentrations of carbon monoxide and a semiconductor based sensor for measuring concentrations of hydrogen gas. A data processor having a computer program for reducing the effects of cross-sensitivities of the sensors to components other than target components of the sensors is also included. Also provided are corresponding or associated methods for real-time, in situ determination of a composition and heating value of a fuel gas.

03 NATURAL GAS↗

Uncertainty Analysis of RELAP5-3D©

As world-wide energy consumption continues to increase, so does the demand for the use of alternative energy sources, such as Nuclear Energy. Nuclear Power Plants currently supply over 370 gigawatts of electricity, and more than 60 new nuclear reactors have been commissioned by 15 different countries. The primary concern for Nuclear Power Plant operation and lisencing has been safety. The safety of the operation of Nuclear Power Plants is no simple matter- it involves the training of operators, design of the reactor, as well as equipment and design upgrades throughout the lifetime of the reactor, etc. To safely design, operate, and understand nuclear power plants, industry and government alike have relied upon the use of best-estimate simulation codes, which allow for an accurate model of any given plant to be created with well-defined margins of safety. The most widely used of these best-estimate simulation codes in the Nuclear Power industry is RELAP5-3D. Our project focused on improving the modeling capabilities of RELAP5-3D by developing uncertainty estimates for its calculations. This work involved analyzing high, medium, and low ranked phenomena from an INL PIRT on a small break Loss-Of-Coolant Accident as wall as an analysis of a large break Loss-Of- Coolant Accident. Statistical analyses were performed using correlation coefficients. To perform the studies, computer programs were written that modify a template RELAP5 input deck to produce one deck for each combination of key input parameters. Python scripting enabled the running of the generated input files with RELAP5-3D on INL’s massively parallel cluster system. Data from the studies was collected and analyzed with SAS. A summary of the results of our studies are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

TEAM project review, year 1

This note summarizes our research activities within the TEAM project between October 2019 and November 2020, funded by the ASCR Advanced Research in Quantum Computing program. Three research efforts have begun during the review period. The first effort, resulting in the JuQBox software, implemented a novel algorithm for optimal control in closed quantum systems based on symplectic time integration and exact calculations of the gradient of the discretized objective function. In our second effort, resulting in a code called Quandary, we consider the optimal control problem for open quantum systems. Here the evolution of the quantum system is modeled by the Lindblad master equation. In this code, we are considering both the realization of quantum gates as well as the optimal reset problem in which a general mixed quantum state is driven towards its ground state. In our third research effort we have begun the development of computational software for what we call Lindblad learning, i.e. the process of characterizing the equations governing our quantum computing hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Automated Blazar Light Curves Using Machine Learning

Every night in a remote clearing called Fenton Hill high in the Jemez Mountains of central New Mexico, a bank of robotically controlled telescopes tilt their lenses to the sky for another round of observation through digital imaging. Los Alamos National Laboratory’s Thinking Telescopes project is watching for celestial transients including high-power cosmic flashes called, and like all science, it can be messy work. To keep the project clicking along, Los Alamos scientists routinely install equipment upgrades, maintain the site, and refine the sophisticated machinelearning computer programs that process those images and extract useful data from them. Each week the system amasses 100,000 digital images of the heavens, some of which are compromised by clouds, wind gusts, focus problems, and so on. For a graduate student at the Lab taking a year’s break between master’s and Ph.D. studies, working with state-of-the-art autonomous telescopes that can make fundamental discoveries feels light years beyond the classroom.

79 ASTRONOMY AND ASTROPHYSICS↗

Layer-by-Layer Analysis of the Thickness Distribution of Silicon Dioxide in the Structure SiO{sub 2}/Si(111) by Inelastic Electron Scattering Cross-Section Spectroscopy

The SiO{sub 2}-concentration profile in the structure SiO{sub 2}/Si(111) is determined from experimental spectra of the cross section for the inelastic scattering of reflected electrons in the primary electron energy range from 300 to 3000 eV. The spectra are analyzed with the use of a proposed algorithm and developed computer program for simulating the spectra of the cross section for the inelastic scattering of reflected electrons for layered structures with an arbitrary number of layers, arbitrary thickness, and variable concentration of components in each layer. The best agreement between the calculated and experimental spectra is attained by varying the silicon dioxide and silicon concentrations in each layer. The results obtained can be used for profiling film-substrate structures with an arbitrary composition.

36 MATERIALS SCIENCE↗

Implementation of new PRODIAG algorithm and simulation-based acceptance test - 219

This paper presents the improvement of an automated-reasoning computer program for nuclear power plant diagnosis, namely PRODIAG. PRODIAG, first developed at the Argonne National Laboratory, uses confluence-based reasoning to perform fault diagnosis in engineering systems that can be described by the conservation laws for mass, energy, and momentum. Therefore, unlike other plant configuration-dependent rule-based diagnosis algorithms, the code is generally applicable to a wide variety of plant configurations by simply replacing the plant instrumentation diagram (PID) database. The need for improvement was raised when the original PRODIAG failed to provide the requisite degree of extensibility and maintainability for new applications. These shortcomings have been addressed through a new code architecture that uses an objected-oriented logic-programming language in place of the original Prolog language and uses an open-source automated reasoning engine in place of the original hardwired rule-evaluation sequence. This paper presents the architecture and implementation details of the code transition. Verification and validation studies were performed for the new PRODIAG, which involved re-running the distinct fault cases defined at the time the code was originally developed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Spectrometer of Synchrotron Radiation Based on Diffraction Focusing a Divergent Beam Formed by a Compound Refractive Lens

The results of the first experimental realization of a spectrometer based on the effect of diffraction focusing of X rays by a flat single crystal are discussed. A secondary X-ray source with a relatively high angular divergence and small sizes was formed at the focus of a compound refractive lens having 50 beryllium biconcave elements with a curvature radius of 50 μm. The silicon spectrometer crystal was cut in the form of a wedge of variable thickness, oriented perpendicular to the diffraction plane. The reflection 111 was used for energies of 8.3 and 12 keV. To simulate the experiment, a computer program was developed, which takes into account accurately and for the first time the focusing of radiation by the lens and its subsequent diffraction in the crystal. A calculated curve for a monochromatic beam has made it possible to determine the monochromator spectrum with high resolution from experimental data for a polychromatic beam. It is shown that monochromator resolution increases with an increase in the distance from the compound refractive lens to the crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗