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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Testing quantum gravity through dumb holes

We propose a method to test the effects of quantum fluctuations on black holes by analyzing the effects of thermal fluctuations on dumb holes, the analogs for black holes. The proposal is based on the Jacobson formalism, where the Einstein field equations are viewed as thermodynamical relations, and so the quantum fluctuations are generated from the thermal fluctuations. It is well known that all approaches to quantum gravity generate logarithmic corrections to the entropy of a black hole and the coefficient of this term varies according to the different approaches to the quantum gravity. It is possible to demonstrate that such logarithmic terms are also generated from thermal fluctuations in dumb holes. In this paper, we claim that it is possible to experimentally test such corrections for dumb holes, and also obtain the correct coefficient for them. This fact can then be used to predict the effects of quantum fluctuations on realistic black holes, and so it can also be used, in principle, to experimentally test the different approaches to quantum gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Monte Carlo chord length sampling for d-dimensional Markov binary mixtures

The Chord Length Sampling (CLS) algorithm is a powerful Monte Carlo method that models the effects of stochastic media on particle transport by generating on-the-fly the material interfaces seen by the random walkers during their trajectories. This annealed disorder approach, which formally consists of solving the approximate Levermore–Pomraning equations for linear particle transport, enables a considerable speed-up with respect to transport in quenched disorder, where ensemble-averaging of the Boltzmann equation with respect to all possible realizations is needed. However, CLS intrinsically neglects the correlations induced by the spatial disorder, so that the accuracy of the solutions obtained by using this algorithm must be carefully verified with respect to reference solutions based on quenched disorder realizations. When the disorder is described by Markov mixing statistics, such comparisons have been attempted so far only for one-dimensional geometries, of the rod or slab type. In this work we extend these results to Markov media in two-dimensional (extruded) and three-dimensional geometries, by revisiting the classical set of benchmark configurations originally proposed by Adams, Larsen and Pomraning and extended by Brantley. In particular, we examine the discrepancies between CLS and reference solutions for scalar particle flux and transmission/reflection coefficients as a function of the material properties of the benchmark specifications and of the system dimensionality.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Mechanism of Catechol Oxidation by the Coupled Bicupric Active Site of Tyrosinase: Completion of the Oxygenase/Oxidase Reaction Cycle

Tyrosinase contains a coupled binuclear copper (CBC) active site, which in its bicuprous form (deoxy-Ty) binds O 2 to form a side-on peroxide [Cu(II) 2 O 2 ] 2+ intermediate (oxy-Ty) that performs the regioselective monooxygenation of monophenols to catechols and their subsequent 2e − oxidation to quinones. Previously, we used spectroscopic, kinetic, and computational methods to elucidate the mechanism of the initial steps in the monooxygenation reaction. Herein, we investigated the final step in catecholate formation and provide experimental and computational results elucidating the catechol oxidation reaction, formally a proton-coupled 2e − process. We employed single-turnover stoppedflow absorption to observe the elusive catecholate-Cu(II) 2 OH intermediate via the anaerobic reaction of resting Cu(II) 2 OH Ty (met-Ty) with the catecholate substrate. This intermediate was cryo-trapped and characterized by electron paramagnetic resonance and X-ray absorption spectroscopies. These experimental results were correlated to quantum mechanics/molecular mechanics (QM/MM) and QM calculations to describe the complete catalytic cycle of tyrosinase, revealing that the 2e − oxidation of catechol to quinone involves two steps: (i) coupled 1H + /1e − transfer from the bound monoanionic catecholate to the bridged hydroxide of the met-Ty active site, followed by (ii) the second 1e − transfer. The latter is calculated to be the rate-limiting step in catechol oxidation, confirmed by experimental solvent kinetic isotope effect studies. Our computational data suggest that quinone release from the protein provides the necessary driving force for this reaction. These results establish a detailed description of this oxidase cycle of the tyrosinase mechanism and more broadly provide molecular-level insights into the diverse reactivity of CBC sites in biology.

QM/MM modeling↗
DOE OSTI2024

The cluster decomposition of the configurational energy of multicomponent alloys

Abstract The cluster expansion method (CEM) is a widely used lattice-based technique in the study of multicomponent alloys. Despite its prevalent use, a clear understanding of expansion terms is lacking. We present a modern mathematical formalism of the CEM and introduce thecluster decomposition—a unique and basis-independent decomposition for functions of the atomic configuration in a crystal. We identify the cluster decomposition as an invariant ANOVA decomposition; and demonstrate how functional analysis of variance and sensitivity analysis can be used to interpret interactions among species. Furthermore, we show how the mathematical structure of the cluster decomposition enables numerical evaluation that scales with the number of clusters and is independent of the number of species. Overall, our work enables rigorous interpretations of interactions among species, provides opportunities to explore parameter estimation beyond linear regression, introduces a numerical efficient implementation, and enables analysis of cluster expansions based on established mathematical and statistical principles.

Chemistry↗
DOE OSTI2016

SU-F-T-659: Nanoparticle-Aided Eye Plaque Radiotherapy

Purpose: Eye plaque brachytherapy is one of the approaches for radiotherapy treatment for ocular cancers: retinoblastoma and choroidal melanoma. This study, investigates the potential benefits of using gold nanoparticles to enhance therapeutic efficacy during eye plaque brachytherapy. Methods: The EYE PHYSICS Inc. Plaque Simulator program distributed by IsoAid, LLC, Port Richey, Florida was used. It is based on the superposition of dose contributions from individual seeds following the TG–43 formalism. Dose enhancement factor (DEF) values for feasible nanoparticle concentrations from previous studies was used to investigate the benefit of using nanoparticles to enhance dose to tumour or reduce dose to healthy tissue. The dose enhancement factor (DEF) represents the ratio of the dose deposited in tumour with nanoparticles divided by dose deposited in the tumour without nanoparticles. The investigation was done for I–125 and Pd–103 typical sources employed for eye plaque brachytherapy. The prescription dose used is 85 Gy. Results: Lower dose enhancement values were obtained for Pd–103. With DEF of 2 due to gold nanoparticles, critical structure doses reduce by a factor of 2. Optic disc dose is 6.69 Gy and 4.571 Gy, opposite retina dose is 4.064 and 2.484 Gy, lens dose is 12.66 Gy and 9.870 Gy, and fovea dose is 9.85 Gy and 7.275 Gy. With DEF of 3 due to gold nanoparticles, critical structure doses reduce by a factor of 3. Optic disc dose is 4.352 Gy and 2.975 Gy, opposite retina dose is 2.644 Gy and 1.618 Gy, lens dose is 8.322 Gy and 6.427 Gy, and fovea dose is 4.815 Gy and 4.737 Gy. Conclusion: The results of this research predict that using gold nanoparticles will lead to major sparing of dose to critical structures. The finding provides more impetus for the development of nanoparticle–aided brachytherapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-201-17: Verification of Dose Distributions From High-Dose-Rate Brachytherapy Ir-192 Source Using a Multiple-Array-Diode-Detector (MapCheck2)

Purpose: To investigate quantitatively the accuracy of dose distributions for the Ir-192 high-dose-rate (HDR) brachytherapy source calculated by the Brachytherapy-Planning system (BPS) and measured using a multiple-array-diode-detector in a heterogeneous medium. Methods: A two-dimensional diode-array-detector system (MapCheck2) was scanned with a catheter and the CT-images were loaded into the Varian-Brachytherapy-Planning which uses TG-43-formalism for dose calculation. Treatment plans were calculated for different combinations of one dwell-position and varying irradiation times and different-dwell positions and fixed irradiation time with the source placed 12mm from the diode-array plane. The calculated dose distributions were compared to the measured doses with MapCheck2 delivered by an Ir-192-source from a Nucletron-Microselectron-V2-remote-after-loader. The linearity of MapCheck2 was tested for a range of dwell-times (2–600 seconds). The angular effect was tested with 30 seconds irradiation delivered to the central-diode and then moving the source away in increments of 10mm. Results: Large differences were found between calculated and measured dose distributions. These differences are mainly due to absence of heterogeneity in the dose calculation and diode-artifacts in the measurements. The dose differences between measured and calculated due to heterogeneity ranged from 5%–12% depending on the position of the source relative to the diodes in MapCheck2 and different heterogeneities in the beam path. The linearity test of the diode-detector showed 3.98%, 2.61%, and 2.27% over-response at short irradiation times of 2, 5, and 10 seconds, respectively, and within 2% for 20 to 600 seconds (p-value=0.05) which depends strongly on MapCheck2 noise. The angular dependency was more pronounced at acute angles ranging up to 34% at 5.7 degrees. Conclusion: Large deviations between measured and calculated dose distributions for HDR-brachytherapy with Ir-192 may be improved when considering medium heterogeneity and dose-artifact of the diodes. This study demonstrates that multiple-array-diode-detectors provide practical and accurate dosimeter to verify doses delivered from the brachytherapy Ir-192-source.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

On the local non-conforming virtual element spaces

The construction of the local non-conforming virtual element spaces is discussed and the isomorphism between their functions and the polynomial moments of such functions is established by a formal argument and two different constructive proofs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

Precise relative location estimation from satellite laser observations

Laser data from the Greenbelt, Bermuda, and Grand Turk Island tec tracking stations, and from observation at Patrick Air Force Base were used in an experiment to determine precise inter-station baselines and relative heights in short orbital arcs of no more than 12 minutes duration. The relative locations obtained are sensitive to reasonable variations in the gravity field to less than 15 cm in inter-station baseline and less than 35 cm in relative height. The formal standard deviations for range measurements at the observed 7 cm noise level are less than 10 cm for either baseline or relative height determinations from the twenty-arc combination solution. The method of data analysis is described, the baseline and relative height results are presented, and gravity model dependence is considered.

Dunn, P. J.↗
NASA NTRS1991

Modeling of electrically thick materials theoretical and experimental aspects

The shielding effectiveness of modern aircraft depends upon the nature and the electromagnetic properties of the materials used. It was shown that thin materials are modeled well by a sheet surface impedance. Unfortunately, this formalism fails in the high frequency domain (when skin effect occurs) and when the materials are physically thick. The purpose is to propose a general formalism that renders the skin effect and can be used with multilayered samples. These electrically thick materials are studied using symmetrical boundary conditions applied to tangential fields. The domain of validity of the method is discussed and experimental validations are presented.

Gobin, V.↗
NASA NTRS1994

Higher order Godunov schemes for isothermal hydrodynamics

In this paper we construct higher order Godunov schemes for isothermal flow. Isothermal hydrodynamics serves as a good representation for several systems of astrophysical interest. The schemes designed here have second-order accuracy in space and time and some are third-order accurate for advection. Moreover, several ingredients of these schemes are essential components of even higher order. The methods designed here have excellent ability to represent smooth flow yet capture shocks with high resolution. Several test problems are presented. The algorithms presented here are compared with other algorithms having a comparable formal order of accuracy.

Balsara, Dinshaw S.↗
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 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↗