Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation types and techniques”

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,063 records · Page 59

Understanding and Exploiting Spatial Properties of System Failures on Extreme-Scale HPC Systems, In: 2015 45th Annual IEEE/IFIP International Conference on Dependable Systems and Networks

As we approach exascale, the scientific simulations are expected to experience more interruptions due to increased system failures. Designing better HPC resilience techniques requires understanding the key characteristics of system failures on these systems. While temporal properties of system failures on HPC systems have been well-investigated, there is limited understanding about the spatial characteristics of system failures and its impact on the resilience mechanisms. Therefore, we examine the spatial characteristics and behavior of system failures. We investigate the interaction between spatial and temporal characteristics of failures and its implications for system operations and resilience mechanisms on large-scale HPC systems. We show that system failures have "spatial locality" at different granularity in the system, study impact of different failure-types, and investigate the correlation among different failure-types. Finally, we propose a novel scheme that exploits the spatial locality in failures to improve application and system performance. Our evaluation shows that the proposed scheme significantly improves the system performance in a dynamic and production-level HPC system.

Computer Science↗

The spectral evolution of the first Galaxies. III. Simulated James Webb Space Telescope spectra of reionization-epoch galaxies with Lyman continuum leakage

Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman-continuum escape fractions (f (esc)) at redshifts $z\approx 7$–9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hβ)–β diagram (previously proposed as a tool for identifying galaxies with very high f (esc)), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman-continuum leakage. Based on our models, we expect that essentially all $z\approx 7\mbox{--}9$ galaxies that exhibit rest-frame $\mathrm{EW}({\rm{H}}\beta )\lesssim 30$ Å to have ${f}_{\mathrm{esc}}\gt 0.5$. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of $z\gt 6$ galaxies can in principle bias the selection, but substantial model deficiencies of this type should at the same time be evident from offsets in the observed distribution of $z\gt 6$ galaxies in the EW(Hβ)–β diagram compared to the simulated distribution. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement.

79 ASTRONOMY AND ASTROPHYSICS↗

Phased array antenna matching: Simulation and optimization of a planar phased array of circular waveguide elements

A computerized simulation of a planar phased array of circular waveguide elements is reported using mutual coupling and wide angle impedance matching in phased arrays. Special emphasis is given to circular polarization. The aforementioned computer program has as variable inputs: frequency, polarization, grid geometry, element size, dielectric waveguide fill, dielectric plugs in the waveguide for impedance matching, and dielectric sheets covering the array surface for the purpose of wide angle impedance matching. Parameter combinations are found which produce reflection peaks interior to grating lobes, while dielectric cover sheets are successfully employed to extend the usable scan range of a phased array. The most exciting results came from the application of computer aided optimization techniques to the design of this type of array.

Dudgeon, J. E.↗

Understanding and Exploiting Spatial Properties of System Failures on Extreme-Scale HPC Systems

As we approach exascale, the scientific simulations are expected to experience more interruptions due to increased system failures. Designing better HPC resilience techniques requires understanding the key characteristics of system failures on these systems. While temporal properties of system failures on HPC systems have been well-investigated, there is limited understanding about the spatial characteristics of system failures and its impact on the resilience mechanisms. Therefore, we examine the spatial characteristics and behavior of system failures. We investigate the interaction between spatial and temporal characteristics of failures and its implications for system operations and resilience mechanisms on large-scale HPC systems. We show that system failures have spatial locality at different granularity in the system; study impact of different failure-types; and investigate the correlation among different failure-types. Finally, we propose a novel scheme that exploits the spatial locality in failures to improve application and system performance. Our evaluation shows that the proposed scheme significantly improves the system performance in a dynamic and production-level HPC system.

Rogers, James H↗

Evaluation of Am–Li neutron spectra data for active well type neutron multiplicity measurements of uranium

Safeguarding nuclear material is an important and challenging task for the international community. One particular safeguards technique commonly used for uranium assay is active neutron correlation counting. This technique involves irradiating unused uranium with ( α,n) neutrons from an Am-Li source and recording the resultant neutron pulse signal which includes induced fission neutrons. Although this non-destructive technique is widely employed in safeguards applications, the neutron energy spectra from an Am-Li sources is not well known. Several measurements over the past few decades have been made to characterize this spectrum; however, little work has been done comparing the measured spectra of various Am-Li sources to each other. This paper examines fourteen different Am-Li spectra, focusing on how these spectra affect simulated neutron multiplicity results using the code Monte Carlo N-Particle eXtended (MCNPX). Two measurement and simulation campaigns were completed using Active Well Coincidence Counter (AWCC) detectors and uranium standards of varying enrichment. The results of this work indicate that for standard AWCC measurements, the fourteen Am-Li spectra produce similar doubles and triples count rates. Finally, the singles count rates varied by as much as 20% between the different spectra, although they are usually not used in quantitative analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Position sensitive alpha detector for an associated particle imaging system

Associated Particle Imaging (API) is a nuclear technique that provides for the nondestructive determination of 3D isotopic distributions. The method is based on the detection of the alpha particles associated with the neutron emitted in the deuterium-tritium (DT) fusion reaction, which provides information regarding the direction and time of the emitted 14 MeV neutron. Inelastic neutron scattering leads to characteristic gamma-ray emission from certain isotopes, for example 12 C, that can be correlated with the neutron interaction location. An API system consisting of a sealed-type neutron generator, gamma detectors, and a position-sensitive alpha detector is under development for the nondestructive quantification of carbon distribution in soils. This paper describes the design of the alpha detector, detector response simulations, and first experimental results. The alpha detector consists of a Yttrium Aluminum Perovskite (YAP) scintillator mounted on the inside of a neutron generator tube. The scintillation light propagates through a sapphire window to a position-sensitive photomultiplier tube mounted on the outside. The 16x16 output signals are connected through a resistive network for a 4-corner readout. The four readout channels are amplified, filtered, and then digitized for the calculation of the alpha position. First test findings demonstrate that an excellent alpha position resolution, better than the 1 mm FWHM required by the application, can be achieved with this detector design.

47 OTHER INSTRUMENTATION↗

An experimental search for near-wall boundary conditions for large eddy simulation

Instantaneous wall shear stress and streamwise velocities have been measured simultaneously in a flat plate, turbulent boundary layer at moderate Reynolds number in an effort to provide experimental support for large eddy simulations. Data were obtained by using a buried-wire wall shear gage and a hot-wire rake positioned in the log region of the flow. All data processing was accomplished with digital data analysis techniques on a minicomputer. Fluctuations of the instantaneous U plus versus Y plus profiles about a mean law of the wall are shown to be significant and complex. Peak cross-correlation values between wall shear stress and the velocities are high and reflect the passage of a large structure inclined at a small angle to the wall. Estimates of this angle are consistent with those made by other investigators. Conditional sampling techniques were used to detect the passage of various sizes and types of flow disturbances (events) and to estimate their mean frequency of occurrence. Events characterized by large and sudden streamwise accelerations were found to be highly coherent throughout the log region and were strongly correlated with large fluctuations in wall shear-stress. Phase randomness between the near-wall quantities and the outer velocities was small. The results suggest that the flow events detected by conditional sampling applied to velocities in the log region may be related to the bursting process.

Robinson, S. K.↗

Project Title: Demonstration High Temperature Superconducting NonPlanar Stellarator Magnet with Advanced Manufactured Assemblies

This is the final report for the project “Demonstration High Temperature Superconducting Non- Planar Stellarator Magnet with Advanced Manufactured Assemblies”, funded by DOE, and performed by Type One Energy from September, 2020, to March 2024 involving the Fusion Technology Institute at the University of Wisconsin–Madison, the Plasma Science and Fusion Center (PSFC) at the Massachusetts Institute of Technology (MIT) and Commonwealth Fusion Systems (CFS) to design and fabricate the first non-planar HTS (REBCO) coil for a high-field stellarator based on the SPARC tokamak’s VIPER cable concept. Stellarators at high fields make high-temperature superconducting magnets necessary for a compact fusion device. But the asymmetric and non-planar nature of its components, especially the magnets make it difficult for scalable producibility. To address these challenges, two promising technologies have emerged: advanced manufacturing (AM) for the supporting plates for forming the magnets, and high-temperature superconducting (HTS) cables inside the plates. AM has advanced enough to produce stellarator components with the necessary geometric complexity, size, and the precision, leading to potentially significant reductions in production time, cost, and waste. The cost of HTS tape has decreased dramatically, and progress in HTS planar magnet development has reached a point where it can be proposed for application to complex 3D non-planar magnets. The main objective of this project is to develop, demonstrate and pre-commercialize a novel, non-planar HTS coil shape that remains superconducting to achieve production scalable reductions in time and cost and performance. The proposed technology is based on the novel concept of a precision sub-scale HTS nonplanar coil assembly. This project focuses on the design, fabrication, material optimization of cable design, and validation and demonstration of the high current carrying capability of superconducting magnets and their support in a complex 3D shape needed for application to stellarator magnetic plasma confinement. The specific objectives of this research program include: (1) The successful application of metal AM to build a precision sub-scale HTS nonplanar coil, (2) An HTS cable and cross-section design that can conform to the required nonplanar coil shape (bend radii as tight as 10-cm) and remains superconducting at an engineering current density of 1.35 kA/cm 2 at 77 K and 1 tesla at the conductor (5 kA in the cable). To achieve the above challenging goals, we have formed a multidisciplinary research team consisting of members from Type One Energy and UW-Madison, MIT PSFC and CFS with complementary skills and strong facilities. The team worked collaboratively on fundamental and applied research on the following three major technical areas: (1) Design, fabrication, and optimization of non-planar HTS Cable The ultimate goal of the project is to determine if commercial REBCO tapes and additive manufacturing can be used to fabricate high field (≥ 10T) non-planar coils with tight bending radii (≃ 100mm) and with a degradation of the critical current (Ic) smaller than 20% with respect to the expected performance. We started with shorter length cable to evaluate the scalability of the production process and eventually reached multiple turns for higher magnetic fields. Our findings suggest that a stellarator coil system of a relevant size, characterized by its asymmetric and non-planar components, can be fabricated using a formed cable in plate method. This system can be simulated using a large-scale modeling approach. The use of hybrid modeling 3 techniques will be pivotal in reducing the complexity of the model and in assessing expected performance in designs. (2) Modeling and simulation of the non-planar HTS Cable Multiphysics simulations are performed using the commercial software and are carried out in self-field conditions, involving 2D and 3D models and twisted around one slot of twist-pitched VIPER cable. Multiphysics simulations are mainly focused on the coil for the critical current evaluation, the magnetic field map, self-Lorentz forces and mechanical, and magnetothermal behavior and the quench dynamics. The detailed model and prediction of the superconducting performance of a stellarator-relevant demonstration cable from numerical simulations supports the results from the actual testing backing the results. A detailed description and results are provided in the later sections. (3) Design, fabrication, and optimization of support for the non-planar HTS Cable The team developed an additive manufactured (AM) coil positioning plate that formed into the required non-planar geometry (with bend radii as tight as 10-cm) and to acceptable tolerances required for a stellarator magnet: (+0.25-mm from ideal on dimensions of coil positioning plates and up to +1-mm from ideal for position of wound coil). The plate materials is also included in this selection process from fabrication and 3D printing perspective and commensurate with eventual application to a fusion reactor. From the cost effectiveness point of view, the HTS coil and plate has the potential to cost less than that made in conventional methods with less waste (<75% waste) reducing time (<50%) and cost (<50%), especially as the AM field matures. The application of advanced manufacturing in the construction of the support plates will also lead to cost reduction, as the cables can be easily replaced, thereby making the assembly modular. With such high primary cost and time savings, high current densities and magnetic field, the funded R&D work has validated the designs, proven the feasibility, and characterized the performance of the HTS coil and plate assembly, paving the way for a relevant-size stellarator coil system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three Dimensional Constraint Effects on the Estimated (Delta)CTOD during the Numerical Simulation of Different Fatigue Threshold Testing Techniques

Variation in constraint through the thickness of a specimen effects the cyclic crack-tip-opening displacement (DELTA CTOD). DELTA CTOD is a valuable measure of crack growth behavior, indicating closure development, constraint variations and load history effects. Fatigue loading with a continual load reduction was used to simulate the load history associated with fatigue crack growth threshold measurements. The constraint effect on the estimated DELTA CTOD is studied by carrying out three-dimensional elastic-plastic finite element simulations. The analysis involves numerical simulation of different standard fatigue threshold test schemes to determine how each test scheme affects DELTA CTOD. The American Society for Testing and Materials (ASTM) prescribes standard load reduction procedures for threshold testing using either the constant stress ratio (R) or constant maximum stress intensity (K(sub max)) methods. Different specimen types defined in the standard, namely the compact tension, C(T), and middle cracked tension, M(T), specimens were used in this simulation. The threshold simulations were conducted with different initial K(sub max) values to study its effect on estimated DELTA CTOD. During each simulation, the DELTA CTOD was estimated at every load increment during the load reduction procedure. Previous numerical simulation results indicate that the constant R load reduction method generates a plastic wake resulting in remote crack closure during unloading. Upon reloading, this remote contact location was observed to remain in contact well after the crack tip was fully open. The final region to open is located at the point at which the load reduction was initiated and at the free surface of the specimen. However, simulations carried out using the constant Kmax load reduction procedure did not indicate remote crack closure. Previous analysis results using various starting K(sub max) values and different load reduction rates have indicated DELTA CTOD is independent of specimen size. A study of the effect of specimen thickness and geometry on the measured DELTA CTOD for various load reduction procedures and its implication in the estimation of fatigue crack growth threshold values is discussed.

Seshadri, Banavara R.↗

Noninvasive liquid level sensing with laser generated ultrasonic waves

This article proposes a noninvasive liquid level sensing technique using laser-generated ultrasound waves for nuclear power plant applications. Liquid level sensors play an important role of managing the coolant system safely and stably in the plant structure. Current sensing techniques are mostly intrusive, performing inside the fluidic structure, which is disadvantageous in terms of the regular maintenance of the plant system. Furthermore, typical intrusive sensors do not perform stably under varying environmental conditions such as temperature and radiation. In this study, sensing units are attached to the outer surface of a liquid vessel to capture guided ultrasound waves in a nonintrusive manner. The signal intensity of the guided wave dissipates when the signal interacts with the internal liquid media. The sensing mechanism is mathematically expressed as an index value to correlate the liquid level with the sensor signal. For the acoustic wave generation, laser-generated ultrasound was adopted instead of using typical contact type transducers. Following the simulation validation of the proposed concept, the performance of the developed sensor was confirmed through experimental results under elevated liquid temperature conditions. Finally, the nonlinear multivariable regression exhibited the best-fit to the datasets measured under the variable liquid level and temperature conditions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Viscosity of α-pinene secondary organic material and implications for particle growth and reactivity

Particles composed of secondary organic material (SOM) are abundant in the lower troposphere and play important roles in climate, air quality, and health. The viscosity of these particles is a fundamental property that is presently poorly quantified for conditions relevant to the lower troposphere. Using two new techniques, namely a bead-mobility technique and a poke-flow technique, in conjunction with simulations of fluid flow, we measure the viscosity of the watersoluble component of SOM produced by α-pinene ozonolysis. The viscosity is comparable to that of honey at 90% relative humidity (RH), comparable to that of peanut butter at 70% RH and greater than or comparable to that of bitumen for ≤ 30% RH, implying that the studied SOM ranges from liquid to semisolid/solid at ambient relative humidities. With the Stokes-Einstein relation, the measured viscosities further imply that the growth and evaporation of SOM by the exchange of organic molecules between the gas and condensed phases may be confined to the surface region when RH ≤ 30%, suggesting the importance of an adsorption-type mechanism for partitioning in this regime. By comparison, for RH ≥ 70% partitioning of organic molecules may effectively occur by an absorption mechanism throughout the bulk of the particle. Finally, the net uptake rates of semi-reactive atmospheric oxidants such as O3 are expected to decrease by two to five orders of magnitude for a change in RH from 90% to ≤ 30% RH, with possible implications for the rates of chemical aging of SOM particles in the atmosphere.

Renbaum-Wolff, Lindsay↗

From solar to stellar corona: The role of wind, rotation, and magnetism

Observations of surface magnetic fields are now within reach for many stellar types thanks to the development of Zeeman–Doppler Imaging. These observations are extremely useful for constraining rotational evolution models of stars, as well as for characterizing the generation of the magnetic field. We recently demonstrated that the impact of coronal magnetic field topology on the rotational braking of a star can be parameterized with a scalar parameter: the open magnetic flux. However, without running costly numerical simulations of the stellar wind, reconstructing the coronal structure of the large-scale magnetic field is not trivial. An alternative—broadly used in solar physics—is to extrapolate the surface magnetic field assuming a potential field in the corona, to describe the opening of the field lines by the magnetized wind. This technique relies on the definition of a so-called source surface radius, which is often fixed to the canonical value of 2.5R{sub ⊙}. However this value likely varies from star to star. To resolve this issue, we use our extended set of 2.5D wind simulations published in 2015 to provide a criterion for the opening of field lines as well as a simple tool to assess the source surface radius and the open magnetic flux. This allows us to derive the magnetic torque applied to the star by the wind from any spectropolarimetric observation. We conclude by discussing some estimations of spin-down timescales made using our technique and compare them to observational requirements.

79 ASTRONOMY AND ASTROPHYSICS↗

TH-AB-209-08: Next Generation Dedicated 3D Breast Imaging with XACT

Purpose: Exposure to radiation increases the risk of cancer. We have designed a new imaging paradigm, X-ray induced acoustic computed tomography (XACT). Applying this innovative technology to breast imaging, an X-ray exposure can generate a 3D acoustic image, which dramatically reduces the radiation dose to patients when compared to conventional breast CT. Methods: Theoretical calculations are done to determine the appropriate X-ray energy and ultrasound frequency in breast XACT imaging. A series of breast CT image along the coronal plane from a patient with calcifications in the breast tissue are used as the source image. HU value based segmentation is done to distinguish the skin, adipose tissue, glandular tissue, breast calcification, and chest bone from each CT image. X-ray dose deposition in each pixel is calculated based on the tissue type by using GEANT4 Monte Carlo toolkits. The initial pressure rise caused by X-ray energy deposition is calculated according to tissue properties. Then, the X-ray induced acoustic wave propagation is simulated by K-WAVE toolkit. Breast XACT images are reconstructed from the recorded time-dependent ultrasound waves. Results: For imaging a breast with large size (16cm in diameter at chest wall), the photon energy of X-ray source and the central frequency of ultrasound detector is determined as 20keV and 5.5MHz. Approximately 10 times contrast between a calcification and the breast tissue can be acquire from XACT image. The calcification can be clearly identified from the reconstructed XACT image. Conclusion: XACT technique takes the advantages of X-ray absorption contrast and high ultrasonic resolution. With the proposed innovative technology, one can potentially reduce radiation dose to patient in 3D breast imaging as compared with current x-ray modalities, while still maintaining high imaging contrast and spatial resolution.

60 APPLIED LIFE SCIENCES↗

Numerical Estimation of the Spent Fuel Ratio

Sabotage of spent nuclear fuel casks remains a concern nearly forty years after attacks against shipment casks were first analyzed and has a renewed relevance in the post-9/11 environment. A limited number of full-scale tests and supporting efforts using surrogate materials, typically depleted uranium dioxide (DUO 2 ), have been conducted in the interim to more definitively determine the source term from these postulated events. However, the validity of these large- scale results remain in question due to the lack of a defensible spent fuel ratio (SFR), defined as the amount of respirable aerosol generated by an attack on a mass of spent fuel compared to that of an otherwise identical surrogate. Previous attempts to define the SFR in the 1980's have resulted in estimates ranging from 0.42 to 12 and include suboptimal experimental techniques and data comparisons. Because of the large uncertainty surrounding the SFR, estimates of releases from security-related events may be unnecessarily conservative. Credible arguments exist that the SFR does not exceed a value of unity. A defensible determination of the SFR in this lower range would greatly reduce the calculated risk associated with the transport and storage of spent nuclear fuel in dry cask systems. In the present work, the shock physics codes CTH and ALE3D were used to simulate spent nuclear fuel (SNF) and DUO 2 targets impacted by a high-velocity jet at an ambient temperature condition. These preliminary results are used to illustrate an approach to estimate the respirable release fraction for each type of material and ultimately, an estimate of the SFR. This page intentionally blank

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Application of multiple objective optimization techniques to finite element model tuning

This report examines tuning a finite element model using vector optimization techniques. Structural models using finite element theory often need to be adjusted so they can accurately simulate the real structure. The goal is to tune the model such that it will reproduce data derived from structural tests. First, the performance indices are extremized using multiple objective optimization theory, producing a set of possible solutions. Next, the solutions are rank ordered according to a decision maker's preferences to select the best answer. The tuning process was applied to a T-38 horizontal stabilizer. Numerous weighted solutions contained a best static deformation model, a best frequency model and three intermediate combinations of these two models. This automated procedure proved to be a versatile method capable of producing solutions for many types of tuning problems.

Devore, Charles R.↗

Climate and atmospheric modeling studies. Climate applications of Earth and planetary observations. Chemistry of Earth and environment

The research conducted during the past year in the climate and atmospheric modeling programs concentrated on the development of appropriate atmospheric and upper ocean models, and preliminary applications of these models. Principal models are a one-dimensional radiative-convective model, a three-dimensional global climate model, and an upper ocean model. Principal applications have been the study of the impact of CO2, aerosols and the solar 'constant' on climate. Progress was made in the 3-D model development towards physically realistic treatment of these processes. In particular, a map of soil classifications on 1 degree x 1 degree resolution has been digitized, and soil properties have been assigned to each soil type. Using this information about soil properties, a method was developed to simulate the hydraulic behavior of soils of the world. This improved treatment of soil hydrology, together with the seasonally varying vegetation cover, will provide a more realistic study of the role of the terrestrial biota in climate change. A new version of the climate model was created which follows the isotopes of water and sources of water (or colored water) throughout the planet. Each isotope or colored water source is a fraction of the climate model's water. It participates in condensation and surface evaporation at different fractionation rates and is transported by the dynamics. A major benefit of this project has been to improve the programming techniques and physical simulation of the water vapor budget of the climate model.

Source record↗

Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes

High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel Benchmarks (NPB). In this paper, we present some interesting performance results of ow OpenMP parallel implementation on different architectures such as the SGI Origin2000, SGI Altix, and Cray MTA-2.

Feng, Huiyu↗

A Review of Basic Energy Reconstruction Techniques in Liquid Xenon and Argon Detectors for Dark Matter and Neutrino Physics Using NEST

Detectors based upon the noble elements, especially liquid xenon as well as liquid argon, as both single- and dual-phase types, require reconstruction of the energies of interacting particles, both in the field of direct detection of dark matter (weakly interacting massive particles WIMPs, axions, etc.) and in neutrino physics. Experimentalists, as well as theorists who reanalyze/reinterpret experimental data, have used a few different techniques over the past few decades. In this paper, we review techniques based on solely the primary scintillation channel, the ionization or secondary channel available at non-zero drift electric fields, and combined techniques that include a simple linear combination and weighted averages, with a brief discussion of the application of profile likelihood, maximum likelihood, and machine learning. Comparing results for electron recoils (beta and gamma interactions) and nuclear recoils (primarily from neutrons) from the Noble Element Simulation Technique (NEST) simulation to available data, we confirm that combining all available information generates higher-precision means, lower widths (energy resolution), and more symmetric shapes (approximately Gaussian) especially at keV-scale energies, with the symmetry even greater when thresholding is addressed. Near thresholds, bias from upward fluctuations matters. For MeV-GeV scales, if only one channel is utilized, an ionization-only-based energy scale outperforms scintillation; channel combination remains beneficial. We discuss here what major collaborations use.

argon↗