Search NASA⌕ Search

SEARCH · Search NASA

Results for “Bounding 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 217 records · Page 12

Non-classical adhesive-bonded joints in practical aerospace construction

Solutions are derived for adhesive-bonded joints of non-classical geometries. Particular attention is given to bonded doublers and to selective reinforcement by unidirectional composites. Non-dimensionalized charts are presented for the efficiency limit imposed on the skin as the result of the eccentricity in the load path through the doubler. It is desirable to employ a relativly large doubler to minimize the effective eccentricity in the load path. The transfer stresses associated with selective reinforcement of metal structures by advanced composites are analyzed. Reinforcement of bolt holes in composites by bonded metal doublers is covered quantitatively. Also included is the adhesive joint analysis for shear flow in a multi-cell torque box, in which the bond on one angle becomes more critical sooner than those on the others, thereby restricting the strength to less than the total of each maximum strength when acting alone. Adhesive plasticity and adherend stiffness and thermal imbalances are included. A simple analysis/design technique of solution in terms of upper and lower bounds on an all-plastic adhesive analysis is introduced.

Hart-Smith, L. J.↗

Modal cost analysis for simple continua

The most popular finite element codes are based upon appealing theories of convergence of modal frequencies. For example, the popularity of cubic elements for beam-like structures is due to the rapid convergence of modal frequencies and stiffness properties. However, for those problems in which the primary consideration is the accuracy of response of the structure at specified locations, it is more important to obtain accuracy in the modal costs than in the modal frequencies. The modal cost represents the contribution of a mode in the norm of the response vector. This paper provides a complete modal cost analysis for simple continua such as beam-like structures. Upper bounds are developed for mode truncation errors in the model reduction process and modal cost analysis dictates which modes to retain in order to reduce the model for control design purposes.

Hu, A.↗

Weight Assessment for Fuselage Shielding on Aircraft With Open-Rotor Engines and Composite Blade Loss

The Federal Aviation Administration (FAA) has been engaged in discussions with airframe and engine manufacturers concerning regulations that would apply to new technology fuel efficient "openrotor" engines. Existing regulations for the engines and airframe did not envision features of these engines that include eliminating the fan blade containment systems and including two rows of counter-rotating blades. Damage to the airframe from a failed blade could potentially be catastrophic. Therefore the feasibility of using aircraft fuselage shielding was investigated. In order to establish the feasibility of this shielding, a study was conducted to provide an estimate for the fuselage shielding weight required to provide protection from an open-rotor blade loss. This estimate was generated using a two-step procedure. First, a trajectory analysis was performed to determine the blade orientation and velocity at the point of impact with the fuselage. The trajectory analysis also showed that a blade dispersion angle of 3deg bounded the probable dispersion pattern and so was used for the weight estimate. Next, a finite element impact analysis was performed to determine the required shielding thickness to prevent fuselage penetration. The impact analysis was conducted using an FAA-provided composite blade geometry. The fuselage geometry was based on a medium-sized passenger composite airframe. In the analysis, both the blade and fuselage were assumed to be constructed from a T700S/PR520 triaxially-braided composite architecture. Sufficient test data on T700S/PR520 is available to enable reliable analysis, and also demonstrate its good impact resistance properties. This system was also used in modeling the surrogate blade. The estimated additional weight required for fuselage shielding for a wing- mounted counterrotating open-rotor blade is 236 lb per aircraft. This estimate is based on the shielding material serving the dual use of shielding and fuselage structure. If the shielding material is not used for dual purpose, and is only used for shielding, then the additional weight per aircraft is estimated to be 428 lb. This weight estimate is based upon a number of assumptions that would need to be revised when applying this concept to an actual airplane design. For example, the weight savings that will result when there is no fan blade containment system, manufacturing limitations which may increase the weight where variable thicknesses was assumed, engine placement on the wing versus aft fuselage, etc.

aircraft safety↗

Performance Analysis on Fault Tolerant Control System

In a fault tolerant control (FTC) system, a parameter varying FTC law is reconfigured based on fault parameters estimated by fault detection and isolation (FDI) modules. FDI modules require some time to detect fault occurrences in aero-vehicle dynamics. In this paper, an FTC analysis framework is provided to calculate the upper bound of an induced-L(sub 2) norm of an FTC system with existence of false identification and detection time delay. The upper bound is written as a function of a fault detection time and exponential decay rates and has been used to determine which FTC law produces less performance degradation (tracking error) due to false identification. The analysis framework is applied for an FTC system of a HiMAT (Highly Maneuverable Aircraft Technology) vehicle. Index Terms fault tolerant control system, linear parameter varying system, HiMAT vehicle.

Shin, Jong-Yeob↗

Automatic Estimation of Verified Floating-Point Round-Off Errors via Static Analysis

This paper introduces a static analysis technique for computing formally verified round-off error bounds of floating-point functional expressions. The technique is based on a denotational semantics that computes a symbolic estimation of floating-point round-o errors along with a proof certificate that ensures its correctness. The symbolic estimation can be evaluated on concrete inputs using rigorous enclosure methods to produce formally verified numerical error bounds. The proposed technique is implemented in the prototype research tool PRECiSA (Program Round-o Error Certifier via Static Analysis) and used in the verification of floating-point programs of interest to NASA.

Moscato, Mariano↗

Mach 10 Stage Separation Analysis for the X43-A

This paper describes the pre-flight stage separation analysis that was conducted in support of the final flight of the X-43A. In that flight, which occurred less than eight months after the successful Mach 7 flight, the X-43A Research Vehicle attained a peak speed of Mach 9.6. Details are provided on how the lessons learned from the Mach 7 flight affected separation modeling and how adjustments were made to account for the increased flight Mach number. Also, the procedure for defining the feedback loop closure and feed-forward parameters employed in the separation control logic are described, and their effect on separation performance is explained. In addition, the range and nominal values of these parameters, which were included in the Mission Data Load, are presented. Once updates were made, the nominal pre-flight trajectory and Monte Carlo statistical results were determined and stress tests were performed to ensure system robustness. During flight the vehicle performed within the uncertainty bounds predicted in the pre-flight analysis and ultimately set the world record for airbreathing powered flight.

Tartabini, Paul V.↗

Vanilloid-dependent TRPV1 opening trajectory from cryo EM ensemble analysis

Single particle cryo-EM often yields multiple protein conformations within a single dataset, but experimentally deducing the temporal relationship of these conformers within a conformational trajectory is not trivial. Here, we use thermal titration methods and cryo-EM in an attempt to obtain temporal resolution of the conformational trajectory of the vanilloid receptor TRPV1 with resiniferatoxin (RTx) bound. Based on our cryo-EM ensemble analysis, RTx binding to TRPV1 appears to induce intracellular gate opening first, followed by selectivity filter dilation, then pore loop rearrangement to reach the final open state. This apparent conformational wave likely arises from the concerted, stepwise, additive structural changes of TRPV1 over many subdomains. Greater understanding of the RTx-mediated long-range allostery of TRPV1 could help further the therapeutic potential of RTx, which is a promising drug candidate for pain relief associated with advanced cancer or knee arthritis.

59 BASIC BIOLOGICAL SCIENCES↗

Near-Earth Phase Risk Comparison of Human Mars Campaign Architectures

A risk analysis of the launch, orbital assembly, and Earth-departure phases of human Mars exploration campaign architectures was completed as an extension of a probabilistic risk assessment (PRA) originally carried out under the NASA Constellation Program Ares V Project. The objective of the updated analysis was to study the sensitivity of loss-of-campaign risk to such architectural factors as composition of the propellant delivery portion of the launch vehicle fleet (Ares V heavy-lift launch vehicle vs. smaller/cheaper commercial launchers) and the degree of launcher or Mars-bound spacecraft element sparing. Both a static PRA analysis and a dynamic, event-based Monte Carlo simulation were developed and used to evaluate the probability of loss of campaign under different sparing options. Results showed that with no sparing, loss-of-campaign risk is strongly driven by launcher count and on-orbit loiter duration, favoring an all-Ares V launch approach. Further, the reliability of the all-Ares V architecture showed significant improvement with the addition of a single spare launcher/payload. Among architectures utilizing a mix of Ares V and commercial launchers, those that minimized the on-orbit loiter duration of Mars-bound elements were found to exceed the reliability of no spare all-Ares V campaign if unlimited commercial vehicle sparing was assumed

Manning, Ted A.↗

mu analysis with real parametric uncertainty

The authors give a broad overview, from a LFT (linear fractional transformation)/mu perspective, of some of the theoretical and practical issues associated with robustness in the presence of real parametric uncertainty, with a focus on computation. Recent results on the properties of mu in the mixed case are reviewed, including issues of NP completeness, continuity, computation of bounds, the equivalence of mu and its bounds, and some direct comparisons with Kharitonov-type analysis methods. In addition, some advances in the computational aspects of the problem, including a novel branch and bound algorithm, are briefly presented together with numerical results. The results suggest that while the mixed mu problem may have inherently combinatoric worst-case behavior, practical algorithms with modest computational requirements can be developed for problems of medium size (less than 100 parameters) that are of engineering interest.

Young, Peter M.↗

Internal Acoustics of a Pintle Valve with Supercritical Helium Flow

Large amplitude flow unsteadiness is a common phenomenon within the high flow rate ducts and valves associated with propulsion systems. Boundary layer noise, shear layers and vortex shedding are a few of the many sources of flow oscillations. The presence of lightly damped acoustic modes can organize and amplify these sources of flow perturbation, causing undesirable loading of internal parts. The present study investigates the self-induced acoustic environment within a pintle valve subject to high Reynolds Number flow of helium gas. Experiments were conducted to measure the internal pressure oscillations of the Ares I Launch Abort System (LAS) Attitude Control Motor (ACM) valve. The AGM consists of a solid propellant gas generator with eight pintle valves attached to the aft end. The pintle valve is designed to deliver variable upstream conditions to an attache( converging diverging nozzle. In order to investigate the full range of operating conditions 28 separate tests were conducted with varying pintle position and upstream pressure. Helium gas was utilized in order to closely mimic the speed of sound of the gas generator exhaust, minimizing required scaling during data analysis. The recordec pressure measurements were interrogated to multiple ends. The development of root mean square (RMS) value! versus Reynolds Number and Pintle position are important to creating bounding unsteady load curves for valve internal parts. Spectral analysis was also performed, helping to identify power spectral densities (PSD) of acoustic natural frequencies and boundary layer noise. An interesting and unexpected result was the identification of an acoustic mode within the valve which does not respond until the valve was over 60% open. Further, the response amplitude around this mode can be as large or larger than those associated with lower frequency modes.

Fishbach, Sean R.↗

An algorithm for the estimation of bounds on the emissivity and temperatures from thermal multispectral airborne remotely sensed data

The effective flux incident upon the detectors of a thermal sensor, after it has been corrected for atmospheric effects, is a function of a non-linear combination of the emissivity of the target for that channel and the temperature of the target. The sensor system cannot separate the contribution from the emissivity and the temperature that constitute the flux value. A method that estimates the bounds on these temperatures and emissivities from thermal data is described. This method is then tested with remotely sensed data obtained from NASA's Thermal Infrared Multispectral Scanner (TIMS) - a 6 channel thermal sensor. Since this is an under-determined set of equations i.e. there are 7 unknowns (6 emissivities and 1 temperature) and 6 equations (corresponding to the 6 channel fluxes), there exist theoretically an infinite combination of values of emissivities and temperature that can satisfy these equations. Using some realistic bounds on the emissivities, bounds on the temperature are calculated. These bounds on the temperature are refined to estimate a tighter bound on the emissivity of the source. An error analysis is also carried out to quantitatively determine the extent of uncertainty introduced in the estimate of these parameters. This method is useful only when a realistic set of bounds can be obtained for the emissivities of the data. In the case of water the lower and upper bounds were set at 0.97 and 1.00 respectively. Five flights were flown in succession at altitudes of 2 km (low), 6 km (mid), 12 km (high), and then back again at 6 km and 2 km. The area selected with the Ross Barnett reservoir near Jackson, Mississippi. The mission was flown during the predawn hours of 1 Feb. 1992. Radiosonde data was collected for that duration to profile the characteristics of the atmosphere. Ground truth temperatures using thermometers and radiometers were also obtained over an area of the reservoir. The results of two independent runs of the radiometer data averaged 7.03 plus or minus .70 for the first run and 7.31 plus or minus .88 for the second run. The results of the algorithm yield a temperature of 7.68 for the low altitude data to 8.73 for the high altitude data.

Jaggi, S.↗

Weighted relaxation for multigrid reduction in time

Current trends in computer architectures now mean that faster computation speed must come primarily from increased concurrency, not faster clock speeds, which are stagnating. Thus, this situation creates bottlenecks for serial algorithms, including the well-known bottleneck for sequential time-integration, where each individual time-value (i.e., time-step) is computed sequentially. One approach to alleviate this and achieve parallelism in time is with multigrid. Here, in this work, we consider multigrid-reduction-in-time (MGRIT), a multilevel method applied to the time dimension that computes multiple time-steps in parallel. Like all multigrid methods, MGRIT relies on the complementary relationship between relaxation on a fine-grid and a correction from the coarse grid to solve the problem. All current MGRIT implementations are based on unweighted-Jacobi relaxation; here we introduce the concept of weighted relaxation to MGRIT. We derive new convergence bounds for weighted relaxation, and use this analysis to guide the selection of relaxation weights. Numerical results then demonstrate that by choosing appropriate non-unitary relaxation weights, one can achieve faster convergence rates and lower iteration counts for MGRIT when compared with unweighted relaxation. In most cases, weighted relaxation yields a 10%–20% saving in iterations, which is significant when using large high-performance computers. For A-stable integration schemes, results also illustrate that under-relaxation can restore convergence in some cases where unweighted relaxation is not convergent.

97 MATHEMATICS AND COMPUTING↗

Toward reducing uncertainty quantification costs in DEM models of particulate flow: Testing simple, sensitivity-based, forward uncertainty propagation techniques

The performance of two conceptually-simple uncertainty quantification techniques are tested against the rigorous nested-loop sampling technique of Roy and Oberkampf (Comput Methods Appl Mech Eng, 200: 2131–2144, 2011) (herein called full-sampling) using two very small-scale DEM-based models of particulate flow (one gas-solid flow and one granular flow). The first simplified forward uncertainty propagation technique, reduced-sampling, uses a sensitivity analysis to eliminate uncertain inputs that have little impact on the model output prior to nested-loop sampling. The second technique, boundary-sampling, uses a sensitivity analysis to inform the selection of two bounding cases for each key model output. In conclusion, the uncertainties in the model outputs obtained via the reduced- and boundary-sampling methods agree well with those from full-sampling for both the gas-solid and granular flow models while yielding computational savings of 65–75% (reduced sampling) and 94–97% (boundary sampling).

42 ENGINEERING↗

Fractionation of Oxygen Isotopes in Uranium Oxides during Peroxide Precipitation and Dry Air Calcination

The interaction between uranium and oxygen-containing substrates is nearly ubiquitous within the nuclear fuel cycle. Given the well-known and predictable oxygen stable isotope compositions of atmospheric oxygen and meteoric waters around the world, the use of these isotopes as a potential geolocation or processing signature of uranium compounds has been of interest within the nuclear safeguards community. This study focuses on measuring the oxygen-stable isotope composition of synthetically produced uranium oxides to determine the mechanism and extent of fractionation during materials processing relevant to the nuclear fuel cycle. Here, metastudtite [UO 2 O 2 (H 2 O) 2 ] samples were first produced using water and hydrogen peroxide of a known oxygen isotope composition. This starting material was calcined in dry air at temperatures ranging from 300 °C to 1000 °C for 20 h producing oxides ranging in composition from UO 3 +x (0 ≤ x ≤ 0.5) to U 3 O 8 . In addition, calcinations between 500 and 800 °C were performed at different time intervals to evaluate the rate of isotopic equilibration. The δ 18 O values of water bound to metastudtite were measured by thermogravimetric analysis coupled to isotope ratio infrared spectroscopy (TGA–IRIS). A redesigned fluorination manifold was constructed at the University of Utah and was used to extract the bulk oxygen from precipitated and calcined products using bromine pentafluoride (BrF 5 ). Mineralization water on metastudtite (measured by TGA–IRIS) was closely associated with the aqueous environment of precipitation with a fractionation of +3.5 ± 0.6‰. In contrast, bulk oxygen from metastudtite (measured by fluorination and gas chromatography combined with isotope ratio mass spectrometry) was found to have an inconsistent fractionation with precipitation water. Samples calcined in dry air exhibited a wide range of oxygen isotope compositions which increased from δ 18 O = +1.4 ± 0.9‰ at 300 °C to +17.3 ± 0.9‰ at 1000 °C. This variation in δ 18 O values was a result of oxygen exchange with atmospheric O 2 and was found to be rapid at calcination temperatures greater than 600°C—with equilibration occurring in less than 2 h. These results provide insights into the incorporation and exchange of oxygen isotopes during precipitation and calcination processes comparable to those employed within the nuclear fuel cycle and establish a foundation for future investigations to build upon.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Structure-guided functional studies of plasmid-encoded dihydrofolate reductases reveal a common mechanism of trimethoprim resistance in Gram-negative pathogens

Two plasmid-encoded dihydrofolate reductase (DHFR) isoforms, DfrA1 and DfrA5, that give rise to high levels of resistance in Gram-negative bacteria were structurally and biochemically characterized to reveal the mechanism of TMP resistance and to support phylogenic groupings for drug development against antibiotic resistant pathogens. Preliminary screening of novel antifolates revealed related chemotypes that showed high levels of inhibitory potency against Escherichia coli chromosomal DHFR (EcDHFR), DfrA1, and DfrA5. Kinetics and biophysical analysis, coupled with crystal structures of trimethoprim bound to EcDHFR, DfrA1 and DfrA5, and two propargyl-linked antifolates (PLA) complexed with EcDHFR, DfrA1 and DfrA5, were determined to define structural features of the substrate binding pocket and guide synthesis of pan-DHFR inhibitors.

59 BASIC BIOLOGICAL SCIENCES↗

Erratum: Planck 2018 results: VI. Cosmological parameters

In the original version, the bounds given in Eqs. (87a) and (87b) on the contribution to the early-time optical depth, (15,30), contained a numerical error in deriving the 95th percentile from the Monte Carlo samples. The corrected 95% upper bounds are: τ(15,30) < 0:018 (lowE, flat τ(15, 30), FlexKnot), (1) τ(15, 30) < 0:023 (lowE, flat knot, FlexKnot): (2) These bounds are a factor of 3 larger than the originally reported results. Consequently, the new bounds do not significantly improve upon previous results from Planck data presented in Millea & Bouchet (2018) as was stated, but are instead comparable. Equations (1) and (2) give results that are now similar to those of Heinrich & Hu (2021), who used the same Planck 2018 data to derive a 95% upper bound of 0.020 using the principal component analysis (PCA) model and uniform priors on the PCA mode amplitudes.

79 ASTRONOMY AND ASTROPHYSICS↗

Toward exascale whole-device modeling of fusion devices: Porting the GENE gyrokinetic microturbulence code to GPU

GENE solves the five-dimensional gyrokinetic equations to simulate the development and evolution of plasma microturbulence in magnetic fusion devices. The plasma model used is close to first principles and computationally very expensive to solve in the relevant physical regimes. In order to use the emerging computational capabilities to gain new physics insights, several new numerical and computational developments are required. Here, we focus on the fact that it is crucial to efficiently utilize GPUs (graphics processing units) that provide the vast majority of the computational power on such systems. In this paper, we describe the various porting approaches considered and given the constraints of the GENE code and its development model, justify the decisions made, and describe the path taken in porting GENE to GPUs. We introduce a novel library called gtensor that was developed along the way to support the process. Performance results are presented for the ported code, which in a single node of the Summit supercomputer achieves a speed-up of almost 15× compared to running on central processing unit (CPU) only. Typical GPU kernels are memory-bound, achieving about 90% of peak. Our analysis shows that there is still room for improvement if we can refactor/fuse kernels to achieve higher arithmetic intensity. We also performed a weak parallel scalability study, which shows that the code runs well on a massively parallel system, but communication costs start becoming a significant bottleneck.

Germaschewski, K. (ORCID:0000000284956354)↗

Accelerating and stabilizing the convergence of vibrational self-consistent field calculations via the direct inversion of the iterative subspace (vDIIS) algorithm

The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.

Chemistry↗