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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2015

Many-body Green’s function theory for electron-phonon interactions: The Kadanoff-Baym approach to spectral properties of the Holstein dimer

We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.

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

CORAL: A framework for rigorous self-validated data modeling and integrative, reproducible data analysis

Abstract Background Many organizations face challenges in managing and analyzing data, especially when relevant datasets arise from multiple sources and methods. Analyzing heterogeneous datasets and additional derived data requires rigorous tracking of their interrelationships and provenance. This task has long been a Grand Challenge of data science and has more recently been formalized in the FAIR principles: that all data objects be Findable, Accessible, Interoperable, and Reusable, both for machines and for people. Adherence to these principles is necessary for proper stewardship of information, for testing regulatory compliance, for measuring the efficiency of processes, and for facilitating reuse of data-analytical frameworks. Findings We present the Contextual Ontology-based Repository Analysis Library (CORAL), a platform that greatly facilitates adherence to all 4 of the FAIR principles, including the especially difficult challenge of making heterogeneous datasets Interoperable and Reusable across all parts of a large, long-lasting organization. To achieve this, CORAL's data model requires that data generators extensively document the context for all data, and our tools maintain that context throughout the entire analysis pipeline. CORAL also features a web interface for data generators to upload and explore data, as well as a Jupyter notebook interface for data analysts, both backed by a common API. Conclusions CORAL enables organizations to build FAIR data types on the fly as they are needed, avoiding the expense of bespoke data modeling. CORAL provides a uniquely powerful platform to enable integrative cross-dataset analyses, generating deeper insights than are possible using traditional analysis tools.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Imposing strong correlated energy spread on relativistic bunches with transverse deflecting cavities

We demonstrate a novel scheme for imposing or removing large energy chirp on relativistic electron bunches. This method relies on the transverse-to-longitudinal mixing accomplished with a set of transverse deflecting cavities. The working principles of the new concept are explained using the first-order matrix formalism in the approximation of the linear single-particle dynamics. The scheme demonstrates versatility regarding the average beam energy and, hence, positioning along the accelerator beam line. The performance of the scheme is numerically investigated with the elegant particle-tracking code. It is shown that the associated nonlinear effects causing emittance deterioration for the extreme quantities of the total energy spread can be effectively compensated by optimization of the Twiss parameters while relying on the eigenemittance analysis. The impact of the longitudinal space charge effects on the beam dynamics in the proposed scheme is also investigated.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Convex optimization of contour deformations

We discuss various formal aspects of contour deformations used to alleviate sign problems; most importantly, relating these contour deformations to a certain convex optimization problem. As a consequence of this connection we describe a general method for proving upper bounds on the average phase achievable by the contour deformation method. Using this method we show that Abelian lattice Yang-Mills in two spacetime dimensions possesses, for many values of the complex coupling, an exponential sign problem that cannot be removed via any contour deformation. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Poster - 11: Radiation barrier thickness calculations for the GammaPod

A consortium of radiotherapy centers in North America is in the process of evaluating a novel new {sup 60}Co teletherapy device, called the GammaPod™ (Xcision Medical Systems, Columbia Maryland), designed specifically for breast SBRT. The GammaPod consists of 36 collimated {sup 60}Co sources with a total activity of 4320 Ci. The sources are housed in a hemispherical source carrier that rotates during treatment to produce a cylindrically symmetric cone of primary beam spanning 16° – 54° degrees from the horizontal. This unique beam geometry presents challenges when designing or evaluating room shielding for the purposes of meeting regulatory requirements, and for ensuring the safety of staff and the public in surrounding areas. Conventional methods for calculating radiation barrier thicknesses have been adapted so that barrier transmission factors for the GammaPod can be determined from a few relevant distances and characteristics of the primary beam. Simple formalisms have been determined for estimating shielding requirements for primary radiation (with a rotating and non-rotating source carrier), patient-scattered radiation, and leakage radiation. When making worst case assumptions, it was found that conventional barrier thicknesses associated with linac treatment suites are sufficient for shielding all sources of radiation from the GammaPod.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1974

Lectures on the scattering of light

The exact (Mie) theory for the scattering of a plane wave by a dielectric sphere is presented. Since this infinite series solution is computationally impractical for large spheres, another formulation is given in terms of an integral equation valid for a bounded, but otherwise general array of scatterers. This equation is applied to the scattering by a single sphere, and several methods are suggested for approximating the scattering cross section in closed form. A tensor scattering matrix is introduced, in terms of which some general scattering theorems are derived. The application of the formalism to multiple scattering is briefly considered.

Saxon, D. S.↗
NASA NTRS1981

Directionality effects in the transfer of X-rays from a magnetized atmosphere: Beam pulse shape

A formalism is presented for radiation transfer in two normal polarization modes in finite and semiinfinite plane parallel uniform atmospheres with a magnetic field perpendicular to the surface and arbitrary propagation angles. This method is based on the coupled integral equations of transfer, including emission, absorption, and scattering. Calculations are performed for atmosphere parameters typical of X-ray pulsars. The directionality of the escaping radiation is investigated for several cases, varying the input distributions. Theoretical pencil beam profiles and X-ray pulse shapes are obtained assuming the radiation is emitted from the polar caps of spinning neutron stars. Implications for realistic models of accreting magnetized X-ray sources are briefly discussed.

Meszaros, P.↗
NASA NTRS1981

Determination of magnetic helicity in the solar wind and implications for cosmic ray propagation

Magnetic helicity (Hm) is the mean value of the correlation between a turbulent magnetic field and the magnetic vector potential. A technique is described for determining Hm and its 'reduced' spectrum from the two point magnetic correlation matrix. The application of the derived formalism to solar wind magnetic fluctuations is discussed, taking into account cases for which only single point measurements are available. The application procedure employs the usual 'frozen in approximation' approach. The considered method is applied to an analysis of several periods of Voyager 2 interplanetary magnetometer data near 2.8 AU. During these periods the correlation length, or energy containing length, was found to be approximately 3 x 10 to the 11th cm

Matthaeus, W. M.↗
DOE OSTI2016

Meson Resonances from Lattice QCD

There has been recent, significant, advances in the determination of the meson spectrum of QCD. Current efforts have focused on the development and application of finite-volume formalisms that allow for the determination of scattering amplitudes as well as resonance behavior in coupled channel systems. I will review some of these recent developments, and demonstrate the viability of the method in meson systems.

Edwards, Robert G.↗