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 127 records · Page 7

A multi-item maintenance center inventory model for low-demand reparable items

In many military and commercial contexts, complex equipment undergoes scheduled maintenance overhauls at regular intervals during which all failed components are replaced. Failure to have replacements on hand for all failed parts requires emergency measures at premium cost. When reparable parts are highly reliable and expensive, both holding and shortage costs are high. This model determines the reparable parts inventory for a maintenance center under three alternative criteria: (1) maximizing job-completion rate subject to constraint on total holding costs, (2) minimizing total holding costs plus expected job noncompletion costs, and (3) minimizing total holding costs subject to a required minimum job-completion rate. Exact solutions may be obtained using dynamic programming. Approximate solutions, found easily by marginal analysis, have readily computed bounds on possible error. The solution methods for the three formulations are illustrated in a simple example.

Schaefer, M. K.↗

An analysis of curvature effects for the control of wall-bounded shear flows

The Reynolds stress transport equations are used to predict the effects of simultaneous and sequential combinations of distortions on turbulent boundary layers. The equations are written in general orthogonal curvilinear coordinates, with the curvature terms expressed in terms of the principal radii of curvature of the respective coordinate surfaces. Results are obtained for the cases of two-dimensional and three-dimensional flows in the limit where production and pressure-strain redistribution dominate over diffusion effects.

Gatski, T. B.↗

Applications and accuracy of the parallel diagonal dominant algorithm

The Parallel Diagonal Dominant (PDD) algorithm is a highly efficient, ideally scalable tridiagonal solver. In this paper, a detailed study of the PDD algorithm is given. First the PDD algorithm is introduced. Then the algorithm is extended to solve periodic tridiagonal systems. A variant, the reduced PDD algorithm, is also proposed. Accuracy analysis is provided for a class of tridiagonal systems, the symmetric, and anti-symmetric Toeplitz tridiagonal systems. Implementation results show that the analysis gives a good bound on the relative error, and the algorithm is a good candidate for the emerging massively parallel machines.

Sun, Xian-He↗

A Feasibility Study to Control Airfoil Shape Using THUNDER

The objective of this study was to assess the capabilities of a new out-of-plane displacement piezoelectric actuator called thin-layer composite-unimorph ferroelectric driver and sensor (THUNDER) to alter the upper surface geometry of a subscale airfoil to enhance performance under aerodynamic loading. Sixty test conditions, consisting of combinations of five angles of attack, four dc applied voltages, and three tunnel velocities, were studied in a tabletop wind tunnel. Results indicated that larger magnitudes of applied voltage produced larger wafer displacements. Wind-off displacements were also consistently larger than wind-on. Higher velocities produced larger displacements than lower velocities because of increased upper surface suction. Increased suction also resulted in larger displacements at higher angles of attack. Creep and hysteresis of the wafer, which were identified at each test condition, contributed to larger negative displacements for all negative applied voltages and larger positive displacements for the smaller positive applied voltage (+102 V). An elastic membrane used to hold the wafer to the upper surface hindered displacements at the larger positive applied voltage (+170 V). Both creep and hysteresis appeared bounded based on the analysis of several displacement cycles. These results show that THUNDER can be used to alter the camber of a small airfoil under aerodynamic loads.

Pinkerton, Jennifer L.↗

Real-Time System Verification by Kappa-Induction

We report the first formal verification of a reintegration protocol for a safety-critical, fault-tolerant, real-time distributed embedded system. A reintegration protocol increases system survivability by allowing a node that has suffered a fault to regain state consistent with the operational nodes. The protocol is verified in the Symbolic Analysis Laboratory (SAL), where bounded model checking and decision procedures are used to verify infinite-state systems by k-induction. The protocol and its environment are modeled as synchronizing timeout automata. Because k-induction is exponential with respect to k, we optimize the formal model to reduce the size of k. Also, the reintegrator's event-triggered behavior is conservatively modeled as time-triggered behavior to further reduce the size of k and to make it invariant to the number of nodes modeled. A corollary is that a clique avoidance property is satisfied.

Pike, Lee S.↗

Traffic Flow Analysis for Package Delivery Drones using a Queueing Model

A key component of the small unmanned aircraft systems traffic management ecosystem is the design of scalable algorithms for strategic deconfliction of drones prior to takeoff. In this work, we focus on efficient flow management of drones on a network of intersecting edges subject to two kinds of spacing constraints: 1) between any two adjacent vehicles on an edge and 2) between any two vehicles on two different edges arriving one after the other at an intersection. The spacing is designed to enable non-intersection of operational volumes corresponding to two different vehicles thereby properly separating the vehicles inside each volume. For simplicity, we assume a constant ground speed for the drones and fixed dimensions for the operational volume blocks. The deconfliction is managed by adjusting the takeoff time of the drones, thereby regulating their arrival time at various crossing waypoints in the network. This framework allows us to study the maximum flow (throughput) of vehicles on a network of edges connecting depots to drop off sites subject to the temporal spacing constraints. The departure scheduling of individual drones results in a combinatorial optimization problem. To alleviate this, we solve a max-flow formulation and use queueing theory to simplify the analysis and provide upper bounds to the underlying optimization problem for individual drone departure scheduling. Our results indicate that throughput drops rapidly after the density of drones in the network passes the max-flow limits.

Alexey A Munishkin↗

Assessment of M2020 Terrain Relative Landing Accuracy: Flight Performance vs Predicts

Terrain Relative Navigation (TRN) was a critical enabling Entry, Descent, and Landing (EDL) technology that enabled Mars 2020 mission Perseverance rover to land at Jezero crater. TRN pro-vides real-time, autonomous, map-relative position determination and generates a landing target based on a priori knowledge of hazards. The required performance for TRN was to land within 60m of the selected target. The required 60m was sub-allocated to various error sources in three major categories: targeting error, knowledge error, and control error. The targeting error is the error in selecting an appropriate landing target and the knowledge of the target on the surface. It includes the Lander Vision System (LVS) position localization with respect the ground, the synchronization between the Lander Vision System measurement and the main Navigation filter, and errors associated with the LVS Reference Map and Safe Target Selec-tion (STS). The knowledge error is the contribution of knowledge growth from the synchronization with LVS to touchdown. The control error encompasses how accurately the system could stay on the desired reference trajectory. The TRN error budget uses a combination of analysis, simulation, and hardware test-ing results to bound the various error contributions obtained during the verification and validation process. This paper first presents a description the TRN system, focusing on the architecture of LVS and STS. The paper then gives detailed overview of the TRN error budget, with a description of the major error contribu-tions in each of the three categories. Next, the paper gives the results for three versions of the error budget, pre-launch, in-flight pre-landing, and post-landing. The paper compares the pre-flight analysis, the pre-landing analysis using in-flight data during cruise, to the post-landing analysis of the TRN performance. Pre-landing analysis best estimate of the landing performance was 33m, compared to the 60m require-ment. Post-landing analysis estimated a landing accuracy of 8.53m or better, much better than the 33m pre-landing estimate. The actual post-landing imagery calculated the distance of the rover to the targeted location to be 5m. The post-landing analysis closely bounds the image-based assessment of landing accu-racy, indicating the success of the error budget architecture in bounding the landing accuracy, as well as the fidelity of the simulations used to model and predict performance.

Chen, Allen↗

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗

HD 143811 AB b: A Directly Imaged Planet Orbiting a Spectroscopic Binary in Sco-Cen

We present confirmation of HD 143811 AB b, a substellar companion to spectroscopic binary HD 143811 AB through direct imaging with the Gemini Planet Imager (GPI) and Keck NIRC2. HD 143811 AB was observed as a part of the GPI Exoplanet Survey in 2016 and 2019 and is a member of the Sco-Cen star formation region. The exoplanet is detected ∼430 mas from the host star by GPI. With two GPI epochs and one from Keck/NIRC2 in 2022, we confirm through common proper motion analysis that the object is bound to its host star. We derive an orbit with a semimajor axis of $64^{+32}_{-14}$ au and eccentricity $0.23^{+0.24}_{-0.16}$. Spectral analysis of the GPI H-band spectrum and NIRC2 L′ photometry provides additional proof that this object is a substellar companion. We compare the spectrum of HD 143811 AB b to PHOENIX stellar models and Exo-Radioactive-Convective Equilibrium Model (REM) exoplanet atmosphere models and find that Exo-REM models provide the best fits to the data. From the Exo-REM models, we derive an effective temperature of $1042^{+178}_{-132}$ K for the planet and translate the derived luminosity of the planet to a mass of 5.6 ± 1.1 M Jup assuming hot-start evolutionary models. HD 143811 AB b is the first directly imaged planet around a binary that is not on an ultrawide orbit. Future characterization of this object will shed light on the formation of planets around binary star systems.

Astronomy and AstroPhysics↗

Mars Data Analysis Program

Grant NAGS12158 addressed a major NASA objective concerning the possibility of a palm ocean or large lake in the northern lowlands of Mars. Our overall approach for this study was an analysis of the graben-bounded giant polygons of Utopia Planitia, but specifically those grabens that define circles rather than open polygons. These circular grabens overlie buried impact craters, and the grabens form because of differential compaction of the overlying material over crater rims and floors. Several years ago, I predicted that the graben circles would bound depressions, and that the depths of these depressions would scale with the diameters of the graben circles. These predictions have been verified by earlier analysis. During this one-year grant (with one-year no-cost extension) we greatly increased the sample size and validated the earlier research robustly. What remained unexplained was why most of the graben circles in Utopia Planitia were double. A new model, involving volumetric compaction rather than simply 2-D compaction, satisfactorily explains the double rings and also provides a measure of relative thickness of the cover material burying the craters as a function of radial distance from the center of the Utopia Basin. Only two materials are likely candidates for the compacting cover material: volcanic ash, or wet sediment. The water in the wet sediment is largely responsible for the volumetric compaction; dry ash will compact vertically but experiences very limited lateral shrinkage. Thus the depressions within the circular grabens and the model explaining the double rings strongly favor wet sediment and thus provide evidence in favor of a past body of standing water in the northern lowlands. Publications supported entirely or in part by this grant are listed below.

McGill, George E.↗

A study of rare earth ion-adsorption clays: The speciation of rare earth elements on kaolinite at basic pH

Synthesized rare earth ion-adsorption clays (La, Ce, Nd, Dy and Y, ~ 400 ppm each) were prepared from kaolinite at various basic pH in open-air conditions. X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), and inductively coupled plasma mass spectrometry (ICP-MS) were used to characterize the effect of pH on the adsorption mechanism of rare earth elements (REEs) on kaolinite. The crystal structure of kaolinite by XRD analysis is not altered after REE adsorption, suggesting that REEs are bound to the surface of kaolinite. Elemental analysis by XPS and ICP-MS show that the surface concentration of REEs on kaolinite is pH dependent, and a local maximum of the surface concentration is achieved at pH 10. Three characteristics of the La 3d 5/2 XPS photoemission feature were measured and compared to those of standard La compounds to probe the speciation of REEs, and electrostatically bound hydroxides such as REE(OH) 2 + are identified as the primary species on the surface of kaolinite. Finally, compared to the typical ion-adsorption clays with intrinsic acidic pH, clays synthesized at basic pH exhibit similarly high rare earth ion-exchange efficiency with ammonium sulfate (AS), suggesting that rare earth ion-adsorption clays can exist in both acidic and basic environments.

36 MATERIALS SCIENCE↗

Uncertainty Quantification of Monolithic Tungsten Target Block Fatigue

The current STS spallation target design consists of 21 water cooled, tantalum clad, monolithic tungsten blocks. The target wheel assembly rotates such that each target block receives a direct proton pulse every 1.4 seconds (21 blocks / 15 Hz pulses). The quasi-steady heating induced by the repetitive beam pulses creates a significant mean stress distribution in the target. On top of this, the rapid energy deposition from a proton pulse induces dynamic stress waves that are of the same order of magnitude as the mean stresses. Over the design goal of 5,000 operational hours per year for a decade, each target must withstand over 130 million repetitions of these dynamic stress events. The fatigue life prediction of the target blocks is a complex sequence of simulations that depend upon a set of input parameters. For the tungsten block, the end goal is the fatigue factor-of-safety (FOS). This study documents a new structural analysis routine using software available to U.S. national laboratories (e.g. CUBIT, Sierra Multiphysics, Dakota) that enables more rapid target fatigue analysis. The increased throughput is amenable to parameter sensitivity and uncertainty quantification (UQ) studies. A set of input parameters (ultimate strength, fatigue strength, density, elastic modulus, Poisson’s ratio, thermal conductivity, coefficient of thermal expansion, specific heat, flow rate) with expected ranges are documented based on literature data, recent experimental testing, and engineering judgement. UQ bounds the resulting fatigue FOS at [0.5, 3.4] using epistemic interval analysis. These bounds are expected to be wider due to model effects (discretization) and beam effects (focus, offset) that were not included. Parameter sensitivity underscores the importance of ultimate strength and fatigue strength of irradiated tungsten on the fatigue calculations.

43 PARTICLE ACCELERATORS↗

Theoretical analysis of key factors achieving reversed magnetic shear q-profiles sustained with lower hybrid waves on EAST

The phase space analysis technique is applied to provide new insights into the fully non-inductive lower hybrid current driven (LHCD) discharges on EAST (Garofalo et al 2017 Nucl. Fusion 57 076037). The analysis shows that there are bounded and unbounded topologies of lower hybrid (LH) waves in phase space. For typical parameters on EAST, the propagation domain for 4.6 GHz LH waves is bounded, while unbounded for 2.45 GHz LH waves and one of the conditions is recognized to be good for achieving an off-axis current profile driven by 4.6 GHz LH waves on EAST. The parametric analysis on the potential power deposition (PPD) region for those experiments demonstrates that the reversed magnetic shear dominates the wave behavior, and confines the LH power absorption to the far off-axis region as long as it occurs in the 4.6 GHz dominated LHCD discharge. GENRAY/CQL3D simulations also confirm this effect of reversed shear on the power absorption of 4.6 GHz waves. A bounded propagation domain and a positive feedback loop between magnetic shear reversal and off-axis LHCD profile could explain the LH-only-sustained strong reversed magnetic shear observed in experiments on multiple machines. In contrast, there is no such effect for 2.45 GHz waves on EAST since the PPD region is less sensitive to the change of reversed shear due to the unbounded propagation domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A new upper bound for the growth factor in Gaussian elimination with complete pivoting

Abstract The growth factor in Gaussian elimination measures how large the entries of an LU factorization can be relative to the entries of the original matrix. It is a key parameter in error estimates, and one of the most fundamental topics in numerical analysis. We produce an upper bound of for the growth factor in Gaussian elimination with complete pivoting — the first improvement upon Wilkinson's original 1961 bound of .

Bisain, Ankit [Department of Mathematics Massachus↗

Reachability Analysis for FollowerStopper: Safety Analysis and Experimental Results

Motivated by earlier work and the developer of a new algorithm, the FollowerStopper, this article uses reachability analysis to verify the safety of the FollowerStopper algorithm, which is a controller designed for dampening stop-and-go traffic waves. With more than 1100 miles of driving data collected by our physical platform, we validate our analysis results by comparing it to human driving behaviors. The FollowerStopper controller has been demonstrated to dampen stop-and-go traffic waves at low speed, but previous analysis on its relative safety has been limited to upper and lower bounds of acceleration. To expand upon previous analysis, reachability analysis is used to investigate the safety at the speeds it was originally tested and also at higher speeds. Two formulations of safety analysis with different criteria are shown: distance-based and time headway-based. The FollowerStopper is considered safe with distance-based criterion. However, simulation results demonstrate that the FollowerStopper is not representative of human drivers - it follows too closely behind vehicles, specifically at a distance human would deem as unsafe. On the other hand, under the time headway-based safety analysis, the FollowerStopper is not considered safe anymore. A modified FollowerStopper is proposed to satisfy time-based safety criterion. Simulation results of the proposed FollowerStopper shows that its response represents human driver behavior better.

Chou, Fang-Chieh↗

Wood Products Analysis

Structural Reliability Consultants' computer program creates graphic plots showing the statistical parameters of glue laminated timbers, or 'glulam.' The company president, Dr. Joseph Murphy, read in NASA Tech Briefs about work related to analysis of Space Shuttle surface tile strength performed for Johnson Space Center by Rockwell International Corporation. Analysis led to a theory of 'consistent tolerance bounds' for statistical distributions, applicable in industrial testing where statistical analysis can influence product development and use. Dr. Murphy then obtained the Tech Support Package that covers the subject in greater detail. The TSP became the basis for Dr. Murphy's computer program PC-DATA, which he is marketing commercially.

Source record↗

Constrained model calibration of grain structure dependent spall dynamics in shock-loaded tantalum

Here, we perform a gas gun experiment by shock loading tantalum samples of varying grain structures to assess the suitability of a numerical model for simulating spall behavior. The observed differences in spall strength, as well spallation and re-compression history, are not captured in uncalibrated hydrodynamic simulations. An optimization is performed on the Johnson spall model to determine the best parameters that fit the observed trends. Linear stability analysis is employed to motivate bounds on those parameters. Herein, optimized simulations agree well with the experimental results, reproducing pullback depth and recompression timescales across the different samples tested. Further, the observed pullback time of the single crystal sample was found to imply, via the stability analysis, a percolation threshold in good agreement with the theoretical value for a body centered cubic lattice. Therefore, the combined linear stability and percolation analysis shows promise and may be applied to other materials with diverse microstructures. Collectively, the findings demonstrate that the model is suitable for reproducing spall-induced free surface behavior across various microstructures, but also points to caution in using model coefficients for uncalibrated microstructures.

58 GEOSCIENCES↗