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,009 records · Page 56

Pushing configuration-interaction to the limit: Towards massively parallel MCSCF calculations

A new large-scale parallel multiconfigurational self-consistent field (MCSCF) implementation in the open-source NWChem computational chemistry code is presented. The generalized active space approach is used to partition large configuration interaction (CI) vectors and generate a sufficient number of batches that can be distributed to the available cores. Massively parallel CI calculations with large active spaces can be performed. The new parallel MCSCF implementation is tested for the chromium trimer and for an active space of 20 electrons in 20 orbitals, which can now routinely be performed. Unprecedented CI calculations with an active space of 22 electrons in 22 orbitals for the pentacene systems were performed and a single CI iteration calculation with an active space of 24 electrons in 24 orbitals for the chromium tetramer was possible. In conclusion, the chromium tetramer corresponds to a CI expansion of one trillion Slater determinants (914 058 513 424) and is the largest conventional CI calculation attempted up to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implications of Voyager 1 observations beyond the heliopause for the local interstellar electron spectrum

Cosmic-ray observations made by the Voyager 1 spacecraft outside the dominant modulating influence of the heliosphere finally allow the comparison of computed galactic spectra with experimental data at lower energies. These computed spectra, based on galactic propagation models, can now be compared with observations at low energies by Voyager 1 and at high energies by the PAMELA space detector at Earth. This improves understanding of basic propagation effects and also provides solar modulation studies with reliable input spectra from 1 MeV to 100 GeV. We set out to reproduce the Voyager 1 electron observations in the energy range of 6-60 MeV, as well as the PAMELA electron spectrum above 10 GeV, using the GALPROP code. By varying the source spectrum and galactic diffusion parameters, specifically the rigidity dependence of spatial diffusion, we find local interstellar spectra that agree with both power-law spectra observed by Voyager 1 beyond the heliopause. The local interstellar spectrum between ∼1 MeV and 100 GeV indicates that it is the combination of two power laws, with E {sup –(1.45} {sup ±} {sup 0.15)} below ∼100 MeV and E {sup –(3.15} {sup ±} {sup 0.05)} above ∼100 MeV. A gradual turn in the spectral shape matching the power laws is found, between 2.0 ± 0.5) GeV and (100 ± 10) MeV. According to our simplified modeling, this transition is caused primarily by galactic propagation effects. We find that the intensity beyond the heliopause at 10 MeV is (350 ± 50) electrons m{sup –2} s{sup –1} sr{sup –1} MeV{sup –1}, decreasing to (50 ± 5) electrons m{sup –2} s{sup –1} sr{sup –1} MeV{sup –1} at 100 MeV.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of an open-source segmented blade design tool

As wind turbines continue to grow ever larger to reduce the cost of energy, their blades follow suit, with the largest commercial offshore blades extending past 100 m. Massive blades such as these raise key transportation and manufacturing challenges, especially for land-based turbines. Segmented blades are one solution and are garnering increased industry and research interest. In this work, a detailed mechanical joint model is integrated into the Wind-Plant Integrated System Design and Engineering Model (WISDEM ® ), which will facilitate future segmented blade research and optimization. WISDEM is used to design a wind turbine with 100-m segmented blades. This wind turbine design is compared to other machines with 100-m monolithic blades designed for rail-transportability. The designs are compared in terms of blade mass and cost, turbine capital cost, annual energy production, and levelized cost of energy, with monolithic designs being the lightest and most economical. However, this result may vary by wind plant location. A variety of segmentation joint types exist, and they will inevitably vary in parameters such as cost, spanwise location, and physical characteristics. This work examines the sensitivity of wind turbine design drivers and annual energy production to a variety of the aforementioned parameters, using the open-source wind turbine design codes OpenFAST and WISDEM, finding that joint mass, stiffness, and location can have significant effects on design drivers.

17 WIND ENERGY↗

ChemEnv : a fast and robust coordination environment identification tool

Coordination or local environments have been used to describe, analyze and understand crystal structures for more than a century. Here, a new tool called ChemEnv, which can identify coordination environments in a fast and robust manner, is presented. In contrast to previous tools, the assessment of the coordination environments is not biased by small distortions of the crystal structure. Its robust and fast implementation enables the analysis of large databases of structures. The code is available open source within the pymatgen package and the software can also be used through a web app available on http://crystaltoolkit.org through the Materials Project.

42 ENGINEERING↗

SU-E-I-66: Radiomics and Image Registration Updates for the Computational Environment for Radiotherapy Research (CERR)

Purpose: To present new tools in CERR for Radiomics, image registration and other software updates and additions. Methods: Radiomics: CERR supports generating 3-D texture metrics based on gray scale co-occurance. Two new ways to calculate texture features were added: (1) Local Texture Averaging: Local texture is calculated around a voxel within the userdefined bounding box. The final texture metrics are the average of local textures for all the voxels. This is useful to detect any local texture patterns within an image. (2) Image Smoothing: A convolution ball of user-defined radius is rolled over an image to smooth out artifacts. The texture metrics are then computed on the smooth image. Image Registration: (1) Support was added to import deformation vector fields as well as non-deformable transformation matrices generated by vendor software and stored in standard DICOM format. (2) Support was added to use image within masks while computing image deformations. CT to MR registration is supported. This registration uses morphological edge information within the images to guide the deformation process. In addition to these features, other noteworthy additions to CERR include (1) Irregularly shaped ROI: This is done by taking intersection between infinitely extended irregular polygons drawn on any of the two views. Such an ROI is more conformal and useful in avoiding any unwanted parts of images that cannot be avoided with the conventional cubic box. The ROI is useful to generate Radiomics metrics. (2) Ability to insert RTDOSE in DICOM format to existing CERR plans. (3) Ability to import multi-frame PET-CT and SPECT-CT while maintaining spatial registration between the two modalities. (4) Ability to compile CERR on Unix-like systems. Results: The new features and updates are available via https://www.github.com/adityaapte/cerr . Conclusion: Features added to CERR increase its utility in Radiomics, Image-Registration and Outcomes modeling.

07 ISOTOPE AND RADIATION SOURCES↗

SU-E-T-212: Influence of the Modulation Index On Daily Quality Assurance in Rapid Arc Treatments

Purpose: At our Institute some measured parameters for daily quality assurance (DQA) of dynamic arc therapy plans showed an unexpected behavior, therefore an investigation of the influence of the magnitude of modulation was conducted. Methods: In our clinical practice all DQAs of dynamic arc therapy plans are measured and analyzed prior to treatments using commercial software. For this study these plans were additionally exported to our in-house software written in MATLAB.The developed software extracted the leaf position, gantry angle, cumulative meterset weight of each control point (CP) and the total number of Monitor Units (MU) of each arc. Based on this information we calculated the leaf travel distance, irradiated segment area, number of MUs and dose rate for each CP. These data allowed us to calculate the modulation indexes (MI) of the plans, applying five different definitions of MI. The results were then correlated to the results of our DQA measurements.To validate the software, additional plans of known MIs were created and analyzed. For confirmation, the calculated parameters were compared to the segmented treatment table (STT) coming from the treatment planning system. Results: All calculated CP-parameters matched the STT by 99% or better. For linac 1, the comparison of the MI evaluation and the DQA results showed a slight tendency: 91.3% failed DQA plans had a MI lower than the average value. For this definition we consider that the lower the MI the higher the modulation. The results of the linac 2 present no significant relevance due to the low sample sizes for each DQA software. Conclusion: Available data and given definitions of the modulation index do not bring conclusive results; one cannot find a clear and distinct correlation with the failure of the DQA. The ongoing analysis with an increased sample size might lead to another conclusion.

07 ISOTOPE AND RADIATION SOURCES↗

SU-E-T-562: Scanned Percent Depth Dose Curve Discrepancy for Photon Beams with Physical Wedge in Place (Varian IX) Using Different Sensitive Volume Ion Chambers

Purpose: To investigate and report the discrepancy of scanned percent depth dose (PDD) for photon beams with physical wedge in place when using ion chambers with different sensitive volumes. Methods/Materials: PDD curves of open fields and physical wedged fields (15, 30, 45, and 60 degree wedge) were scanned for photon beams (6MV and 10MV, Varian iX) with field size of 5x5 and 10x10 cm using three common scanning chambers with different sensitive volumes - PTW30013 (0.6cm3), PTW23323 (0.1cm3) and Exradin A16 (0.007cm3). The scanning system software used was OmniPro version 6.2, and the scanning water tank was the Scanditronix Wellhoffer RFA 300.The PDD curves from the three chambers were compared. Results: Scanned PDD curves of the same energy beams for open fields were almost identical between three chambers, but the wedged fields showed non-trivial differences. The largest differences were observed between chamber PTW30013 and Exradin A16. The differences increased as physical wedge angle increased. The differences also increased with depth, and were more pronounced for 6MV beam. Similar patterns were shown for both 5x5 and 10x10 cm field sizes. For open fields, all PDD values agreed with each other within 1% at 10cm depth and within 1.62% at 20 cm depth. For wedged fields, the difference of PDD values between PTW30013 and A16 reached 4.09% at 10cm depth, and 5.97% at 20 cm depth for 6MV with 60 degree physical wedge. Conclusion: We observed a significant difference in scanned PDD curves of photon beams with physical wedge in place obtained when using different sensitive volume ion chambers. The PDD curves scanned with the smallest sensitive volume ion chamber showed significant difference from larger chamber results, beyond 10cm depth. We believe this to be caused by varying response to beam hardening by the wedges.

07 ISOTOPE AND RADIATION SOURCES↗

TH-AB-207A-07: Radiation Dose Simulation for a Newly Proposed Dynamic Bowtie Filters for CT Using Fast Monte Carlo Methods

Purpose: Dynamic bowtie filter is an innovative design capable of modulating the X-ray and balancing the flux in the detectors, and it introduces a new way of patient-specific CT scan optimizations. This study demonstrates the feasibility of performing fast Monte Carlo dose calculation for a type of dynamic bowtie filter for cone-beam CT (Liu et al. 2014 9(7) PloS one) using MIC coprocessors. Methods: The dynamic bowtie filter in question consists of a highly attenuating bowtie component (HB) and a weakly attenuating bowtie (WB). The HB is filled with CeCl3 solution and its surface is defined by a transcendental equation. The WB is an elliptical cylinder filled with air and immersed in the HB. As the scanner rotates, the orientation of WB remains the same with the static patient. In our Monte Carlo simulation, the HB was approximated by 576 boxes. The phantom was a voxelized elliptical cylinder composed of PMMA and surrounded by air (44cm×44cm×40cm, 1000×1000×1 voxels). The dose to the PMMA phantom was tallied with 0.15% statistical uncertainty under 100 kVp source. Two Monte Carlo codes ARCHER and MCNP-6.1 were compared. Both used double-precision. Compiler flags that may trade accuracy for speed were avoided. Results: The wall time of the simulation was 25.4 seconds by ARCHER on a 5110P MIC, 40 seconds on a X5650 CPU, and 523 seconds by the multithreaded MCNP on the same CPU. The high performance of ARCHER is attributed to the parameterized geometry and vectorization of the program hotspots. Conclusion: The dynamic bowtie filter modeled in this study is able to effectively reduce the dynamic range of the detected signals for the photon-counting detectors. With appropriate software optimization methods, the accelerator-based (MIC and GPU) Monte Carlo dose engines have shown good performance and can contribute to patient-specific CT scan optimizations.

60 APPLIED LIFE SCIENCES↗

Pececillo

Pececillo is a proxy-app for the open source Truchas metal processing code (LA-CC-15-097). It implements many of the physics models used in Truchas: free-surface, incompressible Navier-Stokes fluid dynamics (e.g., water waves); heat transport, material phase change, view factor thermal radiation; species advection-diffusion; quasi-static, elastic/plastic solid mechanics with contact; electomagnetics (Maxwell's equations). The models are simplified versions that retain the fundamental computational complexity of the Truchas models while omitting many non-essential features and modeling capabilities. The purpose is to expose Truchas algorithms in a greatly simplified context where computer science problems related to parallel performance on advanced architectures can be more easily investigated. While Pececillo is capable of performing simulations representative of typical Truchas metal casting, welding, and additive manufacturing simulations, it lacks many of the modeling capabilites needed for real applications.

Jibben, Zechariah↗

Flexible configuration-interaction shell-model many-body solver

BIGSTICK Is a flexible configuration-Interaction open-source shell-model code for the many-fermion problem In a shell model (occupation representation) framework. BIGSTICK can generate energy spectra, static and transition one-body densities, and expectation values of scalar operators. Using the built-in Lanczos algorithm one can compute transition probabflity distributions and decompose wave functions into components defined by group theory.

Johnson, CalvinW.↗

SCALE/DB-Public

This is a collection of public data produced by AMPX or curated from quality sources that the SCALE code system relies on for user problem input.

Brown, Jesse [Oak Ridge National Laboratory (ORNL)↗

Towards exascale for wind energy simulations

We examine large-eddy-simulation modeling approaches and computational performance of two open-source computational fluid dynamics codes for the simulation of atmospheric boundary layer flows that are of direct relevance to wind energy production. The first code, NekRS, is a high-order, unstructured-grid, spectral element code. The second code, AMR-Wind, is a second-order, block-structured, finite-volume code with adaptive mesh refinement capabilities. The objective of this study is to co-develop these codes in order to improve model fidelity and performance for each. These features will be critical for running ABL-based applications such as wind farm analysis on advanced computing architectures. To this end, we investigate the performance of NekRS and AMR-Wind on the Oak Ridge Leadership Facility supercomputers Summit, using 4 to 800 nodes (24 to 4,800 NVIDIA V100 GPUs), and Crusher, the testbed for the Frontier exascale system, using 18 to 384 Graphics Compute Dies on AMD MI250X GPUs. We compare strong- and weak-scaling capabilities, linear solver performance, and time to solution. We also identify leading inhibitors to parallel scaling.

17 WIND ENERGY↗

Final Report - PyVERAIn: An Open Source, Python-based CASL Preprocessing Platform for Reactor Geometry, Code Discretizations, and Parallel Domain Decompositions

The introduction of advanced modeling and simulation practices over the last several decades has proven to produce better products, more quickly, and at the reduced cost. This computational experimentation has enabled the design and development of a variety of innovative products by reducing the need to build expensive prototypes, and by allowing rapid, comprehensive design studies across a wide range of products. Despite these successes, there are significant barriers to widespread adoption of these tools. This project aims to address several of these barriers.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Bison Documentation Update, NEAMS Fuels Product Line Review Meeting

A brief overview of the status of the Bison documentation system is given. The Bison documentation system, which is built using the MOOSEDocs capability, has undergone significant development in the last two years to create a flexible, unified, and comprehensive documentation source for the Bison code. This documentation system, which is both hosted online (INL/MIS-18-50577) and can be built locally, integrates the older separate Bison user and theory manuals into a single format. The development of the documentation system, which is still ongoing, has three main objectives: improved step-by-step tutorial and workshop instructions for new users, implementation of detailed theory and input file information on a single page for each class, and increased accessibility of the documentation information. The progress and future work on these three objectives is reviewed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Developing an Automated Uncertainty Quantification Tool to Improve Watershed-Scale Predictions of Water and Nutrient Cycling

Managing the flow of water, nutrients, and contaminants in watersheds is vital to addressing pressing issues related to water scarcity, access to clean drinking water, energy production, resilience to natural and anthropogenic perturbations, and ecological restoration. Decisions about the management of watersheds critically depend on the accuracy with which the flow of water and chemicals through the watershed can be predicted by computer models. Prediction uncertainty can be reduced by matching the model to data, which are collected in the field at great expense. The contribution of watershed characterization data to reducing uncertainty of relevant model predictions can be evaluated in a so-called data-worth analysis, which provides transparent, quantitative metrics about a data set’s value for the support of relevant watershed management objectives. To achieve this goal, we developed a software package that implements the data-worth analysis approach for use with state-of-the-art watershed models. The purpose of the proposed data-worth analysis is to help decision-makers allocate resources for watershed characterization such that the uncertainty in model predictions can be significantly reduced, which leads to better, more effective management decisions. At the same time, watershed characterization costs can be reduced. The specific technical objectives of this SBIR/STTR Phase II project were to develop a framework and associated software toolsets that implement the uncertainty quantification and data-worth analysis approach for use with state-of-the-art watershed models. This goal was achieved by (A) developing a user-friendly, robust software package that is accessible to a wide audience, including watershed managers, policy-makers, and public stakeholders; (B) by demonstrating application of the prototype on several use cases that are representative of complex watershed management challenges spanning a range of scales and that consider different open-source, DOE-based codes and other modeling platforms; and (C) by gathering information about the needs and requirements from potential users to help guide future developments, ensuring that the final product will be commercially viable. The developed software consists of a graphical user interface that guides the user through a sequence of analysis steps, supported by toolsets that leverage state-of-the-art computational simulation-optimization capabilities. A prototype of the software runs on multiple platforms (PC, Mac, multi-processor Linux environment), is linked to diverse watershed simulators (e.g., ECOSYS, TOUGH2, TOUGHREACT, Amanzi-ATS), performs multiple analysis tasks (predictive simulations, sensitivity analysis, uncertainty analysis, automatic parameter estimation, and data-worth analysis, multicomponent geothermometry), and is readily extensible to include external simulators and analysis tools. The software is being commercialized and will be continually updated to address user needs.

58 GEOSCIENCES↗

Seascape Interface Control Document

This paper serves as the Interface Control Document (ICD) for the Seascape automated test harness developed at Sandia National Laboratories. The primary purposes of the Seascape system are: (1) provide a place for accruing large, curated, labeled data sets useful for developing and evaluating detection and classification algorithms (including, but not limited to, supervised machine learning applications) (2) provide an automated structure for specifying, running and generating reports on algorithm performance. Seascape uses GitLab, Nexus, Solr, and Banana, open source software, together with code written in the Python language, to automatically provision and configure computational nodes, queue up jobs to accomplish algorithms test runs against the stored data sets, gather the results and generate reports which are then stored in the Nexus artifact server.

97 MATHEMATICS AND COMPUTING↗

Use of Kerma for Subzone Dimension Specification in Monte Carlo Based Photon Energy Deposition Calculations

A technique using the photon kerma cross section for a material in combination with the number fraction from a photon energy spectrum has been developed to determine the estimated subzone dimension needed to provide an energy deposition profile in radiation transport calculations. The technique was verified using the ITS code for monoenergetic photon sources and a selection of photon spectra. A Python script was written to use the CEPXS cross-section file with a Rapture calculated transmission spectrum to provide the dimensional estimates in a rapid fashion. The script is available for SNL users through the corporate gitlab server.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗