Search NASA⌕ Search

SEARCH · Search NASA

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 631 records · Page 35

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.↗

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.↗

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,↗

NICMOS Spectroscopy of HD 189733b

Spectral features corresponding to methane and water opacity were reported based on transmission spectroscopy of HD 189733b with Hubble/NICMOS. Recently, these data, and a similar data set for XO-1b, have been reexamined in Gibson et al. 2010, who claim they cannot reliably reproduce prior results. We examine the methods used by the Gibson team and identify two specific issues that could act to increase the formal uncertainties and to create instability in the minimization process. This would also be consistent with the GPA10 finding that they could not identify a way to select among the several instrument models they constructed. In the case of XO-1b, the Gibson team significantly changed the way in which the instrument model is defined (both with respect to the three approaches they used for HD 189733b, and the approach used by previous authors); this change, which omits the effect of the spectrum position on the detector, makes direct intercomparison of results difficult. In the experience of our group, the position of the spectrum on the detector is an important element of the instrument model because of the significant residual structure in the NICMOS spectral flat field. The approach of changing instrument models significantly complicates understanding the data reduction process and interpreting the results. Our team favors establishing a consistent method of handling NICMOS instrument systematic errors and applying it uniformly to data sets.

infrared↗

Qualifying Bulk Metallic Glass Gear Materials for Spacecraft Applications

NASA is evaluating bulk metallic glass (BMG) gears for extreme environment (cryogenic) applications; e.g., Europa Lander. The main purpose of this report is to recommend a quality assurance (QA) protocol for the production of consistent and reliable gear castings. Currently, there are two, separate manufacturers involved; Materion Corporation produces re-melt stock and Visser Precision Cast, Inc. produces BMG castings. Division of the existing alloy specifications document into material supplier-specific documents is proposed. Until a composition-based specification is established, processability- and performance-based methodologies may be an option during manufacturing development. It is recommended that each of the documents has a distinct focus; the processing behavior of crystalline feedstock, and the mechanical behavior of amorphous castings. The report addresses the “gray area” in QA testing between these materials suppliers by outlining an equitable division of responsibilities. The use of a number of dedicated reference dies to produce witness materials for physical and mechanical property evaluation is strongly recommended. Allocation of QA tests between screen testing of gear castings, and formal testing of witness materials, including cryogenic properties, is also proposed. The most effective suite of tests for both physical and mechanical property determination are down-selected from available national and international standards. The contents of this report create a solid basis for selection of QA tests, thereby setting the stage for negotiations between the materials suppliers and the customer.This technical paper:a) provides recommendations to advance materials specifications document JPL D-56223,Rev. A from “draft” to “preliminary” status;b) proposes division of the current document into two supplier-specific documents;c) separates testing methods by “ingots/castings” and “witness materials”;d) classifies test methodology based on “processability” and “performance”;e) assembles candidate ASTM specifications for “screen testing” and “formal testing”;f) identifies standard tests compatible with evaluation of BMG flexspline gears;g) suggests division of responsibilities for processability and performance testing; andh) addresses the QA “gray area” between production of ingots and castings.The report recommends that:a) processability- and performance-based QA methods might be used in lieu ofcompositional specifications during manufacturing development only;b) reference dies be designed for “acceptance” and “qualification” castings;c) witness materials provide physical and mechanical property data;d) formal testing be conducted at ambient and cryogenic temperatures; ande) processing/properties be correlated with compositional variations.

Hales, Stephen J.↗

Granger causal inference for climate change attribution

Abstract Climate change detection and attribution (D&A) is concerned with determining the extent to which anthropogenic activities have influenced specific aspects of the global climate system. D&A fits within the broader field of causal inference, the collection of statistical methods that identify cause and effect relationships. There are a wide variety of methods for making attribution statements, each of which require different types of input data and focus on different types of weather and climate events and each of which are conditional to varying extents. Some methods are based on Pearl causality (direct experimental interference) while others leverage Granger (predictive) causality, and the causal framing provides important context for how the resulting attribution conclusion should be interpreted. However, while Granger-causal attribution analyses have become more common, there is no clear statement of their strengths and weaknesses relative to Pearl-causal attribution and no clear consensus on where and when Granger-causal perspectives are appropriate. In this prospective paper, we provide a formal definition for Granger-based approaches to trend and event attribution and a clear comparison with more traditional methods for assessing the human influence on extreme weather and climate events. Broadly speaking, Granger-causal attribution statements can be constructed quickly from observations and do not require computationally-intesive dynamical experiments. These analyses also enable rapid attribution, which is useful in the aftermath of a severe weather event, and provide multiple lines of evidence for anthropogenic climate change when paired with Pearl-causal attribution. Confidence in attribution statements is increased when different methodologies arrive at similar conclusions. Moving forward, we encourage the D&A community to embrace hybrid approaches to climate change attribution that leverage the strengths of both Granger and Pearl causality.

Risser, Mark D. (ORCID:0000000319561783)↗

Electron-Hydrogen Elastic Scattering

Scattering by single-electron systems is always of interest because the wave function of the target is known exactly. Various approximations have been employed to take into account distortion produced in the target. Among them are the method of polarized orbitals and the close coupling approximation. Recently, e-H and e-He+ S-wave scattering in the elastic region has been studied using the Feshbach projection operator formalism. In this approach, the usual Hartree-Fock and exchange potentials are augmented by an optical potential and the resulting phase shifts have rigorous lower bounds. Now this method is being applied to the e-H P-wave scattering in the elastic region. The number of terms in the Hylleraas-type wave function for the 1,3 P phase shifts is 84 and the resulting phase shifts (preliminary) are given. The results have been given up to five digits because to that accuracy they are rigorous lower bounds. They are in general agreement with the variational (VAR) results of Armstead, and those obtained from the intermediate energy R-matrix method (RM) of Scholz et al., and the finite element method (FEM) of Botero and Shertzer. The later two methods do not provide any bounds on phase shifts.

Bhatia, A. K.↗

Cosmological neutrino mass: a frequentist overview in light of DESI

We derive constraints on the neutrino mass using a variety of recent cosmological datasets, including DESI BAO, the full-shape analysis of the DESI matter power spectrum and the one-dimensional power spectrum of the Lyman-α forest (P1D) from eBOSS quasars as well as the cosmic microwave background (CMB). The constraints are obtained in the frequentist formalism by constructing profile likelihoods and applying the Feldman-Cousins prescription to compute confidence intervals. This method avoids potential prior and volume effects that may arise in a comparable Bayesian analysis. Parabolic fits to the profiles allow one to distinguish changes in the upper limits from variations in the constraining power σ of the different data combinations. We find that all profiles in the ΛCDM model are cut off by the ∑m ν ≥ 0 bound, meaning that the corresponding parabolas reach their minimum in the unphysical sector. The most stringent 95% C.L. upper limit is obtained by the combination of DESI DR2 BAO, Planck PR4 and CMB lensing at 53 meV, below the minimum of 59 meV set by the normal ordering. The corresponding constraining power σ is 43 meV, which highlights the importance of the cut-off by negative values in the determination of the upper limit. Extending ΛCDM to non-zero curvature and w 0 w a CDM relaxes the constraints past 59 meV again, but only w 0 w a CDM exhibits profiles with a minimum at a positive value. Additionally, we extend the formalism to constrain the lightest neutrino mass. For DESI DR2 BAO, Planck PR4 and CMB lensing, we find confidence limits at 20 and 19 meV for normal and inverted ordering, respectively. Using a combination of DESI DR1 full-shape, BBN and eBOSS Lyman-α P1D, we successfully constrain the neutrino mass independently of the CMB. This combination yields m l ≤ 97 and 98 meV in the normal and inverted orderings, and total neutrino mass ∑m ν ≤ 285 meV (95% C.L.). The addition of DESI full-shape or Lyman-α P1D to CMB and DESI BAO results in small but noticeable improvement of the constraining power of the data. Lyman-α free-streaming measurements especially improve the constraint. Since they are based on eBOSS data, this sets a promising precedent for upcoming DESI data.

Frequentist statistics↗

Algorithms for Performance, Dependability, and Performability Evaluation using Stochastic Activity Networks

Modeling tools and technologies are important for aerospace development. At the University of Illinois, we have worked on advancing the state of the art in modeling by Markov reward models in two important areas: reducing the memory necessary to numerically solve systems represented as stochastic activity networks and other stochastic Petri net extensions while still obtaining solutions in a reasonable amount of time, and finding numerically stable and memory-efficient methods to solve for the reward accumulated during a finite mission time. A long standing problem when modeling with high level formalisms such as stochastic activity networks is the so-called state space explosion, where the number of states increases exponentially with size of the high level model. Thus, the corresponding Markov model becomes prohibitively large and solution is constrained by the the size of primary memory. To reduce the memory necessary to numerically solve complex systems, we propose new methods that can tolerate such large state spaces that do not require any special structure in the model (as many other techniques do). First, we develop methods that generate row and columns of the state transition-rate-matrix on-the-fly, eliminating the need to explicitly store the matrix at all. Next, we introduce a new iterative solution method, called modified adaptive Gauss-Seidel, that exhibits locality in its use of data from the state transition-rate-matrix, permitting us to cache portions of the matrix and hence reduce the solution time. Finally, we develop a new memory and computationally efficient technique for Gauss-Seidel based solvers that avoids the need for generating rows of A in order to solve Ax = b. This is a significant performance improvement for on-the-fly methods as well as other recent solution techniques based on Kronecker operators. Taken together, these new results show that one can solve very large models without any special structure.

Deavours, Daniel D.↗

QCD Predictions for Physical Multimeson Scattering Amplitudes

We use lattice QCD calculations of the finite-volume spectra of systems of two and three mesons to determine, for the first time, three-particle scattering amplitudes with physical quark masses. Our results are for combinations of 𝜋 + and 𝐾 + , at a lattice spacing 𝑎 = 0.063 fm, and in the isospin-symmetric limit. We also obtain accurate results for maximal-isospin two-meson amplitudes, with those for 𝜋 + ⁢𝐾 + and 2⁢𝐾 + being the first determinations at the physical point. Dense lattice spectra are obtained using the stochastic Laplacian-Heaviside method, and the analysis leading to scattering amplitudes is done using the relativistic finite-volume formalism. Results are compared to chiral perturbation theory and to phenomenological fits to experimental data, finding good agreement.

hadron-hadron interactions↗

NASA software specification and evaluation system design, part 1

The research to develop methods for reducing the effort expended in software and verification is reported. The development of a formal software requirements methodology, a formal specifications language, a programming language, a language preprocessor, and code analysis tools are discussed.

Source record↗

Lunar and Planetary Science XXXVI, Part 18

Topics discussed include: PoDS: A Powder Delivery System for Mars In-Situ Organic, Mineralogic and Isotopic Analysis Instruments Planetary Differentiation of Accreting Planetesimals with 26Al and 60Fe as the Heat Sources Ground-based Observation of Lunar Surface by Lunar VIS/NIR Spectral Imager Mt. Oikeyama Structure: First Impact Structure in Japan? Central Mounds in Martian Impact Craters: Assessment as Possible Perennial Permafrost Mounds (Pingos) A Further Analysis of Potential Photosynthetic Life on Mars New Insight into Valleys-Ocean Boundary on Mars Using 128 Pixels per Degree MOLA Data: Implication for Martian Ocean and Global Climate Change; Recursive Topography Based Surface Age Computations for Mars: New Insight into Surficial Processes That Influenced Craters Distribution as a Step Toward the Formal Proof of Martian Ocean Recession, Timing and Probability; Laser-induced Breakdown Spectroscopy: A New Method for Stand-Off Quantitative Analysis of Samples on Mars; Milk Spring Channels Provide Further Evidence of Oceanic, >1.7-km-Deep Late Devonian Alamo Crater, Southern Nevada; Exploration of Martian Polar Residual Caps from HEND/ODYSSEY Data; Outflow Channels Influencing Martian Climate: Global Circulation Model Simulations with Emplaced Water; Presence of Nonmethane Hydrocarbons on Pluto; Difference in Degree of Space Weathering on the Newborn Asteroid Karin; Circular Collapsed Features Related to the Chaotic Terrain Formation on Mars; A Search for Live (sup 244)Pu in Deep-Sea Sediments: Preliminary Results of Method Development; Some Peculiarities of Quartz, Biotite and Garnet Transformation in Conditions of Step-like Shock Compression of Crystal Slate; Error Analysis of Remotely-Acquired Mossbauer Spectra; Cloud Activity on Titan During the Cassini Mission; Solar Radiation Pressure and Transient Flows on Asteroid Surfaces; Landing Site Characteristics for Europa 1: Topography; and The Crop Circles of Europa.

Source record↗

On Advanced Estimation Techniques for Exoplanet Detection and Characterization using Ground-based Coronagraphs

The direct imaging of planets around nearby stars is exceedingly difficult. Only about 14 exoplanets have been imaged to date that have masses less than 13 times that of Jupiter. The next generation of planet-finding coronagraphs, including VLT-SPHERE, the Gemini Planet Imager, Palomar P1640, and Subaru HiCIAO have predicted contrast performance of roughly a thousand times less than would be needed to detect Earth-like planets. In this paper we review the state of the art in exoplanet imaging, most notably the method of Locally Optimized Combination of Images (LOCI), and we investigate the potential of improving the detectability of faint exoplanets through the use of advanced statistical methods based on the concepts of the ideal observer and the Hotelling observer. We propose a formal comparison of techniques using a blind data challenge with an evaluation of performance using the Receiver Operating Characteristic (ROC) and Localization ROC (LROC) curves. We place particular emphasis on the understanding and modeling of realistic sources of measurement noise in ground-based AO-corrected coronagraphs. The work reported in this paper is the result of interactions between the co-authors during a week-long workshop on exoplanet imaging that was held in Squaw Valley, California, in March of 2012

coronagraphs↗

On Advanced Estimation Techniques for Exoplanet Detection and Characterization using Ground-Based Coronagraphs

The direct imaging of planets around nearby stars is exceedingly difficult. Only about 14 exoplanets have been imaged to date that have masses less than 13 times that of Jupiter. The next generation of planet-finding coronagraphs, including VLT-SPHERE, the Gemini Planet Imager, Palomar P1640, and Subaru HiCIAO have predicted contrast performance of roughly a thousand times less than would be needed to detect Earth-like planets. In this paper we review the state of the art in exoplanet imaging, most notably the method of Locally Optimized Combination of Images (LOCI), and we investigate the potential of improving the detectability of faint exoplanets through the use of advanced statistical methods based on the concepts of the ideal observer and the Hotelling observer. We provide a formal comparison of techniques through a blind data challenge and evaluate performance using the Receiver Operating Characteristic (ROC) and Localization ROC (LROC) curves. We place particular emphasis on the understanding and modeling of realistic sources of measurement noise in ground-based AO-corrected coronagraphs. The work reported in this paper is the result of interactions between the co-authors during a week-long workshop on exoplanet imaging that was held in Squaw Valley, California, in March of 2012.

instrumentation↗

Lessons Learned for Improving Spacecraft Ground Operations

NASA policy requires each Program or Project to develop a plan for how they will address Lessons Learned. Projects have the flexibility to determine how best to promote and implement lessons learned. A large project might budget for a lessons learned position to coordinate elicitation, documentation and archival of the project lessons. The lessons learned process crosses all NASA Centers and includes the contactor community. o The Office of The Chief Engineer at NASA Headquarters in Washington D.C., is the overall process owner, and field locations manage the local implementation. One tool used to transfer knowledge between program and projects is the Lessons Learned Information System (LLIS). Most lessons come from NASA in partnership with support contractors. A search for lessons that might impact a new design is often performed by a contractor team member. Knowledge is not found with only one person, one project team, or one organization. Sometimes, another project team, or person, knows something that can help your project or your task. Knowledge sharing is an everyday activity at the Kennedy Space Center through storytelling, Kennedy Engineering Academy presentations and through searching the Lessons Learned Information system. o Project teams search the lessons repository to ensure the best possible results are delivered. o The ideas from the past are not always directly applicable but usually spark new ideas and innovations. Teams have a great responsibility to collect and disseminate these lessons so that they are shared with future generations of space systems designers. o Leaders should set a goal for themselves to host a set numbers of lesson learned events each year and do more to promote multiple methods of lessons learned activities. o High performing employees are expected to share their lessons, however formal knowledge sharing presentation are not the norm for many employees.

Bell, Michael↗

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry. I. Effects of Line Coupling on Self-Broadened v (sub 1) and Pure Rotational Bands of NH3

The Robert-Bonamy formalism has been commonly used to calculate half-widths and shifts of spectral lines for decades. This formalism is based on several approximations. Among them, two have not been fully addressed: the isolated line approximation and the neglect of coupling between the translational and internal motions. Recently, we have shown that the isolated line approximation is not necessary in developing semi-classical line shape theories. Based on this progress, we have been able to develop a new formalism that enables not only to reduce uncertainties on calculated half-widths and shifts, but also to model line mixing effects on spectra starting from the knowledge of the intermolecular potential. In our previous studies, the new formalism had been applied to linear and asymmetric-top molecules. In the present study, the method has been extended to symmetric-top molecules with inversion symmetry. As expected, the inversion splitting induces a complete failure of the isolated line approximation. We have calculated the complex relaxation matrices of selfbroadened NH3. The half-widths and shifts in the ν1 and the pure rotational bands are reported in the present paper. When compared with measurements, the calculated half-widths match the experimental data very well, since the inapplicable isolated line approximation has been removed. With respect to the shifts, only qualitative results are obtained and discussed. Calculated off-diagonal elements of the relaxation matrix and a comparison with the observed line mixing effects are reported in the companion paper (Paper II).

Trajectory models↗

Automatic Abstraction in Planning

Traditionally, abstraction in planning has been accomplished by either state abstraction or operator abstraction, neither of which has been fully automatic. We present a new method, predicate relaxation, for automatically performing state abstraction. PABLO, a nonlinear hierarchical planner, implements predicate relaxation. Theoretical, as well as empirical results are presented which demonstrate the potential advantages of using predicate relaxation in planning. We also present a new definition of hierarchical operators that allows us to guarantee a limited form of completeness. This new definition is shown to be, in some ways, more flexible than previous definitions of hierarchical operators. Finally, a Classical Truth Criterion is presented that is proven to be sound and complete for a planning formalism that is general enough to include most classical planning formalisms that are based on the STRIPS assumption.

LOGIC PROGRAMMING↗

Parallel diffusion operator for magnetized plasmas with improved spectral fidelity

Diffusive transport processes in magnetized plasmas are highly anisotropic, with fast parallel transport along the magnetic field lines sometimes faster than perpendicular transport by orders of magnitude. This constitutes a major challenge for describing non-grid-aligned magnetic structures in Eulerian (grid-based) simulations. Here, the present paper describes and validates a new method for parallel diffusion in magnetized plasmas based on the anti-symmetry representation [Halpern and Waltz, Phys. Plasmas 25, 060703 (2018)]. In the anti-symmetry formalism, diffusion manifests as a flow operator involving the logarithmic derivative of the transported quantity. Qualitative plane wave analysis shows that the new operator naturally yields better discrete spectral resolution compared to its conventional counterpart. Numerical simulations comparing the new method against existing finite difference methods are carried out, showing significant improvement. In particular, we find that combining anti-symmetry with finite differences in diagonally staggered grids essentially eliminates the so-called “artificial numerical diffusion” that affects conventional finite difference and finite volume methods.

Anisotropic diffusion↗