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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Synthesis from Design Requirements of a Hybrid System for Transport Aircraft Longitudinal Control

Volume I of this report presents a new method for synthesizing hybrid systems directly from desi gn requirements, and applies the method to design of a hybrid system for longitudinal control of transport aircraft. The resulting system satisfies general requirement for safety and effectiveness specified a priori, enabling formal validation to be achieved. Volume II contains seven appendices intended to make the report accessible to readers with backgrounds in human factors, flight dynamics and control, and formal logic. Major design goals are (1) system design integrity based on proof of correctness at the design level, (2) significant simplification and cost reduction in system development and certification, and (3) improved operational efficiency, with significant alleviation of human-factors problems encountered by pilots in current transport aircraft. This report provides for the first time a firm technical basis for criteria governing design and certification of avionic systems for transport aircraft. It should be of primary interest to designers of next-generation avionic systems.

Hynes, Charles S.↗
NASA NTRS2017

Norton-Thevenin Receptance Coupling (NTRC) as a Payload Design Tool

The NASA Engineering and Safety Center (NESC) is funding a study to develop an alternate method for performing coupled loads analysis called Norton-Thevenin Receptance Coupling (NTRC). NTRC combines Receptance Coupling (RC), a frequency-domain synthesis method and Norton-Thevenin (NT) theory, an impedance based approach for simulating the interaction between dynamic systems. The goal of developing the NTRC method is to provide a tool that payload developers can use to reduce the conservatism in defining preliminary design loads, assess the impact of design changes between formal load cycles, and to perform trade studies for design optimization with a minimum amount of data required from the launch vehicle (LV) provider. NTRC also has the ability to perform parametric loads analysis where many different design configurations can be evaluated. This will result in cost and schedule benefits to the payload developer that are currently not possible under the standard coupled loads analysis (CLA) flow where typically only 2-3 official load cycles are performed by the LV provider over the life of a payload program. NTRC is not envisioned as a replacement for the official load cycles performed by the LV provider but rather as a means to address the types of design issues faced by the payload developer before and between official load cycles.The presentation provides an overview of the NTRC methodology and discusses how NTRC can be used to replicate the results from a standard LV CLA. The presentation covers the benchmarking that has been performed as part of the NESC study to demonstrate the accuracy of the technique for both frequency and time domain dynamic analyses. Future plans for benchmarking the NTRC approach against CLA results for NASAs Space Launch System (SLS) and commercial launch vehicles are discussed and the role that NTRC is envisioned to play in the payload development cycle.

Structural Dynamics↗
NASA NTRS2014

Line Interference Effects Using a Refined Robert-Bonamy Formalism: the Test Case of the Isotropic Raman Spectra of Autoperturbed N2

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 N2 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 N2, opening the way to the analysis of more complex molecular systems.

nitrogen↗
DOE OSTI2020

Calculation of the First Moment of Energy Using D-T Reactivity Formalisms Under the Maxwell-Boltzmann Distribution (Pt. 1)

Nuclear fusion science is an example of a scientific field with a rich history of expert involvement and scientific publications, which together, form an expert-knowledge base. One example of this history is the utilization of published reaction rates from a variety of authors. Investigators for Deuterium-Tritium (D-T) ion fusion can choose from using frequently cited methods: the Bosch and Hale reactivity, thermonuclear reaction rates from Caughlan and Fowler, and the reactivity evaluation from Miley, Towner & Ivich which forms the basis of the Naval Research Lab (NRL) formulary. There are other choices available. Each of the reactivity formulations considered here, are based upon the Maxwell- Boltzmann velocity distribution for D-T fusion ion reactants. Numerical methods for computer codes simulating hot, energetic plasmas, include tabulations of the reactivity, and the first moment of energy. The purpose of this study is to answer the question: what is the first moment of energy, and how has it been formulated? This report is part of a series of reports by the authors on D-T fusion-reaction formalism. The present focus is on defining the mathematical relationship of the first moment of D-T fusion ion kinetic energy, $\langle$E$\rangle$, with the fusion crosssection, fusion reactivity and its derivative with respect to ion-temperature. Three variants of the first moment $\langle$E$\rangle$ are analytically developed and explored: 1) constant cross-section, 2) a normalized first moment, and 3) a particular function of the first moment. Much of the mathematical development is relegated to the Appendices, for a concise presentation of topics and results. A research component of this discussion is our comparison of some of the reactivity and its derivative results between evaluations (1) versus (2) and (3), which the authors have not found in any other publication to date. This investigation is useful for validation, verification (V&V) and uncertainty quantification (UQ) as it compares the work of several authors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fock-Space Schrieffer–Wolff Transformation: Classically-Assisted Rank-Reduced Quantum Phase Estimation Algorithm

We present an extension of many-body downfolding methods to reduce the resources required in the quantum phase estimation (QPE) algorithm. In this paper, we focus on the Schrieffer–Wolff (SW) transformation of the electronic Hamiltonians for molecular systems that provides significant simplifications of quantum circuits for simulations of quantum dynamics. We demonstrate that by employing Fock-space variants of the SW transformation (or rank-reducing similarity transformations (RRST)) one can significantly increase the locality of the qubit-mapped similarity-transformed Hamiltonians. The practical utilization of the SW-RRST formalism is associated with a series of approximations discussed in the manuscript. In particular, amplitudes that define RRST can be evaluated using conventional computers and then encoded on quantum computers. The SW-RRST QPE quantum algorithms can also be viewed as an extension of the standard state-specific coupled-cluster downfolding methods to provide a robust alternative to the traditional QPE algorithms to identify the ground and excited states for systems with various numbers of electrons using the same Fock-space representations of the downfolded Hamiltonian. The RRST formalism serves as a design principle for developing new classes of approximate schemes that reduce the complexity of quantum circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Dispersive and nondispersive 𝐾-matrix formalisms

The modeling of coupled-channel effects has become increasingly important due to the availability of highly precise data for a large variety of hadronic (re)scattering processes. The 𝐾-matrix is a powerful, yet comparatively simple, method to describe scattering amplitudes, including coupled-channel effects, with the aim of interpreting experimental data. Throughout the literature, a range of dispersive and nondispersive 𝐾-matrix methods are employed. Here, we compare the dispersive and nondispersive formulations in the context of the N/D method. It is shown that the methods are equivalent in the physical region under 𝐾-matrix reparametrization. Differences away from the physical region are examined. Applications to synthetic data are used to illustrate the effects of model choices concerning form factors and the application of dispersion relations, with the goal of clarifying best practices. We find no clear preference with regard to dispersive modeling. In contrast, we find that interpretational ambiguity of the bare model parameters—and even of the form of the bare model—is endemic, and recommend a thorough sampling of data and model spaces to assess conclusion robustness.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1988

The ijk forms of factorization methods. I - Vector computers

This paper gives a detailed exposition of the 'ijk forms' of LU and Choleski factorization. Several aspects of these different organizations are discussed and their properties on vector computers are compared. Extensions of the ijk formalism to other algorithms is also given.

Ortega, J. M.↗
DOE OSTI2023

Uncovering a CF 3 Effect on X‐ray Absorption Energies of [Cu(CF 3 ) 4 ] − and Related Copper Compounds by Using Resonant Diffraction Anomalous Fine Structure (DAFS) Measurements**

Abstract Understanding the electronic structures of high‐valent metal complexes aids the advancement of metal‐catalyzed cross coupling methodologies. A prototypical complex with formally high valency is [Cu(CF 3 ) 4 ] − (1), which has a formal Cu(III) oxidation state but whose physical analysis has led some to a Cu(I) assignment in an inverted ligand field model. Recent examinations of1by X‐ray spectroscopies have led previous authors to contradictory conclusions, motivating the re‐examination of its X‐ray absorption profile here by a complementary method, resonant diffraction anomalous fine structure (DAFS). From analysis of DAFS measurements for a series of seven mononuclear Cu complexes including1, here it is shown that there is a systematic trifluoromethyl effect on X‐ray absorption that blue shifts the resonant Cu K‐edge energy by 2–3 eV per CF 3 , completely accounting for observed changes in DAFS profiles between formally Cu(III) complexes like1and formally Cu(I) complexes like (Ph 3 P) 3 CuCF 3 (3). Thus, in agreement with the inverted ligand field model, the data presented herein imply that1is best described as containing a Cu(I) ion with d n count approaching 10.

Chemistry↗
DOE OSTI2022

Lattice Monte Carlo simulations with two impurity worldlines

Abstract We develop the impurity lattice Monte Carlo formalism for the case of two distinguishable impurities in a bath of polarized fermions. The majority particles are treated as explicit degrees of freedom, while the impurities are described by worldlines. The latter serve as localized auxiliary fields, which affect the majority particles. We apply the method to non-relativistic three-dimensional systems of two impurities and a number of majority particles where both the impurity–impurity interaction and the impurity–majority interaction have zero range. We consider the case of an attractive impurity–majority interaction, and we study the formation and disintegration of bound states as a function of the impurity–impurity interaction strength. We also discuss the potential applications of this formalism to other quantum many-body systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Communication: Variation after response in quantum Monte Carlo

We present a new method for modeling electronically excited states that overcomes a key failing of linear response theory by allowing the underlying ground state ansatz to relax in the presence of an excitation. The method is variational, has a cost similar to ground state variational Monte Carlo, and admits both open and periodic boundary conditions. We present preliminary numerical results showing that, when paired with the Jastrow antisymmetric geminal power ansatz, the variation-after-response formalism delivers accuracies for valence and charge transfer single excitations on par with equation of motion coupled cluster, while surpassing coupled cluster’s accuracy for excitations with significant doubly excited character.

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

Optical and magneto-optical properties of ferromagnetic monolayer CrBr 3 : A first-principles G W and G W plus Bethe-Salpeter equation study

The discovery of atomically thin two-dimensional (2D) magnetic semiconductors has triggered enormous research interest recently. In this paper, we use first-principles many-body perturbation theory to study a prototypical 2D ferromagnetic semiconductor, monolayer chromium tribromide (CrBr 3 ). With broken time-reversal symmetry, spin-orbit coupling, and excitonic effects included through the full-spinor GW and GW plus Bethe-Salpeter equation (GW-BSE) methods, we compute the frequency-dependent layer polarizability tensor and dielectric function tensor that govern the optical and magneto-optical (MO) properties. In addition, we provide a detailed theoretical formalism for simulating magnetic circular dichroism, MO Kerr effect, and Faraday effect, demonstrating the approach with monolayer CrBr 3 . Due to reduced dielectric screening in 2D and the localized nature of the Cr d orbitals, we find strong self-energy effects on the quasiparticle band structure of monolayer CrBr 3 that give a 3.8 eV indirect bandgap. Also, excitonic effects dominate the low-energy optical and MO responses in monolayer CrBr 3 where a large exciton binding energy of 2.3 eV is found for the lowest bright exciton state with excitation energy at 1.5 eV. We further find that the MO signals demonstrate strong dependence on the excitation frequency and substrate refractive index. Furthermore, our theoretical framework for modeling optical and MO effects could serve as a powerful theoretical tool for future study of optoelectronic and spintronics devices consisting of van der Waals 2D magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gravitational scattering and beyond from extreme mass ratio effective field theory

We explore a recently proposed effective field theory describing electromagnetically or gravitationally interacting massive particles in an expansion about their mass ratio, also known as the self-force (SF) expansion. By integrating out the deviation of the heavy particle about its inertial trajectory, we obtain an effective action whose only degrees of freedom are the lighter particle together with the photon or graviton, all propagating in a Coulomb or Schwarzschild background. The 0SF dynamics are described by the usual background field method, which at 1SF is supplemented by a “recoil operator” that encodes the wobble of the heavy particle, and similarly computable corrections appearing at 2SF and higher. Our formalism exploits the fact that the analytic expressions for classical backgrounds and particle trajectories encode dynamical information to all orders in the couplings, and from them we extract multiloop integrands for perturbative scattering. As a check, we study the two-loop classical scattering of scalar particles in electromagnetism and gravity, verifying known results. We then present new calculations for the two-loop classical scattering of dyons, and of particles interacting with an additional scalar or vector field coupling directly to the lighter particle but only gravitationally to the heavier particle.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Single-cone finite difference scheme for the (2+1)D Dirac von Neumann equation

An explicit finite difference scheme is presented for the von Neumann equation for (2+1)D Dirac fermions. It is founded upon a staggered space–time grid which ensures a single-cone energy dispersion and performs the time-derivative in one sweep using a three-step leap-frog procedure. It enables a space–time-resolved numerical treatment of the mixed-state dynamics of Dirac fermions within the effective single-particle density matrix formalism. Energy-momentum dispersion, stability and convergence properties are derived. Elementary numerical tests to demonstrate stability properties use parameters which pertain to topological insulator surface states. A method for the simulation of charge injection from an electric contact is presented and tested numerically. Potential extensions of the scheme to a Dirac–Lindblad equation, real-space–time Green's function formulations, and higher-order finite-difference schemes are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

SU-E-T-169: Evaluation of Oncentra TPS for Nasopharynx Brachy Using Patient Specific Voxel Phantom and EGSnrc

Purpose: Nasopharnx carcinoma (NPC) treatment is being carried out using Ir-192 HDR seeds in Mehdieh Hospital in Hamadan, Iran. The Oncentra™ TPS is based on optimized TG-43 formalism which disregards heterogeneity in the treatment area. Due to abundant heterogeneity in head and neck, comparison of the Oncentra™ TPS dose evaluation and an accurate dose calculation method in NPC brachytherapy is the objective of this study. Methods: CT DICOMs of a patient with NPC obtained from Mehdieh Hospital used to create 3D voxel phantom with CTCREATE utility of EGSnrc code package. The voxel phantom together with Ir-192 HDR brachytherapy source were the input to DOSXYZnrc to calculate the 3D dose distribution. The sources were incorporate with type 6 source in DOSXYZnrc and their dwell times were taken into account in final dose calculations. Results: The direct comparison between isodoses as well as DVHs for the GTV, PTV and CTV obtained by Oncentra™ and EGSnrc Monte Carlo code are made. EGSnrc results are obtained using 5×10{sup 9} histories to reduce the statistical error below 1% in GTV and 5% in 5% dose areas. The standard ICRP700 cross section library is employed in DOSXYZnrc dose calculation. Conclusion: A direct relationship between increased dose differences and increased material density (hence heterogeneity) is observed when isodoses contours of the TPS and DOSXYZnrc are compared. Regarding the point dose calculations, the differences range from 1.2% in PTV to 5.6% for cavity region and 7.8% for bone regions. While Oncentra™ TPS overestimates the dose in cavities, it tends to underestimate dose depositions within bones.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-201-10: Experimental Determination of Modified TG-43 Dosimetry Parameters for the Xoft Axxent® Electronic Brachytherapy Source

Purpose: The establishment of an air kerma rate standard at NIST for the Xoft Axxent{sup ®} electronic brachytherapy source (Axxent{sup ®} source) motivated the establishment of a modified TG-43 dosimetry formalism. This work measures the modified dosimetry parameters for the Axxent{sup ®} source in the absence of a treatment applicator for implementation in Xoft’s treatment planning system. Methods: The dose-rate conversion coefficient (DRCC), radial dose function (RDF) values, and polar anisotropy (PA) were measured using TLD-100 microcubes with NIST-calibrated sources. The DRCC and RDF measurements were performed in liquid water using an annulus of Virtual Water™ designed to align the TLDs at the height of the anode at fixed radii from the source. The PA was measured at several distances from the source in a PMMA phantom. MCNP-determined absorbed dose energy dependence correction factors were used to convert from dose to TLD to dose to liquid water for the DRCC, RDF, and PA measurements. The intrinsic energy dependence correction factor from the work of Pike was used. The AKR was determined using a NIST-calibrated HDR1000 Plus well-type ionization chamber. Results: The DRCC was determined to be 8.6 (cGy/hr)/(µGy/min). The radial dose values were determined to be 1.00 (1cm), 0.60 (2cm), 0.42 (3cm), and 0.32 (4cm), with agreement ranging from (5.7% to 10.9%) from the work of Hiatt et al. 2015 and agreement from (2.8% to 6.8%) with internal MCNP simulations. Conclusion: This work presents a complete dataset of modified TG-43 dosimetry parameters for the Axxent{sup ®} source in the absence of an applicator. Prior to this study a DRCC had not been measured for the Axxent{sup ®} source. This data will be used for calculating dose distributions for patients receiving treatment with the Axxent{sup ®} source in Xoft’s breast balloon and vaginal applicators, and for intraoperative radiotherapy. Sources and partial funding for this work were provided by Xoft Inc. (a subsidiary of iCAD). This work was also supported by the Radiological Sciences T32 Training Grant through the University of Wisconsin-Madison Medical Physics department (5T32CA009206-37).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

New charged-particle transport computational capability: the SIT code An L-4 milestone

We have developed a new high-fidelity code for direct transport of charged particles using the simple integral transport method. The code can be coupled to the outputs of any hydrodynamical simulation code in 1-D, 2-D or 3-D. In this report we summarize the formalism involved in treating complex transport problems. We present physical examples wherein we have used the code to calculate the transport of alpha particles. Future work is planned to study the sensitivity of hydrodynamical mix to charged-particle radiochemistry and reaction-in-flight neutrons for the complex inertial confinement fusion problems encountered at NIF and at the Z-machine.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1983

Distributed system modeling of a large space antenna

A general approach for distributed parameter modeling of complex dynamical systems is described. The method consists of dividing the system in parts which can be modeled by simple partial differential equations and coupling the equations thus obtained by applying Hamilton's variational formalism to the entire system. The modeling of a large, offset-fed, wrap-rib antenna is presented to illustrate the approach. Although such models are perhaps not as precise as finite element models, they can be useful for initial physical insight and parametric design.

Hamidi, M.↗
NASA NTRS1985

Contribution functions for Zeeman-split lines, and line formation in photospheric faculae

The transfer of polarized light in an inhomogeneous stellar atmosphere, and the formation of magnetically sensitive spectral lines, are discussed. A new method for the solution of the transfer equations is proposed. The method gives a natural definition of the contribution functions for Stokes' parameters, i.e., functions describing the contributions from different parts along the line-of-sight (LOS). The formalism includes all magneto-optical effects, and allows for an arbitrary variation of magnetic field, velocity field, temperature, density, etc., along the LOS. The formation of FeI lambda 5250.2 in photospheric faculae is described. A potential-field model of a facular element is presented, and spectra profiles and contribution functions are computed for the Stokes parameters I, Q, and V.

Vanballegooijen, A. A.↗