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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-F-18C-04: A Combination of Monoenergetic Reconstruction and Stoichiometric Calibration for Tissue Characterization Using Dual Energy Computed Tomography

Purpose: Dual energy computed tomography (DECT) pre-reconstruction methods require the prior knowledge of the X-ray source spectrum to allow extracting physical parameters needed for radiation therapy dose calculation, such as electron density (ED) and the effective atomic number (EAN). While DECT stoichiometric calibration may provide reliable performance for typical radiation therapy clinical conditions, it is yet to be adapted to prereconstruction methods. The presence of noise and inaccurate spectrum description may lead to systematic errors and artifacts which compromise the accuracy of treatment planning. Methods: A new technique is investigated which consists in applying a DECT stoichiometric calibration method to a set of monoenergetic images obtained with a DECT prereconstruction method. To evaluate the performance of this extended method, a simulation environment is developed to generate DECT scans under well controlled conditions, to reconstruct monoenergetic images of a tissue-equivalent phantom from transformed sinograms and to extract ED and EAN maps using a DECT formalism. Result: Under simulated clinical conditions, the accuracy in determining ED with the extended method versus a pre-reconstruction method alone is shown to be better than 0.35% versus 0.5%, respectively. In the presence of a realistic noise level, EAN determination presents a relative mean error that drops from 2.5% to 0.5% once the calibration is applied. Considering a spectrum alteration by a 1 mm Cu layer, EAN errors are up to 30% for the pre-reconstruction method alone versus less than 3% for the extended method. Conclusion: This study shows that combining pre-reconstruction DECT methods with a stoichiometric calibration considerably improves the accuracy and reliability of tissue characterization for radiation therapy in a clinical context. The presented method could potentially be adapted as gold standard for dose calculation methods based on DECT.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2018

An Eigenvector-based Method of Radio Array Calibration and Its Application to the Tianlai Cylinder Pathfinder

We propose an eigenvector-based formalism for the calibration of radio interferometer arrays. In the presence of a strong dominant point source, the complex gains of the array can be obtained by taking the first eigenvector of the visibility matrix. We use the stable principle component analysis (SPCA) method to help separate outliers and noise from the calibrator signal to improve the performance of the method. This method can be applied with poorly known beam model of the antenna, and is insensitive to outliers or imperfections in the data, and has low computational complexity. It thus is particularly suitable for the initial calibration of the array, which can serve as the initial point for more accurate calibrations. We demonstrate this method by applying it to the cylinder pathfinder of the Tianlai experiment, which aims to measure the dark energy equation of state using the baryon acoustic oscillation (BAO) features in the large scale structure by making intensity mapping observation of the redshifted 21~cm emission of the neutral hydrogen (HI). The complex gain of the array elements and the beam profile in the East-West direction (short axis of the cylinder) are successfully obtained by applying this method to the transit data of bright radio sources.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Huzinaga projection embedding for efficient and accurate energies of systems with localized spin-densities

We demonstrate the accuracy and efficiency of the restricted open-shell and unrestricted formulation of the absolutely localized Huzinaga projection operator embedding method. Restricted open-shell and unrestricted Huzinaga projection embedding in the full system basis is formally exact to restricted open-shell and unrestricted Kohn–Sham density functional theory, respectively. By utilizing the absolutely localized basis, we significantly improve the efficiency of the method while maintaining high accuracy. Furthermore, the absolutely localized basis allows for high accuracy open-shell wave function methods to be embedded into a closed-shell density functional theory environment. The open-shell embedding method is shown to calculate electronic energies of a variety of systems to within 1 kcal/mol accuracy of the full system wave function result. For certain highly localized reactions, such as spin transition energies on transition metals, we find that very few atoms are necessary to include in the wave function region in order to achieve the desired accuracy. This extension further broadens the applicability of our absolutely localized Huzinaga level-shift projection operator method to include open-shell species. Here, we apply our method to several representative examples, such as spin splitting energies, catalysis on transition metals, and radical reactions.

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

Solving the nuclear pairing model with neural network quantum states

In this work, we present a variational Monte Carlo method that solves the nuclear many-body problem in the occupation number formalism exploiting an artificial neural network representation of the groundstate wave function. A memory-efficient version of the stochastic reconfiguration algorithm is developed to train the network by minimizing the expectation value of the Hamiltonian. We benchmark this approach against widely used nuclear many-body methods by solving a model used to describe pairing in nuclei for different types of interaction and different values of the interaction strength. Despite its polynomial computational cost, our method outperforms coupled-cluster and provides energies that are in excellent agreement with the numerically-exact full configuration interaction values.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2010

Cosmic microwave background with Brans-Dicke gravity. I. Covariant formulation

In the covariant cosmological perturbation theory, a 1+3 decomposition ensures that all variables in the frame-independent equations are covariant and gauge invariant and that they have clear physical interpretations. We develop this formalism in the case of Brans-Dicke gravity, and apply this method to the calculation of CMB anisotropy and large scale structures. We modify the publicly available Boltzmann code CAMB to calculate numerically the evolution of the background and adiabatic perturbations, and obtain the temperature and polarization spectra of Brans-Dicke theory for both scalar and tensor modes; the tensor mode results for Brans-Dicke gravity are obtained numerically for the first time. We first present our theoretical formalism in detail, and then explicitly describe the techniques used in modifying the CAMB code. These techniques are also very useful for other gravity models. Next we compare the CMB and large scale structure spectra in Brans-Dicke theory with those in the standard general relativity theory. At last, we investigate the integrated Sachs-Wolfe effect and the CMB lensing effect in Brans-Dicke theory. Constraints on the Brans-Dicke model with current observational data are presented in a companion paper [F. Wu and X. Chen, following Article, Phys. Rev. D 82, 083003 (2010)] (paper II).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Mapping valence electron distributions with multipole density formalism using 4D-STEM

Recent advancement in aberration correction and detector technology opened a door to various applications using 4D-STEM, which yields a diffraction pattern for each scanning position within a crystal unit-cell in scanning transmission electron microscopy (STEM) and generates incredible amounts of data in momentum space. Currently 4D-STEM analysis relies on the center-of-mass of the diffraction patterns in electric field and charge density mapping. It only derives the total projected charge density and is limited to phase objects, e.g. extremely thin samples. Here, we propose a new analytical method to accurately map aspherical valence electron distributions with atom-centered multipolar functions formalism using the whole 4D-STEM dataset. We demonstrate that, with the full dynamical calculations for various sample thicknesses, the method is sensitive not only to the miniscule charge transfer, but also to the atomic site symmetry and aspherical electron orbitals. The process of the refinement is much more robust and reliable than quantitative convergent beam electron diffraction.

4D-STEM↗
DOE OSTI2012

Two-Component Relativistic Calculations of Electric-Field Gradients Using Exact Decoupling Methods: Spin–orbit and Picture-Change Effects

Electric field gradients (EFGs) for heavy nuclei in diatomic molecules (hydrogen halides and group-13 iodides), uranyl (UO 2 2+ ), and a uranyl carbonate complex, were computed with different formally exactly (or nearly exactly) decoupled two-component relativistic Hamiltonians: the “exact two-component” (X2C) method, the Barysz–Sadlej–Snijders (BSS) approach, and up to 35th order in the Douglas–Kroll–Hess (DKH) expansion, utilizing a new implementation in the open-source NWChem quantum chemistry package. Results from two-component Hartree–Fock and density functional theory (DFT) calculations at the scalar relativistic approximation as well as including spin–orbit coupling are reported. Picture-change corrected EFGs obtained with X2C, BSS, and high-order DKH are shown to be numerically equivalent. Effects from spin–orbit coupling on the EFGs tend to be moderate but should not be excluded for reliable predictions. Picture-change effects on the EFG of a heavy atom can be as large as the correct picture-change corrected value, as is shown for the uranium EFG in the uranyl systems. The results demonstrate the necessity of simultaneously including spin–orbit and picture-change effects in calculations of heavy atom EFGs.

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

Spin-flip configuration interaction singles with exact spin-projection: Theory and applications to strongly correlated systems

Spin-flip approaches capture static correlation with the same computational scaling as the ordinary single reference methods. Here, we extend spin-flip configuration interaction singles (SFCIS) by projecting out intrinsic spin-contamination to make it spin-complete, rather than by explicitly complementing it with spin-coupled configurations. We give a general formalism of spin-projection for SFCIS, applicable to any spin states. The proposed method is viewed as a natural unification of SFCIS and spin-projected CIS to achieve a better qualitative accuracy at a low computational cost. While our wave function ansatz is more compact than previously proposed spin-complete SF approaches, it successfully offers more general static correlation beyond biradicals without sacrificing good quantum numbers. It is also shown that our method is invariant with respect to open-shell orbital rotations, due to the uniqueness of spin-projection. We will report benchmark calculations to demonstrate its qualitative performance on strongly correlated systems, including conical intersections that appear both in ground-excited and excited-excited degeneracies.

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

Stability of compressible reduced magnetohydrodynamic equilibria-Analogy with magnetorotational instability

Stability analyses for equilibria of the compressible reduced magnetohydrodynamics (CRMHD) model are carried out by means of the Energy-Casimir (EC) method. Stability results are compared with those obtained for ideal magnetohydrodynamics (MHD) from the classical {delta}W criterion. An identification of the terms in the second variation of the free energy functional for CRMHD with those of {delta}W is made: two destabilizing effects present for CRMHD turn out to correspond to the kink and interchange instabilities in usual MHD, while the stabilizing roles of field line bending and compressibility are also identified in the reduced model. Also, using the EC method, stability conditions in the presence of toroidal flow are obtained. A formal analogy between CRMHD and a reduced incompressible model for magnetized rotating disks, due to Julien and Knobloch [EAS Pub. Series, 21, 81 (2006)], is discovered. In light of this analogy, energy stability analysis shows that the condition for magnetorotational instability (MRI) for the latter model corresponds to the condition for interchange instability in CRMHD, with the Coriolis term and shear velocity playing the roles of the curvature term and pressure gradient, respectively. Using the EC method, stability conditions for the rotating disk model, for a large class of equilibria with possible non-uniform magnetic fields, are obtained. In particular, this shows it is possible for the MRI system to undergo, in addition to the MRI, another instability that is analogous to the kink instability. For vanishing magnetic field, the Rayleigh hydrodynamical stability condition is recovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Higher-order temporal integration for the incompressible Navier–Stokes equations in bounded domains

Here, this paper compares and contrasts higher-order, semi-implicit temporal integration strategies for the incompressible Navier–Stokes methods based on spectral deferred corrections applied to certain gauge or auxiliary variable formulations of the equations. Particular focus is placed on the imposition of boundary conditions in the semi-implicit formulation, the accurate treatment of the pressure term, and the smoothness of the numerical solution for different formulations. The main result presented here is the formulation and numerical validation of a single-step, semi-implicit method called spectral deferred pressure corrections, which can in theory obtain arbitrary formal order of accuracy in time. Lastly, numerical results demonstrate up to eighth-order accuracy, scenarios where optimal temporal accuracy is attained, and scenarios where order reduction is observed due to time-dependent boundary conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Three-body resonances in the φ 4 theory

We study the properties of three-body resonances using a lattice complex scalar φ4 theory with two scalars, with parameters chosen such that one heavy particle can decay into three light ones. We determine the two- and three-body spectra for several lattice volumes using variational techniques, and then analyze them with two versions of the three-particle finite-volume formalism: the Relativistic Field Theory approach and the Finite-Volume Unitarity approach. We find that both methods provide an equivalent description of the energy levels, and we are able to fit the spectra using simple parametrizations of the scattering quantities. By solving the integral equations of the corresponding three-particle formalisms, we determine the pole position of the resonance in the complex energy plane and thereby its mass and width. We find very good agreement between the two methods at different values of the coupling of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2003

Single-Vector Calibration of Wind-Tunnel Force Balances

An improved method of calibrating a wind-tunnel force balance involves the use of a unique load application system integrated with formal experimental design methodology. The Single-Vector Force Balance Calibration System (SVS) overcomes the productivity and accuracy limitations of prior calibration methods. A force balance is a complex structural spring element instrumented with strain gauges for measuring three orthogonal components of aerodynamic force (normal, axial, and side force) and three orthogonal components of aerodynamic torque (rolling, pitching, and yawing moments). Force balances remain as the state-of-the-art instrument that provide these measurements on a scale model of an aircraft during wind tunnel testing. Ideally, each electrical channel of the balance would respond only to its respective component of load, and it would have no response to other components of load. This is not entirely possible even though balance designs are optimized to minimize these undesirable interaction effects. Ultimately, a calibration experiment is performed to obtain the necessary data to generate a mathematical model and determine the force measurement accuracy. In order to set the independent variables of applied load for the calibration 24 NASA Tech Briefs, October 2003 experiment, a high-precision mechanical system is required. Manual deadweight systems have been in use at Langley Research Center (LaRC) since the 1940s. These simple methodologies produce high confidence results, but the process is mechanically complex and labor-intensive, requiring three to four weeks to complete. Over the past decade, automated balance calibration systems have been developed. In general, these systems were designed to automate the tedious manual calibration process resulting in an even more complex system which deteriorates load application quality. The current calibration approach relies on a one-factor-at-a-time (OFAT) methodology, where each independent variable is incremented individually throughout its full-scale range, while all other variables are held at a constant magnitude. This OFAT approach has been widely accepted because of its inherent simplicity and intuitive appeal to the balance engineer. LaRC has been conducting research in a "modern design of experiments" (MDOE) approach to force balance calibration. Formal experimental design techniques provide an integrated view to the entire calibration process covering all three major aspects of an experiment; the design of the experiment, the execution of the experiment, and the statistical analyses of the data. In order to overcome the weaknesses in the available mechanical systems and to apply formal experimental techniques, a new mechanical system was required. The SVS enables the complete calibration of a six-component force balance with a series of single force vectors.

Parker, P. A.↗
DOE OSTI2013

Complex-scaled equation-of-motion coupled-cluster method with single and double substitutions for autoionizing excited states: Theory, implementation, and examples

Theory and implementation of complex-scaled variant of equation-of-motion coupled-cluster method for excitation energies with single and double substitutions (EOM-EE-CCSD) is presented. The complex-scaling formalism extends the EOM-EE-CCSD model to resonance states, i.e., excited states that are metastable with respect to electron ejection. The method is applied to Feshbach resonances in atomic systems (He, H{sup -}, and Be). The dependence of the results on one-electron basis set is quantified and analyzed. Energy decomposition and wave function analysis reveal that the origin of the dependence is in electron correlation, which is essential for the lifetime of Feshbach resonances. It is found that one-electron basis should be sufficiently flexible to describe radial and angular electron correlation in a balanced fashion and at different values of the scaling parameter, {theta}. Standard basis sets that are optimized for not-complex-scaled calculations ({theta} = 0) are not sufficiently flexible to describe the {theta}-dependence of the wave functions even when heavily augmented by additional sets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Note on two formulations of Crank-Nicolson method for Navier-Stokes equations

Here, we consider two formulations of the Crank-Nicolson (CN) method for the Navier-Stokes equations (NSE). The “natural” way of implementing CN for NSE is formally second order accurate in time for both velocity and pressure, whereas another formulation approximates pressure with only first order accuracy in time. Both versions of the method are applied to the benchmark problem of computing drag and lift in the flow around a cylinder. We show that the presumably more accurate version of the CN can create a solution with nonphysical oscillations and give incorrect predictions for the maximal drag coefficient, whereas the other formulation of the method predicts the drag and lift coefficients more accurately and does not introduce nonphysical oscillations. We locate the source of the issue and suggest several remedies.

Crank-Nicolson↗
DOE OSTI2013

Symbolic methods for the evaluation of sum rules of Bessel functions

The use of the umbral formalism allows a significant simplification of the derivation of sum rules involving products of special functions and polynomials. We rederive in this way known sum rules and addition theorems for Bessel functions. Furthermore, we obtain a set of new closed form sum rules involving various special polynomials and Bessel functions. The examples we consider are relevant for applications ranging from plasma physics to quantum optics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

Hamilton/Jacobi perturbation methods applied to the rotational motion of a rigid body in a gravitational field

The formalism for studying perturbations of a triaxial rigid body within the Hamilton-Jacobi framework is developed. The motion of a triaxial artificial earth satellite about its center of mass is studied. Variables are found which permit separation, and the Euler angles and associated conjugate momenta are obtained as functions of canonical constants and time.

Fitzpatrick, P. M.↗
DOE OSTI2015

Theoretically informed Monte Carlo simulation of liquid crystals by sampling of alignment-tensor fields

A theoretically informed coarse-grained Monte Carlo method is proposed for studying liquid crystals. The free energy functional of the system is described in the framework of the Landau-de Gennes formalism. The alignment field and its gradients are approximated by finite differences, and the free energy is minimized through a stochastic sampling technique. The validity of the proposed method is established by comparing the results of the proposed approach to those of traditional free energy minimization techniques. Its usefulness is illustrated in the context of three systems, namely, a nematic liquid crystal confined in a slit channel, a nematic liquid crystal droplet, and a chiral liquid crystal in the bulk. It is found that for systems that exhibit multiple metastable morphologies, the proposed Monte Carlo method is generally able to identify lower free energy states that are often missed by traditional approaches. Importantly, the Monte Carlo method identifies such states from random initial configurations, thereby obviating the need for educated initial guesses that can be difficult to formulate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

NWChem: Quantum Chemistry Simulations at Scale

Methods based on quantum mechanics equations have been developed since the 1930's with the purpose of accurately studying the electronic structure of molecules. However, it is only during the last two decades that intense development of new computational algorithms has opened the possibility of performing accurate simulations of challenging molecular processes with high-order many-body methods. A wealth of evidence indicates that the proper inclusion of instantaneous interactions between electrons (or the so-called electron correlation effects) is indispensable for the accurate characterization of chemical reactivity, molecular properties, and interactions of light with matter. The availability of reliable methods for benchmarking of medium-size molecular systems provides also a unique chance to propagate high-level accuracy across spatial scales through the multiscale methodologies. Some of these methods have potential to utilize computational resources in an effi*cient way since they are characterized by high numerical complexity and appropriate level of data granularity, which can be effi*ciently distributed over multi-processor architectures. The broad spectrum of coupled cluster (CC) methods falls into this class of methodologies. Several recent CC implementations clearly demonstrated the scalability of CC formalisms on architectures composed of hundreds thousand computational cores. In this context NWChem provides a collection of Tensor Contraction Engine (TCE) generated parallel implementations of various coupled cluster methods capable of taking advantage of many thousand of cores on leadership class parallel architectures.

chemistry↗