Search NASA⌕ Search

SEARCH · Search NASA

Results for “standard problem”

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 55 records · Page 3

Flight Software Math Library

The flight software (FSW) math library is a collection of reusable math components that provides typical math utilities required by spacecraft flight software. These utilities are intended to increase flight software quality reusability and maintainability by providing a set of consistent, well-documented, and tested math utilities. This library only has dependencies on ANSI C, so it is easily ported. Prior to this library, each mission typically created its own math utilities using ideas/code from previous missions. Part of the reason for this is that math libraries can be written with different strategies in areas like error handling, parameters orders, naming conventions, etc. Changing the utilities for each mission introduces risks and costs. The obvious risks and costs are that the utilities must be coded and revalidated. The hidden risks and costs arise in miscommunication between engineers. These utilities must be understood by both the flight software engineers and other subsystem engineers (primarily guidance navigation and control). The FSW math library is part of a larger goal to produce a library of reusable Guidance Navigation and Control (GN&C) FSW components. A GN&C FSW library cannot be created unless a standardized math basis is created. This library solves the standardization problem by defining a common feature set and establishing policies for the library s design. This allows the libraries to be maintained with the same strategy used in its initial development, which supports a library of reusable GN&C FSW components. The FSW math library is written for an embedded software environment in C. This places restrictions on the language features that can be used by the library. Another advantage of the FSW math library is that it can be used in the FSW as well as other environments like the GN&C analyst s simulators. This helps communication between the teams because they can use the same utilities with the same feature set and syntax.

McComas, David↗

Radial transmission line analysis of multi-layer structures

The analysis of multi-layer beam tubes is a standard problem and involves axially propagating waves. This treatment is ill suited to a short multi-layer structure such as the present example of a ferrite covered ceramic break in the beam tube at the ERL photo-cathode electron gun. This paper demonstrates that such structure can better be treated by radial wave propagation. The theoretical method is presented and numerical results are compared with measured network analyser data and Microwave Studio generated simulations. The results confirm the concept of radial transmission lines as a valid analytical method. An Energy Recovery Linac (ERL) is being constructed at this laboratory for the purpose of research towards an envisioned Electron Ion Collider. One of the pertinent topics is damping of Higher Order Modes (HOM). In this ERL, the damping is provided by ferrite absorbers in the beam tube. A modified version thereof, a ceramic break surrounded by ferrite, is planed for the superconducting electron gun. The damper here is located at room temperature just outside of the gun. If used in a cavity chain, the ceramic break is in the vacuum tube at helium temperature whereas the ferrite is moved into the cryostat insulating vacuum allowing higher temperatures. The general properties of the ferrite HOM dampers have been published but are more detailed in this paper.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Error analysis of householder transformations as applied to the standard and generalized eigenvalue problems

Backward error analyses of the application of Householder transformations to both the standard and the generalized eigenvalue problems are presented. The analysis for the standard eigenvalue problem determines the error from the application of an exact similarity transformation, and the analysis for the generalized eigenvalue problem determines the error from the application of an exact equivalence transformation. Bounds for the norms of the resulting perturbation matrices are presented and compared with existing bounds when known.

Ward, R. C.↗

Investigation of some problems in developing standards for precracked Charpy slow bend tests

The reported investigation was undertaken in connection with an attempt to develop procedures which would be useful in standardizing a test method for the precracked Charpy slow bend specimen. A number of alloys was studied for which valid plane-strain fracture toughness values have been established. The investigation shows that useful relations between precracked Charpy slow bend results and crack size factors can be obtained under some circumstances. However, it is not yet known what factors control these circumstances.

Succop, G.↗

Problems for the standard black hole/accretion disk models in Cygnus X-1?

Archival EXOSAT and HEAO1-A2 data from Cyg X-1 show the 'high energy excess' above 10 keV seen in X-ray observations of AGN. Using a likelihood ratio test, we are for the first time able to distinguish conclusively in favor of Compton reflection rather than partial covering as the origin of the high energy excess. This supports the idea of an X-ray illuminated accretion disk in Cyg X-1, but the line equivalent width is smaller by a factor of 2-3 than that expected from such a disk. While the larger optical depth required for reflection as opposed to line emission admit the possibility of seeing line without reflection, the converse is not possible. To see a reflection spectrum, including the strong iron absorption edge, implies that strong iron emission must be observed as the line and edge are causally linked.

Done, C.↗

Standardized Definitions for Code Verification Test Problems

This document contains standardized definitions for several commonly used code verification test problems. These definitions are intended to contain sufficient information to set up the test problem in a computational physics code. These definitions are intended to be used in conjunction with exact solutions to these problems generated using Exact- Pack, www.github.com/lanl/exactpack.

97 MATHEMATICS AND COMPUTING↗

COMPARISONS OF THE FINITE-ELEMENT-WITH-DISCONTIGUOUS-SUPPORT METHOD TO CONTINUOUS-ENERGY MONTE CARLO FOR PIN-CELL PROBLEMS

The standard multigroup (MG) method for energy discretization of the transport equation can be sensitive to approximations in the weighting spectrum chosen for cross-section averaging. As a result, MG often inaccurately treats important phenomena such as self-shielding variations across a material. From a finite-element viewpoint, MG uses a single fixed basis function (the pre-selected spectrum) within each group, with no mechanism to adapt to local solution behavior. In this work, we introduce the Finite-Element-with-Discontiguous-Support (FEDS) method, whose only approximation with respect to energy is that the angular flux is a linear combination of unknowns multiplied by basis functions. A basis function is non-zero only in the discontiguous set of energy intervals associated with its energy element. Discontiguous energy elements are generalizations of bands and are determined by minimizing a norm of the difference between snapshot spectra and their averages over the energy elements. We begin by presenting the theory of the FEDS method. We then compare to continuous-energy Monte Carlo for one-dimensional slab and two-dimensional pin-cell problem. We find FEDS to be accurate and efficient at producing quantities of interest such as reaction rates and eigenvalues. Results show that FEDS converges at a rate that is approximately first-order in the number of energy elements and that FEDS is less sensitive to weighting spectrum than standard MG.

Cross sections↗

Deterministically estimated fission source distributions for Monte Carlo k-eigenvalue problems

The standard Monte Carlo (MC) k-eigenvalue algorithm involves iteratively converging the fission source distribution using a series of potentially time-consuming inactive cycles before quantities of interest can be tallied. One strategy for reducing the computational time requirements of these inactive cycles is the Sourcerer method, in which a deterministic eigenvalue calculation is performed to obtain an improved initial guess for the fission source distribution. This method has been implemented in the Exnihilo software suite within SCALE using the SPNSPN or SNSN solvers in Denovo and the Shift MC code. The efficacy of this method is assessed with different Denovo solution parameters for a series of typical k-eigenvalue problems including small criticality benchmarks, full-core reactors, and a fuel cask. Here it is found that, in most cases, when a large number of histories per cycle are required to obtain a detailed flux distribution, the Sourcerer method can be used to reduce the computational time requirements of the inactive cycles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

H→γγ as a Triangle Anomaly: Possible Implications for the Hierarchy Problem

The Standard Model calculation of H→γγ has the curious feature of being finite but regulator-dependent. While dimensional regularization yields a result which respects the electromagnetic Ward identities, additional terms which violate gauge invariance arise if the calculation is done setting d = 4. This discrepancy between the d=4 – ϵ and d = 4 results is recognized as a true ambiguity which must be resolved using physics input; as dimensional regularization respects gauge invariance, the d = 4 – ϵ calculation is accepted as the correct SM result. However, here we point out another possibility; working in analogy with the gauge chiral anomaly, we note that it is possible that the individual diagrams do violate the electromagnetic Ward identities, but that the gauge-invariance-violating terms cancel when all contributions to H→γγ, both from the SM and from new physics, are included. We thus examine the consequences of the hypothesis that the d = 4 calculation is valid, but that such a cancellation occurs. We work in general renormalizable gauge, thus avoiding issues with momentum routing ambiguities. We point out that the gauge-invariance-violating terms in d = 4 arise not just for the diagram containing a SM $W^{\pm}$ boson, but also for general fermion and scalar loops, and relate these terms to a lack of shift invariance in Higgs tadpole diagrams. We then derive the analogue of "anomaly cancellation conditions", and find consequences for solutions to the hierarchy problem. In particular, we find that supersymmetry obeys these conditions, even if it is softly broken at an arbitrarily high scale.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Concepts for on-board satellite image registration. Volume 2: IAS prototype performance evaluation standard definition

Problems encountered in testing onboard signal processing hardware designed to achieve radiometric and geometric correction of satellite imaging data are considered. These include obtaining representative image and ancillary data for simulation and the transfer and storage of a large quantity of image data at very high speed. The high resolution, high speed preprocessing of LANDSAT-D imagery is considered.

Daluge, D. R.↗

Dipolar dark matter with massive bigravity

Massive gravity theories have been developed as viable IR modifications of gravity motivated by dark energy and the problem of the cosmological constant. On the other hand, modified gravity and modified dark matter theories were developed with the aim of solving the problems of standard cold dark matter at galactic scales. Here we propose to adapt the framework of ghost-free massive bigravity theories to reformulate the problem of dark matter at galactic scales. We investigate a promising alternative to dark matter called dipolar dark matter (DDM) in which two different species of dark matter are separately coupled to the two metrics of bigravity and are linked together by an internal vector field. We show that this model successfully reproduces the phenomenology of dark matter at galactic scales (i.e. MOND) as a result of a mechanism of gravitational polarisation. The model is safe in the gravitational sector, but because of the particular couplings of the matter fields and vector field to the metrics, a ghost in the decoupling limit is present in the dark matter sector. However, it might be possible to push the mass of the ghost beyond the strong coupling scale by an appropriate choice of the parameters of the model. Crucial questions to address in future work are the exact mass of the ghost, and the cosmological implications of the model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗

Solving the Trivial Crossing Problem While Preserving the Nodal Symmetry of the Wave Function

In an adiabatic mixed quantum-classical simulation, the avoided crossing of weakly coupled eigenstates can lead to unphysical discontinuities in wave function dynamics, otherwise known as the trivial crossing problem. A standard solution to the trivial crossing problem eliminates spatial discontinuities in wave function dynamics by imposing changes to the eigenstate of the wave function. In this paper, we show that this solution has the side effect of introducing transient discontinuities in the nodal symmetry of the wave function. We present an alternative solution to the trivial crossing problem that preserves both the spatial and nodal structure of the adiabatic wave function. By considering a model of exciton dynamics on conjugated polymer systems, we show that failure to preserve wave function symmetry yields exciton dynamics that depends unphysically on polymer system size. We demonstrate that our symmetry-preserving solution to the trivial crossing problem yields more realistic dynamics and can thus improve the accuracy of simulations of larger systems that are prone to the trivial crossing problem.

74 ATOMIC AND MOLECULAR PHYSICS↗

Robust Resilient Signal Reconstruction under Adversarial Attacks

We consider the problem of signal reconstruction for a system under sparse signal corruption by a malicious agent. The reconstruction problem follows the standard error coding problem that has been studied extensively in the literature. We include a new challenge of robust estimation of the attack support. The problem is then cast as a constrained optimization problem merging promising techniques in the area of deep learning and estimation theory. A pruning algorithm is developed to reduce the "false positive" uncertainty of data-driven attack localization results, thereby improving the probability of correct signal reconstruction. Sufficient conditions for the correct reconstruction and the associated reconstruction error bounds are obtained for both exact and inexact attack support estimation. Moreover, a simulation of a water distribution system is presented to validate the proposed techniques.

Robust, Signal reconstruction, Resilient estimator↗

Verification of technical elements of the advanced spacecraft based upon the CCSDS recommendation

We are going to meet the era when advanced spacecraft such as space stations are developed and operated. The current system of the satellite operations control will need to undergo many changes. We consider that the future system will require the following functions: the function for interchanging data between international agencies, processing the various kinds of space data, and distributing data as many unspecified users require. However, we have to solve the following problems in order to satisfy these requirements: the problem of standardization of space data communication protocol, establishment of multimedia data management method, and standardization of the user interface. This paper describes three techniques to solve the above mentioned problems. That is, standardization of the data communication protocol between space and ground by AOS (Advanced Orbiting System) protocol of CCSDS (Consultative Committee for Space Data Systems) Recommendation, management of multimedia data by catalog reference, standardization of user interface by SFDU(Standard Formatted Data Unit) of CCSDS Recommendation.

Hara, Hideo↗

Quantum phase estimation for a class of generalized eigenvalue problems

Quantum phase estimation provides a path to quantum computation of solutions to Hermitian eigenvalue problems Hv = λv , such as those occurring in quantum chemistry. It is natural to ask whether the same technique can be applied to generalized eigenvalue problems Av = λBv , which arise in many areas of science and engineering. Here, we answer this question affirmatively. A restricted class of generalized eigenvalue problems could be solved as efficiently as standard eigenvalue problems. A paradigmatic example is provided by Sturm-Liouville problems. Another example comes from linear ideal magnetohydrodynamics, where phase estimation could be used to determine the stability of magnetically confined plasmas in fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SEMI-ACTIVE CONTROL FOR TWO-BODY OCEAN WAVE ENERGY CONVERTER BY USING HYBRID MODEL PREDICTIVE CONTROL

Model predictive control (MPC) has been considered as one important feed-forward optimal control strategy for ocean wave energy converter (WEC) targeted on power maximization. The capability of MPC to handle system constraints (ex. stroke, velocity, actuator limitations), and the availability to provide optimal solution for linear system provide potential for the implementation of such algorithm in the WEC control. However, currently, only active MPC control has been introduced for single and two-body WECs. Such control strategy may introduce negative power during the optimization process, since the power take-off (PTO) damping has no constraint. In this paper, we proposed a hybrid MPC strategy in limiting both the PTO damping force and PTO damping to avoid negative power generation during cost function minimization (negative power minimization) for the two-body WEC. The problem is formulated into a quadratic programming (QP) problem targeted at power maximization. However, the standard QP problem formulation cannot be directly applied to the semi-active control problem due to the PTO damping constraints. Therefore, the problem is reformulated as a Mixed-integer Quadratic Programming (MIQP) problem, which contains logical switch to select constraint matrices based on the sign of the relative velocity between the buoy and submerged body. The optimal solution is compared with those of the active MPC control strategy and the passive model with the same irregular wave input.

16 TIDAL AND WAVE POWER↗