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,333 records · Page 74

Parser for Sabin-to-Mahoney Transition Model of Quasispecies Replication

This code is a data parse for preparing output from the Qspp agent-based stochastic simulation model for plotting in Excel. This code is specific to a set of simulations that were run for the purpose of preparing data for a publication. It is necessary to make this code open-source in order to publish the model code (Qspp), which has already been released. There is a necessity of assuring that results from using Qspp for a publication

Ecale Zhou, Carol↗

Parser for Sabin-to-Mahoney Transition Model of Quasispecies Replication

This code is a data parse for preparing output from the Qspp agent-based stochastic simulation model for plotting in Excel. This code is specific to a set of simulations that were run for the purpose of preparing data for a publication. It is necessary to make this code open-source in order to publish the model code (Qspp), which has already been released. There is a necessity of assuring that results from using Qspp for a publication

Ecale Zhou, Carol↗

An efficient system for reliably transmitting image and video data over low bit rate noisy channels

This research project is intended to develop an efficient system for reliably transmitting image and video data over low bit rate noisy channels. The basic ideas behind the proposed approach are the following: employ statistical-based image modeling to facilitate pre- and post-processing and error detection, use spare redundancy that the source compression did not remove to add robustness, and implement coded modulation to improve bandwidth efficiency and noise rejection. Over the last six months, progress has been made on various aspects of the project. Through our studies of the integrated system, a list-based iterative Trellis decoder has been developed. The decoder accepts feedback from a post-processor which can detect channel errors in the reconstructed image. The error detection is based on the Huber Markov random field image model for the compressed image. The compression scheme used here is that of JPEG (Joint Photographic Experts Group). Experiments were performed and the results are quite encouraging. The principal ideas here are extendable to other compression techniques. In addition, research was also performed on unequal error protection channel coding, subband vector quantization as a means of source coding, and post processing for reducing coding artifacts. Our studies on unequal error protection (UEP) coding for image transmission focused on examining the properties of the UEP capabilities of convolutional codes. The investigation of subband vector quantization employed a wavelet transform with special emphasis on exploiting interband redundancy. The outcome of this investigation included the development of three algorithms for subband vector quantization. The reduction of transform coding artifacts was studied with the aid of a non-Gaussian Markov random field model. This results in improved image decompression. These studies are summarized and the technical papers included in the appendices.

Costello, Daniel J., Jr.↗

Development of the computer code for multi-calculations of source term of spent fuel - 15351

Many PWR spent fuels have been transferred between the neighboring nuclear power plants at Kori, Hanbit, Hanwul site, and from nuclear power plant to the KAERI(Korea Atomic Energy Research Institute) for the research purpose in Korea. In case of the PHWR type spent fuels, these have been transferred from the spent fuel pool to an on-site dry storage facility. Total radioactivity of spent fuels has to be calculated before transport in advance by the regulation requirement. In case of PWR, KN-12 transport cask can transfer 12 spent fuel assemblies at one time and in case of PHWR, Hi-STAR 63 transport cask can transfer 120 spent fuel assemblies at one time. Therefore, total radioactivity of all transported spent fuel assemblies that are transported generally are calculated by using the ORIGEN-ARP to meet this regulatory requirement. Basically, ORIGEN-ARP calculates the one fuel assembly at one time. SCALE code including this code provides the batch mode to perform the multi-calculations of the source term. In this study, in order to use the batch mode of ORIGEN-ARP, ARP-auto program was developed to generate input files and batch file of ORIGEN-ARP automatically and to make the summary file extracting the result from output files of ORIGEN-ARP. This code enhances the efficiency at the field. The working time to complete the calculation of source terms of 120 spent fuel assemblies take only 20 minutes. In order to confirm the credibility of this code, total radioactivity of 12 spent fuel assemblies of 17x17 Westinghouse type were calculated and compared to the previous result. The initial enrichments of these spent fuels are 3.46 ∼ 3.50 w/o and burnups are 29 ∼ 43 GWD/MTU. As the results, the relative errors are about 0.9 ∼ 3.0%. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Model for Simulating a Spiral Software-Development Process

A discrete-event simulation model, and a computer program that implements the model, have been developed as means of analyzing a spiral software-development process. This model can be tailored to specific development environments for use by software project managers in making quantitative cases for deciding among different software-development processes, courses of action, and cost estimates. A spiral process can be contrasted with a waterfall process, which is a traditional process that consists of a sequence of activities that include analysis of requirements, design, coding, testing, and support. A spiral process is an iterative process that can be regarded as a repeating modified waterfall process. Each iteration includes assessment of risk, analysis of requirements, design, coding, testing, delivery, and evaluation. A key difference between a spiral and a waterfall process is that a spiral process can accommodate changes in requirements at each iteration, whereas in a waterfall process, requirements are considered to be fixed from the beginning and, therefore, a waterfall process is not flexible enough for some projects, especially those in which requirements are not known at the beginning or may change during development. For a given project, a spiral process may cost more and take more time than does a waterfall process, but may better satisfy a customer's expectations and needs. Models for simulating various waterfall processes have been developed previously, but until now, there have been no models for simulating spiral processes. The present spiral-process-simulating model and the software that implements it were developed by extending a discrete-event simulation process model of the IEEE 12207 Software Development Process, which was built using commercially available software known as the Process Analysis Tradeoff Tool (PATT). Typical inputs to PATT models include industry-average values of product size (expressed as number of lines of code), productivity (number of lines of code per hour), and number of defects per source line of code. The user provides the number of resources, the overall percent of effort that should be allocated to each process step, and the number of desired staff members for each step. The output of PATT includes the size of the product, a measure of effort, a measure of rework effort, the duration of the entire process, and the numbers of injected, detected, and corrected defects as well as a number of other interesting features. In the development of the present model, steps were added to the IEEE 12207 waterfall process, and this model and its implementing software were made to run repeatedly through the sequence of steps, each repetition representing an iteration in a spiral process. Because the IEEE 12207 model is founded on a waterfall paradigm, it enables direct comparison of spiral and waterfall processes. The model can be used throughout a software-development project to analyze the project as more information becomes available. For instance, data from early iterations can be used as inputs to the model, and the model can be used to estimate the time and cost of carrying the project to completion.

Mizell, Carolyn↗

Data and Scripts associated with “Lambda-PFLOTRAN: Workflow for Incorporating Organic Matter Chemistry Informed by Ultra High Resolution Mass Spectrometry into Biogeochemical Modeling.”

This data package is associated with the publication “Lambda-PFLOTRAN: Workflow for Incorporating Organic Matter Chemistry Informed by Ultra High Resolution Mass Spectrometry into Biogeochemical Modeling” submitted to Geoscientific Model Development (Muller et al., 2024). In this manuscript, organic matter chemistry and thermodynamics are directly connected to reactive transport simulators through the newly developed Lambda-PFLOTRAN (Parallel Reactive Flow and Transport model) workflow tool that succinctly incorporates organic matter chemistry data generated from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) into reaction networks to simulate aerobic respiration of the organic matter and the resulting biogeochemistry. Lambda-PFLOTRAN is a python-based workflow, executed through a Jupyter Notebook interface, that digests raw FTICR-MS data, develops a representative reaction network based on substrate-explicit thermodynamic modeling (also termed lambda modeling due to its key thermodynamic parameter λ used therein), and completes a biogeochemical simulation with the open source, reactive flow, and transport code PFLOTRAN. This data package contains Jupyter Notebook based workflows for two test cases for running biogeochemical simulations of organic matter oxidation identified by FTICR-MS. It contains four primary folders (workflow, data, src, and analysis), a file-level metadata file (Muller_2024_Lambda_PFLOTRAN_Manuscript_Data_Package_flmd.csv) that lists all the files contained in this data package with a short description of each, and a data dictionary (Muller_2024_Lambda_PFLOTRAN_Manuscript_Data_Package_dd.csv) file that describes the tabular column headers. The ‘workflow’ folder contains the Jupyter Notebook based workflows for running the lambda analysis, PFLOTRAN simulation, sensitivity analysis and parameter estimation. The ‘data’ folder contains the FTICR-MS data, initial conditions, and incubation data for test cases 1 and 2 in folders titled ‘WHONDRS’ and ‘Colloids’, respectively. The data folder also has a ‘Database’ folder containing a reaction network for bulk organic matter (assumed to be CH2O) and a general database for PFLOTRAN (hanford_rxn_network). The CH2O reaction network defines bulk organic matter oxidation. Biogeochemical simulations are completed for both the lambda binned organic matter and bulk organic matter reaction networks. The ‘hanford_rxn_network’ database includes information required for PFLTORAN simulations including ion size, molar mass, and charge of the aqueous species, gases, and minerals phases. The ‘src’ folder contains python source codes for performing lambda analysis, PFLOTRAN simulation, sensitivity analysis and parameter estimation. The ‘analysis’ folder contains outputs from the test cases 1 and 2 including lambda analysis, PFLOTRAN runs and the calibration results.

54 ENVIRONMENTAL SCIENCES↗

SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD

We present a new capability of the ice sheet model SICOPOLIS that enables flexible adjoint code generation via source transformation using the open-source algorithmic differentiation (AD) tool OpenAD.The adjoint code enables efficient calculation of the sensitivities of a scalar-valued objective function or quantity of interest (QoI) to a range of important, often spatially varying and uncertain model input variables, including initial and boundary conditions, as well as model parameters.Compared to earlier work on the adjoint code generation of SICOPOLIS, our work makes several important advances:(i) it is embedded within the up-to-date trunk of the SICOPOLIS repository – accounting for 1.5 decades of code development and improvements – and is readily available to the wider community;(ii) the AD tool used, OpenAD, is an open-source tool;(iii) the adjoint code developed is applicable to both Greenland and Antarctica, including grounded ice as well as floating ice shelves, with an extended choice of thermodynamical representations.A number of code refactorization steps were required. They are discussed in detail in an Appendix as they hold lessons for the application of AD to legacy codes at large.As an example application, we examine the sensitivity of the total Antarctic Ice Sheet volume to changes in initial ice thickness, austral summer precipitation, and basal and surface temperatures across the ice sheet.Simulations of Antarctica with floating ice shelves show that over 100 years of simulation the sensitivity of total ice sheet volume to the initial ice thickness and precipitation is almost uniformly positive, while the sensitivities to surface and basal temperature are almost uniformly negative. Sensitivity to austral summer precipitation is largest on floating ice shelves from Queen Maud to Queen Mary Land. The largest sensitivity to initial ice thickness is at outlet glaciers around Antarctica. Comparison between total ice sheet volume sensitivities to surface and basal temperature shows that surface temperature sensitivities are higher broadly across the floating ice shelves, while basal temperature sensitivities are highest at the grounding lines of floating ice shelves and outlet glaciers. A uniformly perturbed region of East Antarctica reveals that, among the four control variables tested here, total ice sheet volume is the most sensitive to variations in austral summer precipitation as formulated in SICOPOLIS.Comparison between adjoint- and finite-difference-derived sensitivities shows good agreement, lending confidence that the AD tool is producing correct adjoint code.The new modeling infrastructure is freely available at http://www.sicopolis.net (last access: 2 April 2020) under the development trunk.

58 GEOSCIENCES↗

Titanium alloy response sensitivity to variations in spectral reconstructions of National Ignition Facility xenon line-emission x-ray sources

Thermomechanical shock experiments on the National Ignition Facility (NIF) aim to study high strain rate dynamic material response. In such experiments, the NIF laser is used to generate high fluence x-ray emission sources, which irradiate material samples of interest. Under sufficiently high x-ray energy deposition, thermomechanical impulses are generated in the materials. While it is known that the characteristics of x-ray generated impulses vary as a function of incident x-ray spectra, it remains unclear how spectral assumptions and uncertainties in NIF spectral reconstructions affect our interpretation of impulsive loading. Here, in this paper, we simulate the response of a standard titanium alloy baseline sample to synthetic analytically derived and measured NIF xenon line-emission x-ray sources with a radiation hydrodynamics code. We vary the source spectral characteristics based on different source reconstruction techniques to understand the resulting variation in baseline sample response and compare the simulated response with experimental results. We find that the response is highly sensitive to assumptions made about the spectral contents and that knowledge of spectral uncertainties bounds our understanding of the resulting material response. The results of this effort help to extend our ability to use baseline material samples to extract quantitative properties from x-ray experiments on the NIF.

Alloys↗

SU-F-T-28: Evaluation of BEBIG HDR Co-60 After-Loading System for Skin Cancer Treatment Using Conical Surface Applicator

Purpose: To evaluate the possibility of utilizing the BEBIG HDR 60Co remote after-loading system for malignant skin surface treatment using Monte Carlo (MC) simulation technique. Methods: First TG-43 parameters of BEBIG-Co-60 and Nucletron Ir-192-mHDR-V2 brachytherapy sources were simulated using MCNP6 code to benchmark the sources against the literature. Second a conical tungsten-alloy with 3-cm diameter of Planning-Target-Volume (PTV) at surface for use with a single stepping HDR source is designed. The HDR source is modeled parallel to treatment plane at the center of the conical applicator with a source surface distance (SSD) of 1.5-cm and a removable plastic end-cap with a 1-mm thickness. Third, MC calculated dose distributions from HDR Co-60 for conical surface applicator were compared with the simulated data using HDR Ir-192 source. The initial calculations were made with the same conical surface applicator (standard-applicator) dimensions as the ones used with the Ir-192 system. Fourth, the applicator wall-thickness for the Co-60 system was increased (doubled) to diminish leakage dose to levels received when using the Ir-192 system. With this geometry, percentage depth dose (PDD), and relative 2D-dose profiles in transverse/coronal planes were normalized at 3-mm prescription-depth evaluated along the central axis. Results: PDD for Ir-192 and Co-60 were similar with standard and thick-walled applicator. 2D-relative dose distribution of Co-60, inside the standard-conical-applicator, generated higher penumbra (7.6%). For thick-walled applicator, it created smaller penumbra (<4%) compared to Ir-192 source in the standard-conicalapplicator. Dose leakage outside of thick-walled applicator with Co-60 source was approximately equal (≤3%) with standard applicator using Ir-192 source. Conclusion: Skin cancer treatment with equal quality can be performed with Co-60 source and thick-walled conical applicators instead of Ir-192 with standard applicators. These conical surface applicator must be used with a protective plastic end-cap to eliminate electron contamination and over-dosage of the skin.

60 APPLIED LIFE SCIENCES↗

shiftNMFk 1.1: Robust Nonnegative matrix factorization with kmeans clustering and signal shift, for allocation of unknown physical sources, toy version for open sourcing with publications

This code is a toy (short) version of CODE-2016-83. From a general perspective, the code represents an unsupervised adaptive machine learning algorithm that allows efficient and high performance de-mixing and feature extraction of a multitude of non-negative signals mixed and recorded by a network of uncorrelated sensor arrays. The code identifies the number of the mixed original signals and their locations. Further, the code also allows deciphering of signals that have been delayed in regards to the mixing process in each sensor. This code is high customizable and it can be efficiently used for a fast macro-analyses of data. The code is applicable to a plethora of distinct problems: chemical decomposition, pressure transient decomposition, unknown sources/signal allocation, EM signal decomposition. An additional procedure for allocation of the unknown sources is incorporated in the code.

Alexandrov, BoianS.↗

shiftNMFk 1.1: Robust Nonnegative matrix factorization with kmeans clustering and signal shift, for allocation of unknown physical sources, toy version for open sourcing with publications

This code is a toy (short) version of CODE-2016-83. From a general perspective, the code represents an unsupervised adaptive machine learning algorithm that allows efficient and high performance de-mixing and feature extraction of a multitude of non-negative signals mixed and recorded by a network of uncorrelated sensor arrays. The code identifies the number of the mixed original signals and their locations. Further, the code also allows deciphering of signals that have been delayed in regards to the mixing process in each sensor. This code is high customizable and it can be efficiently used for a fast macro-analyses of data. The code is applicable to a plethora of distinct problems: chemical decomposition, pressure transient decomposition, unknown sources/signal allocation, EM signal decomposition. An additional procedure for allocation of the unknown sources is incorporated in the code.

Alexandrov, Boian S.↗

An observation of the Galactic center hard X-ray source, 1E 1740.7-2942, with the Caltech coded-aperture telescope

The Galactic center region hard X-ray source IE 1740.7-2942 has been observed with the Caltech Gamma-Ray Imaging Payload (GRIP) from Alice Springs, Australia, on 1988 April 12 and on 1989 April 3 and 4. We report here results from the 1989 measurements based on 14 hr of observation of the Galactic center region. The observations showed IE 1740.7-2942 to be in its normal state, having a spectrum between 35 and 200 keV characterized by a power law with an exponent of -2.2 +/- 0.3 and flux at 100 keV of (7.0 +/- 0.7) x 10 exp -5 sq cm s keV. No flux was detected above 200 keV. A search for time variability in the spectrum of IE 1740.7-2942 on one hour time scales showed no evidence for variability.

Heindl, William A.↗

PFLOTRAN-RepoTREND Source Term Comparison Summary.

Code inter-comparison studies are useful exercises to verify and benchmark independently developed software to ensure proper function, especially when the software is used to model high-consequence systems which cannot be physically tested in a fully representative environment. This summary describes the results of the first portion of the code inter-comparison between PFLOTRAN and RepoTREND, which compares the radionuclide source term used in a typical performance assessment.

97 MATHEMATICS AND COMPUTING↗

Development of an Integrated Mechanistic Source Term Assessment Capability for Lead- and Sodium-Cooled Fast Reactors

The scope of this project entails development of mechanistic source term assessment capabilities to inform and guide plant design activities and to support license applications. The main objective of this project is to couple NL’s SAS4A safety analysis code with the FAI’s FATE facility modeling code for mechanistic source term assessments in liquid-metal-cooled fast reactors. The end product is a coupled SAS4A-FATE code system for integrated mechanistic source term assessments and Level-2 PRA as the precursor for off-site dispersion analyses (Level-3 PRA). The SAS4A code is used to model system transients with fuel failures and its results are linked to the FATE code to predict radionuclide transport through the reactor vessel cover gas and containment system. The integrated code package offers a state-of-the-art radionuclide tracking capability for a spectrum of accident scenarios consistent with the plant dynamic response including the reactivity feedback to support Probabilistic Risk Assessments (PRA).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Image coding by adaptive block quantization.

A new source encoder called the adaptive block quantizer is proposed for coding data sources that emit a sequence of correlated real numbers with known first- and second-order statistics. Blocks of source output symbols are first classified and then block quantized in a manner that depends on their classification. The system is optimized relative to both the mean square error and the subjective quality of the reconstructed data for a certain class of pictorial data, and the resulting system performance demonstrated. Some interesting relationships between mean square error and subjective picture quality are presented.

Tasto, M.↗

Nonlinear Optimization in "Synchrotron Radiation Workshop" Code

"Synchrotron Radiation Workshop" (SRW) is an open source physical optics computer code for calculation of detailed characteristics of Synchrotron Radiation (SR) generated by relativistic electrons in magnetic fields of arbitrary configuration. SRW also supports high-accuracy simulations of fully- and partially-coherent radiation propagation through X-ray optical beamlines, facilitated by so-called “Virtual Beamline” module. In this present tech note, we discuss the nonlinear optimization features of SRW with emphasis on the optimization library in SRW. This optimization library not only can be used to optimize the simulation beamline in SRW, it also can be integrated with real beamline operations and optimize the parameters based on real measurements. We also show some demos of the optimization library recently developed in SRW.

36 MATERIALS SCIENCE↗

Recent Advances in the Use of TRIPOLI-4 Monte Carlo Code for Radiation Protection and Shielding Calculations - Paper 47

TRIPOLI-4 is a general purpose radiation transport code. It uses the continuous-energy Monte Carlo method to simulate neutron, photon, electron and positron transport in three-dimensional geometries. The main application areas of TRIPOLI-4 code include radiation protection and shielding, fission reactor physics, nuclear criticality safety, fusion reactor design, and nuclear instrumentation. For the radiation protection and shielding calculations with TRIPOLI-4 code, miscellaneous radiation sources options, various variance reduction options, and many useful tally options are available. Perturbation calculations on material densities and nuclide concentrations are routinely used in radiation protection applications. Parallel calculations with TRIPOLI-4 are also easy to perform with a PC under the linux/unix OS. To help the users to get the feedback from their test runs during the input data preparation and adjustment and to analyze the intermediate and final results from their production runs, several user-interface tools have been introduced. By using some application cases in the fields of radiation shielding, radiation dosimetry, and gamma spectrometry, this paper summarizes and demonstrates some essential improvements and possibilities in the use of TRIPOLI-4 code on fixed source calculations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗