Search NASA⌕ Search

SEARCH · Search NASA

Results for “version”

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 811 records · Page 45

Symplectic neural network and its application to charged particle dynamics in electromagnetic fields

Recently, machine learning models have shown many successes in various applications in science and technology. In this work, we focus on the charged particle dynamics, with the development of a class of symplectic neural networks, including a linear version, SympMat, and a nonlinear version, HénonNet. Both are designed to preserve the structure of Hamiltonian systems. We show that they can be used to model relevant Hamiltonian systems of interest in plasma physics and astrophysics, for linear and nonlinear charged particle dynamics, with the potential to bridge multi-scale simulations. These symplectic neural networks are adapted to the applications in plasma simulations and particle-wave interaction with parametric dependence and periodicity, where we have investigated their performance and accuracy. In particular, SympMat is shown to outperform the traditional Boris particle pusher down to the sub-gyroperiod scale in the case of charged particles in uniform magnetic fields. HénonNet successfully predicts the hot electron distribution, which is validated against theoretical results. These results highlight the potential of symplectic neural networks as a trajectory integrator for particle-in-cell simulations or a fast surrogate to replace conventional numerical schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Restoring permutational invariance in the Jordan–Wigner transformation

The Jordan–Wigner transformation is a powerful tool for converting systems of spins into systems of fermions, or vice versa. While this mapping is exact, the transformation itself depends on the labelling of the spins. One consequence of this dependence is that approximate solutions of a Jordan–Wigner-transformed Hamiltonian may depend on the (physically inconsequential) labelling of the spins. In this work, we turn to an extended Jordan–Wigner transformation which remedies this problem and which may also introduce some correlation atop the Hartree–Fock solution of a transformed spin Hamiltonian. We demonstrate that this extended Jordan–Wigner transformation can be thought of as arising from a unitary version of the Lie algebraic similarity transformation (LAST) theory. Here, we show how these ideas, particularly in combination with the standard (non-unitary) version of LAST, can provide a potentially powerful tool for the treatment of the XXZ and J 1 –J 2 Heisenberg Hamiltonians.

74 ATOMIC AND MOLECULAR PHYSICS↗

Comparison of state-interaction and spinor-representation calculations of spin-orbit coupling within exact two-component coupled-cluster theories

Here a benchmark study of state-interaction and spinor-representation calculations of spin-orbit coupling using the exact two-component Hamiltonians with atomic mean-field integrals (the X2CAMF schemes) at the equation-of-motion coupled-cluster singles and doubles level is reported. We adopt a version of the X2CAMF scheme with spin-orbit integrals correct to first order for the state-interaction calculations and a version correct to infinite order for the spinor-representation calculations. The differences between the state-interaction calculations with minimum active spaces to account for (quasi-)degeneracy and the spinor-representation calculations thus correspond to higher-order spin-orbit contributions. This state-interaction approach with scalar-relativistic effects accurately included in the reference functions and in the spin-orbit integrals is shown to exhibit robust performance for elements across the periodic table. On the other hand, the more rigorous spinor representation shows more rapid convergence with respect to the number of correlated electrons and is the preferred choice for accurate calculations for heavy elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quasiclassical sampling and Wigner sampling of initial vibrational coordinates and momenta for polyatomic molecules in Monte Carlo molecular dynamics simulations

In a quasiclassical trajectory simulation, the vibrational modes are initialised with quantised vibrational energies, but vibrational phases are sampled by Monte Carlo. This requires an algorithm to assign coordinates and momenta to the various atoms. In this work, we present two methods for implementing this for nonrotating polyatomic molecules, namely, fixed-energy vibrational-state-selected initial conditions and thermal initial conditions. We also present a method for initiating classical trajectories with a ground-state Wigner distribution. These vibrational treatments are sufficient to initialise trajectories for unimolecular processes, and we also show how they can be applied to simulate bimolecular collision processes. The treatments of unimolecular and bimolecular collision processes are available in two Python codes called wigner_state_selected.py and bimolecular_collision.py, respectively, which will generate initial condition files that are recognisable by the SHARC and SHARC-MN computer programs for dynamics calculations. Both codes are available as standalone programs, as well as being included in SHARC-MN, and they will be included in future versions of SHARC. Here, the methods implemented in these codes are mostly also available in the ANT computer program, and those that are not available in ANT will be incorporated in future versions of ANT.

Wigner distribution↗

Analysis of AP1000 Small-Break Loss-of-Coolant Accident Using Reactor Transient Simulator

The Westinghouse Electric Company’s Advanced Passive Reactor (AP1000) is characterized by the incorporation of passive safety systems (PSSs) designed to ensure core cooling during transient events. The assessment of PSSs requires evaluation of their performance through a combination of experiments and simulations employing various thermal-hydraulic codes. In addition, detailed evaluation of PSSs for a specific reactor system transient analysis such as loss-of-coolant-accident analysis supports understanding representative integral effects test facility development and the further evolution model development and assessment process. Developing a reactor system code is a complex and time-consuming process that requires significant engineering expertise and effort. It can take several months to even years to complete in the early stages of reactor system design and analysis. However, this process can be expedited through the use of transient simulator models for similar reactor systems, which can be used for lesson learning and training purposes. This study uses the Personal Computer Transient Analyzer (PCTRAN) code. The main advantage of PCTRAN is its ease of use and ability to run faster than real time. This study presents the results obtained for a small-break loss-of-coolant accident (SBLOCA) for two breaks using the full version (licensed) of PCTRAN. The purpose of this investigation is to evaluate the overall system behavior during the postulated SBLOCA event as well as assess the capability of the PCTRAN code to reproduce the system response during transient events. The obtained results were compared with the Westinghouse NOTRUMP system code. The PCTRAN code proved to be reliable in predicting the qualitative behavior of the system in both transient cases. As for the system response, it was found that it is contingent on the activation time of the PSSs. The differences in reactor coolant system pressure between the two codes were attributed to the critical flow model and simplification of mass and energy balance. Despite PCTRAN’s limitations, it can still provide a reasonable prediction of various reactor parameters such as pressure, mass flow rate, and void fraction during a SBLOCA scenario. It is worth noting that PCTRAN currently employs a bulk approach similar to that of the Modular Accident Analysis Program (MAAP) and MELCOR codes. However, the upcoming version of PCTRAN will include an artificial intelligence–based detection and accident prevention system, as well as different models for different reactor components. Consequently, PCTRAN has the potential to be upgraded to match the system thermal-hydraulic codes of the U.S. Nuclear Regulatory Commission and become more widely used in cybersecurity to safeguard nuclear power plants from cyberattacks.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validation of Fast Reactor Depletion Tools Using EBR-II Measured Data

The validation of simulation tools for calculating fuel depletion and evolution in fast reactors is vital for design, licensing, deployment, operations, and material accountancy. The Physics Analysis Database (PADB) and Analytical Laboratory (AL) database contain measured data collected from Experimental Breeder Reactor II (EBR-II) and were used to validate the most recent versions of the Argonne Reactor Computation (ARC) tool suite and ORIGEN-S for calculating isotopic compositions in irradiated fast reactor fuel. The PADB contains important modeling and operational information about the EBR-II core design, fuel cycle, and analytical results from the legacy versions of the ARC tool suite. The AL database contains the measured isotopic compositions of irradiated samples taken from core subassemblies. A new procedure was developed for the ARC tool suite to perform the EBR-II depletion simulation, as well as to perform more detailed isotopic calculations using ORIGEN-S calculations by coupling it with the ARC suite. Both the ARC and ARC-ORIGEN results were compared with the AL measured data for all relevant samples and showed good agreement for the major actinides. Good agreement with measured data was also achieved using the ARC-ORIGEN approach for several fission products.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A New Shutdown Dose Rate Benchmark Problem for Representative Fusion Applications

Here, this work introduces a new benchmark problem for calculating shutdown dose rates (SDDRs) aimed at fusion reactor applications. The model is designed to represent a simplified version of a typical ITER port plug. The responses of interest include neutron flux, gamma flux, and gamma SDDR at 12 different locations scattered throughout the port. This article outlines the geometry specifications of the problem, provides material definitions for the components, specifies the required responses to be calculated, and presents the source definition information. The need for this benchmark arises from the limited availability of publicly accessible references, with only one benchmark representing the typical dimensions and materials found in fusion systems. This existing benchmark has been cited extensively, reflecting the demand within the scientific community to test both established and novel workflows for SDDR calculations. However, since its presentation at a conference in 2011, the results have become increasingly well known. Moreover, the absence of formal publication and peer review has led to the details of this benchmark being extracted from secondary sources, such as subsequent studies that reference it. As a result, analysts are left with significant flexibility in interpreting the key parameters, which can be adjusted to account for unknown systematic errors, ultimately reproducing the already well-known responses. This new benchmark serves as an updated version of that earlier work, with the aim of providing a more reliable description of the materials and their impurities, which is crucial for assessing activation and subsequent gamma emission. Additionally, it seeks to provide a geometry that more closely represents an ITER port plug. The improvements in the problem definition will lead to a more reproducible benchmark problem, while also presenting the radiation transport community with a completely new challenge. The results will be published in a future article to allow analysts adequate time to analyze this problem independently.

Benchmark↗

Isodrastic magnetic fields for suppressing transitions in guiding-centre motion

In a magnetic field, transitions between classes of guiding-centre motion can lead to cross-field diffusion and escape. We say a magnetic field is isodrastic if guiding centres make no transitions between classes of motion. This is an important ideal for enhancing confinement. First, we present a weak formulation, based on the longitudinal adiabatic invariant, generalising omnigenity. To demonstrate that isodrasticity is strictly more general than omnigenity, we construct weakly isodrastic mirror fields that are not omnigenous. Then we present a strong formulation that is exact for guiding-centre motion. We develop a first-order treatment of the strong version via a Melnikov function and show that it recovers the weak version. The theory provides quantification of deviations from isodrasticity that can be used as objective functions in optimal design. The theory is illustrated with some simple examples.

97 MATHEMATICS AND COMPUTING↗

Neural network model of the multi-mode anomalous transport module for accelerated transport simulations

Here, a neural network version of the Multi-Mode Anomalous Transport Module, known as MMMnet, has been developed to calculate plasma turbulent diffusivities in DIII-D with a calculation time suitable for control applications. MMMnet uses a simple artificial neural network structure to predict the ion thermal, electron thermal, and toroidal momentum diffusivities while reproducing Multi-Mode Model (MMM) data with good accuracy and keeping the calculation time as a fraction of that associated with MMM. Model-based control techniques require models with fast calculation times, making many existing physics-oriented predictive codes unsuitable. The control-oriented predictive code COTSIM (Control Oriented Transport SIMulator) calculates the most significant plasma dynamics in response to the different actuators while running at a speed useful for control design. In order to achieve this calculation speed, COTSIM often relies on scaling laws and control-level models. Replacing some of these scaling laws and control-level models with neural network versions of more complex physics-level models has the potential of increasing the range of validity and the level of accuracy of COTSIM without compromising its computational speed. In this work, MMMnet is integrated into COTSIM to improve the turbulent diffusivity predictions, which will in turn improve the prediction accuracy associated with the dynamics of many plasma properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neural net modeling of equilibria in NSTX-U

Neural networks (NNs) offer a path towards synthesizing and interpreting data on faster timescales than traditional physics-informed computational models. In this work we develop two NNs relevant to equilibrium and shape control modeling, which are part of a suite of tools being developed for the National Spherical Torus Experiment-Upgrade for fast prediction, optimization, and visualization of plasma scenarios. The networks include Eqnet, a free-boundary equilibrium solver trained on the EFIT01 (Equilibrium FITtting 01) reconstruction algorithm, and Pertnet, which is trained on the Gspert code and predicts the non-rigid plasma response, a nonlinear term that arises in shape control modeling. The NNs are trained with different combinations of inputs and outputs in order to offer flexibility in use cases. In particular, Eqnet can use magnetic diagnostics as inputs and act as an EFIT-like reconstruction algorithm, or, by using pressure and current profile information the NN can act as a forward Grad–Shafranov equilibrium solver. This forward-mode version is envisioned to be implemented in the suite of tools for simulation of plasma scenarios. The reconstruction-mode version gives some performance improvements compared to the online reconstruction code real-time EFIT, especially when vessel eddy currents are significant. Here, we report strong performance for all NNs indicating that the models could reliably be used within closed-loop simulations or other applications. Some limitations are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling and analysis of the tritium fuel cycle for ARC- and STEP-class D-T fusion power plants

Abstract The limited tritium resources available for the first fusion power plants (FPPs) make fuel self-sufficiency and tritium inventory minimization leading issues in FPP design. This work builds on the model proposed by Abdou et al (2020 Nucl. Fusion 61 013001), which analyzed the fuel cycle (FC) of a DEMOnstration nuclear FPP-class FPP with a time-dependent system-level model. Here, we use a modified version of their model to analyze the FC of an Affordable, Robust, Compact (ARC)-class tokamak and two versions of a Spherical Tokamak for Energy Production (STEP)-class tokamak. The ARC-class tokamak breeds tritium in a 2LiF + BeF 2 liquid immersion blanket, while the STEP-class tokamak breeds tritium utilizing either a liquid-lithium blanket design or an encapsulated breeding blanket. A time-dependent system-level model is developed in Matlab Simulink ® to simulate the evolution of tritium flows and tritium inventories in the FC. The main goals of this work are to assess tritium self-sufficiency of the ARC- and STEP-class designs and to determine quantitative design requirements that can be used to analyze the adequacy of a proposed FC system. These design requirements are aimed at achieving a low tritium inventory doubling time ( t d ) and a low start-up inventory ( I s t a r t u p ) while keeping the required tritium breeding ratio (TBR r ) as low as possible. We also consider how improvements in FC technology and POs affect TBR r and I s t a r t u p . The model results show that TBR r for ARC- and STEP-class FPPs should be achievable if the tritium burn efficiency (TBE) reaches 0.5%–1% (TBR r < 1.2). This assumes significant, but attainable, improvements over current abilities. However, the model results indicate that an FPP must achieve ambitious performance targets, including FPP availability > 70%, tritium processing time < 4 h, and the implementation of direct internal recycling (DIR). If future research yields major improvements to achievable TBE, it may be possible to achieve tritium self-sufficiency while operating at lower availability and without implementing DIR.

Physics↗

OpenFAST simulation of floating wind turbines with large heading change

Previous versions of OpenFAST, the physics-based wind turbine engineering and design tool developed by the National Renewable Energy Laboratory, were limited to small rotations of the floating platform. This prevented OpenFAST from being used to simulate important events, such as the loss of a mooring line, or specific floater concepts that might experience large platform yaw motion. To overcome this limitation, we modify the structural dynamics and hydrodynamics modules of OpenFAST to allow unrestricted platform yaw motion. We apply the improved version of OpenFAST to simulate the drifting of a floating wind turbine system after the loss of a mooring line. The results, including both global motion and the internal structural loads at selected locations, appear credible and consistent with expectations.

17 WIND ENERGY↗

Developments in Performance and Portability of BlockGen

For more than a decade Monte Carlo event generators with the current matrix element algorithms have been used for generating hard scattering events on CPU platforms, with excellent flexibility and good efficiency. While the HL-LHC is approaching and precision requirements are becoming more demanding, many studies have been made to solve the bottleneck in the current Monte Carlo event generator tool chains. The novel BlockGen family of fast matrix element algorithms shown in this report, is one of the new developments that are more suitable for GPU acceleration. We report the development experience of porting BlockGen using Kokkos. Moreover, we discuss the performance of the Kokkos version in comparison with the dedicated GPU version in CUDA.

Bothmann, E. [Gottingen U.]↗

Optimizing the hit finding algorithm for liquid argon TPC neutrino detectors using parallel architectures

Neutrinos are particles that interact rarely, so identifying them requires large detectors which produce lots of data. Processing this data with the computing power available is becoming even more difficult as the detectors increase in size to reach their physics goals. Liquid argon time projection chamber (LArTPC) neutrino experiments are expected to grow in the next decade to have 100 times more wires than in currently operating experiments, and modernization of LArTPC reconstruction code, including parallelization both at data- and instruction-level, will help to mitigate this challenge. The LArTPC hit finding algorithm is used across multiple experiments through a common software framework. In this paper we discuss a parallel implementation of this algorithm. Using a standalone setup we find speedup factors of two times from vectorization and 30–100 times from multi-threading on Intel architectures. The new version has been incorporated back into the framework so that it can be used by experiments. On a serial execution, the integrated version is about 10 times faster than the previous one and, once parallelization is enabled, more speedups comparable to the standalone program are achieved.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DUNE Photon Detection System

The Deep Underground Neutrino Experiment (DUNE) aims to provide a broad physics program primarily addressed to probing CP violation in the neutrino sector and identifying the neutrino mass hierarchy. The search for proton decay, the observation of supernova neutrino bursts, and the investigation of solar neutrinos represent other additional goals of DUNE experiment, which can be enhanced by the use and the high performances of the Photon Detection System (PDS). Using the technology based on liquid argon Time Projection Chamber (LArTPC), the experiment plans to observe neutrino interactions inside detectors located 1300 km away from the Long Baseline Neutrino Facility (LBNF) at Fermi National Accelerator Laboratory (FNAL), where the neutrinos are produced. The experiment consists of two main parts, respectively the far and near detectors. The far site will comprise four detector modules. The first module is constituted by a vertical drift single-phase LArTPC, while the second module provides a horizontal drift single-phase LArTPC. The configurations of the other modules are still under definition. Neutrino detection in LArTPCs is achieved by identifying the charge and light generated from its interactions with liquid argon. The wire planes of the instrumented anode and the PDS detect these signals, respectively. The PDS in the first two modules, detailed in the present document, uses a modified version of the so-called ARAPUCA technology, named X-ARAPUCA. This system consisting of a highly reflecting box with an entrance window made by dichroic filters and wavelength shifters, creates a trap to detect the VUV (128 nm) scintillation photons. The X-ARAPUCA of the first module is called Supercell, and it has dimensions of 488×100 mm 2 , while Megacell is the second module version, with an active area of 60×60 cm 2 . This latter configuration also represents a significant technological advancement. Since half of the modules are placed on the cathode at high voltage, they are powered and read out using innovative power-over-fiber (PoF) and signal-over-fiber (SoF) technologies. Meanwhile, the other half are installed in a membrane behind the field cage, with a total transparency of around 70%.

Neutrino detectors↗

Identification of tau leptons using a convolutional neural network with domain adaptation

A tau lepton identification algorithm,DeepTau, based on convolutional neural network techniques, has been developed in the CMS experiment to discriminate reconstructed hadronic decays of tau leptons (τ h ) from quark or gluon jets and electrons and muons that are misreconstructed as τ h candidates. The latest version of this algorithm, v2.5, includes domain adaptation by backpropagation, a technique that reduces discrepancies between collision data and simulation in the region with the highest purity of genuine τh candidates. Additionally, a refined training workflow improves classification performance with respect to the previous version of the algorithm, with a reduction of 30–50% in the probability for quark and gluon jets to be misidentified as τ h candidates for given reconstruction and identification efficiencies. This paper presents the novel improvements introduced in theDeepTau algorithm and evaluates its performance in LHC proton-proton collision data at √(s) = 13 and 13.6 TeV collected in 2018 and 2022 with integrated luminosities of 60 and 35 fb -1 , respectively. Techniques to calibrate the performance of the τ h identification algorithm in simulation with respect to its measured performance in real data are presented, together with a subset of results among those measured for use in CMS physics analyses.

Large detector-systems performance↗

Arctic extreme precipitation changes from 1980 to 2022 in response to sea ice decline and enhanced atmospheric rivers

Arctic extreme precipitation (EP) broadly impacts permafrost degradation, glacier and snow cover changes, and ice sheet mass balance as well as ecosystems. However, investigation of EP spatiotemporal variations over the Arctic remains challenging, and their primary drivers are still poorly understood. Performance estimation of three state-of-the-art reanalysis products (Climate Forecast System Reanalysis (CFSR), European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA-5), and Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2)) against gauge-based precipitation observations reveals that MERRA-2 outperforms other reanalysis for Arctic EP changes. Based on MERRA-2 data, both annual EP amount and occurrence days averaged over the Arctic show statistically significant positive trends during 1980–2022 (3.37 ± 1.03 mm dec −1 and 0.42 ± 0.17 d dec −1 , respectively), with the most pronounced increase in the autumn. Spatial heterogeneity in annual and seasonal EP trends is found across the Arctic, with the largest positive annual trends of 30 mm dec −1 over the Bering Sea and the Denmark Strait. The significant EP increase is closely associated with intensified atmospheric rivers (ARs) and widespread decline in sea ice concentration (SIC). Specially, SIC and ARs are responsible for 12% and 50% of Arctic EP inter-annual variance, respectively, while ARs directly contribute 28.3% of the total annual EP amounts. These findings explain the mechanistic controls on Arctic EP, providing critical insights for projecting polar weather and climate extremes and their impacts on the Arctic environment.

atmospheric rivers↗

Kicking the can down the road: understanding the effects of delaying the deployment of stratospheric aerosol injection

Abstract Climate change is a prevalent threat, and it is unlikely that current mitigation efforts will be enough to avoid unwanted impacts. One potential option to reduce climate change impacts is the use of stratospheric aerosol injection (SAI). Even if SAI is ultimately deployed, it might be initiated only after some temperature target is exceeded. The consequences of such a delay are assessed herein. This study compares two cases, with the same target global mean temperature of ∼1.5° C above preindustrial, but start dates of 2035 or a ‘delayed’ start in 2045. We make use of simulations in the Community Earth System Model version 2 with the Whole Atmosphere Coupled Chemistry Model version 6 (CESM2-WACCM6), using SAI under the SSP2-4.5 emissions pathway. We find that delaying the start of deployment (relative to the target temperature) necessitates lower net radiative forcing (−30%) and thus larger sulfur dioxide injection rates (+20%), even after surface temperatures converge, to compensate for the extra energy absorbed by the Earth system. Southern hemisphere ozone is higher from 2035 to 2050 in the delayed start scenario, but converges to the same value later in the century. However, many of the surface climate differences between the 2035 and 2045 start simulations appear to be small during the 10–25 years following the delayed SAI start, although longer simulations would be needed to assess any longer-term impacts in this model. In addition, irreversibilities and tipping points that might be triggered during the period of increased warming may not be adequately represented in the model but could change this conclusion in the real world.

Brody, Ezra (ORCID:000900030008681X)↗