Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel systems”

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,369 records · Page 76

Big Data Seismology

The discipline of seismology is based on observations of ground motion that are inherently undersampled in space and time. Our basic understanding of earthquake processes and our ability to resolve 4D Earth structure are fundamentally limited by data volume. Currently, Big Data Seismology is an emergent revolution involving the use of large, data-dense inquiries that is providing new opportunities to make fundamental advances in these areas. This article reviews recent scientific advances enabled by Big Data Seismology through the context of three major drivers: the development of new data-dense sensor systems, improvements in computing, and the development of new types of techniques and algorithms. Each driver is explored in the context of both global and exploration seismology, alongside collaborative opportunities that combine the features of long-duration data collections (common to global seismology) with dense networks of sensors (common to exploration seismology). The review explores some of the unique challenges and opportunities that Big Data Seismology presents, drawing on parallels from other fields facing similar issues. Finally, recent scientific findings enabled by dense seismic data sets are discussed, and we assess the opportunities for significant advances made possible with Big Data Seismology. This review is designed to be a primer for seismologists who are interested in getting up-to-speed with how the Big Data revolution is advancing the field of seismology.

58 GEOSCIENCES↗

Polarization-resolved core exciton dynamics in LiF using attosecond transient absorption spectroscopy

The ability to control absorption by modifying the polarization of light presents an exciting opportunity to experimentally determine the orbital alignment of absorption features. Here, attosecond extreme ultraviolet (XUV) transient absorption spectroscopy is used to investigate the polarization dependence of core exciton dynamics in LiF thin films at the Li + K edge. XUV pulses excite electrons from the Li 1s core level into the conduction band, allowing for the formation of a p-orbital-like core exciton, aligned along the XUV light polarization axis. A sub-5 fs near-infrared (NIR) probe pulse then arrives at variable time delays, perturbing the XUV-excited states and allowing the coherence decay of the core exciton to be mapped. The coherence lifetimes are found to be ≈ 2.4 ± 0.4 fs, which is attributed to a phonon-mediated dephasing mechanism as in previous core exciton studies. The differential absorption features are also shown to be sensitive to the relative polarization of the XUV and NIR fields. The parallel NIR probe induces couplings between the initial XUV-excited p-like bright exciton and s-like dark excitons. When crossed pump and probe polarizations are used, the coupling between the bright and dark states is no longer dipole-allowed, and the transient absorption signal associated with the coupling is suppressed by approximately 90%. This interpretation is supported by simulations of a few-level model system, as well as analysis of the calculated band structure. Here, the results indicate that laser polarization can serve as a powerful experimental tool for exploring the orbital alignment of core excitonic states in solid-state materials.

Attosecond laser spectroscopy↗

Efficient Low-Order Refined Preconditioners for High-Order Matrix-Free Continuous and Discontinuous Galerkin Methods

In this paper, we design preconditioners for the matrix-free solution of high-order continuous and discontinuous Galerkin discretizations of elliptic problems based on finite element method--spectral element method (FEM-SEM) equivalence and additive Schwarz methods. The high-order operators are applied without forming the system matrix, making use of sum factorization for efficient evaluation. The system is preconditioned using a spectrally equivalent low-order (p = 1) finite element operator discretization on a refined mesh. The low-order refined mesh is anisotropic and not shape regular in the polynomial degree of the high-order operator, requiring specialized solvers to treat the anisotropy. We make use of an element-structured, geometric multigrid V-cycle with ordered ILU(0) smoothing. The preconditioner is parallelized through an overlapping additive Schwarz method that is robust in h and p. The method is extended to interior penalty and Bassi and Rebay (BR2) discontinuous Galerkin discretizations, for which it is also robust in the size of the penalty parameter. Finally, numerical results are presented on a variety of examples, verifying the uniformity of the preconditioner.

97 MATHEMATICS AND COMPUTING↗

Fuel Fabrication Capability Assessment in Support of Advanced Reactor Deployments

More than 30 U.S. companies are designing a variety of advanced reactor concepts, and several companies are planning to demonstrate their reactor designs in the mid-2020s to late 2030s time frame. In 2020, the U.S. Department of Energy (DOE) announced a series of awards under the Advanced Reactor Demonstration Program (ARDP) to accelerate the successful deployment of 10 of these reactors under three pathways. TerraPower and X-energy were awarded grants under the Advanced Reactor Demonstration Program to deploy their respective Natrium reactor and Xe-100 reactor designs in the next 7–10 years. These demonstrations are in addition to several parallel programs, including the U.S. Department of Defense’s (DoD’s) interest in the development of microreactors, and interest of the National Aeronautics and Space Administration in space nuclear power and propulsion. The National Reactor Innovation Center’s (NRIC’s) mission is to accelerate the demonstration and deployment of advanced reactors; NRIC is partnering with several reactor developers and harnessing the world-class capabilities of the U.S. National Laboratory system to deliver on its mission. Several of these reactor designs will require advanced fuel forms that are not commercially available today, including metal fuel, molten salt fuel, TRi-structural ISOtropic (TRISO) particle fuel, and uranium nitride fuel. Recognizing that there may be potential gaps in the laboratory-scale process development and pilot-scale first-of-a-kind (FOAK) production of these fuel forms leading to delivery of the FOAK cores, NRIC commissioned this study to look at the challenges that need to be overcome for successful deliveries, including the evaluation of existing facilities and the potential need for a new fuel fabrication facility.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Astronomy Commons Platform: A Deployable Cloud-based Analysis Platform for Astronomy

Abstract We present a scalable, cloud-based science platform solution designed to enable next-to-the-data analyses of terabyte-scale astronomical tabular data sets. The presented platform is built on Amazon Web Services (over Kubernetes and S3 abstraction layers), utilizes Apache Spark and the Astronomy eXtensions for Spark for parallel data analysis and manipulation, and provides the familiar JupyterHub web-accessible front end for user access. We outline the architecture of the analysis platform, provide implementation details and rationale for (and against) technology choices, verify scalability through strong and weak scaling tests, and demonstrate usability through an example science analysis of data from the Zwicky Transient Facility’s 1Bn+ light-curve catalog. Furthermore, we show how this system enables an end user to iteratively build analyses (in Python) that transparently scale processing with no need for end-user interaction. The system is designed to be deployable by astronomers with moderate cloud engineering knowledge, or (ideally) IT groups. Over the past 3 yr, it has been utilized to build science platforms for the DiRAC Institute, the ZTF partnership, the LSST Solar System Science Collaboration, and the LSST Interdisciplinary Network for Collaboration and Computing, as well as for numerous short-term events (with over 100 simultaneous users). In a live demo instance, the deployment scripts, source code, and cost calculators are accessible. 4 4 http://hub.astronomycommons.org/

79 ASTRONOMY AND ASTROPHYSICS↗

AXEAP (ARGONNE X-RAY EMISSION PACKAGE)

Argonne X-ray Emission Package (AXEAP), a singular purpose software package for processing X-ray emission (XES) images collected with a 2-dimensional position sensitive pixel array detector, has been developed. AXEAP can rapidly convert XES image files into a spectral form by applying parallel computation and unsupervised machine learning to compute vast amount of image data. Special focus has been placed on designing user-friendly-interface for processing multiple edges, non-resonant and resonant x-ray emission image analysis, in order to make data processing quick and easy. AXEAP is free software and is written in MATLAB, armed with powerful libraries and toolboxes. The software runs on all common operating systems such as Linux, Window, and Mac.

SUN, CHENGJUN↗

Two-Tower Quantum Matrix Chain Multiplication: Trading Qubits for Depth

Matrix chain multiplication -- computing $\mathcal{W} = M^{(0)}\cdots M^{(K-1)}$ where $M^{(k)} \in \mathbb{R}^{P_k \times P_{k+1}}$-- arises in scientific computing, machine learning, and graph analysis. Despite the importance of this problem, for chains of distinct matrices, the classical number of operations grows linearly with the chain length $K$ and polynomially in the matrix dimensions. We present \emph{Two-Tower Matrix Multiplication}, a quantum subroutine that encodes the product $\mathcal{W}$ of the $K$ matrices into a quantum state in circuit depth $\mathcal{O}(\max_{k} \mathrm{polylog} (P_k P_{k+1}))$, which is independent of~$K$ within the QRAM-based state-preparation model, whereas the qubit count is $\mathcal{O}\bigl(\sum_{k} \log P_k \bigr)$; the total gate count remains linear in $K$, so the gain is in the circuit depth. The construction interleaves state-preparation operators across two layers; within each layer, all operators act on disjoint registers and execute in parallel. This subroutine can be specialized for the chain-vector case, which computes the product of $K-1$ matrices applied to a vector. We prove the correctness of the subroutine for all $K$ and provide two implementations using the Qiskit and QCLAB frameworks. The subroutine is applicable to any downstream quantum algorithm that operates on a matrix encoded in the statevector, including norm estimation, graph-matrix powers, linear system solving, and quantum machine learning kernels.

Antonioli, Giacomo [Pisa U.] (ORCID:00090000668703↗

Accelerating Biomimetic Solar - Energy Harvesting: Mapping the Interaction Landscape of Plasmonic-Excitonic Hybrid Nanosystems (Final Report)

In general, excitonic and plasmonic nanoscale materials in close proximity show high potential for significant breakthroughs in energy related materials research. The interactions between these two kinds of materials result in coupled optical transitions (plexcitons), distinct from those of both the individual exciton and plasmon as well as from those of the sum of their constituents (synergistic effects). By linking together materials-research and physical-research approaches, this project contributes to a concerted approach on nanomaterials energy research. The project’s overall goal is to accelerate the development of well-defined plexcitonic model systems consisting of carefully engineered plasmonic and excitonic nanomaterial— essential for both gaining a fundamental understanding of plexcitonic nanomaterials and the development of novel design principles for biomimetic solar energy harvesting. During the 3-year project period and the terminal renewal with limited support for a 12-month period, we successfully synthesized and characterized (1) a robust excitonic nanomaterial and (2) a library of plasmonic nanoparticles as well as developed (3) a microfluidic platform for homogenous nanosynthesis as summarized below: (1) Robust Excitonic Nanomaterial. Supramolecular assemblies are Nature’s most successful material system for solar energy harvesting. However, photovoltaic devices based on artificial supramolecular assemblies continue to be stymied by disappointing efficiencies and poor stability. The conceptual failure may lie in current solar cell architectures, which rely on solidifying supramolecular assemblies as an ensemble into a solid matrix, neglecting the intrinsic fragility of the assemblies’ internal structure, thus disrupting or even destroying their delicate optoelectronic properties, that is, delicate Frenkel excitonic properties. Supramolecular assemblies may finally serve as usable light harvesting material systems for solar energy conversion technologies, only if they meet the following criteria: (a) Stability, that is, the fragile structure including its delicate Frenkel excitonic character needs to be stable, (b) Robustness, that is, resistant against elevated and fluctuating temperatures, and (c) Viability for device integration, that is, capable of being immobilized onto solid substrates. Here, by developing a nanocomposite via a tunable, cage-like scaffold design, we successfully provided stable supramolecular nanocomposites, that inhabit robust Frenkel excitons despite harming environmental conditions such as extreme heat stress. (2) Library of Plasmonic Nanoparticles. Naturally, current models describing plasmonic hybrid quantum states—plasmonic hybridizations—parallel those developed for molecular orbitals, equating individual plasmonic nanostructures with “atoms” and the plasmonic nanoassemblies with “molecules.” In analogy to organic synthesis, a suitably robust fabrication method would allow for “atom-like” manipulation of “molecule-like” plasmonic nanoassemblies; of high value for next-generation energy nanotechnologies. Despite this frequent comparison, current plasmonic nanoassembly fabrication methods favor top-down templating over wet-chemical synthesis, however, achieving precise control over nanostructure’s geometry and surface characteristics remain an art and a scientific challenge. The conceptual failure may lie in the current wet-chemical synthesis paradigm, as it relies on the accessibility of a multi-dimensional synthesis parameter space through limited, rather one-dimensional synthesis procedures by employing step-by-step approaches. Solution-based nanoarchitectonics for rational design of precisely built plasmonic nanoassemblies via solution-based fabrication may finally be possible only if multi-dimensional syntheses approaches are available that allow for comprehensive control over the plasmonic nanomaterials’ (a) Structural Properties and (b) Surface Properties. Here, by developing an innovative multidimensional 1,3-propanediol based polyol synthesis, we successfully provided control over the plasmonic building-block’s geometry (size and shape) together with its surface characteristics. Our results present a critical step toward the vision of a “periodic table-like” system for plasmonic materials based on straightforward wet-chemical syntheses for energy nanotechnologies. Developing deliberate modifications on this synthesis, we generated a library of plasmonic nanostructures covering the vast parameter space—opening the door for fundamental investigation of plexcitonic model systems. (3) Microfluidic Platform for Homogenous Nanosynthesis. Control over structural properties of plexcitonic nanocomposites remains a challenge due to current limitations in nanosynthesis techniques. Slight variations in nanostructure’s geometry impact their optoelectronic properties, demanding precise synthesis beyond the capabilities of solution-based (batch) synthesis processes. In contrast, the small, confined liquid volumes used in microfluidics—a reaction technique where the manipulation of fluids takes place in channels with dimensions of tens of micrometers—allows for homogenous synthesis conditions, providing excellent control of the reaction and, as a result, of the materials’ geopmetry and composition. However, thus far, the majority of plexcitonic systems has been developed via batch synthesis. Here, by successfully developing a two-channel microreactor, our microfluidic-supported synthesis approach combines the advantages of both microfluidics and batch platforms, allowing for precise spatio-temporal control over all synthesis parameters opening the possibility for homogenous nanosythnesis of well-defined plexcitonic model systems.

14 SOLAR ENERGY↗

New opportunities in produced water management: A market-based approach to produced water trading

Produced water (PW) is a byproduct of oil and gas (O&G) production. Obtained alongside the more valuable energy products, PW is usually characterized by high levels of salinity and often contains many contaminants (chemicals, soluble and insoluble oil, organics, etc.) making it unsuitable for release without substantial treatment. Couple this with the fact that PW is typically obtained at multiple times the rate of oil or gas, and the added transport, treatment, and disposal costs become a serious challenge for operators. These realities have led to ad-hoc practices including cooperation between industry competitors to recycle, share, or otherwise mitigate PW costs. The National Energy Technology Laboratory (NETL) in partnership with the Ground Water Protection Council (GWPC) is pursuing novel technology solutions to address PW issues that complement or improve ad-hoc practices adopted by operators. In this paper, we observe that well-established market management practices used in electrical power generation have natural analogues in the PW supply chain. These parallels open up a new line of research where we view PW management as a market equilibrium problem, and explore solutions that foster active and data-based collaboration among operators through market structures similar to power markets, with the ultimate objective of improving PW management costs and recycling rates. Here, we make a case for our observations, present a PW market clearing optimization model that shows how such a market system could operate in the O&G space, and provide an illustrative case study for demonstration.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

The continuous readout stream of the MicroBooNE liquid argon time projection chamber for detection of supernova burst neutrinos

The MicroBooNE continuous readout stream is a parallel readout of the MicroBooNE liquid argon time projection chamber (LArTPC) which enables detection of non-beam events such as those from a supernova neutrino burst. The low energies of the supernova neutrinos and the intense cosmic-ray background flux due to the near-surface detector location makes triggering on these events very challenging. Instead, MicroBooNE relies on a delayed trigger generated by SNEWS (the Supernova Early Warning System) for detecting supernova neutrinos. The continuous readout of the LArTPC generates large data volumes, and requires the use of real-time compression algorithms (zero suppression and Huffman compression) implemented in an FPGA (field-programmable gate array) in the readout electronics. In this paper we present the results of the optimization of the data reduction algorithms, and their operational performance. To demonstrate the capability of the continuous stream to detect low-energy electrons, a sample of Michel electrons from stopping cosmic-ray muons is reconstructed and compared to a similar sample from the lossless triggered readout stream.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Nanosecond machine learning regression with deep boosted decision trees in FPGA for high energy physics

We present a novel application of the machine learning / artificial intelligence method called boosted decision trees to estimate physical quantities on field programmable gate arrays (FPGA). The software package fwXmachina features a new architecture called parallel decision paths that allows for deep decision trees with arbitrary number of input variables. It also features a new optimization scheme to use different numbers of bits for each input variable, which produces optimal physics results and ultraefficient FPGA resource utilization. Problems in high energy physics of proton collisions at the Large Hadron Collider (LHC) are considered. Estimation of missing transverse momentum (E T miss ) at the first level trigger system at the High Luminosity LHC (HL-LHC) experiments, with a simplified detector modeled by Delphes, is used to benchmark and characterize the firmware performance. The firmware implementation with a maximum depth of up to 10 using eight input variables of 16-bit precision gives a latency value of $\mathcal{O}$(10) ns, independent of the clock speed, and $\mathcal{O}$(0.1)% of the available FPGA resources without using digital signal processors.

Instruments & Instrumentation↗

Evolution of the SLATE linear algebra library

SLATE (Software for Linear Algebra Targeting Exascale) is a distributed, dense linear algebra library targeting both CPU-only and GPU-accelerated systems, developed over the course of the Exascale Computing Project (ECP). While it began with several documents setting out its initial design, significant design changes occurred throughout its development. In some cases, these were anticipated: an early version used a simple consistency flag that was later replaced with a full-featured consistency protocol. In other cases, performance limitations and software and hardware changes prompted a redesign. Sequential communication tasks were parallelized; host-to-host MPI calls were replaced with GPU device-to-device MPI calls; more advanced algorithms such as Communication Avoiding LU and the Random Butterfly Transform (RBT) were introduced. Early choices that turned out to be cumbersome, error prone, or inflexible have been replaced with simpler, more intuitive, or more flexible designs. Applications have been a driving force, prompting a lighter weight queue class, nonuniform tile sizes, and more flexible MPI process grids. Of paramount importance has been building a portable library that works across several different GPU architectures – AMD, Intel, and NVIDIA – while keeping a clean and maintainable codebase. Here we explore the evolving design choices and their effects, both in terms of performance and software sustainability.

Gates, Mark↗

Chromatin Changes in Phytochrome Interacting Factor-Regulated Genes Parallel Their Rapid Transcriptional Response to Light

As sessile organisms, plants must adapt to a changing environment, sensing variations in resource availability and modifying their development in response. Light is one of the most important resources for plants, and its perception by sensory photoreceptors (e.g., phytochromes) and subsequent transduction into long-term transcriptional reprogramming have been well characterized. Chromatin changes have been shown to be involved in photomorphogenesis. However, the initial short-term transcriptional changes produced by light and what factors enable these rapid changes are not well studied. Here, we define rapidly light-responsive, Phytochrome Interacting Factor (PIF) direct-target genes (LRP-DTGs). We found that a majority of these genes also show rapid changes in Histone 3 Lysine-9 acetylation (H3K9ac) in response to the light signal. Detailed time-course analysis of transcript and chromatin changes showed that, for light-repressed genes, H3K9 deacetylation parallels light-triggered transcriptional repression, while for light-induced genes, H3K9 acetylation appeared to somewhat precede light-activated transcript accumulation. However, direct, real-time imaging of transcript elongation in the nucleus revealed that, in fact, transcriptional induction actually parallels H3K9 acetylation. Collectively, the data raise the possibility that light-induced transcriptional and chromatin-remodeling processes are mechanistically intertwined. Histone modifying proteins involved in long term light responses do not seem to have a role in this fast response, indicating that different factors might act at different stages of the light response. This work not only advances our understanding of plant responses to light, but also unveils a system in which rapid chromatin changes in reaction to an external signal can be studied under natural conditions.

59 BASIC BIOLOGICAL SCIENCES↗

Probing Particle Impingement in Boilers Using High-Performance Computing with Parallel CPUs and GPUs

The major goals of the project are to calculate and analyze particle impingement within boilers, quantify effects of particulates in boilers, and predict damage rates of boilers under different cycling modes. Collectively, these initiatives develop insight into existing coal plant challenges using advanced modeling tools, particularly those leveraging high-performance computing resources. High-performance CFD computing forms a central theme in this project that will employ a high degree of coordination and communication between these initiatives to realize a final, rigorously sound, and validated computational capability upon completion. These results will create a holistic, comprehensive, systems-level assessment of damage rates under different cycling modes. Together, these objectives will develop critical insight into damage mechanisms in existing coal plant challenges for accurately and efficiently assessing operating performance in fossil energy power plants.

20 FOSSIL-FUELED POWER PLANTS↗

Motion of Molecules in Supramolecular Scaffolds Enhances Bone Regeneration

The regeneration of human tissues is a great scientific challenge and a critical factor to achieve a long healthspan and prevent disabilities due to injury or disease. Materials chemistry can contribute to this goal with the development of bioactive supramolecular systems that can signal cells for regeneration. Recent work in our laboratory using in vivo models of spinal cord injury and cartilage regeneration has demonstrated that the motion of bioactive molecules in supramolecular scaffolds enhances receptor signaling. We report here on a novel molecular strategy to control supramolecular motion in filamentous assemblies using bone regeneration as a functional target. The supramolecular assemblies are composed of monomers that arrange, by design, with either parallel or antiparallel β-sheets, and some of them contain a terminal peptide sequence that binds BMP-2. We found that parallel β-sheet supramolecular assemblies promote greater osteogenic differentiation of progenitor cells in vitro relative to antiparallel assemblies, as well as superior quality of newly regenerated bone in a rat model of spinal fusion. Furthermore, these assemblies drastically reduce the dangerous supraphysiological dose of BMP-2 used clinically for spinal fusion. Here, we attribute the enhanced bioactivity to the weaker nature of hydrogen bonds in parallel relative to antiparallel β-sheet assemblies, which in turn allows greater supramolecular motion and cell signaling of the growth factor-binding molecules.

Anatomy↗

Probing Quantum Telecloning on Superconducting Quantum Processors

Quantum information cannot be perfectly cloned, but approximate copies of quantum information can be generated. Quantum telecloning combines approximate quantum cloning, more typically referred to as quantum cloning, and quantum teleportation. Quantum telecloning allows approximate copies of quantum information to be constructed by separate parties, using the classical results of a Bell measurement made on a prepared quantum telecloning state. Quantum telecloning can be implemented as a circuit on quantum computers using a classical coprocessor to compute classical feedforward instructions using if statements based on the results of a midcircuit Bell measurement in real time. We present universal symmetric optimal 1 → M telecloning circuits and experimentally demonstrate these quantum telecloning circuits for M = 2 up to M = 10, natively executed with real-time classical control systems on IBM Quantum superconducting processors, known as dynamic circuits. We perform the cloning procedure on many different message states across the Bloch sphere, on seven IBM Quantum processors, optionally using the error suppression technique X–X sequence digital dynamical decoupling. Two circuit optimizations are utilized: one that removes ancilla qubits for M = 2,3 , and one that reduces the total number of gates in the circuit but still uses ancilla qubits. Parallel single-qubit tomography with maximum likelihood estimation density matrix reconstruction is used in order to compute the mixed-state density matrices of the clone qubits, and clone quality is measured using quantum fidelity. These results present one of the largest and most comprehensive noisy intermediate-scale quantum computer experimental analyses on (single qubit) quantum telecloning to date. The clone fidelity sharply decreases to 0.5 for M > 5 , but for M = 2 , we are able to achieve a mean clone fidelity of up to 0.79 using dynamical decoupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A scalable framework for efficient coupling of thermal and microstructural simulations in additive manufacturing

Predicting microstructure evolution in metal additive manufacturing (AM) is important for process optimization, but spatiotemporal scale disparities between thermal transport and microstructure evolution create significant challenges for efficient data transfer between simulation codes. To address this, we present Stork, a scalable framework for coupling thermal and microstructural simulations. Stork uses a sparse data representation to identify and store active solidification sub-volumes, enabling highly parallel quad-linear interpolation from coarse thermal grids to fine microstructure grids without large intermediate storage. We demonstrate the framework by coupling the semi-analytic heat transfer code 3DThesis with the time-parallel cellular automata code Toucan. This approach achieves over two orders of magnitude reduction in data generation time and file size compared to prior workflows. Numerical studies show that quad-linear interpolation preserves grain morphology and crystallographic texture in laser powder bed fusion (LPBF) simulations for coarsening ratios up to 16. Overall, Stork provides a scalable pathway for high-throughput, component-scale AM simulations on modern high-performance computing systems.

36 MATERIALS SCIENCE↗

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗