Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coding standards”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A comprehensive guide to the physics and usage of PYTHIA 8.3

This manual describes the PYTHIA 8.3 event generator, the most recent version of an evolving physics tool used to answer fundamental questions in particle physics. The program is most often used to generate high-energy-physics collision "events", i.e. sets of particles produced in association with the collision of two incoming high-energy particles, but has several uses beyond that. The guiding philosophy is to produce and reproduce properties of experimentally obtained collisions as accurately as possible. The program includes a wide ranges of reactions within and beyond the Standard Model, and extending to heavy ion physics. Emphasis is put on phenomena where strong interactions play a major role. The manual contains both pedagogical and practical components. All included physics models are described in enough detail to allow the user to obtain a cursory overview of used assumptions and approximations, enabling an informed evaluation of the program output. A number of the most central algorithms are described in enough detail that the main results of the program can be reproduced independently, allowing further development of existing models or the addition of new ones. Finally, a chapter dedicated fully to the user is included towards the end, providing pedagogical examples of standard use cases, and a detailed description of a number of external interfaces. The program code, the online manual, and the latest version of this print manual can be found on the PYTHIA web page: https://www.pythia.org/

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

THE AGORA HIGH-RESOLUTION GALAXY SIMULATIONS COMPARISON PROJECT. II. ISOLATED DISK TEST

Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt–Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of ∼3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2018

Optimization of Elastodynamic Finite Integration Technique on Intel Xeon Phi Knights Landing Processors

This work describes the development and optimization of an implementation of an isotropic elastodynamic finite integration technique (EFIT) code for parallelized computation on Intel Knights Landing (KNL) hardware. EFIT is a numerical approach resulting in standard staggered-grid finite difference equations for the elastodynamic equations of motion to simulate bulk waves is solids. The computationally efficient simulation of elastodynamic wave propagation and interactions in aerospace materials is of high-interest in the fields of nondestructive evaluation (NDE) and structural health monitoring (SHM). Ultrasonic inspection uses an ultrasonic signal, generated at the surface of the material/structure via use of a piezoelectric transducer, to propagate sound waves into the material where it interacts with any existing defects, as well as with structural boundaries and any material inhomogeneity. Reflections from defects and boundaries are then measured by a transducer. Realistic ultrasound simulation tools can significantly aid the development and optimization of inspection techniques and can assist in the interpretation of experimental data. The optimization of an elastodynamics simulation code for the KNL Many Integrated Core processor was performed. The optimization focused on data locality and vectorization. Results show that tiling of the data to exploit the cache behavior and allow for significant utilization of the KNL hardware. The MPI implementation allows for a scalable implementation enabling large problems to be simulated. The model results were validated against theoretical dispersion curves to within 2% of the group velocity, and within 0.5% of the phase velocity of the A0 mode. Aggressive use of tiling, threading, and vectorization techniques allowed for dramatically improved time to solution.

William C. Schneck↗
DOE OSTI2016

THE RELATION BETWEEN SOLAR ERUPTION TOPOLOGIES AND OBSERVED FLARE FEATURES. II. DYNAMICAL EVOLUTION

A long-established goal of solar physics is to build understanding of solar eruptions and develop flare and coronal mass ejection (CME) forecasting models. In this paper, we continue our investigation of nonlinear forces free field (NLFFF) models by comparing topological properties of the solutions to the evolution of the flare ribbons. In particular, we show that data-constrained NLFFF models of three erupting sigmoid regions (SOL2010-04-08, SOL2010-08-07, and SOL2012-05-12) built to reproduce the active region magnetic field in the pre-flare state can be rendered unstable and the subsequent sequence of unstable solutions produces quasi-separatrix layers that match the flare ribbon evolution as observed by SDO/AIA. We begin with a best-fit equilibrium model for the pre-flare active region. We then add axial flux to the flux rope in the model to move it across the stability boundary. At this point, the magnetofrictional code no longer converges to an equilibrium solution. The flux rope rises as the solutions are iterated. We interpret the sequence of magnetofrictional steps as an evolution of the active region as the flare/CME begins. The magnetic field solutions at different steps are compared with the flare ribbons. The results are fully consistent with the three-dimensional extension of the standard flare/CME model. Our ability to capture essential topological features of flaring active regions with a non-dynamic magnetofrictional code strongly suggests that the pre-flare, large-scale topological structures are preserved as the flux rope becomes unstable and lifts off.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Field Testing of Compartmentalization Methods for Multifamily Construction

The 2012 International Energy Conservation Code (IECC) has an airtightness requirement of 3 air changes per hour at 50 Pascals test pressure (3 ACH50) for single-family and multifamily construction (in climate zones 3–8). The Leadership in Energy & Environmental Design certification program and ASHRAE Standard 189 have comparable compartmentalization requirements. ASHRAE Standard 62.2 will soon be responsible for all multifamily ventilation requirements (low rise and high rise); it has an exceptionally stringent compartmentalization requirement. These code and program requirements are driving the need for easier and more effective methods of compartmentalization in multifamily buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

A fast Monte Carlo cell-by-cell simulation for radiobiological effects in targeted radionuclide therapy using pre-calculated single-particle track standard DNA damage data

Introduction: We developed a new method that drastically speeds up radiobiological Monte Carlo radiation-track-structure (MC-RTS) calculations on a cell-by-cell basis. Methods: The technique is based on random sampling and superposition of single-particle track (SPT) standard DNA damage (SDD) files from a “pre-calculated” data library, constructed using the RTS code TOPAS-nBio, with “time stamps” manually added to incorporate dose-rate effects. This time-stamped SDD file can then be input into MEDRAS, a mechanistic kinetic model that calculates various radiation-induced biological endpoints, such as DNA double-strand breaks (DSBs), misrepairs and chromosomal aberrations, and cell death. As a benchmark validation of the approach, we calculated the predicted energy-dependent DSB yield and the ratio of direct-to-total DNA damage, both of which agreed with published in vitro experimental data. We subsequently applied the method to perform a superfast cell-by-cell simulation of an experimental in vitro system consisting of neuroendocrine tumor cells uniformly incubated with 177 Lu. Results and discussion: The results for residual DSBs, both at 24 and 48 h post-irradiation, are in line with the published literature values. Our work serves as a proof-of-concept demonstration of the feasibility of a cost-effective “in silico clonogenic cell survival assay” for the computational design and development of radiopharmaceuticals and novel radiotherapy treatments more generally.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS2010

Use of MODIS Terra Imagery to Estimate Surface Water Quality Standards, Using Lake Thonotosassa, Florida, as a Case Study

Lake Thonotosassa is a highly eutrophied lake located in an area with rapidly growing population in the Tampa Bay watershed, Florida. The Florida Administrative Code has designated its use for "recreation, propagation and maintenance of a healthy, well-balanced population of fish and wildlife." Although this lake has been the subject of efforts to improve water quality since 1970, overall water quality has remained below the acceptable state standards, and has a high concentration of nutrients. This condition is of great concern to public health since it has favored episodic blooms of Cyanobacteria. Some Cyanobacterial species release toxins that can reach humans through drinking water, fish consumption, and direct contact with contaminated water. The lake has been historically popular for fishing and water sports, and its overflow water drains into the Hillsborough River, the main supply of municipal water for the City of Tampa, this explains why it has being constantly monitored in situ for water quality by the Environmental Protection Commission of Hillsborough County (EPC). Advances in remote sensing technology, however, open the possibility of facilitating similar types of monitoring in this and similar lakes, further contributing to the implementation of surveillance systems that would benefit not just public health, but also tourism and ecosystems. Although traditional application of this technology to water quality has been focused on much larger coastal water bodies like bays and estuaries, this study evaluates the feasibility of its application on a 46.6 km2 freshwater lake. Using surface reflectance products from Moderate-Resolution Imaging Spectroradiometer (MODIS) Terra, this study evaluates associations between remotely sensed data and in situ data from the EPC. The parameters analyzed are the surface water quality standards used by the State of Florida and general indicators of trophic status.

Moreno, Max J.↗
DOE OSTI2019

Development of a Standard Test Method for Reducing the Uncertainties in Measuring the Capture Efficiency of Range Hoods

Cooking activities are a major source of indoor air pollutants. To control pollutants generated from cooking activities, a range hood is commonly used in residential kitchens. Several building codes require that a range hood be installed in new homes to control pollutants from cooking, and the required airflow rates for range hoods are specified by indoor air quality (IAQ) standards. However, air flow alone does not tell us how much of the cooking pollutants are exhausted by the range hood. A better metric to evaluate range hood IAQ performance is capture efficiency (CE)—the fraction of contaminants emitted during cooking that are exhausted directly to the outside via the range hood. This paper summarizes the development of a range hood CE test method for use in laboratory testing and equipment rating.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2001

Comparison of the NASCAP/GEO, POLAR, SEE Charging Handbook, and NASCAP-2K.1 Spacecraft Charging Codes

The NASA Charging Analyzer Program (NASCAP) spacecraft charging software developed by Maxwell Technologies has been widely used for the past fifteen to twenty years in satellite design and investigation of spacecraft charging related anomalies. Individual versions of the NASCAP software are available for use in low inclination, low Earth orbit environments (NASCAP[LEO) and geostationary orbit environments (NASCAP/GEO). In addition, the Potentials of Large objects in the Auroral Region (POLAR) code is available for use in LEO polar orbit environments. NASCAP/GEO and POLAR were both written in the 1980's using algorithms appropriate for the computers of the time. They solve the Poisson-Vlasov system for currents and densities assuming limited speed and memory of computer systems standard for the day. In addition, use of the charging models requires individual input files that are not readily transported into the various codes to facilitate comparison of results by the user.

Neergaard, Linda E.↗
DOE OSTI2026

fp-tools

The output of a molecular dynamics simulation is a data file containing properties of the motion of the simulated particles as a function of time (e.g., positions, velocities). From this output, scientists can extract various properties which give insight into the physics of the system; for example, the radial distribution function can be computed from the particle positions, which gives insight into how close the system is to melting or freezing. Many publicly available codes to perform molecular dynamics simulations exist and have been well adopted by the scientific community (e.g., LAMMPS, VASP). However, it remains standard practice for scientists to write their own post-processing scripts to extract various properties of interest from the molecular dynamics simulation data. Since in many cases scientists are interested in the same set of “textbook” properties, there is a large duplication of effort in writing these post-processing scripts. Our proposed code, FP-Tools, is a C++ toolkit designed to help scientists extract and analyze quantities of interest from the output of a molecular dynamics simulation, eliminating the need for them to write these codes themselves and thus reducing this duplication of effort. The quantities which are computed by our code are well documented in the literature (either in publications or textbooks), and the algorithms we implement to compute these properties are also well known in the field. We are not introducing new science or methods here; rather, our primary goal is to provide a useful tool to the community in the form of a well-documented, easy to use package.

Hartman, Leah↗
DOE OSTI2016

Errors from approximation of ODE systems with reduced order models

This is a code to calculate the error from approximation of systems of ordinary differential equations (ODEs) by using Proper Orthogonal Decomposition (POD) Reduced Order Models (ROM) methods and to compare and analyze the errors for two POD ROM variants. The first variant is the standard POD ROM, the second variant is a modification of the method using the values of the time derivatives (a.k.a. time-derivative snapshots). The code compares the errors from the two variants under different conditions.

Vassilevska, Tanya↗
DOE OSTI2022

Deep Learning of Dark Energy Spectroscopic Instrument Mock Spectra to Find Damped Lyα Systems

We have updated and applied a convolutional neural network (CNN) machine-learning model to discover and characterize damped Ly α systems (DLAs) based on Dark Energy Spectroscopic Instrument (DESI) mock spectra. We have optimized the training process and constructed a CNN model that yields a DLA classification accuracy above 99% for spectra that have signal-to-noise ratios (S/N) above 5 per pixel. The classification accuracy is the rate of correct classifications. This accuracy remains above 97% for lower S/N ≈1 spectra. This CNN model provides estimations for redshift and H i column density with standard deviations of 0.002 and 0.17 dex for spectra with S/N above 3 pixel -1 . Also, this DLA finder is able to identify overlapping DLAs and sub-DLAs. Further, the impact of different DLA catalogs on the measurement of baryon acoustic oscillations (BAO) is investigated. The cosmological fitting parameter result for BAO has less than 0.61% difference compared to analysis of the mock results with perfect knowledge of DLAs. This difference is lower than the statistical error for the first year estimated from the mock spectra: above 1.7%. We also compared the performances of the CNN and Gaussian Process (GP) models. Our improved CNN model has moderately 14% higher purity and 7% higher completeness than an older version of the GP code, for S/N > 3. Both codes provide good DLA redshift estimates, but the GP produces a better column density estimate by 24% less standard deviation. A credible DLA catalog for the DESI main survey can be provided by combining these two algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Coupling of CTF and TRANSFORM using the Functional Mockup Interface

An in-memory coupling between the sub-channel thermal hydraulics code COBRA-TF (CTF), which is included in the Virtual Environment for Reactor Applications (VERA), and the systems code Transient Simulation Framework of Reconfigurable Models (TRANSFORM) was developed in this work. Data exchange is accomplished by using the Functional Mock-Up Interface (FMI), an open standard for coupling dynamic systems models together. The FMI-based coupling necessitated the development of a novel FORTRAN wrapper for communicating with Functional Mock-Up Units. This wrapper will facilitate future FMI-based code couplings with other FORTRAN-based programs. CTF-TRANSFORM is exercised on a simplified Molten Salt Reactor Experiment (MSRE) model in steady-state and transient configurations. The coupled CTF-TRANSFORM model is shown to predict core temperature deltas similar to those available in historical MSRE operational data and is shown to be robust to fast power transients.The coupling paves the way to performing sensitivity and design studies impossible with VERA/CTF alone, in which secondary and tertiary loop operating conditions and design parameters may be perturbed and their impact on the core operating conditions studied.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Propagation Analysis of the 3rd Sonic Boom Prediction Workshop Cases using sBOOM

The 3rd Sonic Boom Prediction Workshop was held on January 4-5, 2020 in conjunction with the AIAA SciTech conference. The workshop had 23 participants who represented 6 different countries and 18 unique organizations that included government, industry, and academia. The workshop was attended by over 50 individuals. The motivation for the workshop stems from the goal of obtaining supersonic commercial overland flight. In order to replace the current prohibition with a certification standard, an international effort is required to quantify the accuracy and reliability of prediction methods. The workshops also identify deficiencies in existing methods where further research should be focused. The workshop consisted of two days with the first day focused on near-field computational fluid dynamics (CFD) and the second day on propagation techniques. For the atmospheric propagation portion of the workshop, two required cases and one optional case were prescribed for the participants to exercise their propagation implementations. This paper presents and discusses the results that were obtained using the sonic boom propagation code sBOOM. Results show that sonic boom carpets could be much different under arbitrary atmospheric conditions as opposed to standard atmospheric conditions. Additionally, several unique features when using sBOOM are presented.

Sonic Boom↗
DOE OSTI2014

A Code of Radiation Damage by Gamma - Paper 62

In recent years, the study of radiation damage of materials has seen a thrust because of the expectation of the high rate of radiation damage in fusion reactors, accelerator driven sub critical systems and damage of a control system of space vehicles by the high energy neutrons and gamma showers. Embrittlement of reactor vessels, corrosion of steel and swelling of fuel rods are the well known empirical examples of severe radiation damage. It is also realized that damage by high energy gammas is in no way insignificant, and it is equivalent to the damage by the neutrons. Secondly, value addition of certain materials like enhancement of electrical resistance, mechanical strength, color enhancement of precious stones on irradiation by neutrons and gammas and polymerization in the presence of gamma radiation add to the necessity of a Monte Carlo (MC) code of radiation damage in a wide range of energy. The JAIPU code, for the study of radiation damage by neutrons up to 10 MeV energy, incorporates features like development of secondary atomic cascade and PKA energy corrected for the ionization lose and production of Frenkel pair. MC estimation of both the primary and secondary displacements has been compared with the Frenkel pairs calculated using the NRT model. In the present work, JAIPU code is extended for irradiation of materials by high energy gammas. Here cascade generated by the incident gamma produces electrons by way of photo-electric, Compton and pair production until it is absorbed. The recoiled Frenkel pairs produced in electron + atom collision have been considered for the purpose of atomic cascade. Standard cross section libraries for gamma and electron interaction are used while for atom + atom collision we have made use of the IOTA code. In addition, in the code, the lattice structure of the irradiated material is simulated to know the positions of the photon interactions with atoms and the energy of electron it adds an extra feature of the position and size of the voids and the clusters produced on irradiation. The code has the possibility of calculation of 'absorbed radiation dose' and 'KERMA' and the interatomic distances of the displaced atoms. It may be pointed out that the code uses only the MC simulation of the number of displacements unlike several other approaches where the NRT model with PKA energy is used for gross estimation of the number of the displaced atoms. Thus, the code keeps an account of the positions of displaced atoms for the in-depth modification of the irradiated material. 1 MeV gamma seems to damage Fe up to several millimeters of its depth and this is mostly by the recoiled Compton electrons than photo electrons. The code is in the developmental stage and it has the capability of applications in the developing field of Radiology besides the shielding and other studies of radiation damage. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2018

Modeling of Time-Dependent Strength Degradation of SiC/SiC Ceramic Matrix Composites via MAC/GMC Computer Code

Silicon carbide fiber reinforced silicon carbide (SiC/SiC) ceramic matrix composites (CMCs) display time-dependent strength degradation at intermediate temperatures (600 to 900 degrees Centigrade). This is generally believed to be an oxidation induced phenomenon. The understanding of the effect of temperature with environment (oxidation) is key towards development of SiC/SiC CMCs with a reliable load carrying capacity. Various theories have been proposed to explain the strength degradation. One suggests that the boron nitride (BN) coating deposited on the fibers oxidizes causing fusion of fibers. Another theory proposes that the SiC fibers are oxidized forming a silica scale leading to premature fiber failure. A more recent theory suggests that SiC fiber strength is intrinsically time-dependent due to slow crack growth of flaws in the fibers. An empirical model, termed as a "fiber classic model," which is based on a standard slow crack growth type power-law, has been implemented within NASA's micromechanics-based MAC/GMC computer code as a user routine. Model parameters for this "classic model" were calibrated from stress-rupture data of Hi-Nicalon (TM) monofilaments using the maximum likelihood estimation (MLE) technique. This new capability in the MAC/GMC computer code was then used to predict the stress-rupture behavior of Hi-Nicalon (TM) tows as well as 2-D SiC/SiC composites reinforced with Hi-Nicalon (TM) fibers. Results demonstrate that the MAC/GMC with this new capability successfully predicts the time-to-failure vs. applied stress within the intermediate temperature range at various scales as well as laminated composites in an oxidizing environment.

Environmental degradation↗
DOE OSTI2016

Zoomed high-resolution simulations of Multi-coupled Dark Energy: cored galaxy density profiles at high redshift

We perform for the first time high-resolution zoom-in re-simulations of individual halos in the context of the Multi-coupled Dark Energy (McDE) scenario, which is characterised by the existence of two distinct Dark Matter particle species with opposite couplings to a Dark Energy scalar field. We compare the structural properties of the simulated halos to the standard ΛCDM results. The zoomed-in initial conditions are set up using a specifically designed code called ZInCo that we publicly release along with the present paper. Our numerical results allow to investigate in detail and with unprecedented resolution the halo segregation process that characterises McDE cosmologies from its very early stages. In particular, we find that in contrast to what could be inferred from previous numerical analysis at lower resolution, the segregation process is already in place at redshifts as high as z ∼ 7. Most remarkably, we find that the subsequent evolution of the segregation leads to the formation of cored total matter density profiles with a core size that progressively increases in time. The shape of the cored profiles can be accurately predicted as the superposition of two NFW profiles with an increasing offset, thereby confirming the interpretation of the simulations results in terms of the segregation of the two dark matter components of the halo as a consequence of their different coupling to the Dark Energy field.

79 ASTRONOMY AND ASTROPHYSICS↗