Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods: data 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,675 records · Page 93

Bayesian Estimation of Thermonuclear Reaction Rates for Deuterium+Deuterium Reactions

The study of d+d reactions is of major interest since their reaction rates affect the predicted abundances of D, 3 He, and 7 Li. In particular, recent measurements of primordial D/H ratios call for reduced uncertainties in the theoretical abundances predicted by Big Bang nucleosynthesis (BBN). Different authors have studied reactions involved in BBN by incorporating new experimental data and a careful treatment of systematic and probabilistic uncertainties. To analyze the experimental data, Coc et al. used results of ab initio models for the theoretical calculation of the energy dependence of S -factors in conjunction with traditional statistical methods based on χ 2 minimization. Bayesian methods have now spread to many scientific fields and provide numerous advantages in data analysis. Astrophysical S -factors and reaction rates using Bayesian statistics were calculated by Iliadis et al. In this paper we present a similar analysis for two d+d reactions, d(d, n) 3 He and d(d, p) 3 H, that has been translated into a total decrease of the predicted D/H value by 0.16%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

When will low-contrast features be visible in a STEM X-ray spectrum image?

When will a small or low-contrast feature, such as an embedded second-phase particle, be visible in a scanning transmission electron microscopy (STEM) X-ray map? This work illustrates a computationally inexpensive method to simulate X-ray maps and spectrum images (SIs), based upon the equations of X-ray generation and detection. To particularize the general procedure, an example of nanostructured ferritic alloy (NFA) containing nm-sized Y 2 Ti 2 O 7 embedded precipitates in ferritic stainless steel matrix is chosen. The proposed model produces physically appearing simulated SI data sets, which can either be reduced to X-ray dot maps or analyzed via multivariate statistical analysis. Comparison to NFA X-ray maps acquired using three different STEM instruments match the generated simulations quite well, despite the large number of simplifying assumptions used. A figure of merit of electron dose multiplied by X-ray collection solid angle is proposed to compare feature detectability from one data set (simulated or experimental) to another. The proposed method can scope experiments that are feasible under specific analysis conditions on a given microscope. As a result, future applications, such as spallation proton–neutron irradiations, core-shell nanoparticles, or dopants in polycrystalline photovoltaic solar cells, are proposed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Variational inference of effective range parameters for 3 He− 4 He scattering

We use two different methods, Monte Carlo sampling and variational inference (VI), to perform a Bayesian calibration of the effective-range parameters in 3 He– 4 He elastic scattering. The parameters are calibrated to data from a recent set of 3 He– 4 He elastic scattering differential cross section measurements. Analysis of these data for E lab ≤ 4.3 MeV yields a unimodal posterior for which both methods obtain the same structure. However, the effective-range expansion amplitude does not account for the 7/2 − state of 7 Be so, even after calibration, the description of data at the upper end of this energy range is poor. The data up to E lab = 2.6 MeV can be well described, but calibration to this lower-energy subset of the data yields a bimodal posterior. After adapting VI to treat such a multi-modal posterior we find good agreement between the VI results and those obtained with parallel-tempered Monte Carlo sampling.

effective field theory↗

SU-E-T-153: Establish a Comprehensive Patient-Specific Plan QA Database for Instituitional Quality Control Program

Purpose: Since the publication of TG-119, several new treatment techniques requiring new QA devices have been implemented in the clinic. To monitor and analyze our institutional QA performance, we have created a comprehensive QA database using Research Electronic Data Capture (REDCap). The database will also assist us in creating IMRT QA Analysis practice guidelines. Methods: Since Oct. 2013, 336 patient plan QA results were entered into the database. Plan parameters such as plan ID, treatment site, technique, energy, optimization constraints, modulation factor, leaf speed, leaf opening.etc were automatically extracted from the Varian Eclipse database to allow us to refine our evaluation and analysis method. Specific QA device, LINAC-related information, and measurement and analysis results were manually entered by the QA team. IMRT plans were measured using MapCHECK2 while RapidArc plans were measured using ArcCHECK. Distance-To-Agreement 3%/3mm without global maximum normalization was used. Results: The data indicates that different treatment techniques might benefit from a different site-specific action level(AL) depending on the complexity of the plan and optimization parameter used (e.g., breast IMRT QA= 97.8% and pelvis IMRT QA=93.1%). Different QA devices may also benefit from a different AL (MapCHECK2 = 94.1% while ArcCHECK = 83.0%). The relationship between the parameters and passing rate suggests that the complexity of each plan, characterized by leaf travel, leaf opening and modulation factor, affect the passing rate significantly. The database is reviewed regularly, and any abnormal point of the QA result or a trend of lower QA passing rate on a specific LINAC is further investigated. Conclusion: Establishing a comprehensive QA database provides an overview of the quality assurance program. It not only helps in answering the question “what is a reasonable and achievable standard for each institution”, but also saves time in monitoring and investigating the daily machine and QA performance.

07 ISOTOPE AND RADIATION SOURCES↗

Unveiling acoustic physics of the CMB using nonparametric estimation of the temperature angular power spectrum for Planck

Estimation of the angular power spectrum is one of the important steps in Cosmic Microwave Background (CMB) data analysis. Here, we present a nonparametric estimate of the temperature angular power spectrum for the Planck 2013 CMB data. The method implemented in this work is model-independent, and allows the data, rather than the model, to dictate the fit. Since one of the main targets of our analysis is to test the consistency of the ΛCDM model with Planck 2013 data, we use the nuisance parameters associated with the best-fit ΛCDM angular power spectrum to remove foreground contributions from the data at multipoles ℓ ≥50. We thus obtain a combined angular power spectrum data set together with the full covariance matrix, appropriately weighted over frequency channels. Our subsequent nonparametric analysis resolves six peaks (and five dips) up to ℓ ∼1850 in the temperature angular power spectrum. We present uncertainties in the peak/dip locations and heights at the 95% confidence level. We further show how these reflect the harmonicity of acoustic peaks, and can be used for acoustic scale estimation. Based on this nonparametric formalism, we found the best-fit ΛCDM model to be at 36% confidence distance from the center of the nonparametric confidence set—this is considerably larger than the confidence distance (9%) derived earlier from a similar analysis of the WMAP 7-year data. Another interesting result of our analysis is that at low multipoles, the Planck data do not suggest any upturn, contrary to the expectation based on the integrated Sachs-Wolfe contribution in the best-fit ΛCDM cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectral tomographic analysis of Bremsstrahlung X-rays generated in a laser-produced plasma

A new approach is proposed to analyze Bremsstrahlung X-rays that are emitted from laser-produced plasmas (LPP) and are measured by a stack type spectrometer. This new method is based on a spectral tomographic reconstruction concept with the variational principle for optimization, without referring to the electron energy distribution of a plasma. This approach is applied to the analysis of some experimental data obtained at a few major laser facilities to demonstrate the applicability of the method. Furthermore, slope temperatures of X-rays from LPP are determined with a two-temperature model, showing different spectral characteristics of X-rays depending on laser properties used in the experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tomographic analysis of neutron and gamma pulse shape distributions from liquid scintillation detectors at Joint European Torus

The Joint European Torus (JET, Culham, UK) is the largest tokamak in the world devoted to nuclear fusion experiments of magnetic confined Deuterium (D)/Deuterium-Tritium (DT) plasmas. Neutrons produced in these plasmas are measured using various types of neutron detectors and spectrometers. Two of these instruments on JET make use of organic liquid scintillator detectors. The neutron emission profile monitor implements 19 liquid scintillation counters to detect the 2.45 MeV neutron emission from D plasmas. A new compact neutron spectrometer is operational at JET since 2010 to measure the neutron energy spectra from both D and DT plasmas. Liquid scintillation detectors are sensitive to both neutron and gamma radiation but give light responses of different decay time such that pulse shape discrimination techniques can be applied to identify the neutron contribution of interest from the data. The most common technique consists of integrating the radiation pulse shapes within different ranges of their rising and/or trailing edges. In this article, a step forward in this type of analysis is presented. The method applies a tomographic analysis of the 3-dimensional neutron and gamma pulse shape and pulse height distribution data obtained from liquid scintillation detectors such that n/γ discrimination can be improved to lower energies and additional information can be gained on neutron contributions to the gamma events and vice versa.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Round Robin Summary for ASTM Method C1733

This report was prepared to summarize the results of an inter- and intra-laboratory round robin of ASTM Method C1733-17 (now C1733-20a) Distribution Coefficients of Inorganic Species by the Batch Method. The testing used the sorption of cesium (Cs) to a standard reference soil from the San Joaquin Valley (Standard Reference Material 2709a). The data collected from the nine participants was used for a statistical analysis to assess the precision of the method and generate a precision and bias statement to update the C1733-20a standard.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Monte Carlo sampling and correlation coefficients of critical experiments

We study the influence of calculation parameters on Monte-Carlo generated correlation coefficients of k{sub eff} values for a data set of critical experiments. The correlation coefficients can be used e.g. to reduce conservatism in bias estimations, but generally must be treated correctly if a data set includes correlations. We find that the actual value of the correlation coefficient for a given data set depends on the number of samples and the ratio of the sampling uncertainty to the Monte-Carlo uncertainty of the individual k{sub eff} calculation. We also discuss that the number of input parameters of the model to calculate the k{sub eff} values do not necessarily impact the convergence of the correlation coefficient with increasing sample number. However, we explain how varying the aforementioned parameters lead to different converged results for the correlation coefficients. The conclusion is that a correlation coefficient should not be discussed isolated but with the underlying uncertainties of both, the k{sub eff} estimation and the Monte Carlo sampling approach. The underlying modelling approach of the experimental data has to be described in detail. This should include a sensitivity analysis of the impact of input parameters on the k{sub eff} uncertainty with respect to the modelling approach.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An In Situ , Automated High-Explosives Aging Method Utilizing Two-Dimensional Gas Chromatography–Mass Spectrometry

Understanding chemical changes that occur in high explosives as they age is of great importance to the safe employment and storage of these compounds. Traditional methods of aging high explosives even under accelerated aging conditions are time intensive with durations on the order of months to years. The nature of traditional aging analyses reduces each sample to a snapshot data point often separated widely in time, requiring many assumptions as to how the degradation products develop. Further complicating matters, several analytical techniques are typically employed for each sample analysis in order to ascertain an entire picture of the decomposition pathways. To address these shortcomings with existing methods, a new method of accelerated aging of high explosives utilizing comprehensive two-dimensional gas chromatography coupled to high-resolution mass spectrometry (GC × GC-HRMS) was developed using 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) as a model compound for method development. This in situ automated method reduces the time scale of aging to a matter of hours using the inlet of the GC × GC as the aging vessel. GC × GC in combination with HRMS allowed for the collection of both evolved gases and other decomposition products produced during the entire aging process in real time with HRMS providing far greater certainty in identification of explosives aging products. Additionally, this method allowed for a higher throughput of samples with greatly simplified sample preparation. Chemometric analysis of the GC × GC-HRMS data set via the alteration analysis (ALA) enabled discovery of statistically significant chemical changes providing insight into the variation of decomposition pathways with varying aging temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Block-diagonalization as a tool for the robust diabatization of high-dimensional potential energy surfaces

Robust diabatization techniques are key for the development of high-dimensional coupled potential energy surfaces (PESs) to be used in multi-state quantum dynamics simulations. In the present study we demonstrate that, besides the actual diabatization technique, common problems with the underlying electronic structure calculations can be the reason why a diabatization fails. After giving a short review of the theoretical background of diabatization, we propose a method based on the block-diagonalization to analyse the electronic structure data. This analysis tool can be used in three different ways: First, it allows to detect issues with the ab initio reference data and is used to optimize the setup of the electronic structure calculations. Second, the data from the block-diagonalization are utilized for the development of optimal parametrized diabatic model matrices by identifying the most significant couplings. Third, the block-diagonalization data are used to fit the parameters of the diabatic model, which yields an optimal initial guess for the non-linear fitting required by standard or more advanced energy based diabatization methods. The new approach is demonstrated by the diabatization of 9 electronic states of the propargyl radical, yielding fully coupled full-dimensional (12D) PESs in closed form.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗

Performance and safety parameters for the high flux isotope reactor

A Monte Carlo depletion model for the High Flux Isotope Reactor (HFIR) Cycle 400 and its use in calculating parameters of relevance to the reactor performance and safety during the reactor cycle are presented in this paper. This depletion model was developed to serve as a reference for the design of a low-enriched uranium (LEU) fuel for an ongoing study to convert HFIR from high-enriched uranium (HEU) to LEU fuel; both HEU and LEU depletion models use the same methodology and ENDF/B-VII nuclear data as discussed in this paper. The calculated HFIR Cycle 400 parameters, which are compared with measurement data from critical experiments performed at HFIR, data included in the HFIR Safety Analysis Report (SAR), or data reported by previous calculations, provide a basis for verification or updating of the corresponding SAR data. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Second derivative analysis and alternative data filters for multi-dimensional spectroscopies: A Fourier-space perspective

The second derivative image (SDI) method is uniformly applied to sharpen dispersive data features in multi-dimensional spectroscopies such as angle resolved photoemission spectroscopy (ARPES). In this work, the SDI function is represented in Fourier space, where it has the form of a multi-band pass filter. The interplay of the SDI procedure with undesirable noise and background features in ARPES data sets is reviewed, and it is shown that final image quality can be improved by eliminating higher Fourier harmonics of the SDI filter. Furthermore, we discuss extensions of SDI-like band pass filters to higher dimensional data sets, and how one can create even more effective filters with some a priori knowledge of the spectral features.

36 MATERIALS SCIENCE↗

Extended maximum likelihood halo-independent analysis of dark matter direct detection data

We extend and correct a recently proposed maximum-likelihood halo-independent method to analyze unbinned direct dark matter detection data. Instead of the recoil energy as independent variable we use the minimum speed a dark matter particle must have to impart a given recoil energy to a nucleus. This has the advantage of allowing us to apply the method to any type of target composition and interaction, e.g. with general momentum and velocity dependence, and with elastic or inelastic scattering. We prove the method and provide a rigorous statistical interpretation of the results. As first applications, we find that for dark matter particles with elastic spin-independent interactions and neutron to proton coupling ratio f{sub n}/f{sub p}=−0.7, the WIMP interpretation of the signal observed by CDMS-II-Si is compatible with the constraints imposed by all other experiments with null results. We also find a similar compatibility for exothermic inelastic spin-independent interactions with f{sub n}/f{sub p}=−0.8.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Covariance of PFNS Results from the Chi-Nu Experiment

The prompt fission neutron spectrum (PFNS) from neutron-induced fission is a fundamental quantity for the behavior of nuclear reactors, and has been measured many times on a wide variety of nuclei and covering different ranges of incident and emitted neutron energies. However, results from past measurements are frequently called into question in modern nuclear data evaluations because of a lack of thorough experimental documentation and incomplete uncertainty analyses. The Chi-Nu experiment at Los Alamos National Laboratory was designed to produce high-precision measurements of the PFNS of major actinides over a wide range of incident and emitted neutron energies, and with the documentation and covariance analysis required to ensure that the results of this experiment maintain their impact long into the future, thereby avoiding this pitfall of past measurements. In this work we describe the Chi-Nu experiment along with summaries of the treatment of and methods developed to address two important components of the analysis of Chi-Nu data: random-coincidence backgrounds and MCNP simulations. Furthermore, we describe the first results for correlations not just between all data points collected on a single target nucleus, but also between all data points from separate Chi-Nu measurements on 235 U and 239 Pu. These correlations are important for accurately calculating ratios of the PFNS from one actinide to another, which are rare and can be informative for nuclear data evaluation efforts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessment of Accelerated Stress Testing Data for Silicon Photovoltaics Using Tensor Decomposition Methods

In this work, we examine the use of high-order tensor decompositions to analyze degradation pathways emerging from accelerated stress testing of silicon photovoltaic (PV) modules. Matrix-based decompositions are powerful tools for studying two-dimensional data arrays and form the foundation of a host of classical data analysis techniques. Tensors are high-order extrapolations of matrices that are able to account for more parameter dimensions, and a variety of tensor decomposition methods have been developed that similarly seek to extend insights from matrix decompositions to higher dimensions. Applying and interpreting tensor decomposition methods to sequences of PV module image data, we seek to uncover and isolate different degradation modes occurring from accelerated stress testing procedures. Further, we consider the contributions of different modes to PV module performance degradations.

data analysis↗

Scalable, In-situ Data Clustering Data Analysis for Extreme Scale Scientific Computing (Final Report)

The objective of this project is to address challenges in the design and development of scalable in-situ data clustering and analytics algorithms and software. Our goal is to develop parallel software consisting of a set of spatio-temporal data clustering and anomaly detection functions, both of which are very important for large-scale analysis and have wide applicability for in-situ runs as well as post-processing analysis. Our design principles for in-situ analysis consider the following: (1) identify parts of the computation can be done close to the data within the nodes, while it is still in memory; (2) extract analysis components can (and should) be performed in remote staging and analysis nodes; (3) develop error-bound approximation methods for applications tolerable for small errors; (4) identify the type of derived distributions and statistics, for spatio-temporal data, that can be kept locally in order to both accelerate computations and meet energy constraints in subsequent iterations and phases; (5) use a self-describing data format so that data can be consistent and understood among local storage (memory and SSDs) and at staging and analysis nodes, thereby providing portability and flexibility; (6) develop service-oriented functions that can schedule in-situ and post-hoc analysis tasks based on the dynamic requirements of applications. Our development focus is to produce the parallel data analysis software/library that will be scalable, reusable, extensible, and generic for applications in different disciplines. The software will be able to run in-situ with the simulations as well as post-hoc analysis. This approach will satisfy many synergistic requirements for data intensive applications executed on data coming from instruments and experiments. In particular, the proposed multilevel approach is directly applicable to perform design tradeoffs for running part of the algorithms near the instruments and the rest on remote (analysis) systems.

97 MATHEMATICS AND COMPUTING↗