Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time parallelization”

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

The Pointing Self-calibration Algorithm for Aperture Synthesis Radio Telescopes

This paper is concerned with algorithms for calibration of direction-dependent effects (DDE) in aperture synthesis radio telescopes (ASRT). After correction of direction-independent effects (DIE) using self-calibration, imaging performance can be limited by the imprecise knowledge of the forward gain of the elements in the array. In general, the forward gain pattern is directionally dependent and varies with time due to a number of reasons. Some factors, such as rotation of the primary beam with Parallactic Angle for Azimuth–Elevation mount antennas are known a priori. Some, such as antenna pointing errors and structural deformation/projection effects for aperture-array elements cannot be measured a priori. Thus, in addition to algorithms to correct for DD effects known a priori, algorithms to solve for DD gains are required for high dynamic range imaging. Here, we discuss a mathematical framework for antenna-based DDE calibration algorithms and show that this framework leads to computationally efficient optimal algorithms that scale well in a parallel computing environment. As an example of an antenna-based DD calibration algorithm, we demonstrate the Pointing SelfCal (PSC) algorithm to solve for the antenna pointing errors. Our analysis show that the sensitivity of modern ASRT is sufficient to solve for antenna pointing errors and other DD effects. We also discuss the use of the PSC algorithm in real-time calibration systems and extensions for antenna Shape SelfCal algorithm for real-time tracking and corrections for pointing offsets and changes in antenna shape.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Machine Learning for Big Data: A Study to Understand Limits at Scale

This report aims to empirically understand the limits of machine learning when applied to Big Data. We observe that recent innovations in being able to collect, access, organize, integrate, and query massive amounts of data from a wide variety of data sources have brought statistical data mining and machine learning under more scrutiny, evaluation and application for gleaning insights from the data than ever before. Much is expected from algorithms without understanding their limitations at scale while dealing with massive datasets. In that context, we pose and address the following questions How does a machine learning algorithm perform on measures such as accuracy and execution time with increasing sample size and feature dimensionality? Does training with more samples guarantee better accuracy? How many features to compute for a given problem? Do more features guarantee better accuracy? Do efforts to derive and calculate more features and train on larger samples worth the effort? As problems become more complex and traditional binary classification algorithms are replaced with multi-task, multi-class categorization algorithms do parallel learners perform better? What happens to the accuracy of the learning algorithm when trained to categorize multiple classes within the same feature space? Towards finding answers to these questions, we describe the design of an empirical study and present the results. We conclude with the following observations (i) accuracy of the learning algorithm increases with increasing sample size but saturates at a point, beyond which more samples do not contribute to better accuracy/learning, (ii) the richness of the feature space dictates performance - both accuracy and training time, (iii) increased dimensionality often reflected in better performance (higher accuracy in spite of longer training times) but the improvements are not commensurate the efforts for feature computation and training and (iv) accuracy of the learning algorithms drop significantly with multi-class learners training on the same feature matrix and (v) learning algorithms perform well when categories in labeled data are independent (i.e., no relationship or hierarchy exists among categories).

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2015

SU-E-T-143: Automatic PBS Treatment Planning for Prostate Cancer

Purpose: This study presents a novel idea of using automated planning technique for parallel-opposed prostate plans with pencil beam scanning proton therapy. Methods: A random 5 clinical prostate patients treated with parallel-opposed modulated-scanning proton beam in our proton center were selected in this study. All of the 5 cases were re-planned with in-house developed python scripts integrated into RayStation 4.0 clinical version. This automatic tool was designed to perform patient DICOM data import, external structure contours, CT density assignment, prescription to PTV assignment, beam placement, optimization, and dose calculation. Results: The mean time to generate a complete treatment plan was 1 minute per plan. For the automatically generated plans, 5 out of 5 plans (100%) were deemed clinically acceptable by one experienced radiation oncologist. In addition, the automated plans were overall dosimetrically equivalent to the clinical plans when scored for target coverage, rectum and bladder dose. Mann-Whitney U tests did not show significant dosimetric differences for normal tissue structures between the clinical and automated plans. Conclusion: We have developed a robust and automated method for fully inversed planned paralleled-opposed prostate proton treatment planning. This efficient method can be readily integrated into clinical practice. The tool produces clinically acceptable plans using only manually contoured anatomic structures and simulation CT data set. We anticipate that this tool will not only improve patient access to high-quality proton treatment by simplifying the planning process and will also reduce the effort and cost of incorporating more advanced planning into clinical practice in the near future.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

VeloxChem: A Python-driven density-functional theory program for spectroscopy simulations in high-performance computing environments

An open-source program named VeloxChem has been developed for the calculation of electronic real and complex linear response functions at the levels of Hartree–Fock and Kohn–Sham density functional theories. With an object-oriented program structure written in a Python/C++ layered fashion, VeloxChem enables time-efficient prototyping of novel scientific approaches without sacrificing computational efficiency, so that molecular systems involving up to and beyond 500 second-row atoms (or some 10,000 contracted and in part diffuse Gaussian basis functions) can be routinely addressed. In addition, VeloxChem is equipped with a polarizable embedding scheme for the treatment of the classical electrostatic interactions with an environment that in turn is modeled by atomic site charges and polarizabilities. The underlying hybrid message passing interface (MPI)/open multiprocessing (OpenMP) parallelization scheme makes VeloxChem suitable for execution in high-performance computing cluster environments, showing even slightly beyond linear scaling for the Fock matrix construction with use of up to 16,384 central processing unit (CPU) cores. An efficient—with respect to convergence rate and overall computational cost—multifrequency/gradient complex linear response equation solver enables calculations not only of conventional spectra, such as visible/ultraviolet/X-ray electronic absorption and circular dichroism spectra, but also time-resolved linear response signals as due to ultra-short weak laser pulses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Analysis the configuration of earthing system based on high-low and low-high soil structure

Each TNB transmission tower requires a tower footing resistance (TFR) with a lower grounding resistance value that depends on the transmission line voltage. For 132kV and 275kV tower, the TFR must less than 10Ω and 500kV tower must less than 5Ω. The TFR is changeable with variable factors such as soil resistivity. Low TFR provides essential protection to the fault such as lightning strike that may occur at any time. The fault current flow to the lowest resistance path and easily disperses to earth. Back flashover voltage across the insulator of transmission lines may occur when the TFR is high. The TFR is influenced by soil resistivity. There are three parameters affecting the soil resistivity; moisture content, salt content and temperature of the soil. High moisture content in soil will reduce the soil resistivity and resultant low TFR. Small scale moisture control by using Micro Reservoir (MR) irrigation with semi-permeable membranes have the power to offer the stable moisture in soil. By using osmosis concept, it is the process of net movement of water molecules from high potential water to lower potential water though a semi permeable membrane. The MR can withstand for 3 to 5 days without continuous water supply. The MR installed in the centre of the tower that contains a multiple parallel of electrode rods. The concentrated of electrode rods grounding configuration with a combination of MR will improve the TFR even at multilayer soil. As a result, MR gives a little improvement to TFR. The MR in area of concentrated electrode rod configuration to ensure the soil always wet and moist at all times. The changes in soil affect the tower-footing-resistance. The tower-footing-resistance measurement at afternoon is higher than at evening because of the temperature and moisture content in soil is change due to sun radiation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Enabling Modular Autonomous Feedback‐Loops in Materials Science through Hierarchical Experimental Laboratory Automation and Orchestration

Abstract Materials acceleration platforms (MAPs) operate on the paradigm of integrating combinatorial synthesis, high‐throughput characterization, automatic analysis, and machine learning. Within a MAP, one or multiple autonomous feedback loops may aim to optimize materials for certain functional properties or to generate new insights. The scope of a given experiment campaign is defined by the range of experiment and analysis actions that are integrated into the experiment framework. Herein, the authors present a method for integrating many actions within a hierarchical experimental laboratory automation and orchestration (HELAO) framework. They demonstrate the capability of orchestrating distributed research instruments that can incorporate data from experiments, simulations, and databases. HELAO interfaces laboratory hardware and software distributed across several computers and operating systems for executing experiments, data analysis, provenance tracking, and autonomous planning. Parallelization is an effective approach for accelerating knowledge generation provided that multiple instruments can be effectively coordinated, which the authors demonstrate with parallel electrochemistry experiments orchestrated by HELAO. Efficient implementation of autonomous research strategies requires device sharing, asynchronous multithreading, and full integration of data management in experimental orchestration, which to the best of the authors’ knowledge, is demonstrated for the first time herein.

36 MATERIALS SCIENCE↗
DOE OSTI2017

The influence of seawater carbonate chemistry, mineralogy, and diagenesis on calcium isotope variations in Lower-Middle Triassic carbonate rocks

The geological calcium cycle is linked to the geological carbon cycle through the weathering and burial of carbonate rocks. As a result, calcium (Ca) isotope ratios ( 44 Ca/ 40 Ca, expressed as δ 44/40 Ca) can help to constrain ancient carbon cycle dynamics if Ca cycle behavior can be reconstructed. But, the δ 44/40 Ca of carbonate rocks is influenced not only by the δ 44/40 Ca of seawater but also by diagenetic processes and fractionation associated with carbonate precipitation. In this study, we investigate the dominant controls on carbonate δ 44/40 Ca in Upper Permian to Middle Triassic limestones (ca. 253 to 244 Ma) from south China and Turkey. This time interval is ideal for assessing controls on Ca isotope ratios in carbonate rocks because fluctuations in seawater δ 44/40 Ca may be expected based on several large carbon isotope (δ 13 C) excursions ranging from -2 to + 8‰. Parallel negative δ 13 C and δ 44/40 Ca excursions were previously identified across the end-Permian extinction horizon. Here, we find a second negative excursion in δ 44/40 Ca of ~ 0.2‰ within Lower Triassic strata in both south China and Turkey; however, this excursion is not synchronous between regions and thus cannot be interpreted to reflect secular change in the δ 44/40 Ca of global seawater. Additionally, δ 44/40 Ca values from Turkey are consistently 0.3‰ lower than contemporaneous samples from south China, providing further support for local or regional influences. By measuring δ 44/40 Ca and Sr concentrations ([Sr]) in two stratigraphic sections located at opposite margins of the Paleo-Tethys Ocean, we can determine whether the data represent global conditions (e.g., secular variations in the δ 44/40 Ca of seawater) versus local controls (e.g., original mineralogy or diagenetic alteration). The [Sr] and δ 44/40 Ca data from this study are best described statistically by a log-linear correlation that also exists in many previously published datasets of various geological ages. By using a model of early marine diagenetic water-rock interaction, we illustrate that this general correlation can be explained by the chemical evolution of bulk carbonate sediment samples with different initial mineralogical compositions that subsequently underwent recrystallization. Although early diagenetic resetting and carbonate mineralogy strongly influence the carbonate δ 44/40 Ca values, the relationship between [Sr] and δ 44/40 Ca holds potential for reconstructing first-order secular changes in seawater δ 44/40 Ca composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Numerical and experimental analysis of mechanically induced failure in electric vehicle battery modules

Mitigating thermal runaway and cell-to-cell propagation is essential for improving the safety of electric and hybrid vehicles. Enhancing digital twin capabilities to predict battery mechanical abuse is particularly critical for automotive and aerospace applications, where crashworthiness is a key concern. Understanding failure conditions and propagation in battery modules during mechanical abuse is complex due to interactions between structural deformation, heat transfer, electrochemical processes, exothermic reactions and mechanical fracture. While prior studies have focused on modeling cell-level behavior, extending these models to module or pack level is necessary for a system level understating of electric vehicle safety. This study develops coupled large deformation finite element models that simultaneously solve for electrochemistry, material failure, internal short circuit and thermal runaway propagation. The models account for mechanical and thermal interactions between lithium-ion cells and other battery components while the contact interfaces are evolving with time. Model-predicted voltage, temperature and force responses are compared with experimental data for validation. The results demonstrate that the approach captures key failure mechanisms, including thermal propagation through heat transfer, electrical propagation from short circuits in parallel-connected cells, and mechanical propagation via penetration and crack formation. These findings show that computational models are valuable tools for understanding battery module failure and providing insight that can reduce the need for extensive experimental testing.

25 ENERGY STORAGE↗
DOE OSTI2017

Structure and decomposition of the silver formate Ag(HCO{sub 2})

Crystal structure of the silver formate Ag(HCO{sub 2}) has been determined (orthorhombic, sp.gr. Pccn, a=7.1199(5), b=10.3737(4), c=6.4701(3)Å, V=477.88(4) Å{sup 3}, Z=8). The structure contains isolated formate ions and the pairs Ag{sub 2}{sup 2+} which form the layers in (001) planes (the shortest Ag–Ag distances is 2.919 in the pair and 3.421 and 3.716 Å between the nearest Ag atoms of adjacent pairs). Silver formate is unstable compound which decompose spontaneously vs time. Decomposition was studied using Rietveld analysis of the powder diffraction patterns. It was concluded that the diffusion of Ag atoms leads to the formation of plate-like metal particles as nuclei in the (100) planes which settle parallel to (001) planes of the silver formate matrix. - Highlights: • Silver formate Ag(HCO{sub 2}) was synthesized and characterized. • Layered packing of Ag-Ag pairs in the structure was found. • Decomposition of Ag(HCO{sub 2}) and formation of metal phase were studied. • Rietveld-refined micro-structural characteristics during decomposition reveal the space relationship between the matrix structure and forming Ag phase REPLACE with: Space relationship between the matrix structure and forming Ag phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Dynamics of the vaginal wall dose in HDR interstitial brachytherapy for gynecological cancer: Systematic analysis of phantom vs patient case

This study aimed to investigate the dynamics of the vaginal wall dose for interstitial brachytherapy (ISBT). A patient undergoing ISBT was selected as the patient case. The phantom case was generated to simulate the patient case in all regards with the exception of parallel needle positions. The vaginal wall was contoured as a 0.5-cm expansion around the vaginal surface of the obturator. The prescribed ISBT dose was 20 Gy in 4 fractions. Six treatment plans were generated by modifying relative dwell times and needle positions (DTNP). The volume of the vaginal wall receiving > 150% of prescription dose (V{sub > 150%}) and D{sub 2cc} of the vaginal wall were compared among plans. The V{sub > 150%} was much larger in the patient case (49.3%) due to unparallel needles compared with the phantom case (14.3%) without modification (plan 1). Among the 6 plans, reduced dwell time (plan 3) and no dwell time (plans 5 and 6) on the vaginal surface needles had the lowest vaginal wall doses with the use of a central obturator needle in both cases. In comparison of patient case plans 1, 3, 5, and 6, V{sub 150%} was 49.2%, 19.0%, 21.3%, and 28.7%, respectively, and D{sub 2cc} was 41.15 Gy, 33.10 Gy, 36.51 Gy, and 34.37 Gy, respectively, which was limited around each loaded needle. Modification of DTNP is able to reduce the vaginal wall volume exceeding 150% of the prescription dose in the patient case. Understanding these dynamics of the vaginal wall dose will improve dose optimization of ISBT and may reduce vaginal morbidities.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Oceanic Heat Delivery to the Antarctic Continental Shelf: Large-Scale, Low-Frequency Variability

According to the article, onshore penetration of oceanic water across the Antarctic continental slope (ACS) plays a major role in global sea level rise by delivering heat to the Antarctic marginal seas, thus contributing to the basal melting of ice shelves. Here the time-mean (Φ mean ) and eddy (Φ eddy ) components of the heat transport (Φ) across the 1,000-m isobath along the entire ACS are investigated using a 0.1° global coupled ocean/sea ice simulation based on the Los Alamos Parallel Ocean Program (POP) and sea ice (Community Ice CodE) models. In this work, comparison with in situ hydrography shows that the model successfully represents the basic water mass structure, with a warm bias in the Circumpolar Deep Water layer. Segments of on-shelf Φ, with lengths of O(100–1,000 km), are found along the ACS. The circumpolar integral of the annually averaged Φ is O(20 TW), with Φ eddy always on-shelf, while Φ mean fluctuates between on-shelf and off-shelf. Stirring along isoneutral surfaces is often the dominant process by which eddies transport heat across the ACS, but advection of heat by both mean flow-topography interactions and eddies can also be significant depending on the along- and across-slope location. The seasonal and interannual variability of the circumpolarly integrated Φ mean is controlled by convergence of Ekman transport within the ACS. Prominent warming features at the bottom of the continental shelf (consistent with observed temperature trends) are found both during high-Southern Annular Mode and high-Niño 3.4 periods, suggesting that climate modes can modulate the heat transfer from the Southern Ocean to the ACS across the entire Antarctic margin.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Rapid discovery and evolution of nanosensors containing fluorogenic amino acids

Binding-activated optical sensors are powerful tools for imaging, diagnostics, and biomolecular sensing. However, biosensor discovery is slow and requires tedious steps in rational design, screening, and characterization. Here we report on a platform that streamlines biosensor discovery and unlocks directed nanosensor evolution through genetically encodable fluorogenic amino acids (FgAAs). Building on the classical knowledge-based semisynthetic approach, we engineer ~15 kDa nanosensors that recognize specific proteins, peptides, and small molecules with up to 100-fold fluorescence increases and subsecond kinetics, allowing real-time and wash-free target sensing and live-cell bioimaging. An optimized genetic code expansion chemistry with FgAAs further enables rapid (~3 h) ribosomal nanosensor discovery via the cell-free translation of hundreds of candidates in parallel and directed nanosensor evolution with improved variant-specific sensitivities (up to ~250-fold) for SARS-CoV-2 antigens. Altogether, this platform could accelerate the discovery of fluorogenic nanosensors and pave the way to modify proteins with other non-standard functionalities for diverse applications.

Biosensors↗
DOE OSTI2015

Todorokite-type manganese oxide nanowires as an intercalation cathode for Li-ion and Na-ion batteries

Extended hydrothermal treatment at an elevated temperature of 220 °C allowed high yield synthesis of manganese oxide nanowires with a todorokite crystal structure suitable for ions intercalation. The flexible, high aspect ratio nanowires are 50–100 nm in diameter and up to several microns long, with 3 × 3 structural tunnels running parallel to the nanowire longitudinal axis. Moreover, the tunnels are occupied by magnesium ions and water molecules, with the chemical composition found to be Mg 0.2 MnO 2 ·0.5H 2 O. The todorokite nanowires were, for the first time, electrochemically tested in both Li-ion and Na-ion cells. A first discharge capacity of 158 mA h g -1 was achieved in a Na-ion system, which was found to be greater than the first discharge capacity in a Li-ion system (133 mA h g -1 ). In spite of the large structural tunnel dimensions, todorokite showed a significant first cycle capacity loss in a Na-ion battery. After 20 cycles, the capacity was found to stabilize around 50 mA h g -1 and remained at this level for 100 cycles. In a Li-ion system, todorokite nanowires showed significantly better capacity retention with 78% of its initial capacity remaining after 100 cycles. Rate capability tests also showed superior performance of todorokite nanowires in Li-ion cells compared to Na-ion cells at higher current rates. Finally, these results highlight the difference in electrochemical cycling behavior of Li-ion and Na-ion batteries for a host material with spacious 3 × 3 tunnels tailored for large Na + ion intercalation.

25 ENERGY STORAGE↗
DOE OSTI2025

Quadrupolar density structures in driven magnetic reconnection experiments with a guide field

Magnetic reconnection is a ubiquitous process in plasma physics, driving rapid and energetic events such as coronal mass ejections. Reconnection between magnetic fields with arbitrary shear can be decomposed into an anti-parallel reconnecting component and a non-reconnecting guide-field component, which is parallel to the reconnecting electric field. This guide field modifies the structure of the reconnection layer and the reconnection rate. We present results from experiments on the MAIZE pulsed-power generator (500 kA peak current, 200 ns rise time), which use two exploding wire arrays, tilted in opposite directions, to embed a guide field in the plasma flows with a relative strength b≡B g /B rec =0, 0.4, or 1. The reconnection layers in these experiments have widths that are less than the ion skin depth, d i =c/ω pi , indicating the importance of the Hall term, which generates a distinctive quadrupolar magnetic field structure along the separatrices of the reconnection layer. Using laser imaging interferometry, we observe quadrupolar structures in the line-integrated electron density, consistent with the interaction of the embedded guide field with the quadrupolar Hall field. Our measurements extend over much larger length scales (40d i ) at higher β (∼1) than previous experiments, providing an insight into the global structure of the reconnection layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

High-resolution hybrid simulations of kinetic plasma turbulence at proton scales

We investigate properties of plasma turbulence from magnetohydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations with energy equipartition and vanishing correlation. Once the turbulence is fully developed, we observe an MHD inertial range, where the spectra of the perpendicular magnetic field and the perpendicular proton bulk velocity fluctuations exhibit power-law scaling with spectral indices of −5/3 and −3/2, respectively. This behavior is extended over a full decade in wavevectors and is very stable in time. A transition is observed around proton scales. At sub-ion scales, both spectra steepen, with the former still following a power law with a spectral index of ∼−3. A −2.8 slope is observed in the density and parallel magnetic fluctuations, highlighting the presence of compressive effects at kinetic scales. The spectrum of the perpendicular electric fluctuations follows that of the proton bulk velocity at MHD scales, and flattens at small scales. All these features, which we carefully tested against variations of many parameters, are in good agreement with solar wind observations. The turbulent cascade leads to on overall proton energization with similar heating rates in the parallel and perpendicular directions. While the parallel proton heating is found to be independent on the resistivity, the number of particles per cell, and the resolution employed, the perpendicular proton temperature strongly depends on these parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Large-scale atomistic calculations of clusters in intense x-ray pulses

Here, we present the methodology of our recently developed Monte-Carlo/ Molecular-Dynamics method for studying the fundamental ultrafast dynamics induced by high-fluence, high-intensity x-ray free electron laser (XFEL) pulses in clusters. The quantum nature of the initiating ionization process is accounted for by a Monte Carlo method to calculate probabilities of electronic transitions, including photo absorption, inner-shell relaxation, photon scattering, electron collision and recombination dynamics, and thus track the transient electronic configurations explicitly. The freed electrons and ions are followed by classical particle trajectories using a molecular dynamics algorithm. These calculations reveal the surprising role of electron-ion recombination processes that lead to the development of nonuniform spatial charge density profiles in x-ray excited clusters over femtosecond timescales. In the high-intensity limit, it is important to include the recombination dynamics in the calculated scattering response even for a 2- fs pulse. We also demonstrate that our numerical codes and algorithms can make e!cient use of the computational power of massively parallel supercomputers to investigate the intense-field dynamics in systems with increasing complexity and size at the ultrafast timescale and in non-linear x-ray interaction regimes. In particular, picosecond trajectories of XFEL clusters with attosecond time resolution containing millions of particles can be e!ciently computed on upwards of 262,144 processes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Doppler coherence imaging of divertor and SOL flows in ASDEX upgrade and Wendelstein 7-X

Measurement of ion dynamics is of great importance to the study of plasma exhaust and impurity particle flows in the scrape-off-layer (SOL) and divertor in magnetic confinement experiments. Doppler coherence imaging spectroscopy (CIS) is a passive optical diagnostic that produces 2D images of line-integrated ion flow velocity. Doppler CIS flow measurements of neutral deuterium and impurities were conducted for the first time in the divertor of the medium-sized tokamak experiment ASDEX Upgrade. A detailed sightline, emission and magnetic field analysis was undertaken to identify location and direction of flows. Under the assumption of toroidal axisymmetry, they revealed mainly parallel impurity flows in the proximity of the X-point and target plates on the order of 20 km s -1 . Two Doppler CIS systems are currently set up for the optimized stellarator Wendelstein 7-X (W7-X). With this being said, the main task of this diagnostic will be the measurement of impurity ion flows in one of the island divertors of W7-X. EMC3-EIRENE simulations have been carried out to estimate the flow behaviour in the W7-X SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A hybrid 3D/2D field response calculation for liquid argon detectors with PCB based anode plane

Liquid Argon Time Projection Chamber (LArTPC) technology is commonly utilized in neutrino detector designs. It enables detailed reconstruction of neutrino events with high spatial precision and low energy threshold. Its field response (FR) model describes the time-dependent electric currents induced in the anode-plane electrodes when ionization electrons drift nearby. An accurate and precise FR is a crucial input to LArTPC detector simulations and charge reconstruction. Established LArTPC designs have been based on parallel wire planes. It allows accurate and computationally economic two-dimensional (2D) FR models utilizing the translational symmetry along the direction of the wires. Recently, novel LArTPC designs utilize electrodes formed on printed circuit board (PCB) in the shape of strips with through holes. The translational symmetry is no longer a good approximation near the electrodes and a new FR calculation that employs regions with three dimensions (3D) has been developed. Extending the 2D models to 3D would be computationally expensive. Fortuitously, the nature of strips with through holes allows for a computationally economic approach based on the finite-difference method (FDM). In this paper, we present a new software package pochoir that calculates LArTPC field response for these new strip-based anode designs. This package combines 3D calculations in the volume near the electrodes with 2D far-field solutions to achieve fast and precise field response computation. We apply the resulting FR to simulate and reconstruct samples of cosmic-ray muons and 39 Ar decays from a Vertical Drift (VD) detector prototype operated at CERN. We find the difference between real and simulated data within 5%. Current state-of-the-art LArTPC software requires a 2D FR which we provide by averaging over one dimension and estimate that variations lost in this average are smaller than 7%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗