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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

The phase delay and its complex time: From stationary states up to wave packets

Complex time is often invoked about tunneling effect where the classical phase delay is completed with a crucial filter effect. Usually the complex times are obtained by considering the flux-flux correlation function, but this can be obtained by a very simple approach using the search of the maximum of the generalized complex phase function, including the amplitude of the wave function. Various aspects of the phase delay are presented in the case of wave packets impinging on simple or resonant quantum barriers. Formal links with the classical mechanics give birth to quasi-trajectories of the quantum particle, totally compatible with the quantum mechanics. - Highlights: Black-Right-Pointing-Pointer The stationary phase method is extended in including the variations of the spectra. Black-Right-Pointing-Pointer The complex phase delay leads to a complex trajectory inside and out-side the barrier. Black-Right-Pointing-Pointer Examples of quasi-trajectories are given in case of different quantum barriers. Black-Right-Pointing-Pointer Phase delays are specified for resonant tunneling or above-barrier wave-packets. Black-Right-Pointing-Pointer The coherence between the quasi-trajectories and quantum mechanics is shown.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Analysis of calibration data for the uranium active neutron coincidence counting collar with attention to errors in the measured neutron coincidence rate

The declared linear density of 238 U and 235 U in fresh low enriched uranium light water reactor fuel assemblies can be verified for nuclear safeguards purposes using a neutron coincidence counter collar in passive and active mode, respectively. The active mode calibration of the Uranium Neutron Collar - Light water reactor fuel (UNCL) instrument is normally performed using a non-linear fitting technique. The fitting technique relates the measured neutron coincidence rate (the predictor) to the linear density of 235 U (the response) in order to estimate model parameters of the nonlinear Parte equation, which traditionally is used to model the calibration data. Alternatively, following a simple data transformation, the fitting can also be performed using standard linear fitting methods. This article compares performance of the nonlinear technique to the linear technique, using a range of possible error variance magnitudes in the measured neutron coincidence rate. We develop the required formalism and then apply the traditional (nonlinear) and alternative approaches (linear) to the same experimental and corresponding simulated representative datasets. Here, we find that, in this context, because of the magnitude of the errors in the predictor, it is preferable not to transform to a linear model, and it is preferable not to adjust for the errors in the predictor when inferring the model parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Combined Effects of Hemicolligation and Ion Pairing on Reduction Potentials of Biphenyl Radical Cations

Here, formal reduction potentials of highly oxidizing and short-lived radical cations of substituted biphenyls generated by pulse radiolysis in 1,2-dichloroethane (DCE) were measured using a redox equilibrium ladder method. The effect of halide ion–radical interactions on reduction potentials of biphenyls was examined by utilizing the ability of DCE to release Cl – in the vicinity of the radical cation. The Hammett correlation of measured potentials across a range of over 700 mV shows saturation at high Hammett sigma values. This effect has been explained by both ion-pairing and hemicolligation interactions between biphenyl radical cations and Cl – and appears to modulate reduction potentials by as much as 400 mV. This finding offers a convenient way to manipulate the energetics of electron transfer involving organic redox species.

14 SOLAR ENERGY↗
DOE OSTI2020

Coupled Cluster Green's function formulations based on the effective Hamiltonians

In this work, we demonstrate that the effective Hamiltonians obtained with the downfolding procedure based on double unitary coupled cluster (DUCC) ansatz can be used in the context of Green’s function coupled cluster (GFCC) formalism to calculate spectral functions of molecular systems. This combined approach (DUCC-GFCC) provides a significant reduction of numerical effort and good agreement with the corresponding all-orbital GFCC methods in energy windows that are consistent with the choice of active space. These features are demonstrated on the example of two benchmark systems: H 2 O and N 2 , where DUCC-GFCC calculations were performed for active spaces of various sizes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

{chi}{sup 2} DISCRIMINATORS FOR TRANSITING PLANET DETECTION IN KEPLER DATA

The Kepler spacecraft observes a host of target stars to detect transiting planets. Requiring a 7.1{sigma} detection in three years of data yields over 100,000 detections, many of which are false alarms. After a second cut is made on a robust detection statistic, some 50,000 or more targets still remain. These false alarms waste resources as they propagate through the remainder of the software pipeline and so a method to discriminate against them is crucial in maintaining the desired sensitivity to true events. This paper describes a {chi}{sup 2} test which represents a novel application of an existing formalism developed for false alarm mitigation in searches for gravitational waves. Using this technique, the false alarm rate can be lowered to {approx}5%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-E-T-447: Electronic Brachytherapy (EBT) Treatment of Cervical Cancer - First Clinical Experience

Purpose: To study the first trial patient in which an electronic brachytherapy (EBT) x-ray source is utilized for treatment of cervical cancer. Methods: During patient treatment, a miniaturized x-ray source was used in combination with a customized titanium tandem and ovoid applicator set. The semi-specialized source was modeled with formalisms outlined by AAMP Task Group 43. Multiple models were used to compensate for variable attenuation conditions as a function of source positions. Varian Brachyvision treatment planning software was utilized on CT data sets for dose calculations prior to treatment delivery. The dose was prescribed to “point A” as defined by American Brachytherapy society. Additional treatments plans were created from those clinically utilized in patient care and were recalculated for an existing Ir-192 source model. Dose volume histograms (DVH) and point dose calculations were compared between the modalities for the clinical condition present in patients treated with EBT. Results: Clinical treatment times, though longer than those typically experienced by Ir-192 users, were manageable. Instantaneous dose rates at personal positions within the treatment vault were lower than those measured during intra operative radiation therapy and breast EBT treatments. Due to lower average photon energy in EBT, dose gradients within the treatment plans were as expected steeper than those observed in Ir-192 based brachytherapy. DVH comparisons between Ir-192 and EBT treatments showed an expected decrease in the integral dose to normal tissues of interest for EBT. In comparing plans created for EBT delivery with those calculated for Ir-192, average dose values for EBT were more than 4%, 11%, and 9% lower at predefined bladder, rectum and “point B” positions, respectively. Conclusion: For the first time, we have demonstrated that the utilizing electronic brachytherapy system for tandem and ovoid based treatment of cancer of the cervix is feasible, and potentially beneficial, in clinical environment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A generic high-dose rate {sup 192}Ir brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism

Purpose: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) {sup 192}Ir source and a virtual water phantom were designed, which can be imported into a TPS. Methods: A hypothetical, generic HDR {sup 192}Ir source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic {sup 192}Ir source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra{sup ®} Brachy with advanced collapsed-cone engine (ACE) and BrachyVision ACUROS{sup TM}]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and PENELOPE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201){sup 3} voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR {sup 192}Ir source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. Results: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ACE algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 ± 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ACE at clinically relevant distances. Conclusions: A hypothetical, generic HDR {sup 192}Ir source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

P–V–T Equation of State of Iridium Up to 80 GPa and 3100 K

In the present study, the high-pressure high-temperature equation of the state of iridium has been determined through a combination of in situ synchrotron X-ray diffraction experiments using laser-heating diamond-anvil cells (up to 48 GPa and 3100 K) and density-functional theory calculations (up to 80 GPa and 3000 K). The melting temperature of iridium at 40 GPa was also determined experimentally as being 4260 (200) K. The results obtained with the two different methods are fully consistent and agree with previous thermal expansion studies performed at ambient pressure. The resulting thermal equation of state can be described using a third-order Birch–Murnaghan formalism with a Berman thermal-expansion model. The present equation of the state of iridium can be used as a reliable primary pressure standard for static experiments up to 80 GPa and 3100 K. A comparison with gold, copper, platinum, niobium, rhenium, tantalum, and osmium is also presented. On top of that, the radial-distribution function of liquid iridium has been determined from experiments and calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Formal Solutions for Polarized Radiative Transfer. I. The DELO Family

The discussion regarding the numerical integration of the polarized radiative transfer equation is still open and the comparison between the different numerical schemes proposed by different authors in the past is not fully clear. Aiming at facilitating the comprehension of the advantages and drawbacks of the different formal solvers, this work presents a reference paradigm for their characterization based on the concepts of order of accuracy , stability , and computational cost . Special attention is paid to understand the numerical methods belonging to the Diagonal Element Lambda Operator family, in an attempt to highlight their specificities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Recent Progress on the Stochastic Objective Decision Aide (SODA) Application

Idaho State University, in collaboration with Idaho National Laboratory, is engaged in development of a simple to use software application that utilizes stochastic techniques for the calculation of radiation doses associated with hypothetical radioactive material releases. This work is being performed through a grant from the U.S. Department of Energy (DOE) Nuclear Safety Research and Development (NSR and D) Program. Initial progress on the application was reported previously. Recent additions to the application are reported here. The Stochastic Objective Decision Aide (SODA) provides an easy to use software application to help decision makers support their judgment about the need or lack of need for safety structures, systems, and components. SODA is not intended to replace traditional design basis accident dose calculation methodology, nor is it intended to replace existing software used to perform formal dose consequence calculations. Rather, SODA is intended to supplement information obtained from existing codes to give a decision maker more information about an accident scenario. SODA also provides a simple method to perform parametric or sensitivity studies on individual input parameters. Non-reactor nuclear facilities operated under the approval authority of the US DOE use unmitigated hazard evaluations to determine if potential radiological doses associated with design basis events challenge dose evaluation guidelines. Unmitigated events that sufficiently challenge dose evaluation guidelines merit the selection of safety structures, systems, or components (SSCs) to prevent or mitigate the hazard. Dose consequence calculations traditionally involve use of the 'five-factor' formula: ST = MAR.DR.ARF.RF.LPF where ST is the source term (Bq), MAR is the total available material-at-risk (Bq), DR is the damage ratio (no units), ARF is the airborne release fraction (no units), RF is the respirable fraction (no units), and LPF is the leak path factor (no units). Potential radiation doses are then calculated using: CED = χ/Q.BR.ST.DCF (2) where CED is the committed effective dose (Sv), χ/Q is the plume dispersion (s/m{sup 3}), BR is the breathing rate (m{sup 3}/s), ST is the source term (Bq), and DCF is the dose conversion factor (Sv/Bq). Conservative single value input parameters are typically used to represent ARF, RF, LPF and BR. The traditional methodology, while conservative, can lead to skewed conclusions in the balance between cost and risk reduction resulting in over engineered systems with greater design, construction and operating costs. Rather than using a bounding single point value for each parameter in the dose consequence calculation, SODA uses distributions for some or all of the parameters. Each input parameter distribution is stochastically sampled and the resulting dose consequence calculation is repeated many times to develop a dose consequence distribution. The resulting distribution provides more information for the decision maker to rely upon when determining the need for safety SSCs. The SODA application was developed using MATLAB and it incorporates the use of Monte Carlo techniques as well as a graphical user interface (GUI). The application includes user selection of the governing distribution for parameters such as the MAR, DR, ARF, RF, LPF,χ/Q, BR, and DCF. While MATLAB was used for the development work, the application is distributed as a self-contained executable program. Features of the application include pull down menus with available distributions for the various parameters. The user has the option to plot the input parameter distribution that results from the random sampling process. The resulting dose consequence distribution is automatically analyzed to determine its similarity to well know distributions. This was done using the Bayesian Information Criterion method. If the resulting distribution is sufficiently similar to a well-known distribution, the application provides information about the distributions such as the mean value and other parameters. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

SU-D-201-07: A Survey of Radiation Oncology Residents’ Training and Preparedness to Lead Patient Safety Programs in Clinics

Purpose: Safety and quality has garnered increased attention in radiation oncology, and physicians and physicists are ideal leaders of clinical patient safety programs. However, it is not clear whether residency programs incorporate formal patient safety training and adequately equip residents to assume this leadership role. A national survey was conducted to evaluate medical and physics residents’ exposure to safety topics and their confidence with the skills required to lead clinical safety programs. Methods: Radiation oncology residents were identified in collaboration with ARRO and AAPM. The survey was released in February 2016 via email using REDCap. This included questions about exposure to safety topics, confidence leading safety programs, and interest in training opportunities (i.e. workshops). Residents rated their exposure, skills, and confidence on 4 or 5-point scales. Medical and physics residents responses were compared using chi-square tests. Results: Responses were collected from 56 of 248 (22%) physics and 139 of 690 (20%) medical residents. More than two thirds of all residents had no or only informal exposure to incident learning systems (ILS), root cause analysis (RCA), failure mode and effects analysis (FMEA), and the concept of human factors engineering (HFE). Likewise, 63% of residents had not heard of RO-ILS. Response distributions were similar, however more physics residents had formal exposure to FMEA (p<0.0001) and felt they were adequately trained to lead FMEAs in clinic (p<0.001) than medical residents. Only 36% of residents felt their patient safety training was adequate, and 58% felt more training would benefit their education. Conclusion: These results demonstrate that, despite increasing desire for patient safety training, medical and physics residents’ exposure to relevant concepts is low. Physics residents had more exposure to FMEA than medical residents, and were more confident in leading FMEA. This suggests that increasing resident exposure to specific topics may increase their confidence.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

A Structured Quasi-Newton Algorithm for Optimizing with Incomplete Hessian Information

We present a structured quasi-Newton algorithm for unconstrained optimization problems that have unavailable second-order derivatives or Hessian terms. We provide a formal derivation of the well-known Broyden--Fletcher--Goldfarb--Shanno (BFGS) secant update formula that approximates only the missing Hessian terms, and we propose a linesearch quasi-Newton algorithm based on a modification of Wolfe conditions that converges to first-order optimality conditions. We also analyze the local convergence properties of the structured BFGS algorithm and show that it achieves superlinear convergence under the standard assumptions used by quasi-Newton methods using secant updates. In conclusion, we provide a thorough study of the practical performance of the algorithm on the CUTEr suite of test problems and show that our structured BFGS-based quasi-Newton algorithm outperforms the unstructured counterpart(s).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Kernel reconstruction methods for Doppler broadening — Temperature interpolation by linear combination of reference cross sections at optimally chosen temperatures

This paper establishes a new family of methods to perform temperature interpolation of nuclear interactions cross sections, reaction rates, or cross sections times the energy. One of these quantities at temperature T is approximated as a linear combination of quantities at reference temperatures (T j ). The problem is formalized in a cross section independent fashion by considering the kernels of the different operators that convert cross section related quantities from a temperature T 0 to a higher temperature T — namely the Doppler broadening operation. Doppler broadening interpolation of nuclear cross sections is thus here performed by reconstructing the kernel of the operation at a given temperature T by means of linear combination of kernels at reference temperatures (T j ). The choice of the L 2 metric yields optimal linear interpolation coefficients in the form of the solutions of a linear algebraic system inversion. The optimization of the choice of reference temperatures (T j ) is then undertaken so as to best reconstruct, in the L∞ sense, the kernels over a given temperature range [T min ,T max ]. The performance of these kernel reconstruction methods is then assessed in light of previous temperature interpolation methods by testing them upon isotope 238 U. Temperature-optimized free Doppler kernel reconstruction significantly outperforms all previous interpolation-based methods, achieving 0.1% relative error on temperature interpolation of 238 U total cross section over the temperature range [300 K,3000 K] with only 9 reference temperatures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

DETECTION OF LENSING SUBSTRUCTURE USING ALMA OBSERVATIONS OF THE DUSTY GALAXY SDP.81

We study the abundance of substructure in the matter density near galaxies using ALMA Science Verification observations of the strong lensing system SDP.81. We present a method to measure the abundance of subhalos around galaxies using interferometric observations of gravitational lenses. Using simulated ALMA observations we explore the effects of various systematics, including antenna phase errors and source priors, and show how such errors may be measured or marginalized. We apply our formalism to ALMA observations of SDP.81. We find evidence for the presence of a M = 10{sup 8.96±0.12} M {sub ⊙} subhalo near one of the images, with a significance of 6.9 σ in a joint fit to data from bands 6 and 7; the effect of the subhalo is also detected in both bands individually. We also derive constraints on the abundance of dark matter (DM) subhalos down to M ∼ 2 × 10{sup 7} M {sub ⊙}, pushing down to the mass regime of the smallest detected satellites in the Local Group, where there are significant discrepancies between the observed population of luminous galaxies and predicted DM subhalos. We find hints of additional substructure, warranting further study using the full SDP.81 data set (including, for example, the spectroscopic imaging of the lensed carbon monoxide emission). We compare the results of this search to the predictions of ΛCDM halos, and find that given current uncertainties in the host halo properties of SDP.81, our measurements of substructure are consistent with theoretical expectations. Observations of larger samples of gravitational lenses with ALMA should be able to improve the constraints on the abundance of galactic substructure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A comparative study of Rosenbrock-type and implicit Runge-Kutta time integration for discontinuous Galerkin method for unsteady 3D compressible Navier-Stokes equations

A comparative study of two classes of third-order implicit time integration schemes is presented for a third-order hierarchical WENO reconstructed discontinuous Galerkin (rDG) method to solve the 3D unsteady compressible Navier-Stokes equations: — 1) the explicit first stage, single diagonally implicit Runge-Kutta (ESDIRK3) scheme, and 2) the Rosenbrock-Wanner (ROW) schemes based on the differential algebraic equations (DAEs) of Index-2. Compared with the ESDIRK3 scheme, a remarkable feature of the ROW schemes is that, they only require one approximate Jacobian matrix calculation every time step, thus considerably reducing the overall computational cost. A variety of test cases, ranging from inviscid flows to DNS of turbulent flows, are presented to assess the performance of these schemes. Here, numerical experiments demonstrate that the third-order ROW scheme for the DAEs of index-2 can not only achieve the designed formal order of temporal convergence accuracy in a benchmark test, but also require significantly less computing time than its ESDIRK3 counterpart to converge to the same level of discretization errors in all of the flow simulations in this study, indicating that the ROW methods provide an attractive alternative for the higher-order time-accurate integration of the unsteady compressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A fourth-order accurate finite volume method for ideal MHD via upwind constrained transport

Here we present a fourth-order accurate finite volume method for the solution of ideal magnetohydrodynamics (MHD). The numerical method combines high-order quadrature rules in the solution of semi-discrete formulations of hyperbolic conservation laws with the upwind constrained transport (UCT) framework to ensure that the divergence-free constraint of the magnetic field is satisfied. A novel implementation of UCT that uses the piecewise parabolic method (PPM) for the reconstruction of magnetic fields at cell corners in 2D is introduced. The resulting scheme can be expressed as the extension of the second-order accurate constrained transport (CT) Godunov-type scheme that is currently used in the Athena astrophysics code. After validating the base algorithm on a series of hydrodynamics test problems, we present the results of multidimensional MHD test problems which demonstrate formal fourth-order convergence for smooth problems, robustness for discontinuous problems, and improved accuracy relative to the second-order scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

A method for the construction of stable Galileon models consistent with the Planck data results

The reconstruction procedure, which has proven quite useful to obtain viable models of the universe evolution, is here employed in order to construct inflation models. It has the advantages that it ensures full consistency with astronomical observations and that it allows to evaluate the stability of the resulting cosmological model in a rather simple way. The reconstruction is carried out in detail for two different types of Lagrangian, included in the frame of G-inflation, and explicit models for each Lagrangian are obtained. As a bonus for having used this reconstruction formalism, the final models are easily adjusted to satisfy the observational constraints—imposed by the most recent data releases of the Planck mission—on the spectral index, the tensor to scalar ratio, and the running of the spectral index. Further, it turns also to be not difficult to impose the models to be stable. Thus, the method here developed provides a general and very efficient tool, a quite natural procedure in order to construct models consistent with very precise observations. It can be also applied to other models, besides the ones here considered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Classical nuclear motion coupled to electronic non-adiabatic transitions

Based on the exact factorization of the electron-nuclear wave function, we have recently proposed a mixed quantum-classical scheme [A. Abedi, F. Agostini, and E. K. U. Gross, Europhys. Lett. 106, 33001 (2014)] to deal with non-adiabatic processes. Here we present a comprehensive description of the formalism, including the full derivation of the equations of motion. Numerical results are presented for a model system for non-adiabatic charge transfer in order to test the performance of the method and to validate the underlying approximations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗