SEARCH · Engineering Papers

Results for “Simulation framework”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,675 records · Page 93

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Human-Building Collaboration: Toward lighting enabled collaborative system design

The profession of lighting design is evolving as contemporary lighting systems increasingly rely on cutting-edge computational technologies, sensors, and IOT systems. This trend requires designers to incorporate ideas of automation, dynamic controls, and user-system interaction into their design logic. However, real-time lighting simulations are constrained by inherent limitations arising from the reductive assumptions inevitably introduced in simulated lighting environments. This raises the question of how can designers account for the discrepancies between simulated and real lighting environments, and how can collaboration between humans and autonomous lighting systems bridge this gap. To address this question, we propose a protocol for designing collaborative interactions between humans and systems. Furthermore, this protocol builds on the Human-Lighting System Interaction Framework and demonstrates how human and system intelligence can be combined to fine-tune lighting qualities in a given space. Our paper shows how interactive lighting systems can customize lighting based on user preferences in real-time and how global lighting configurations can be adjusted over time. Specifically, we demonstrate: a) Human-system collaboration assumptions and goals, as well as how the protocol can be integrated into digitally programmable lighting systems. b) Implementations of collaboration that reveal how system autonomy, performance, and user experience are improved over short and long-term timeframes. c) How lighting design can be enhanced beyond simulation-driven design optimization capacities. The associated affordances and limitations are discussed with respect to existing lighting simulation design frameworks.

autonomous systems↗

NASA NTRSDate not supplied

Jet Noise Prediction for Chevron Nozzle Flows with Wall-Modeled Large-Eddy Simulation

This paper presents results from ongoing research on jet noise prediction with wall-modeled large eddy simulations (WMLES) performed with the LAVA computational framework. In particular the focus of this study is on mixing enhancements from a single stream chevron nozzle at Reynolds number of 1×10(exp 6). Although the concept of chevron nozzles to reduce jet noise is not new, our understanding of its impact on the overall noise is still not well understood. As a first step towards predicting noise reduction due to mixing enhancement concepts from first principles with WMLES, we simulate the noise generated by a single stream chevron nozzle SMC001 as well as its equivalent axisymmetric round jet nozzle SMC000. Detailed comparisons are made with a dedicated experiment conducted at NASA’s Glenn Research Center and good agreement was achieved. Two different approaches to introduce a turbulent boundary layer were compared but show no major impact on the results. This is especially important given future work were multi-stream nozzles are considered and extended costs of resolving the internal BL would have a bigger cost impact. A permeable Ffowcs Williams Hawkings (FWH) surface enclosing the jet is used to predict far-field noise from the simulated flow-field and excellent comparison to microphone array measurements is achieved within the resolved frequency band. Sensitivity of far-field noise predictions to grid resolution is systematically documented. Near-field comparisons to PIV data shows great agreement for both velocity and normal stresses, however a systematic TKE overshoot at the nozzle exit is seen in the shear-layer. The paper shows a cost reduction of an order of magnitude compared to an earlier study of this configuration due to algorithmic and software improvements and demonstrates that WMLES can be used as a cost-competitive approach for jet noise predictions.

CST↗

DOE OSTI2021

Simulations of divertor heat flux width using transport code with cross-field drifts under the BOUT++ framework

The fluid transport code [trans-electric field (Er) module] under the BOUT++ framework has been used to simulatedivertor heat flux width and boundary Er with all drifts and the sheath potential in the scrape-off layer. The calculated steady state radial Er in the pedestal region has been compared with that of experimental measurements from the Alcator C-Mod tokamak. The magnitude and shape of Er are similar to those of the experimental data. In order to understand the relative role of cross-field drifts vs turbulent transport in setting the heat flux width, four C-Mod enhanced Dα H-mode discharges with a lower single null divertor configuration should be simulated. BOUT++ transport simulations with cross-field drifts included yield similar heat flux width λq to that of experimental measurements (within a factor of 2) from both the probe and the surface thermocouple diagnostics and show a similar trend with plasma current to that of the Eich experimental scaling. The simulations show that both drifts and turbulent transport compete to determine the heat flux width. The magnetic drifts play a dominant role in setting the divertor heat-flux width, while the E × B drift decreases the heat flux width by 10%–25%, leading to improved agreement with the experiment relative to Goldston’s model. A turbulence diffusivity scan (χ) identifies two distinct regimes: a drift dominant regime when χ is small and a turbulence dominant regime when χ is large. The Goldston heuristic drift model yields a lower limit of the width λq.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Optimal Power Hardware-in-the-Loop Interfacing: Applying Modern Control For Design and Verification of High-Accuracy Interfaces

Here, we develop a novel approach to designing power hardware-in-the-loop (PHIL) simulation interfaces that maximizes simulation bandwidth and accuracy by leveraging a modern control framework that explicitly considers objectives on accuracy and can automatically synthesize an optimal controller that meets these objectives while stabilizing the closed-loop system. The method developed improves upon common approaches to PHIL interface design that typically involve multiple design steps for manual compensation and stabilization design that can result in interfaces that are stable but have suboptimal bandwidth and accuracy. The approach developed is general and can be applied to most PHIL system configurations. The modeling framework allows for the inclusion of scaling factors and hardware-under-test current injection models that are common to practical PHIL simulations and have significant impacts on stability. We also present practical methods and metrics for verifying the absolute accuracy of PHIL simulation results without relying on relative comparisons to potentially inaccurate models or previous simulation results. We demonstrate the accuracy evaluation method and show the improved performance achievable when using the optimal PHIL interface design approach in an experimental case study involving a 100-kVA battery inverter.

42 ENGINEERING↗

DOE OSTI2023

As a Matter of Dynamical Range – Scale Dependent Energy Dynamics in MHD Turbulence

Abstract Magnetized turbulence is ubiquitous in many astrophysical and terrestrial plasmas but no universal theory exists. Even the detailed energy dynamics in magnetohydrodynamic (MHD) turbulence are still not well understood. We present a suite of subsonic, super-Alfvénic, high plasma beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scale forcing and dissipation scales, and their dissipation mechanism (implicit large eddy simulation, ILES, and direct numerical simulation (DNS)). Using an energy transfer analysis framework we calculate the effective numerical viscosities and resistivities, and demonstrate that all ILES calculations of MHD turbulence are resolved and correspond to an equivalent visco-resistive MHD turbulence calculation. Increasing the number of grid points used in an ILES corresponds to lowering the dissipation coefficients, i.e., larger (kinetic and magnetic) Reynolds numbers for a constant forcing scale. Independently, we use this same framework to demonstrate that—contrary to hydrodynamic turbulence—the cross-scale energy fluxes are not constant in MHD turbulence. This applies both to different mediators (such as cascade processes or magnetic tension) for a given dynamical range as well as to a dependence on the dynamical range itself, which determines the physical properties of the flow. We do not observe any indication of convergence even at the highest resolution (largest Reynolds numbers) simulation at 2048 3 cells, calling into question whether an asymptotic regime in MHD turbulence exists, and, if so, what it looks like.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

A Fast Framework for Generating Radioactive Mixture Spectra and Its Application to Remote High-Performance Mixture Identification

Remote detection of radioactive materials in mixtures using handheld or portal detectors remains a challenge because of factors such as low concentration, environmental interference, sensor noise, and other complications. This work introduces a fast framework for generating realistic mixture spectra. Moreover, we present mixture isotope identification using data generated by the fast framework. Researchers have examined a range of conventional and recent algorithms within the fields of machine learning and deep learning. An application to uranium enrichment-level prediction has been included. Extensive simulation experiments validated the efficacy of the proposed framework.

GADRAS↗

DOE OSTI2020

Mean Climate and Tropical Rainfall Variability in Aquaplanet Simulations Using the Model for Prediction Across Scales‐Atmosphere

Abstract Aquaplanet experiments are important tools for understanding and improving physical processes simulated by global models; yet, previous aquaplanet experiments largely differ in their representation of subseasonal tropical rainfall variability. This study presents results from aquaplanet experiments produced with the Model for Prediction Across Scales‐Atmosphere (MPAS‐A)—a community model specifically designed to study weather and climate in a common framework. The mean climate and tropical rainfall variability simulated by MPAS‐A with varying horizontal resolution were compared against results from a recent suite of aquaplanet experiments. This comparison shows that, regardless of horizontal resolution, MPAS‐A produces the expected mean climate of an aquaplanet framework with zonally symmetric but meridionally varying sea‐surface temperature. MPAS‐A, however, has a stronger signal of tropical rainfall variability driven by convectively coupled equatorial waves. Sensitivity experiments with different cumulus parameterizations, physics packages, and vertical grids consistently show the presence of those waves, especially equatorial Kelvin waves, in phase with lower‐tropospheric convergence. Other models do not capture such rainfall‐kinematics phasing. These results suggest that simulated tropical rainfall variability depends not only on the cumulus parameterization (as suggested by previous studies) but also on the coupling between physics and dynamics of climate and weather prediction models.

Rios‐Berrios, R.↗

DOE OSTI2022

HI-STORM Overpack and MPC-32 Thermal-Hydraulic Model with MOOSE Framework

Nuclear power is a significant source of electricity in the United States, but the average age of nuclear power plants is around 40 years old. The safe management of the spent nuclear fuel (SNF) is a key aspect of the back-end of the nuclear fuel cycle. Spent fuel dry storage systems are becoming a popular and effective solution in this area, given the absence of a final disposal system. The spent fuel cask system (dry cask method) provides a feasible solution to maintain spent fuel for 60 years before final disposal. Dry cask storage has many characteristics that make it attractive. It fulfills the safety requirements of the Nuclear Regulatory Commission (NRC) while providing modularity and flexibility to contractors. The HI-STORM overpack and MPC-32 canister are the main parts of the HI-STORM 100 dry cask storage system. These components remove heat from the system using natural circulation, requiring no human intervention. This is the characteristic that provides passive heat removal and low maintenance features in dry cask storage systems. To develop a thermal model for a dry cask storage system, the physics behind the system should be defined clearly. There are two natural circulation loops in the system; circulation of helium cools down the nuclear assemblies in the MPC, while circulation of air cools down the walls of the MPC. This work aims to develop a thermal model of the MPC-32 canister and HI-STORM overpack using the Multiphysics Object-Oriented Simulation Environment (MOOSE). MOOSE is an open-source framework developed by Idaho National Laboratory (INL) for multiscale, multiphysics simulations. In this study, we will investigate and demonstrate the thermal-hydraulics modeling capabilities of the MOOSE framework, including natural circulation, heat transfer, and porous flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Time projection chamber for GADGET II

The established Gaseous Detector with Germanium Tagging (GADGET) detection system is used to measure weak, low-energy 𝛽-delayed proton decays. It consists of the Gaseous Proton Detector equipped with a MICROMEGAS (MM) readout to detect protons and other charged particles calorimetrically, surrounded by the Segmented Germanium Array (SeGA) for high-resolution detection of prompt 𝛾 rays. To upgrade GADGET's Proton Detector to operate as a compact time projection chamber (TPC) for the detection, three-dimensional imaging and identification of low-energy 𝛽-delayed single- and multiparticle emissions mainly of interest to astrophysical studies. A new high granularity MM board with 1024 pads has been designed, fabricated, installed, and tested. A high-density data acquisition system based on generic electronics for TPCs (GET) has been installed and optimized to record and process the gas avalanche signals collected on the readout pads. The TPC's performance has been tested using a 220 Rn 𝛼-particle source and cosmic-ray muons. In addition, decay events in the TPC have been simulated by adapting the attpcroot data analysis framework. Furthermore, a novel application of two-dimensional convolutional neural networks for GADGET II event classification is introduced. The optimization of data throughput is also addressed. The GADGET II TPC is capable of detecting and identifying 𝛼 particles as well as measuring their track direction, range, and energy. The extracted energy resolution of the GADGET II TPC using P10 gas is about 5.4% at 6.288 MeV ( 220 Rn 𝛼 events), computed using charge integration. Based on a systematic simulation study, we estimated the detection efficiency of the GADGET II TPC for protons and 𝛼 particles, respectively. It has also been demonstrated that the GADGET II TPC is capable of tracking minimum-ionizing particles (i.e., cosmic-ray muons). From these measurements, the electron drift velocity was measured under typical operating conditions. In addition to being one of the first generation of micropattern gaseous detectors (MPGDs) to utilize a resistive anode applied to low-energy nuclear physics, the GADGET II TPC will also be the first TPC surrounded by a high-efficiency array of high-purity germanium 𝛾-ray detectors. As a result, the TPC of GADGET II has been designed, fabricated, and tested and is ready for operation at the Facility for Rare Isotope Beams for radioactive-beam-line experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Ecosystems and Networks Integrated with Genes and Molecular Assemblies (ENIGMA) (Final Scientific/Technical Report, Subcontract Award No. 6953691)

The overarching objective the Baliga lab is to develop a mechanistic understanding of the field relevant metabolic processes underlying community partitioning of respiration pathways and physiological state shifts. In order to dissect how environmental perturbations influence microbial regulatory networks, we have constructed data-driven models that capture the dynamic interplay between abiotic and biotic factors during laboratory simulations. Specifically, we have developed multiple computational models (i.e., cMonkey2, EGRIN 2.0, and multiple genome-scale metabolic models) connecting environmental influence to genome-encoded regulatory programs. In addition, we have generated new tools such as Live Anaerobic Cell Sorting (LAnCS) and Fluidized Bed Reactors (FBRs) to investigate the dynamics of community partitioning of metabolic processes supporting sulfate or nitrate respiration, which are two critical activities observed from the field site at Oak Ridge National Labs (ORNL). This framework of connecting field phenomena with laboratory simulations and computational modeling has created the Environmental Simulations and Modelling (EnvSim) campaign within ENIGMA. The campaign, which spans multiple laboratories across ENIGMA and is led by the Baliga group, has established synthetic communities (SynComs) to discover, characterize, and dissect key microbial processes relevant to field observations, such as emissions of the greenhouse gas nitrous oxide (N 2 O). For example, the EnvSim campaign so far has deduced four potential mechanisms that may account for the N 2 O emissions at the FRC and are currently being investigated by labs across ENIGMA. Denitrification may be driven by complete denitrifiers, however, their NosZ enzymes, which catalyze the final step in denitrification by converting N 2 O to N 2 , may be sensitive at a lower pH. Additionally, the denitrification process could be partitioned among organisms and some may have pH-sensitive NosZ genes. Another factor contributing to variable N 2 O emissions relates to the metal co-factors involved in the denitrification pathway. For instance, excess Cu, Al, Mn, U, Ni, Co, Cu, and/or Cd may have inhibitory effects on multiple enzymatic steps during denitrification, while the enzymatic production of nitrite, the precursor of N 2 O, can be limited by the essential metal Mo. Lastly, abiotic production of N 2 O may occur as a result of chemodenitrification, in which metals like Fe, Mn, and some organic compounds can drive redox reactions that convert nitrogen cycle intermediates to N 2 O under the right conditions. Current investigations and analyses have generated multiple transcriptomic profiles aiding the refinement of new metabolic and gene regulatory network model for our denitrifying SynCom. In sum, the previous funding cycle has generated 33 peer-reviewed publications spanning predictive network biology to co-operativity of mutualistic interspecies interactions of evolved communities. Our integration of multiple omics datasets and the development of new network modelling tools has revealed important insights that can explain field observed phenomenon and generated new hypotheses to be investigated in the field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Geologic framework model for the Dry Alluvium Geology (DAG) experiment testbed, Yucca Flat, Nevada National Security Site

Geologic framework models (GFMs) provide a methodology for integrating geology into, and thus geologically informing, other modeling and simulation activities. GFMs provide a three-dimensional (3-D), geology-based digital framework for developing and parametrizing meshes and evaluating model and simulation results. This report describes a 3-D GFM constructed for the Source Physics Experiment Phase II Dry Alluvium Geology test series located in Yucca Flat at the Nevada National Security Site. The geology in the Yucca Flat region is complex and diverse, which creates challenges to modeling seismic wave propagation from SPE tests. The Yucca Flat GFM helps address these challenges by providing the 3-D distribution of relevant geologic features and physical properties necessary to more effectively model seismic wave propagation. The GFM includes 7 model layers and 48 faults that cut and offset the layers. An appendix is included that provides quantitative data on physical properties for each model layer.

58 GEOSCIENCES↗

NASA NTRS2018

Flight Demonstration of Novel Atmospheric Satellite Concept

The major focus of the Phase II effort described herein is to develop and demonstrate an aircraft capable of autonomously sailing (i.e., to cruise without propulsion or external assistance), and thereby prove that the dual-aircraft platform (DAP) atmospheric satellite concept is potentially viable. This sailing mode of flight was identified as the number-1 enabling technology required for the stratospheric DAP concept (also known as Stratosat) in the NIAC (NASA Innovative Advanced Concept) Phase I effort. No scientific demonstration of this technology has ever been done or documented to our knowledge. This report describes efforts to take a major step towards the sailing mode of flight capability using a single aircraft connected by cable to a moving ground vehicle which uses sufficient crosswind to cruise without propulsion while "pulling" the ground vehicle forward (i.e., without external assistance). The development of a prototype aircraft is described in terms of novel and key hardware and software elements. A specialized prototype aircraft is described, including a novel cable release mechanism, novel "lateron" control surfaces, and a highly-accurate onboard wind measurement system. Additionally, a novel means to safely connect the aircraft to the moving ground vehicle is described involving a fishing rod/reel and integrated load cell. All of these devices were designed and developed in-house and validated in flight testing. Software is developed to provide look-up tables that give the flight condition targets (i.e., 3-D position relative to ground vehicle, forward speed, aircraft orientation, etc.), based on current wind speed and direction. These tables are successfully validated in flight simulation and used onboard the aircraft. High fidelity analysis of the aircraft aerodynamics are described - required to produce accurate target sailing flight conditions. A novel wind tunnel measurement technique is developed to accurately assess the aerodynamics of the ultra-thin cable. A new specialized flight simulator is described which is utilized to develop and verify the flight software required onboard the aircraft, and to support training of pilots for flying the aircraft while tethered to a ground vehicle. The DAP flight simulator was developed within the Matlab-Simulink framework and included detailed treatment of aircraft/cable aerodynamics, cable dynamics, experimentally-derived propeller-motor thrust curves, actuator responsiveness, and realistic air turbulence. The specialized formation flight controller algorithm, developed using this flight simulator, and onboard the aircraft is described. Finally, a novel auto-tuning software is described and verified within the flight simulator that is shown to refine the sailing flight condition targets during flight using an optimization technique involving doublet maneuvers. Virtual flights using the auto-tuning software indicate that the prototype aircraft should be able to reach and hold sailing conditions despite moderate levels of turbulence provided there is sufficient mean wind available. An overview of the flight testing program is provided. Hundreds of short flights were conducted, primarily using a “dead” short runway at Deland Municipal Airport which permitted use of a moving ground vehicle. Additional flight tests at Space Florida’s Shuttle Landing Facility are also described. First year results from these tests in which the aircraft is controlled manually, demonstrated that excessive flight testing would be required for a pilot to learn to sail with visual cues. However, second year results from autonomous flight these tests included successful demonstration of the closed-loop autonomous formation flight capability (i.e., autonomously determine, reach, and hold the required 3-D location relative to the ground vehicle required for sailing). The next step of using the auto-tune software to autonomously refine the aircraft orientation targets to finally achieve sailing remains the primary goal of future work.

Flight↗

DOE OSTI2020

Real-World Distribution System Modeling Framework for Transmission-and-Distribution Cosimulation: Preprint

This paper presents a modeling methodology for realistic distribution system simulation and analysis. The methodology involves three major approaches: utility model conversion, feeder load modeling, and feeder model validation. The feeder models obtained from the utility are converted to the format that is more flexible for analysis and algorithm development. The load profiles down to each node are modeled in detail using advanced metering infrastructure data and supervisory control and data acquisition (SCADA) system measured load data. Then the distribution system models are validated by comparing the simulated feeder-head voltage results and SCADA measured voltage data. To better understand the bulk system operations and the interactions between transmission and distribution systems, the modeled realistic distribution systems are integrated into a transmission-and-distribution cosimulation framework to perform the system simulation from the bulk system down to each node in the distribution system.

load modeling↗

DOE OSTI2024

Assessment of Current MACCS Capabilities for Modeling Atmospheric Physical and Chemical Transformations

The physical and chemical transformation during atmospheric transport of radionuclides released into the environment has the possibility of impacting consequence modeling results. Accordingly, this report identifies physical and chemical transformations that may occur following release of chemically reactive radioactive species, how those transformations may affect modeling of consequences of release to the atmosphere and identifies current capabilities – in both MACCS and other state-of-practice atmospheric transport and dispersions models– to model those transformations. It was found that the inclusion of physical and chemical transformations is currently very limited in current state-of-practice codes for atmospheric dispersion of radionuclides. State-of-practice atmospheric dispersion codes appear to be typically limited to simulating either physical-chemical transformations or radioactive transformations, but not both. A state-of-practice atmospheric dispersion code capable of performing parallel physical, chemical, and radioactive transformation was not identified. A few atmospheric dispersion codes capable of modeling physical and chemical transport of specific species such as tritium or uranium hexafluoride were identified. Consequently, there is currently no information available that clearly suggests updates to the MACCS code are needed to bring it up to state-of-practice. However, investigations concluded that the MACCS computational framework can currently accommodate multiple physical/chemical forms in one simulation. Additionally, with some major assumptions, the computational framework in MACCS can accommodate parallel physical-chemical and radioactive transformations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2020

In-Situ Inference: Bringing Advanced Data Science Into Exascale Simulations

As simulations generate ever-increasing amounts of data, there are correspondingly richer opportunities for analysis and scientific discovery - discoveries that will be missed if most of the data must be discarded before it is analyzed. Because future exascale architectures will be increasingly storage-limited, it will not be possible to save the vast majority of simulation data for later analysis, requiring analysis to occur “in-situ” within the simulation. However, existing in-situ data analysis frameworks provide little or no support for one of the most sophisticated forms of data science: probabilistic statistical modeling or uncertainty quantification (UQ), and the accompanying challenge of inference - fitting those statistical models to massive simulation output. Our goal is to develop the fundamental statistical algorithms and computer science needed to perform statistical inference in-situ (in HPC simulations) to the full stream of data those simulations generate.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Assessment of Physical Security Modeling and Simulation in the Vulnerability Assessment Process

This report provides a comprehensive assessment of physical security modeling and simulation tools available for use in the vulnerability assessment (VA) process for nuclear facilities. It outlines the historical evolution of VA methodologies, emphasizing the transition from traditional layer-based approaches to a more holistic framework that integrates detection probabilities directly into combat simulations. The document details the critical components of the VA process, including the characterization of targets, threats, and protective measures, as well as the development of adversary scenarios that reflect both insider and outsider threats. It highlights the importance of performance assurance programs, emphasizing the need for continuous evaluation and testing of security systems to ensure their effectiveness against evolving threats. Additionally, the report discusses the significance of utilizing accredited modeling and simulation tools in accredited areas to accurately represent adversary actions and the corresponding responses of protective forces. By establishing a systematic approach to VA, this document aims to enhance the overall security posture of nuclear facilities, ensuring compliance with regulatory standards while effectively mitigating risks associated with potential adversarial actions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2026

Understanding carbon sourcing and transport originating from the helicon antenna surfaces during high-power helicon discharge in DIII-D Tokamak

The high-power helicon wave system in the DIII-D tokamak could potentially introduce new plasma–material interaction (PMI) challenges owing to rectified RF sheath potentials that develop near the antenna and surrounding plasma-facing components. We present the first application of the STRIPE (Simulated Transport of RF Impurity Production and Emission) framework to helicon-induced PMIs, extending previous STRIPE studies of ICRH antennas by incorporating net erosion, local re-deposition, and three-dimensional global impurity transport. The integrated workflow couples SOLPS, COMSOL, RustBCA, GITR, and GITRm to simulate carbon erosion, re-deposition, and impurity transport for two experimentally constrained DIII-D H-mode helicon operating scenarios with different antenna–plasma gaps, coupled RF powers, and edge plasma conditions. COMSOL predicts rectified RF sheath potentials of 1–5 kV localized near the lower portion of the antenna, where the magnetic field intersects the surface at grazing incidence. Carbon self-sputtering dominates the erosion source, whereas RF-accelerated D + ions contribute approximately 1% of the total gross erosion. The smaller-gap operating scenario exhibits substantially stronger gross erosion, enhanced local re-deposition (∼12%), and a larger confined carbon inventory owing to increased plasma accessibility and broader RF sheath coverage. Comparison with available DIII-D measurements shows no distinct change in the global carbon signal that correlates with the helicon RF pulse, consistent with the simulations indicating that the helicon-generated carbon source remains small compared with the existing background carbon inventory under the present graphite-wall operating conditions. These results demonstrate the capability of STRIPE to integrate RF sheath modeling, plasma transport, surface interaction physics, and three-dimensional impurity transport for the interpretation of helicon-induced PMIs. The study further identifies the principal sources of modeling uncertainty, including grazing-angle RF sheath physics, slow-wave resolution, plasma-background extrapolation, and trace-impurity transport assumptions, providing a framework for future validation and model development.

Kumar, Atul [Oak Ridge National Laboratory (ORNL),↗

DOE OSTI2021

Evaluation of Portable Acceleration Solutions for LArTPC Simulation Using Wire-Cell Toolkit

The Liquid Argon Time Projection Chamber (LArTPC) technology plays an essential role in many current and future neutrino experiments. Accurate and fast simulation is critical to developing efficient analysis algorithms and precise physics model projections. The speed of simulation becomes more important as Deep Learning algorithms are getting more widely used in LArTPC analysis and their training requires a large simulated dataset. Heterogeneous computing is an efficient way to delegate computationally intensive tasks to specialized hardware. However, as the landscape of compute accelerators quickly evolves, it becomes increasingly difficult to manually adapt the code to the latest hardware or software environments. A solution which is portable to multiple hardware architectures without substantially compromising performance would thus be very beneficial, especially for long-term projects such as the LArTPC simulations. In search of a portable, scalable and maintainable software solution for LArTPC simulations, we have started to explore high-level portable programming frameworks that support several hardware backends. In this paper, we present our experience porting the LArTPC simulation code in the Wire-Cell Toolkit to NVIDIA GPUs, first with the CUDA programming model and then with a portable library called Kokkos. Preliminary performance results on NVIDIA V100 GPUs and multi-core CPUs are presented, followed by a discussion of the factors affiecting the performance and plans for future improvements.

Yu, Haiwang↗