Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal Methods”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Computational Aspects of Single-Molecule Kinetics for Coupled Catalytic Cycles: A Spectral Analysis

Catalysis from single active sites is analyzed using methods developed from single molecule kinetics. Using a stochastic Markov state description, the observable properties of general catalytic networks of reactions are expressed using an eigenvalue decomposition of the transition matrix for the Markov process. By the use of a sensitivity analysis, the necessary eigenvalues and eigenvectors are related to the energies of controlling barriers and wells located along the reaction routes. A generalization of the energetic span theory allows the eigenvalues to be computed from several activation energies corresponding to distinct barrier-well pairings. The formalism is demonstrated accurately for model problems for a physically realistic mechanism for an alkene hydrogenation reaction on a single atom catalyst.

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

A summary of the Hanford single-shell tank structural analysis of record - 15526

Washington River Protection Solutions, LLC (WRPS), the US DOE Hanford Site tank contractor, developed an enhanced Single-Shell Tank Integrity Project (SSTIP) in 2009 to better understand the integrity of the single-shell tanks (SSTs). A panel of experts from industry and academia was created to recommend project activities in four main areas: structural integrity, liner degradation, leak integrity and prevention, and contamination migration. One of the panel's primary recommendations was to perform a detailed structural integrity analysis of the SSTs. WRPS and the U.S. Department of Energy, Office of River Protection commissioned the Pacific Northwest National Laboratory (PNNL) to perform this structural integrity analysis of record for the Hanford single shell tanks. This is one of seventeen panel recommendations that form the M-45-10-1 Change Package for the Hanford Federal Agreement and Compliance Order; also known as the Tri-Party Agreement. Finite element analysis was used to predict the structural response of the SSTs to the historical thermal and operating loads, plus design basis seismic loads. Bounding thermal histories were established from waste temperature records and applied to the model to include the thermal degradation of concrete modulus and strength, plus cracking due to differential thermal expansion while under in situ static loading conditions. PNNL contracted with Becht Engineering to evaluate the dynamic response of the tanks to a design basis earthquake. The degraded stiffness of the reinforced concrete was also incorporated in the seismic analysis. The combined response to static and seismic loads was then evaluated against the design requirements of the American Concrete Institute standard, ACI-349-06, for nuclear safety-related concrete structures. Specific analyses were performed for the four major tank designs.This paper presents an overview of the analysis methods used and it summarizes the major analytical findings that are common to the tanks and others that are specific to the individual tank designs. The analyses of record determined that including the bounding effects of their operating histories, each tank design currently satisfies the ACI-349-06 structural integrity design requirements. Although leak tightness is not directly addressed in the structural integrity analysis, this analysis will be a key input in the formal tank integrity assessment due in 2018. (authors)

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

SU-E-T-757: TMRs Calculated From PDDs Versus the Direct Measurements for Small Field SRS Cones

Purpose: To investigate the variation of TMR for SRS cones obtained by TMR scanning, calculation from PDDs, and point measurements. The obtained TMRs were also compared to the representative data from the vendor. Methods: TMRs for conical cones of 4, 5, 7.5, 10, 12.5, 15, and 17.5 mm diameter (jaws set to 5×5 cm) were obtained for 6X FFF and 10X FFF energies on a Varian Edge linac. TMR scanning was performed with a Sun Nuclear 3D scanner and Edge detector at 100 cm SDD. TMR point measurements were measured with a Wellhofer tank and Edge detector, at multiple depths from 0.5 to 20 cm and 100 cm SDD. PDDs for converting to TMR were scanned with a Wellhofer system and SFD detector. The formalism of converting PDD to TMR, given in Khan’s book (4th Edition, p.161) was applied. Sp values at dmax were obtained by measuring Scp and Sc of the cones (jaws set to 5×5 cm) using the Edge detector, and normalized to the 10×10 cm field. Results: Along the central axis beyond dmax, the RMS and maximum percent difference of TMRs obtained with different methods were as follows: (a) 1.3% (max=3.5%) for the calculated TMRs from PDDs versus direct scanning; (b) 1.2% (max=3.3%) for direct scanning versus point measurement; (c) 1.8% (max=5.1%) for the calculated versus point measurements; (d) 1.0% (max=3.6%) for direct scanning versus vendor data; (e) 1.6% (max=7.2%) for the calculated versus vendor data. Conclusion: The overall accuracy of TMRs calculated from PDDs was comparable with that of direct scanning. However, the uncertainty at depths greater than 20 cm, increased up to 5% when compared to point measurements. This issue must be considered when developing a beam model for small field SRS planning using cones.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Pesticide risk perceptions and the differences between farmers and extensionists: Towards a knowledge-in-context model

A growing body of literature analyzes farmer perceptions of pesticide risk, but much less attention has been given to differences in risk perception between farmers and technical experts. Furthermore, inconsistencies in knowledge have too easily been explained in terms of lack of knowledge rather than exploring the underlying reasons for particular forms of thinking about pesticide risks. By doing this, the division between expert and lay knowledge has been deepened rather than transcended. Objective: This study aims to understand differences and similarities among the perceptions of pesticide risks of farmers, farm workers, and technical experts such as extensionists, by applying a social science approach towards knowledge and risk attitudes. Methods: Semi-structured interviews and field observations were conducted to smallholders, farm workers, extensionists, health professionals and scientists involved in the use and handling of pesticides. Subsequently, a survey was carried out to quantify the farmers and extensionists' acceptance or rejection of typical assertions expressed previously in the semi-structured interviews. Results: Smallholders showed to gain knowledge from their own experiences and to adapt pesticides practices, which is a potential basis for transforming notions of pesticide safety and risk reduction strategies. Though extensionists have received formal education, they sometimes develop ideas deviating from the technical perspective. The risk perception of the studied actors appeared to vary according to their role in the agricultural labor process; they varied much less than expected according to their schooling level. Conclusions: Commitment to the technical perspective is not dramatically different for extensionists on the one hand and farmers as well as farm workers on the other hand. Ideas about a supposed lack of knowledge by farmers and the need of formal training are too much driven by a deficit model of knowledge. Further research on risk perceptions of pesticides and training of rural people will benefit from the development of a knowledge-in-context model. -- Highlights: • Researching perceptions of farmers' extensionists and other professionals. • Experts as well as farmers deviate from the technical perspective. • Blaming who is responsible for pesticide problems creates expert-lay division. • Qualitative and quantitative methods, not as complementary but integrated. • Knowledge-in-context model as an alternative to the knowledge-deficit model.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Global Sensitivity Measures from Given Data

Simulation models support managers in the solution of complex problems. International agencies recommend uncertainty and global sensitivity methods as best practice in the audit, validation and application of scientific codes. However, numerical complexity, especially in the presence of a high number of factors, induces analysts to employ less informative but numerically cheaper methods. This work introduces a design for estimating global sensitivity indices from given data (including simulation input–output data), at the minimum computational cost. We address the problem starting with a statistic based on the L1-norm. A formal definition of the estimators is provided and corresponding consistency theorems are proved. The determination of confidence intervals through a bias-reducing bootstrap estimator is investigated. The strategy is applied in the identification of the key drivers of uncertainty for the complex computer code developed at the National Aeronautics and Space Administration (NASA) assessing the risk of lunar space missions. We also introduce a symmetry result that enables the estimation of global sensitivity measures to datasets produced outside a conventional input–output functional framework.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2015

Real-time and imaginary-time quantum hierarchal Fokker-Planck equations

We consider a quantum mechanical system represented in phase space (referred to hereafter as “Wigner space”), coupled to a harmonic oscillator bath. We derive quantum hierarchal Fokker-Planck (QHFP) equations not only in real time but also in imaginary time, which represents an inverse temperature. This is an extension of a previous work, in which we studied a spin-boson system, to a Brownian system. It is shown that the QHFP in real time obtained from a correlated thermal equilibrium state of the total system possesses the same form as those obtained from a factorized initial state. A modified terminator for the hierarchal equations of motion is introduced to treat the non-Markovian case more efficiently. Using the imaginary-time QHFP, numerous thermodynamic quantities, including the free energy, entropy, internal energy, heat capacity, and susceptibility, can be evaluated for any potential. These equations allow us to treat non-Markovian, non-perturbative system-bath interactions at finite temperature. Through numerical integration of the real-time QHFP for a harmonic system, we obtain the equilibrium distributions, the auto-correlation function, and the first- and second-order response functions. These results are compared with analytically exact results for the same quantities. This provides a critical test of the formalism for a non-factorized thermal state and elucidates the roles of fluctuation, dissipation, non-Markovian effects, and system-bath coherence. Employing numerical solutions of the imaginary-time QHFP, we demonstrate the capability of this method to obtain thermodynamic quantities for any potential surface. It is shown that both types of QHFP equations can produce numerical results of any desired accuracy. The FORTRAN source codes that we developed, which allow for the treatment of Wigner space dynamics with any potential form (TanimuranFP15 and ImTanimuranFP15), are provided as the supplementary material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Source size of x rays from self-modulated laser wakefield accelerators

A comparative study of x-ray sources generated with different mechanisms from self-modulated laser wakefield acceleration (SM-LWFA) electrons was performed to compare the source size or spatial resolution for use in high energy density science applications. We examine the source size of betatron, inverse Compton scattering, and bremsstrahlung radiation with a Fresnel diffraction based formalism and a modified x-ray ray tracing model. We observe the dependence of source size on the radiation generation process, laser parameters, and compare to what is possible in other regimes of LWFA, as well as current methods. This information is significant as we begin to explore the use of light sources driven by SM-LWFA for use as a diagnostic at large-scale laser facilities where blowout regime LWFA is not possible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Thermal conductivity in large-J two-dimensional antiferromagnets: Role of phonon scattering

Motivated by the recent heat transport experiments in 2D antiferromagnets, such as La 2 CuO 4 , where the exchange coupling J is larger than the Debye energy Θ D , we discuss different types of relaxation processes for magnon heat current with a particular focus on coupling to 3D phonons. In this work, we study thermal conductivity by these in-plane magnetic excitations using two distinct techniques, Boltzmann formalism within the relaxation-time approximation and memory-function approach. Within these approaches, a close consideration is given to the scattering of magnons by both acoustic and optical branches of phonons. A remarkable accord between the two methods with regards to the asymptotic behavior of the effective relaxation rates is demonstrated. Additional scattering mechanisms, due to grain boundaries, impurities, and finite correlation length in the paramagnetic phase, are discussed and included in the calculations of the thermal conductivity κ(T). Again, we demonstrate a close similarity of the results from the two techniques of calculating κ(T). Our complementary approach strongly suggests that scattering from optical or zone-boundary phonons is important for magnon heat current relaxation in a high temperature window of Θ D ≲T$\ll$J.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Convolution Lagrangian perturbation theory for biased tracers beyond general relativity

We compare analytic predictions for real and Fourier space two-point statistics for biased tracers from a variety of Lagrangian perturbation theory approaches against those from state of the art N-body simulations in f(R) Hu-Sawicki and the nDGP braneworld modified gravity theories. We show that the novel physics of gravitational collapse in scalar tensor theories with the chameleon or the Vainshtein screening mechanism can be effectively factored in with bias parameters analytically predicted using the peak-background split formalism when updated to include the environmental sensitivity of modified gravity theories as well as changes to the halo mass function. Here, we demonstrate that convolution Lagrangian perturbation theory (CLPT) and standard perturbation theory (SPT) approaches provide accurate analytic methods to predict the correlation function and power spectra, respectively, for biased tracers in modified gravity models and are able to characterize both the baryon acoustic oscillation, power-law, and small scale regimes needed for upcoming galaxy surveys such as DESI, Euclid, LSST and WFIRST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Summary of LANL Critical Benchmark Comparison Study and Revisions for Cases Involving HEU, LEU, MIX, and Pu

This report documents results obtained for revisions made to cases involving Highly Enriched Uranium (HEU), Intermediate Enriched Uranium (IEU), a mixture of Pu and Uranium (MIX), as well as Pu cases. A previous summary of revisions for HEU an Pu cases was reported and additional investigations into four cases originally presented therein uncovered further revisions which led to better agreement with other transport codes, those cases are updated in this report. The summary of all cases reported in Reference 2 is updated in this report. In addition, a previous summary of revisions for LEU and MIX was reported, a summary of those revisions in reproduced in this report for a comprehensive summary of changes to benchmarks beginning in fiscal year 2020 to current date. The report focuses on the changes made to LANL benchmarks modeled with MCNP6 using ENDF/B-VII.1 nuclear data that appeared to have discrepant results when compared with results of other codes. Feedback was used to pinpoint review of benchmark input files and to revise them when necessary. This report documents the results of review and revision of specific benchmarks highlighted as possibly discrepant in the comparison study. In addition, there is an effort tied to this work involving collaboration between LANL XCP and NCS Divisions in the development of a shared review/revision procedure and use of a new benchmark repository. LANL has a benchmark library of critical experiments from the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook modeled for use with MCNP. This collection is now over 1100 benchmarks, referred to as the Whisper-1.1 library because it is used with the sensitivity/uncertainty package, Whisper, which supports nuclear criticality safety validation and is released with MCNP6.2. The collection, originally created several decades ago, is a combination of smaller collections, which has been revised and expanded, by various groups at LANL over the years. The original authors are no longer at the laboratory and little formal documentation of review and revision of these benchmarks exists today. A branch of the benchmark collection was already the subject of a formal review undertaken by the LANL NCS Division and expanded to include XCP Division.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Effect of Impurities on the Formation Energy of Point Defects in the γ-TiAl Alloy

The projector augmented-wave method within density functional theory is used to study the effect of 4d elements and Group IIIA and IVA elements on the energetics of point defect formation and elastic moduli of the γ-TiAl alloy. To calculate the point defect formation energy, the grand canonical formalism was used. It is shown that the formation energy of aluminum vacancies decreases by approximately 1.3 eV with increasing its content in the series of alloys Ti{sub 3}Al–TiAl–TiAl{sub 3}, whereas the formation energy of titanium vacancies changes insignificantly. In general, the formation of titanium vacancies in these alloys is preferred to the formation of aluminum vacancies. It has been found that Nb, Mo, Tc, Ru, Rh, and Pd impurities on the aluminum sublattice contribute to an increase in the formation energy of an aluminum vacancies, and Mo and Tc also lead to an increase in the formation energy of titanium vacancies. All 4d impurities, if they substitute for titanium, reduce the formation energy of aluminum vacancies, and Nb and Mo are favorable for increasing the formation energy of titanium vacancies. The influence of impurities on the chemical bond in the γ-TiAl alloy and its elastic moduli and characteristics based on them is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Parallelized domain decomposition for multi-dimensional Lagrangian random walk mass-transfer particle tracking schemes

Lagrangian particle tracking schemes allow a wide range of flow and transport processes to be simulated accurately, but a major challenge is numerically implementing the inter-particle interactions in an efficient manner. This article develops a multi-dimensional, parallelized domain decomposition (DDC) strategy for mass-transfer particle tracking (MTPT) methods in which particles exchange mass dynamically. We show that this can be efficiently parallelized by employing large numbers of CPU cores to accelerate run times. In order to validate the approach and our theoretical predictions we focus our efforts on a well-known benchmark problem with pure diffusion, where analytical solutions in any number of dimensions are well established. In this work, we investigate different procedures for “tiling” the domain in two and three dimensions (2-D and 3-D), as this type of formal DDC construction is currently limited to 1-D. An optimal tiling is prescribed based on physical problem parameters and the number of available CPU cores, as each tiling provides distinct results in both accuracy and run time. We further extend the most efficient technique to 3-D for comparison, leading to an analytical discussion of the effect of dimensionality on strategies for implementing DDC schemes. Increasing computational resources (cores) within the DDC method produces a trade-off between inter-node communication and on-node work. For an optimally subdivided diffusion problem, the 2-D parallelized algorithm achieves nearly perfect linear speedup in comparison with the serial run-up to around 2700 cores, reducing a 5 h simulation to 8 s, while the 3-D algorithm maintains appreciable speedup up to 1700 cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Few-Body Bound States and Resonances in Finite Volume

Since the pioneering work of Lüscher in the 1980s it is well known that considering quantum systems in finite volume, specifically, finite periodic boxes, can be used as a powerful computational tool to extract physical observables. While this formalism has been worked out in great detail in the two-body sector, much effort is currently being invested into deriving analogous relations for systems with more constituents. This work is relevant not only for nuclear physics, where lattice methods are now able to calculate few- and many-nucleon states, but also for other fields such as simulations of cold atoms. This article discusses recent progress regarding the extraction of few-body bound-state and resonance properties from finite-volume calculations of systems with an arbitrary number of constituents.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quantification of human contribution to soil moisture-based terrestrial aridity

Abstract Current knowledge of the spatiotemporal patterns of changes in soil moisture-based terrestrial aridity has considerable uncertainty. Using Standardized Soil Moisture Index (SSI) calculated from multi-source merged data sets, we find widespread drying in the global midlatitudes, and wetting in the northern subtropics and in spring between 45°N–65°N, during 1971–2016. Formal detection and attribution analysis shows that human forcings, especially greenhouse gases, contribute significantly to the changes in 0–10 cm SSI during August–November, and 0–100 cm during September–April. We further develop and apply an emergent constraint method on the future SSI’s signal-to-noise (S/N) ratios and trends under the Shared Socioeconomic Pathway 5-8.5. The results show continued significant presence of human forcings and more rapid drying in 0–10 cm than 0–100 cm. Our findings highlight the predominant human contributions to spatiotemporally heterogenous terrestrial aridification, providing a basis for drought and flood risk management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Estimating statistical isotropy violation in CMB due to non-circular beam and complex scan in minutes

Mild, unavoidable deviations from circular-symmetry of instrumental beams along with scan strategy can give rise to measurable Statistical Isotropy (SI) violation in Cosmic Microwave Background (CMB) experiments. If not accounted properly, this spurious signal can complicate the extraction of other SI violation signals (if any) in the data. However, estimation of this effect through exact numerical simulation is computationally intensive and time consuming. A generalized analytical formalism not only provides a quick way of estimating this signal, but also gives a detailed understanding connecting the leading beam anisotropy components to a measurable BipoSH characterisation of SI violation. In this paper, we provide an approximate generic analytical method for estimating the SI violation generated due to a non-circular (NC) beam and arbitrary scan strategy, in terms of the Bipolar Spherical Harmonic (BipoSH) spectra. Our analytical method can predict almost all the features introduced by a NC beam in a complex scan and thus reduces the need for extensive numerical simulation worth tens of thousands of CPU hours into minutes long calculations. As an illustrative example, we use WMAP beams and scanning strategy to demonstrate the easability, usability and efficiency of our method. We test all our analytical results against that from exact numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Cosmic shear power spectra in practice

Cosmic shear is one of the most powerful probes of Dark Energy, targeted by several current and future galaxy surveys. Lensing shear, however, is only sampled at the positions of galaxies with measured shapes in the catalog, making its associated sky window function one of the most complicated amongst all projected cosmological probes of inhomogeneities, as well as giving rise to inhomogeneous noise. Partly for this reason, cosmic shear analyses have been mostly carried out in real-space, making use of correlation functions, as opposed to Fourier-space power spectra. Since the use of power spectra can yield complementary information and has numerical advantages over real-space pipelines, it is important to develop a complete formalism describing the standard unbiased power spectrum estimators as well as their associated uncertainties. Building on previous work, this paper contains a study of the main complications associated with estimating and interpreting shear power spectra, and presents fast and accurate methods to estimate two key quantities needed for their practical usage: the noise bias and the Gaussian covariance matrix, fully accounting for survey geometry, with some of these results also applicable to other cosmological probes. Here, we demonstrate the performance of these methods by applying them to the latest public data releases of the Hyper Suprime-Cam and the Dark Energy Survey collaborations, quantifying the presence of systematics in our measurements and the validity of the covariance matrix estimate. We make the resulting power spectra, covariance matrices, null tests and all associated data necessary for a full cosmological analysis publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-F-19A-10: Recalculation and Reporting Clinical HDR 192-Ir Head and Neck Dose Distributions Using Model Based Dose Calculation

Purpose: To retrospectively re-calculate dose distributions for selected head and neck cancer patients, earlier treated with HDR 192Ir brachytherapy, using Monte Carlo (MC) simulations and compare results to distributions from the planning system derived using TG43 formalism. To study differences between dose to medium (as obtained with the MC code) and dose to water in medium as obtained through (1) ratios of stopping powers and (2) ratios of mass energy absorption coefficients between water and medium. Methods: The MC code Algebra was used to calculate dose distributions according to earlier actual treatment plans using anonymized plan data and CT images in DICOM format. Ratios of stopping power and mass energy absorption coefficients for water with various media obtained from 192-Ir spectra were used in toggling between dose to water and dose to media. Results: Differences between initial planned TG43 dose distributions and the doses to media calculated by MC are insignificant in the target volume. Differences are moderate (within 4–5 % at distances of 3–4 cm) but increase with distance and are most notable in bone and at the patient surface. Differences between dose to water and dose to medium are within 1-2% when using mass energy absorption coefficients to toggle between the two quantities but increase to above 10% for bone using stopping power ratios. Conclusion: MC predicts target doses for head and neck cancer patients in close agreement with TG43. MC yields improved dose estimations outside the target where a larger fraction of dose is from scattered photons. It is important with awareness and a clear reporting of absorbed dose values in using model based algorithms. Differences in bone media can exceed 10% depending on how dose to water in medium is defined.

60 APPLIED LIFE SCIENCES↗