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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

The Density Matrix of H20 - N2 In the Coordinate Representation: A Monte Carlo Calculation of the Far-Wing Line Shape

The far-wing line shape theory within the binary collision and quasistatic framework has been developed using the coordinate representation. Within this formalism, the main computational task is the evaluation of multidimensional integrals whose variables are the orientational angles needed to specify the initial and final positions of the system during transition processes. Using standard methods, one is able to evaluate the 7-dimensional integrations required for linear molecular systems, or the 7-dimensional integrations for more complicated asymmetric-top (or symmetric-top) molecular systems whose interaction potential contains cyclic coordinates. In order to obviate this latter restriction on the form of the interaction potential, a Monte Carlo method is used to evaluate the 9-dimensional integrations required for systems consisting of one asymmetric-top (or symmetric-top) and one linear molecule, such as H20-N2. Combined with techniques developed previously to deal with sophisticated potential models, one is able to implement realistic potentials for these systems and derive accurate, converged results for the far-wing line shapes and the corresponding absorption coefficients. Conversely, comparison of the far-wing absorption with experimental data can serve as a sensitive diagnostic tool in order to obtain detailed information on the short-range anisotropic dependence of interaction potentials.

Ma, Q.↗
NASA NTRS2011

Effect of Liquid Penetrant Sensitivity on Probability of Detection

The objective of the task is to investigate the effect of liquid penetrant sensitivity level on probability of detection (POD) of cracks in various metals. NASA-STD-5009 currently requires the use of only sensitivity level 4 liquid penetrants for NASA Standard Level inspections. This requirement is based on the fact that the data used to establish the reliably detectable flaw sizes penetrant inspection was from studies performed in the 1970s using penetrant deemed to be equivalent only to modern day sensitivity level 4 penetrants. However, many NDE contractors supporting NASA Centers routinely use sensitivity level 3 penetrants. Because of the new NASA-STD-5009 requirement, these contractors will have to either shift to sensitivity level 4 penetrants or perform formal POD demonstration tests to qualify their existing process. We propose a study to compare the POD generated for two penetrant manufactures, Sherwin and Magnaflux, and for the two most common penetrant inspection methods, water washable and post emulsifiable, hydrophilic. NDE vendors local to GSFC will be employed. A total of six inspectors will inspect a set of crack panels with a broad range of fatigue crack sizes. Each inspector will perform eight inspections of the panel set using the combination of methods and sensitivity levels described above. At least one inspector will also perform multiple inspections using a fixed technique to investigate repeatability. The hit/miss data sets will be evaluated using both the NASA generated DOEPOD software and the MIL-STD-1823 software.

Parker, Bradford H.↗
NASA NTRS2019

Simultaenous Retrieval of Surface Roughness Parameters from Combined Active-Passive SMAP Observations

Soil roughness strongly influences processes like erosion, infiltration, moisture and evaporation of soils as well as growth of agricultural plants. An approach to soil roughness based on active-passive microwave covariation is proposed in order to simultaneously retrieve the vertical RMS height (s) and horizontal correlation length (l) of soil surfaces from simultaneously measured radar and radiometer microwave signatures. The approach is based on a retrieval algorithm for active-passive covariation including the improved Integral Equation Method (I2EM). The algorithm is tested with the global active-passive microwave observations of the SMAP mission. The developed roughness retrieval algorithm shows independence of permittivity for > 10 [-] due to the covariation formalism. Results reveal that s and l can be estimated simultaneously by the proposed approach since surface patterns of non-vegetated areas become evident on global scale. In regions with sandy deserts, like the Sahara or the outback in Australia, determined and confirm rather smooth to semi-rough surface roughness patterns with small vertical RMS heights and corresponding higher horizontal correlation lengths.

correlation length↗
NASA NTRS2023

Radiative interaction of atmosphere and surface: write up with elements of code

In passive satellite remote sensing of the Earth, separation of the path radiance (atmosphere-only contribution) from the surface reflection remains a “significant challenge”. Recent literature names it among the gaps in radiative transfer (RT) topics that “require continued research in the near future”. The challenge comes from multiple reflections (bouncing) between the atmosphere and surface – radiative interaction. In this paper we use a known RT technique, the matrix-operator method (MOM), and a new modification of the monochromatic vector RT (vRT) code IPOL (Intensity and POLarization) to simulate the interaction of a plane-parallel atmosphere and a few widely used surface reflection models. Following the idea of the Green’s function method, IPOL no longer takes the surface model parameters on input. Instead, it provides the path radiance, and the atmospheric reflection and transmission matrices as output. Despite many RT codes use the MOM formalism, this output does not seem common. The surface reflection matrix is computed externally. Therefore, this paper extends the Green’s function atmospheric correction technique to the case of polarized light. Aiming clarity rather than performance, we explain in Python the structure of the surface matrices for the isotropic (Lambertian), directional unpolarized, and polarized ocean reflection models. We then combine these surface matrices and the precomputed IPOL output to get numerically accurate signal at the top of atmosphere (TOA) and test it vs. published benchmarks. Then, for each benchmark scenario we show how to get the surface from the TOA signal, i.e. perform the RT-based atmospheric correction.

radiative transfer↗
DOE OSTI2013

Classifying with confidence from incomplete information.

For this paper, we consider the problem of classifying a test sample given incomplete information. This problem arises naturally when data about a test sample is collected over time, or when costs must be incurred to compute the classification features. For example, in a distributed sensor network only a fraction of the sensors may have reported measurements at a certain time, and additional time, power, and bandwidth is needed to collect the complete data to classify. A practical goal is to assign a class label as soon as enough data is available to make a good decision. We formalize this goal through the notion of reliability—the probability that a label assigned given incomplete data would be the same as the label assigned given the complete data, and we propose a method to classify incomplete data only if some reliability threshold is met. Our approach models the complete data as a random variable whose distribution is dependent on the current incomplete data and the (complete) training data. The method differs from standard imputation strategies in that our focus is on determining the reliability of the classification decision, rather than just the class label. We show that the method provides useful reliability estimates of the correctness of the imputed class labels on a set of experiments on time-series data sets, where the goal is to classify the time-series as early as possible while still guaranteeing that the reliability threshold is met.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2019

An Overview of the United States Department of Energy's Nuclear Criticality Safety Program and Future Challenges

The United States (US) Department of Energy (DOE)/National Nuclear Security Administration (NNSA) Nuclear Criticality Safety Program (NCSP) was formally established in January 1998 as part of the Defense Nuclear Facility Safety Board (DNFSB) Recommendation 97-2 implementation plan. The NCSP was initially focused on seven NCSP tasks: Critical Experiments, Benchmarking, Analytical Methods, Nuclear Data, Training and Qualification, Information Preservation and Dissemination, and Applicable Ranges of Bounding Curves and Data. The first NCSP 5-year plan was published in August 1999. Now, more than 20 years later, NCSP continues to serve the United States and the nuclear criticality safety (NCS) community around five technical program elements (TPEs): Analytical Methods, Information Preservation and Dissemination, Integral Experiments, Nuclear Data, and Training and Education. NCSP has three chartered support groups, the Criticality Safety Support Group, the Nuclear Data Advisory Group, and the DOE Criticality Safety Coordinating Team, all of which provide technical advisement to the NCS community, to DOE, and to the NCSP manager. NCSP has a comprehensive mission and vision document that provides technical and budget priorities for TPE attributes and goals to ensure that the needs of the NCS community and DOE NCS are prioritized appropriately to meet the goals defined by the 1998 DNFSB recommendation. Like any DOE/NNSA program, NCSP faces significant challenges in the future that involves program funding, rising nuclear facility costs, aging critical assembly infrastructure, knowledge retention, and other issues. This paper discusses NCSP and provides some insights into these and other challenges moving into the future.

Bowen, Douglas G.↗
DOE OSTI2016

Electrochemical Electron Transfer and Proton-Coupled Electron Transfer: Effects of Double Layer and Ionic Environment on Solvent Reorganization Energies

Electron transfer and proton coupled electron transfer (PCET) reactions at electrochemical interfaces play an essential role in a broad range of energy conversion processes. The reorganization energy, which is a measure of the free energy change associated with solute and solvent rearrangements, is a key quantity for calculating rate constants for these reactions. We present a computational method for including the effects of the double layer and ionic environment of the diffuse layer in calculations of electrochemical solvent reorganization energies. This approach incorporates an accurate electronic charge distribution of the solute within a molecular-shaped cavity in conjunction with a dielectric continuum treatment of the solvent, ions, and electrode using the integral equations formalism polarizable continuum model. The molecule-solvent boundary is treated explicitly, but the effects of the electrode-double layer and double layer-diffuse layer boundaries, as well as the effects of the ionic strength of the solvent, are included through an external Green’s function. The calculated total reorganization energies agree well with experimentally measured values for a series of electrochemical systems, and the effects of including both the double layer and ionic environment are found to be very small. This general approach was also extended to electrochemical PCET and produced total reorganization energies in close agreement with experimental values for two experimentally studied PCET systems. This research was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center, funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences.

Ghosh, Soumya↗
DOE OSTI2025

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

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

Hydrogen monitoring requirements in the global technical regulation on hydrogen and fuel cell vehicles

Here, the United Nations Economic Commission for Europe Global Technical Regulation (GTR) Number 13 ( Global Technical Regulation on Hydrogen and Fuel Cell Vehicles ) is the defining document regulating safety requirements in hydrogen vehicles, and in particular, fuel cell electric vehicles (FCEVs). GTR Number 13 has been formally adopted and will serve as the basis for the national regulatory standards for FCEV safety in North America (led by the United States), Japan, Korea, and the European Union. The GTR defines safety requirements for these vehicles, including specifications on the allowable hydrogen levels in vehicle enclosures during in-use and post-crash conditions and on the allowable hydrogen emissions levels in vehicle exhaust during certain modes of normal operation. However, in order to be incorporated into national regulations, that is, to be legally binding, methods to verify compliance with the specific requirements must exist. In a collaborative program, the Sensor Laboratories at the National Renewable Energy Laboratory in the United States and the Joint Research Centre, Institute for Energy and Transport in the Netherlands have been evaluating and developing analytical methods that can be used to verify compliance with the hydrogen release requirements as specified in the GTR.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2017

A New Method for Generating Probability Tables in the Unresolved Resonance Region

One new method for constructing probability tables in the unresolved resonance region (URR) has been developed. This new methodology is an extensive modification of the single-level Breit-Wigner (SLBW) pseudo-resonance pair sequence method commonly used to generate probability tables in the URR. The new method uses a Monte Carlo process to generate many pseudo-resonance sequences by first sampling the average resonance parameter data in the URR and then converting the sampled resonance parameters to the more robust R-matrix limited (RML) format. Furthermore, for each sampled set of pseudo-resonance sequences, the temperature-dependent cross sections are reconstructed on a small grid around the energy of reference using the Reich-Moore formalism and the Leal-Hwang Doppler broadening methodology. We then use the effective cross sections calculated at the energies of reference to construct probability tables in the URR. The RML cross-section reconstruction algorithm has been rigorously tested for a variety of isotopes, including 16 O, 19 F, 35 Cl, 56 Fe, 63 Cu, and 65 Cu. The new URR method also produced normalized cross-section factor probability tables for 238 U that were found to be in agreement with current standards. The modified 238 U probability tables were shown to produce results in excellent agreement with several standard benchmarks, including the IEU-MET-FAST-007 (BIG TEN), IEU-MET-FAST-003, and IEU-COMP-FAST-004 benchmarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

X-ray Thomson scattering absolute intensity from the f-sum rule in the imaginary-time domain

We present a formally exact and simulation-free approach for the normalization of X-ray Thomson scattering (XRTS) spectra based on the f-sum rule of the imaginary-time correlation function (ITCF). Our method works for any degree of collectivity, over a broad range of temperatures, and is applicable even in nonequilibrium situations. In addition to giving us model-free access to electronic correlations, this new approach opens up the intriguing possibility to extract a plethora of physical properties from the ITCF based on XRTS experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Line interference effects using a refined Robert-Bonamy formalism: The test case of the isotropic Raman spectra of autoperturbed N{sub 2}

A symmetrized version of the recently developed refined Robert-Bonamy formalism [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)] is proposed. This model takes into account line coupling effects and hence allows the calculation of the off-diagonal elements of the relaxation matrix, without neglecting the rotational structure of the perturbing molecule. The formalism is applied to the isotropic Raman spectra of autoperturbed N{sub 2} for which a benchmark quantum relaxation matrix has recently been proposed. The consequences of the classical path approximation are carefully analyzed. Methods correcting for effects of inelasticity are considered. While in the right direction, these corrections appear to be too crude to provide off diagonal elements which would yield, via the sum rule, diagonal elements in good agreement with the quantum results. In order to overcome this difficulty, a re-normalization procedure is applied, which ensures that the off-diagonal elements do lead to the exact quantum diagonal elements. The agreement between the (re-normalized) semi-classical and quantum relaxation matrices is excellent, at least for the Raman spectra of N{sub 2}, opening the way to the analysis of more complex molecular systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Multi-reference approach to the calculation of photoelectron spectra including spin-orbit coupling

X-ray photoelectron spectra provide a wealth of information on the electronic structure. The extraction of molecular details requires adequate theoretical methods, which in case of transition metal complexes has to account for effects due to the multi-configurational and spin-mixed nature of the many-electron wave function. Here, the restricted active space self-consistent field method including spin-orbit coupling is used to cope with this challenge and to calculate valence- and core-level photoelectron spectra. The intensities are estimated within the frameworks of the Dyson orbital formalism and the sudden approximation. Thereby, we utilize an efficient computational algorithm that is based on a biorthonormal basis transformation. The approach is applied to the valence photoionization of the gas phase water molecule and to the core ionization spectrum of the [Fe(H{sub 2}O){sub 6}]{sup 2+} complex. The results show good agreement with the experimental data obtained in this work, whereas the sudden approximation demonstrates distinct deviations from experiments.

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

cclib 2.0: An updated architecture for interoperable computational chemistry

Interoperability in computational chemistry is elusive, impeded by the independent development of software packages and idiosyncratic nature of their output files. The cclib library was introduced in 2006 as an attempt to improve this situation by providing a consistent interface to the results of various quantum chemistry programs. The shared API across programs enabled by cclib has allowed users to focus on results as opposed to output and to combine data from multiple programs or develop generic downstream tools. Initial development, however, did not anticipate the rapid progress of computational capabilities, novel methods, and new programs; nor did it foresee the growing need for customizability. Here, we recount this history and present cclib 2, focused on extensibility and modularity. We also introduce recent design pivots—the formalization of cclib’s intermediate data representation as a tree-based structure, a new combinator-based parser organization, and parsed chemical properties as extensible objects.

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

EMBEDDED LENSING TIME DELAYS, THE FERMAT POTENTIAL, AND THE INTEGRATED SACHS–WOLFE EFFECT

We derive the Fermat potential for a spherically symmetric lens embedded in a Friedman–Lemaître–Robertson–Walker cosmology and use it to investigate the late-time integrated Sachs–Wolfe (ISW) effect, i.e., secondary temperature fluctuations in the cosmic microwave background (CMB) caused by individual large-scale clusters and voids. We present a simple analytical expression for the temperature fluctuation in the CMB across such a lens as a derivative of the lens’ Fermat potential. This formalism is applicable to both linear and nonlinear density evolution scenarios, to arbitrarily large density contrasts, and to all open and closed background cosmologies. It is much simpler to use and makes the same predictions as conventional approaches. In this approach the total temperature fluctuation can be split into a time-delay part and an evolutionary part. Both parts must be included for cosmic structures that evolve and both can be equally important. We present very simple ISW models for cosmic voids and galaxy clusters to illustrate the ease of use of our formalism. We use the Fermat potentials of simple cosmic void models to compare predicted ISW effects with those recently extracted from WMAP and Planck data by stacking large cosmic voids using the aperture photometry method. If voids in the local universe with large density contrasts are no longer evolving we find that the time delay contribution alone predicts values consistent with the measurements. However, we find that for voids still evolving linearly, the evolutionary contribution cancels a significant part of the time delay contribution and results in predicted signals that are much smaller than recently observed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

COSMIC RAY TRANSPORT WITH MAGNETIC FOCUSING AND THE “TELEGRAPH” MODEL

Cosmic rays (CR), constrained by scattering on magnetic irregularities, are believed to propagate diffusively. However, a well-known defect of diffusive approximation, whereby some of the particles propagate unrealistically fast, has directed interest toward an alternative CR transport model based on the “telegraph” equation. Though, its derivations often lack rigor and transparency leading to inconsistent results. We apply the classic Chapman–Enskog method to the CR transport problem. We show that no “telegraph” (second order time derivative) term emerges in any order of a proper asymptotic expansion with systematically eliminated short timescales. Nevertheless, this term may formally be converted from the fourth order hyper-diffusive term of the expansion. However, both the telegraph and hyperdiffusive terms may only be important for a short relaxation period associated with either strong pitch-angle anisotropy or spatial inhomogeneity of the initial CR distribution. Beyond this period the system evolves diffusively in both cases. The term conversion, that makes the telegraph and Chapman–Enskog approaches reasonably equivalent, is possible only after this relaxation period. During this period, the telegraph solution is argued to be unphysical. Unlike the hyperdiffusion correction, it is not uniformly valid and introduces implausible singular components to the solution. These dominate the solution during the relaxation period. Because they are shown not to be inherent in the underlying scattering problem, we argue that the telegraph term is involuntarily acquired in an asymptotic reduction of the problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Many-body theory of phonon-induced spin relaxation and decoherence

First-principles calculations enable accurate predictions of electronic interactions and dynamics. However, computing the electron spin dynamics remains challenging. The spin-orbit interaction causes various dynamical phenomena that couple with phonons, such as spin precession and spin-flip e-ph scattering, which are difficult to describe with current first-principles calculations. Here, in this work, we show a rigorous framework to study phonon-induced spin relaxation and decoherence, by computing the spin-spin correlation function and its vertex corrections due to e-ph interactions. We apply this approach to a model system and develop corresponding first-principles calculations of spin relaxation in GaAs. Our vertex-correction formalism is shown to capture the Elliott-Yafet, Dyakonov-Perel, and strong-precession mechanisms—three independent spin decoherence regimes with distinct physical origins—thereby unifying their theoretical treatment and calculation. Our method is general and enables quantitative studies of spin relaxation, decoherence, and transport in a wide range of materials and devices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Graphic Processor Realization of the Formal Solution of the Problem of Radiation Transfer in the Neutron Star Atmosphere

The realization (on the CUDA platform) of the search for the formal solution to the problem of radiation transfer in the strongly magnetized plasma in the atmosphere of a neutron star on graphic processors with the CUDA architecture. The solution is obtained using the Galerkin method with finite elements. The realization of graphic processors makes it possible to substantially accelerate computations in constructing models for the atmospheres of neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗