Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel optimization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measurements of Magnetic Field Penetration of Materials for Superconducting Radiofrequency Cavities

Superconducting Radio Frequency (SRF) cavities used in particle accelerators are typically formed from or coated with superconducting materials. Currently high purity niobium is the material of choice for SRF cavities which have been optimized to operate near their theoretical field limits. This brings about the need for significant R&D efforts to develop next generation superconducting materials which could outperform Nb and keep up with the demands of new accelerator facilities. To achieve high quality factors and accelerating gradients, the cavity material should be able to remain in the superconducting Meissner state under high RF magnetic field without penetration of quantized magnetic vortices through the cavity wall. Therefore, the magnetic field at which vortices penetrate in a supercon- ductor is one of the key parameters of merit of SRF cavities. Techniques to measure the onset of magnetic field penetration on thin film samples need to be developed to mitigate the issues with the conventional magnetometry measurements which are strongly influenced by the film orientation and shape and edge effects. The applied magnetic field is also needed to be parallel to the one side of the superconductor to resemble the magnetic field profile at the surface of the SRF cavities operating at fundamental accelerating mode. In this work we report the development of an experimental setup called Magnetic Field Penetration (MFP) magnetometer to measure the field of full flux penetration through bulk, thin films and multilayered superconductors. Our system combines a small superconducting solenoid which can generate the magnetic field up to 500 mT at the sample surface and three Hall probes to detect the full flux penetration through the planner superconductor with 2 inch diameter. This setup was used to study alternative materials which could potentially outperform niobium, as well as SIS multilayer coatings on niobium.

Senevirathne, Iresha↗
DOE OSTI2022

A critical comparison of CH∙∙∙π versus π∙∙∙π interactions in the benzene dimer: obtaining benchmarks at the CCSD(T) level and assessing the accuracy of lower scaling methods

In this study, we have established CCSD(T)/CBS (Complete Basis Set) limits for 3 stationary points on the benzene dimer potential energy surface, corresponding to the π∙∙∙π (parallel displaced or PD(C 2h ), minimum) and CH∙∙∙π (T-shaped or T(C 2v ), transition state) and tilted T-shaped (or TT(Cs), minimum) bonding scenarios considering both the structure and binding energy. The CCSD(T)/CBS binding energies are -2.65 ± 0.02 (PD), -2.74 ± 0.03 (T), and -2.83 ± 0.01 kcal mol -1 (TT). To this end, the CH∙∙∙π is ~0.2 kcal mol -1 stronger than the π∙∙∙π interaction, whereas the tilting of the CH donating benzene molecule with respect to the other benzene is worth 0.1 kcal mol -1 . As previously discussed in the literature, the MP2 level of theory does not provide a close match for either the energy or structure, yet the SCS-MP2 yields structures in excellent agreement with respect to the CCSD(T) result. It is found that the SCS-MI-MP2 also gives optimized structures very close to SCS-MP2 (within ~0.01 Å of the benchmark). Despite the closer match in structure, the spin-biased MP2 methods (SCS-, SCS-MI-, and SOS-MP2) incorrectly predict the relative stabilities of the isomers. That said, none of the spin biased MP2 methods offers a good compromise between energy and structure for the systems examined. Finally, the CCSD(T)/CBS benchmarks were used to assess the performance of 13 DFT functionals selected from different rungs of Jacob's ladder. Several functionals such as TPSS-D3, B3LYP-D3, B97-D, B97-D3, and B2PLYP-D3 provided a good description of the binding energies for both CH∙∙∙π and π∙∙∙π interactions, yielding values within 6% of the CCSD(T)/CBS benchmark values. Unlike the MP2 methods, these functionals correctly predict the relative stability of the PD(C 2h ) and T(C 2v ) dimers. Further, we find that there is no systematic improvement as Jacob's ladder is ascended (increased complexity of functional). The best functionals that result in a good compromise between structure and energy accuracy are B97-D3 and B2PLYP-D3 for both the CH∙∙∙π and π∙∙∙π interaction. Despite the impressive performance of these functionals, a challenge that remains is ensuring the transferability of these density functionals in accurately describing the interaction between dimers of larger aromatic molecules, the latter requiring high-level benchmarks for these systems.

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

Novel Geometric Operations for Linear Programming

This report summarizes the work performed under the project "Linear Programming in Strongly Polynomial Time." Linear programming (LP) is a classic combinatorial optimization problem heavily used directly and as an enabling subroutine in integer programming (IP). Specifically IP is the same as LP except that some solution variables must take integer values (e.g. to represent yes/no decisions). Together LP and IP have many applications in resource allocation including general logistics, and infrastructure design and vulnerability analysis. The project was motivated by the PI's recent success developing methods to efficiently sample Voronoi vertices (essentially finding nearest neighbors in high-dimensional point sets) in arbitrary dimension. His method seems applicable to exploring the high-dimensional convex feasible space of an LP problem. Although the project did not provably find a strongly-polynomial algorithm, it explored multiple algorithm classes. The new medial simplex algorithms may still lead to solvers with improved provable complexity. We describe medial simplex algorithms and some relevant structural/complexity results. We also designed a novel parallel LP algorithm based on our geometric insights and implemented it in the Spoke-LP code. A major part of the computational step is many independent vector dot products. Our parallel algorithm distributes the problem constraints across processors. Current commercial and high-quality free LP solvers require all problem details to fit onto a single processor or multicore. Our new algorithm might enable the solution of problems too large for any current LP solvers. We describe our new algorithm, give preliminary proof-of-concept experiments, and describe a new generator for arbitrarily large LP instances.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

High-rate heavy-ion tracking using MWPCs and PPACs with an Anti-Discharge Unit

This work reports on the development of two robust, heavy-ion beam tracking concepts operating at low pressure (< 15 Torr) for high rate applications (> 200 kHz). The first concept consists of a Multi-Wire Proportional Counter (MWPC) with a central anode consisting of 12 μm Au-plated Tungsten wires spaced 1 mm from each other. The anode grid is sandwiched between two segmented cathodes aligned orthogonally in the two dimensions for (x,y) particle localization. The second detector is a Parallel-Plate Avalanche Counter (PPAC) that uses the same readout geometry as the MWPC, but replaces the wire anode with a 150 nm silver layer deposited on both sides of a thin (< 1 mg/cm 2 ) polypropylene foil. Additionally, the bias circuitry for the PPAC anode central foil is equipped with an Anti-Discharge Unit (ADU) to prevent transitions from proportional operation to streamer formation, thereby avoiding damaging discharges. The localization capability of both detectors was tested with a low-rate alpha-particle source (241-Am). A position resolution of < 1 mm (FWHM) was achieved under stable, high-gas-gain (> 1000) operating conditions. Their performance in terms of detection efficiency as a function of the isotope charge (Z) was determined by irradiating the detectors with a cocktail beam (Z ≤ 15) with energy of ∼ 100 MeV/u. Full detection efficiency is maintained for all available fragments under optimal operational conditions (i.e., voltage bias). Full detection efficiency was achieved at rates above 200 kHz by irradiating the detectors with a 1 cm diameter 238 U beam at an energy of 143 MeV/u.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Flow Loop Study of ECT-Based Volume Fraction Monitoring in Oil–Water Two-Phase Flows

Measurement of phase volume fractions in water-containing multiphase flows is necessary for the optimization of a host of industrial flow processes. Many water-containing multiphase flows can be classified as either water-dispersed or water-continuous mixtures. A recently developed approach based on Hanai’s mixture formula and utilizing electrical capacitance tomography (ECT) sensors have shown good potential for obtaining water volume fraction estimates in two-phase water-containing flows with different water salinity levels. However, the proposed approach was investigated via controlled experiments restricted to static configurations while, in practice, multiphase flows can be dynamic and unpredictable. In this work, we perform a flow loop study of the proposed ECT-based method for volume fraction estimation in oil–water two-phase flows. Here, we evaluate the performance of the proposed method in both water-dispersed and water-continuous flow regimes by employing different types of capacitive sensors in cylindrical arrangements and parallel-plate rectangular arrangements.

02 PETROLEUM↗
DOE OSTI2025

Deep Reinforcement Learning-Based Control of Energy Storage for Interarea Oscillation Damping

With the increasing electricity consumption and lack of transmission investment, today's power systems are operated much closer to their limits, raising concerns of inter-area oscillations that deteriorate the system stability. Here, this article presents a novel energy storage placement and control approach for enhanced damping of interarea oscillations. Combining the residual analysis and dominant mode analysis, we are able to identify the advantageous locations for placing energy storage that achieve improved damping performance. To overcome the challenges, such as fixed control parameters and insufficient damping, we propose to use a deep reinforcement learning-based approach for energy storage control. A state-of-the-art guided surrogate-gradient-based evolutionary strategy is used to train a learning agent in a robust, efficient, and reproducible manner. Parallel computing is also adopted to speed up the training process. The proposed strategy has been tested on both medium and large-scale systems. The proposed methods have demonstrated their effectiveness in mitigating various interarea oscillations within a timeframe of 20 s, thereby averting system collapse and enhancing power grid stability effectively.

25 ENERGY STORAGE↗
DOE OSTI2023

Distributed out-of-memory NMF on CPU/GPU architectures

We propose an efficient distributed out-of-memory implementation of the non-negative matrix factorization (NMF) algorithm for heterogeneous high-performance-computing systems. The proposed implementation is based on prior work on NMFk, which can perform automatic model selection and extract latent variables and patterns from data. In this work, we extend NMFk by adding support for dense and sparse matrix operation on multi-node, multi-GPU systems. The resulting algorithm is optimized for out-of-memory problems where the memory required to factorize a given matrix is greater than the available GPU memory. Memory complexity is reduced by batching/tiling strategies, and sparse and dense matrix operations are significantly accelerated with GPU cores (or tensor cores when available). Input/output latency associated with batch copies between host and device is hidden using CUDA streams to overlap data transfers and compute asynchronously, and latency associated with collective communications (both intra-node and inter-node) is reduced using optimized NVIDIA Collective Communication Library (NCCL) based communicators. Benchmark results show significant improvement, from 32X to 76x speedup, with the new implementation using GPUs over the CPU-based NMFk. Good weak scaling was demonstrated on up to 4096 multi-GPU cluster nodes with approximately 25,000 GPUs when decomposing a dense 340 Terabyte-size matrix and an 11 Exabyte-size sparse matrix of density 10 -6 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Capacity Investment under Bayesian Information Updates at Reporting Periods: Model and Application

We consider capacity addition decisions by a new product manufacturer faced with uncertain technology alternatives. The manufacturing capacity addition and technology development occurs in parallel, with preliminary results from a technology project's success providing valuable information to the manufacturer in adding capacity. We solve a stochastic dynamic program with Bayesian updates to obtain the manufacturer's expected profit‐maximizing capacity investment decision. Our model and applications are motivated by the Critical Materials Institute (CMI) (funded by the Department of Energy), which manages research projects focused on mitigating critical material constraints, vital to renewable energy technologies such as direct‐drive wind turbines, electric vehicles, and energy‐efficient lighting. We capture three unique aspects of the problem: first, the learning from project progress depends on task‐based stochastic outcomes and a project's percent‐done. Second, the underlying technology's profitability is based on a model of competition. Third, we evaluate the impact of progress across a portfolio of projects based on a manufacturer's capacity addition. We develop a heuristic that produces results that are close to optimal and can thus be used for large problem sizes. The managerial insights from an application of our model to CMI projects include: (i) technology projects that report the percent‐done of a project earlier increase expected manufacturer profit; (ii) careful choice of “safe bets,” that is, technologies with low profitability but a high probability of success, can increase expected manufacturer profit; (iii) a portfolio of projects can increase profits significantly over separate project evaluation; and (iv) dynamic management of project resources can increase overall profit.

Vedantam, Aditya↗
DOE OSTI2020

Giant momentum-dependent spin splitting in centrosymmetric low-$\textit{Z}$ antiferromagnets

The energy vs crystal momentum $\textit{E(k)}$ diagram for a solid (band structure) constitutes the road map for navigating its optical, magnetic, and transport properties. By selecting crystals with specific atom types, composition, and symmetries, one could design a target band structure and thus desired properties. A particularly attractive outcome would be to design energy bands that are split into spin components with a momentum-dependent splitting, as envisioned by Pekar and Rashba, enabling spintronic application. The current paper provides “design principles” for wave-vector dependent spin splitting (SS) of energy bands that parallel the traditional Dresselhaus and Rashba spin-orbit coupling (SOC)–induced splitting, but originates from a fundamentally different source—antiferromagnetism. Here, we identify a few generic antiferromagnetic (AFM) prototypes with distinct SS patterns using magnetic symmetry design principles. These tools allow also the identification of specific AFM compounds with SS belonging to different prototypes. A specific compound—centrosymmetric tetragonal MnF 2 —is used via density functional band-structure calculations to quantitatively illustrate one type of AFM SS. Unlike the traditional SOC-induced effects restricted to noncentrosymmetric crystals, we show that antiferromagnetic-induced spin splitting broadens the playing field to include even centrosymmetric compounds, and gives SS comparable in magnitude to the best known (“giant’” SOC effects, even without traditional SOC, and consequently does not rely on the often-unstable high atomic number elements required for high SOC. We envision that use of the current design principles to identify an optimal antiferromagnet with spin-split energy bands would be beneficial for efficient spin-charge conversion and spin-orbit torque applications without the burden of requiring compounds containing heavy elements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Liquid–Liquid Microfluidic Flows for Ultrafast 5-Hydroxymethyl Furfural Extraction

We study the effects of liquid–liquid flow patterns on the extraction of 5-hydroxymethyl furfural from the aqueous phase into an organic phase (methyl isobutyl ketone, ethyl acetate, and 2-pentanol) at ultrashort residence times of 0.2–20 s in a coiled capillary microchannel. Slug, droplet, slug-droplet, dispersed, parallel, annular, and irregular two-phase flow patterns are generated. We evaluate the extraction efficiency and the volumetric mass transfer coefficient via experiments and computational fluid dynamics (CFD). The CFD predictions are in quantitative agreement with the experimental data. The convective and diffusive contributions to the mass transfer of different flow patterns are analyzed. For all solvent pairs, the extraction efficiency varies irregularly with increasing flow rate due to changes in the mass transfer mechanism caused by a transition between flow regimes. Overall, we demonstrate >90% extraction efficiencies with high volumetric mass transfer coefficients ranging from 0.006 to 2.17 s–1. Furthermore, a merit index is defined to assess the trade-off between the fast mass transfer and the pressure drop in biphasic microchannels. In general, the mass transfer rate increases significantly at high flow rates at the expense of a high-pressure drop. Furthermore, the ethyl acetate/water irregular flow pattern is optimal due to a large throughput, extraction efficiency, and mass transfer coefficient, with a minimal pressure drop.

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

Specifications and pre-production of n + -in-p large-format strip sensors fabricated in 6-inch silicon wafers, ATLAS18, for the Inner Tracker of the ATLAS Detector for High-Luminosity Large Hadron Collider

Abstract The ATLAS experiment is constructing new all-silicon inner tracking system for HL-LHC. The strip detectors cover the radial extent of 40 to 100 cm. A new approach is adopted to use p-type silicon material, making the readout in n + -strips, so-called n + -in-p sensors. This allows for enhanced radiation tolerance against an order of magnitude higher particle fluence compared to the LHC. To cope with varying hit rates and occupancies as a function of radial distance, there are two barrel sensor types, the short strips (SS) for the inner 2 and the long strips (LS) for the outer 2 barrel cylinders, respectively. The barrel sensors exhibit a square, 9.8 × 9.8 cm 2 , geometry, the largest possible sensor area from a 6-inch wafer. The strips are laid out in parallel with a strip pitch of 75.5 μm and 4 or 2 rows of strip segments. The strips are AC-coupled and biased via polysilicon resistors. The endcap sensors employ a “stereo-annulus” geometry exhibiting a skewed-trapezoid shapes with circular edges. They are designed in 6 unique shapes, R0 to R5, corresponding to progressively increasing radial extents and which allows them to fit within the petal geometry and the 6-inch wafer maximally. The strips are in fan-out geometry with an in-built rotation angle, with a mean pitch of approximately 75 μm and 4 or 2 rows of strip segments. The eight sensor types are labeled as ATLAS18xx where xx stands for SS, LS, and R0 to R5. According to the mechanical and electrical specifications, CAD files for wafer processing were laid out, following the successful designs of prototype barrel and endcap sensors, together with a number of optimizations. A pre-production was carried out prior to the full production of the wafers. The quality of the sensors is reviewed and judged excellent through the test results carried out by vendor. These sensors are used for establishing acceptance procedures and to evaluate their performance in the ATLAS collaboration, and subsequently for pre-production of strip modules and stave and petal structures.

Instruments & Instrumentation↗
DOE OSTI2024

Early experiences on the OLCF Frontier system with AthenaPK and Parthenon–Hydro

The Oak Ridge Leadership Computing Facility (OLCF) has been preparing the nation's first exascale system, Frontier, for production and end users. Frontier is based on HPE Cray's new EX architecture and Slingshot interconnect and features 74 cabinets of optimized 3rd Gen AMD EPYC CPUs for HPC and AI and AMD Instinct 250X accelerators. As a part of this preparation, “real-world” user codes have been selected to help assess the functionality, performance, and usability of the system. This article describes early experiences using the system in collaboration with the Hamburg Observatory for two selected codes, which have since been adopted in the OLCF test harness. Experiences discussed include efforts to resolve performance variability and per-cycle slowdowns. Results are shown for a performance portable astrophysical magnetohydronamics code, AthenaPK, and a mini-application stressing the core functionality of a performance portable block-structured adaptive mesh refinement framework, Parthenon-Hydro. Here, these results show good scaling characteristics to the full system. At the largest scale, the Parthenon-Hydro miniapp reaches a total of $1.7$ $\times$ $10^{13}$ zone-cycles/s on 9216 nodes (73,728 logical GPUs) at ≈92% weak scaling parallel efficiency (starting from a single node using a second-order, finite-volume method).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Study of the Alfven eigenmodes stability in CFQS plasma using a Landau closure model

The aim of this study is to analyze the stability of the Alfven eigenmodes (AE) in the Chinese First Quasi-axisymmetric Stellarator (CFQS). The AE stability is calculated using the code FAR3d that solves the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moment for the energetic particles (EP) species including the effect of the helical couplings and acoustic modes. The Landau damping and resonant destabilization effects are added in the model by a given closure relation. The simulation results indicate the destabilization of n = 1 to 4 AEs by EP during the slowing down process, particularly n = 1 and n = 2 toroidal AEs (TAE), n = 3 elliptical AE (EAE) and n = 4 non circular AE (NAE). If the resonance is caused by EPs with an energy above 17 keV (weakly thermalized EP), n = 2 EAEs and n = 3 NAEs are unstable. On the other hand, EPs with an energy below 17 keV (late thermalization stage) lead to the destabilization of n = 3 and n = 4 TAEs. The simulations for an off-axis NBI injection indicate the further destabilization of n = 2 to 4 AEs although the growth rate of the n = 1 AEs slightly decreases, so no clear optimization trend with respect to the NBI deposition region is identified. In addition, n = 2, 4 helical AE (HAE) are unstable above an EP β threshold. Also, if the thermal β of the simulation increases (higher thermal plasma density) the AE stability of the plasma improves. Overall, the simulations including the effect of the finite Larmor radius and electron-ion Landau damping show the stabilization of the n = 1 to 4 EAE/NAEs as well as a decrease of the growth rate and frequency of the n = 1 to 4 BAE/TAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Roles of Carburization and Oxidation on Transition-Metal Catalyzed Surface Reactions

Many industrial chemical transformation processes are based on reactions catalyzed by transition metals. In some cases the nano-size particles used as catalysts are synthesized simultaneously or right before the initiation of the main catalytic reaction. Synthesis treatments involve complex processes such as reductions followed by calcinations in oxygenated environments coupled to the presence of gases used for the main reaction. Synthesis details are crucial because they have important consequences on the catalyst shape, size, and chemical composition, and may also affect the structure and composition of the catalyst support. All these effects may alter the main catalyzed reaction in unexpected ways. This work focuses on the synthesis of the overall catalytic system that includes catalyst and support taking place in parallel with the initial stages of the main catalyzed reaction. Two important processes are characterized: carburization and oxidation of the transition metal catalyst during synthesis or initial stages of the main reaction. The objective is to understand how the coking mechanism coupled to surface oxidation and changes in the substrate influence the dry reforming of methane taken as a model reaction. The approach uses first principles computational methods to investigate the catalyst synthesis and characterization and catalytic reaction mechanisms. These studies are complemented with experimental characterization done by the group of Dr. Renu Sharma at the National Institute of Standards and Technology (unfunded collaborator) using high resolution transmission electron microscopy. Specific goals are to identify the role of the synthesis process on the initial nanocatalyst structure and composition, and to determine how that structure and composition may evolve after being exposed to the reactants atmosphere. These studies facilitate a systematic analysis of conditions where the nanoparticle composition can be optimized regarding activity and stability. Moreover, the fundamental study of catalyst/support interactions yields valuable insights for other catalytic processes and is a first step toward controlled nanocatalyst synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Scalable Deep Learning-Based Microarchitecture Simulation on GPUs

Cycle-accurate microarchitecture simulators are essential tools for designers to architect, estimate, optimize, and manufacture new processors that meet specific design expectations. However, conventional simulators based on discrete-event methods often require an exceedingly long time-to-solution for the simulation of applications and architectures at full complexity and scale. Given the excitement around wielding the machine learning (ML) hammer to tackle various architecture problems, there have been attempts to employ ML to perform architecture simulations, such as Ithemal and SimNet. However, the direct application of existing ML approaches to architecture simulation may be even slower due to overwhelming memory traffic and stringent sequential computation logic. This work proposes the first graphics processing unit (GPU)-based microarchitecture simulator that fully unleashes the potential of GPUs to accelerate state-of-the-art ML-based simulators. First, considering the application traces are loaded from central processing unit (CPU) to GPU for simulation, we introduce various designs to reduce the data movement cost between CPUs and GPUs. Second, we propose a parallel simulation paradigm that partitions the application trace into sub-traces to simulate them in parallel with rigorous error analysis and effective error correction mechanisms. Combined, this scalable GPU-based simulator outperforms by orders of magnitude the traditional CPU-based simulators and the state-of-the-art ML-based simulators, i.e., SimNet and Ithemal.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

MPACT Verification With Magnox Reactor Neutronics Progression Problems

MPACT is a state-of-the-art core simulator designed to perform high-fidelity analysis using whole-core, three-dimensional, pin-resolved neutron transport calculations on modern parallel computing hardware. MPACT was originally developed to model light water reactors, and its capabilities are being extended to simulate gas-cooled, graphite-moderated cores such as Magnox reactors. To verify MPACT’s performance in this new application, the code is being formally benchmarked using representative problems. Progression problems are a series of example models that increase in complexity designed to test a code’s performance. The progression problems include both beginning-of-cycle and depletion calculations. Reference solutions for each progression problem have been generated using Serpent 2, a continuous-energy Monte Carlo reactor physics burnup calculation code.Using the neutron multiplication eigenvalue ke as a metric, MPACT’s performance is assessed on each of the progression problems. Initial results showed that MPACT’s multi-group cross section libraries, originally developed for pressurized water reactor problems, were not sufficient to accurately solve Magnox problems. MPACT’s improved performance on the progression problems is demonstrated using this new optimized cross section library.

Luciano, Nicholas↗
DOE OSTI2023

Threat Agnostic Virulence Assessment of Pathogens

Virulence assessment of new, emerging, and engineered pathogens is critical to mounting an appropriate response to a biothreat agent. The capacity of the pathogen to colonize human and harm tissues must be characterized to understand pathogenicity pathways and optimize diagnosis and treatment of resulting disease. Respiratory pathogens are of interest because they can have high transmissibility rates, as observed with the SARS-CoV-2 virus, the causative agent of Covid-19. Current technologies are insufficient to assess threats due to their reliance on systems with only one cell type and on sequencing the pathogen. However, it is known that sequence is not an accurate predictor of function, and sequencing can be unreliable for newly emerged or engineered pathogens. An ideal system would consist of relevant epithelial cell types and an assay sensitive enough to detect changes in host responses that do not rely on DNA sequencing. We chose a system consisting of host lung epithelial cells that can be used to assess the virulence of unknown respiratory pathogens. We interrogated pathogens using this model and assess features of pathogenicity. Our objective is to leverage PNNLs strengths in tissue engineering and proteomics capabilities to build a multiple reaction monitoring (MRM) or parallel reaction monitoring (PRM) liquid chromatography-tandem mass spectrometry assay for human host cell proteins whose abundance is influenced by infection. These responses can were then assessed for relative virulence using pathogen agnostic signatures. When confronted with a pathogen, cells activate dedicated signaling pathways, typically through phosphorylation of regulatory proteins and downstream activation of host cell networks.

59 BASIC BIOLOGICAL SCIENCES↗