Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical Testing Approach”

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 595 records · Page 33

Hardware Verification and Validation for a Navigation Sensor Software Model in Support of Flight Vehicle Performance Analysis

… or, “It’s in the details, how to make complicated software perform like complicated hardware.” In attempts to minimize development time and quickly build an operational vehicle, NASA’s Space Launch System (SLS) has had to be intentional about integrated testing. Constraints on budget and schedule have required balance between testing needs and the desire for an integrated flight vehicle as soon as possible. To provide key insights early in design and analysis cycles, a large amount of effort has shifted into maturing and validating models at the component level with integrated testing as a means to validate their integration. In terms of SLS Navigation, this, and the model-based design approach have pushed explicit requirements for sensor models to be validated against flight hardware to high precision. This paper covers the approach taken to verify and validate the models for the two key navigation sensors on the SLS vehicle, the Redundant Inertial Navigation Sensor and the Rate Gyro Assembly. These models are used in performance evaluation, fault detection, and operations development extensively. Using a mix of data from hardware vendor documentation and testing reports, limited in-house testing, and integration activities, these models were able to be validated against flight hardware at multiple levels, from the internal software design to statistical behavior at the raw sensor and integrated box levels. The high level of insight into the hardware elements is instrumental to support flight certification activities and building confidence in SLS Navigation capability. Focused testing enabled additional insight and validation that proved invaluable and the resulting insights were used to focus and mature models. Additionally, of having validated performance-based hardware models enables a wide breadth of activities including detailed fault detection studies and integration into future vehicle frameworks, such as an upper stage and provide a valuable asset to continued SLS analysis and design.

Thomas Park↗

Temperature-dependent mechanical properties and crystal plasticity parameters for additively manufactured Haynes-214 alloy: Experiments and numerical modeling

Our experimental mechanical testing data demonstrated that the additively manufactured (AM) laser powder bed fusion (L-PBF) Haynes-214 alloy exhibits non-linear mechanical properties as the temperature rises from ambient to 870 °C. Crystal plasticity (CP) simulations provide an effective approach to gaining deeper insights into microstructure-property linkages under thermomechanical loading. This method can reduce the need for costly high-temperature mechanical testing while accounting for the effects of crystallographic texture and grain morphology on the mechanical behavior of AM materials. However, calibrating a CP model is time-consuming because individual simulations are computationally expensive and hundreds (or more) of iterations over parameter sets may be required. To address this issue, we have designed a machine learning-differential evolution (ML-DE) CP framework that can accurately interpolate the tensile properties of AM L-PBF Haynes-214 alloy across a wide temperature range from ambient to 870 °C, with minimal reliance on experimental data. The framework uses electron backscatter diffraction (EBSD) measurements to generate statistically equivalent microstructural volume elements to serve as inputs to the CP modeling framework. Stress–strain curves were generated from 1000 CP simulations, which serve as the training data set for the three ML regression algorithms explored: linear, extra-trees, and multi-layer perceptron. These three regression models were independently evaluated to compare their efficiency and identify the most suitable algorithm for the given problem. Results revealed that the extra-trees ML regressor outperforms the other models in both qualitative and quantitative aspects with an R 2 of 0.98. Subsequently, the differential evolution optimization approach is employed to calibrate the ML-based CP material parameters with experimental results obtained at various temperatures. Finally, temperature-dependent CP material parameters are formulated. The effectiveness and efficiency of the designed framework are validated through comparison with experimental results, demonstrating a high degree of agreement. These calibrated parametric constitutive equations enable further use of the CP model to study the deformation behavior of this alloy under a wide range of thermo-mechanical loading conditions.

36 MATERIALS SCIENCE↗

Advanced Signal Decomposition Analysis and Anomaly Detection in Photovoltaic Systems

With the rapid expansion of large-scale photovoltaic (PV) plants, it is paramount for solar stakeholders to understand the reliability and efficiency of their plants to inform maintenance decisions, increase production, and understand the design factors that impact performance. Diagnosing underperformance in PV plants is challenging due to the relatively few monitoring points with respect to the large geographic footprint of the plant. This work introduces a cutting-edge method that transforms the analysis and management of key factors influencing PV plant performance, including performance loss rate (PLR), recoverable soiling, and major system changes. Identifying these factors is critical for deriving actionable insights. Leveraging advanced analytical techniques such as wavelet transformation, robust regression, and extreme point analysis, this approach provides a nuanced understanding of these factors. This method has been tested across two synthetic datasets and one real dataset, consistently surpassing existing benchmarks by achieving a lower median mean absolute error and reduced error variability across all comparable components.

14 SOLAR ENERGY↗

Sensitive neutron transverse polarization analysis using a 3 He spin filter

Here we report an experimental implementation for neutron transverse polarization analysis that is capable of detecting a small angular change ($\ll$10 –3 rad) in neutron spin orientation. This approach is demonstrated for monochromatic beams, and we show that it could be extended to polychromatic neutron beams. Our approach employs a 3 He spin filter inside a solenoid with an analyzing direction perpendicular to the incident neutron polarization direction. The method was tested with polarized neutron beams and a spin rotator placed inside a μ -metal shield just upstream of the analyzer. No cryogenic superconducting shields or additional neutron spin manipulations are needed. With a counting detector, we experimentally show that the angular resolution δ θ = 1 / ( P n A N ) rad is only determined by the counting statistics for the total counts N and the product of the neutron polarization P n and the analyzing power A . With a high-flux neutron beam, 10 –6 rad angular sensitivity is feasible within a day. This simple, classical-quantum-limited transverse polarization analysis scheme may reduce the overall complexity of experimental implementation for applications requiring sensitive neutron polarimetry and improve the precision in fundamental science studies and polarized neutron imaging.

47 OTHER INSTRUMENTATION↗

Angular systematics-free cosmological analysis of galaxy clustering in configuration space

ABSTRACT Galaxy redshift surveys are subject to incompleteness and inhomogeneous sampling due to the various constraints inherent to spectroscopic observations. This can introduce systematic errors on the summary statistics of interest, which need to be mitigated in cosmological analysis to achieve high accuracy. Standard practices involve applying weighting schemes based on completeness estimates across the survey footprint, possibly supplemented with additional weighting schemes accounting for density-dependent effects. In this work, we concentrate on pure angular systematics and describe an alternative approach consisting in analysing the galaxy two-point correlation function where angular modes are nulled. By construction, this procedure removes all possible known and unknown sources of angular observational systematics, but also part of the cosmological signal. We use a modified Landy–Szalay estimator for the two-point correlation function that relies on an additional random catalogue where angular positions are randomly drawn from the galaxy catalogue, and provide an analytical model to describe this modified statistic. We test the model by performing an analysis of the full anisotropic clustering in mock catalogues of luminous red and emission-line galaxies at 0.43 < z < 1.1. We find that the model fully accounts for the modified correlation function in redshift space, without introducing new nuisance parameters. The derived cosmological parameters from the analysis of baryon acoustic oscillations and redshift-space distortions display slightly larger statistical uncertainties, mostly for the growth rate of structure parameter fσ8 that exhibits a $50{{\ \rm per\ cent}}$ statistical error increase, but free from angular systematic error.

79 ASTRONOMY AND ASTROPHYSICS↗

Novel Data Reduction Based on Statistical Similarity

Applications such as scientific simulations and power grid monitoring are generating so much data quickly that compression is essential to reduce storage requirement or transmission capacity. To achieve better compression, one is often willing to discard some repeated information. These lossy compression methods are primarily designed to minimize the Euclidean distance between the original data and the compressed data. But this measure of distance severely limits either reconstruction quality or compression performance. In this paper, we propose a new class of compression method by redefining the distance measure with a statistical concept known as exchangeability. This approach reduces the storage requirement and captures essential features, while reducing the storage requirement. In this paper, we report our design and implementation of such a compression method named IDEALEM. To demonstrate its effectiveness, we apply it on a set of power grid monitoring data, and show that it can reduce the volume of data much more than the best known compression method while maintaining the quality of the compressed data. Finally, in these tests, IDEALEM captures extraordinary events in the data, while its compression ratios can far exceed 100.

97 MATHEMATICS AND COMPUTING↗

The Role of Adjuvant Radiotherapy in the Treatment of Pleural Mesothelioma

Background: Pleural mesothelioma is a rare but aggressive form of cancer. Local recurrence represents the majority of treatment failures and overall survival (OS) outcomes remain dismal. Adding locoregional treatment with radiotherapy after surgical resection has been considered but its role remains uncertain. Objective: The purpose of this study was to evaluate the outcomes of adjuvant radiation therapy (RT) for patients with malignant pleural mesothelioma. Methods: The National Cancer Data Base (NCDB) was queried (2004–2013) for patients with malignant mesothelioma. Patients were divided into three groups: observation, surgery alone, and surgery followed by adjuvant RT. Statistics included Fisher’s exact or Chi square tests to analyze categorical proportions between groups, Kaplan–Meier analysis to evaluate OS, and Cox proportional hazards modeling to determine variables associated with OS. Propensity matching was performed to make comparisons between homogenous groups. Results: Overall, the surgery plus radiotherapy group had a higher median survival (21.4 months) compared with surgery alone (16.59 months) [p < 0.001]. RT was more likely to be delivered after extrapleural pneumonectomy than with lung-sparing surgical approaches. On multivariable analysis, receipt of surgery plus radiotherapy, chemotherapy administration, and higher socioeconomic status were associated with improved OS (p < 0.0001). After propensity matching, receipt of surgery plus radiotherapy and chemotherapy administration were still associated with improved OS (p < 0.05). Conclusions: In the treatment of malignant pleural mesothelioma, adjuvant radiotherapy after surgical intervention was associated with improved OS. This study is the largest study of adjuvant radiotherapy to date, and our findings highlight the need for additional prospective data.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Next Generation Hydrogen Storage Vessels Enabled by Carbon Fiber Infusion with a Low Viscosity, High Toughness System

Onboard hydrogen storage for fuel cell vehicles holds considerable commercial promise which could be accelerated by lower cost storage vessels. In past years, the use of lower cost reinforcing fiber had been studied as an intuitive avenue to reduce vessel costs. Although the matrix resin does not play a large role in bearing stress loads in many areas of a filament wound pressure vessel, the matrix is relied upon for stress transfer between fibers and appears to play an important role in certain regions of the vessel construction such as the “tangent region” where the cylinder section transitions to the dome. Also, the matrix resin can play a role in the damage tolerance of vessels after impact events. Furthermore, the effect of dry spots or voids in the laminate can affect stress transfer between fibers in a manner that is challenging to quantify. Therefore, excess fiber must be used to account for this variability in performance. Unfortunately, in the case hydrogen pressure vessels the cost of the vessel is strongly dependent on the amount of carbon-fiber used in the design and manufacturing. An extensive investigation on the use of non-traditional fabrication processes enabled by an ultra-low viscosity thermoset resin (<10 cP) with high fracture toughness was conducted to determine the opportunity to reduce carbon-fiber usage and the associated high cost. The hypothesis was based on enabling a lower variation in vessel performance from a reduction in void defects (often found with traditional wet-winding) and a more defect-tolerant vessel resulting from use of a matrix resin with high fracture toughness. These two features were anticipated to justify a more aggressive design approach to the vessel winding pattern. From the work in this project, the potential use of the vacuum infusion approach to produce high quality vessels was validated. However, the team recognized the need to re-optimize the design, the dry-process, and the vacuum infusion process at each scale-up stage. During the planning phase, the critical need to re-optimize at this frequency was under-estimated. Also, the team identified the lack of industry testing standards on subscale parts for providing predictive value of vessel performance. Therefore, significant effort is needed to produce statistically significant numbers of vessels. To produce such numbers of vessels the project needed in retrospect to partner closely with a series vessel manufactures to develop a proof of concept production line to make the vessels. However, without the appropriate data, a manufacturing partner was challenging to fully engage during this timeline. New technologies relying on changing raw materials and manufacturing processes carry inherent risks. This project did not reach its ultimate objectives of demonstrated a vessel with lower carbon-fiber. But some companies expressed interest in learning from the partial successes of this project for future challenges in pressure vessel applications.

08 HYDROGEN↗

The construction of life prediction models for the design of Stirling engine heater components

The service life of Stirling-engine heater structures of Fe-based high-temperature alloys is predicted using a numerical model based on a linear-damage approach and published test data (engine test data for a Co-based alloy and tensile-test results for both the Co-based and the Fe-based alloys). The operating principle of the automotive Stirling engine is reviewed; the economic and technical factors affecting the choice of heater material are surveyed; the test results are summarized in tables and graphs; the engine environment and automotive duty cycle are characterized; and the modeling procedure is explained. It is found that the statistical scatter of the fatigue properties of the heater components needs to be reduced (by decreasing the porosity of the cast material or employing wrought material in fatigue-prone locations) before the accuracy of life predictions can be improved.

Petrovich, A.↗

Damage detection in elastic structures using vibratory residual forces and weighted sensitivity

A methodology is presented for detecting structural damage in elastic structures by nondestructive means. Measured modal test data along with a correlated analytical structural model are used to locate potentially damaged regions using residual modal force vectors and to conduct a weighted sensitivity analysis to assess the extent of mass and/or stiffness variations, where damage is characterized as a stiffness reduction. The current approach is unique among other approaches in that it accounts for (1) variations in system mass, system stiffness, and mass center (locations), (2) perturbations of both the natural frequencies and modal vectors, and (3) statistical confidence factors for the structural parameters and potential experimental instrumentation error. Moreover, this procedure can be used with either full or reduced models. A wide variety of numerical examples are presented that show that the current method provides a precise indication of both the location and the extent of structural damage.

Ricles, J. M.↗

A GPS-Based Pitot-Static Calibration Method Using Global Output-Error Optimization

Pressure-based airspeed and altitude measurements for aircraft typically require calibration of the installed system to account for pressure sensing errors such as those due to local flow field effects. In some cases, calibration is used to meet requirements such as those specified in Federal Aviation Regulation Part 25. Several methods are used for in-flight pitot-static calibration including tower fly-by, pacer aircraft, and trailing cone methods. In the 1990 s, the introduction of satellite-based positioning systems to the civilian market enabled new inflight calibration methods based on accurate ground speed measurements provided by Global Positioning Systems (GPS). Use of GPS for airspeed calibration has many advantages such as accuracy, ease of portability (e.g. hand-held) and the flexibility of operating in airspace without the limitations of test range boundaries or ground telemetry support. The current research was motivated by the need for a rapid and statistically accurate method for in-flight calibration of pitot-static systems for remotely piloted, dynamically-scaled research aircraft. Current calibration methods were deemed not practical for this application because of confined test range size and limited flight time available for each sortie. A method was developed that uses high data rate measurements of static and total pressure, and GPSbased ground speed measurements to compute the pressure errors over a range of airspeed. The novel application of this approach is the use of system identification methods that rapidly compute optimal pressure error models with defined confidence intervals in nearreal time. This method has been demonstrated in flight tests and has shown 2- bounds of approximately 0.2 kts with an order of magnitude reduction in test time over other methods. As part of this experiment, a unique database of wind measurements was acquired concurrently with the flight experiments, for the purpose of experimental validation of the optimization method. This paper describes the GPS-based pitot-static calibration method developed for the AirSTAR research test-bed operated as part of the Integrated Resilient Aircraft Controls (IRAC) project in the NASA Aviation Safety Program (AvSP). A description of the method will be provided and results from recent flight tests will be shown to illustrate the performance and advantages of this approach. Discussion of maneuver requirements and data reduction will be included as well as potential applications.

Foster, John V.↗

Wind-Tunnel Balance Characterization for Hypersonic Research Applications

Wind-tunnel research was recently conducted at the NASA Langley Research Center s 31-Inch Mach 10 Hypersonic Facility in support of the Mars Science Laboratory s aerodynamic program. Researchers were interested in understanding the interaction between the freestream flow and the reaction control system onboard the entry vehicle. A five-component balance, designed for hypersonic testing with pressurized flow-through capability, was used. In addition to the aerodynamic forces, the balance was exposed to both thermal gradients and varying internal cavity pressures. Historically, the effect of these environmental conditions on the response of the balance have not been fully characterized due to the limitations in the calibration facilities. Through statistical design of experiments, thermal and pressure effects were strategically and efficiently integrated into the calibration of the balance. As a result of this new approach, researchers were able to use the balance continuously throughout the wide range of temperatures and pressures and obtain real-time results. Although this work focused on a specific application, the methodology shown can be applied more generally to any force measurement system calibration.

Lynn, Keith C.↗

Covariance Structure of Pu Isotopic Compositions Obtained Using FRAM 5.2

In the nondestructive assay of special nuclear materials, isotopic determination based on the analysis of high-purity Ge high-resolution gamma-ray spectra plays a pivotal role. The isotopic results, which are based on a relative measurement of peak intensities, are used to interpret data acquired using quantitative (mass) nondestructive assay methods. In the case of plutonium measurements, this includes both passive and active neutron techniques, and also calorimetry (heat). Historically, precision estimates on the Pu- isotopic composition have been generated, on a case-by-case basis, by least squares fitting algorithms embedded in the analysis codes, but all other sources of uncertainty have generally been estimated using a top-down approach by item category. A largely neglected consideration is the covariance between the reported abundances and derived quantities – notably Pu-239-effective mass, Pu-240-effective mass and specific thermal power. In this work we study the covariance structure by statistically evaluating repeat spectrum measurements of two known low burn-up heavily filtered Pu items reduced using the spectrum analysis code FRAM 5.2 and alternate standard parameter sets. For routine measurements, however, only a single spectrum from an item will typically be acquired. We therefore discuss how to generate the sought after covariance information from an individual spectrum by spectrum re-sampling, an approach which could be added to provide item specific reports, and compare the two approaches to test its reliability.

Croft, Stephen↗

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗

Comprehensive Measurement and Analysis of the User-Perceived I/O Performance in a Production Leadership-Class Storage System

With the increase of the scale and intensity of the parallel I/O workloads generated by those scientific applications running on high performance computing facilities, understanding the I/O dynamics, especially the root cause of the I/O performance variability and degradation in HPC environment, have become extremely critical to the HPC community. In this paper, we run extensive I/O measuring tests on a production leadership-class storage system to capture the performance variabilities of large-scale parallel I/O. Analyzing these results and its statistic correlation revealed some valuable insights into the characteristics of the storage system and the root cause of I/O performance variability. Further, we leverage these findings and propose an I/O middleware design refactoring which can improve the performance of the parallel I/O by optimizing the data striping and placement. Our preliminary evaluation results demonstrate the proposed approach can reduce the average per-process write latency by at least 80% and the maximum per-process write latency by at least 20%.

Wan, Lipeng↗

Development and in-flight performance of the Mariner 9 spacecraft propulsion system

On November 14, 1971, Mariner 9 was decelerated into orbit about Mars by a 1334 N (300 lbf) liquid bipropellant propulsion system. This paper describes and summarizes the development and in-flight performance of this pressure-fed, nitrogen tetroxide/monomethyl hydrazine bipropellant system. The design of all Mariner propulsion subsystems has been predicted upon the premise that simplicity of approach, coupled with thorough qualification and margin-limits testing, is the key to cost-effective reliability. The qualification test program and analytical modeling are also discussed. Since the propulsion subsystem is modular in nature, it was completely checked, serviced, and tested independent of the spacecraft. Proper prediction of in-flight performance required the development of three significant modeling tools to predict and account for nitrogen saturation of the propellant during the six-month coast period and to predict and statistically analyze in-flight data.

Evans, D. D.↗

Explicit 3D continuum fracture modeling with smooth particle hydrodynamics

Impact phenomena shaped our solar system. As usual for most solar system processes, the scales are far different than we can address directly in the laboratory. Impact velocities are often much higher than we can achieve, sizes are often vastly larger, and most impacts take place in an environment where the only gravitational force is the mutual pull of the impactors. The Smooth Particle Hydrodynamics (SPH) technique has been applied in the past to the simulations of giant impacts. In these simulations, the colliding objects were so massive (at least a sizeable fraction of the Earth's mass) that material strength was negligible compared to gravity. This assumption can no longer be made when the bodies are much smaller. To this end, we have developed a 3D SPH code that includes a strength model to which we have added a von Mises yielding relation for stresses beyond the Hugoniot Elastic Limit. At the lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution. Following Grady and Kipp and Melosh et al., we introduce a state variable D ('damage'), 0 less than D less than 1, which expresses the local reduction in strength due to crack growth under tensile loading. Unfortunately for the hydrodynamics, Grady and Kipp's model predicts which fragments are the most probable ones and not the ones that are really formed. This means, for example, that if a given laboratory experiment is modeled, the fragment distribution obtained from the Grady-Kipp theory would be equivalent to a ensemble average over many realizations of the experiment. On the other hand, the hydrodynamics itself is explicit and evolves not an ensemble average but very specific fragments. Hence, there is a clear incompatibility with the deterministic nature of the hydrodynamics equations and the statistical approach of the Grady-Kipp dynamical fracture model. We remedy these shortcomings by making the incipient flaw distribution explicit, i.e., particles carry activation strains which are distributed at random with a probability of occurrence given by the Weibull distribution. If the local principal axis strain exceeds this limit, damage starts to grow. By growing explicit cracks together with statistical cracks (damage) at the sub-particle scale, we ensure that material strength and fragmentation is independent of model resolution. We tested our scheme by simulating laboratory impact experiments on basalt spheres.

Benz, W.↗

Software Reliability 2002

In FY01 we learned that hardware reliability models need substantial changes to account for differences in software, thus making software reliability measurements more effective, accurate, and easier to apply. These reliability models are generally based on familiar distributions or parametric methods. An obvious question is 'What new statistical and probability models can be developed using non-parametric and distribution-free methods instead of the traditional parametric method?" Two approaches to software reliability engineering appear somewhat promising. The first study, begin in FY01, is based in hardware reliability, a very well established science that has many aspects that can be applied to software. This research effort has investigated mathematical aspects of hardware reliability and has identified those applicable to software. Currently the research effort is applying and testing these approaches to software reliability measurement, These parametric models require much project data that may be difficult to apply and interpret. Projects at GSFC are often complex in both technology and schedules. Assessing and estimating reliability of the final system is extremely difficult when various subsystems are tested and completed long before others. Parametric and distribution free techniques may offer a new and accurate way of modeling failure time and other project data to provide earlier and more accurate estimates of system reliability.

Wallace, Dolores R.↗