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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Spatial correlations in driven-dissipative photonic lattices

In this paper, we study the nonequilibrium steady-state of interacting photons in cavity arrays as described by the driven-dissipative Bose–Hubbard and spin-1/2 XY model. For this purpose, we develop a self-consistent expansion in the inverse coordination number of the array ($\sim 1/z$) to solve the Lindblad master equation of these systems beyond the mean-field approximation. Our formalism is compared and benchmarked with exact numerical methods for small systems based on an exact diagonalization of the Liouvillian and a recently developed corner-space renormalization technique. We then apply this method to obtain insights beyond mean-field in two particular settings: (i) we show that the gas–liquid transition in the driven-dissipative Bose–Hubbard model is characterized by large density fluctuations and bunched photon statistics. (ii) We study the antibunching–bunching transition of the nearest-neighbor correlator in the driven-dissipative spin-1/2 XY model and provide a simple explanation of this phenomenon.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

Efficient optimal reconstruction of linear fields and band-powers from cosmological data

We present an efficient implementation of Wiener filtering of real-space linear field and optimal quadratic estimator of its power spectrum Band-powers. We first recast the field reconstruction into an optimization problem, which we solve using quasi-Newton optimization. We then recast the power spectrum estimation into the field marginalization problem, from which we obtain an expression that depends on the field reconstruction solution and a determinant term. We develop a novel simulation based method for the latter. We extend the simulations formalism to provide the covariance matrix for the power spectrum. We develop a flexible framework that can be used on a variety of cosmological fields and present results for a variety of test cases, using simulated examples of projected density fields, projected shear maps from galaxy lensing, and observed Cosmic Microwave Background (CMB) temperature anisotropies, with a wide range of map incompleteness and variable noise. For smaller cases where direct numerical inversion is possible, we show that our solution matches that created by direct Wiener Filtering at a fraction of the overall computation cost. Even more significant reduction of computational is achieved by this implementation of optimal quadratic estimator due to the fast evaluation of the Hessian matrix. This technique allows for accurate map and power spectrum reconstruction with complex masks and nontrivial noise properties.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces – part II

The modeling of the near-field interaction in the scattering-type scanning near-field optical microscope (s-SNOM) is rapidly advancing, although an accurate yet versatile modeling framework that can be easily adapted to various complex situations is still lacking. In this work, we propose a time-efficient numerical scheme in the quasi-electrostatic limit to capture the tip-sample interaction in the near field. This method considers an extended tip geometry, which is a significant advantage compared to the previously reported method based on the point-dipole approximation. Using this formalism, we investigate, among others, nontrivial questions such as uniaxial and biaxial anisotropy in the near-field interaction, the relationship between various experimental parameters (e.g. tip radius, tapping amplitude, etc.), and the tip-dependent spatial resolution. The demonstrated method further sheds light on the understanding of the contrast mechanism in s-SNOM imaging and spectroscopy, while also representing a valuable platform for future quantitative analysis of the experimental observations.

Chen, Xinzhong (ORCID:0000000261033848)↗

NASA NTRS1972

Solution of non-isoenergetic supersonic flows by method of characteristics, volume 3

The calculation of supersonic flow fields by the method of characteristics. The theoretical approach to the solution of these flow fields and a computer program to implement the numerical solution of the flow equations are discussed. This versatile program has a flexible set of boundary conditions enabling the calculation of nozzles, plumes and many other complex flow fields. A complete derivation of the equations of motion for reacting gas systems is presented. An important consequence of this derivation is that, for the reaction assumptions which were made, the thermochemistry was shown to be uncoupled from the flow solution and as such could be solved separately. The methods of characteristics equations are shown to be formally the same for ideal, frozen, and equilibrium reacting gas mixtures.

Prozan, R. J.↗

NASA NTRS1982

Numerical computation of exponential matrices using the Cayley-Hamilton theorem

A method for computing exponential matrices, which often arise naturally in the solution of systems of linear differential equations, is developed. An exponential matrix is generated as a linear combination of a finite number (equal to the matrix order) of matrices, the coefficients of which are scalar infinite sums. The method can be generalized to apply to any formal power series of matrices. Attention is focused upon the exponential function, and the matrix exponent is assumed tri-diagonal in form. In such cases, the terms in the coefficient infinite sums can be extracted, as recursion relations, from the characteristic polynomial of the matrix exponent. Two numerical examples are presented in some detail: (1) the three dimensional infinitesimal rotation rate matrix, which is skew symmetric, and (2) an N-dimensional tri-diagonal and symmetric finite difference matrix which arises in the numerical solution of the heat conduction partial differential equation. In the second example, the known eigenvalues and eigenvectors of the finite difference matrix permit an analytical solution for the exponential matrix, through the theory of diagonalization and similarity transformations, which is used for independent verification. The convergence properties of the scalar infinite summations are investigated for finite difference matrices of various orders up to ten, and it is found that the number of terms required for convergence increases slowly with the order of the matrix.

Walden, H.↗

NASA NTRS1984

Configuration study for a 30 GHz monolithic receive array, volume 2

The formalism of the sidelobe suppression algorithm and the method used to calculate the system noise figure for a 30 GHz monolithic receive array are presented. Results of array element weight determination and performance studies of a Gregorian aperture image system are also given.

Nester, W. H.↗

NASA NTRS1994

Program Helps Design Tests Of Developmental Software

Computer program called "A Formal Test Representation Language and Tool for Functional Test Designs" (TRL) provides automatic software tool and formal language used to implement category-partition method and produce specification of test cases in testing phase of development of software. Category-partition method useful in defining input, outputs, and purpose of test-design phase of development and combines benefits of choosing normal cases having error-exposing properties. Traceability maintained quite easily by creating test design for each objective in test plan. Effort to transform test cases into procedures simplified by use of automatic software tool to create cases based on test design. Method enables rapid elimination of undesired test cases from consideration and facilitates review of test designs by peer groups. Written in C language.

Hops, Jonathan↗

NASA NTRS1995

Uncertainties for two-dimensional models of solar rotation from helioseismic eigenfrequency splitting

Observed solar p-mode frequency splittings can be used to estimate angular velocity as a function of position in the solar interior. Formal uncertainties of such estimates depend on the method of estimation (e.g., least-squares), the distribution of errors in the observations, and the parameterization imposed on the angular velocity. We obtain lower bounds on the uncertainties that do not depend on the method of estimation; the bounds depend on an assumed parameterization, but the fact that they are lower bounds for the 'true' uncertainty does not. Ninety-five percent confidence intervals for estimates of the angular velocity from 1986 Big Bear Solar Observatory (BBSO) data, based on a 3659 element tensor-product cubic-spline parameterization, are everywhere wider than 120 nHz, and exceed 60,000 nHz near the core. When compared with estimates of the solar rotation, these bounds reveal that useful inferences based on pointwise estimates of the angular velocity using 1986 BBSO splitting data are not feasible over most of the Sun's volume. The discouraging size of the uncertainties is due principally to the fact that helioseismic measurements are insensitive to changes in the angular velocity at individual points, so estimates of point values based on splittings are extremely uncertain. Functionals that measure distributed 'smooth' properties are, in general, better constrained than estimates of the rotation at a point. For example, the uncertainties in estimated differences of average rotation between adjacent blocks of about 0.001 solar volumes across the base of the convective zone are much smaller, and one of several estimated differences we compute appears significant at the 95% level.

Genovese, Christopher R.↗

NASA NTRS2001

Science@NASA: Direct to People Via the Internet

NASA's founding charter includes the requirement for reporting all scientific results to the public. This requirement is based on the principal that the exploration of space results in real benefits to humanity and that those benefits are to be shared as widely as practical. When NASA was founded, the traditional education and outreach methods were through the news media and the formal and informal (museums, planetariums exhibits, etc.) educational communities. With the nearly ubiquitous availability of the Internet, a third choice presents itself: communicating directly with individuals in their homes. This powerful approach offers benefits and pitfalls that must be addressed to be effective. This paper covers an integrated approach to providing high quality NASA research information to multiple audiences via a family of websites. The paper discuss the content generation, review, and production process and provide metrics on evaluating the results.

Koczor, R. J.↗

NASA NTRS2015

Kodiak: An Implementation Framework for Branch and Bound Algorithms

Recursive branch and bound algorithms are often used to refine and isolate solutions to several classes of global optimization problems. A rigorous computation framework for the solution of systems of equations and inequalities involving nonlinear real arithmetic over hyper-rectangular variable and parameter domains is presented. It is derived from a generic branch and bound algorithm that has been formally verified, and utilizes self-validating enclosure methods, namely interval arithmetic and, for polynomials and rational functions, Bernstein expansion. Since bounds computed by these enclosure methods are sound, this approach may be used reliably in software verification tools. Advantage is taken of the partial derivatives of the constraint functions involved in the system, firstly to reduce the branching factor by the use of bisection heuristics and secondly to permit the computation of bifurcation sets for systems of ordinary differential equations. The associated software development, Kodiak, is presented, along with examples of three different branch and bound problem types it implements.

Smith, Andrew P.↗

NASA NTRS2021

Hydrocarbon Cold-Flow Simulation of LOX/LH2 Rocket Aerodynamics

In FY19, this project developed a test method which employed heated ethane in the physical simulation of rocket plume aerodynamics. Based on this work, EUS partnered with significant leveraged funding. A portable, self-contained test bed was also fabricated to support validation of the new method. In FY20, this project focused on the formal validation of ethane’s performance as an aerodynamic simulant as well as the development of the analytical and hardware infrastructure required to carry out that task. The test bed was activated, a series of heated-ethane tests on one or more rocket/diffuser configurations were conducted, and characterizations of the performance of the ethane simulation compared to hot-fire were performed. A series of 127 tests were conducted on 13 aerodynamic configurations to evaluate ethane’s performance as a simulant. It was found that the new methodology can reproduce rocket plume flow fields with +/- 0.5% pressure error at ~3% physical scale. The result is a cost-per-test reduction of ~99% compared to the hot-fire subscale testing typically used to obtain equivalent data. Some initial efforts have taken place, and future efforts are being planned with LaRC supersonic retro propulsion and MSFC SLS aerodynamics teams, which are interested in using the ethane system to wind-tunnel test human-rated Mars lander retrorockets and SLS Booster. Details of the prototypical experiments were published as NASA TM 2020-5009122,and have been presented in JANNAF Virtual Conference paper, "New Developments in Retro propulsion Testing for Mars Entry, Descent and Landing".

Daniel Jones↗

DOE OSTI2015

Joint neutron measurements on the flat-top HEU critical assembly machine - 14388

Flat-Top is a critical assembly with a spherical core (either ∼17 kg of HEU or ∼6 kg of δ-phase Pu) surrounded by a spherical 1000 kg natural uranium reflector. In 2013, joint LANL and CEA measurements were conducted to measure Flat-Top configurations with the HEU core using He-3 detector systems. Similar experiments were previously documented in which the Caliban and Godiva critical assemblies were measured. The core and reflector of Flat-Top have a glory hole in which nuclear material or non-nuclear material may be placed. Five configurations were measured and included subcritical configurations at various reactivity states as well as a configuration at delayed critical. The various reactivity states were achieved by varying the placement of the safety blocks and control rods (composed of natural uranium). A description of the assembly and the specific configurations measured will be presented. The subcritical and delayed critical measurements were analyzed using the Hage-Cifarelli formalism of the Feynman Variance-to-Mean method. This analysis can be used to calculate the multiplication of each configuration, from which the multiplication factor and reactivity are determined. These results are compared to the reactivity states based upon control rod worth curves and MCNP simulations. For this work, a new MCNP model was made to accurately represent the measured configurations. The ICSBEP handbook contains a benchmark of Flat-Top but all models are very simple and do not have the required fidelity to accurately represent the measured configurations. This work also examines the number of neutrons emitted per fission for this assembly (which has a significant contribution of induced fissions in U-238 due to the large natural uranium reflector). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Beyond quantum cluster theories: multiscale approaches for strongly correlated systems

The degrees of freedom that confer to strongly correlated systems their many intriguing properties also render them fairly intractable through typical perturbative treatments. For this reason, the mechanisms responsible for their technologically promising properties remain mostly elusive. Computational approaches have played a major role in efforts to fill this void. In particular, dynamical mean field theory and its cluster extension, the dynamical cluster approximation have allowed significant progress. However, despite all the insightful results of these embedding schemes, computational constraints, such as the minus sign problem in quantum Monte Carlo (QMC), and the exponential growth of the Hilbert space in exact diagonalization (ED) methods, still limit the length scale within which correlations can be treated exactly in the formalism. A recent advance aiming to overcome these difficulties is the development of multiscale many body approaches whereby this challenge is addressed by introducing an intermediate length scale between the short length scale where correlations are treated exactly using a cluster solver such QMC or ED, and the long length scale where correlations are treated in a mean field manner. At this intermediate length scale correlations can be treated perturbatively. This is the essence of multiscale many-body methods. Furthermore, we will review various implementations of these multiscale many-body approaches, the results they have produced, and the outstanding challenges that should be addressed for further advances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1995

Hybrid Upwind Splitting (HUS) by a Field-by-Field Decomposition

We introduce and develop a new approach for upwind biasing: the hybrid upwind splitting (HUS) method. This original procedure is based on a suitable hybridization of current prominent flux vector splitting (FVS) and flux difference splitting (FDS) methods. The HUS method is designed to naturally combine the respective strengths of the above methods while excluding their main deficiencies. Specifically, the HUS strategy yields a family of upwind methods that exhibit the robustness of FVS schemes in the capture of nonlinear waves and the accuracy of some FDS schemes in the resolution of linear waves. We give a detailed construction of the HUS methods following a general and systematic procedure directly performed at the basic level of the field by field (i.e. waves) decomposition involved in FDS methods. For such a given decomposition, each field is endowed either with FVS or FDS numerical fluxes, depending on the nonlinear nature of the field under consideration. Such a design principle is made possible thanks to the introduction of a convenient formalism that provides us with a unified framework for upwind methods. The HUS methods we propose bring significant improvements over current methods in terms of accuracy and robustness. They yield entropy-satisfying approximate solutions as they are strongly supported in numerical experiments. Field by field hybrid numerical fluxes also achieve fairly simple and explicit expressions and hence require a computational effort between that of the FVS and FDS. Several numerical experiments ranging from stiff 1D shock-tube to high speed viscous flows problems are displayed, intending to illustrate the benefits of the present approach. We assess in particular the relevance of our HUS schemes to viscous flow calculations.

Coquel, Frederic↗

NASA NTRS1995

Software Safety Progress in NASA

NASA has developed guidelines for development and analysis of safety-critical software. These guidelines have been documented in a Guidebook for Safety Critical Software Development and Analysis. The guidelines represent a practical 'how to' approach, to assist software developers and safety analysts in cost effective methods for software safety. They provide guidance in the implementation of the recent NASA Software Safety Standard NSS-1740.13 which was released as 'Interim' version in June 1994, scheduled for formal adoption late 1995. This paper is a survey of the methods in general use, resulting in the NASA guidelines for safety critical software development and analysis.

Radley, Charles F.↗

NASA NTRS2008

Robustness Analysis and Reliable Flight Regime Estimation of an Integrated Resilent Control System for a Transport Aircraft

Formal robustness analysis of aircraft control upset prevention and recovery systems could play an important role in their validation and ultimate certification. As a part of the validation process, this paper describes an analysis method for determining a reliable flight regime in the flight envelope within which an integrated resilent control system can achieve the desired performance of tracking command signals and detecting additive faults in the presence of parameter uncertainty and unmodeled dynamics. To calculate a reliable flight regime, a structured singular value analysis method is applied to analyze the closed-loop system over the entire flight envelope. To use the structured singular value analysis method, a linear fractional transform (LFT) model of a transport aircraft longitudinal dynamics is developed over the flight envelope by using a preliminary LFT modeling software tool developed at the NASA Langley Research Center, which utilizes a matrix-based computational approach. The developed LFT model can capture original nonlinear dynamics over the flight envelope with the ! block which contains key varying parameters: angle of attack and velocity, and real parameter uncertainty: aerodynamic coefficient uncertainty and moment of inertia uncertainty. Using the developed LFT model and a formal robustness analysis method, a reliable flight regime is calculated for a transport aircraft closed-loop system.

Shin, Jong-Yeob↗

DOE OSTI2017

Nonperturbative comparison of clover and highly improved staggered quarks in lattice QCD and the properties of the Φ meson

Here, we compare correlators for pseudoscalar and vector mesons made from valence strange quarks using the clover quark and highly improved staggered quark (HISQ) formalisms in full lattice QCD. We use fully nonperturbative methods to normalise vector and axial vector current operators made from HISQ quarks, clover quarks and from combining HISQ and clover fields. This allows us to test expectations for the renormalisation factors based on perturbative QCD, with implications for the error budget of lattice QCD calculations of the matrix elements of clover-staggered $b$-light weak currents, as well as further HISQ calculations of the hadronic vacuum polarisation. We also compare the approach to the (same) continuum limit in clover and HISQ formalisms for the mass and decay constant of the $\phi$ meson. Our final results for these parameters, using single-meson correlators and neglecting quark-line disconnected diagrams are: $m_{\phi} =$ 1.023(5) GeV and $f_{\phi} = $ 0.238(3) GeV in good agreement with experiment. These results come from calculations in the HISQ formalism using gluon fields that include the effect of $u$, $d$, $s$ and $c$ quarks in the sea with three lattice spacing values and $m_{u/d}$ values going down to the physical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2014

A Bayesian method for using simulator data to enhance human error probabilities assigned by existing HRA methods

In the past several years, several international organizations have begun to collect data on human performance in nuclear power plant simulators. The data collected provide a valuable opportunity to improve human reliability analysis (HRA), but these improvements will not be realized without implementation of Bayesian methods. Bayesian methods are widely used to incorporate sparse data into models in many parts of probabilistic risk assessment (PRA), but Bayesian methods have not been adopted by the HRA community. In this paper, we provide a Bayesian methodology to formally use simulator data to refine the human error probabilities (HEPs) assigned by existing HRA methods. We demonstrate the methodology with a case study, wherein we use simulator data from the Halden Reactor Project to update the probability assignments from the SPAR-H method. The case study demonstrates the ability to use performance data, even sparse data, to improve existing HRA methods. Furthermore, this paper also serves as a demonstration of the value of Bayesian methods to improve the technical basis of HRA.

99 GENERAL AND MISCELLANEOUS↗