Search NASA⌕ Search

SEARCH · Search NASA

Results for “data intensive computing”

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,657 records · Page 92

Fast equilibration of coarse-grained polymeric liquids

The study of macromolecular systems may require large computer simulations that are too time consuming and resource intensive to execute in full atomic detail. The integral equation coarse-graining approach by Guenza and co-workers enables the exploration of longer time and spatial scales without sacrificing thermodynamic consistency, by approximating collections of atoms using analytically-derived soft-sphere potentials. Because coarse-grained (CG) characterizations evolve polymer systems far more efficiently than the corresponding united atom (UA) descriptions, we can feasibly equilibrate a CG system to a reasonable geometry, then transform back to the UA description for a more complete equilibration. Automating the transformation between the two different representations simultaneously exploits CG efficiency and UA accuracy. By iteratively mapping back and forth between CG and UA, we can quickly guide the simulation towards a configuration that would have taken many more time steps within the UA representation alone. Accomplishing this feat requires a diligent workflow for managing input/output coordinate data between the different steps, deriving the potential at runtime, and inspecting convergence. Here in this article, we present a lightweight workflow environment that accomplishes such fast equilibration without user intervention. The workflow supports automated mapping between the CG and UA descriptions in an iterative, scalable, and customizable manner. We describe this technique, examine its feasibility, and analyze its correctness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detecting and Removing Data Artifacts in Hadamard Transform Ion Mobility-Mass Spectrometry Measurements

Applying Hadamard transform multiplexing to ion mobility separations (IMS) can significantly improve the signal-to-noise ratio and throughput for IMS coupled mass spectrometry (MS) measurements by increasing the ion utilization efficiency. However, it has been determined that both fluctuations in ion intensity as well as spatial shifts in the multiplexed data lower the signal-to-noise ratios and appear as noise in downstream processing of the data. To address this problem, we have developed a novel algorithm that discovers and eliminates data artifacts. The algorithm uses knowledge of the true signal peaks derived from the encoded data and allows for both artifacts and noise to be removed with high confidence, decreasing the likelihood of false identifications in subsequent data processing. The result is that IMS-MS can be applied to increase measurement sensitivity while avoiding artifacts that have previously limited its utility.

47 OTHER INSTRUMENTATION↗

Characterization of an imploding cylindrical plasma for electron transport studies using x-ray emission spectroscopy

We report here on the characterization of the conditions of an imploding cylindrical plasma by time-resolved x-ray emission spectroscopy. Knowledge about this implosion platform can be applied to studies of particle transport for inertial confinement fusion schemes or to astrophysical plasmas. A cylindrical Cl-doped CH foam within a tube of solid CH was irradiated by 36 beams (I total ~5×10 14 W/cm 2 , 1.5ns square pulse, and E total ~ 16.2kJ) of the OMEGA-60 laser to radially compress the CH toward the axis. The analysis of the time-resolved spectra showed that the compression can be described by four distinct phases, each presenting different plasma conditions. First the ablation of the cylinder is dominant; second, the foam is heated and induces a significant jump in emission intensities; third, the temperature and density of the foam reaches a maximum; and finally, the plasma expands. Ranges for the plasma temperature were inferred with the atomic physics code SCRAM (Spectroscopic Collisional-Radiative Atomic Model) and the experimental data have been compared to hydrodynamic simulations performed with the 2D code FLASH, which showed a similar implosion dynamic over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Implementing High-Performance Geometric Multigrid Solver with Naturally Grained Messages

Structured-grid linear solvers often require manually packing and unpacking of communication data to achieve high performance.Orchestrating this process efficiently is challenging, labor-intensive, and potentially error-prone.In this paper, we explore an alternative approach that communicates the data with naturally grained messagesizes without manual packing and unpacking. This approach is the distributed analogue of shared-memory programming, taking advantage of the global addressspace in PGAS languages to provide substantial programming ease. However, its performance may suffer from the large number of small messages. We investigate theruntime support required in the UPC ++ library for this naturally grained version to close the performance gap between the two approaches and attain comparable performance at scale using the High-Performance Geometric Multgrid (HPGMG-FV) benchmark as a driver.

97 MATHEMATICS AND COMPUTING↗

D-Guaranteed Discrimination of Statistical Hypotheses: a Review of Results and Unsolved Problems

We compare two sequential d-guaranteed tests and an optimal d-guaranteed test based on a fixed number of observations with respect to the average number of observations within the most accepted practical applications of Bayesian probabilistic models. We consider the sequential “first skipping” test, the sequential locally efficient test based on the score statistic, and the test based on a fixed number of observations that minimizes the necessary sample size. We study various characteristics of these tests connected with the number of observations within three probabilistic models, namely, the normal (ϑ, σ{sup 2}) distribution of the observed random variable and the normal a priori distribution of ϑ with fixed σ{sup 2}; the exponential distribution with the intensity parameter ϑ and the a priori gamma distribution of ϑ; and the Bernoulli sampling with the success probability ϑ and the a priori beta distribution of ϑ. We discuss the connection of the d-posterior approach with the compound decision problem as applied to the analysis of data provided by microchips (when the false discovery rate, or FDR, for short, is treated as the d-risk of the first kind). We present the vast data on characteristics of the mentioned tests obtained by the method of mathematical modeling in several tables. We discuss unsolved problems of the d-guaranteed discrimination of hypotheses with the minimal number of observations and approaches to their solution.

97 MATHEMATICS AND COMPUTING↗

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗

Porting the WAVEWATCH III (v6.07) wave action source terms to GPU

Abstract. Surface gravity waves play a critical role in several processes, including mixing, coastal inundation, and surface fluxes. Despite the growing literature on the importance of ocean surface waves, wind–wave processes have traditionally been excluded from Earth system models (ESMs) due to the high computational costs of running spectral wave models. The development of the Next Generation Ocean Model for the DOE’s (Department of Energy) E3SM (Energy Exascale Earth System Model) Project partly focuses on the inclusion of a wave model, WAVEWATCH III (WW3), into E3SM. WW3, which was originally developed for operational wave forecasting, needs to be computationally less expensive before it can be integrated into ESMs. To accomplish this, we take advantage of heterogeneous architectures at DOE leadership computing facilities and the increasing computing power of general-purpose graphics processing units (GPUs). This paper identifies the wave action source terms, W3SRCEMD, as the most computationally intensive module in WW3 and then accelerates them via GPU. Our experiments on two computing platforms, Kodiak (P100 GPU and Intel(R) Xeon(R) central processing unit, CPU, E5-2695 v4) and Summit (V100 GPU and IBM POWER9 CPU) show respective average speedups of 2× and 4× when mapping one Message Passing Interface (MPI) per GPU. An average speedup of 1.4× was achieved using all 42 CPU cores and 6 GPUs on a Summit node (with 7 MPI ranks per GPU). However, the GPU speedup over the 42 CPU cores remains relatively unchanged (∼ 1.3×) even when using 4 MPI ranks per GPU (24 ranks in total) and 3 MPI ranks per GPU (18 ranks in total). This corresponds to a 35 %–40 % decrease in both simulation time and usage of resources. Due to too many local scalars and arrays in the W3SRCEMD subroutine and the huge WW3 memory requirement, GPU performance is currently limited by the data transfer bandwidth between the CPU and the GPU. Ideally, OpenACC routine directives could be used to further improve performance. However, W3SRCEMD would require significant code refactoring to make this possible. We also discuss how the trade-off between the occupancy, register, and latency affects the GPU performance of WW3.

58 GEOSCIENCES↗

A stochastic precipitation disaggregation scheme for GCM applications

In the surface hydrologic parameterization of general circulation models (GCMs), it is commonly assumed that the precipitation processes are homogeneous over a GCM grid square and that the precipitation intensity is uniformly distributed. Based on evidence that the spatial distribution of precipitation within a GCM grid square is crucial for the land surface hydrology parameterization, a few researchers have explored the impacts of assuming that the precipitation is exponentially distributed. This paper explores the suitability of the aforementioned assumptions. First, a statistical analysis is conducted of historical precipitation data for three GCM grids in different regions of the United States. The analysis suggests that neither the uniform nor the exponential distribution assumption may be suitable at the GCM grid scale and, that instead, the spatial variability in precipitation is characterized by statistical patterns that are inhomogeneous. These patterns vary from grid to grid and are induced by the interaction between atmospheric conditions and various land surface characteristics, such as topographical features, surface properties, etc. Within the same grid square, however, the statistical patterns are generally constant from year to year. Based on this analysis, a computationally viable (i.e., usable with GCMs) stochastic precipitation disaggregation scheme that utilizes these stable statistical patterns is proposed. The method was used to generate spatially distributed hourly rainfall for a summer season in the southwestern region of the continental United States. Analysis of the results shows that the methodology preserves the seasonal characteristics of spatial variability in precipitation that is observed in the long-term historical data.

Gao, Xiaogang↗

Textural analysis on contrast-enhanced CT in pancreatic neuroendocrine neoplasms: association with WHO grade

Purpose: Grades of pancreatic neuroendocrine neoplasms (PNENs) are associated with the choice of treatment strategies. Texture analysis has been used in tumor diagnosis and staging evaluation. In this study, we aim to evaluate the potential ability of texture parameters in differentiation of PNENs grades. Materials and methods: 37 patients with histologically proven PNENs and underwent pretreatment dynamic contrast-enhanced computed tomography examinations were retrospectively analyzed. Imaging features and texture features at contrast-enhanced images were evaluated. Receiver operating characteristic curves were used to determine the cut-off values and the sensitivity and specificity of prediction. Results: There were significant differences in tumor margin, pancreatic duct dilatation, lymph nodes invasion, size, portal enhancement ratio (PER), arterial enhancement ratio (AER), mean grey-level intensity, kurtosis, entropy, and uniformity among G1, G2, and pancreatic neuroendocrine carcinoma (PNEC) G3 (p < 0.01). Similar results were found between pancreatic neuroendocrine tumors (PNETs) G1/G2 and PNEC G3. AER and PER showed the best sensitivity (0.86–0.94) and specificity (0.92–1.0) for differentiating PNEC G3 from PNETs G1/G2. Mean grey-level intensity, entropy, and uniformity also showed acceptable sensitivity (0.73–0.91) and specificity (0.85–1.0). Mean grey-level intensity was also showed acceptable sensitivity (91% to 100%) and specificity (82% to 91%) in differentiating PNET G1 from PNET G2. Conclusions: Our data indicated that texture parameters have potential in grading PNENs, in particular in differentiating PNEC G3 from PNETs G1/G2.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Desalination metamodels and a framework for cross-comparative performance simulations

There is an opportunity to save energy and reduce operational expenses when choosing a suitable desalination method aided by computational modeling. Existing models are not conducive to generalized comparisons between different desalination methods. Therefore, the work in this study developed metamodels for six desalination methods, grouped them into thermal and molecular transport families, and validated their predictive performance within 9% difference from published data. This validated framework allowed comparisons of desalination methods at their prescribed ranges of operational conditions that they were designed for. These conditions specify feed salinity ranges of 1.6 to 2.4 g/kg for Capacitive Deionization and Reverse Osmosis (RO), 2.8 to 4.2 g/kg for Electrodialysis, 28 to 42 g/kg for Thermovapor Compression and Humidification-Dehumidification, and 37 to 55 g/kg for Multi-Effect Distillation (MED). Despite different operational conditions, all models exhibit non-linear, positive correlation between energy consumption and system size in response to feed salinity and production rate. The framework is also employed in a cross-comparative analysis between MED and RO whose results suggest that energy intensity for MED is an order of magnitude greater than RO for the same operational conditions, but actual operational costs are comparable. Overall, the framework is ready for deployment in case studies of actual desalination plants.

42 ENGINEERING↗

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗

Convective flow effects on protein crystal growth

A high-resolution microscopic interferometric setup for the monitoring of protein morphologies has been developed. Growth or dissolution of a crystal can be resolved with a long-term depth resolution of 200 A and a lateral resolution of 2 microns. This capability of simultaneously monitoring the interfacial displacement with high local depth resolution has yielded several novel results. We have found with lysozyme that (1) the normal growth rate is oscillatory, and (2) depending on the impurity content of the solution, the growth step density is either greater or lower at the periphery of a facet than in its center. The repartitioning of Na plus and Cl minus ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed, to interpret the large body of data in unified way. The results strongly suggest that (1) the ion to lysozyne ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter electron microscopy results appear to confirm this finding, which could have far-reaching consequences for x-ray diffraction studies. A computational model for diffusive-convective transport in protein crystallization has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies and our kinetics data for the growth of lysozyme. The results show that even in the small cell employed, protein concentration nonuniformities and gravity-driven solutal convection can be significant. The calculated convection velocities are of the same order to magnitude as those found in earlier experiments. As expected, convective transport, i.e., at Og, lysozyme crystal growth remains kinetically limited. The salt distribution in the crystal is predicted to be non-uniform at both 1g and 0g, as a consequence of protein depletion in the solution. Static and dynamic light scattering studies in undersaturated and supersaturated solutions have been performed. Diffusivities in undersaturated solutions, were found to vary with lysozyme concentrations. Depending on the salt concentration, the diffusivities either increase or decrease. Interestingly, the corresponding static scattering intensities behave oppositely, Our current analysis indicates that these changes are inconsistent with aggregation in undersaturated solutions. However, the data are compatible with concentration-dependent changes of the interactions between protein and salt.

Rosenberger, Franz↗

A Statistical Study of the Global Structure of the Ring Current

In this paper we derive the average configuration of the ring current as a function of the state of the magnetosphere as indicated by the Dst index. We sort magnetic field data from the Combined Release and Radiation Effects Satellite (CRRES) by spatial location and by the Dst index in order to produce magnetic field maps. From these maps we calculate local current systems by taking the curl of the magnetic field. NN7e find both the westward (outer) and the eastward (inner) components of the ring current. We find that the ring current intensity varies linearly with D.St as expected, and that the ring current is asymmetric for all Dst values. The azimuthal peak of the ring current is located in the afternoon sector for quiet conditions, and near midnight for disturbed conditions. The ring current also moves closer to the Earth during disturbed conditions. We are able to recreate the Dst index by integrating the magnetic perturbations caused by the ring current. We find that we needed to apply a 20 nT offset to Dst, and assume a perfectly conducting Earth to obtain an optimal agreement between the computed and the observed Dst. We interpret the 20 nT offset as the magnetic field generated by the quiet time ring current used as baseline in computing Dst.

Jorgensen, A. M.↗

Particulate Emission and Optical Measurements in n-Heptane Low-Swirl Non-Premixed Flames at Elevated Pressures

A series of ground combustion tests conducted at elevated pressures demonstrates in-chamber optical measurements of soot emissions, combined with conventional gaseous and particle emissions sampling techniques. Utilizing a micro-radial-entry counter-swirl (MRX) burner, we successfully stabilized an n-heptane low-swirl, non-premixed spray flame across a wide range of equivalence ratios at pressures reaching 250 psia. Our optical diagnostics feature a temporally gated diffused-backlight illumination extinction imaging (DBI-EI) system, equipped with a 470-nm strobe high-intensity LED array as the illumination source, along with a line-of-sight telecentric camera. By accessing the flame tube, the pulsed DBI-EI provides illuminating optical extinction images, enabling experimental estimation of soot volume fractions within the reaction zone. This imaging technique will allow for real-time observation of soot formation dynamics occurring above the air-blast fuel injector. To further enhance our analysis, we employed simultaneous gaseous and particulate extractive diagnostics, interpreting in-flame particulate information through CO2 sampling, exhaust particle sizing, and particle counting. The processed dataset reveals a promising application of DBI-EI for propulsion emission analysis in high-pressure environments, where optical beam-steering effects can significantly impact results. This report is essential to bridging the divide between in-flame optical diagnostics and exhaust extractive measurements. By achieving an understanding of non-premixed combustion emission characteristics, we attempt to provide invaluable data for computer code validation, advancing the field of combustion research and supporting the development of cleaner, more efficient engines.

high-pressure flame↗

Particulate Emission and Optical Measurements in n-Heptane Low-Swirl Non-Premixed Flames at Elevated Pressures

A series of ground combustion tests conducted at elevated pressures demonstrates in-chamber optical measurements of soot emissions, combined with conventional gaseous and particle emissions sampling techniques. Utilizing a micro-radial-entry counter-swirl (MRX) burner, we successfully stabilized an n-heptane low-swirl, non-premixed spray flame across a wide range of equivalence ratios at pressures reaching 250 psia. Our optical diagnostics feature a temporally gated diffused-backlight illumination extinction imaging (DBI-EI) system, equipped with a 470-nm strobe high-intensity LED array as the illumination source, along with a line-of-sight telecentric camera. By accessing the flame tube, the pulsed DBI-EI provides illuminating optical extinction images, enabling experimental estimation of soot volume fractions within the reaction zone. This imaging technique will allow for real-time observation of soot formation dynamics occurring above the air-blast fuel injector. To further enhance our analysis, we employed simultaneous gaseous and particulate extractive diagnostics, interpreting in-flame particulate information through CO2 sampling, exhaust particle sizing, and particle counting. The processed dataset reveals a promising application of DBI-EI for propulsion emission analysis in high-pressure environments, where optical beam-steering effects can significantly impact results. This report is essential to bridging the divide between in-flame optical diagnostics and exhaust extractive measurements. By achieving an understanding of non-premixed combustion emission characteristics, we attempt to provide invaluable data for computer code validation, advancing the field of combustion research and supporting the development of cleaner, more efficient engines.

high-pressure flame↗

Resilience Metrics Framework for Solar Photovoltaics

This presentation was given at the Photovoltaic Specialist Conference (PVSC) 54 in New Orleans, Louisiana. Photovoltaic (PV) systems are routinely exposed to extreme weather, including wind and hail storms. Historically, most systems have proven to be resilient to such events, but some storms have damaged PV systems, leading to physical and financial loss. Storm hardening measures and specific system attributes can reduce this risk. This work introduces a set of resilience metrics and a framework for quantifying, comparing, and predicting PV system resilience. The framework is divided into two parts: 1) predictive, attribute metrics based on site and component characteristics, and 2) impact metrics that assess post-storm performance. Metrics are weighted and aggregated, producing hazard-specific resilience scores. We derive damage functions from storm-impacted PV systems, establishing a baseline against which post-storm performance can be compared. This damage was widely variable across hail and wind intensities, and field hail damage was less than predicted by laboratory tests, suggesting that system features - in addition to storm conditions - influence damage likelihood. Finally, the metrics framework is demonstrated using three case studies of storm damaged PV systems. Although additional data are needed to create attribute specific damage functions and establish metric weights, this study presents a methodology for evaluating PV resilience and contributes new damage functions to the literature.

14 SOLAR ENERGY↗

Benchmark Generation and Simulation at Extreme Scale

The path to extreme scale high-performance computing (HPC) poses several challenges related to power, performance, resilience, productivity, programmability, data movement, and data management. Investigating the performance of parallel applications at scale on future architectures and the performance impact of different architectural choices is an important component of HPC hardware/software co-design. Simulations using models of future HPC systems and communication traces from applications running on existing HPC systems can offer an insight into the performance of future architectures. This work targets technology developed for scalable application tracing of communication events. It focuses on extreme-scale simulation of HPC applications and their communication behavior via lightweight parallel discrete event simulation for performance estimation and evaluation. Instead of simply replaying a trace within a simulator, this work promotes the generation of a benchmark from traces. This benchmark is subsequently exposed to simulation using models to reflect the performance characteristics of future-generation HPC systems. This technique provides a number of benefits, such as eliminating the data intensive trace replay and enabling simulations at different scales. The presented work features novel software co-design aspects, combining the ScalaTrace tool to generate scalable trace files, the ScalaBenchGen tool to generate the benchmark, and the xSim tool to assess the benchmark characteristics within a simulator.

Lagadapati, Mahesh↗

A method for normalizing pathology images to improve feature extraction for quantitative pathology

Purpose: With the advent of digital slide scanning technologies and the potential proliferation of large repositories of digital pathology images, many research studies can leverage these data for biomedical discovery and to develop clinical applications. However, quantitative analysis of digital pathology images is impeded by batch effects generated by varied staining protocols and staining conditions of pathological slides. Methods: To overcome this problem, this paper proposes a novel, fully automated stain normalization method to reduce batch effects and thus aid research in digital pathology applications. Their method, intensity centering and histogram equalization (ICHE), normalizes a diverse set of pathology images by first scaling the centroids of the intensity histograms to a common point and then applying a modified version of contrast-limited adaptive histogram equalization. Normalization was performed on two datasets of digitized hematoxylin and eosin (H&E) slides of different tissue slices from the same lung tumor, and one immunohistochemistry dataset of digitized slides created by restaining one of the H&E datasets. Results: The ICHE method was evaluated based on image intensity values, quantitative features, and the effect on downstream applications, such as a computer aided diagnosis. For comparison, three methods from the literature were reimplemented and evaluated using the same criteria. The authors found that ICHE not only improved performance compared with un-normalized images, but in most cases showed improvement compared with previous methods for correcting batch effects in the literature. Conclusions: ICHE may be a useful preprocessing step a digital pathology image processing pipeline.

61 RADIATION PROTECTION AND DOSIMETRY↗