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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

ON IMPROVING ANALYTICAL MODELS OF COSMIC REIONIZATION FOR MATCHING NUMERICAL SIMULATION

The methods for studying the epoch of cosmic reionization vary from full radiative transfer simulations to purely analytical models. While numerical approaches are computationally expensive and are not suitable for generating many mock catalogs, analytical methods are based on assumptions and approximations. We explore the interconnection between both methods. First, we ask how the analytical framework of excursion set formalism can be used for statistical analysis of numerical simulations and visual representation of the morphology of ionization fronts. Second, we explore the methods of training the analytical model on a given numerical simulation. We present a new code which emerged from this study. Its main application is to match the analytical model with a numerical simulation. Then, it allows one to generate mock reionization catalogs with volumes exceeding the original simulation quickly and computationally inexpensively, meanwhile reproducing large-scale statistical properties. These mock catalogs are particularly useful for cosmic microwave background polarization and 21 cm experiments, where large volumes are required to simulate the observed signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The many-body Wigner Monte Carlo method for time-dependent ab-initio quantum simulations

The aim of ab-initio approaches is the simulation of many-body quantum systems from the first principles of quantum mechanics. These methods are traditionally based on the many-body Schrödinger equation which represents an incredible mathematical challenge. In this paper, we introduce the many-body Wigner Monte Carlo method in the context of distinguishable particles and in the absence of spin-dependent effects. Despite these restrictions, the method has several advantages. First of all, the Wigner formalism is intuitive, as it is based on the concept of a quasi-distribution function. Secondly, the Monte Carlo numerical approach allows scalability on parallel machines that is practically unachievable by means of other techniques based on finite difference or finite element methods. Finally, this method allows time-dependent ab-initio simulations of strongly correlated quantum systems. In order to validate our many-body Wigner Monte Carlo method, as a case study we simulate a relatively simple system consisting of two particles in several different situations. We first start from two non-interacting free Gaussian wave packets. We, then, proceed with the inclusion of an external potential barrier, and we conclude by simulating two entangled (i.e. correlated) particles. The results show how, in the case of negligible spin-dependent effects, the many-body Wigner Monte Carlo method provides an efficient and reliable tool to study the time-dependent evolution of quantum systems composed of distinguishable particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Sparse maps—A systematic infrastructure for reduced-scaling electronic structure methods. I. An efficient and simple linear scaling local MP2 method that uses an intermediate basis of pair natural orbitals

In this work, a systematic infrastructure is described that formalizes concepts implicit in previous work and greatly simplifies computer implementation of reduced-scaling electronic structure methods. The key concept is sparse representation of tensors using chains of sparse maps between two index sets. Sparse map representation can be viewed as a generalization of compressed sparse row, a common representation of a sparse matrix, to tensor data. By combining few elementary operations on sparse maps (inversion, chaining, intersection, etc.), complex algorithms can be developed, illustrated here by a linear-scaling transformation of three-center Coulomb integrals based on our compact code library that implements sparse maps and operations on them. The sparsity of the three-center integrals arises from spatial locality of the basis functions and domain density fitting approximation. A novel feature of our approach is the use of differential overlap integrals computed in linear-scaling fashion for screening products of basis functions. Finally, a robust linear scaling domain based local pair natural orbital second-order Möller-Plesset (DLPNO-MP2) method is described based on the sparse map infrastructure that only depends on a minimal number of cutoff parameters that can be systematically tightened to approach 100% of the canonical MP2 correlation energy. With default truncation thresholds, DLPNO-MP2 recovers more than 99.9% of the canonical resolution of the identity MP2 (RI-MP2) energy while still showing a very early crossover with respect to the computational effort. Based on extensive benchmark calculations, relative energies are reproduced with an error of typically <0.2 kcal/mol. The efficiency of the local MP2 (LMP2) method can be drastically improved by carrying out the LMP2 iterations in a basis of pair natural orbitals. While the present work focuses on local electron correlation, it is of much broader applicability to computation with sparse tensors in quantum chemistry and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Systems, methods and apparatus for generation and verification of policies in autonomic computing systems

Described herein is a method that produces fully (mathematically) tractable development of policies for autonomic systems from requirements through to code generation. This method is illustrated through an example showing how user formulated policies can be translated into a formal mode which can then be converted to code. The requirements-based programming method described provides faster, higher quality development and maintenance of autonomic systems based on user formulation of policies.Further, the systems, methods and apparatus described herein provide a way of analyzing policies for autonomic systems and facilities the generation of provably correct implementations automatically, which in turn provides reduced development time, reduced testing requirements, guarantees of correctness of the implementation with respect to the policies specified at the outset, and provides a higher degree of confidence that the policies are both complete and reasonable. The ability to specify the policy for the management of a system and then automatically generate an equivalent implementation greatly improves the quality of software, the survivability of future missions, in particular when the system will operate untended in very remote environments, and greatly reduces development lead times and costs.

Hinchey, Michael G.↗
NASA NTRS1987

Heating of the solar corona by the resonant absorption of Alfven waves

An improved method for calculating the resonance absorption heating rate is discussed and the results are compared with observations in the solar corona. To accomplish this, the wave equation for a dissipative, compressible plasma is derived from the linearized magnetohydrodynamic equations for a plasma with transverse Alfven speed gradients. For parameters representative of the solar corona, it is found that a two-scale description of the wave motion is appropriate. The large-scale motion, which can be approximated as nearly ideal, has a scale which is on the order of the width of the loop. The small-scale wave, however, has a transverse scale much smaller than the width of the loop, with a width of about 0.3-250 km, and is highly dissipative. These two wave motions are coupled in a narrow resonance region in the loop where the global wave frequency equals the local Alfven wave frequency. Formally, this coupling comes about from using the method of matched asymptotic expansions to match the inner and outer (small and large scale) solutions. The resultant heating rate can be calculated from either of these solutions. A formula derived using the outer (ideal) solution is presented, and shown to be consistent with observations of heating and line broadening in the solar corona.

Davila, Joseph M.↗
DOE OSTI2016

A feasibility study of dynamic adaptive radiotherapy for nonsmall cell lung cancer

Purpose: The final state of the tumor at the end of a radiotherapy course is dependent on the doses given in each fraction during the treatment course. This study investigates the feasibility of using dynamic adaptive radiotherapy (DART) in treating lung cancers assuming CBCT is available to observe midtreatment tumor states. DART adapts treatment plans using a dynamic programming technique to consider the expected changes of the tumor in the optimization process. Methods: DART is constructed using a stochastic control formalism framework. It minimizes the total expected number of tumor cells at the end of a treatment course, which is equivalent to maximizing tumor control probability, subject to the uncertainty inherent in the tumor response. This formulation allows for nonstationary dose distributions as well as nonstationary fractional doses as needed to achieve a series of optimal plans that are conformal to the tumor over time, i.e., spatiotemporally optimal plans. Sixteen phantom cases with various sizes and locations of tumors and organs-at-risk (OAR) were generated using in-house software. Each case was planned with DART and conventional IMRT prescribing 60 Gy in 30 fractions. The observations of the change in the tumor volume over a treatment course were simulated using a two-level cell population model. Monte Carlo simulations of the treatment course for each case were run to account for uncertainty in the tumor response. The same OAR dose constraints were applied for both methods. The frequency of replanning was varied between 1, 2, 5 (weekly), and 29 times (daily). The final average tumor dose and OAR doses have been compared to quantify the potential dosimetric benefits of DART. Results: The average tumor max, min, mean, and D95 doses using DART relative to these using conventional IMRT were 124.0%–125.2%, 102.1%–114.7%, 113.7%–123.4%, and 102.0%–115.9% (range dependent on the frequency of replanning). The average relative maximum doses for the cord and esophagus, mean doses for the heart and lungs, and D05 for the unspecified tissue resulting 84%–102.4%, 99.8%–106.9%, 66.9%–85.6%, 58.2%–78.8%, and 85.2%–94.0%, respectively. Conclusions: It is feasible to apply DART to the treatment of NSCLC using CBCT to observe the midtreatment tumor state. Potential increases in the tumor dose and reductions in the OAR dose, particularly for parallel OARs with mean or dose–volume constraints, could be achieved using DART compared to nonadaptive IMRT.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2014

Gradient Calculation Methods on Arbitrary Polyhedral Unstructured Meshes for Cell-Centered CFD Solvers

A survey of gradient reconstruction methods for cell-centered data on unstructured meshes is conducted within the scope of accuracy assessment. Formal order of accuracy, as well as error magnitudes for each of the studied methods, are evaluated on a complex mesh of various cell types through consecutive local scaling of an analytical test function. The tests highlighted several gradient operator choices that can consistently achieve 1st order accuracy regardless of cell type and shape. The tests further offered error comparisons for given cell types, leading to the observation that the "ideal" gradient operator choice is not universal. Practical implications of the results are explored via CFD solutions of a 2D inviscid standing vortex, portraying the discretization error properties. A relatively naive, yet largely unexplored, approach of local curvilinear stencil transformation exhibited surprisingly favorable properties

Meshes↗
DOE OSTI2017

First-principles Green-Kubo method for thermal conductivity calculations

We present a first-principles approach to calculate the phonon thermal conductivity based on the Green-Kubo formalism. In this approach, the density functional theory energy is distributed to each atom, and a two-step method in the molecular dynamics is introduced to avoid the atomic position R wrapping problem in a periodic system when the heat current is calculated. We show that this first-principles Green-Kubo approach is particularly suitable for disordered systems like amorphous and liquid, where the thermal conductivities are small due to strong phonon scattering but difficult to be calculated using anharmonic interaction energy. We have applied our method to liquid Ar, liquid Si, and amorphous Si. In this work, the calculated thermal conductivities agree well with previous theoretical and experimental results. We have also compared our method to previous works combining first-principles simulations with the Green-Kubo formalism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Quantum mechanics of a constrained particle on an ellipsoid: Bein formalism and Geometric momentum

In this work we apply the Dirac method in order to obtain the classical relations for a particle on an ellipsoid. We also determine the quantum mechanical form of these relations by using Dirac quantization. Then by considering the canonical commutation relations between the position and momentum operators in terms of curved coordinates, we try to propose the suitable representations for momentum operator that satisfy the obtained commutators between position and momentum in Euclidean space. We see that our representations for momentum operators are the same as geometric one.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

Use of a minimum rate of change formalism to quantify variability of extragalactic X-ray sources

A method to obtain rigorous, quantitative constraints on the rate of variability is suggested. The method is motivated by the case when the existence of time variability is unequivocal, but the statistical uncertainties are too large to apply more direct methods such as power spectrum analysis. Conceptually the data are fitted to a function of time and a set of free parameters. The method of Lagrange multipliers is used to solve for that parameter set which minimizes the value of the derivative at time t(m). The result can be related physically to the minimum efficiency with which rest mass must be converted to radiation energy. Slightly different formulations of the general principal are used to estimate maximum source size scales associated with the variability.

Schwartz, Daniel A.↗
DOE OSTI2025

A Methodology for the Analysis of Water Oxidation Electrocatalysts in the Absence of Limiting Current that Avoids the Pitfalls of Existing Methods

Water oxidation is an important reaction studied as a way to generate electrons from water, to promote water splitting and the formation of green hydrogen. When using electrodes to drive homogeneous water oxidation catalysis, cyclic voltammograms are analyzed to provide catalytic rate constants. There are two main methods, foot-of-the-wave analysis (FOWA) and limiting current analysis. FOWA relies on approximations inherent to analyzing water oxidation catalysis, such as determining the formal potential of the catalytic intermediate, E 0 cat . Limiting current methods are the optimal way to analyze catalyst performance but rely on observable limiting current, which is virtually never seen in water oxidation. To avoid those issues, a method is proposed for analyzing nonideal cyclic voltammetry waveshapes in water oxidation: by analyzing rate data across a large range of potentials, an optimal potential, E 0 cat , can be obtained, where catalytic current, i cat , is nearly independent of scan rate and has a linear dependency on buffer concentration. Here, the method is applied to four homogeneous water oxidation catalysts with prior extensive electrochemical elucidation, all of which lack an ideal, purely kinetic waveshape in cyclic voltammetry. Application of the method avoids the biases of the other methods cited for the kinetic analyses of water oxidation catalysts.

14 SOLAR ENERGY↗
NASA NTRS1992

Design and Analysis Techniques for Concurrent Blackboard Systems

Blackboard systems are a natural progression of knowledge-based systems into a more powerful problem solving technique. They provide a way for several highly specialized knowledge sources to cooperate to solve large, complex problems. Blackboard systems incorporate the concepts developed by rule-based and expert systems programmers and include the ability to add conventionally coded knowledge sources. The small and specialized knowledge sources are easier to develop and test, and can be hosted on hardware specifically suited to the task that they are solving. The Formal Model for Blackboard Systems was developed to provide a consistent method for describing a blackboard system. A set of blackboard system design tools has been developed and validated for implementing systems that are expressed using the Formal Model. The tools are used to test and refine a proposed blackboard system design before the design is implemented. My research has shown that the level of independence and specialization of the knowledge sources directly affects the performance of blackboard systems. Using the design, simulation, and analysis tools, I developed a concurrent object-oriented blackboard system that is faster, more efficient, and more powerful than existing systems. The use of the design and analysis tools provided the highly specialized and independent knowledge sources required for my concurrent blackboard system to achieve its design goals.

Mcmanus, John William↗
DOE OSTI2016

Wave function continuity and the diagonal Born-Oppenheimer correction at conical intersections

We demonstrate that though exact in principle, the expansion of the total molecular wave function as a sum over adiabatic Born-Oppenheimer (BO) vibronic states makes inclusion of the second-derivative nonadiabatic energy term near conical intersections practically problematic. In order to construct a well-behaved molecular wave function that has density at a conical intersection, the individual BO vibronic states in the summation must be discontinuous. When the second-derivative nonadiabatic terms are added to the Hamiltonian, singularities in the diagonal BO corrections (DBOCs) of the individual BO states arise from these discontinuities. In contrast to the well-known singularities in the first-derivative couplings at conical intersections, these singularities are non-integrable, resulting in undefined DBOC matrix elements. Though these singularities suggest that the exact molecular wave function may not have density at the conical intersection point, there is no physical basis for this constraint. Instead, the singularities are artifacts of the chosen basis of discontinuous functions. We also demonstrate that continuity of the total molecular wave function does not require continuity of the individual adiabatic nuclear wave functions. We classify nonadiabatic molecular dynamics methods according to the constraints placed on wave function continuity and analyze their formal properties. Based on our analysis, it is recommended that the DBOC be neglected when employing mixed quantum-classical methods and certain approximate quantum dynamical methods in the adiabatic representation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1982

Discrete Fourier transforms of nonuniformly spaced data

Time series or spatial series of measurements taken with nonuniform spacings have failed to yield fully to analysis using the Discrete Fourier Transform (DFT). This is due to the fact that the formal DFT is the convolution of the transform of the signal with the transform of the nonuniform spacings. Two original methods are presented for deconvolving such transforms for signals containing significant noise. The first method solves a set of linear equations relating the observed data to values defined at uniform grid points, and then obtains the desired transform as the DFT of the uniform interpolates. The second method solves a set of linear equations relating the real and imaginary components of the formal DFT directly to those of the desired transform. The results of numerical experiments with noisy data are presented in order to demonstrate the capabilities and limitations of the methods.

Swan, P. R.↗
NASA NTRS2008

Applications of the Hybrid Theory to the Scattering of Electrons from HE+ and Li++ and Resonances in these Systems

Applications of the hybrid theory to the scattering of electrons from Ile+ and Li++ and resonances in these systems, A. K. Bhatia, NASA/Goddard Space Flight Center- The Hybrid theory of electron-hydrogen elastic scattering [I] is applied to the S-wave scattering of electrons from He+ and Li++. In this method, both short-range and long-range correlations are included in the Schrodinger equation at the same time. Phase shifts obtained in this calculation have rigorous lower bounds to the exact phase shifts and they are compared with those obtained using the Feshbach projection operator formalism [2], the close-coupling approach [3], and Harris-Nesbet method [4]. The agreement among all the calculations is very good. These systems have doubly-excited or Feshbach resonances embedded in the continuum. The resonance parameters for the lowest ' S resonances in He and Li+ are calculated and they are compared with the results obtained using the Feshbach projection operator formalism [5,6]. It is concluded that accurate resonance parameters can be obtained by the present method, which has the advantage of including corrections due to neighboring resonances and the continuum in which these resonances are embedded.

Bhatia, Anand K.↗
DOE OSTI2020

Bismuth Detection in Alkaline Electrolyte via Anodic Stripping Voltammetry for Battery Separator Evaluation

Abstract Anodic stripping voltammetry (ASV) has been widely used for the detection of several heavy metal ions in neutral and acidic solution, in many cases employing electrodes and/or solutions incorporating Bi. In this work we demonstrate that Bi(OH) 4 − ion concentration can be measured in highly alkaline 8.5 M KOH solution using ASV. The addition of Pb in similar concentrations to the Bi(OH) 4 − being measured is shown to improve both the sensitivity and precision of the method. When the Pb additive is employed, a formal limit of detection of 8.5 ppb is achieved, compared to 17.3 ppb when the Pb additive is not used. Due to the use of Bi additives in alkaline battery chemistries, it follows that separators which limit Bi(OH) 4 − diffusion into the bulk electrolyte and away from the electrodes are of interest. For this purpose, we utilize ASV to determine Bi(OH) 4 − diffusion rates through Celgard 3501, cellophane 350P00, and Nafion 211. Bi(OH) 4 − crossover rates, as determined by ASV, are shown to be repeatable and consistent with expectations from the known separator structure.

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

The photon content of the neutron

In this work, we complete our CT18qed study with the neutron’s photon parton distribution function (PDF), which is essential for the nucleus scattering phenomenology. Two methods, CT18lux and CT18qed, based on the LUXqed formalism and the DGLAP evolution, respectively, to determine the neutron’s photon PDF have been presented. Various low-Q 2 non-perturbative variations have been carefully examined, which are treated as additional uncertainties on top of those induced by quark and gluon PDFs. The impacts of the momentum sum rule as well as isospin symmetry violation have been explored and turned out to be negligible. A detailed comparison with other neutron’s photon PDF sets has been performed, which shows a great improvement in the precision and a reasonable uncertainty estimation. Finally, two phenomenological implications are demonstrated with photon-initiated processes: neutrino-nucleus W-boson production, which is important for the near-future TeV–PeV neutrino observations, and the axion-like particle production at a high-energy muon beam-dump experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Higher-order tails and RG flows due to scattering of gravitational radiation from binary inspirals

Abstract We establish and develop a novel methodology to treat higher-order non-linear effects of gravitational radiation that is scattered from binary inspirals, which employs modern scattering-amplitudes methods on the effective picture of the binary as a composite particle. We spell out our procedure to study such effects: assembling tree amplitudes via generalized-unitarity methods and employing the closed-time-path formalism to derive the causal effective actions, which encompass the full conservative and dissipative dynamics. We push through to a new state of the art for these higher-order effects, up to the third subleading tail effect, at order$$ {G}_N^5 $$ G N 5 and the 5-loop level, which corresponds to the 8.5PN order. We formulate the consequent dissipated energy for these higher-order corrections, and carry out a renormalization analysis, where we uncover new subleading RG flow of the quadrupole coupling. For all higher-order tail effects we find perfect agreement with partial observable results in PN and self-force theories, where available.

Physics↗