Search NASA⌕ Search

SEARCH · Search NASA

Results for “cluster counts”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Parallel performance of algebraic multigrid domain decomposition

Algebraic multigrid (AMG) is a widely used scalable solver and preconditioner for large-scale linear systems resulting from the discretization of a wide class of elliptic PDEs. While AMG has optimal computational complexity, the cost of communication has become a significant bottleneck that limits its scalability as processor counts continue to grow on modern machines. This article examines the design, implementation, and parallel performance of a novel algorithm, algebraic multigrid domain decomposition (AMG-DD), designed specifically to limit communication. The goal of AMG-DD is to provide a low-communication alternative to standard AMG V-cycles by trading some additional computational overhead for a significant reduction in communication cost. Numerical results show that AMG-DD achieves superior accuracy per communication cost compared with AMG, and speedup over AMG is demonstrated on a large GPU cluster.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Towards Superior Software Portability with SHAD and HPX C++ Libraries

As hardware architectures and software stacks complexity grows, development productivity, performance and software portability, quickly evolve from desirable features to actual needs. SHAD, the Scalable High-performance Algorithms and Data-structures C++ library is designed to mitigate these issues: it provides general purpose building blocks as well as high-level custom utilities, and offers a shared-memory programming abstraction which facilitates the programming of complex systems, scaling up to High Performance Computing clusters. SHAD’s portability is achieved through an abstract runtime interface, which decouples the upper layers of the library and hides the low level details of the underlying architecture. This layer enables SHAD to interface with different runtime/threading systems, e.g. Intel TBB and Global Memory and Threading (GMT). However, current backends targeting distributed systems, rely on a centralized controller which may possibly limit scalability up to hundreds of nodes and creates a network hot spot due to all to one communication for synchronization, and possibly resulting in degraded performance at high process counts. In this research, we explore HPX, the C++ standard library for parallelism and concurrency, as an additional backend in support of the SHAD library, and present the methodologies in support of local and remote task executions in SHAD with respect to HPX. Finally, we evaluate the proposed system by comparing against existing backends of SHAD and analyzing their performance on C++ Standard Template Library algorithms.

Wu, Nanmiao↗
NASA NTRS2004

Detecting the Sun in 2-Gyr Old Open Clusters

The aim of this observation was to take advantage of the large collecting area and sensitivity of the European Space Agency's XMM-Newton telescope to survey the X ray brightness of solar-mass stars in the 2 Gyr old NGC 752 cluster in order to chart out the dependence of coronal X ray emission of Sun like stars as a function of age. Dr. Simon's stated role as one of the co-investigators in the project was limited to assisting the project's European Principal Investigator with the analysis and interpretation of the X ray measurements, and to assist in summarizing the results for publication. The observation was executed successfully, and approximately 130 X ray sources were detected in the XMM field of view and subsequently measured to obtain X ray count rates. To date neither the measurements nor the raw data have been distributed by the European P.I. to the co-investigators, It is expected that the measurements will be distributed to the co-investigators by the end of 2004, and once that is done, the scientific analysis by Dr. Simon can then begin. At present, since the data and measurements have not been shared with the co-investigators, there are no results to report.

Simon, Theodore↗
DOE OSTI2021

Comparative genome analyses suggest a hemibiotrophic lifestyle and virulence differences for the beech bark disease fungal pathogens Neonectria faginata and Neonectria coccinea

Abstract Neonectria faginata and Neonectria coccinea are the causal agents of the insect-fungus disease complex known as beech bark disease (BBD), known to cause mortality in beech forest stands in North America and Europe. These fungal species have been the focus of extensive ecological and disease management studies, yet less progress has been made toward generating genomic resources for both micro- and macro-evolutionary studies. Here, we report a 42.1 and 42.7 mb highly contiguous genome assemblies of N. faginata and N. coccinea, respectively, obtained using Illumina technology. These species share similar gene number counts (12,941 and 12,991) and percentages of predicted genes with assigned functional categories (64 and 65%). Approximately 32% of the predicted proteomes of both species are homologous to proteins involved in pathogenicity, yet N. coccinea shows a higher number of predicted mitogen-activated protein kinase genes, virulence determinants possibly contributing to differences in disease severity between N. faginata and N. coccinea. A wide range of genes encoding for carbohydrate-active enzymes capable of degradation of complex plant polysaccharides and a small number of predicted secretory effector proteins, secondary metabolite biosynthesis clusters and cytochrome oxidase P450 genes were also found. This arsenal of enzymes and effectors correlates with, and reflects, the hemibiotrophic lifestyle of these two fungal pathogens. Phylogenomic analysis and timetree estimations indicated that the N. faginata and N. coccinea species divergence may have occurred at ∼4.1 million years ago. Differences were also observed in the annotated mitochondrial genomes as they were found to be 81.7 kb (N. faginata) and 43.2 kb (N. coccinea) in size. The mitochondrial DNA expansion observed in N. faginata is attributed to the invasion of introns into diverse intra- and intergenic locations. These first draft genomes of N. faginata and N. coccinea serve as valuable tools to increase our understanding of basic genetics, evolutionary mechanisms and molecular physiology of these two nectriaceous plant pathogenic species.

Salgado-Salazar, Catalina↗
DOE OSTI2020

Fisher for complements: extracting cosmology and neutrino mass from the counts-in-cells PDF

In this work, we comprehensively analyse the cosmology dependence of counts-in-cells statistics. Furthermore, we focus on the shape of the one-point probability distribution function (PDF) of the matter density field at mildly non-linear scales. Based on large-deviation statistics, we parametrize the cosmology dependence of the matter PDF in terms of the linear power spectrum, the growth factor, the spherical collapse dynamics, and the non-linear variance. We extend our formalism to include massive neutrinos, finding that the total matter PDF is highly sensitive to the total neutrino mass Mν and can disentangle it from the clustering amplitude σ 8 . Additionally, using more than a million PDFs extracted from the Quijote simulations, we determine the response of the matter PDF to changing parameters in the νΛCDM model and successfully cross-validate the theoretical model and the simulation measurements. We present the first νΛCDM Fisher forecast for the matter PDF at multiple scales and redshifts, and its combination with the matter power spectrum. We establish that the matter PDF and the matter power spectrum are highly complementary at mildly non-linear scales. The matter PDF is particularly powerful for constraining the matter density Ωm, clustering amplitude σ 8 and the total neutrino mass Mν. Adding the mildly non-linear matter PDF to the mildly non-linear matter power spectrum improves constraints on Ω m by a factor of 5 and σ 8 by a factor of 2 when considering the three lowest redshifts. In our joint analysis of the matter PDF and matter power spectrum at three redshifts, the total neutrino mass is constrained to better than 0.01 eV with a total volume of 6 (Gpc h -1 ) 3 . We discuss how density-split statistics can be used to translate those encouraging results for the matter PDF into realistic observables in galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Implications of crater distributions on Venus

The horizontal locations of craters on Venus are consistent with randomness. However, (1) randomness does not make crater counts useless for age indications; (2) consistency does not imply necessity or optimality; and (3) horizontal location is not the only reference frame against which to test models. Re (1), the apparent smallness of resurfacing areas means that a region on the order of one percent of the planet with a typical number of craters, 5-15, will have a range of feature ages of several 100 My. Re (2), models of resurfacing somewhat similar to Earth's can be found that are also consistent and more optimal than random: i.e., resurfacing occurring in clusters, that arise and die away in lime intervals on the order of 50 My. These agree with the observation that there are more areas of high crater density, and fewer of moderate density, than optimal for random. Re (3), 799 crater elevations were tested; there are more at low elevations and fewer at high elevations than optimal for random: i.e., 54.6 percent below the median. Only one of 40 random sets of 799 was as extreme.

Kaula, W. M.↗
NASA NTRS2003

Coronal Activity in the R CrA T Association

For this project we intend to extract net counts and variability information for the X-ray sources found in 9 archival ROSAT PSPC and 6 archival ROSAT HRI images in the region of the R CrA T association. These data will be merged with an extensive database of optical photometry and spectroscopy plus published near-infrared photometry to (1) identify new association members and (2) derive, for the first time, an accurate coronal luminosity function for the T Tauri members of this T association. These data will be used to make direct comparisons between the coronal luminosity functions for other T associations and those of large clusters to assess (a) the importance of the star-formation environment in initial coronal activity levels, (b) the effects of PMS evolution on dynamo activity as a function of mass and age, and (c) the level of contamination by field post-T Tauri stars on association membership surveys.

Patten, Brian M.↗
DOE OSTI2023

Comparing the Performance of Julia on CPUs versus GPUs and Julia-MPI versus Fortran-MPI: a case study with MPAS-Ocean (Version 7.1)

Abstract. Some programming languages are easy to develop at the cost of slow execution, while others are fast at runtime but much more difficult to write. Julia is a programming language that aims to be the best of both worlds – a development and production language at the same time. To test Julia's utility in scientific high-performance computing (HPC), we built an unstructured-mesh shallow water model in Julia and compared it against an established Fortran-MPI ocean model, the Model for Prediction Across Scales–Ocean (MPAS-Ocean), as well as a Python shallow water code. Three versions of the Julia shallow water code were created: for single-core CPU, graphics processing unit (GPU), and Message Passing Interface (MPI) CPU clusters. Comparing identical simulations revealed that our first version of the Julia model was 13 times faster than Python using NumPy, where both used an unthreaded single-core CPU. Further Julia optimizations, including static typing and removing implicit memory allocations, provided an additional 10–20× speed-up of the single-core CPU Julia model. The GPU-accelerated Julia code was almost identical in terms of performance to the MPI parallelized code on 64 processes, an unexpected result for such different architectures. Parallelized Julia-MPI performance was identical to Fortran-MPI MPAS-Ocean for low processor counts and ranges from 2× faster to 2× slower for higher processor counts. Our experience is that Julia development is fast and convenient for prototyping but that Julia requires further investment and expertise to be competitive with compiled codes. We provide advice on Julia code optimization for HPC systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

The Analysis of the Patterns of Radiation-Induced DNA Damage Foci by a Stochastic Monte Carlo Model of DNA Double Strand Breaks Induction by Heavy Ions and Image Segmentation Software

To create a generalized mechanistic model of DNA damage in human cells that will generate analytical and image data corresponding to experimentally observed DNA damage foci and will help to improve the experimental foci yields by simulating spatial foci patterns and resolving problems with quantitative image analysis. Material and Methods: The analysis of patterns of RIFs (radiation-induced foci) produced by low- and high-LET (linear energy transfer) radiation was conducted by using a Monte Carlo model that combines the heavy ion track structure with characteristics of the human genome on the level of chromosomes. The foci patterns were also simulated in the maximum projection plane for flat nuclei. Some data analysis was done with the help of image segmentation software that identifies individual classes of RIFs and colocolized RIFs, which is of importance to some experimental assays that assign DNA damage a dual phosphorescent signal. Results: The model predicts the spatial and genomic distributions of DNA DSBs (double strand breaks) and associated RIFs in a human cell nucleus for a particular dose of either low- or high-LET radiation. We used the model to do analyses for different irradiation scenarios. In the beam-parallel-to-the-disk-of-a-flattened-nucleus scenario we found that the foci appeared to be merged due to their high density, while, in the perpendicular-beam scenario, the foci appeared as one bright spot per hit. The statistics and spatial distribution of regions of densely arranged foci, termed DNA foci chains, were predicted numerically using this model. Another analysis was done to evaluate the number of ion hits per nucleus, which were visible from streaks of closely located foci. In another analysis, our image segmentaiton software determined foci yields directly from images with single-class or colocolized foci. Conclusions: We showed that DSB clustering needs to be taken into account to determine the true DNA damage foci yield, which helps to determine the DSB yield. Using the model analysis, a researcher can refine the DSB yield per nucleus per particle. We showed that purely geometric artifacts, present in the experimental images, can be analytically resolved with the model, and that the quantization of track hits and DSB yields can be provided to the experimentalists who use enumeration of radiation-induced foci in immunofluorescence experiments using proteins that detect DNA damage. An automated image segmentaiton software can prove useful in a faster and more precise object counting for colocolized foci images.

Ponomarev, Artem↗
DOE OSTI2020

Development of Gamma Background Radiation Digital Twin with Machine Learning Algorithms: Application of Unsupervised Machine Learning to Detection of Anomalies and Nuisances in Gamma Background Radiation Environmental Screening Data

Environmental screening of gamma radiation consists of detecting weak nuisance and anomaly signal in the presence of strong and highly varying background. In a typical scenario, a mobile detector-spectrometer continuously measures gamma radiation spectra in short, e.g., one-second, signal acquisition intervals. The measurement data is a 2D matrix, where one dimension is gamma ray energy, and the other dimension is the number of measurements or total time. In principle, gamma radiation sources can be detected and identified from the measured data by their unique spectral lines. Detecting sources from data measured in a search scenario is difficult due to the highly varying background because of naturally occurring radioactive material (NORM), and low signal-to-noise ratio (S/N) of spectral signal measured during one-second acquisition intervals. The objective of this work is to explore unsupervised machine learning (ML) algorithms for development of a digital twin of gamma radiation background, and for detection and identification of weak nuisances and anomalies events in the presence of highly fluctuating background. In one segment of work, we developed a gamma background estimation model using a Longshort term memory (LSTM) network for one-step CPS time series prediction. The LSTM model was validated with two data sets of measurements from two independent NaI detectors positioned on a mobile platform. The data sets contained background radiation only and no orphan isotope sources. The LSTM model was constructed and tested using data from one of the detectors. Performance of the LSTM model was validate through one-step prediction of CPS time series of another NaI detector without re-training. This approach allows to create a digital twin for nuclear background estimation. Using LSTM, it could be possible to detect a source through subtraction of the estimated counts from the measured background. In another segment of work, we investigated detection of gamma emitting sources in the presence of complex background using unsupervised machine learning. Spectral lines of isotopes are difficult to observe in one-second measurements. Averaging over the entire measurement campaign data set reveals spectral lines of most common background isotopes. Spectral lines of orphan sources, which might appear only in a few measurements during the campaign, will be washed out if averaging is performed over the entire measurement data set. The approach we have explored consists of extracting one-second measurements containing weak spectral features through data clustering. Averaging one-second spectra in a cluster should reveal the presence of anomaly sources. We created two ML models using K-means clustering and Neural Network Self-organizing Map (SOM). Performance of these ML models was benchmarked using search data. One data set contained 137 Cs source, and another dataset contained 131 I source.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

The stellar population and luminosity function in M31 bulge and Inner Disk Fields

We report infrared photometry and stellar identifications for stars in five fields in the M31 bulge located from 2 to 11 arcmin from the nucleus. These fields have been chosen such that the bulge/disk star ratio predicted from Kent's (1989) small bulge model varies from 7:1 to 1:5, allowing a study of near pure disk and near pure bulge stellar populations. We reject the hypothesis of Davies et al. (1991) that luminous stars found within 500 pc of the nucleus are due to a contaminating disk population. We find that the bulge contains stars in excess of M(sub bol) = -5 mag and that the bulge luminosity function has a distinct shape different from the disk fields. We find many stars redder than (J-K) = 2 mag, and suggest that these stars may be the counterparts of the IRAS-selected Galactic bulge Miras studied by Whitelock et at. (1991). The number of bright stars (M(sub bol) is less than -5 mag) falls off more rapidly than the r band surface brightness. By building model fields out of a bulge luminosity function and artificial stars, we are able to show that the change in the luminosity function toward the center cannot be explained simply by the mismeasurement of overcrowded star images. However, these tests also raise the possibility that the asymptotic giant branch (AGB) tip may be approximately equal to 1 mag fainter than actually measured in our most crowded field, reaching only M(sub bol) = -5. We compare observed counts of AGB stars with those predicted from theoretical lifetimes using a technique of general interest for this problem, the Fuel Consumption Theorem of Renzini & Buzzoni (1986) Spectral Evolution of Galaxies (Reidel, Dordrecht). Our methodology is generally applicable to the study of other resolved extragalactic stellar populations. The number of observed stars per magnitude up to a luminosity of M(bol) = -5.5 mag is consistent with AGB evolution of the whole population of the innermost bulge field with the standard lifetime on the AGB of 1.3 Myr/mag. We advance the possibility that the bulge of M31 may be younger than the oldest Galactic globular clusters. We note that M33 has recently been found to have an r(exp 1/4)-law spheroid consisting of intermediate-age stars; bulges can form later than the old halo population.

Rich, R. Michael↗
DOE OSTI2022

Chromatin structure predicts survival in glioma patients

The pathological changes in epigenetics and gene regulation that accompany the progression of low-grade to high-grade gliomas are under-studied. The authors use a large set of paired atac-seq and RNA-seq data from surgically resected glioma specimens to infer gene regulatory relationships in glioma. Thirty-eight glioma patient samples underwent atac-seq sequencing and 16 samples underwent additional RNA-seq analysis. Using an atac-seq/RNA-seq correlation matrix, atac-seq peaks were paired with genes based on high correlation values (|r 2 | > 0.6). Samples clustered by IDH1 status but not by grade. Surprisingly there was a trend for IDH1 mutant samples to have more peaks. The majority of peaks are positively correlated with survival and positively correlated with gene expression. Constructing a model of the top six atac-seq peaks created a highly accurate survival prediction model (r 2 = 0.68). Four of these peaks were still significant after controlling for age, grade, pathology, IDH1 status and gender. Grade II, III, and IV (primary) samples have similar transcription factors and gene modules. However, grade IV (recurrent) samples have strikingly few peaks. Patient-derived glioma cultures showed decreased peak counts following radiation indicating that this may be radiation-induced. This study supports the notion that IDH1 mutant and IDH1 wildtype gliomas have different epigenetic landscapes and that accessible chromatin sites mapped by atac-seq peaks tend to be positively correlated with expression. The data in this study leads to a new model of treatment response wherein glioma cells respond to radiation therapy by closing open regions of DNA.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2005

Coronal Activity in the R CrA T Association

Brian Patten is the Principal Investigator of the NASA ROSS-ADP project Coronal Activity in the R CrA T Association. For this project we have extracted net counts and variability information for all of the X-ray sources found in 23 archival ROSAT PSPC and HRI images in the region of the R CrA T association. These data have been merged with an extensive database of optical and near-infrared photometry, optical spectroscopy, and parallax data. These data have been used to (1) identify new association members and clarify the membership status of a number of previously suspected members of the association, and (2) derive, for the first time, an accurate coronal luminosity function for the T Tauri members of this T association and make direct comparisons between the coronal luminosity functions for other T associations and those of large clusters. We have used our survey data to assess (a) the importance of the star-formation environment in initial coronal activity levels, (b) the effects of PMS evolution on dynamo activity as a function of mass and age, and (c) the level of contamination by field post-T Tauri stars on association membership surveys.

Patten, Brian M.↗
DOE OSTI2022

Role of chemical disorder on radiation-induced defect production and damage evolution in NiFeCoCr

Understanding chemical disorder in many concentrated solid solution alloys (CSAs) at the levels of electrons and atoms has attracted increasing attention as a path forward to reveal and identify underlying mechanisms for extraordinary mechanical properties and improved radiation tolerance. Single-phase NiFeCoCr CSA is a common base for many high-entropy alloys (HEAs) that have shown improved mechanical strength and radiation tolerance. In this study, defect production and damage evolution in NiFeCoCr under ion irradiation at room temperature to dose over 20 dpa are determined using ion channeling technique along both <100> and <110> directions utilizing multiple probing beam energies. The results obtained from the multi-axial and multi-energy channeling analysis are compared with those previously obtained for Ni crystals irradiated under similar conditions. The influence of chemical complexity on defect production and clustering at early-stage under room temperature irradiation up to dose of 1 dpa is discussed based on positron annihilation spectroscopy results. Defect structure evaluation in Ni and NiFeCoCr is also discussed based on transmission electron microscopy results over a prolonged irradiation at both room and elevated temperatures. Compared with chemically complex NiFeCoCr, larger dislocation loops thus less lattice strain are expected to form in pure Ni. Moreover, the role of chemical disorder in this CSA is also investigated based on ab initio calculations using large supercells. Finally, to understand the impact of chemical complexity effect on defect structure evolution, this integrated research effort attempts to link the relatively large charge redistribution due to difference in valence electron counts resulting from alloying different 3d transition metal elements, moderate lattice distortion arising from similar adaptable atomic size, and notable suppressed or delayed damage evolution in NiFeCoCr.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Detection of trace quantities of dimethylacetamide during the curing process of a polyimide film

An experimental study of positive ion clustering of dimethylacetamide (DMAC) with water molecules using high pressure mass spectrometry techniques is reported. The possibilities of using these techniques to determine the concentration levels of DMAC evolved from a polyimide film during the curing process, are examined. The test sample of gas contained DMAC of the order of a few micrograms, water, nitrogen, oxygen, and other compounds that might evolve during the curing cycle of a polyimide film. It is found that PPM concentrations such as DMAC lost from polyimide during curing stages can be detected with this technique. Improvement of the detection of lower concentrations (50 picograms) can be made by using electron multipliers with a gain in the 10 to the 8th power range and pulse counting techniques.

Smith, A. C.↗
DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
DOE OSTI2022

Revisiting the Definition of Rapid Intensification of Tropical Cyclones by Clustering the Initial Intensity and Inner‐Core Size

Abstract Rapid intensification (RI) of tropical cyclones (TCs) provides a great challenge in operational forecasting and contributes significantly to the development of major TCs. RI is commonly defined as an increase in the maximum sustained surface wind speed of at least a certain threshold within 24 hr. The most widely used threshold is 30 kt (15.4 m/s), which was determined statistically. Here we propose a new definition for RI by objectively clustering TCs using the intensification rate, initial intensity, and radius of the maximum wind speed (RMW). A group of 770 samples is separated at a threshold of 45 kt (23.2 m/s). The threshold is 40 kt (20.6 m/s) for the western North Atlantic, where TC size measurements are more reliable. Monte Carlo experiments demonstrate that the proposed threshold is robust even considering the uncertainty in RMW of as high as 30 km. We show that, when a TC undergoes RI, its maximum wind speed is approximately 60 ± 15 kt (30.9 ± 7.7 m/s) and the RMW is 45 ± 20 km. The new threshold outperforms the conventional threshold of 30 kt/24 hr in describing the bimodal distribution of lifetime maximum intensity and explaining the annual count of Category 5 TCs. This new definition provides a more physically based threshold and describes a more reliable representation of extreme events. Although more comparisons are needed for operational application, it is likely to be desirable for case‐based process studies and could provide a more valuable metric for TC intensification classification and research.

Li, Yi↗
DOE OSTI2021

ResNet and CycleGAN for pulse shape discrimination of He-4 detector pulses: Recovering pulses conventional algorithms fail to label unanimously

Pulse shape discrimination (PSD) capable detectors, such as He-4, that respond to neutron and gamma-ray 7 interactions have a threshold deposited energy value below which n/γ discrimination vanishes when using 8 conventional PSD algorithms. Recent attempts in applying supervised learning based artificial neural 9 networks for PSD use the pulses in the separated regions to train the networks so they can be used to classify 10 another set of separated pulses. In doing so, pulses previously indistinguishable are not recovered for 11 classification, which would have increased the number of neutron and gamma-ray pulses that could be used 12 for further analysis. Assuming the reason why conventional PSD algorithms have unseparated regions is 13 because the parameter space of the algorithms fail to capture the intrinsic (but subtle) distinguishing 14 behavior of some of the neutron and gamma-ray pulses, a cycle-consistent generative adversarial network 15 (CycleGAN) was trained to amplify those differences and extract well separated neutron and gamma-ray 16 clusters. Results show that, once the network is trained with pulses from separated and unseparated regions, 17 it was able to transform the pulses in the unseparated region to improve the PSD. Subsequent n/γ 18 classification was performed using deep residual network (ResNet) that takes pulses with 512 data points 19 as an input. Two different ResNets were explored – simple ResNet and modified ResNet which takes 20 segmented pulse inputs in the first layer and the corresponding time axis values in the last hidden layer. 21 The later approach enables the network to extract time correlated pulse features to enhance its ability to 22 capture the pulse behaviors relevant for PSD. Although it achieves slightly lower accuracy, 99.41% versus 23 99.89%, based on simply counting the number of correct n/γ labels assigned, compared to the simple 24 ResNet, the modified ResNets architecture was able to decreases the cross-entropy loss function by half, 25 which implies that the correct n/γ labels assigned are less likely to be accidental. PSD parameter 26 distributions based on n/γ classification by ResNet before and after transforming unseparated pulses using 27 CycleGAN show that by enhancing the separation between neutrons and gamma-rays, the transformation 28 helps improve the performance of classifier networks that are trained using labeled dataset. The 29 enhancement of neutron and gamma-ray separation by the CycleGAN increased the PSD figure of merit 30 (FOM) by up to 70% in some regions. Here, the results show that, if a given detector achieves clear separation 31 between neutron and gamma-ray pulses in any energy region, such neural network approaches can help 32 lower the energy threshold for the separation and increasing the number of neutron and gamma-ray pulses 33 that can be used for further analysis.

4He↗