Search NASA⌕ Search

SEARCH · Search NASA

Results for “Formal analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,081 records · Page 60

ALDF Data Retrieval Algorithms for Validating the Optical Transient Detector (OTD) and the Lightning Imaging Sensor (LIS)

A linear algebraic solution is provided for the problem of retrieving the location and time of occurrence of lightning ground strikes from in Advanced Lightning Direction Finder (ALDF) network. The ALDF network measures field strength, magnetic bearing, and arrival time of lightning radio emissions and solutions for the plane (i.e.. no Earth curvature) are provided that implement all of these measurements. The accuracy of the retrieval method is tested using computer-simulated data sets and the relative influence of bearing and arrival time data on the outcome of the final solution is formally demonstrated. The algorithm is sufficiently accurate to validate NASA's Optical Transient Detector (OTD) and Lightning Imaging System (LIS). We also introduce a quadratic planar solution that is useful when only three arrival time measurements are available. The algebra of the quadratic root results are examined in detail to clarify what portions of the analysis region lead to fundamental ambiguities in source location. Complex root results are shown to be associated with the presence of measurement errors when the lightning source lies near an outer sensor baseline of the ALDF network. For arbitrary noncollinear network geometries and in the absence of measurement errors, it is shown that the two quadratic roots are equivalent (no source location ambiguity) on the outer sensor baselines. The accuracy of the quadratic planar method is tested with computer-generated data sets and the results are generally better than those obtained from the three station linear planar method when bearing errors are about 2 degrees.

Koshak, W. J.↗

Fusing Quantitative Requirements Analysis with Model-based Systems Engineering

A vision is presented for fusing quantitative requirements analysis with model-based systems engineering. This vision draws upon and combines emergent themes in the engineering milieu. "Requirements engineering" provides means to explicitly represent requirements (both functional and non-functional) as constraints and preferences on acceptable solutions, and emphasizes early-lifecycle review, analysis and verification of design and development plans. "Design by shopping" emphasizes revealing the space of options available from which to choose (without presuming that all selection criteria have previously been elicited), and provides means to make understandable the range of choices and their ramifications. "Model-based engineering" emphasizes the goal of utilizing a formal representation of all aspects of system design, from development through operations, and provides powerful tool suites that support the practical application of these principles. A first step prototype towards this vision is described, embodying the key capabilities. Illustrations, implications, further challenges and opportunities are outlined.

optimization↗

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry: Line Mixing Effects in the V1 Band of NH3 - II

Line mixing effects have been calculated in the ν1 parallel band of self-broadened NH3. The theoretical approach is an extension of a semi-classical model to symmetric-top molecules with inversion symmetry developed in the companion paper [Q. Ma and C. Boulet, J. Chem. Phys. 144, 224303 (2016)]. This model takes into account line coupling effects and hence enables the calculation of the entire relaxation matrix. A detailed analysis of the various coupling mechanisms is carried out for Q and R inversion doublets. The model has been applied to the calculation of the shape of the Q branch and of some R manifolds for which an obvious signature of line mixing effects has been experimentally demonstrated. Comparisons with measurements show that the present formalism leads to an accurate prediction of the available experimental line shapes. Discrepancies between the experimental and theoretical sets of first order mixing parameters are discussed as well as some extensions of both theory and experiment.

Boulet, C.↗

Analysis of the long-term drift rates and oscillations of Jupiter’s largest vortices

Jupiter’s troposphere is dynamically characterized by a series of alternating eastward and westward jets. Embedded between these alternating jets, large vortices, characterized by their high speed closed circulations can be found in several latitudinal bands. The Great Red Spot (GRS) and the White Ovals (WOS) located at ≈22 and ≈30 degrees south respectively1 are the largest of these vortices. The GRS has undergone significant changes in size, coloration and maximum circumferential wind speed over the last two decades. The three WOS also experienced changes in size, and they merged over time until there was only one remaining, which since then has been formally known as oval BA. Here we used over four decades worth of data from the Database for Object Positions on Jupiter (JUPOS), to investigate the long-term mean longitudinal drift of these largest vortices, and to characterize the deviations from this mean drift. A linear fit to the GRS drift speed for the last two decades gives a rate of change of (3.1 ± 0.4)^(-3°)/day/year. For oval BA we find that its annual drift rate changed from -0.39°/day in 2016 to -0.58°/day in 2020 following a steep trend that is best fitted by a second order polynomial. Regarding oscillations, we find that for the last two decades the GRS has continued to oscillate longitudinally with the same period (90 days) and peak-to-peak amplitude (1.2 degrees longitude) as reported historically. For the WOS, the data prior to 2000 seems to indicate that they presented similar oscillations with periods of ≈158 days and peak-to-peak amplitudes of ≈1.6 degrees. The analysis of the data for oval BA since the year 2000 reveals evidence of an oscillation with a period of ≈138 days and peak-to-peak amplitude of ≈1.6 degrees longitude.

Raúl Morales-Juberías↗

Sensitivity Analysis for Atmospheric Infrared Sounder (AIRS) CO2 Retrieval

The Atmospheric Infrared Sounder (AIRS) is a thermal infrared sensor able to retrieve the daily atmospheric state globally for clear as well as partially cloudy field-of-views. The AIRS spectrometer has 2378 channels sensing from 15.4 micrometers to 3.7 micrometers, of which a small subset in the 15 micrometers region has been selected, to date, for CO2 retrieval. To improve upon the current retrieval method, we extended the retrieval calculations to include a prior estimate component and developed a channel ranking system to optimize the channels and number of channels used. The channel ranking system uses a mathematical formalism to rapidly process and assess the retrieval potential of large numbers of channels. Implementing this system, we identifed a larger optimized subset of AIRS channels that can decrease retrieval errors and minimize the overall sensitivity to other iridescent contributors, such as water vapor, ozone, and atmospheric temperature. This methodology selects channels globally by accounting for the latitudinal, longitudinal, and seasonal dependencies of the subset. The new methodology increases accuracy in AIRS CO2 as well as other retrievals and enables the extension of retrieved CO2 vertical profiles to altitudes ranging from the lower troposphere to upper stratosphere. The extended retrieval method for CO2 vertical profile estimation using a maximum-likelihood estimation method. We use model data to demonstrate the beneficial impact of the extended retrieval method using the new channel ranking system on CO2 retrieval.

sensitivity analysis↗

Near-axis expansion of stellarator equilibrium at arbitrary order in the distance to the axis

A direct construction of equilibrium magnetic fields with toroidal topology at arbitrary order in the distance from the magnetic axis is carried out, yielding an analytical framework able to explore the landscape of possible magnetic flux surfaces in the vicinity of the axis. This framework can provide meaningful analytical insight into the character of high-aspect-ratio stellarator shapes, such as the dependence of the rotational transform and the plasma beta limit on geometrical properties of the resulting flux surfaces. The approach developed here is based on an asymptotic expansion on the inverse aspect ratio of the ideal magnetohydrodynamics equation. The analysis is simplified by using an orthogonal coordinate system relative to the Frenet–Serret frame at the magnetic axis. The magnetic field vector, the toroidal magnetic flux, the current density, the field line label and the rotational transform are derived at arbitrary order in the expansion parameter. Moreover, a comparison with a near-axis expansion formalism employing an inverse coordinate method based on Boozer coordinates (the so-called Garren–Boozer construction) is made, where both methods are shown to agree at lowest order. Finally, as a practical example, a numerical solution using a W7-X equilibrium is presented, and a comparison between the lowest-order solution and the W7-X magnetic field is performed.

Physics↗

Bi-Kappa Proton Mirror and Cyclotron Instabilities in the Solar Wind

Abstract The charged particles in the solar wind are often observed to possess a nonthermal tail in the velocity distribution function, a feature that can be fitted with the Kappa model. The anisotropic, or bi-Kappa, model of protons, electrons, and other charged particles is thus adopted in the literature for interpreting the data as well as in the context of the analysis of wave–particle interactions. The present paper develops an approximate but efficient theory of the mirror and cyclotron instabilities excited by the bi-Kappa protons in the solar wind. A velocity moment-based quasi-linear theory of these instabilities is also formulated in order to investigate the saturation behavior. Applications of the formalism are made for instabilities close to the marginally unstable state, which is typical of the solar wind near 1 au.

79 ASTRONOMY AND ASTROPHYSICS↗

On making things the best - Aeronautical uses of optimization /Wright Bros. lecture/

The paper's purpose is to summarize and evaluate the results of an investigation into the degree to which formal optimization methods have contributed practically to the design and operation of atmospheric flight vehicles. The nature of this technology is reviewed and illustrated with simple structural examples. A series of published successful applications is described, from the fields of aerodynamics, structures, guidance and control, optimal trajectories and vehicle configuration optimization. The corresponding improvements over conventional analysis are assessed. Speculations are offered as to why these tools have made such little headway toward acceptance by designers. The growing need for their use in the future is explained; they hold out an unparalleled opportunity for improved efficiencies.

Ashley, H.↗

Discrete-time Markovian stochastic Petri nets

We revisit and extend the original definition of discrete-time stochastic Petri nets, by allowing the firing times to have a 'defective discrete phase distribution'. We show that this formalism still corresponds to an underlying discrete-time Markov chain. The structure of the state for this process describes both the marking of the Petri net and the phase of the firing time for each transition, resulting in a large state space. We then modify the well-known power method to perform a transient analysis even when the state space is infinite, subject to the condition that only a finite number of states can be reached in a finite amount of time. Since the memory requirements might still be excessive, we suggest a bounding technique based on truncation.

Ciardo, Gianfranco↗

Statistical Methods for Effective Spatial and Temporal Scaling in Support of Probabilistic Performance Assessments - 17340

Performance assessments (PAs) for the disposal of low-level radioactive waste have traditionally been performed in a deterministic manner, with inputs specified as single values as if they carry no uncertainty. Modern PAs, however, explicitly incorporate uncertainty through use of probability distributions for inputs instead of single values. While this is clearly an improvement, challenges remain in specifying appropriate distributions, many of which relate to capturing the appropriate spatial and/or temporal scale given the PA model, and matching the scale(s) of data and information to that specified by the model. More careful attention to scaling will lead to more defensible characterizations of physical variability and uncertainty in model inputs, thus leading to more defensible information regarding outcomes of interest to be used in decision making, and hence more defensible decisions. This work introduces general strategies and considerations for explicitly incorporating considerations of scale into distribution development. The strategies are explained as a series of steps and then presented with an example that starts to illustrate the new statistical concepts. This work frames the problem in terms of hierarchical statistical modeling, providing important connections to current statistical approaches to formally and rigorously account for temporal and spatial scaling. The larger goal is to build a more formal and unified approach to distribution development to help teams of modelers, subject matter experts, and statisticians work together to develop distributions and justify the choices made in the process, ultimately leading to more informed decision-making. Without this rigor, statistical distributions are not specified correctly and can have unexpected consequences for modeling and sensitivity analysis. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Tuning valley degeneracy with band inversion

Valley degeneracy is a key feature of the electronic structure that benefits the thermoelectric performance of a material. Despite recent studies which claim that high valley degeneracy can be achieved with inverted bands, our analysis of rock-salt IV–VI compounds using first-principles calculations and k•p perturbation theory demonstrates that mere band inversion is an insufficient condition for high valley degeneracy; rather, there is a critical degree to which the bands must be inverted to induce multiple carrier pockets. The so-called “band inversion parameter” is formalized as a chemically-tunable property, offering a design route to achieving high valley degeneracy in compounds with inverted bands. In conclusion, we predict that the valley degeneracy of rock-salt IV–VI compounds can be increased from N V = 4 to N V = 24, which could result in a corresponding increase in the thermoelectric figure of merit zT .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Permutation-invariant distance between atomic configurations

We present a permutation-invariant distance between atomic configurations, defined through a functional representation of atomic positions. This distance enables us to directly compare different atomic environments with an arbitrary number of particles, without going through a space of reduced dimensionality (i.e., fingerprints) as an intermediate step. Moreover, this distance is naturally invariant through permutations of atoms, avoiding the time consuming associated minimization required by other common criteria (like the root mean square distance). Finally, the invariance through global rotations is accounted for by a minimization procedure in the space of rotations solved by Monte Carlo simulated annealing. A formal framework is also introduced, showing that the distance we propose verifies the property of a metric on the space of atomic configurations. Two examples of applications are proposed. The first one consists in evaluating faithfulness of some fingerprints (or descriptors), i.e., their capacity to represent the structural information of a configuration. The second application concerns structural analysis, where our distance proves to be efficient in discriminating different local structures and even classifying their degree of similarity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analytic Gaussian covariance matrices for galaxy N-point correlation functions

Here, we derive analytic covariance matrices for the N-point correlation functions (NPCFs) of galaxies in the Gaussian limit. Our results are given for arbitrary N and projected onto the isotropic basis functions given by spherical harmonics and Wigner 3j symbols. A numerical implementation of the 4PCF covariance is compared to the sample covariance obtained from a set of lognormal simulations, Quijote dark matter halo catalogues, and MultiDark-Patchy galaxy mocks, with the latter including realistic survey geometry. The analytic formalism gives reasonable predictions for the covariances estimated from mock simulations with a periodic-box geometry. Furthermore, fitting for an effective volume and number density by maximizing a likelihood based on Kullback-Leibler divergence is shown to partially compensate for the effects of a nonuniform window function. Our result is recently shown to facilitate NPCF analysis on a realistic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-G-TeP2-07: Dosimetric Characterization of a New HDR Multi-Channel Esophageal Applicator for Brachytherapy

Purpose: To characterize the dose distribution of a new multi-channel esophageal applicator for brachytherapy HDR treatment, and particularly the effect of the presence of air or water in the applicator’s expansion balloon. Methods: A new multi-channel (6) inflatable applicator for esophageal HDR has been developed in house and tested in a simple water phantom. CT image sets were obtained under several balloon expansions (80ml of air, 50 cc of water), and channel loadings and used with the Oncentra (Elekta) planning system based on TG43 formalism. 400 cGy was prescribed to a plane 1cm away from the applicator. Planar dose distributions were measured for that plane and one next to the applicator using Gafchromic EBT3 film and scanned by a Vidar VXR-12 film digitizer. Film and TPS generated dose distributions of film were sent to OmniPro I’mRT (iba DOSIMETRY) for analysis. 2D dose profiles in both X and Y directions were compared and gamma analysis performed. Results: Film dose measurement of the air-inflated applicator is lower than the TPS calculated dose by as much as 60%. Only 80.8% of the pixels passed the gamma criteria (3%/3mm). For the water-inflated applicator, the measured film dose is fairly close to the TPS calculated dose (typically within <3%). 99.84% of the pixels passed the gamma criteria (3%/3mm). Conclusion: TG43 based calculations worked well when water was used in the expansion balloon. However, when air is present in that balloon, the neglect of heterogeneity corrections in the TG43 calculation results in large differences between calculated and measured doses. This could result in severe underdosing when used in a patient. This study illustrates the need for a TPS with an advanced algorithm which can account for heterogeneity. Supported by Innovations Department, Cleveland Clinic.

60 APPLIED LIFE SCIENCES↗

A note on sound radiation into a uniformly flowing fluid

Sound radiation generated by mechanical or vibratory source arrangements and, in some instances of planar or piston types, encompasses all magnitudes of the wavelength or frequency as well as disposition of the observation site. The comparative increase in average total radiated power from infinitesimally small or point sources of periodic strength, on passing from rest to steady rectilinear motion, depends only on Mach number in the noncompact case; and a formally analogous rise in output is linked, to the presence of a steady parallel flow past a rigid plane wall in which a compact circular piston executes normal oscillations. The perspective of effects connected with such a background is flow is widened by regard for an elongated or strip piston, which prototype permits a straightforward and general analysis.

Levine, H.↗

Topological variation in the magnetic field line at the dayside magnetopause

This paper uses normal mode analysis to examine topological changes in magnetic field lines at the magnetopause. For a modified Harris sheet, the normal modes include a symmetric solution (tearing mode) and an antisymmetric solution (twisting mode). Both modes form flux ropes. The tearing mode in the presence of shear flow is studied. By modeling the external region plasma with drift kinetic formalism, the growth rate of the tearing mode is found to be reduced by shear flow. On the other hand, the shear flow generates a Kelvin-Helmholtz instability which makes the magnetic field perturbation a twisting mode structure. When magnetopause parameters are used, it is difficult to induce a tearing mode in the presence of shear flow, though a shear flow twisting mode did grow under these conditions.

Wang, Zhi↗

Lagrangian methods in the analysis of nonlinear wave interactions in plasma

An averaged-Lagrangian method is developed for obtaining the equations which describe the nonlinear interactions of the wave (oscillatory) and background (nonoscillatory) components which comprise a continuous medium. The method applies to monochromatic waves in any continuous medium that can be described by a Lagrangian density, but is demonstrated in the context of plasma physics. The theory is presented in a more general and unified form by way of a new averaged-Lagrangian formalism which simplifies the perturbation ordering procedure. Earlier theory is extended to deal with a medium distributed in velocity space and to account for the interaction of the background with the waves. The analytic steps are systematized, so as to maximize calculational efficiency. An assessment of the applicability and limitations of the method shows that it has some definite advantages over other approaches in efficiency and versatility.

Galloway, J. J.↗

High-order gradients with the shifted boundary method: An embedded enriched mixed formulation for elliptic PDEs

Here, we propose an extension of the embedded boundary method known as “shifted boundary method” to elliptic diffusion equations in mixed form (e.g., Darcy flow, heat diffusion problems with rough coefficients, etc.). Our aim is to obtain an improved formulation that, for linear finite elements, is at least second-order accurate for both flux and primary variable, when either Dirichlet or Neumann boundary conditions are applied. Following previous work of Nishikawa and Mazaheri in the context of residual distribution methods, we consider the mixed form of the diffusion equation (i.e., with Darcy-type operators), and introduce an enrichment of the primary variable. This enrichment is obtained exploiting the relation between the primary variable and the flux variable, which is explicitly available at nodes in the mixed formulation. The proposed enrichment mimics a formally quadratic pressure approximation, although only nodal unknowns are stored, similar to a linear finite element approximation. We consider both continuous and discontinuous finite element approximations and present two approaches: a non-symmetric enrichment, which, as in the original references, only improves the consistency of the overall method; and a symmetric enrichment, which enables a full error analysis in the classical finite element context. Combined with the shifted boundary method, these two approaches are extended to high-order embedded computations, and enable the approximation of both primary and flux (gradient) variables with second-order accuracy, independently on the type of boundary conditions applied. We also show that the primary variable is third-order accurate, when pure Dirichlet boundary conditions are embedded.

97 MATHEMATICS AND COMPUTING↗