Search NASA⌕ Search

SEARCH · Search NASA

Results for “default parameters”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-E-T-587: Monte Carlo Versus Ray-Tracing for Treatment Planning Involving CNS Tumors On the MultiPlan System for CyberKnife Radiosurgery

Purpose: MultiPlan, the treatment planning system for the CyberKnife Robotic Radiosurgery system offers two approaches to dose computation, namely Ray-Tracing (RT), the default technique and Monte Carlo (MC), an option. RT is deterministic, however it accounts for primary heterogeneity only. MC on the other hand has an uncertainty associated with the calculation results. The advantage is that in addition, it accounts for heterogeneity effects on the scattered dose. Not all sites will benefit from MC. The goal of this work was to focus on central nervous system (CNS) tumors and compare dosimetrically, treatment plans computed with RT versus MC. Methods: Treatment plans were computed using both RT and MC for sites covering (a) the brain (b) C-spine (c) upper T-spine (d) lower T-spine (e) L-spine and (f) sacrum. RT was first used to compute clinically valid treatment plans. Then the same treatment parameters, monitor units, beam weights, etc., were used in the MC algorithm to compute the dose distribution. The plans were then compared for tumor coverage to illustrate the difference if any. All MC calculations were performed at a 1% uncertainty. Results: Using the RT technique, the tumor coverage for the brain, C-spine (C3–C7), upper T-spine (T4–T6), lower T-spine (T10), Lspine (L2) and sacrum were 96.8%, 93.1%, 97.2%, 87.3%, 91.1%, and 95.3%. The corresponding tumor coverage based on the MC approach was 98.2%, 95.3%, 87.55%, 88.2%, 92.5%, and 95.3%. It should be noted that the acceptable planning target coverage for our clinical practice is >95%. The coverage can be compromised for spine tumors to spare normal tissues such as the spinal cord. Conclusion: For treatment planning involving the CNS, RT and MC appear to be similar for most sites but for the T-spine area where most of the beams traverse lung tissue. In this case, MC is highly recommended.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2010

Energy dependence of elliptic flow from heavy-ion collision models

We have compared the experimental data on charged-particle elliptic flow parameter (v{sub 2}) in Au + Au collisions at midrapidity for {radical}(s{sub NN})=9.2, 19.6, 62.4, and 200 GeV with results from various models in heavy-ion collisions such as the ultrarelativistic quantum molecular dynamics (UrQMD) model, a multiphase transport model (AMPT), and heavy-ion jet interaction generator (HIJING). We observe that the average v{sub 2} (<v{sub 2}>) from the transport model UrQMD agrees well with the measurements at {radical}(s{sub NN})=9.2 GeV but increasingly falls short of the experimental <v{sub 2}> values as the beam energy increases. The difference in <v{sub 2}> is of the order of 60% at {radical}(s{sub NN})=200 GeV. The <v{sub 2}> results from HIJING are consistent with zero, while those from AMPT with default settings, a model based on HIJING with additional initial- and final-state rescattering effects included, give a <v{sub 2}> value of about 4% for all the beam energies studied. This is in contrast to an increase in <v{sub 2}> with beam energy for the experimental data. A different version of the AMPT model, which includes partonic effects and quark coalescence as a mechanism of hadronization, gives higher values of <v{sub 2}> among the models studied and is in agreement with the measured <v{sub 2}> values at {radical}(s{sub NN})=200 GeV. These studies show that the experimental <v{sub 2}> has substantial contribution from partonic interactions at {radical}(s{sub NN})=200 GeV, whose magnitude reduces with decrease in beam energy. We also compare the available data on the transverse momentum and pseudorapidity dependence of v{sub 2} to those from these models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

ELM model simulations of Plum Island Ecosystems LTER low marsh site 2018-2020

Model simulations using the Department of Energy's Energy Exascale Earth System Model (E3SM) land model (ELM) with improved capabilities to represent vegetation response to salinity and inundation. The simulations were conducted for a tidal salt marsh at Plum Island Ecosystems Long Term Ecological Research (LTER) site near Rowley, Massachusetts, USA; the site is a low marsh dominated by Spartina alterniflora. The model was forced with site-specific meteorology, salinity and tidal cycles from 2018-2020. Four sets of model simulations are included and described below:1. Parameterization of the salinity response function. These simulations tested different combinations of values for optimal salinity and salinity tolerance.2. Model evaluation. This comparison conducted simulations using the default model, the salinity function only, the submergence function only, and both the salinity and submergence functions. 3. Salinity scenarios. These simulations used the 2018 salinity input data varied by -5 to +10 ppt salinity.4. Water level scenarios. These simulations used the tide height varied by -10 to +50 cm. These simulations were used to demonstrate how the salinity and submergence functions better represent carbon uptake by tidal salt marshes.The data package includes netCDF files used as forcing files for tide height and salinity, one for each year 2018-2020 at observed salinity concentrations, and an additional three forcing files in which salinity concentrations were varied 5 ppt lower, 5 ppt higher, and 10 ppt higher than the measured 2018 time series. Also included are python scripts for creating forcing files, plain text parameter and command files for running simulations, model outputs in netCDF format, and python scripts for visualizing outputs. Code for the modified E3SM model is archived in Sulman et al 2023 at doi:10.15485/1991625. More detail about files is provided in the README.md file.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Performance Assessment of a New Gamma Spectroscopy Toolkit: Peak Doctor and the Spectral Nondestructive Assay Platform (SNAP) - 17397

High resolution gamma spectroscopy (HRGS) is a widely used technique for non-destructive characterization of nuclear materials, and assay of radioactive waste. The ability to accurately evaluate raw spectral data is an important step in producing high-quality results. A spectral analysis application known as Peak Doctor has been developed under a Cooperative Research and Development Agreement (CRADA) between LANL and Pajarito Scientific Corporation (PSC). The software has undergone several years of testing and validation by the Radioactive Waste Characterization Team of the LANL Waste Management Services Group (WM-SVS). PeakDoctor is designed to work with the Spectral Nondestructive Assay Platform (SNAP{sup TM}) to provide a complete HRGS analysis engine. SNAP{sup TM} performs model specific corrections for attenuation and geometry. These applications are the major components of a new toolkit designed primarily for portable far-field non-destructive assay. Peak Doctor performs a very accurate peak fit analysis with the goal of improving upon the capabilities of existing commercially available routines. It is designed to be user-configurable and can be customized to any given type of detector and any given type of spectrum. Peak Doctor was developed to resolve many of the inherent complexities in gamma ray spectra including overlapping peaks, complex multiplet regions, large step peaks, backscatter regions and Compton edges. Unlike most spectral analysis platforms, Peak Doctor does not rely on the energy calibration and detector resolution parameters that are saved by the operator in the spectrum file. Instead, key parameters are pre-defined for individual detectors. The software fits a continuum using a series of splines between a series of junctions or 'knots'. Fitting techniques are used to join the splines together to create a near-continuous multi-element curve. Optimizing the location of the knots is one of Peak Doctor's most powerful features. The knots are initially configured using a pre-defined list and are removed or added automatically during the fitting process. Knots are also created when new features are identified such as Step Peaks and Compton Edges. The user can also manually specify knots in order to improve the continuum fit. The initial fit is improved on by removing any outlier regions (which are converted to peaks) and the continuum fitting process is repeated. At this stage, an experienced user can work interactively with the fit by manually converting peak regions to continuum or vice versa. The final step is to fit a series of peaks above the continuum based on detector parameters. Single and multiplet peaks are then processed by either the tailing (piecewise) or non-tailing (Gaussian) Levenberg Marquardt non-linear peak fitting routine depending on the user defined selection. Only a rough knowledge of the detector characteristics (energy calibration and resolution) is required for the fitting process to proceed with good accuracy. Any peaks that are missed can be manually added by the user and will also be fitted. Discontinuities in the peak / continuum fit can be rapidly identified by reviewing the 'residuals scale' plot. For waste streams that are well characterized, it is often the case that routinely encountered spectra will have similar peaks and continuum features. For these cases a handy 'One Step' method is available that proceeds in an automated fashion without user interaction. Peak Doctor has been configured with a set of default spectrum types including plutonium / transuranic waste (including Weapons Grade Pu, Heat Source Pu, Am-241, Am-243 and Np-237), depleted or enriched uranium, nuclear power plant / intermediate level waste (with fission and activation products such as Cs-137 and Co-60) as well as Naturally Occurring Radioactive Materials (NORM). The user can easily add new spectrum types for specific waste streams or other types of radiological / special nuclear materials. The output from PeakDoctor (including net peak area, centroid energy and a list of minimum detectable activity levels at pre-defined regions) is used by SNAP{sup TM} to complete the data analysis process. Isotope activities are calculated using configurable template-driven geometry and attenuation correction methods. The software accounts for matrix and container wall effects, collimator / angular corrections, detector efficiency (based on calibration curves) and stand-off distance. SNAP{sup TM} uses published algorithms incorporating a series of decision points during its operation to allow the analyst to provide logical input to guide the software during the analysis of each unique spectrum. The Peak Doctor / SNAP{sup TM} toolkit has been extensively tested and validated under a variety of conditions using 'benchmark' reference spectra. Test scenarios reported include HPGe-based measurements of well characterized uranium and plutonium standards and radioactive sources with a variety of different geometries and matrix materials including surrogate drums and boxes. (authors)

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

Automation and optimization of stopping and range of ions in matter simulation runtime

Prior to every ion implantation experiment a simulation of the ion range and other relevant parameters is performed using Monte-Carlo based codes. Although increasing computing power has improved the speed of these calculations, the demands on Monte-Carlo codes are also increasing, requiring evaluation of the optimal number of simulations while ensuring accuracy within threshold bounds. We evaluate the “Stopping and Range of Ions in Matter” (SRIM) code due to its widespread usage. We show how dividing simulations into multiple parallel simulations with different random seeds can lead to calculation speedup and find lower bounds for the required number of ion traces simulated based on an exemplar system of a Ga focused ion beam and a high energy C beam as used in high linear energy transfer testing. Here our results indicate simulations can yield results within the underlying data accuracy of SRIM at 10X and 100X shorter simulation time than the SRIM default values.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2007

Full Life-Cycle Defect Management Assessment: Initial Inspection Data Collection Results and Research Questions for Further Study

It is often the case in software projects that when schedule and budget resources are limited, the Verification and Validation (V&V) activities suffer. Fewer V&V activities can be afforded and moreover, short-term challenges can result in V&V activities being scaled back or dropped altogether. As a result, too often the default solution is to save activities for improving software quality until too late in the life-cycle, relying on late-term code inspections followed by thorough testing activities to reduce defect counts to acceptable levels. As many project managers realize, however, this is a resource-intensive way of achieving the required quality for software. The Full Life-cycle Defect Management Assessment Initiative, funded by NASA s Office of Safety and Mission Assurance under the Software Assurance Research Program, aims to address these problems by: Improving the effectiveness of early life-cycle V&V activities to make their benefits more attractive to team leads. Specifically, we focus on software inspection, a proven method that can be applied to any software work product, long before executable code has been developed; Better communicating this effectiveness to software development teams, along with suggestions for parameters to improve in the future to increase effectiveness; Analyzing the impact of early life-cycle V&V on the effectiveness and cost required for late life-cycle V&V activities, such as testing, in order to make the tradeoffs more apparent. This white paper reports on an initial milestone in this work, the development of a preliminary model of inspection effectiveness across multiple NASA Centers. This model contributes toward reaching our project goals by: Allowing an examination of inspection parameters, across different types of projects and different work products, for an analysis of factors that impact defect detection effectiveness. Allowing a comparison of this NASA-specific model to existing recommendations in the literature regarding how to plan effective inspections. Forming a baseline model which can be extended to incorporate factors describing: the numbers and types of defects that are missed by inspections; how such defects flow downstream through software development phases; how effectively they can be caught by testing activities in the late stages of development. The model has been implemented in a prototype web-enabled decision-support tool which allows developers to enter their inspection data and receive feedback based on a comparison against the model. The tool also allows users to access reusable materials (such as checklists) from projects included in the baseline. Both the tool itself and the model underlying it will continue to be extended throughout the remainder of this initiative. As results of analyzing inspection effectiveness for defect containment are determined, they can be shared via the tool and also via updates to existing training courses on metrics and software inspections. Moreover, the tool will help satisfy key CMMI requirements for the NASA Centers, as it will enable NASA to take a global view across peer review results for various types of projects to identify systemic problems. This analysis can result in continuous improvements to the approach to verification.

Shull, Forrest↗
DOE OSTI2024

MITgcm-AD v2: Open source tangent linear and adjoint modeling framework for the oceans and atmosphere enabled by the Automatic Differentiation tool Tapenade

The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.

Adjoints↗
DOE OSTI2014

SU-E-J-271: A Feasibility Study to Utilize XR-QA2 Radiochromic Films as An Alternative Tool for Measuring Beam Parameters of Kilovoltage CBCT Systems

Purpose: The beam parameters of the x-ray beams used in image guidance are difficult to measure in water due to the low exposure. This study is to investigate if XR-QA2 radiochromic films can be used in obtaining beam parameters, such as beam profiles and outputs for kV-CBCT systems. Methods: The kV-CBCT beams were from a Varian OBI system. The films used were XR-QA2-GAFCHROMICTM films placed in the middle of two Plastic -Water slab phantoms at isocenter. The scanner used to analyze the films was an Epson Expression 1680 flat-bed scanner. Pre and 24 h-post-irradiated scans of films were obtained by using the reflective mode. The corresponding pixel-by-pixel dose in the films was obtained using a dose calibration equation as a function of net reflectance from the literature. Beam profiles of different scan protocols with and without bow-tie filters were measured. Monte Carlo simulations were used to compare the accuracy of the measured dose profiles. Results: The beam profiles measured from XR-QA2 films for different clinical default kV-CBCT scan protocols were benchmarked against Monte Carlo simulated profiles which have been validated by measurements performed in water phantom. We confirmed that a minimum post-irradiated time of 24 hours is necessary to pass to scan the films to account for light scatter stability. We found that the uncertainties of the results depend on the scan orientation of the film when analyzed in different directions. Conclusion: The XR-QA2 radiochromic films are sensitive to measure the x-ray beam profiles with irradiated dose ranging 1–5 cGy. The film measurements can be performed with an easy setup compared to ion chamber measurements in a water phantom. These films provide an alternative method to obtain the beam profiles for kV-CBCT beams.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Description of Sensor Assignment Optimization Method as Deployed on a Multi-Node Cluster

Data analytic methods are being developed to address the problem of how to assign a sensor set in a nuclear facility such that a requisite level of process monitoring capability is realized and that the sensor set is sufficiently rich to determine the status of the individual sensors with respect to need for calibration. There is an awareness in the nuclear industry that data analytics combined with rich sensor sets represent a means to improve operations and reduce costs. In the industry the calibration problem has been previously approached as an empirical data-driven problem with several methods having been developed. However, the experience of the utilities over the past ten years with these methods indicates that the absence of physics-based information renders the data-driven approach less reliable. Complicating factors such as the inherent variability of operation (both equipment alignment and operating condition) can confound a pure data-driven approach while there are no rigorous guidelines for determining what constitutes an adequate sensor set. The solution under development to overcome these shortcomings supplements the data analytic method with process information in a so-called process-constrained data-analytic approach. Simple balance equations are written for generic components (e.g., mechanical pump, valve, and heat exchanger). These do not require a priori knowledge of process parameters, such as heat transfer coefficients or friction factors. All that is needed on the part of the utility user is to identify the components and how they are connected. This report describes the development of a parallel computing capability for determining the optimal sensor set. The optimal sensor set problem suffers from the curse of dimensionality. Computation time increases exponentially as the size of the system grows. To overcome this difficulty a pre-conditioner algorithm is developed to find an approximate solution close the actual solution. This serves as a seed for the full-blown algorithm and acts to constrain the space that must searched. The optimization algorithms are described and the implementation on a parallel computing platform is described. The application of the method to a use case we are solving in collaboration with our utility partner served to illustrate how the default sensor set in a nuclear plant may not provide sufficient coverage to infer sensor calibration status.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

PREDICTIONS OF ION PRODUCTION RATES AND ION NUMBER DENSITIES WITHIN THE DIAMAGNETIC CAVITY OF COMET 67P/CHURYUMOV-GERASIMENKO AT PERIHELION

We present a one-dimensional ion chemistry model of the diamagnetic cavity of comet 67P/Churyumov-Gerasimenko, the target comet for the ESA Rosetta mission. We solve the continuity equations for ionospheric species and predict number densities of electrons and selected ions considering only gas-phase reactions. We apply the model to the subsolar direction and consider conditions expected to be encountered by Rosetta at perihelion (1.29 AU) in 2015 August. Our default simulation predicts a maximum electron number density of {approx}8 Multiplication-Sign 10{sup 4} cm{sup -3} near the surface of the comet, while the electron number densities for cometocentric distances r > 10 km are approximately proportional to 1/r {sup 1.23} assuming that the electron temperature is equal to the neutral temperature. We show that even a small mixing ratio ({approx}0.3%-1%) of molecules having higher proton affinity than water is sufficient for the proton transfer from H{sub 3}O{sup +} to occur so readily that other ions than H{sub 3}O{sup +}, such as NH{sub 4} {sup +} or CH{sub 3}OH{sub 2} {sup +}, become dominant in terms of volume mixing ratio in part of, if not throughout, the diamagnetic cavity. Finally, we test how the predicted electron and ion densities are influenced by changes of model input parameters, including the neutral background, the impinging EUV solar spectrum, the solar zenith angle, the cross sections for photo- and electron-impact processes, the electron temperature profile, and the temperature dependence of ion-neutral reactions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Comparing Simulation Results with Traditional PRA Model on a Boiling Water Reactor Station Blackout Case Study

A previous study used RELAP and RAVEN to conduct a boiling water reactor station black-out (SBO) case study in a simulation based environment to show the capabilities of the risk-informed safety margin characterization methodology. This report compares the RELAP/RAVEN simulation results with traditional PRA model results. The RELAP/RAVEN simulation run results were reviewed for their input parameters and output results. The input parameters for each simulation run include various timing information such as diesel generator or offsite power recovery time, Safety Relief Valve stuck open time, High Pressure Core Injection or Reactor Core Isolation Cooling fail to run time, extended core cooling operation time, depressurization delay time, and firewater injection time. The output results include the maximum fuel clad temperature, the outcome, and the simulation end time. A traditional SBO PRA model in this report contains four event trees that are linked together with the transferring feature in SAPHIRE software. Unlike the usual Level 1 PRA quantification process in which only core damage sequences are quantified, this report quantifies all SBO sequences, whether they are core damage sequences or success (i.e., non core damage) sequences, in order to provide a full comparison with the simulation results. Three different approaches were used to solve event tree top events and quantify the SBO sequences: “W” process flag, default process flag without proper adjustment, and default process flag with adjustment to account for the success branch probabilities. Without post-processing, the first two approaches yield incorrect results with a total conditional probability greater than 1.0. The last approach accounts for the success branch probabilities and provides correct conditional sequence probabilities that are to be used for comparison. To better compare the results from the PRA model and the simulation runs, a simplified SBO event tree was developed with only four top events and eighteen SBO sequences (versus fifty-four SBO sequences in the original SBO model). The estimated SBO sequence conditional probabilities from the original SBO model were integrated to the corresponding sequences in the simplified SBO event tree. These results were then compared with the simulation run results.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2014

SU-E-T-247: Multi-Leaf Collimator Model Adjustments Improve Small Field Dosimetry in VMAT Plans

Purpose: The Elekta beam modulator linac employs a 4-mm micro multileaf collimator (MLC) backed by a fixed jaw. Out-of-field dose discrepancies between treatment planning system (TPS) calculations and output water phantom measurements are caused by the 1-mm leaf gap required for all moving MLCs in a VMAT arc. In this study, MLC parameters are optimized to improve TPS out-of-field dose approximations. Methods: Static 2.4 cm square fields were created with a 1-mm leaf gap for MLCs that would normally park behind the jaw. Doses in the open field and leaf gap were measured with an A16 micro ion chamber and EDR2 film for comparison with corresponding point doses in the Pinnacle TPS. The MLC offset table and tip radius were adjusted until TPS point doses agreed with photon measurements. Improvements to the beam models were tested using static arcs consisting of square fields ranging from 1.6 to 14.0 cm, with 45° collimator rotation, and 1-mm leaf gap to replicate VMAT conditions. Gamma values for the 3-mm distance, 3% dose difference criteria were evaluated using standard QA procedures with a cylindrical detector array. Results: The best agreement in point doses within the leaf gap and open field was achieved by offsetting the default rounded leaf end table by 0.1 cm and adjusting the leaf tip radius to 13 cm. Improvements in TPS models for 6 and 10 MV photon beams were more significant for smaller field sizes 3.6 cm or less where the initial gamma factors progressively increased as field size decreased, i.e. for a 1.6cm field size, the Gamma increased from 56.1% to 98.8%. Conclusion: The MLC optimization techniques developed will achieve greater dosimetric accuracy in small field VMAT treatment plans for fixed jaw linear accelerators. Accurate predictions of dose to organs at risk may reduce adverse effects of radiotherapy.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

MBARI-WEC September and October 2022 Field Data

This data is needed to simulate a model of the MBARI-WEC (Monterey Bay Aquarium Research Institute, Wave Energy Converter device) in a simulation environment (e.g. Gazebo) for 56 observation dates in the time between September and October 2022, and to compare the simulation outputs to the corresponding field data of the physical MBARI-WEC. There were 50 observations chosen in Sept and 6 observations in Oct. To help understand terms below, a summary of the system can be found at the github link in the downloads section below. The Gazebo MBARI-WEC model is also provided, should users wish to simulate using this platform. There are 4 *mat files included. ................................................................................................................................................................................................................................... Spectrum and Simulation Inputs: September2022_spectrum_siminputs.mat and October2022_spectrum_siminputs.mat has data needed for simulation inputs in table format. These include the ocean spectrum for an observation and operating parameters of the MBARI-WEC during that observation. They are organized as rows representing an observation and columns representing data. For example, for the September *mat there are 50 rows. The first 7 columns are Datetime, sig_waveheight, peak_period, mean_period, heaveconedoor_status, pistonpos_mean, and scale_factor: - Datetime is the date and time the observation occurred in PST - sig_waveheight is the significant wave height of the ocean spectrum during that observation in meters - peak_period is the peak period of the ocean spectrum during that observation in seconds - mean_period is the mean period of the ocean spectrum during that observation in seconds - heaveconedoor_status is the status of the heave cone doors where 0 represents the doors are open and 1 represents they are closed - pistonpos_mean is the mean position of the PTO ram (piston) in meters - scale_factor is an additional factor of 0.5 --1.4 applied to a default damping relationship The next columns are data needed to represent the ocean spectrum. First are the frequencies [Hz] labeled as "f0-f38", then the variance density [m2/Hz] labeled as "vardens0-vardens38". October2022_spectrum_siminputs.mat follows as a similar format as above, but includes a larger amount of ocean spectrum frequencies and variance density elements. ................................................................................................................................................................................................................................... Field data: The field data is found in MBARIWEC_septdata.mat and MBARIWEC_octdata.mat for the observations of September and October, respectively. These contain data in a struct format. The struct contains the following fields for each observation: PC_BattCurr, PC_LoadCurr, PC_RPM, PC_Voltage, SC_Range, SC_Velocity, DateTime, where: - PC_BattCurr is the current flowing to or from the onboard batteries in Amps - PC_LoadCurr is the current flowing to the load dump in Amps - PC_RPM is the electric/hydraulic motor shaft speed (directly coupled) in RPM - PC_Voltage is the bus voltage at the power converter in Volts - SC_Range is the PTO ram (piston) position in meters where 0 is fully retracted and 2.03 is fully extended - SC_Velocity is the PTO ram (piston) velocity in meters/sec - DateTime is the date and time of the sampled field data in each observation in PST - Electric Power is equal to: PC_Voltage*(PC_BattCurr + PC_LoadCurr) in Watts For example, upon loading MBARIWEC_octdata.mat, the aforementioned fields would be loaded, each with {6x1} cells for the 6 observations chosen in October. Within the first cell of e.g. SC_Range would be sampled data representing the field data of the MBARIWEC PTO piston position for say, one hour, of the first October observation. The corresponding field DateTime would...

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

IDAES-PSE 2.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.0.0 Release Highlights Removal of deprecated features from IDAES v1 Update to Pyomo v6.5 – this required a number of updates to support the new NL solver writer and to address some changes in Pyomo Creation of new testing suite for backward compatibility, model robustness and verification More general implementation of the Helmholtz EoS. This brings some new features like standard property diagrams, choice of mass or mole basis, and new state variable options Standardizing names in Heat Exchanger models (breaking change from v2.0.0a2): Control Volumes named hot_side and cold_side Ports names hot_side_inlet, hot_side_outlet, cold_side_inlet and cold_side_outlet Config Blocks names hot_side_config and cold_side_config Config arguments for user provided names for each side: hot_side_name and cold_side_name. Updating Keras surrogate tool to use v1.1 of OMLT New prototype API for model initialization (idaes.core.initialization) The new API uses "Model Initializer" objects instead of class methods, allowing for the definition of multiple initialization routines for a single model A number of common, model agnostic initialization routines have also been defined, including initialization from data, block-decomposition and a general hierarchical approach equivalent to the existing method for common unit models New metadata for thermophysical properties – valid_range This can be used to record the range of values over which a property value can be trusted, such as the range of experimental data used to regress parameters A number of new utility functions have been added to check for properties with values outside the valid range and to set bounds based on this metadata Updated construction of balance expressions in Control Volumes to remove unneeded terms In the past, unneeded terms were added as a constant 0 term, however they will now be dropped entirely from the expression This was necessary due to more strict unit checking in the new Pyomo solver writer which no longer ignores 0 terms Updates to metadata for thermophysical properties to better define known properties and units of measurement This results in more strict enforcement of standard naming for thermophysical and reaction properties Users can still define custom properties, but these must be done explicitly using the define_custom_properties() method instead of being implicitly created by add_property() Updated convergence tester utility tool to support definition of benchmark files (JSON format) and comparison of performance to benchmarks Set default iteration limit for IPOPT in IDAES config to 200 iterations Update scaling of example models to work with new Pyomo NL solver writer Improve testing of extensions and examples infrastructure to avoid need for downloading files Updated distillation column to centralize common functionality and remove a number of Pyomo warnings

IDAES↗
DOE OSTI2016

Redesign of the DFT/MRCI Hamiltonian

The combined density functional theory and multireference configuration interaction (DFT/MRCI) method of Grimme and Waletzke [J. Chem. Phys. 111, 5645 (1999)] is a well-established semi-empirical quantum chemical method for efficiently computing excited-state properties of organic molecules. As it turns out, the method fails to treat bi-chromophores owing to the strong dependence of the parameters on the excitation class. In this work, we present an alternative form of correcting the matrix elements of a MRCI Hamiltonian which is built from a Kohn-Sham set of orbitals. It is based on the idea of constructing individual energy shifts for each of the state functions of a configuration. The new parameterization is spin-invariant and incorporates less empirism compared to the original formulation. By utilizing damping techniques together with an algorithm of selecting important configurations for treating static electron correlation, the high computational efficiency has been preserved. The robustness of the original and redesigned Hamiltonians has been tested on experimentally known vertical excitation energies of organic molecules yielding similar statistics for the two parameterizations. Besides that, our new formulation is free from artificially low-lying doubly excited states, producing qualitatively correct and consistent results for excimers. The way of modifying matrix elements of the MRCI Hamiltonian presented here shall be considered as default choice when investigating photophysical processes of bi-chromophoric systems such as singlet fission or triplet-triplet upconversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Differentiable programming for online training of a neural artificial viscosity function within a staggered grid Lagrangian hydrodynamics scheme

Lagrangian methods to solve the inviscid Euler equations produce numerical oscillations near shock waves. A common approach to reducing these oscillations is to add artificial viscosity (AV) to the discrete equations. The AV term acts as a dissipative mechanism that attenuates oscillations by smearing the shock across a finite number of computational cells. However, AV introduces several control parameters that are not determined by the underlying physical model, and hence, in practice are tuned to the characteristics of a given problem. We seek to improve the standard quadratic-linear AV form by replacing it with a learned neural function that reduces oscillations relative to exact solutions of the Euler equations, resulting in a hybrid numerical-neural hydrodynamic solver. Because AV is an artificial construct that exists solely to improve the numerical properties of a hydrodynamic code, there is no offline ‘viscosity data’ against which a neural network can be trained before inserting into a numerical simulation, thus requiring online training. We achieve this via differentiable programming, i.e. end-to-end backpropagation or adjoint solution through both the neural and differential equation code, using automatic differentiation of the hybrid code in the Julia programming language to calculate the necessary loss function gradients. A novel offline pre-training step accelerates training by initializing the neural network to the default numerical AV scheme, which can be learned rapidly by space-filling sampling over the AV input space. We find that online training over early time steps of simulation is sufficient to learn a neural AV function that reduces numerical oscillations in long-term hydrodynamic shock simulations. These results offer an early proof-of-principle that online differentiable training of hybrid numerical schemes with novel neural network components can improve certain performance aspects existing in purely numerical schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

IDAES-PSE 2.6.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.6.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New Intersphinx extension automatically linking Jupyter notebook examples to project documentation New end-to-end diagnostics example demonstrated on a real problem New complementarity formulation for VLE with cubic equations of state, backward compatibility for old formulation New solver interface with presolve (ipopt_v2) in support of upcoming changes to the initialization and APIs methods, with default set to ipopt to maintain backwards compatibility; this will deprecate once all examples have been updated New forecaster and parameterized bidder methods within grid integration library Updated surrogates API and examples to support Keras 3, with backwards compatibility for older formats such as TensorFlow SavedModel (TFSM) Updated costing base dictionary to include the 2023 cost year index value Updated ProcessBlock to include information on the constructing block class Updated Flowsheet Visualizer to allow visualize() method to return value and functions Bug Fixes Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Fixed typos flagged by June update to crate-ci/typos and removed DMF-related exceptions Minor corrections of units of measurement handling in power plant waste/transport costing expressions, control volume material holdup expressions, and BTX property package parameters Fixed throwing >7500 numpy deprecation warnings by replacing scalar value assignment with element extraction and item iteration calls Testing and Robustness Migrated slow tests (>10s) to integration, impacting test coverage but also yielding a nearly 30% decrease in local test runtime Pinned pint to avoid issues with older supported Python versions Pinned codecov versions to avoid tokenless upload behavior with latest version Bumped extensions to version 3.4.2 to allow pointing to non-standard install location Deprecations and Removals Python 3.8 is no longer supported. The supported Python versions are 3.9 through 3.12 The Data Management Framework (DMF) is no longer supported. Importing idaes.core.dmf will cause a deprecation warning to be displayed until the next release The SOFC Keras surrogates have been removed. The current version of the SOFC surrogate model in the examples repository is a PySMO Kriging model.

AS↗
DOE OSTI2022

Gas puff imaging on the TCV tokamak

We present the design and operation of a suite of Gas Puff Imaging (GPI) diagnostic systems installed on the Tokamak à Configuration Variable (TCV) for the study of turbulence in the plasma edge and Scrape-Off-Layer (SOL). These systems provide the unique ability to simultaneously collect poloidal 2D images of plasma dynamics at the outboard midplane, around the X-point, in both the High-Field Side (HFS) and Low-Field Side (LFS) SOL, and in the divertor region. We describe and characterize an innovative control system for deuterium and helium gas injection, which is becoming the default standard for the other gas injections at TCV. Extensive pre-design studies and the different detection systems are presented, including an array of avalanche photodiodes and a high-speed CMOS camera. First results with spatial and time resolutions of up to ≈2 mm and 0.5 µs, respectively, are described, and future upgrades of the GPI diagnostics for TCV are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗