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.

SourceResearch connectionFrom yearThrough year

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

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

DOE OSTI2013

Full dimensional potential energy surface for the ground state of H{sub 4}{sup +} system based on triatomic-in-molecules formalism

In this work, we present a global potential energy surface for the ground electronic state of the H{sub 4}{sup +} based on ab initio calculations. The final fit is based on triatomics-in-molecules (TRIM) approximation and it includes extra four-body terms for the better description of some discrepancies found on the TRIM model. The TRIM method itself allows a very accurate description of the asymptotic regions. The global fit uses more than 19 000 multireference configuration interaction ab initio points. The global potential energy surface has an overall root mean square error of 0.013 eV for energies up to 2 eV above the global minimum. This work presents an analysis of the stationary points, reactant and product channels, and crossing between the two lowest TRIM adiabatic states. It is as well included a brief description of the two first excited states of the TRIM matrix, concluding that TRIM method is a very good approximation not only for the ground state but also for at least two of the excited states of H{sub 4}{sup +} system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Application of vibrational correlation formalism to internal conversion rate: Case study of Cu{sub n} (n = 3, 6, and 9) and H{sub 2}/Cu{sub 3}

This work reports non-radiative internal conversion (IC) rate constants obtained for Cu{sub n} with n = 3, 6, and 9 and H{sub 2} on Cu{sub 3}. The Time-Dependent Density Functional Theory (TDDFT) method was employed with three different functionals in order to investigate the electronic structures and the absorption spectra. The performance of the generalized gradient approximation of Perdew, Burke and Ernzerhof (PBE) and the hybrid B3LYP and PBE0 exchange correlation functionals in combination with the SVP and the def2-TZVP basis sets was examined. TDDFT results were used as input data to compute internal conversion rate constants. For this purpose, we have developed a program package. A description of the theoretical background used in our numerical implementation and the program input file is presented. In view of future applications of this program package in photoinduced catalysis, we present the analysis of the IC rate processes for the photodissociation of H{sub 2} on Cu{sub 3}. These results showed the applicability of the method and the computational program to identify the vibrational modes in transition metal clusters giving rise to the largest IC rate constant due to their interactions with the excited electronic states occurring in the hot-electron induced dissociation phenomena.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2016

SU-F-J-195: On the Performance of Four Dual Energy CT Formalisms for Extracting Proton Stopping Powers

Purpose: Dual energy CT can predict stopping power ratios (SPR) for ion therapy treatment planning. Several approaches have been proposed recently, however accuracy and practicability in a clinical workflow are unaddressed. The aim of this work is to provide a fair comparison of available approaches in a human-like phantom to find the optimal method for tissue characterization in a clinical situation. Methods: The SPR determination accuracy is investigated using simulated DECT images. A virtual human-like phantom is created containing 14 different standard human tissues. SECT (120 kV) and DECT images (100 kV and 140 kV Sn) are simulated using the software ImaSim. The single energy CT (SECT) stoichiometric calibration method and four recently published calibration-based DECT methods are implemented and used to predict the SPRs from simulated images. The difference between SPR predictions and theoretical SPR are compared pixelwize. Mean, standard deviation and skewness of the SPR difference distributions are used as measures for bias, dispersion and symmetry. Results: The average SPR differences and standard deviations are (0.22 ± 1.27)% for SECT, and A) (−0.26 ± 1.30)%, B) (0.08 ± 1.12)%, C) (0.06 ± 1.15)% and D) (−0.05 ± 1.05)% for the four DECT methods. While SPR prediction using SECT is showing a systematic error on SPR, the DECT methods B, C and D are unbiased. The skewness of the SECT distribution is 0.57%, and A) −0.19%, B) −0.56%, C) −0.29% and D) −0.07% for DECT methods respectively. Conclusion: The here presented DECT methods B, C and D outperform the commonly used SECT stoichiometric calibration. These methods predict SPR accurately without a bias and within ± 1.2% (68th percentile). This indicates that DECT potentially improves accuracy of range predictions in proton therapy. A validation of these findings using clinical CT images of real tissues is necessary.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

Investigation of surface termination by directional elastic peak electron spectroscopy: Experiment and theory

Experimental and theoretical details concerning the directional elastic peak electron spectroscopy are presented. The application of this experimental method to the identification of a crystalline structure of surface layers is shown for Cu(111) and Ru(0001), which enables the analysis of different surface terminations associated with different sequences of atoms along the surface normal. Theoretical data are obtained by applying the multiple scattering theory and different sources of phase shifts. The quantitative analysis performed by an R-factor reveals almost the same populations of A and B terminated (0001) terraces and proves explicitly the presence of steps of atomic height on the clean Ru(0001) surface. This is not the case for the Cu(111) surface, where the terrace termination does not affect the distribution of atomic directions within the first few atomic layers. The averaging of theoretical data versus scattering volume defined by the sphere radius R{sub max} around the emitter site is discussed in view of the computation time optimization. - Highlights: • We examined the structure and surface termination of fcc and hcp monocrystals. • Crystalline structure within the first few atomic layers was determined by DEPES. • Multiple scattering formalism was used to obtain theoretical results. • Populations of terraces with different terminations were determined. • Scattering parameters are discussed in the context of computation time optimization.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Nonlinear dynamics from linear quantum evolutions

Linear dynamics restricted to invariant submanifolds generally gives rise to nonlinear dynamics. Submanifolds in the quantum framework may emerge for several reasons: one could be interested in specific properties possessed by a given family of states, either as a consequence of experimental constraints or inside an approximation scheme. In this work we investigate such issues in connection with a one parameter group ϕ{sub t} of transformations on a Hilbert space, H, defining the unitary evolutions of a chosen quantum system. Two procedures will be presented: the first one consists in the restriction of the vector field associated with the Schrödinger equation to a submanifold invariant under the flow ϕ{sub t}. The second one makes use of the Lagrangian formalism and can be extended also to non-invariant submanifolds, even if in such a case the resulting dynamics is only an approximation of the flow ϕ{sub t}. Such a result, therefore, should be conceived as a generalization of the variational method already employed for stationary problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

What can be learned from optical two-color diffusion and thermodiffusion experiments on ternary fluid mixtures?

A number of optical techniques have been developed during the recent years for the investigation of diffusion and thermodiffusion in ternary fluid mixtures, both on ground and on-board the International Space Station. All these methods are based on the simultaneous measurement of refractive index changes at two different wavelengths. Here, we discuss and compare different techniques with the emphasis on optical beam deflection (OBD), optical digital interferometry, and thermal diffusion forced Rayleigh scattering (TDFRS). We suggest to formally split the data evaluation into a phenomenological parameterization of the measured transients and a subsequent transformation from the refractive index into the concentration space. In all experiments, the transients measured at two different detection wavelengths can be described by four amplitudes and two eigenvalues of the diffusion coefficient matrix. It turns out that these six parameters are subjected to large errors and cannot be determined reliably. Five good quantities, which can be determined with a high accuracy, are the stationary amplitudes, the initial slopes as defined in TDFRS experiments and by application of a heuristic criterion for similar curves, a certain mean diffusion coefficient. These amplitudes and slopes are directly linked to the Soret and thermodiffusion coefficients after transformation with the inverse contrast factor matrix, which is frequently ill-conditioned. Since only five out of six free parameters are reliably determined, including the single mean diffusion coefficient, the determination of the four entries of the diffusion matrix is not possible. We apply our results to new OBD measurements of the symmetric (mass fractions 0.33/0.33/0.33) ternary benchmark mixture n-dodecane/isobutylbenzene/1,2,3,4-tetrahydronaphthalene and existing literature data for the same system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Unstable states in quantum theory

Various approaches to the problem of describing unstable particles are reviewed. Fundamental problems that arise in quantum field description of these particles and the ways of their solution are considered. Among them, there is an approach related to the notion of the smeared (continuous) mass, which originates from the finite lifetime of unstable particles. The quantum field model of unstable particles with smeared mass, which is built upon two basic axiomatic elements, is considered in detail. The basic processes with unstable particles (decay and scattering) are considered within the framework of the model and the formalism for describing physical characteristics of those processes is developed. The model is successfully applied to describing the processes of pair and triple boson production at the linear collider, top quark pair production, and certain hadronic decays. Based on this model, the factorization method is developed, which allows a description of complicated and multistep scattering and decay processes with unstable particles to be considerably simplified.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Certifying the Performance of Fixative Technologies under Open Air Demolition Activities for D and D - 20096

Department of Energy (DOE) facilities undergoing deactivation and decommissioning (D and D) activities struggle with safety concerns of residual radioactive contamination leftover after gross decontamination efforts have concluded. The positive effects of implementing fixative technologies for decommissioning and maintenance is known across the DOE complex. However, there is no standardized metric in place that quantifies the performance of these technologies such as incombustible fixative platforms that immobilize the contamination in a solid polymer material under normal operating conditions or when exposed to stressors. The primary focus of this research is to provide the empirical data necessary to properly characterize the effects of implementing fixative technologies on mitigating the release of contamination. In addition, this effort has been used to identify and establish uniform testing protocols to quantify operational parameters of fixating platforms and properly credit such systems under a variety of operational conditions and stressors that can arise during open air demolition efforts (e.g. impact, water, thermal, etc.). The established protocols will be provided to the ASTM International E10.03 Subcommittee on Radiological Protection for Decontamination and Decommissioning of Nuclear Facilities and Components for the potential of formal standardization. Attributes to be addressed during testing are material compromise from impact stress, thermal stress for incombustible fixative materials, and immobilization factors (how well an immobilizing platform is capable of retaining fixated material during thermal and impact stressors). In developing these test methods, a greater understanding of the material's behavior under anticipated and unanticipated events, such as decommissioning activities and contingency events as outlined in the Basis of Interim Operations (BIO), can be successfully characterized. There are several key components of the experimental methodology that is essential in empirically certifying fixative technologies. A surrogate contaminant with a unique signature has been utilized with the implementation of a uniform contamination process of test coupons that is quantifiable and replicable for direct fixative comparisons. A modular test chamber has also been used to apply the various stressors of interest. Collection of released contamination that includes airborne and resettled contaminant particles and analyzed using mass spectrometry. This research directly supports DOE complex-wide concerns in terms of final disposition of nuclear facilities and the considerations made during their time between ceasing facility operation and final disposition. The results of this effort aims to provide a standardized road map and logical decision flow to better facilitate fixative technologies into D and D activities required to achieve the desired facility end state. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2014

An asymptotic-preserving Lagrangian algorithm for the time-dependent anisotropic heat transport equation

We propose a Lagrangian numerical algorithm for a time-dependent, anisotropic temperature transport equation in magnetized plasmas in the large guide field regime. The approach is based on an analytical integral formal solution of the parallel (i.e., along the magnetic field) transport equation with sources, and it is able to accommodate both local and non-local parallel heat flux closures. The numerical implementation is based on an operator-split formulation, with two straightforward steps: a perpendicular transport step (including sources), and a Lagrangian (field-line integral) parallel transport step. Algorithmically, the first step is amenable to the use of modern iterative methods, while the second step has a fixed cost per degree of freedom (and is therefore scalable). Accuracy-wise, the approach is free from the numerical pollution introduced by the discrete parallel transport term when the perpendicular to parallel transport coefficient ratio X ⊥ /X ∥ becomes arbitrarily small, and is shown to capture the correct limiting solution when ε = X⊥L 2 ∥ /X1L 2 ⊥ → 0 (with L∥∙ L⊥ , the parallel and perpendicular diffusion length scales, respectively). Therefore, the approach is asymptotic-preserving. We demonstrate the capabilities of the scheme with several numerical experiments with varying magnetic field complexity in two dimensions, including the case of transport across a magnetic island.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

A new dipolar potential for numerical simulations of polar fluids on the 4D hypersphere

We present a new method for Monte Carlo or Molecular Dynamics numerical simulations of three-dimensional polar fluids. The simulation cell is defined to be the surface of the northern hemisphere of a four-dimensional (hyper)sphere. The point dipoles are constrained to remain tangent to the sphere and their interactions are derived from the basic laws of electrostatics in this geometry. The dipole-dipole potential has two singularities which correspond to the following boundary conditions: when a dipole leaves the northern hemisphere at some point of the equator, it reappears at the antipodal point bearing the same dipole moment. We derive all the formal expressions needed to obtain the thermodynamic and structural properties of a polar liquid at thermal equilibrium in actual numerical simulation. We notably establish the expression of the static dielectric constant of the fluid as well as the behavior of the pair correlation at large distances. We report and discuss the results of extensive numerical Monte Carlo simulations for two reference states of a fluid of dipolar hard spheres and compare these results with previous methods with a special emphasis on finite size effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Real time quantum dynamics preaveraged over imaginary time path integral: A formal basis for both Centroid Molecular Dynamics and Ring Polymer Molecular Dynamics

An exact real time quantum dynamics preaveraged over imaginary time path integral is formulated for general condensed phase equilibrium ensemble. This formulation results in the well-known centroid dynamics approach upon filtering of centroid constraint, and provides a rigorous framework to understand and analyze a related quantum dynamics approximation method called ring polymer molecular dynamics. The formulation also serves as the basis for developing new kinds of quantum dynamics that utilize the cyclic nature of the imaginary time path integral.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Effective long wavelength scalar dynamics in de Sitter

We discuss the effective infrared theory governing a light scalar's long wavelength dynamics in de Sitter spacetime. We show how the separation of scales around the physical curvature radius k / a ∼ H can be performed consistently with a window function and how short wavelengths can be integrated out in the Schwinger-Keldysh path integral formalism. At leading order, and for time scales Δ t >> H {sup −1}, this results in the well-known Starobinsky stochastic evolution. However, our approach allows for the computation of quantum UV corrections, generating an effective potential on which the stochastic dynamics takes place. The long wavelength stochastic dynamical equations are now second order in time, incorporating temporal scales Δ t ∼ H {sup −1} and resulting in a Kramers equation for the probability distribution—more precisely the Wigner function—in contrast to the more usual Fokker-Planck equation. This feature allows us to non-perturbatively evaluate, within the stochastic formalism, not only expectation values of field correlators, but also the stress-energy tensor of φ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2015

Low-power, ultrafast, and dynamic all-optical tunable plasmonic analog to electromagnetically induced transparency in two resonators side-coupled with a waveguide system

We theoretically and numerically investigate a low-power, ultrafast, and dynamic all-optical tunable plasmonic analog to electromagnetically induced transparency (EIT) in two nanodisk resonators side-coupled to a metal-insulator-metal plasmonic waveguide system. The optical Kerr effect is enhanced by the slow light effect of the plasmonic EIT-like effect and the plasmonic waveguide based on graphene-Ag composite material structures with giant effective Kerr nonlinear coefficient. The optical Kerr effect modulation method is applied to improve tuning rate with response time of subpicoseconds or even femtoseconds. With dynamically tuning the propagation phase of the plasmonic waveguide, π-phase shift of the transmission spectrum in the plasmonic EIT-like system is achieved under excitation of a pump light with an intensity as low as 5.85 MW/cm{sup 2}. The group delay is controlled between 0.09 and 0.4 ps. All observed schemes are analyzed rigorously through finite-difference time-domain simulations and coupled-mode formalism. Results show a new direction toward the low power consumption and ultrafast responses of integration plasmonic photonic devices and all-optical dynamical storage of light devices in optical communication and quantum information processing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

The Epstein–Glaser causal approach to the light-front QED{sub 4}. I: Free theory

In this work we present the study of light-front field theories in the realm of the axiomatic theory. It is known that when one uses the light-cone gauge pathological poles (k{sup +}){sup −n} arises, demanding a prescription to be employed in order to tame these ill-defined poles and to have the correct Feynman integrals due to the lack of Wick rotation in such theories. In order to shed a new light on this long standing problem we present here a discussion based on the use of rigorous mathematical machinery of the distributional theory combined with physical concepts, such as causality, to show how to deal with these singular propagators in a general fashion without making use of any prescription. The first step of our development will consist in showing how the analytic representation for propagators arises by requiring general physical properties within the framework of Wightman’s formalism. From that we shall determine the equal-time (anti)commutation relations in the light-front form for the scalar and fermionic fields, as well as for the dynamical components of the electromagnetic field. In conclusion, we introduce the Epstein–Glaser causal method in order to have a mathematical rigorous description of the free propagators of the theory, allowing us to discuss a general treatment for propagators of the type (k{sup +}){sup −n}. Afterwards, we show that at given conditions our results reproduce known prescriptions in the literature. - Highlights: • We develop the analytic representation for propagators in Wightman’s framework. • We make use of the analytic representation to obtain equal-time (anti)commutation relations in the light-front. • We derive the free Feynman propagators for the light-front quantum electrodynamics in the Epstein–Glaser approach. • We determine a general expression for the propagator associated to the light-cone poles (k{sup +}){sup −n} in the causal approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

A Review of the Lawrence Livermore Nuclear Accident Dosimeter 1980s-present

A Nuclear Accident Dosimetry program is a federal requirement for all facilities that have the potential to have a criticality accident. Personnel Nuclear Accident Dosimeter (PNAD) theory and analytical procedures are driven by various scientific needs and interacting regulations. A brief history of the status of USA Department of Energy (DOE) nuclear accident dosimetry regulations, recommendations, and performance testing criteria are given. Then, the history of the Lawrence Livermore National Laboratory (LLNL) PNAD is explored, including changes in the physical dosimeter and adjustments of the analysis method through the last four decades. Finally, the performance of LLNL’s PNAD at criticality accident intercomparison training exercises since 2009 is explored. In general, reported neutron doses have been within or close to DOE-STD-1098 performance criteria while reported gamma doses have been outside of DOE-STD-1098 performance criteria. Reported total absorbed doses have varied in meeting ANSI/HPS N13.3 and ANSI/HPS N13.3 (R2019) performance criteria. Dosimetry staff retirement and turnover have left historical knowledge gaps, yet provided opportunities within the NAD program at LLNL. This review paper serves as an overview of the history and status of the NAD program. Brief technical, procedural and programmatic recommendations to improve LLNL’s NAD program are given. Technical recommendations include investigating orientation factors through modeling or empirical experimentation, investigating gamma dosimetry methods for high-dose scenarios, and exploring other dosimetric methods for simpler, quicker NAD analysis. Procedural recommendations include better documentation of conversion factor (activity-to-fluence and fluence-to-dose) derivations and spectrum uses, and updated analysis spreadsheets or simple Graphic User Interfaces for dose calculations. In conclusion, programmatic recommendations include formalized training for NAD analysts, and having multiple SMEs trained on the NAD program.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Growth phase estimation for abundant bacterial populations sampled longitudinally from human stool metagenomes

Longitudinal sampling of the stool has yielded important insights into the ecological dynamics of the human gut microbiome. However, human stool samples are available approximately once per day, while commensal population doubling times are likely on the order of minutes-to-hours. Despite this mismatch in timescales, much of the prior work on human gut microbiome time series modeling has assumed that day-to-day fluctuations in taxon abundances are related to population growth or death rates, which is likely not the case. Here, we propose an alternative model of the human gut as a stationary system, where population dynamics occur internally and the bacterial population sizes measured in a bolus of stool represent a steady-state endpoint of these dynamics. We formalize this idea as stochastic logistic growth. We show how this model provides a path toward estimating the growth phases of gut bacterial populations in situ. We validate our model predictions using an in vitro Escherichia coli growth experiment. Finally, we show how this method can be applied to densely-sampled human stool metagenomic time series data. We discuss how these growth phase estimates may be used to better inform metabolic modeling in flow-through ecosystems, like animal guts or industrial bioreactors.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2014

Phase transformation of ZnMoO{sub 4} by localized thermal spike

We show that ZnMoO{sub 4} remains in stable phase under thermal annealing up to 1000 °C, whereas it decomposes to ZnO and MoO{sub 3} under transient thermal spike induced by 100 MeV Ag irradiation. The transformation is evidenced by X-ray diffraction (XRD), Raman spectroscopy, and X-ray photoelectron spectroscopy (XPS). Thin films of ZnMoO{sub 4} were synthesized by thermal evaporation and subsequent annealing in oxygen ambient at 600 °C for 4 h. XRD results show that as the irradiation fluence increases, the peak related to ZnMoO{sub 4} decreases gradually and eventually disappear, whereas peaks related to ZnO grow steadily up to fluence of 3 × 10{sup 12} ions/cm{sup 2} and thereafter remain stable till highest fluence. This indicates that polycrystalline ZnMoO{sub 4} film has transformed to polycrystalline ZnO thin film. The Raman lines related to ZnMoO{sub 4} are observed to have disappeared with increasing irradiation fluence. XPS results show modification in bonding and depletion of Mo from near surface region after the ion irradiation. Cross-sectional transmission electron microscopy result shows the formation of ion track of diameter 12–16 nm. These results demonstrate that ion beam methods provide the means to control phase splitting of ZnMoO{sub 4} to ZnO and MoO{sub 3} within nanometric dimension along the ion track. The observation of phase splitting and Mo loss are explained in the framework of ion beam induced thermal spike formalism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Parametrizations of coupled betatron motion for strongly coupled lattices

The coupling of transverse motion is a natural occurrence in particle accelerators, either in the form of a residual coupling arising from imperfections or originating by design from strong systematic coupling fields. While the first can be treated perturbatively, the latter requires a robust approach adapted to strongly coupled optics, and a parametrization of the linear optics must be performed to explore beam dynamics in such peculiar lattices. This work highlights the key physical interpretations of the main parametrization formalisms to describe linear coupled optics, along with explicit links and comparisons of these parametrizations. Concepts rarely illustrated in other works, such as forced mode flips and local coupling, are explored in detail, clarifying some anomalies that can arise in lattice functions. The analytical methods have been implemented in a reference Python package and connected with ray-tracing and integration codes to explore examples of strongly coupled lattices, which are discussed in detail to highlight the key physical interpretations of the parametrizations and characteristics of the lattices.

47 OTHER INSTRUMENTATION↗