Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel processing”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Homomolecular Triplet–Triplet Annihilation in Metalloporphyrin Photosensitizers

Metalloporphyrins are ubiquitous in their applications as triplet photosensitizers, particularly for promoting sensitized photochemical upconversion processes. In this study, bimolecular excited state triplet–triplet quenching kinetics, termed homomolecular triplet–triplet annihilation (HTTA), exhibited by the traditional triplet photosensitizers–zinc(II) tetraphenylporphyrin (ZnTPP), palladium(II) octaethylporphyrin (PdOEP), platinum(II) octaethylporphyrin (PtOEP), and platinum(II) tetraphenyltetrabenzoporphyrin (PtTPBP)–were revealed using conventional transient absorption spectroscopy. Nickel(II) tetraphenylporphyrin was used as a control sample as it is known to be rapidly quenched intramolecularly through ligand-field state deactivation and, therefore, cannot result in triplet–triplet annihilation (TTA). The single wavelength transients associated with the metalloporphyrin triplet excited state decay–measured as a function of incident laser pulse energy in toluene–were well modeled using parallel first- and second-order kinetics, consistent with HTTA being operable. The combined transient kinetic data enabled the determination of the first-order rate constants (k T ) for excited triplet decay in ZnTPP (4.0 × 10 3 s –1 ), PdOEP (3.6 × 10 3 s –1 ), PtOEP (1.2 × 10 4 s –1 ), and PtTPBP (2.1 × 10 4 s –1 ) as well as the second-order rate constant (k TT ) for HTTA in ZnTPP (5.5 × 10 9 M –1 s –1 ), PdOEP (1.1 × 10 10 M –1 s –1 ), PtOEP (7.1 × 10 9 M –1 s –1 ), and PtTPBP (1.6 × 10 10 M –1 s –1 ). In most instances, triplet excited state extinction coefficients are either reported for the first time or have been revised using ultrafast transient absorption spectroscopy and singlet depletion: ZnTPP (78,000 M –1 cm –1 ) at 470 nm, PdOEP (67,000 M –1 cm –1 ) at 430 nm, PtOEP (51,000 M –1 cm –1 ) at 418 nm, and PtTPBP (100,000 M –1 cm –1 ) at 460 nm. Furthermore, the combined experimental results establish competitive time scales for homo- and heteromolecular TTA rate constants, implying the significance of considering HTTA processes in future research endeavors harnessing TTA photochemistry using common metalloporphyrin photosensitizers.

14 SOLAR ENERGY↗
DOE OSTI2021

Scalable balanced training of conditional generative adversarial neural networks on image data

Here, we propose a distributed approach to train deep convolutional generative adversarial neural network (DC-CGANs) models. Our method reduces the imbalance between generator and discriminator by partitioning the training data according to data labels, and enhances scalability by performing a parallel training where multiple generators are concurrently trained, each one of them focusing on a single data label. Performance is assessed in terms of inception score, Fréchet inception distance, and image quality on MNIST, CIFAR10, CIFAR100, and ImageNet1k datasets, showing a significant improvement in comparison to state-of-the-art techniques to training DC-CGANs. Weak scaling is attained on all the four datasets using up to 1000 processes and 2000 NVIDIA V100 GPUs on the OLCF supercomputer Summit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Design for Recycling Principles Applicable to Selected Clean Energy Technologies: Crystalline-Silicon Photovoltaic Modules, Electric Vehicle Batteries, and Wind Turbine Blades

Abstract The global growth of clean energy technology deployment will be followed by parallel growth in end-of-life (EOL) products, bringing both challenges and opportunities. Cumulatively, by 2050, estimates project 78 million tonnes of raw materials embodied in the mass of EOL photovoltaic (PV) modules, 12 billion tonnes of wind turbine blades, and by 2030, 11 million tonnes of lithium-ion batteries. Owing partly to concern that the projected growth of these technologies could become constrained by raw material availability, processes for recycling them at EOL continue to be developed. However, none of these technologies are typically designed with recycling in mind, and all of them present challenges to efficient recycling. This article synthesizes and extends design for recycling (DfR) principles based on a review of published industrial and academic best practices as well as consultation with experts in the field. Specific principles developed herein apply to crystalline-silicon PV modules, batteries like those used in electric vehicles, and wind turbine blades, while a set of broader principles applies to all three of these technologies and potentially others. These principles are meant to be useful for stakeholders—such as research and development managers, analysts, and policymakers—in informing and promoting decisions that facilitate DfR and, ultimately, increase recycling rates as a way to enhance the circularity of the clean energy economy. The article also discusses some commercial implications of DfR. Graphical Abstract

36 MATERIALS SCIENCE↗
DOE OSTI2021

Caffeine v0.1.0

Caffeine is the CoArray Fortran Framework of Efficient Interfaces to Network Environments. Caffeine aims to produce a parallel runtime library that will support Fortran compilers with a programming-model-agnostic application binary interface (ABI) to various lower-level communication libraries. The current version of Caffeine uses the GASNet-EX networking middleware, also developed at Berkeley Lab. On many combinations of applications and platforms, GASNet-EX outperforms the widely used Message Passing Interface (MPI). Through GASNet-EX's support for communicating between graphics processing units (GPU), GASNet-EX has features that specifically target the emerging, leading-edge exascale computing platforms.

Rouson, Damian↗
DOE OSTI2024

Influence of Antarctic and Greenland Continental Shelf Circulation on High‐Latitude Oceans in E3SM

The science objectives of this project are to simulate and understand the impacts of both deep-basin warm-water intrusions and land-ice melt on the continental shelf circulations and sea-ice distributions around the margins of Greenland and Antarctica. As well, the role of subsurface ocean heat from the Atlantic on declining sea-ice cover in the Arctic is explored. Mesoscale processes and fine bathymetry are implicated in cross-shelf property transports around both Greenland and Antarctica. Therefore, we configured and ran an atmospheric reanalysis-forced global ocean/sea-ice simulation on a grid that reduces from 8 km at the Equator to 2 km at the poles (UH8to2) with 60 vertical levels. It was produced using the Energy Exascale Earth System Model ‘‘HiLAT’’ code (E3SMv0-HiLAT) that uses the Parallel Ocean Program (POP) and CICE5 as its ocean and sea-ice components, respectively. Two main UH8to2 simulations were carried out: one for 1975-2009 and the other for July 2016-2020 after it was initialized from a 1/25° data-assimilative ocean/sea-ice prediction system ocean/sea-ice state. The UH8to2 is not coupled to an active land-ice model. Rather, land-ice melt is represented by observationally informed freshwater fluxes (FWFs). Short (multi-year) UH8to2 simulations were conducted to understand sensitivities when Greenland ice sheet (GrIS) melt is released only at the ocean surface or when it is distributed over the upper water column in accordance with fjord melt plume behavior; these cases were compared with a no GrIS melt case. West Greenland continental shelf currents were fastest in the vertical distribution case and an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy was found relative to the surface release case. Further, salinity is lower and meltwater volume greater in the eastern Labrador Sea in the vertical distribution case. For the Arctic, the veracity of the UH8to2 was evaluated for 2017-2020 using available observations. Simulated seasonal sea-ice thickness and concentration are realistic, but the ice is unrealistically thin in the central and eastern Arctic in the fall. Comparisons of vertical sections of ocean temperature, salinity, and buoyancy collected from Ice-Tethered Profilers (ITPs) in the eastern Arctic in the fall and winter of 2019/2020 and co-located/concurrent UH8to2 fields show the stratification over the top 100 m of the water column is too low in the model, the simulated mixed layer too deep, and the simulated subsurface Atlantic Water (AW) too warm; these biases may contribute to the sea-ice biases. A model intercomparison study using the UH8to2 and a forced 1/25° regional Arctic ocean/sea-ice (uses the HYbrid Coordinate Ocean Model and CICE5) simulation further investigates the relationship between AW and sea-ice in the eastern Arctic. The models show a mesoscale-rich pulse of Atlantic Water extending into the eastern basin that reaches maximum intensity in late winter of 2018, after which it decreases in strength. Concurrent and co-located sea-ice melt or the inhibition of sea-ice growth is seen and is attributed to halocline mesoscale eddies doming into the mixed layer with convection bringing this heat into the vicinity of the sea-ice.

58 GEOSCIENCES↗
DOE OSTI2025

ACHILLES-GENIE Interface for Neutrino Simulations, ADRIANO2 Tile Prototype for High-Granularity Dual-Readout Calorimetry

ICARUS (Imaging Cosmic And Rare Underground Signals) is a liquid argon time projection chamber (LArTPC) detector that pursues the sterile neutrino, which relies on accurate simulations of neutrino-argon interactions. REDTOP (Rare Eta Decays To Observe new Physics) is a proposed low-energy, high-intensity meson factory designed to explore rare $\eta$/$\eta'$ meson decays and probe physics beyond the Standard Model. As a next-generation experiment, this requires both accurate simulations and innovative detector technologies. This project contributes to both ICARUS, from a simulation perspective, and REDTOP, from both a simulation and detection perspective, through the event generation of lepton-nucleon interactions and the physical enhancement of the calorimeter technology within the REDTOP detector. We developed an interface between ACHILLES (A CHIcago Land Lepton Event Simulator), a theory-driven lepton-level event generator, and GENIE, a robust event generator framework used for neutrino physics. By incorporating the precise theoretical cross-section calculations of ACHILLES into the experimental realism of GENIE, the interface allows for improved accuracy of neutrino-nucleon simulations, which can be adapted for the proton beam specifications of the REDTOP meson factory as well as for the ICARUS experiment. In parallel, we developed an improved prototype for the ADRIANO2 (A Dual Readout Integrally Active Non-segmented Option) dual-readout calorimeter tiles for the REDTOP detector. To improve the efficiency of the lead-glass tiles trapping Cherenkov light for energy reconstruction and particle identification, we optimized the application of a highly reflective coating. Through viscosity and thickness control, masking, and a custom spray technique, we refined the coating process to reduce surface defects and improve light yield. Together, these efforts strengthen the ICARUS neutrino program and REDTOP's capability of detecting rare decay events.

Visser, Erin [Michigan State U.]↗
DOE OSTI2026

Quality Control of Silicon Sensor Modules for Particle Detectors

The High-Luminosity Large Hadron Collider (HL-LHC) will produce a higher rate of particle collisions than the current Large Hadron Collider (LHC), requiring important upgrades to the Compact Muon Solenoid (CMS) to handle an increased amount of data. An important upgrade is the Phase-2 Outer Tracker Upgrade, which consists of 13,000 silicon sensor modules made of two parallel silicon sensors and readout electronics. These modules undergo careful quality control checks both during and after module assembly to ensure precise and reliable detector performance. This project focuses on precision testing for quality control of silicon sensor modules at Fermilab. Hands-on work includes visual inspection, current-voltage testing, module testing, and ultraviolet (UV) light exposure of modules showing abnormal current-voltage behavior. The ultraviolet exposure process improves the abnormal sensor readout data by placing the selected sensor side of the module directly under the UV light inside a controlled box. In addition to laboratory testing and ultraviolet experiments, I developed a Python-based data tool that connects to a module database and allows selected testing conditions and module information to be retrieved and displayed efficiently. These different testing procedures, experimental processes, and computational tools support the broader goal of identifying module issues and improving modules that will be used in the CMS Outer Tracker Phase-2 Upgrade.

Siddiqui, Hooriya [DuPage Coll.] (ORCID:0009000151↗
DOE OSTI2025

Enhancing ICARUS and REDTOP Software and Hardware: Event Generator Interface Development and Calorimeter Tile Prototype

ICARUS (Imaging Cosmic And Rare Underground Signals) is a liquid argon time projection chamber (LArTPC) detector that pursues the sterile neutrino, which relies on accurate simulations of neutrino-argon interactions. REDTOP (Rare Eta Decays To Observe new Physics) is a proposed low-energy, high-intensity meson factory designed to explore rare $\eta$/$\eta'$ meson decays and probe physics beyond the Standard Model. As a next-generation experiment, this requires both accurate simulations and innovative detector technologies. This project contributes to both ICARUS, from a simulation perspective, and REDTOP, from both a simulation and detection perspective, through the event generation of lepton-nucleon interactions and the physical enhancement of the calorimeter technology within the REDTOP detector. We developed an interface between ACHILLES (A CHIcago Land Lepton Event Simulator), a theory-driven lepton-level event generator, and GENIE, a robust event generator framework used for neutrino physics. By incorporating the precise theoretical cross-section calculations of ACHILLES into the experimental realism of GENIE, the interface allows for improved accuracy of neutrino-nucleon simulations, which can be adapted for the proton beam specifications of the REDTOP meson factory as well as for the ICARUS experiment. In parallel, we developed an improved prototype for the ADRIANO2 (A Dual Readout Integrally Active Non-segmented Option) dual-readout calorimeter tiles for the REDTOP detector. To improve the efficiency of the lead-glass tiles trapping Cherenkov light for energy reconstruction and particle identification, we optimized the application of a highly reflective coating. Through viscosity and thickness control, masking, and a custom spray technique, we refined the coating process to reduce surface defects and improve light yield. Together, these efforts strengthen the ICARUS neutrino program and REDTOP's capability of detecting rare decay events.

Visser, Erin [Michigan State U.] (ORCID:0009000184↗
DOE OSTI2026

Quality Control of Silicon Sensor Modules for Particle Detectors

The High-Luminosity Large Hadron Collider (HL-LHC) will produce a higher rate of particle collisions than the current Large Hadron Collider (LHC), requiring significant upgrades to the Compact Muon Solenoid (CMS) to handle the increased amount of data. An important upgrade is the Phase-2 Outer Tracker Upgrade, which consists of 13,000 silicon sensor modules made of two parallel silicon sensors and readout electronics. These modules undergo careful quality control checks both during and after module assembly to ensure precise and reliable detector performance. This project focuses on precision testing for quality control of silicon sensor modules at Fermilab. Hands-on work includes visual inspection, current-voltage testing, module testing, and ultraviolet (UV) light exposure of modules showing abnormal current-voltage behavior. The ultraviolet exposure process improves the abnormal sensor readout data by placing the selected sensor side of the module directly under the UV light inside a controlled box. In addition to laboratory testing and ultraviolet experiments, I developed a Python-based data tool that connects to a module database and allows selected testing conditions and module information to be retrieved and displayed efficiently. These different testing procedures, experimental processes, and computational tools support the broader goal of identifying module issues and improving modules that will be used in the CMS Outer Tracker Phase-2 Upgrade.

Siddiqui, Hooriya [DuPage Coll.; Fermilab] (ORCID:↗
DOE OSTI2023

Large-Signal Stability Improvement of Parallel Grid-Forming Inverter-Driven Black Start: Preprint

With the rapid increase of inverter-based resources in modern grids, advanced grid-forming (GFM) inverter capa- bilities, such as system restoration and operation under faults, are urgently needed to realize power electronics-dominant grids at scale. One such capability is inverter-driven black start using GFM inverters. This paper analyzes the ability of two recently proposed advanced GFM controls to help GFM inverters sustain system-wide, off-nominal conditions and remain synchronized until they can overcome the momentary overloading as more GFMs join the process without generator sequence coordination or communications and finally stabilize the grid. Through an extensive set of 1,200 full-order electromagnetic transient simu- lations, we evaluate the black-start process while employing the various GFM inverter controls. The results show that the GFM current limiter and primary control have a significant impact on the stability of the system during dynamic operating conditions and thus impact the success of inverter-driven system restoration.

current limit↗
DOE OSTI2022

Editorial for the Virtual Special Issue of Journal of Power Sources “Low temperature fuel cells and electrolysers” - Science and engineering: Let's play this game hand in hand!

Electrochemistry is already and will increasingly be amongst the major drivers of our energetic mix decarbonization, so to fight anthropogenic global warming. The hydrogen economy vector (ideally produced in water electrolysers using renewable electricity and transformed back into electricity in fuel cells), in parallel with the conversion of carbon dioxide and nitrogen gas into value-added products are seen as pillars of the so-called green and sustainable Energy Revolution. The hydrogen technologies presently attract increasing interest, their momentum growing irremediably in Europe, Asia and Americas. Here, this Virtual Special Issue (VSI) gathers the latest scientific results of major research laboratories worldwide, all addressing various aspects of electrochemical energy and material conversion processes.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Dissecting protein tyrosine phosphatase signaling by engineered chemogenetic control of its activity

Protein tyrosine phosphatases (PTPases) are critical mediators of dynamic cell signaling. A tool capable of identifying transient signaling events downstream of PTPases is essential to understand phosphatase function on a physiological time scale. We report a broadly applicable protein engineering method for allosteric regulation of PTPases. This method enables dissection of transient events and reconstruction of individual signaling pathways. Implementation of this approach for Shp2 phosphatase revealed parallel MAPK and ROCK II dependent pathways downstream of Shp2, mediating transient cell spreading and migration. Furthermore, we show that the N-SH2 domain of Shp2 regulates MAPK-independent, ROCK II-dependent cell migration. Engineered targeting of Shp2 activity to different protein complexes revealed that Shp2-FAK signaling induces cell spreading whereas Shp2-Gab1 or Shp2-Gab2 mediates cell migration. We identified specific transient morphodynamic processes induced by Shp2 and determined the role of individual signaling pathways downstream of Shp2 in regulating these events. Broad application of this approach is demonstrated by regulating PTP1B and PTP-PEST phosphatases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Spinbox: tools for many-body quantum systems in a Monte Carlo context

Spinbox is a piece of software that facilitates quantum mechanical calculations relevant to Monte Carlo simulation of atomic nuclei. At the front lines of research on the nuclear many-body problem are a large number of supercomputer-scale simulation codes. These codes produce valuable results but can be hard to understand, especially for those without intimate knowledge of the relevant theoretical methods. Thus, tools that fill pedagogical roles are extremely valuable. Spinbox makes it easy for one to replicate and analyze the computational processes relevant to a Quantum Monte Carlo (QMC) simulation that may be difficult to understand/debug/analyze due to the scale of the corresponding simulation software. Spinbox is written in Python using other state-of-the-art Python modules for numerical calculations. While a number of Python libraries exist that are suited to general quantum many-body calculations, the motivation of Spinbox is quite particular. In Diffusion Monte Carlo methods (DMC, GFMC, AFDMC), the central calculation is the imaginary-time propagation of individual samples of the many-body wavefunction. Although quantum wavefunctions generally must be described by a probability distribution over a basis, DMC imbues particles (within one sample) with classical spatial coordinates. This method is unusual, so other Python packages are typically not set up to do this easily. Furthermore, the software has built-in options for nuclear systems assuming isospin symmetry, which can be set up with other libraries but is a nontrivial process to do so. Features: - numerical representation of samples of the many-body wavefunctions, including tensor-product states (used in AFDMC) - numerical representation of many-body operators, including tensor-product operators: general, spin, imaginary-time propagation, etc. - the correct associated arithmetic and algebra, implemented as class methods - classes for representing realistic nuclear two- and three-body Hamiltonians (e.g. Argonne V18, Illinois NNN) - large-scale parallel integration over random variables, crucial for the AFDMC method My goal is to make this package open source so that anyone may use it and contribute to it, particularly other researchers doing AFDMC calculations

Fox, Jordan↗
DOE OSTI2024

Doppler Lidar Motion-Correction Wind Profiles (DLMCPROF-WIND) Value-Added Product Report

Wind speed and direction, together with pressure, temperature, and relative humidity, are the most fundamental atmospheric state parameters. Accurate measurement of these parameters is crucial for numerical weather prediction. Vertically resolved wind measurements in the atmospheric boundary layer are particularly important for modeling pollutant and aerosol transport. The U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility currently operates several scanning coherent Doppler lidar (DL) systems that provide accurate height-resolved measurements of wind speed and direction. These instruments operate in the near infrared (IR;1.5 microns) and provide range-resolved measurements of radial velocity, attenuated aerosol backscatter, and signal-to-noise ratio (SNR). The systems are operated using a fixed scan schedule consisting of plan position indicator (PPI) or Doppler Beam-Swing (DBS) scans. PPI scans are performed by scanning the beam in azimuth while maintaining a fixed elevation angle, and DBS scans are similar but typically also include a vertical beam. Radial velocity data from these scans are processed to yield profiles of wind speed direction. For stationary ground-based deployments, the Doppler Lidar Horizontal Wind Profiles (DLPROF-WIND) Value-Added-Product (VAP) provides height- and time-resolved measurements of the winds (Newsom and Krishnamurthy 2022). For operation on a moving platform, modifications to the existing DLPROF-WIND VAP are required to compensate for the effects of the platform motion. This report describes a parallel VAP, DLMCPROF-WIND, that was developed for computing motion-compensated wind profiles from ARM Doppler lidar data acquired during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) campaign, when the lidar was deployed on the German icebreaker Polarstern.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Summary of Savannah River Site FY23 Salt Waste Qualification Data

The Savannah River National Laboratory (SRNL) analyzed samples from Savannah River Site (SRS) Waste Tanks 41H and 21H to support qualification of Salt Waste Processing Facility (SWPF) Waste Batches 8 and 9 for processing (the FY23 Salt Batch Qualification samples). These Tanks (i.e. 41H and 21H) are blend tanks for feed to SWPF. None of the samples displayed any unusual or unexpected characteristics such as large amounts of solids, floating solids, or unusual color. Characterization of these samples confirmed similar chemical composition and characteristics to previous salt waste batches. The results for Batches 8 and 9 were provided by SRNL to Savannah River Mission Completion (SRMC), the Liquid Waste Operations subcontractor at SRS, as External Sample Results (Laboratory Information Management System (LIMS)) Reports. Additionally, a separate technical memo was issued by SRNL to report re-test data for Batch 8 for Cs-137 for filtered samples only which were run at the request of SWPF. For Batch 9, a set of samples was also analyzed in parallel by the SWPF-Analytical Laboratory (SWPF-AL). The SWPF-AL data is included herein for comparison with the SRNL data where applicable. The analytical results (both rapid, typically 4 weeks, and long term, typically 8 weeks) for Batches 8 and 9 are now summarized and discussed in this technical report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Influence of energy density on the microstructure, growth orientation, and anisotropy of magnetic properties in additively manufactured Fe-3.8wt%Si transformer steels

Fe-3.8wt%Si transformer steels were processed using two different additive manufacturing (AM) techniques, laser powder bed fusion (LPBF) and directed energy deposition (DED). While the LPBF processed samples exhibited a strong <001> orientation of the BCC grains along the build axis, the DED processed samples exhibited a randomized texture along the build axis. DED processed samples showed substantially coarser columnar grains as compared to their LPBF counterparts. Here, the columnar grains exhibited a substantial number of low-angle sub-grain boundaries. All samples exhibited very good soft magnetic properties, with saturation magnetization (M s ) values ranging from 205 - 232 emu/gm, and coercivity (H c ) values ranging from 1.2 – 4.2 Oe. The Coercivity (H c ) values were significantly lower when the magnetic field was applied parallel to the build axis, as compared to being perpendicular, which can be rationalized based on the columnar nature of the grains, resulting in a higher number density of grain boundaries in case of the field applied perpendicular to the build axis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Study of the structure of exploding flat foils at superhigh current density

We have investigated the features produced in flat metal foils as a result of nanosecond explosions and phase transitions at current densities of (0.1–2) • 10 9 A/cm 2 . Thin foils made of aluminum, copper, nickel, and titanium, with thicknesses of 1–15 μm, were tested. Here, the exploded foil structure was studied using X-pinch radiography and laser shadow imaging. Al, Ti, and Cu (5 and 7 μm thick) foils had an initial two-dimensional structure. At the same time, no apparent structure was registered in 1 μm Cu and 5 μm Ni foils. Experiments on generators with different output parameters have shown that the dominant structures that developed in the explosion are either cracks or bubbles and may depend on the orientation of the initial structure in the thin foil with respect to the current direction. In addition, the energy deposited in the foil differs by a factor of about 1.5 for the orthogonal vs parallel current and initial foil structure directions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Pore-Scale Simulation of Spectral Induced Polarization

A new model was developed to simulate spectral induced polarization (SIP) in porous media. The model is based on a pore-network approximation of the soil or rock, in which bonds between pores are represented by equivalent circuits. In the circuit model, electrolytic and surface conduction occur in parallel, with the electrolyte behaving as a single resistor, and the mineral surface behaving as a resistor and capacitor in series. Based on the pore geometry, fluid conductivity, and mineral surface conductivity, an impedance network is generated. SIP spectra are calculated by solving the sinusoidal steady-state electrical problem for a range of frequencies. The new model provides a tool to bridge the gap between (1) our fundamental understanding of polarization processes at the scale of a single pore or grains and (2) existing simulation capabilities for macroscale systems that do not account for pore-scale properties. The model is implemented in Matlab and can solve two-dimensional and three-dimensional (2D and 3D) networks. Examples for regular cubic lattices are presented, but the method and code are flexible with respect to pore geometries and allow for future leveraging of existing open-source capabilities for generation of pore networks and extraction of network geometries from micro imaging of core samples.

58 GEOSCIENCES↗