Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Distributed State Estimation Over Time-Varying Graphs: Exploiting the Age-of-Information

Here, we study the problem of designing a distributed observer for an LTI system over a time-varying communication graph. The limited existing work on this topic imposes various restrictions either on the observation model or on the sequence of communication graphs. In contrast, we propose a single-time-scale distributed observer that works under mild assumptions. Specifically, our communication model only requires strong-connectivity to be preserved over non-overlapping, contiguous intervals that are even allowed to grow unbounded over time. We show that under suitable conditions that bound the growth of such intervals, joint observability is sufficient to track the state of any discrete-time LTI system exponentially fast, at any desired rate. We also develop a variant of our algorithm that is provably robust to worst-case adversarial attacks, provided the sequence of graphs is sufficiently connected over time. The key to our approach is the notion of a "freshness-index" that keeps track of the age-of-information being diffused across the network. Such indices enable nodes to reject stale estimates of the state, and, in turn, contribute to stability of the error dynamics.

42 ENGINEERING↗
DOE OSTI2014

Development and Benchmarking of a Hybrid PIC Code For Dense Plasmas and Fast Ignition (Final Report)

Radiation processes play an important role in the study of both fast ignition and other inertial confinement schemes, such as plasma jet driven magneto-inertial fusion, both in their effect on energy balance, and in generating diagnostic signals. In the latter case, warm and hot dense matter may be produced by the convergence of a plasma shell formed by the merging of an assembly of high Mach number plasma jets. This innovative approach has the potential advantage of creating matter of high energy densities in voluminous amount compared with high power lasers or particle beams. An important application of this technology is as a plasma liner for the flux compression of magnetized plasma to create ultra-high magnetic fields and burning plasmas. HyperV Technologies Corp. has been developing plasma jet accelerator technology in both coaxial and linear railgun geometries to produce plasma jets of sufficient mass, density, and velocity to create such imploding plasma liners. An enabling tool for the development of this technology is the ability to model the plasma dynamics, not only in the accelerators themselves, but also in the resulting magnetized target plasma and within the merging/interacting plasma jets during transport to the target. Welch pioneered numerical modeling of such plasmas (including for fast ignition) using the LSP simulation code. Lsp is an electromagnetic, parallelized, plasma simulation code under development since 1995. It has a number of innovative features making it uniquely suitable for modeling high energy density plasmas including a hybrid fluid model for electrons that allows electrons in dense plasmas to be modeled with a kinetic or fluid treatment as appropriate. In addition to in-house use at Voss Scientific, several groups carrying out research in Fast Ignition (LLNL, SNL, UCSD, AWE (UK), and Imperial College (UK)) also use LSP. A collaborative team consisting of HyperV Technologies Corp., Voss Scientific LLC, FAR-TECH, Inc., Prism Computational Sciences, Inc. and Advanced Energy Systems Inc. joined efforts to develop new physics and numerical models for LSP in several key areas to enhance the ability of LSP to model high energy density plasmas (HEDP). This final report details those efforts. Areas addressed in this research effort include: adding radiation transport to LSP, first in 2D and then fully 3D, extending the EMHD model to 3D, implementing more advanced radiation and electrode plasma boundary conditions, and installing more efficient implicit numerical algorithms to speed complex 2-D and 3-D computations. The new capabilities allow modeling of the dominant processes in high energy density plasmas, and further assist the development and optimization of plasma jet accelerators, with particular attention to MHD instabilities and plasma/wall interaction (based on physical models for ion drag friction and ablation/erosion of the electrodes). In the first funding cycle we implemented a solver for the radiation diffusion equation. To solve this equation in 2-D, we used finite-differencing and applied the parallelized sparse-matrix solvers in the PETSc library (Argonne National Laboratory) to the resulting system of equations. A database of the necessary coefficients for materials of interest was assembled using the PROPACEOS and ATBASE codes from Prism. The model was benchmarked against Prism's 1-D radiation hydrodynamics code HELIOS, and against experimental data obtained from HyperV's separately funded plasma jet accelerator development program. Work in the second funding cycle focused on extending the radiation diffusion model to full 3-D, continued development of the EMHD model, optimizing the direct-implicit model to speed up calculations, add in multiply ionized atoms, and improved the way boundary conditions are handled in LSP. These new LSP capabilities were then used, along with analytic calculations and Mach2 runs, to investigate plasma jet merging, plasma detachment and transport, restrike and advanced jet accelerator design. In addition, a strong linkage to diagnostic measurements was made by modeling plasma jet experiments on PLX to support benchmarking of the code. A large number of upgrades and improvements advancing hybrid PIC algorithms were implemented in LSP during the second funding cycle. These include development of fully 3D radiation transport algorithms, new boundary conditions for plasma-electrode interactions, and a charge conserving equation of state that permits multiply ionized high-Z ions. The final funding cycle focused on 1) mitigating the effects of a slow-growing grid instability which is most pronounced in plasma jet frame expansion problems using the two-fluid Eulerian remap algorithm, 2) extension of the Eulerian Smoothing Algorithm to allow EOS/Radiation modeling, 3) simulations of collisionless shocks formed by jet merging, 4) simulations of merging jets using high-Z gases, 5) generation of PROPACEOS EOS/Opacity databases, 6) simulations of plasma jet transport experiments, 7) simulations of plasma jet penetration through transverse magnetic fields, and 8) GPU PIC code development. The tools developed during this project are applicable not only to the study of plasma jets, but also to a wide variety of HEDP plasmas of interest to DOE, including plasmas created in short-pulse laser experiments performed to study fast ignition concepts for inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The mystery of missing species in atom probe tomography of composite materials

There is a serious problem in atom probe tomography of composite materials such as oxides that even from stoichiometric samples one observes non-stoichiometric ion yields. We present a quantitative model that explains the non-stoichiometry allowing a fit to experimental data of ion yields as a function of applied field to extract activation barriers and prefactors. The numbers are confirmed by density functional theory. We also show that for oxides the missing oxygen is thermally desorbed as neutral O{sub 2}, either directly or associatively. Finally, we suggest methods to improve the experimental setup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Classical transport equations for burning gas-metal plasmas

Thermonuclear inertial confinement fusion plasmas confined by a heavy metal shell may be subject to the mixing of metal into the gas with a resulting degradation of fusion yield. Classical plasma diffusion driven by a number of gradients can provide a physical mechanism to produce atomic mix, possibly in concert with complex hydrodynamic structures and/or turbulence. This paper gives a derivation of the complete dissipative plasma hydrodynamics equations from kinetic theory, for a binary ionic mixture plasma consisting of electrons, e, a light (hydrogenic gas) ion species, i, and a heavy, high ZI plasma metal species, I. A single mean ionization state for the heavy metal, ZI, is assumed to be provided by some independent thermodynamic model of the heavy metal ZI=ZI(ni,nI,Te). The kinetic equations are solved by a generalized Chapman-Enskog expansion that assumes small Knudsen numbers for all species: NKe≡λe/L≪1,NKi≡λi/L≪1. The small electron to ion mass ratio, me/mi≪1, is utilized to account for electron-ion temperature separation, Te≠Ti, and to decouple the electron and ion transport coefficient calculations. This produces a well ordered perturbation theory for the electrons, resulting in the well known “Spitzer” problem of Spitzer and collaborators and solved independently by Braginskii. The formulation in this paper makes clear the inherent symmetry of the transport and gives an analytic solution for all values of the effective charge Zeff, including Zeff<1. The electron problem also determines the ambipolar electric field and the “thermal forces” on both ion species that are needed for the ion kinetic solution. The ion transport problem makes use of the small mass ratio between ion species, mi/mI≪1, to identify an “ion Spitzer problem” that is mathematically identical to that for the electrons but with different thermodynamic forces. The ionic scattering parameter, ΔI≡nIZI2/ni, replaces the Zeff of the electron problem, but has an extended domain, 0≤ΔI<∞, to cover all mixture fractions from the pure gas to the pure metal plasma. The extension of the Spitzer problem to include this extended domain is given in this work. The resulting transport equations for the binary gas-metal plasma mixture are complete and accurate through second order. All transport coefficients are provided in analytic form.

Molvig, Kim↗
DOE OSTI2020

PPINN: Parareal physics-informed neural network for time-dependent PDEs

Physics-informed neural networks (PINNs) encode physical conservation laws and prior physical knowledge into the neural networks, ensuring the correct physics is represented accurately while alleviating the need for supervised learning to a great degree. While effective for relatively short-term time integration, when long time integration of the time-dependent PDEs is sought, the time–space domain may become arbitrarily large and hence training of the neural network may become prohibitively expensive. To this end, we develop a parareal physics-informed neural network (PPINN), hence decomposing a long-time problem into many independent short-time problems supervised by an inexpensive/fast coarse-grained (CG) solver. In particular, the serial CG solver is designed to provide approximate predictions of the solution at discrete times, while initiate many fine PINNs simultaneously to correct the solution iteratively. There is a two-fold benefit from training PINNs with small-data sets rather than working on a large-data set directly, i.e., training of individual PINNs with small-data is much faster, while training the fine PINNs can be readily parallelized. Consequently, compared to the original PINN approach, the proposed PPINN approach may achieve a significant speed-up for long-time integration of PDEs, assuming that the CG solver is fast and can provide reasonable predictions of the solution, hence aiding the PPINN solution to converge in just a few iterations. To investigate the PPINN performance on solving time-dependent PDEs, we first apply the PPINN to solve the Burgers equation, and subsequently we apply the PPINN to solve a two-dimensional nonlinear diffusion–reaction equation. Furthermore, our results demonstrate that PPINNs converge in a few iterations with significant speed-ups proportional to the number of time-subdomains employed.

42 ENGINEERING↗
DOE OSTI2021

Mitigating the Sign Problem through Basis Rotations

Quantum Monte Carlo simulations of quantum many-body systems are plagued by the Fermion sign problem. The computational complexity of simulating Fermions scales exponentially in the projection time β and system size. The sign problem is basis dependent and an improved basis, for fixed errors, leads to exponentially quicker simulations. Herein, we show how to use sign-free quantum Monte Carlo simulations to optimize over the choice of basis on large two-dimensional systems. We numerically illustrate these techniques decreasing the “badness” of the sign problem by optimizing over single-particle basis rotations on one- and two-dimensional Hubbard systems. We find a generic rotation which improves the average sign of the Hubbard model for a wide range of U and densities for L×4 systems. In one example improvement, the average sign (and hence simulation cost at fixed accuracy) for the 16×4 Hubbard model at U/t=4 and n=0.75 increases by exp[8.64(6)β]. For typical projection times of β ≳ 100, this accelerates such simulation by many orders of magnitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Confined Ni-In intermetallic alloy nanocatalyst with excellent coking resistance for methane dry reforming

Carbon dioxide and methane are two main greenhouse gases which are contributed to serious global warming. Fortunately, dry reforming of methane (DRM), a very important reaction developed decades ago, can convert these two major greenhouse gases into value-added syngas or hydrogen. The main problem retarding its industrialization is the seriously coking formation upon the nickel-based catalysts. Herein, a series of confined indium-nickel (In-Ni) intermetallic alloy nanocatalysts (In x Ni@SiO 2 ) have been prepared and displayed superior coking resistance for DRM reaction. The sample containing 0.5 wt.% of In loading (In 0.5 Ni@SiO 2 ) shows the best balance of carbon deposition resistance and DRM reactivity even after 430 h long term stability test. The boosted carbon resistance can be ascribed to the confinement of core–shell structure and to the transfer of electrons from Indium to Nickel in In-Ni intermetallic alloys due to the smaller electronegativity of In. Additionally, both the silica shell and the increase of electron cloud density on metallic Ni can weaken the ability of Ni to activate C–H bond and decrease the deep cracking process of methane. The reaction over the confined InNi intermetallic alloy nanocatalyst was conformed to the Langmuir-Hinshelwood (L-H) mechanism revealed by in situ diffuse reflectance infrared Fourier transform spectroscopy (in-situ DRIFTS). This work provides a guidance to design high performance coking resistance catalysts for methane dry reforming to efficiently utilize these two main greenhouse gases.

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

Integrable symplectic maps with a polygon tessellation

Identifying integrable dynamics remains a formidable challenge, and despite centuries of research, only a handful of examples are known to date. In this article, we explore a distinct form of area-preserving (symplectic) mappings derived from the stroboscopic Poincaré cross section of a kicked rotator—an oscillator subjected to an external force periodically switched on in short pulses. The significance of this class of problems extends to various applications in physics and mathematics, including particle accelerators, crystallography, and studies of chaos. Notably, Suris's theorem constrains the integrability within this category of mappings, outlining potential scenarios with analytic invariants of motion. In this paper, we challenge the assumption of the analyticity of the invariant by exploring piecewise linear transformations on a torus ( T 2 ) and associated systems on the plane ( R 2 ), incorporating arithmetic quasiperiodicity and discontinuities. Introducing a new automated technique, we discovered previously unknown scenarios featuring polygonal invariants that form perfect tessellations and, moreover, fibrations of the plane or torus. This work reveals a novel category of planar tilings characterized by discrete symmetries that emerge from the invertibility of transformations and are intrinsically linked to the presence of integrability. Our algorithm relies on the analysis of the Poincaré rotation number and its piecewise monotonic nature for integrable cases, contrasting with the noisy behavior in the case of chaos, thereby allowing for clear separation. Some of the newly discovered systems exhibit the peculiar behavior of “integrable diffusion,” characterized by infinite and quasirandom hopping between tiles while being confined to a set of invariant segments. Finally, through the implementation of a smoothening procedure, all mappings can be generalized to quasi-integrable scenarios with suppressed volume occupied by chaotic trajectories, thereby opening doors to potential practical applications. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Experiments and numerical modeling of surrogate droplet combustion for a 93 octane gasoline calibration fuel

Gasoline is the world’s most prolific liquid transportation fuel for ground transportation vehicles. Yet, a fundamental understanding of its burning characteristics is still lacking. Complicating this matter is the highly multicomponent nature of gasoline (comprised of hundreds of chemical species) and the seasonable variability of its composition due to variations in the crude oil feedstocks used in refining. The need to bring order to this complex problem has led to development of calibration fuels (Tier II in presented study) which on average will provide repeatable metrics that replicate performance of refined gasoline for various engine designs. This study presents results of experiments using the configuration of an isolated droplet to study the combustion dynamics of calibration fuel and three different surrogates (S4, S5, S6) for the calibration fuel. Droplet burning is under microgravity conditions, where gas transport is determined entirely by fuel evaporation in order to promote spherical symmetry in the burning process. The initial droplet diameters ranged between 0.55 mm and 0.60 mm. The droplets were spark-ignited and their burning histories were recorded by digital video cameras to obtain quantitative data. For each fuel, the evolutions of droplet, flame, and soot-shell diameters were measured to calculate the burn rates, flame-droplet diameter ratios (FSR) and soot-droplet diameter ratios. Additionally, a detailed numerical modeling is presented. This simulation includes multicomponent phase equilibrium, unsteady gas and liquid transport, radiative heat transfer, internal liquid diffusion and detailed combustion chemistry of different fuels during their burning process. Experiment results show that the calibration fuel and surrogates have a close burn rates and FSR, and S5 presents a relatively similar results with calibration fuel. Predictions from the numerical simulation using the surrogate are compared with experiment measurements, which showed a good agreement with the experiment results, including the droplet burning rate and FSR.

Guo, Songtao↗
DOE OSTI2023

Development of Spectroscopic Tools for Nuclear Forensics Signatures Development

Diversion of nuclear materialsfrom legitimate commercial processes for use by nefarious actors is a frightening and surreptitious problem. Smuggled nuclear material can be used to manufacture a simple radiological dispersion device (RDD) and potentially incorporated into the development of a nuclear weapon. With some frequency, smuggled nuclear material is interdicted in foreign countries as actors cross border checkpoints. Nuclear material interdicted in this way can be analyzed using a variety of tools. Plutonium-239 is commonly used as the fissile material in nuclear weapon production as an alternative to uranium-235. Plutonium in a production, refining, or finishing facility will exist in many forms including oxide precursors (PuF4, PuF3, Pu oxalate, etc.), oxide, and metal. The ability to identify plutonium in each of these forms and determine properties, characteristics, or history of these materials is crucial for determining if the processes occurring in a nuclear facility are legitimate or not. This presentation will focus on the development and application ofspectroscopic tools developed at SRNL using doubled-walled cells (DWC) to characterize Pu-bearing compounds. Advances in the use of vibrational spectroscopy techniques (Raman and IR), UV-VIS-SWIR diffuse reflectance spectroscopy, gamma emission spectroscopy, luminescence spectroscopy, laser-induced breakdown spectroscopy in conjunction with XRD and SEM in the characterization of Pu-bearing compounds will be presented with emphasis in material processing history.

Villa-Aleman, Eliel↗
DOE OSTI2021

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗
DOE OSTI2016

High-fidelity simulation-driven model development for coarse-grained computational fluid dynamics

Nuclear reactor safety analysis requires analysis of a broad range of accident scenarios and determining their consequences. For a nuclear power plant behavior, it is impossible to obtain experimental data at sufficiently large scales to support calibration and validation of various system models. In single-phase flow convective problems, high-resolution methods of Computational Fluid Dynamics (CFD) such as Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) can provide high fidelity results under conditions where physical measurements are unavailable. However, such high-resolution simulations are computationally expensive and not suitable for simulation of long transient scenarios in nuclear reactor accidents. In this work, we investigate the use of a high fidelity simulation-driven approach to model sub-grid scale (SGS) effect in Coarse Grained Computational Fluid Dynamics CG-CFD. Notably, fine-mesh simulations are used to construct physics-informed statistical surrogate of coarse-mesh SGS model. For an initial analysis, we consider a case of turbulent natural convection in a volumetrically heated fluid layer with a thermally insulated lower boundary and isothermal upper boundary. This scenario of unstable stratification is relevant to turbulent natural convection in a severe nuclear reactor accident, as well as in containment mixing and passive cooling. For this type of flow, the SGS effect is modeled by an added turbulent diffusivity in the energy equation. It is shown that a global correction for the energy equation is sufficient to achieve a significant improvement to the CG-CFD prediction of thermal mixing in the fluid layer. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Patch2Self2: Self-supervised Denoising on Coresets via Matrix Sketching

Diffusion MRI (dMRI) non-invasively maps brain white matter yet necessitates denoising due to low signal-to-noise ratios. Patch2Self (P2S) employing self-supervised techniques and regression on a Casorati matrix effectively denoises dMRI images and has become the new de-facto standard in this field. P2S however is resource intensive both in terms of running time and memory usage as it uses all voxels (n) from all-but-one held-in volumes (d-1) to learn a linear mapping Phi : \mathbb R ^ n x(d-1) \mapsto \mathbb R ^ n for denoising the held-out volume. The increasing size and dimensionality of higher resolution dMRI acquisitions can make P2S infeasible for large-scale analyses. This work exploits the redundancy imposed by P2S to alleviate its performance issues and inspect regions that influence the noise disproportionately. Specifically this study makes a three-fold contribution: (1) We present Patch2Self2 (P2S2) a method that uses matrix sketching to perform self-supervised denoising. By solving a sub-problem on a smaller sub-space so called coreset we show how P2S2 can yield a significant speedup in training time while using less memory. (2) We present a theoretical analysis of P2S2 focusing on determining the optimal sketch size through rank estimation a key step in achieving a balance between denoising accuracy and computational efficiency. (3) We show how the so-called statistical leverage scores can be used to interpret the denoising of dMRI data a process that was traditionally treated as a black-box. Experimental results on both simulated and real data affirm that P2S2 maintains denoising quality while significantly enhancing speed and memory efficiency achieved by training on a reduced data subset.

Fadnavis, Shreyas↗
DOE OSTI2015

First and second order operator splitting methods for the phase field crystal equation

In this paper, we present operator splitting methods for solving the phase field crystal equation which is a model for the microstructural evolution of two-phase systems on atomic length and diffusive time scales. A core idea of the methods is to decompose the original equation into linear and nonlinear subequations, in which the linear subequation has a closed-form solution in the Fourier space. We apply a nonlinear Newton-type iterative method to solve the nonlinear subequation at the implicit time level and thus a considerably large time step can be used. By combining these subequations, we achieve the first- and second-order accuracy in time. We present numerical experiments to show the accuracy and efficiency of the proposed methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Experimental study on magnetically insulated transmission line electrode surface evolution process under MA/cm current density

The design of high-current density magnetically insulated transmission line (MITL) is a difficult problem of current large-scale Z-pinch device. In particular, a thorough understanding of the MITL electrode surface evolution process under high current density is lacking. On the “QiangGuang-I” accelerator, the load area possesses a low inductance short-circuit structure with a diameter of 2.85 mm at the cathode, and three reflux columns with a diameter of 3 mm and uniformly distributed circumference at the anode. The length of the high density MITL area is 20 mm. A laser interferometer is used to assess and analyze the state of the MITL cathode and anode gap, and their evolution process under high current density. Experimental results indicate that evident current loss is not observed in the current density area at pulse leading edge, and peak when the surface current density reaches MA/cm. Analysis on electrode surface working conditions indicates that when the current leading edge is at 71.5% of the peak, the total evaporation of MITL cathode structure can be realized by energy deposition caused by ohmic heating. The electrode state changes, and diffusion conditions are reflected in the laser interferometer image. The MITL cathode area mainly exists in metal vapor form. The metal vapor density in the cathode central region is higher than the upper limit of laser penetration density (∼4 × 10{sup 21}/cm{sup 3}), with an expansion velocity of ∼0.96 km/s. The metal vapor density in the electrode outer area may lead to evident distortion of fringes, and its expansion velocity is faster than that in the center area (1.53 km/s).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗
DOE OSTI2019

Algebraic Hybridization and Static Condensation with Application to Scalable $H$(div) Preconditioning

We propose here a unified algebraic approach for the practical application and preconditioning of static condensation and hybridization, two popular techniques in finite element discretizations. We demonstrate the use of this algebraic framework for the construction of scalable solvers for problems involving $H$(div)-spaces discretized by conforming (Raviart--Thomas) elements of arbitrary order. We illustrate through numerical experiments the relative performance of the two (in some sense dual) techniques in comparison with a state-of-the-art parallel solver, ADS [T. V. Kolev and P. S. Vassilevski, SIAM J. Sci. Comput., 34 (2012), pp. A3079--A3098], available at http://www.llnl.gov/casc/hypre and http://mfem.org. Based on these results, we recommend the use of the hybridization technique in practice, due to its clearly demonstrated superior performance with increased benefit for higher-order elements.

97 MATHEMATICS AND COMPUTING↗