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 1,495 records · Page 83

Intershell correlations in nonresonant Compton scattering of an X-ray photon by an atom

The role of intershell correlations in nonresonant Compton scattering of an X-ray photon by a free multielectron atom is studied theoretically for the Ar atom. The results of calculation are of a predictive nature. The developed mathematical formalism is general in nature and can be applied to a wide set of elements from the Periodic Table, for which the description of the wavefunctions of scattering states in the nonrelativistic Hartree-Fock approximation remains correct.

74 ATOMIC AND MOLECULAR PHYSICS↗

Dose Uncertainty Quantification in AGC Experiment

This memorandum formally documents completion of the Advanced Reactor Technologies Level 3 Milestone (M3NT-18IN060504037), “Interoffice Memorandum detailing status of the AGC irradiation uncertainty analysis” due SeptemberAugust 30, 2018. This memorandum details the progress made toward the issuance of a final report detailing the AGC uncertainty analysis including the proposed methodology, progress on activities, and explanation of how the data will be used. Radiation damage estimation is an important component of the post-irradiation analysis of the AGC experiments. It depends primarily on the fast fluence, which is determined using well-established methods of spectral adjustments. These are based on best estimates from models (e.g. MCNP). Input cross-sections, and measured activities from flux wires in the experiment, propagate uncertainties in each of these. While the methods of propagating uncertainty are well-established, the final uncertainty estimates they provide are only as good as the estimates of uncertainty in the inputs on which they are based. The purpose of this work is to implement a method that eliminates some deficiencies in the ways input uncertainties are presently estimated. The rRadiation damage estimates for the AGC experiments are ostensibly based on experimental data in the form of flux wire activation measurements; each AGC experiment has been instrumented with flux wires at various locations along the axial length to provide this data. However, the estimatesy also depend on the results of MCNP calculations. To understand why, consider that the objective of the analysis is to estimate some number of parameters (the energy group fluences ?_i) that best reproduce the observations (the measured activities of the flux wires a_i).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High-Level Model Articulation with BuildingSync and OpenStudio

The use of the BuildingSync schema for describing the contents of buildings is becoming more common with its recent integration into the Audit Template tool as well as ASHRAE’s Building Energy Quotient (bEQ) web portal. Although BuildingSync was initially created to store and transfer data related to building energy audits (as defined by ASHRAE Standard 211), it has since been expanded to store the data needed to articulate fully defined physics-based building energy models. BuildingSync combined with abstracted high-level input methods defined in OpenStudio’s Standards project and the newly developed BuildingSync gem allows BuildingSync eXtensible Markup Language (XML) (Bray, Paoli, Sperberg-McQueen, Maler, Eve (Sun Microsystems, & Francois, 2008) document to be converted into OpenStudio models. Automatic generation of a building energy model from data collected during an audit 1) eliminates the need to separately generate an energy model, 2) improves consistency between model generation, and 3) simplifies evaluation of various energy efficiency measures. This manuscript will discuss this open-source project, including: development of the validation infrastructure necessary to provide formalized expectations of informational requirements for BuildingSync documents; development of the BuildingSync gem for translation of BuildingSync documents to OpenStudio models as well as example implementations; and finally, the manuscript will elaborate on the advantages and disadvantages of using high-level models generated from BuildingSync as surrogates to detailed models.

30 DIRECT ENERGY CONVERSION↗

ON GALACTIC DENSITY MODELING IN THE PRESENCE OF DUST EXTINCTION

Inferences about the spatial density or phase-space structure of stellar populations in the Milky Way require a precise determination of the effective survey volume. The volume observed by surveys such as Gaia or near-infrared spectroscopic surveys, which have good coverage of the Galactic midplane region, is highly complex because of the abundant small-scale structure in the three-dimensional interstellar dust extinction. We introduce a novel framework for analyzing the importance of small-scale structure in the extinction. This formalism demonstrates that the spatially complex effect of extinction on the selection function of a pencil-beam or contiguous sky survey is equivalent to a low-pass filtering of the extinction-affected selection function with the smooth density field. We find that the angular resolution of current 3D extinction maps is sufficient for analyzing Gaia sub-samples of millions of stars. However, the current distance resolution is inadequate and needs to be improved by an order of magnitude, especially in the inner Galaxy. We also present a practical and efficient method for properly taking the effect of extinction into account in analyses of Galactic structure through an effective selection function. We illustrate its use with the selection function of red-clump stars in APOGEE using and comparing a variety of current 3D extinction maps.

79 ASTRONOMY AND ASTROPHYSICS↗

Two-component hybrid time-dependent density functional theory within the Tamm-Dancoff approximation

We report the implementation of a two-component variant of time-dependent density functional theory (TDDFT) for hybrid functionals that accounts for spin-orbit effects within the Tamm-Dancoff approximation (TDA) for closed-shell systems. The influence of the admixture of Hartree-Fock exchange on excitation energies is investigated for several atoms and diatomic molecules by comparison to numbers for pure density functionals obtained previously [M. Kühn and F. Weigend, J. Chem. Theory Comput. 9, 5341 (2013)]. It is further related to changes upon switching to the local density approximation or using the full TDDFT formalism instead of TDA. Efficiency is demonstrated for a comparably large system, Ir(ppy){sub 3} (61 atoms, 1501 basis functions, lowest 10 excited states), which is a prototype molecule for organic light-emitting diodes, due to its “spin-forbidden” triplet-singlet transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Foundations of variational discrete action theory

Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

SU-E-P-08: Alarming Range of Seed Activities Ordered for I-125 Plaque Brachytherapy

Purpose: To investigate the variation in I-125 seed activities ordered by various clinics for their plaque brachytherapy cases under a standardized set of assumptions. Methods: A majority of the plaque programs in North America were contacted and a survey was designed to give a few standardized cases to allow inter-comparison of seed activities ordered. Tumor dose, treatment duration, number of seeds, plaque, and tumor apex were held constant in order to reveal differences in prescription point, seed type, and seed activity. Results: While the survey is presently underway, preliminary results show alarmingly wide variations between centers. Differences up to 45% have been found with 15% differences being common. Conclusion: Though knowledge of the TG-43 dose calculation formalism is common, a number of factors in the field of plaque brachytherapy lead to alarming differences in activity of I-125 seeds being ordered for a given tumor. Knowledge of the present reality of widely varying treatment activities, and thus doses to tumor and normal structures, should serve as motivation for centers involved in this modality to review their programs with others in the community and share their experiences.

60 APPLIED LIFE SCIENCES↗

Investigation of scattering processes in quantum few-body systems involving long-range interaction by the complex-rotation method

The complex-rotation method adapted to solving the multichannel scattering problem in the two-body system where the interaction potential contains the long-range Coulomb components is described. The scattering problem is reformulated as the problem of solving a nonhomogeneous Schroedinger equation in which the nonhomogeneous term involves a Coulomb potential cut off at large distances. The incident wave appearing in the nonhomogeneous term is a solution of the Schroedinger equation with longrange Coulomb interaction. This formulation is free from approximations associated with a direct cutoff of Coulomb interaction at large distances. The efficiency of this formalism is demonstrated by considering the example of solving scattering problems in the {alpha}-{alpha} and p-p systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Perioperative Gemcitabine + Erlotinib Plus Pancreaticoduodenectomy for Resectable Pancreatic Adenocarcinoma: ACOSOG Z5041 (Alliance) Phase II Trial

Background: There is considerable interest in a neoadjuvant approach for resectable pancreatic ductal adenocarcinoma (PDAC). This study evaluated perioperative gemcitabine + erlotinib (G+E) for resectable PDAC. Methods: A multicenter, cooperative group, single-arm, phase II trial was conducted between April 2009 and November 2013 (ACOSOG Z5041). Patients with biopsy-confirmed PDAC in the pancreatic head without evidence of involvement of major mesenteric vessels (resectable) were eligible. Patients (n = 123) received an 8-week cycle of G+E before and after surgery. The primary endpoint was 2-year overall survival (OS), and secondary endpoints included toxicity, response, resection rate, and time to progression. Resectability was assessed retrospectively by central review. The study closed early due to slow accrual, and no formal hypothesis testing was performed. Results: Overall, 114 patients were eligible, consented, and initiated protocol treatment. By central radiologic review, 97 (85%) of the 114 patients met the protocol-defined resectability criteria. Grade 3+ toxicity was reported in 60% and 79% of patients during the neoadjuvant phase and overall, respectively. Twenty-two of 114 (19%) patients did not proceed to surgery; 83 patients (73%) were successfully resected. R0 and R1 margins were obtained in 67 (81%) and 16 (19%) resected patients, respectively, and 54 patients completed postoperative G+E (65%). The 2-year OS rate for the entire cohort (n = 114) was 40% (95% confidence interval [CI] 31–50), with a median OS of 21.3 months (95% CI 17.2–25.9). The 2-year OS rate for resected patients (n = 83) was 52% (95% CI 41–63), with a median OS of 25.4 months (95% CI 21.8–29.6). Conclusions: For resectable PDAC, perioperative G+E is feasible. Further evaluation of neoadjuvant strategies in resectable PDAC is warranted with more active systemic regimens.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Water Splitting by Thin Film Metal-Oxo Catalysts

The dropping price of silicon photovoltaics in the United States is causing load defection to solar supply at an accelerated pace. This conversion to solar and, more generally, other renewable energy sources has accordingly turned the energy research focus from generation to one of storage. Truly disruptive improvements in energy storage technologies are limited by energy density. This limitation, however, does not apply to fuels, which possess the energy density needed for large-scale energy storage. The first step of the basic science needed to drive such historic restructuring of the U.S. energy infrastructure begins with the solar-driven generation of hydrogen and oxygen from water. The solar-produced hydrogen may then be combined with carbon dioxide to deliver any number of fuels. Obviously, light does not directly act on water to engender its splitting into its elemental components. Hence, catalysts are needed to drive the oxygen evolution reaction (OER) and hydrogen evolution reaction (HER). Of these two reactions, the four-electron, four-proton oxidation of OER is the more kinetically challenging reaction, and therefore the development of energy efficient solar fuels processes demands that OER be accomplished at a minimal overpotential. The research completed in this program developed catalysts that drive OER and at the same time meet the important criteria of (1) using non-critical materials that (2) are easy to assemble and (3) accomplish OER under simple conditions. Research was designed to uncover the chemical principles that underlie the self-assembly of metal oxide oxygen evolving catalysts (M-OEC) from the metals of M = Mn, Co, and Ni. For example, a dogma of heterogeneous catalysis of any sort is that “edges” matter in promoting catalytic transformations. We provided a rationale for such dogma by showing that the OER in Co-OEC occurred at a dimensionally reduced dicobalt edge site. Edge site reactivity was clearly revealed analyzing 18O labeled OER with differential electrochemical mass spectrometry (DEMS) of Co-OEC. The OER mechanism of M-OECs was examined with complementary studies of model dicobalt compounds that captured the critical steps of the OER reaction. Additionally, the role of activating M-OECS with metal ion dopants was defined by developing structure–function relationships, guided by the principles of inorganic chemistry. We found that the M(IV) oxidation state in oxidic OER frameworks was correlated to the presence of the dopant metal, as assessed by coulometric titration and ICP-MS analysis. To investigate why greater M(IV) valence is beneficial to greater catalytic OER activity, we probed the influence of formal M valence on the electronic structure of oxygen ions in M-OECs by undertaking O and Ni K-edge spectroscopy, which revealed greater M-O covalency and hence M-oxyl radical character with M(IV) formation. Such oxyl radical character is consistent with increasing evidence for the role of oxygen radicals in O–O bond formation by a proton-coupled electron transfer mechanism involving water, to generate a hydroperoxide intermediate from which oxygen is generate. In accomplishing this science, the DOE program leveraged its expertise in spectroscopy and structural methods, inorganic and materials synthesis, and electrochemical characterization. The knowledge garnered from this proposed program enables the design of next generation catalysts with improved OER kinetics that operate over a wide range of conditions and environments.

25 ENERGY STORAGE↗

Applying the Gyrokinetic Formulation of Magnetized Plasma Dynamics to Z-pinch plasmas

This project adapts and applies the gyrokinetic Eulerian code COGENT to an emerging area of highly-magnetized Z-pinch plasma systems. Simulation efforts for these plasmas currently use either magneto-hydrodynamic (MHD) fluid codes, which suffer from fidelity being restricted to narrow plasma conditions, or fully kinetic particle-in-cell (PIC) methods, which suffer from the run-time needed to resolve the very small spatial and temporal scales required for the plasma conditions. The progress with achieving detailed theoretical understanding of Z-pinch plasma systems can be greatly expedited using some combination of gyrokinetic and drift-ordered MHD models, developed in the tokamak community and implemented in the COGENT code. These formulations are applicable to a much larger parameter space than standard MHD models, while providing a substantial speed-up in run-time over fully kinetic codes. In particular, the gyrokinetic formalism is successfully applied to a flow-stabilized Z pinch (FSZP) system – a promising fusion concept where a sheared axial flow is used to provide stability against large-scale ideal-MHD modes. Results of gyrokinetic simulations with COGENT agree well with those from fully kinetic simulations, while requiring ~500x less run-time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optical Constants Determination of Samarium Holmium and Erbium in the 1.5-850 eV Spectral Range using a Transmittance Method

The optical constants {beta} and {sigma} of the complex refractive index {tilde n} = 1 - {delta} + i{beta} of Sm, Ho, and Er were obtained in the 1.5-850 eV energy range using a transmittance method. Thin films of Sm, Ho, and Dr were deposited by magnetron sputtering, and transmittance was measured using synchrotron radiation using a high vacuum condition. All films were directly coated on Si photodiodes, which were used as coating substrates, as well as photon detectors. Si was used as capping layer while a thin W layer was used as barrier against interface diffusion between Si and the highly reactive rare earth elements. The constants {beta} were extracted from transmittance results, and the constants {sigma} were calculated based on measured {beta} values using the Kramers-Kronig formalism. Small deficiencies determined from the present data using the partial sum rules were partially attributed to the sputtered film densities that could be slightly lower than the bulk values.

36 MATERIALS SCIENCE↗

A Comparative Study of the Efficacy of Local/Global and Parametric/Nonparametric Machine Learning Methods for Establishing Structure–Property Linkages in High-Contrast 3D Elastic Composites

Reduced-order structure–property (S-P) linkages play a pivotal role in the tailored design of materials for advanced engineering components. There is a critical need to distill these from the simulation datasets aggregated using sophisticated, computationally expensive, physics-based numerical tools (e.g., finite element methods). The recent emergence of materials data science approaches has opened new avenues for addressing this challenge. In this paper, we critically compare the relative merits of the application of four distinct machine learning approaches for their efficacy in extracting microstructure-property linkages from the finite element simulation data aggregated on high-contrast elastic composites with different microstructures. The machine learning approaches selected for the study have included different combinations of local/global and parametric/nonparametric approaches. Furthermore, the nonparametric approaches selected for this study are based on Gaussian Process (GP) models that allow for a formal treatment of uncertainty quantification in the predicted values. The predictive performances of these different approaches have been compared against each other using rigorous cross-validation error metrics. Furthermore, their sensitivity to both the dataset size and dimensionality has been investigated.

36 MATERIALS SCIENCE↗

Structure, stability, and dissociation of small ionic silicon oxide clusters [SiO{sub n}{sup +}(n = 3, 4)]: Insight from density functional and topological exploration

The structures, energies, isomerization, and decomposition pathways of small ionic silicon oxide clusters, SiO{sub n}{sup +} (n = 3, 4), on doublet and quartet energy surfaces are investigated by density functional theory. New structural isomers of these ionic clusters have been obtained with this systematic study. The energy ordering of the isomeric cluster ions on doublet spin surface is found to follow the same general trend as that of the neutral ones, while it differs on the quartet surface. Our computational results reveal the energetically most preferred decomposition pathways of the ionic clusters on both spin surfaces. To comprehend the reaction mechanism, bonding evolution theory has also been employed using atoms in molecules formalism. The possible reasons behind the structural deformation of some isomers on quartet surface have also been addressed. Our results are expected to provide important insight into the decomposition mechanism and relative stability of the SiO{sub n}{sup +} clusters on both the energy surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mutual diffusion of binary liquid mixtures containing methanol, ethanol, acetone, benzene, cyclohexane, toluene, and carbon tetrachloride

Mutual diffusion coefficients of all 20 binary liquid mixtures that can be formed out of methanol, ethanol, acetone, benzene, cyclohexane, toluene, and carbon tetrachloride without a miscibility gap are studied at ambient conditions of temperature and pressure in the entire composition range. The considered mixtures show a varying mixing behavior from almost ideal to strongly non-ideal. Predictive molecular dynamics simulations employing the Green-Kubo formalism are carried out. Radial distribution functions are analyzed to gain an understanding of the liquid structure influencing the diffusion processes. It is shown that cluster formation in mixtures containing one alcoholic component has a significant impact on the diffusion process. The estimation of the thermodynamic factor from experimental vapor-liquid equilibrium data is investigated, considering three excess Gibbs energy models, i.e., Wilson, NRTL, and UNIQUAC. It is found that the Wilson model yields the thermodynamic factor that best suits the simulation results for the prediction of the Fick diffusion coefficient. Four semi-empirical methods for the prediction of the self-diffusion coefficients and nine predictive equations for the Fick diffusion coefficient are assessed and it is found that methods based on local composition models are more reliable. Finally, the shear viscosity and thermal conductivity are predicted and in most cases favorably compared with experimental literature values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stabilization of a locally minimal forest

The method of partial stabilization of locally minimal networks, which was invented by Ivanov and Tuzhilin to construct examples of shortest trees with given topology, is developed. According to this method, boundary vertices of degree 2 are not added to all edges of the original locally minimal tree, but only to some of them. The problem of partial stabilization of locally minimal trees in a finite-dimensional Euclidean space is solved completely in the paper, that is, without any restrictions imposed on the number of edges remaining free of subdivision. A criterion for the realizability of such stabilization is established. In addition, the general problem of searching for the shortest forest connecting a finite family of boundary compact sets in an arbitrary metric space is formalized; it is shown that such forests exist for any family of compact sets if and only if for any finite subset of the ambient space there exists a shortest tree connecting it. The theory developed here allows us to establish further generalizations of the stabilization theorem both for arbitrary metric spaces and for metric spaces with some special properties. Bibliography: 10 titles.

97 MATHEMATICS AND COMPUTING↗

Automated Vulnerability Detection for Compiled Smart Grid Software

While testing performed with proper experimental controls can provide scientifically quantifiable evidence that software does not contain unintentional vulnerabilities (bugs), it is insufficient to show that intentional vulnerabilities exist, and impractical to certify devices for the expected long lifetimes of use. For both of these needs, rigorous analysis of the software itself is essential. Automated software behavior computation applies rigorous static software analysis methods based on function extraction (FX) to compiled software to detect vulnerabilities, intentional or unintentional, and to verify critical functionality. This analysis is based on the compiled firmware, takes into account machine precision, and does not rely on heuristics or approximations early in the analysis.

99 GENERAL AND MISCELLANEOUS↗

Approximate but accurate quantum dynamics from the Mori formalism: I. Nonequilibrium dynamics

We present a formalism that explicitly unifies the commonly used Nakajima-Zwanzig approach for reduced density matrix dynamics with the more versatile Mori theory in the context of nonequilibrium dynamics. Employing a Dyson-type expansion to circumvent the difficulty of projected dynamics, we obtain a self-consistent equation for the memory kernel which requires only knowledge of normally evolved auxiliary kernels. To illustrate the properties of the current approach, we focus on the spin-boson model and limit our attention to the use of a simple and inexpensive quasi-classical dynamics, given by the Ehrenfest method, for the calculation of the auxiliary kernels. For the first time, we provide a detailed analysis of the dependence of the properties of the memory kernels obtained via different projection operators, namely, the thermal (Redfield-type) and population based (NIBA-type) projection operators. We further elucidate the conditions that lead to short-lived memory kernels and the regions of parameter space to which this program is best suited. Via a thorough analysis of the different closures available for the auxiliary kernels and the convergence properties of the self-consistently extracted memory kernel, we identify the mechanisms whereby the current approach leads to a significant improvement over the direct usage of standard semi- and quasi-classical dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗