Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model checking”

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 415 records · Page 23

First results of electron temperature measurements by the use of multi-pass Thomson scattering system in GAMMA 10

A multi-pass Thomson scattering (TS) has the advantage of enhancing scattered signals. We constructed a multi-pass TS system for a polarisation-based system and an image relaying system modelled on the GAMMA 10 TS system. We undertook Raman scattering experiments both for the multi-pass setting and for checking the optical components. Moreover, we applied the system to the electron temperature measurements in the GAMMA 10 plasma for the first time. The integrated scattering signal was magnified by approximately three times by using the multi-pass TS system with four passes. The electron temperature measurement accuracy is improved by using this multi-pass system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Final Report on XStack: Software Synthesis for High Productivity ExaScale Computing

The goal of the project was to develop a programming model that would significantly improve productivity in the high-performance computing domain by bringing together three components: a) Automated equivalence checking, b) Sketch-based program synthesis, and c) Autotuning. The report provides an executive summary of the research accomplished through this project. At the end of the report is appended a paper that describes in more detail the key technical accomplishments from this project, and which was published in SC 2014.

97 MATHEMATICS AND COMPUTING↗

Modeling Transverse Space Charge effects in IOTA with pyORBIT

The role and mitigation of space charge effects are important aspects of the beam physics research to be performed at the Integrable Optics Test Accelerator (IOTA) at Fermilab. The impact of nonlinear integrability (partial and complete) on space charge driven incoherent and coherent resonances will be a new feature of this accelerator and has the potential to influence the design of future low energy proton synchrotrons. In this report we model transverse space charge effects using the PIC code pyORBIT. First we benchmark the single particle tracking results against MADX with checks on symplecticity, tune footprints, and dynamic aperture in a partially integrable lattice realized with a special octupole string insert. Our space charge calculations begin with an examination of the 4D symplecticity. Short term tracking is done first with the initial bare lattice and then with a matching of the rms values with space charge. Next, we explore slow initialization of charge as a technique to establish steady state and reduce emittance growth and beam loss following injection into a mismatched lattice. We establish values of space charge simulation parameters so as to ensure numerical convergence. Finally, we compare the simulated space charge induced tune shifts and footprints against theory.

43 PARTICLE ACCELERATORS↗

Lorentz symmetry and IR structure of the BFSS matrix model

The BFSS matrix model relates flat space M-theory to a large N limit of matrix quantum mechanics describing N non-relativistic D0-branes. M-theory, being a theory of gravity in flat space, has a rich infrared structure that includes various soft theorems and an infinite set of conserved charges associated to asymptotic symmetries. In this work, we ask: to what extent is this infrared structure present in BFSS? We find that all the salient features concerning the infrared structure of M-theory carry over naturally to the quantum mechanics dual. Moreover, we demonstrate that the dual statement of the soft graviton theorem in the matrix model implies that D0-brane scattering amplitudes in BFSS enjoy the full 11d Lorentz symmetry of M-theory, a claim which has been long anticipated. We also offer several first-principle consistency checks for our findings, including a computation of the soft theorem which does not presuppose the BFSS duality and a non-trivial match between several known symmetries of M-theory and BFSS that appear naturally in this formalism. These calculations give non-perturbative evidence in support of the BFSS duality as a model of flat space holography.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Relativistic strange stars in Tolman–Kuchowicz spacetime

Highlights: • We study relativistic strange star using Tolman–Kuchowicz metric. • Bag constant (55–75) MeV/fm{sup 3} remains within the predicted range for stable strange star. • Using observed values for mass and radius, we calculate different physical parameters. • This model satisfies energy conditions, TOV equations and all other stability criteria. • Solutions are completely free from any singularity and give stable stellar system. - Abstract: In this article we propose a relativistic model of a static spherically symmetric anisotropic strange star with the help of Tolman–Kuchowicz (TK) metric potentials (Tolman, 1939 and Kuchowicz, 1968). The form of the potentials are λ(r)=ln(1+ar{sup 2}+br{sup 4}) and ν(r)=Br{sup 2}+2lnC where a, b, B and C are constants which we have to evaluate using boundary conditions. We also consider the simplest form of the phenomenological MIT bag equation of state (EOS) to represent the strange quark matter (SQM) distribution inside the stellar system. Here, the radial pressure p{sub r} relates with the density profile ρ as follows, p{sub r}(r)=(1/3)[ρ(r)−4B{sub g}], where B{sub g} is the Bag constant. To check the physical acceptability and stability of the stellar system based on the obtained solutions, we have performed various physical tests. It is shown that the model satisfies all the stability criteria, including nonsingular nature of the density and pressure, implies stable nature. Here, the Bag constant for different strange star candidates are found to be (68–70) MeV/fm{sup 3} which satisfies all the acceptability criteria and remains in the experimental range.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Investigation of methodology for uncertainty quantification of model parameters

Quantifying the uncertainty contributors of best-estimate thermal hydraulic code is getting more and more attention in safety analysis of nuclear industry during recent decades. Yet for uncertainties of certain model parameters which may not be directly measured, the quantification process is usually subjective and inaccurate. This paper investigates the statistical methodologies in order to get the probability density function (pdf) of model parameters more objectively based on observed experimental responses. The simplification of mathematical model is described for the parameter estimation, and solutions based on Bayesian theorem and Maximum Likelihood Estimation (MLE) are classified and demonstrated for both deterministic approach and stochastic approach. As an application, the uncertainties of model parameters related to reflood phenomena implemented in RELAP5 code are quantified. It is indicated that the statistical method which is independent of best-estimate code is feasible and easy to apply. And through the check of uncertainty propagation, it demonstrates that the uncertainty bands can envelope most of the experiment measurement with advantage of accuracy. (authors)

42 ENGINEERING↗

Two-loop free energy of three-dimensional antiferromagnets in external magnetic and staggered fields

Using a model-independent low-energy effective field theory, we calculate the free energy of three-dimensional antiferromagnets in a combination of mutually perpendicular external magnetic and staggered fields at the next-to-next-to-leading, two-loop order. Renormalization is carried out analytically, and the renormalization group invariance of the result is checked explicitly. The free energy is thus expressed solely in terms of temperature, the external fields, and a set of low-energy coupling constants, to be determined by experiment or by matching to the microscopic model of a given concrete material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Atacama Cosmology Telescope: DR6 gravitational lensing and SDSS BOSS cross-correlation measurement and constraints on gravity with the 𝐸 𝐺 statistic

We derive new constraints on the 𝐸 𝐺 statistic as a test of gravity, combining the cosmic microwave background (CMB) lensing map estimated from Data Release 6 (DR6) of the Atacama Cosmology Telescope with Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (SDSS BOSS) CMASS and LOWZ galaxy data. We develop an analysis pipeline to measure the cross-correlation between CMB lensing maps and galaxy data, following a blinding policy and testing the approach through null and consistency checks. By testing the equivalence of the spatial and temporal gravitational potentials, the 𝐸 𝐺 statistic can distinguish Λ⁢ CDM from alternative models of gravity. We find 𝐸 𝐺 ⁡(𝑧 eff = 0.555) = 0.3⁢1$^{+0.06}_{−0.05}$ for Atacama Cosmology Telescope (ACT) and CMASS data at 68.28% confidence level, and 𝐸 𝐺 ⁡(𝑧 eff = 0.316) = 0.4⁢9$^{+0.14}_{−0.11}$ for the ACT and LOWZ. Systematic errors are estimated to be 3% and 4%, respectively. Including CMB lensing information from Planck PR4 results in 𝐸 𝐺 ⁡(𝑧 eff = 0.555) = 0.3⁢4$^{+0.05}_{−0.05}$ with CMASS and 𝐸 𝐺 ⁡(𝑧 eff = 0.316) = 0.4⁢3$^{+0.11}_{−0.09}$ with LOWZ. These are consistent with predictions for the Λ⁢ CDM model that best fits the Planck CMB anisotropy and SDSS BOSS baryon acoustic oscillations (BAO), where 𝐸$^{GR}_{𝐺⁡}$(𝑧 eff =0.555) =0.401 ± 0.005 for CMB lensing combined with CMASS and 𝐸$^{GR}_{𝐺}$⁡(𝑧 eff = 0.316) = 0.452 ± 0.005 combined with LOWZ. We also find 𝐸 𝐺 to be scale independent, with probability to exceed >5%, as predicted by general relativity. The methods developed in this work are also applicable to improved future analyses with upcoming spectroscopic galaxy samples and CMB lensing measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

TU-C-17A-03: An Integrated Contour Evaluation Software Tool Using Supervised Pattern Recognition for Radiotherapy

Purpose: Radiotherapy (RT) contours delineated either manually or semiautomatically require verification before clinical usage. Manual evaluation is very time consuming. A new integrated software tool using supervised pattern contour recognition was thus developed to facilitate this process. Methods: The contouring tool was developed using an object-oriented programming language C# and application programming interfaces, e.g. visualization toolkit (VTK). The C# language served as the tool design basis. The Accord.Net scientific computing libraries were utilized for the required statistical data processing and pattern recognition, while the VTK was used to build and render 3-D mesh models from critical RT structures in real-time and 360° visualization. Principal component analysis (PCA) was used for system self-updating geometry variations of normal structures based on physician-approved RT contours as a training dataset. The inhouse design of supervised PCA-based contour recognition method was used for automatically evaluating contour normality/abnormality. The function for reporting the contour evaluation results was implemented by using C# and Windows Form Designer. Results: The software input was RT simulation images and RT structures from commercial clinical treatment planning systems. Several abilities were demonstrated: automatic assessment of RT contours, file loading/saving of various modality medical images and RT contours, and generation/visualization of 3-D images and anatomical models. Moreover, it supported the 360° rendering of the RT structures in a multi-slice view, which allows physicians to visually check and edit abnormally contoured structures. Conclusion: This new software integrates the supervised learning framework with image processing and graphical visualization modules for RT contour verification. This tool has great potential for facilitating treatment planning with the assistance of an automatic contour evaluation module in avoiding unnecessary manual verification for physicians/dosimetrists. In addition, its nature as a compact and stand-alone tool allows for future extensibility to include additional functions for physicians’ clinical needs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

ECP-copa/ExaMiniMD

ExaMiniMD is a proxy application for Molecular Dynamics (MD) codes such as LAMMPS. Its intend is to facilitate research and development on MD algorithms in particular in the context of CoDesign efforts with hardware vendors. Currently, ExaMiniMD implements a basic Lennard Jones two body potential as well as the SNAP many body potential, both of which are available in LAMMPS. It also includes basic neighbor list builds as well as constant energy integration as found in LAMMPS and other MD codes. A feature not generally available in other MD codes is a set of statistical correctness checks which will help other developers judge whether changes to the code still give correct answers. ExaMiniMD is using the Kokkos programming model to achieve performance portability and is able to run across all major HPC platforms (GPUs, ManyCore, CPUs).

Thompson, Aidan↗

Three-dimensional Magnetohydrodynamical Simulations of the Morphology of Head–Tail Radio Galaxies Based on the Magnetic Tower Jet Model

The distinctive morphology of head–tail radio galaxies reveals strong interactions between the radio jets and their intra-cluster environment, the general consensus on the morphology origin of head–tail sources is that radio jets are bent by violent intra-cluster weather. We demonstrate in this paper that such strong interactions provide a great opportunity to study the jet properties and also the dynamics of the intra-cluster medium (ICM). By three-dimensional magnetohydrodynamical simulations, we analyze the detailed bending process of a magnetically dominated jet, based on the magnetic tower jet model. We use stratified atmospheres modulated by wind/shock to mimic the violent intra-cluster weather. Core sloshing is found to be inevitable during the wind-cluster core interaction, which induces significant shear motion and could finally drive ICM turbulence around the jet, making it difficult for the jet to survive. We perform a detailed comparison between the behavior of pure hydrodynamical jets and the magnetic tower jet and find that the jet-lobe morphology could not survive against the violent disruption in all of our pure hydrodynamical jet models. On the other hand, the head–tail morphology is well reproduced by using a magnetic tower jet model bent by wind, in which hydrodynamical instabilities are naturally suppressed and the jet could always keep its integrity under the protection of its internal magnetic fields. Finally, we also check the possibility for jet bending by shock only. We find that shock could not bend the jet significantly, and thus could not be expected to explain the observed long tails in head–tail radio galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of the Nuclear Equation of State on Relativistic Turbulence-induced Core-collapse Supernovae

Abstract The nuclear equation of state (EOS) is an important component in the evolution of core-collapse supernovae. In this paper we make a survey of various EOSs in the literature and analyze their effect on spherical core-collapse models in which the effects of three-dimensional turbulence is modeled by a general relativistic formulation of Supernova Turbulence In Reduced-dimensionality (STIR). We show that the viability of the explosion is quite EOS dependent and that it best correlates with the early-time interior entropy density of the proto–neutron star. We check that this result is not progenitor dependent, although the lowest-mass progenitors show different explosion properties, due to the different pre-collapse nuclear composition. Larger central entropies also induce more vigorous proto–neutron star convection in our one-dimensional turbulence model, as well as a wider convective layer.

79 ASTRONOMY AND ASTROPHYSICS↗

SALT3: An Improved Type Ia Supernova Model for Measuring Cosmic Distances

Abstract A spectral-energy distribution (SED) model for Type Ia supernovae (SNe Ia) is a critical tool for measuring precise and accurate distances across a large redshift range and constraining cosmological parameters. We present an improved model framework, SALT3, which has several advantages over current models—including the leading SALT2 model (SALT2.4). While SALT3 has a similar philosophy, it differs from SALT2 by having improved estimation of uncertainties, better separation of color and light-curve stretch, and a publicly available training code. We present the application of our training method on a cross-calibrated compilation of 1083 SNe with 1207 spectra. Our compilation is 2.5× larger than the SALT2 training sample and has greatly reduced calibration uncertainties. The resulting trained SALT3.K21 model has an extended wavelength range 2000–11,000 Å (1800 Å redder) and reduced uncertainties compared to SALT2, enabling accurate use of low- z I and iz photometric bands. Including these previously discarded bands, SALT3.K21 reduces the Hubble scatter of the low- z Foundation and CfA3 samples by 15% and 10%, respectively. To check for potential systematic uncertainties, we compare distances of low (0.01 < z < 0.2) and high (0.4 < z < 0.6) redshift SNe in the training compilation, finding an insignificant 3 ± 14 mmag shift between SALT2.4 and SALT3.K21. While the SALT3.K21 model was trained on optical data, our method can be used to build a model for rest-frame NIR samples from the Roman Space Telescope. Our open-source training code, public training data, model, and documentation are available at https://saltshaker.readthedocs.io/en/latest/ , and the model is integrated into the sncosmo and SNANA software packages.

79 ASTRONOMY AND ASTROPHYSICS↗

2021 Summer SPE Project

SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.

58 GEOSCIENCES↗

Re-evaluation of Spent Nuclear Fuel Assay Data for the Three Mile Island Unit 1 Reactor and Application to Code Validation

Destructive radiochemical assay measurements of spent nuclear fuel rod segments from an assembly irradiated in the Three Mile Island unit 1 (TMI-1) pressurized water reactor have been performed at Oak Ridge National Laboratory (ORNL). Assay data are reported for five samples from two fuel rods of the same assembly. The TMI-1 assembly was a 15 X 15 design with an initial enrichment of 4.013 wt% 235 U, and the measured samples achieved burnups between 45.5 and 54.5 gigawatt days per metric ton of initial uranium (GWd/t). Measurements were performed mainly using inductively coupled plasma mass spectrometry after elemental separation via high performance liquid chromatography. High precision measurements were achieved using isotope dilution techniques for many of the lanthanides, uranium, and plutonium isotopes. Measurements are reported for more than 50 different isotopes and 16 elements. One of the two TMI-1 fuel rods measured in this work had been measured previously by Argonne National Laboratory (ANL), and these data have been widely used to support code and nuclear data validation. Recently, ORNL provided an important opportunity to independently cross check results against previous measurements performed at ANL. The measured nuclide concentrations are used to validate burnup calculations using the SCALE nuclear systems modeling and simulation code suite. These results show that the new measurements provide reliable benchmark data for computer code validation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SU-G-TeP4-06: An Integrated Application for Radiation Therapy Treatment Plan Directives, Management, and Reporting

Purpose: With electronic medical records, patient information for the treatment planning process has become disseminated across multiple applications with limited quality control and many associated failure modes. We present the development of a single application with a centralized database to manage the planning process. Methods: The system was designed to replace current functionalities of (i) static directives representing the physician intent for the prescription and planning goals, localization information for delivery, and other information, (ii) planning objective reports, (iii) localization and image guidance documents and (iv) the official radiation therapy prescription in the medical record. Using the Eclipse Scripting Application Programming Interface, a plug-in script with an associated domain-specific SQL Server database was created to manage the information in (i)–(iv). The system’s user interface and database were designed by a team of physicians, clinical physicists, database experts, and software engineers to ensure usability and robustness for clinical use. Results: The resulting system has been fully integrated within the TPS via a custom script and database. Planning scenario templates, version control, approvals, and logic-based quality control allow this system to fully track and document the planning process as well as physician approval of tradeoffs while improving the consistency of the data. Multiple plans and prescriptions are supported along with non-traditional dose objectives and evaluation such as biologically corrected models, composite dose limits, and management of localization goals. User-specific custom views were developed for the attending physician review, physicist plan checks, treating therapists, and peer review in chart rounds. Conclusion: A method was developed to maintain cohesive information throughout the planning process within one integrated system by using a custom treatment planning management application that interfaces directly with the TPS. Future work includes quantifying the improvements in quality, safety and efficiency that are possible with the routine clinical use of this system. Supported in part by NIH-P01-CA-059827.

60 APPLIED LIFE SCIENCES↗

pnnl/ruleset-checking-tool

The Ruleset Checking Tool (RCT) will be a tool for verifying implementation of a ruleset (such as ASHRAE Standard 90.1 Appendix G) in building performance modeling (BPM) tools.

McNeill, James↗

Efforts to enhance reproducibility in a human performance research project

Background: Ensuring the validity of results from funded programs is a critical concern for agencies that sponsor biological research. In recent years, the open science movement has sought to promote reproducibility by encouraging sharing not only of finished manuscripts but also of data and code supporting their findings. While these innovations have lent support to third-party efforts to replicate calculations underlying key results in the scientific literature, fields of inquiry where privacy considerations or other sensitivities preclude the broad distribution of raw data or analysis may require a more targeted approach to promote the quality of research output. Methods: We describe efforts oriented toward this goal that were implemented in one human performance research program, Measuring Biological Aptitude, organized by the Defense Advanced Research Project Agency's Biological Technologies Office. Our team implemented a four-pronged independent verification and validation (IV&V) strategy including 1) a centralized data storage and exchange platform, 2) quality assurance and quality control (QA/QC) of data collection, 3) test and evaluation of performer models, and 4) an archival software and data repository. Results: Our IV&V plan was carried out with assistance from both the funding agency and participating teams of researchers. QA/QC of data acquisition aided in process improvement and the flagging of experimental errors. Holdout validation set tests provided an independent gauge of model performance. Conclusions: In circumstances that do not support a fully open approach to scientific criticism, standing up independent teams to cross-check and validate the results generated by primary investigators can be an important tool to promote reproducibility of results.

59 BASIC BIOLOGICAL SCIENCES↗