Search NASA⌕ Search

SEARCH · Search NASA

Results for “core optimization”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

High Performance Data Transfer for Distributed Data Intensive Sciences

We report on the development of ZX software providing high performance data transfer and encryption. The design scales in: computation power, network interfaces, and IOPS while carefully balancing the available resources. Two U.S. patent-pending algorithms help tackle data sets containing lots of small files and very large files, and provide insensitivity to network latency. It has a cluster-oriented architecture, using peer-to-peer technologies to ease deployment, operation, usage, and resource discovery. Its unique optimizations enable effective use of flash memory. Using a pair of existing data transfer nodes at SLAC and NERSC, we compared its performance to that of bbcp and GridFTP and determined that they were comparable. With a proof of concept created using two four-node clusters with multiple distributed multi-core CPUs, network interfaces and flash memory, we achieved 155Gbps memory-to-memory over a 2x100Gbps link aggregated channel and 70Gbps file-to-file with encryption over a 5000 mile 100Gbps link.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2009

Application of Air Coupled Acoustic Thermography (ACAT) for Inspection of Honeycomb Sandwich Structures

The application of a noncontact air coupled acoustic heating technique is investigated for the inspection of advanced honeycomb composite structures. A weakness in the out of plane stiffness of the structure, caused by a delamination or core damage, allows for the coupling of acoustic energy and thus this area will have a higher temperature than the surrounding area. Air coupled acoustic thermography (ACAT) measurements were made on composite sandwich structures with damage and were compared to conventional flash thermography. A vibrating plate model is presented to predict the optimal acoustic source frequency. Improvements to the measurement technique are also discussed.

Winfree, William P.↗
DOE OSTI2026

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

Risk reduction using DDP (Defect Detection and Prevention): Software support and software applications

Risk assessment and mitigation is the focus of the Defect Detection and Prevention (DDP) process, which has been applied to spacecraft technology assessments and planning, both hardware and software. DDP's major elements and their relevance to core requirement engineering concerns are summarized. The accompanying research demonstration illustrates DDP's tool support, and further customizations for application to software.

Software-specific customizations↗
DOE OSTI2025

ON-OFF neuromorphic ISING machines using Fowler-Nordheim annealers

We introduce NeuroSA, a neuromorphic architecture specifically designed to ensure asymptotic convergence to the ground state of an Ising problem using a Fowler-Nordheim quantum mechanical tunneling based threshold-annealing process. The core component of NeuroSA consists of a pair of asynchronous ON-OFF neurons, which effectively map classical simulated annealing dynamics onto a network of integrate-and-fire neurons. The threshold of each ON-OFF neuron pair is adaptively adjusted by an FN annealer and the resulting spiking dynamics replicates the optimal escape mechanism and convergence of SA, particularly at low-temperatures. To validate the effectiveness of our neuromorphic Ising machine, we systematically solved benchmark combinatorial optimization problems such as MAX-CUT and Max Independent Set. Across multiple runs, NeuroSA consistently generates distribution of solutions that are concentrated around the state-of-the-art results (within 99%) or surpass the current state-of-the-art solutions for Max Independent Set benchmarks. Furthermore, NeuroSA is able to achieve these superior distributions without any graph-specific hyperparameter tuning. For practical illustration, we present results from an implementation of NeuroSA on the SpiNNaker2 platform, highlighting the feasibility of mapping our proposed architecture onto a standard neuromorphic accelerator platform.

42 ENGINEERING↗
DOE OSTI2020

Designing a parallel Feel-the-Way clustering algorithm on HPC systems

This paper introduces a new parallel clustering algorithm, named Feel-the-Way clustering algorithm, that provides better or equivalent convergence rate than the traditional clustering methods by optimizing the synchronization and communication costs. Our algorithm design centers on how to optimize three factors simultaneously: reduced synchronizations, improved convergence rate, and retained same or comparable optimization cost. To compare the optimization cost, we use the Sum of Square Error (SSE) cost as the metric, which is the sum of the square distance between each data point and its assigned clusters. Compared with the traditional MPI k-means algorithm, the new Feel-the-Way algorithm requires less communications among participating processes. As for the convergence rate, the new algorithm requires fewer number of iterations to converge. As for the optimization cost, it obtains the SSE costs that are close to the k-means algorithm. In the paper, we first design the full-step Feel-the-Way k-means clustering algorithm that can significantly reduce the number of iterations that are required by the original k-means clustering method. Next, we improve the performance of the full-step algorithm by adopting an optimized sampling-based approach, named reassignment-history-aware sampling. Our experimental results show that the optimized sampling-based Feel-the-Way method is significantly faster than the widely used k-means clustering method, and can provide comparable optimization costs. More extensive experiments with several synthetic datasets and real-world datasets (e.g., MNIST, CIFAR-10, ENRON, and PLACES-2) show that the new parallel algorithm can outperform the open source MPI k-means library by up to 110% on a high-performance computing system using 4,096 CPU cores. In addition, the new algorithm can take up to 51% fewer iterations to converge than the k-means clustering algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

RAFT: Reconfigurable Array of High-Efficiency Ducted Turbines for Hydrokinetic Energy Harvesting

Diversifying the energy harvesting portfolio is crucial to achieving the ambitious goal of transitioning to clean energy by 2030. Marine hydrokinetic energy has garnered renewed interest due to its high harvesting potential in the U.S., and the resource's reliability and predictability—remaining relatively constant on a daily basis and available 24/7. However, there are currently few commercial devices capable of harnessing the energy from flowing water. This project aims to bridge that gap by designing and evaluating a novel hydrokinetic turbine concept that can efficiently harvest energy from both rivers and tidal streams. The RAFT (Reconfigurable Array of High-Efficiency Ducted Turbines) concept introduces a duct surrounding the turbine rotor and creates an array of small 5-kW units. The duct serves two primary purposes: (1) it enhances hydrodynamic efficiency by accelerating flow to the rotor, and (2) it functions as a structural component, facilitating the formation of modular arrays that lower costs. This project focuses on demonstrating this concept and validating these benefits through simulations and scaled prototype testing. The project team includes 8 faculty members and over 20 students from 3 universities, organized into three core areas: hydrodynamics, electrical systems, and structural analysis, with additional teams dedicated to system integration, environmental assessment and risk management, and tech-to-market strategy. The team successfully demonstrated the increased hydrodynamic efficiency of a ducted turbine compared to an unducted version using high-fidelity simulations and prototype tests. Moreover, design optimization efforts led to surpassing the SHARKS program's goal of 60% reduction of the levelized cost of energy with a significant margin.

13 HYDRO ENERGY↗
DOE OSTI2019

Synchrotron Radiation Research and Analysis of the Particulate Matter in Deep Ice Cores: An Overview of the Technical Challenges

Airborne dust extracted from deep ice core perforations can provide chemical and mineralogical insight into the history of the climate and atmospheric conditions, with unrivalled temporal resolution, time span and richness of information. The availability of material for research and the natural complexity of the particulate, however, pose significant challenges to analytical methods. We present the developments undertaken to optimize the experimental techniques, materials and protocols for synchrotron radiation-based analysis, in particular for the acquisition of combined Synchrotron Radiation X-Ray Fluorescence and X-ray Absorption Spectroscopy data.

Cibin, Giannantonio↗
NASA NTRS2000

Development of Polymer Coatings for the ProSEDS Tether

The ProSEDS mission is designed to provide an on-orbit demonstration of the electrodynamic propulsion capabilities of tethers in space. The ProSEDS experiment will be a secondary payload on a Delta 11 unmanned, expendable booster. A 5 km conductive tether is attached to the deployer baseplate on the Delta 11 second stage and collects current from the low Earth orbit (LEO) plasma to facilitate de-orbit of the Delta II second stage. The conductive tether is attached to a 10-15 km non-conductive tether, which in turn is attached to an endmass. A bare metal tether would have the best conductivity but thermal concerns preclude this design. A conductive polymer developed by Triton Systems has been optimized for optimum conductivity and thermo-optical properties. The current design for the ProSEDS conductive tether is seven individually coated strands of 28 AWG aluminum wire, coated with 12.7 micrometers (0.5 mil) atomic oxygen-resistant conductive polymer composed of a mixture of COR and PANi, wrapped around a braided Kevlar 29 core. Extensive testing has been performed at the Marshall Space Flight Center to qualify this material for flight on ProSEDS. Atomic oxygen exposure has been performed, with solar absorptance and infrared emittance measured before and after exposure. Plasma chamber tests have been completed, as well as tether deployment tests. Also developed for the ProSEDS mission was the insulating polymer TOR-BP. Approximately 200 meters of the conductive tether closest to the Delta II second stage is insulated to prevent any electron reconnection to the tether from the plasma contactor. The insulating material is TOR-BP with a dielectric strength of TBD.

Vaughn, Jason A.↗
NASA NTRS2015

Aluminum 2195 T8 Gore Development for Space Launch System Core and Upper Stage

Gores are pie-shaped panels that are welded together to form the dome ends of rocket fuel tanks as shown in figure 1. Replacing aluminum alloy 2219 with aluminum (Al)-lithium (Li) alloy 2195 as the Space Launch System (SLS) cryogenic tank material would save enormous amounts of weight. In fact, it has been calculated that simply replacing Al 2219 gores with Al 2195 gores on the SLS core stage domes could save approximately 3,800 pound-mass. This is because the Al-Li 2195 alloy exhibits both higher mechanical properties and lower density than the SLS baseline Al 2219 alloy. Indeed, the known advantages of Al 2195 led to its use as a replacement for Al 2219 in the shuttle external tank program. The required thicknesses of Al 2195 gores for either SLS core stage tanks or upper stage tanks will depend on the specific design configurations. The required thicknesses or widths may exceed the current experience base in the manufacture of such gores by the stretch-forming process. Accordingly, the primary objective of this project was to enhance the formability of Al 2195 by optimizing the heat treatment and stretch-forming process for gore thicknesses up to 0.75 inches, which envelop the maximum expected gore thicknesses for SLS tank configurations.

Volz, Martin↗
DOE OSTI2017

A Methodology for Dynamic Security Risk Quantification and Optimal Resource Allocation of Security Assets

In this presentation we overview a methodology for dynamic security risk quantification and optimal resource allocation of security assets for high profile venues. This methodology is especially applicable to venues that require security screening operations such as mass transit (e.g., train or airport terminals), critical infrastructure protection (e.g., government buildings), and largescale public events (e.g., concerts or professional sports). The method starts by decomposing the three core components of risk -- threat, vulnerability, and consequence -- into their various subcomponents. For instance, vulnerability can be decomposed into availability, accessibility, organic security, and target hardness and each of these can be evaluated against the potential threats of interest for the given venue. Once evaluated, these subcomponents are rolled back up to compute the specific value for the vulnerability core risk component. Likewise, the same is done for consequence and threat, and then risk is computed as the product of these three components. A key aspect of our methodology is dynamically quantifying risk. That is, we incorporate the ability to uniquely allow the subcomponents and core components, and in turn, risk, to be quantified as a continuous function of time throughout the day, week, month, or year as appropriate.

Brigantic, Robert T.↗
DOE OSTI2017

Neutronics Analysis of TREAT Multi-SERTTA Calibration Test Vehicle (Multi-SERTTA CAL)

Transient testing is a necessary step in the performance evaluation of Accident Tolerant Fuels (ATF) developed as part of the Fuel Cycle Research and Development Advanced Fuels Campaign initiated and supported by the United States Department of Energy (DOE-NE). These advanced fuels are expected to be used in future reactors or replace the current fuel design used in existing power reactors for enhanced performance and safety. Upon completion of a process proscribed by the National Environmental Policy Act, the Transient Reactor Test (TREAT) facility, located at Idaho National Laboratory (INL), was chosen by DOE-NE to provide the desired specialized testing. TREAT has successfully operated from start-up in 1959 to 1994 when it was placed on operational stand-by. Careful preparations are being made for operational readiness and TREAT is projected to resume testing no later than 2018. Past TREAT testing primarily focused on sodium-cooled reactor fuels (SRF). Test vehicle designs varied to accommodate testing needs ranging from separate effects tests to integral tests approaching simulation of actual reactor environment conditions including forced convection cooling. The current focus on light water reactor (LWR) ATF requires design of new experiment vehicles to accommodate the water/steam environment and associated complexities regarding desired energy deposition into the specimen and proper containment. The first set of experiments in TREAT restart is expected to be conducted with the Multi-SERTTA (Static Environment Rodlet Transient Test Apparatuses). Difficult to quantify and calculate integral and time dependent effects inherent in transient testing limit predictive modeling capability; therefore, calibration testing in TREAT prior to the actual experiment is necessary to empirically determine the effectiveness of core-to-specimen power coupling in order to achieve the desired energy deposition into the test specimen during a nuclear power transient. This has been done in the past, in part, by evaluating the difference in energy deposition into dosimeter wires subjected to various TREAT time dependent power conditions. These calibration tests were typically conducted in a test vehicle neutronically similar to the actual test vehicle, but more simplistic in design in order to quickly retrieve and analyze dosimeter wires and minimize containment structure costs due to several trials. Prior calibration testing for SRF experiments could be conducted without coolant as sodium has no notable neutronic effect on the core-to-specimen power coupling factor (PCF). Such is not the case with LWR fuel where water serves as both moderator and coolant. This presents a design challenge as the calibration test vehicle must retain heated, high-pressure water to preserve neutronic equivalence representative of near prototypic reactor conditions while allowing for quick retrieval of dosimeters. The solution to this challenge comes in the form of a thin walled dry well positioned through the center of the basic Multi-SERTTA design. Test vehicle reactivity worth and PCF are two important factors with respect to ascertainment of neutronic equivalence. Similarity in test vehicle reactivity worth is desired to avoid the need to change operational core rod positions in TREAT which could significantly alter neutron flux distribution and consequently power distribution throughout the core and into the test specimen. The geometry and composition of various materials within the experiment containment structure can also significantly alter the quality of neutron signal delivery from the core to the test specimen. Introduction of an Inconel tube, albeit very thin, is expected to impact PCF. Simple modifications to the CAL design may compensate for the anticipated change. Evaluation of these factors and proposed modifications through neutronics analysis serves to inform experiment design and is a valuable tool towards achieving the optimal CAL design.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

National Spherical Torus Experiment Upgrade Collaborative Research on Configuration Optimization of Advanced Operating Scenarios and Control Including Macroscopy Stability EOG for Period 3/1/14 through 2/28/18

General Atomics has collaborated with PPPL since 1999 in support of NSTX and NSXT-U by leveraging DIII-D program expertise. John Ferron was the principal investigator (PI) for the first nine years (to 2008) and Rob La Haye since then. The research effort evolved over a number of three-year grants and originally included both the core (control, stability, resonant error fields) and the boundary (ELMs, divertor, error fields). In 2014 the current and last grant extension was stretched to four years and split into a core grant (advanced scenarios and control, macroscopic plasma stability and extending 3D magnetics diagnostics) and a boundary grant. Todd Evans is the PI for the new boundary grant which while also for four years expires later than the core grant; a separate end of grant (EOG) report for that is due in July 2018.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Risk-informed Physical Security Using Dynamic Force-on-Force Modelling Tools (ABSTRACT)

Physical security requirement of nuclear power plants (NPPs) in the US is regulated in a mainly prescriptive manner. This “one size fits all” approach may consequently overlook plant-specific characteristics which are beneficial for the inherent safety and security of the plant. Although a performance-based security measure has been promulgated in 10 CFR 73.55(r), current evaluation methods generally emphasize static and conservative assumptions on physical security uncertainties. This conservatism contributes to the overall operational and maintenance costs, which account to 70% of total operating expenditures for an NPP. This paper describes an ongoing work at the Idaho National Laboratory within the Light Water Reactor Sustainability pathway, which aims to optimize physical protection system with regards to system effectiveness and dollar costs. To achieve this goal, a risk-informed approach to physical security is proposed. This approach leverages the Event Modeling Risk Assessment using Linked Diagrams (EMRALD) tool to model facility infiltration and Force-on-Force (FoF) engagement in a dynamic manner. The dynamic FoF analysis was conducted by coupling human performance uncertainties modeled in EMRALD to existing commercial 3-dimensional FoF software, i.e. SCRIBE3D, AVERT and Simajin. Meanwhile, risk analysis was accomplished by coupling the attack outcomes to existing plant Probabilistic Risk Assessment (PRA) model with FLEX strategy. This approach reduces the conservatisms in security modeling and further acknowledges the possibility for plant operators to mitigate the attack outcomes. A sabotage attack on a hypothetical NPP facility was devised to demonstrate the proposed methodology. Dynamic FoF rules were implemented in EMRALD coupled to SCRIBE3D. This dynamic evaluation revealed an additional 21% system effectiveness in protecting the target set relative to the static modelling approach. Furthermore, the conditional core damage probability could be reduced by up to 46% by crediting operator actions in mitigating the sabotage outcomes. We evaluated the relative importance of each component of the physical security posture included in the model and assessed tradeoffs between components to identify potential opportunities for cost reduction while maintaining a specific level of system effectiveness. Results of the cost analysis will be presented in the conference. This study is expected to help NPP operators in allocating resources optimally to satisfy the security requirements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Carbon fiber synthesis from pitch: Insights from ReaxFF based molecular dynamics simulations

Despite the potential to lowering fabrication costs of pitch-based carbon fiber (CF) compared to polyacrylonitrile-based CF, large scale utilization is currently limited by the complexity in manufacturing. This is mainly due to the lack of consistency between the pitch structure/chemistry and CF properties under different processing conditions. Here, we employed full atomistic simulations based on reactive force fields to investigate the conversion of pitch precursors to CFs. Additionally, it was found that at the carbonization stage, the consumption rates of precursors, gas formations, and solid yields (carbonized CFs) are all dependent on the composition and arrangement of molecular moieties in the precursor. Specifically, precursors with sp 3 carbons in the cata-condensed core area have limited decomposition pathways (i.e., limited reaction sites), which leads to the formation of carbonized CFs with similar compositions. Furthermore, carbonized CFs possess distinct structural characteristics, and finally lead to the formation of CFs with different moduli. By analyzing CF crystallinity and elemental composition, general rationales are proposed for the selection of optimal precursors and processing conditions, and thus shed lights on parameter screening during industrial CF synthesis. For instance, inclusion of sp 3 carbons in cata- condensed rings can be preferred over the presence of sp 3 carbons on side chains.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2018

Structural basis for MTR4–ZCCHC8 interactions that stimulate the MTR4 helicase in the nuclear exosome-targeting complex

The nuclear exosome-targeting (NEXT) complex functions as an RNA exosome cofactor and is involved in surveillance and turnover of aberrant transcripts and noncoding RNAs. NEXT is a ternary complex composed of the RNA-binding protein RBM7, the scaffold zinc-knuckle protein ZCCHC8, and the helicase MTR4. While RNA interactions with RBM7 are known, it remains unclear how NEXT subunits collaborate to recognize and prepare substrates for degradation. Here, we show that MTR4 helicase activity is enhanced when associated with RBM7 and ZCCHC8. While uridine-rich substrates interact with RBM7 and are preferred, optimal activity is observed when substrates include a polyadenylated 3' end. We identify a bipartite interaction of ZCCHC8 with MTR4 and uncover a role for the conserved C-terminal domain of ZCCHC8 in stimulating MTR4 helicase and ATPase activities. A crystal structure reveals that the ZCCHC8 C-terminal domain binds the helicase core in a manner that is distinct from that observed for Saccharomyces cerevisiae exosome cofactors Trf4p and Air2p. Furthermore, our results are consistent with a model whereby effective targeting of substrates by NEXT entails recognition of elements within the substrate and activation of MTR4 helicase activity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A.I. for nuclear physics

Nuclear science is concerned with the understanding of the nature of matter, its basic constituents and their interaction to form the elements and the properties we observe. This includes the forms of matter we see around us and also exotic forms such as those that existed in the first moments after the Big Bang and that exist today inside neutron stars. The techniques, tools, and expertise needed for nuclear physics (NP) research are therefore diverse in nature. State-of-the art accelerators are being developed to illuminate the dynamical basis of the core of the atom in terms of the fundamental constituents called quarks and gluons and to increase the number of isotopes with known properties. This scientific infrastructure is reaching scales and complexities that require computational methods for tasks such as anomaly detection in operational data. New methodologies are needed to detect anomalies and to optimize operating parameters, predict failures as well as to discover new optimization algorithms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗