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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Autogenic-Feedback Training (AFT) as a preventive method for space motion sickness: Background and experimental design

Finding an effective treatment for the motion sickness-like symptoms that occur in space has become a high priority for NASA. The background research is reviewed and the experimental design of a formal life sciences shuttle flight experiment designed to prevent space motion sickness in shuttle crew members is presented. This experiment utilizes a behavioral medicine approach to solving this problem. This method, Autogenic-Feedback Training (AFT), involves training subjects to voluntarily control several of their own physiological responses to environmental stressors. AFT has been used reliably to increase tolerance to motion sickness during ground-based tests in over 200 men and women under a variety of conditions that induce motion sickness, and preliminary evidence from space suggests that AFT may be an effective treatment for space motion sickness as well. Proposed changes to this experiment for future manifests are included.

Cowings, Patricia S.↗
DOE OSTI2015

Conditioning of Pu-containing radioactive waste generated in the Hot-lab: where are we now - 15400

To achieve an optimal strategy for interim storage and final disposal, solid radioactive waste has to be sorted and treated. This strategy is based on the requirement to reduce the volume of waste designated for long term disposal, mainly for economical and safety reasons. In addition, the radiological limitations for shipping of the waste packages to a disposal site have to be taken into account as well as the improved regulatory requirements for material under IAEA control. In case of the operational Pu-contaminated and -containing waste from the PSI-Hot laboratory in Switzerland, the existing and over years practiced conditioning method needed improvement. The existing authorized specifications of the conditioning process had to be revised for an optimization of the interim-storage and the deep geological disposal. This was done by the PSI Section for Dismantling and Waste Management. The improvement focused on technical aspects as well as improvement of the documentation for the certification by final depositor (Nagra) and the safety authority (ENSI). This work started in 2008. In 2012 ENSI issued an authorization to realize the first conditioning according the new specification, formally as a certification test ('Type Test') prior to final approval by the ENSI. But, although a licensed method by the competent authority is available, modifications could be necessary and new tasks come up. In 2013 the PSI has produced 13 waste-packages as a Type Test. To obtain the legal authorization ('Type Approval'), some open issues ought to be answered. Waste is continuously generated and so the question arises: When will the next conditioning be done and applying which method? (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Voltage Restoration in Microgrids using Temporal Logic Specifications

This paper proposes an energy storage controller synthesis method for voltage restoration in microgrids with respect to temporal logic specifications (TLSs). TLSs is introduced in this paper as a formalism to control the voltage variation of a critical bus against an operational bounds over time. The power system with synchronous generator (SG) connected to a critical load bus is modeled as a set of differential-algebraic equations and a simplified analytical model is derived to describe the voltage variation of this critical bus. The control objective is to schedule an optimal control input signal from a supportive energy storage system (ESS) connected to the critical bus, such that the voltage variation of the latter satisfies the TLSs, such as a finite-time restoration. The proposed control is verified on a lumped distribution system model. With this control diagram, supportive controllers can be designed to make voltage behaviors comply with grid codes and avoid unnecessary relay actions.

Taousser, Fatima Zohra↗
DOE OSTI2023

Loop-by-loop differential equations for dual (elliptic) Feynman integrals

We present a loop-by-loop method for computing the differential equations of Feynman integrals using the recently developed dual form formalism. We give explicit prescriptions for the loop-by-loop fibration of multi-loop dual forms. Then, we test our formalism on a simple, but non-trivial, example: the two-loop three-mass elliptic sunrise family of integrals. We obtain an ε-form differential equation within the correct function space in a sequence of relatively simple algebraic steps. In particular, none of these steps relies on the analysis of q-series. Then, we discuss interesting properties satisfied by our dual basis as well as its simple relation to the known ε-form basis of Feynman integrands. The underlying K3-geometry of the three-loop four-mass sunrise integral is also discussed. Finally, we speculate on how to construct a “good” loop-by-loop basis at three-loop.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Moment expansion of polarized dust SED: A new path towards capturing the CMB B -modes with LiteBIRD

Accurate characterization of the polarized dust emission from our Galaxy will be decisive in the quest for the cosmic microwave background (CMB) primordial B-modes. An incomplete modeling of its potentially complex spectral properties could lead to biases in the CMB polarization analyses and to a spurious measurement of the tensor-to-scalar ratio r. It is particularly crucial for future surveys like the LiteBIRD satellite, the goal of which is to constrain the faint primordial signal leftover by inflation with an accuracy on the tensor-to-scalar ratio r of the order of 10 –3 . Variations of the dust properties along and between lines of sight lead to unavoidable distortions of the spectral energy distribution (SED) that cannot be easily anticipated by standard component-separation methods. This issue can be tackled using a moment expansion of the dust SED, an innovative parametrization method imposing minimal assumptions on the sky complexity. In the present paper, we apply this formalism to the B-mode cross-angular power spectra computed from simulated LiteBIRD polarization data at frequencies between 100 and 402 GHz that contain CMB, dust, and instrumental noise. The spatial variation of the dust spectral parameters (spectral index β and temperature T) in our simulations lead to significant biases on r (~21 σ r ) if not properly taken into account. Performing the moment expansion in β, as in previous studies, reduces the bias but does not lead to sufficiently reliable estimates of r. We introduce, for the first time, the expansion of the cross-angular power spectra SED in both β and T, showing that, at the sensitivity of LiteBIRD, the SED complexity due to temperature variations needs to be taken into account in order to prevent analysis biases on r. Thanks to this expansion, and despite the existing correlations between some of the dust moments and the CMB signal responsible for a rise in the error on r, we can measure an unbiased value of the tensor-to-scalar ratio with a dispersion as low as σ r = 8.8 × 10 –4 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Far-infrared collision-induced absorption in rare gas mixtures: Quantum and semi-classical calculations

We compare calculations of the translational collision-induced spectra and their integrated intensities of both He–Ar and Ne–Ar collisional complexes, using the quantum mechanical and a semiclassical formalism. Advanced potential energy and induced dipole functions are used for the calculations. The quantum method used is as described previously [L. Frommhold, Collision-induced Absorption in Gases (Cambridge University Press, 1993 and 2006)]. The semiclassical method is based on repeated classical atom-atom scattering calculations to simulate an ensemble average; subsequent Fourier transform then renders the binary absorption coefficient as a function of frequency. The problem of classical calculations is the violation of the principle of detailed balance, which may be introduced only artificially in classical calculations. Nevertheless, it is shown that the use of classical trajectories permits a fairly accurate reproduction of the experimental spectra, comparable to the quantum mechanical results at not too low temperatures and for collisional pairs of not too small reduced mass. Inexpensive classical calculations may thus be promising to compute spectra also of molecular pairs, or even of polyatomic collisional pairs with anisotropic intermolecular interactions, for which the quantum approach is still inefficient or impractical.

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

Dielectric response and electric modulus of Y{sub 2}CrCoO{sub 6} perovskite

Y{sub 2}CrCoO{sub 6} perovskite prepared by solution combustion method and sintered at 1073 K has been characterized by dielectric spectroscopy and electric modulus formalism. Temperature and frequency dependent measurements of permitivitty reveal that observed relaxation might be related to the hopping conductivity, i.e. universal dielectric response. The presence of electrode polarization is dominant at low frequencies. The electric modulus dependencies enable us to distinguish and separate the relaxation processes connected to the conduction processes in the material. The presences of both grain and grain boundary effects are established, each dominant in different frequency and temperature range. The conductivity through grain boundaries obeys metalic behavior, while conductivity through grains shows semiconductor behavior. The electrical behavior of this material depends on the differences in (Cr-O) and (Co-O) bond lenghts, Co{sup 3+} being in the low-spin state, resulting in shorter Co-O and thus stronger π bonding e.g. more efficient overlapping of the Co{sup 3+} d-orbitals with oxygen p{sub π} orbitals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

DL_POLY Quantum 2.0: A modular general-purpose software for advanced path integral simulations

DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman’s path integral formalism of quantum mechanics. Here, the flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Thermodynamic assessment of the Pr-O system

We found that the Calphad method was used to perform a thermodynamic assessment of the Pr–O system. Compound energy formalism representations were developed for the fluorite α-PrO 2–x and bixbyite σ-Pr 3 O 5 ± x solid solutions while the two-sublattice liquid model was used to describe the binary melt. The series of phases between Pr 2 O 3 and PrO 2 were taken to be stoichiometric. Moreover, the equilibrium oxygen pressure, phase equilibria, and enthalpy data were used to optimize the adjustable parameters of the models for a self-consistent representation of the thermodynamic behavior of the Pr–O system from 298 K to melting.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

QECC-Synth: A Layout Synthesizer for Quantum Error Correction Codes on Sparse Architectures

Quantum Error Correction (QEC) codes are essential for achieving fault-tolerant quantum computing (FTQC). However, their implementation faces significant challenges due to disparity between required dense qubit connectivity and sparse hardware architectures. Current approaches often either underutilize QEC circuit features or focus on manual designs tailored to specific codes and architectures, limiting their capability and generality. In response, we introduce QECC-Synth, an automated compiler for QEC code implementation that addresses these challenges. We leverage the ancilla bridge technique tailored to the requirements of QEC circuits and introduces a systematic classification of its design space flexibilities. We then formalize this problem using the MaxSAT framework to optimize these flexibilities. Evaluation shows that our method significantly outperforms existing methods while demonstrating broader applicability across diverse QEC codes and hardware architectures.

Yin, Keyi [University of California, San Diego]↗
NASA NTRS1984

Water vapor profile retrieval simulation studies for the HIRS/MSU and AMTS/MSU sounders

A combined physical statistical method for the retrieval of atmospheric water vapor column densities has been developed at the Goddard Laboratory for Atmospheric Sciences. In this method differences between observed and calculated clear column brightness temperatures are used in a regression formalism to obtain estimates of integrated column densities in given atmospheric layers. The physics of the radiative transfer problem is contained in the calculated clear column brightness temperatures and in accounting for the effects of clouds on the observed radiances. The regression matrix is equivalent to a set of empirical sensitivity factors.

Reuter, D.↗
NASA NTRS1987

SHARP simulation of discontinuities in highly convective steady flow

For steady multidimesional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗
NASA NTRS1988

Simple high-accuracy resolution program for convective modelling of discontinuities

For steady multidimensional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗
NASA NTRS2008

Safer Systems: A NextGen Aviation Safety Strategic Goal

The Joint Planning and Development Office (JPDO), is charged by Congress with developing the concepts and plans for the Next Generation Air Transportation System (NextGen). The National Aviation Safety Strategic Plan (NASSP), developed by the Safety Working Group of the JPDO, focuses on establishing the goals, objectives, and strategies needed to realize the safety objectives of the NextGen Integrated Plan. The three goal areas of the NASSP are Safer Practices, Safer Systems, and Safer Worldwide. Safer Practices emphasizes an integrated, systematic approach to safety risk management through implementation of formalized Safety Management Systems (SMS) that incorporate safety data analysis processes, and the enhancement of methods for ensuring safety is an inherent characteristic of NextGen. Safer Systems emphasizes implementation of safety-enhancing technologies, which will improve safety for human-centered interfaces and enhance the safety of airborne and ground-based systems. Safer Worldwide encourages coordinating the adoption of the safer practices and safer systems technologies, policies and procedures worldwide, such that the maximum level of safety is achieved across air transportation system boundaries. This paper introduces the NASSP and its development, and focuses on the Safer Systems elements of the NASSP, which incorporates three objectives for NextGen systems: 1) provide risk reducing system interfaces, 2) provide safety enhancements for airborne systems, and 3) provide safety enhancements for ground-based systems. The goal of this paper is to expose avionics and air traffic management system developers to NASSP objectives and Safer Systems strategies.

Darr, Stephen T.↗
DOE OSTI2019

Radioactive Decay Computation with a Graph Approach for Solving a Linear Dynamic System

ABSTRACT: This article gives two graph-based algorithms for radioactive decay computation. The first algorithm identifies the connected components of the graph induced from the given radioactive decay dynamics to reduce the size of the problem. The solutions are derived over the precalculated connected components, respectively and independently. The second algorithm utilizes acyclic structure of radioactive decay dynamics. The algorithm evaluates the reachable vertices of the induced system graph from the initially activated vertices and finds the minimal set of starting vertices populating the entire reachable vertices. Then, the decay calculations are performed over the reachable vertices from the identified minimal starting vertices, respectively, with the partitioned initial value over the reachable vertices. A formal argument is given to show that the proposed graph inspired divide and conquer calculation methods perform the intended radioactive decay calculation. Empirical efforts comparing the proposed radioactive decay calculation algorithms are presented.

42 ENGINEERING↗
DOE OSTI2024

On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH

We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1999

A Spherical Earth Solution for TOA Lightning Location Retrieval

With the advent of high-speed digital computers, the method of chi square minimization is a highly practical means for analyzing a wide variety of (otherwise intractable) nonlinear inversion problems in applied mathematical physics. Little thought or effort is required to apply the chi square method to obtain quick and reasonable estimates of a solution, and the method offers a means to assess retrieval errors. Because the method is simple and practical it is sometimes hastily applied to problems that can be solved by formal analytic or quasi-analytic means. Presently, Global Atmospherics Inc. (GAI) finds the minimum of a chi square function to analyze time-of-arrival (TOA) and magnetic bearing data derived from the National Lightning Detection Network (NLDN); ellipsoidal Earth geometry is assumed. An analytic solution to this problem has not yet been attained, but the consideration and solving of less general problem statements might eventually lead to a final solution. In the present study, the problem of retrieving lightning ground-strike location on a spherical Earth surface using a network of 4 or more time-of-arrival (TOA) sensors is considered. It is shown that this problem has an analytic solution and therefore does not require the use of nonlinear estimation theory (such as the chi square method mentioned above). The mathematical robustness of the analytic solution is tested using computer-generated lightning sources and simulated TOA measurement errors. A quasi-analytic extension of the spherical Earth solution for an oblate spheroidal Earth geometry is considered in a related study. The incorporation of magnetic bearing information into these analytic solutions would lead to a general and elegant analytic retrieval scheme that would most likely replace the chi square estimation theory currently employed by Global Atmospherics Inc. (GAI).

Koshak, William J,↗
DOE OSTI2023

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗