Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source code”

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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Orb5: A global electromagnetic gyrokinetic code using the PIC approach in toroidal geometry

This paper presents the current state of the global gyrokinetic code Orb5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The Orb5 code solves the electromagnetic VlasovMaxwell system of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization density and drift-kinetic electrons. Variants of the physical model can be selected for electrons such as assuming an adiabatic response or a “hybrid” model in which passing electrons are assumed adiabatic and trapped electrons are drift-kinetic. A Fourier filter as well as various control variates and noise reduction techniques enable simulations with good signal-to-noise ratios at a limited numerical cost. They are completed with different momentum and zonal flow-conserving heat sources allowing for temperature-gradient and flux-driven simulations. The code, which runs on both CPUs and GPUs, is well benchmarked against other similar codes and analytical predictions, and shows good scalability up to thousands of nodes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A NUMERICAL SCHEME FOR SPECIAL RELATIVISTIC RADIATION MAGNETOHYDRODYNAMICS BASED ON SOLVING THE TIME-DEPENDENT RADIATIVE TRANSFER EQUATION

We develop a numerical scheme for solving the equations of fully special relativistic, radiation magnetohydrodynamics (MHDs), in which the frequency-integrated, time-dependent radiation transfer equation is solved to calculate the specific intensity. The radiation energy density, the radiation flux, and the radiation stress tensor are obtained by the angular quadrature of the intensity. In the present method, conservation of total mass, momentum, and energy of the radiation magnetofluids is guaranteed. We treat not only the isotropic scattering but also the Thomson scattering. The numerical method of MHDs is the same as that of our previous work. The advection terms are explicitly solved, and the source terms, which describe the gas–radiation interaction, are implicitly integrated. Our code is suitable for massive parallel computing. We present that our code shows reasonable results in some numerical tests for propagating radiation and radiation hydrodynamics. Particularly, the correct solution is given even in the optically very thin or moderately thin regimes, and the special relativistic effects are nicely reproduced.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Validation of the intrinsic spatial efficiency method for non cylindrical homogeneous sources using MC simulation

The Monte Carlo simulation of the gamma spectroscopy systems is a common practice in these days. The most popular softwares to do this are MCNP and Geant4 codes. The intrinsic spatial efficiency method is a general and absolute method to determine the absolute efficiency of a spectroscopy system for any extended sources, but this was only demonstrated experimentally for cylindrical sources. Due to the difficulty that the preparation of sources with any shape represents, the simplest way to do this is by the simulation of the spectroscopy system and the source. In this work we present the validation of the intrinsic spatial efficiency method for sources with different geometries and for photons with an energy of 661.65 keV. In the simulation the matrix effects (the auto-attenuation effect) are not considered, therefore these results are only preliminaries. The MC simulation is carried out using the FLUKA code and the absolute efficiency of the detector is determined using two methods: the statistical count of Full Energy Peak (FEP) area (traditional method) and the intrinsic spatial efficiency method. The obtained results show total agreement between the absolute efficiencies determined by the traditional method and the intrinsic spatial efficiency method. The relative bias is lesser than 1% in all cases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Numerical study of atomic production rate in hydrogen negative ion sources with the effect of non-equilibrium electron energy distribution function

Spatial non-uniformity of the dissociative hydrogen atom (H{sup 0}) production has been investigated in a large negative ion source (JAEA 10 A source) with the electron energy distribution function (EEDF) obtained by a Monte-Carlo simulation code for electron transport in 3D3V (three dimensional real and velocity) space. It has been shown that the H{sup 0} production rate becomes larger in the upper region (one side in the longitudinal direction) of the source chamber. This spatial non-uniformity of the H{sup 0} production profile is mainly explained by the non-equilibrium features of the EEDF in the upper region, i.e., the EEDF consists of thermal electron component with kinetic energy ε < 25 eV and fast electron component with energy ε > 25 eV in the upper region, while the EEDF mainly consists of only thermal electrons in the bottom region. These characteristics for the EEDF and the energy dependence of cross-sections for dissociation and dissociative ionization processes lead to the non-uniform profile of the H{sup 0} production. The above numerical results of the spatial H{sup 0} non-uniformity are validated and confirmed by comparisons with those by spectroscopic measurement. It has been clarified that the non-equilibrium (fast electron) component of the EEDF has a large contribution to the non-uniformity of the H{sup 0} production rate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

The resolved outflow from 3C 48

We investigate the properties of the high-velocity outflow driven by the young radio jet of 3C 48, a compact-steep-spectrum source. We use the Space Telescope Imaging Spectrograph on board the Hubble Space Telecope to obtain (1) low-resolution UV and optical spectra and (2) multi-slit medium-resolution spectra of the ionized outflow. With supporting data from ground-based spectrographs, we are able to accurately measure the ratios of diagnostic emission lines such as [O III] λ5007, [O III] λ3727, [N II] λ6548, Hα, Hβ, [Ne V] λ3425, and [Ne III] λ3869. We fit the observed emission-line ratios using a range of ionization models, powered by active galactic nucleus (AGN) radiation and shocks, produced by the MAPPINGS code. We have determined that AGN radiation is likely the dominant ionization source. The outflow's density is estimated to be in the range n = 10{sup 3}-10{sup 4} cm{sup –3}, the mass is ∼6 × 10{sup 6} M {sub ☉}, and the metallicity is likely equal to or higher than solar. Compared with the typical outflows associated with more evolved radio jets, this young outflow is denser, less massive, and more metal rich. Multi-slit observations allow us to construct a two-dimensional velocity map of the outflow that shows a wide range of velocities with distinct velocity components, suggesting a wide-angle clumpy outflow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

TU-AB-201-09: Calibration of An Element of a New Directional Pd-103 Planar Source Array

Purpose: The CivaSheet™ is a new directional Pd-103 planar source array, with a variable number of discrete source elements referred to as “dots”. Each dot consists of a polymer capsule containing {sup 103}Pd and a gold shield that attenuates radiation on one side of the device to define hot and cold dose regions. Fluorescence from the gold shield is observed in the dot spectrum. Since CivaSheet™ is a planar directional source, conventional methods used for calibration of azimuthally symmetric sources are not applicable. The purpose of this work is to establish an air-kerma-strength standard and a transfer to a well chamber for clinical calibration. Methods: Primary air-kerma strength measurement of the dots was performed using a variable-aperture free-air chamber (VAFAC). Charge measurements were recorded using a well chamber with a custom insert. Anisotropy measurements were performed using a Sodium-Iodide detector. Spectral measurements were performed using a low-energy germanium detector and compared to a source without gold. The dot geometry was modeled using the MCNP6 radiation transport code. Results: Air-kerma strength measurements of a batch of four dots performed with the VAFAC were within ±1.5% of the average measured value and the measurement precision was within ±0.5%. Anisotropy measurements indicated uniform emission within the measurement uncertainty for the solid angle defining the VAFAC aperture used. Charge measurements of each dot using the well chamber in four cardinal angle source orientations were within ±1.5% of the average measured values. The spectral study of a dot resulted in identification of fluorescence from the gold shield and primary spectral energies that were compared to MCNP6 simulations. Conclusion: Calibration procedures for the new directional Pd-103 source dot were established for future clinical use, based on the results of experimental and Monte Carlo investigations. This work was partially supported by NCI contract (HHSN261201200052C) through CivaTech Oncology Inc.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Integrated, Interoperable Software Environment for Fusion Simulation and Data Analysis Tools SBIR Phase II

This SBIR effort was focused on developing a production ready system to address the integration and interoperability challenges with analysis and visualization in fusion simulations. Our overarching technical objective was to minimize the code development simulation scientists incur when coupling their simulation codes with different analysis frameworks. To this end, we developed an open-source software library to make data exchange between application easier and a Web application to manage and display analysis extracts from simulations. We have also augmented existing libraries funded by DOE such as ADIOS and VTK-m. When used together, these make it significantly easier to integrate simulation and analysis capability. We demonstrated the flexibility of our approach using two common simulation codes in the fusion community, XGC1 and GTC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Thermochimica

Thermochimica is a software library that determines a unique combination of phases and their compositions at thermochemical equilibrium. Thermochimica can be used for stand-alone calculations or it can be directly coupled to other codes. This release of the software does not have a graphical user interface (GUI) and it can be executed from the command line or from an Application Programming Interface (API). Also, it is not intended for thermodynamic model development or for constructing phase diagrams. The main purpose of the software is to be directly coupled with a multi-physics code to provide material properties and boundary conditions for various physical phenomena. Significant research efforts have been dedicated to enhance computational performance through advanced algorithm development, such as improved estimation techniques and non-linear solvers. Various useful parameters can be provided as output from Thermochimica, such as: determination of which phases are stable at equilibrium, the mass of solution species and phases at equilibrium, mole fractions of solution phase constituents, thermochemical activities (which are related to partial pressures for gaseous species), chemical potentials of solution species and phases, and integral Gibbs energy (referenced relative to standard state). The overall goal is to provide an open source computational tool to enhance the predictive capability of multi-physics codes without significantly impeding computational performance.

Simunovic, Srdjan↗
DOE OSTI2024

Shielding Benchmark Comparison - MCNP6.2

This report documents the calculations performed for a shielding code comparison between various Department of Energy (DOE) sites. These shielding calculations were performed using an experiment drawn from the International Handbook of Evaluated Criticality Safety Benchmark Experiments, published by the Organisation for Economic Cooperation and Development/Nuclear Energy Agency (OECD/NEA). The benchmark selected for comparison is ALARM-CF-AIR-LAB-001 (“Neutron Fields in the Three-Section Concrete Labyrinth from Cf-252 Source, Benchmark ALARM CF AIR LAB-001”). The Y-12 submission for this code comparison was performed with Monte Carlo N Particle (MCNP) Transport Code System, Version 6.2 and Automated Variance Reduction Generator (ADVANTG).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Model-assisted decision analyses related to a chromium plume at Los Alamos national Laboratory - 15449

?he chromium plume at Los Alamos National Laboratory is a complex contaminant-remediation site with spatial extent of several square kilometers. It is located in a regional aquifer near water-supply wells and other points of compliance. The plume originated on the ground surface and the contaminant migrated along complex pathways through a thick vadose zone (∼300 m) that includes intermediate lateral zones of saturation and vertical preferential flow paths. Here we discuss the development of a modeling computational framework for the site. The framework work incorporates data assimilation and model-results analysis tools including decision support techniques. The decision tools are demonstrated using synthetic problems consistent with the interpretation of existing site data and site modeling results. The analyses aim at robust and efficient environmental management of the plume and explore a series of alternative remedial options. The remedial options include natural attenuation (NA), enhanced attenuation (EA), contaminant source removal (in the vadose zone), contaminant extraction (in the regional aquifer at the plume centroid and peripheries), biogeochemical remediation (injection of fluids stimulating growth of organisms in the aquifer impacting chromium concentrations), as wells as hydraulic controls on the groundwater flow and transport in the vadose zone and the regional aquifer. The applied model analyses and decision-support algorithms are implemented in code MADS (Model Analyses for Decision Support; http://mads.lanl.gov). MADS is an open-source framework for model-based decision support employing system and physics simulation models. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Design of the new extraction grid for the NIO1 negative ion source

NIO1 is a compact source of negative ions jointly developed by RFX and INFN, to study the physics of production and acceleration of H{sup −} beams. Negative ions, up to 120 mA of current, are extracted from a radiofrequency driven plasma, by means of a gridded electrode (plasma grid, PG) featuring 9 apertures arranged in a 3x3 square lattice. The same aperture pattern is replicated in the following electrodes, allowing ion acceleration up to 60 keV. All electrodes are realized in copper, by electro-deposition technique, leaving empty slots in the metal to place magnets and to flow water for the grid cooling. The first set of electrodes was completed, installed in the source and tested. At the same time, an upgrade of the extraction system was carried out, in order to optimize the beam optics and to explore alternative electrostatic configurations. In particular, the accelerator will be modified by completely replacing the EG grid, exploiting the modularity of NIO1. The new electrode will feature other slots in between apertures, to place additional magnets. This allows testing different magnetic configurations, to optimize electron filtering and residual ion deflection. The present paper describes the theoretical activities driving the design of these new extractors, carried out with most updated numerical codes, and exploiting the synergy with the refined modeling of the 40 A ITER negative ion sources, under development at Consorzio RFX. Beam simulations are performed both with tracing codes (SLACCAD and OPERA) and with particle in cell codes (ACCPIC)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Preparatory Verification of a New Monte Carlo Code cosRMC based on VERA Core Physics Benchmark

To promote nuclear power plant system design and simulation software self-reliance development, DOE of China approved a R and D key project to develop an integrated engine for pressurized water reactor core and system design and analysis, referred as COSINE hereafter. The scope of this code package includes reactor physics, thermal hydraulics, safety analysis, severe accidents, probability safety assessment and etc. Many universities, research institutes and nuclear related corporations in China are involved in the code development, among which the State Nuclear Power Software Development Center (SNPSDC) plays the leading role. With the requirement of accurate three-dimensional modeling of the advanced, special or new complex core in reactor neutronics and shielding analysis, Monte Carlo method is becoming a more powerful tool for reactor analysis. Thus, a new Monte Carlo transport code named cosRMC (abbreviation for COSINE Reactor Monte Carlo) has been developed jointly by SNPSDC and Department of Engineering Physics, Tsinghua University, the main developer of RMC code. The cosRMC now has many functions as criticality calculation, fixed-source calculation, burnup calculation, and etc. Some techniques for geometry treatment, new burnup algorithm, source convergence acceleration, massive tally and parallel calculation have been implemented in cosRMC to improve the efficiency. In order to ensure its reliability, several types of V and V work have been and are being carried out for cosRMC. In this paper, some preparatory verification results based on VERA Core Physics Benchmark are presented. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

SU-E-T-81: Comparison of Microdosimetric Quantities Calculated Using the Track Structure Monte Carlo Algorithms Geant4-DNA and NOREC

Purpose: Microdosimetric quantities, such as the lineal energy, have been shown to correlate with the biological response to radiation and the relative biological effect of different radiation types. Track-structure Monte Carlo simulations are an important tool for investigating these responses and for developing mechanistic models to explain them. However, some of the cross-sectional data used in these algorithms has large uncertainties; thus, it is important to investigate how the implementation of the different codes affects the quantities of interest. Methods: Two of the most widely-used publicly available track-structure Monte Carlo codes, Geant4-DNA and NOREC, were used generate electron tracks for two particle sources. One source was a mono-energetic parallel beam of electrons with energies from 5 to 500-keV, and the lineal energy for each track was calculated in 1-mm-spheres arranged in planar arrays at multiple distances from the source. The second source was mono-energetic, uniformly-distributed, and isotropic source, and the lineal energy was scored in a single 30-mm-sphere for energies between 300-eV and 5-keV. Results: The dose-mean lineal energy for the parallel-beam simulations almost all agreed within 5%. For the uniformly-distributed source, at the lowest energies there was strong agreement between the algorithms, but the Geant4-DNA simulations showed slightly more high-energy events for more energetic electrons, but the dose-mean lineal energy agreed to within 4% for all energies. Conclusion: While there were slight differences in the results between the codes, these were consistent with previous studies of the stopping power and angular scattering distributions. Importantly, the computation time for Geant4-DNA was larger than for NOREC, largely due to approximations used in the NOREC for energies below 10-eV. This study shows that these approximation does not have a major impact on the microdosimetry on the energy and length scales investigated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Source scaling comparison and validation in Central Italy: data intensive direct S waves versus the sparse data coda envelope methodology

SUMMARY Robustness of source parameter estimates is a fundamental issue in understanding the relationships between small and large events; however, it is difficult to assess how much of the variability of the source parameters can be attributed to the physical source characteristics or to the uncertainties of the methods and data used to estimate the values. In this study, we apply the coda method by Mayeda et al. using the coda calibration tool (CCT), a freely available Java-based code (https://github.com/LLNL/coda-calibration-tool) to obtain a regional calibration for Central Italy for estimating stable source parameters. We demonstrate the power of the coda technique in this region and show that it provides the same robustness in source parameter estimation as a data-driven methodology [generalized inversion technique (GIT)], but with much fewer calibration events and stations. The Central Italy region is ideal for both GIT and coda approaches as it is characterized by high-quality data, including recent well-recorded seismic sequences such as L'Aquila (2009) and Amatrice–Norcia–Visso (2016–2017). This allows us to apply data-driven methods such as GIT and coda-based methods that require few, but high-quality data. The data set for GIT analysis includes ∼5000 earthquakes and more than 600 stations, while for coda analysis we used a small subset of 39 events spanning 3.5 < Mw < 6.33 and 14 well-distributed broad-band stations. For the common calibration events, as well as an additional 247 events (∼1.7 < Mw < ∼5.0) not used in either calibration, we find excellent agreement between GIT-derived and CCT-derived source spectra. This confirms the ability of the coda approach to obtain stable source parameters even with few calibration events and stations. Even reducing the coda calibration data set by 75 per cent, we found no appreciable degradation in performance. This validation of the coda calibration approach over a broad range of event size demonstrates that this procedure, once extended to other regions, represents a powerful tool for future routine applications to homogeneously evaluate robust source parameters on a national scale. Furthermore, the coda calibration procedure can homogenize the Mw estimates for small and large events without the necessity of introducing any conversion scale between narrow-band measures such as local magnitude (ML) and Mw, which has been shown to introduce significant bias.

Morasca, Paola (ORCID:0000000265254867)↗
DOE OSTI2020

SPISEA: A Python-based Simple Stellar Population Synthesis Code for Star Clusters

We present Stellar Population Interface for Stellar Evolution and Atmospheres (SPISEA), an open-source Python package that simulates simple stellar populations. The strength of SPISEA is its modular interface which offers the user control of 13 input properties including (but not limited to) the initial mass function, stellar multiplicity, extinction law, and the metallicity-dependent stellar evolution and atmosphere model grids used. The user also has control over the initial–final mass relation in order to produce compact stellar remnants (black holes, neutron stars, and white dwarfs). We demonstrate several outputs produced by the code, including color–magnitude diagrams, HR-diagrams, luminosity functions, and mass functions. SPISEA is object-oriented and extensible, and we welcome contributions from the community. The code and documentation are available on GitHub (https://github.com/astropy/SPISEA) and ReadtheDocs (https://spisea.readthedocs.io/en/latest/), respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Radiation shielding analysis of the concrete hutch of a long beamline of the Advanced Photon Source Upgrade project

The radiation shielding analysis of the external concrete enclosure of a feature long beamline of the Advanced Photon Source Upgrade project is described in this work. The beamline has 2 undulators that can operate in two different configurations and the source terms are estimated for the worst configuration using the SRW code in the OASYS package as well as the STAC8 code. The shielding analysis is carried out using the STAC8 code, and for the case of pink beam, with the FLUKA code also. Effective dose rates on contact with the outside surface of the concrete walls for pink, reflected full spectrum and monochromatic beams are calculated. For the monochromatic beams, 2 sets of discrete bandwidths (BW) from the Double Crystal Monochromator (DCM) and the full reflectance of the Double Multilayer Monochromator (DMM) are used. Further, due to the size of the direct beam and the energy domain of interest in the direct and scattered fields, the importance of the equivalent dose to the lens of the eye is discussed with respect to the shielding calculations of such facilities. The dose rates outside the concrete hutch in the lateral directions calculated using STAC8 is about 50% higher than the FLUKA results for the case of pink beam but may be considered as a reasonable agreement due to the differences in the calculational methodologies. Those dose rates in the lateral and forward directions are higher for the pink beam to allow full occupancy from the radiation protection point of view. For the full spectrum after 3 reflections, the effective dose rates are less than 0.5 Sv/h in the lateral direction but is higher in the forward region below about 22° with respect to the beam direction. For the DCM beams, except for the scenario with one mirror and 0.1%BW, all dose rates are less than 0.5 Sv/h outside the concrete hutch in the lateral directions. Here, the higher harmonics are what contributes to the dose outside the shield and the bandwidths assumed to derive the photon flux thus becomes important. The dose rates with the DMM are generally lower compared to the 0.1% BW DCM beam but are higher than the results obtained with XOP calculated bandwidths.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

The prototype of a detector for monitoring the cosmic radiation neutron flux on ground

This work presents a comparison between the results of experimental tests and Monte Carlo simulations of the efficiency of a detector prototype for on-ground monitoring the cosmic radiation neutron flux. The experimental tests were made using one conventional {sup 241}Am-Be neutron source in several incidence angles and the results were compared to that ones obtained with a Monte Carlo simulation made with MCNPX Code.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

ON THE FERMI LARGE AREA TELESCOPE SURPLUS OF DIFFUSE GALACTIC GAMMA-RAY EMISSION

Recent observations of diffuse Galactic γ-ray emission (DGE) by the Fermi Large Area Telescope (Fermi-LAT) have shown significant deviations, above a few GeV to about 100 GeV, from DGE models that use the GALPROP code for the propagation of cosmic ray (CR) particles outside their sources in the Galaxy and their interaction with the target distributions of the interstellar gas and radiation fields. The surplus of radiation observed is most pronounced in the inner Galaxy, where the concentration of CR sources is strongest. The present study investigates this 'Fermi-LAT Galactic Plane Surplus' by estimating the γ-ray emission from the sources themselves, which is disregarded in the above DGE models. It is shown that the expected hard spectrum of CRs, still confined in their sources (source cosmic rays, SCRs), can indeed explain this surplus. The method is based on earlier studies regarding the so-called EGRET GeV excess, which by now is generally interpreted as an instrumental effect. The contribution from SCRs is also predicted to increasingly exceed the DGE models above 100 GeV, up to γ-ray energies of about 10 TeV, where the corresponding surplus exceeds the hadronic part of the DGE by about one order of magnitude. Above such energies, the emission surplus should decrease again with energy due to the finite lifetime of the assumed supernova remnant sources. Observations of the DGE in the inner Galaxy at 15 TeV with the ground-based Milagro γ-ray detector and, at TeV energies, with the ARGO-YBJ detector are interpreted to provide confirmation of a significant SCR contribution to the DGE.

79 ASTRONOMY AND ASTROPHYSICS↗