Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solver”

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

Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

A realistic parametrization of convection and convective boundary mixing in conventional stellar evolution codes is still the subject of ongoing research. Furthermore, to improve the current situation, multidimensional hydrodynamic simulations are used to study convection in stellar interiors. Such simulations are numerically challenging, especially for flows at low Mach numbers which are typical for convection during early evolutionary stages. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. Simulations of convection for a realistic stellar profile are analyzed regarding the properties of convective boundary mixing. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M ⊙ star model. The results obtained with the low-Mach AUSM + -up solver were compared to results when using its non low-Mach variant AUSM B + -up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from 180 × 90 2 to 810 × 540 2 cells and the nuclear energy release was boosted by factors of 3 × 10 3 , 1 × 10 4 , and 3 × 10 4 to study the dependence of the results on these parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
DOE OSTI2021

Modeling of a chain of three plasma accelerator stages with the WarpX electromagnetic PIC code on GPUs

The fully electromagnetic particle-in-cell code WarpX is being developed by a team of the U.S. DOE Exascale Computing Project (with additional non-U.S. collaborators on part of the code) to enable the modeling of chains of tens to hundreds of plasma accelerator stages on exascale supercomputers, for future collider designs. The code is combining the latest algorithmic advances (e.g., Lorentz boosted frame and pseudo-spectral Maxwell solvers) with mesh refinement and runs on the latest computer processing unit and graphical processing unit (GPU) architectures. In this paper, we summarize the strategy that was adopted to port WarpX to GPUs, report on the weak parallel scaling of the pseudo-spectral electromagnetic solver, and then present solutions for decreasing the time spent in data exchanges from guard regions between subdomains. In Sec. IV, we demonstrate the simulations of a chain of three consecutive multi-GeV laser-driven plasma accelerator stages.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Observation of synergy between lower hybrid waves at two frequencies in EAST

Synergistic effects between two frequencies of lower hybrid (LH) waves—operating at 2.45 and 4.6 GHz—were observed in experiment on EAST for the first time. At low density (ne,lin≈2.0×1019 m−3), simultaneous injection of a 65/35 mix of 2.45/4.6 GHz power achieved an lower hybrid current drive efficiency that was 25% higher than what should be expected from the linear combination of the two sources. The experiment was interpreted with time-dependent simulations, using the equilibrium and transport solver TRANSP, coupled with the ray-tracing code GENRAY and the Fokker–Planck solver CQL3D. For each discharge, profiles of current and hard x-rays from simulation and measurement agree within uncertainties. An examination of the electron distribution function indicates that the LH synergy is supported by the increased width of the LH resonance plateau in the simultaneous injection case compared to independent injection.

Choi, W. (ORCID:0000000268912300)↗
DOE OSTI2021

Spatial profiles of photon chemical potential in near-field thermophotovoltaic cells

Emitted photons stemming from the radiative recombination of electron–hole pairs carry chemical potential in radiative energy converters. This luminescent effect can substantially alter the local net photogeneration in near-field thermophotovoltaic cells. Several assumptions involving the luminescent effect are commonly made in modeling photovoltaic devices; in particular, the photon chemical potential is assumed to be zero or a constant prescribed by the bias voltage. The significance of photon chemical potential depends upon the emitter temperature, the semiconductor properties, and the injection level. Hence, these assumptions are questionable in near-field thermophotovoltaic devices with nanoscale separation distances between the emitter and the cell. In the present work, an iterative solver that combines fluctuational electrodynamics with the drift-diffusion model is developed to tackle the coupled photon and charge transport problem, enabling the determination of the spatial profile of photon chemical potential beyond the detailed balance approach. Here, the difference between the results obtained by allowing the photon chemical potential to vary spatially and by assuming a constant value demonstrates the limitations of the conventional approaches. This study is critically important for performance evaluation of near-field thermophotovoltaic systems with nanoscale vacuum gaps.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Numerical study of $\delta$-function current sheets arising from resonant magnetic perturbations

We report general three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially regions of stochastic field lines when islands overlap. Numerically resolving singular current sheets in the ideal magnetohydrodynamics (MHD) limit has been a significant challenge. This work presents numerical solutions of the Hahm–Kulsrud–Taylor (HKT) problem, which is a prototype for resonant singular current sheet formation. The HKT problem is solved by two codes: a Grad–Shafranov (GS) solver and the Stepped Pressure Equilibrium Code (SPEC) code. The GS solver has built-in nested flux surfaces with prescribed magnetic fluxes. The SPEC code implements multi-region relaxed magnetohydrodynamics (MRxMHD), whereby the solution relaxes to a Taylor state in each region while maintaining force balance across the interfaces between regions. As the number of regions increases, the MRxMHD solution appears to approach the ideal MHD solution assuming a continuum of nested flux surfaces. We demonstrate agreement between the numerical solutions obtained from the two codes through a convergence study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hybrid-VPIC: An open-source kinetic/fluid hybrid particle-in-cell code

Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Finally, special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Accurate simulation of direct laser acceleration in a laser wakefield accelerator

In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The traditional particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Furthermore, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. Furthermore, the azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A gas dynamic perspective on particle lifting in electrostatic discharge-like devices using multiphase particle-in-cell approach

Although particle–laden electrostatic discharges are widely used in laboratories as well as in industrial applications, the mechanism of particle lifting for particles initially at rest in such highly unsteady systems is not well understood. A multiphase gas–particle solver is developed using the multiphase particle-in-cell (MP-PIC) approach to emulate the interaction of a compressible shock-dominated gas phase with the dense particle phase. First, the two-way coupled solver is initially used to study the interaction of a planar traveling shock with a vertical curtain of particulates. The gas and particle phase evolution was found to be in good agreement with a similar experimental study in Ling et al. [Phys. Fluids 24, 113301 (2012)]. Second, the MP-PIC code is used to study the interaction of an expanding blast wave with a thick bed of particles. The simulation considered forces such as quasi-steady drag, pressure-gradient, added-mass, Saffman, and Magnus forces. We observe that the vertical liftoff particles close to the shock impingement point in this configuration are associated with the quasi-steady drag, pressure gradient, and added-mass forces. Also, the Saffman lift and Magnus forces contribute to lifting particles located radially farther away from the shock impingement point. In addition, the study finds a decrease in particle lifting efficiency with decreasing plasma kernel length and shock strength.

42 ENGINEERING↗
DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

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

Sparsified time-dependent Fourier neural operators for fusion simulations

This paper presents a sparsified Fourier neural operator for coupled time-dependent partial differential equations (ST-FNO) as an efficient machine learning surrogate for fluid and particle-based fusion codes such as NIMROD (Non-Ideal Magnetohydrodynamics with Rotation - Open Discussion) and GTC (Gyrokinetic Toroidal Code). ST-FNO leverages the structures in the governing equations and utilizes neural operators to represent Green's function-like numerical operators in the corresponding numerical solvers. Once trained, ST-FNO can rapidly and accurately predict dynamics in fusion devices compared with first-principle numerical algorithms. In general, ST-FNO represents an efficient and accurate machine learning surrogate for numerical simulators for multi-variable nonlinear time-dependent partial differential equations, with the proposed architectures and loss functions. The efficacy of ST-FNO has been demonstrated using quiescent H-mode simulation data from NIMROD and kink-mode simulation data from GTC. The ST-FNO H-mode results show orders of magnitude reduction in memory and central processing unit usage in comparison with the numerical solvers in NIMROD when computing fields over a selected poloidal plane. The ST-FNO kink-mode results achieve a factor of 2 reduction in the number of parameters compared to baseline FNO models without accuracy loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Feasibility Study on Implementing a Staggered-Grid Finite Volume Method for System Analysis Code Development Under the MOOSE Framework

Here, this work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key to the test bed is the implementation of high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. The test bed utilized a more flexible code structure to enable the finite volume method implementation and direct interacting with the solver package, instead of using the natively supported finite element method by the framework. Using a suite of selected test problems with different problem sizes and levels of complexity, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. For a complex reactor model, transient simulation was performed using the newly developed finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development.

MOOSE↗
DOE OSTI2020

Diffusion Synthetic Acceleration for Heterogeneous Domains, Compatible with Voids

A standard approach to solving the S N transport equations is to use source iteration with diffusion synthetic acceleration (DSA). Although this approach is widely used and effective on many problems, there remain some practical issues with DSA preconditioning, particularly on highly heterogeneous domains. For large-scale parallel simulation, it is critical that both (a) preconditioned source iteration converges rapidly and (b) the action of the DSA preconditioner can be applied using fast, scalable solvers, such as algebraic multigrid (AMG). For heterogeneous domains, these two interests can be at odds. In particular, there exist DSA diffusion discretizations that can be solved rapidly using AMG, but they do not always yield robust/fast convergence of the larger source iteration. Conversely, there exist robust DSA discretizations where source iteration converges rapidly on difficult heterogeneous problems, but fast parallel solvers like AMG tend to struggle applying the action of such operators. Moreover, very few current methods for the solution of deterministic transport are compatible with voids. This paper develops a new heterogeneous DSA preconditioner based on only preconditioning the optically thick subdomains. The resulting method proves robust on a variety of heterogeneous transport problems, including a linearized hohlraum mesh related to inertial confinement fusion. Moreover, the action of the preconditioner is easily computed using O(1) AMG iterations, convergence of the transport iteration typically requires 2 to 5× fewer iterations than current state-of-the-art “full” DSA, and the proposed method is trivially compatible with voids. On the hohlraum problem, rapid convergence is obtained by preconditioning less than 3% of the mesh elements with five to ten AMG iterations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part I: Flux and Eigenvalue Solutions

Computing in the energy dimension is one of the greatest challenges confronting present-day deterministic neutron transport solvers. Accurately resolving the neutron flux as neutrons downscatter across resonances in the nuclear cross sections currently requires considerable computing power and suffers from approximation errors. Flux uncertainty resulting from the uncertainty of the resonance structure is the single-largest cause of reactivity uncertainty. Any additional reference solution for the critical neutron downscattering problem with resonance phenomena would be a boon to verification and validation of neutronics codes. This paper establishes a benchmark to verify the accuracy of neutron transport criticality solvers along the energy dimension. For the first time, the analytic solution of the flux amplitude is derived in the particular case of an infinite homogeneous medium with isotropic scattering in the center of mass and an arbitrary number of no-threshold, neutral particle reaction resonances (e.g., radiative capture, fission, and resonance scattering). Furthermore, original analytic expressions are established to quantify the discrepancy between the ψ k (E) and ψ α (E) flux amplitudes, respective solutions of the multiplication factor k, or the exponential time-evolution frequency α eigenproblems. The physical study of these relations led to analysis of their first-order relative difference near the criticality condition α=0. Finally, numerical solutions are provided to a benchmark problem constituted of the first resonance of 239 Pu, the 6.67-eV resonance of 238 U, and a scattering isotope with a flat cross section, allowing for the computational verification of the energy resolution of current neutron transport criticality codes. Through these novel results, this analytic benchmark can serve as a reference to verify the energy resolution and sensitivity analysis of neutron transport criticality calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The Random Ray Method Versus Multigroup Monte Carlo: The Method of Characteristics in OpenMC and SCONE

The Random Ray Method (TRRM) is a recently developed approach to solving neutral particle transport problems based on the Method of Characteristics. While the method previously has been implemented only in closed-source or limited-functionality codes, this work describes its implementation in two open-source Monte Carlo codes: OpenMC and SCONE. The random ray implementations required small modifications to the existing Multigroup Monte Carlo (MGMC) solvers, offering a rare venue for redundant, fine-grained, "apples-to-apples" speed and accuracy comparisons between transport methods. To this end, TRRM and MGMC solvers are evaluated against each other using each code's native capabilities on reactor eigenvalue problems with different degrees of energy discretization. On the C5G7 benchmark (featuring only seven energy groups), TRRM achieves a maximum pin power error comparable to or lower than that of MGMC for a given run time. On a problem with 69 energy groups, MGMC is found to scale more efficiently, obtaining a lower pin power error for a given run time. However, the defining difference between the two transport methods is found to be their vastly different uncertainty distributions. Specifically, TRRM is found to maintain similar levels of accuracy and uncertainty throughout the simulation domain whereas MGMC can exhibit orders-of-magnitude greater errors in areas of the problem that feature low neutron flux. For instance, TRRM provided an up to 373 times speed advantage compared with MGMC for computing the flux in low-flux regions in the moderator surrounding the C5G7 core.

42 ENGINEERING↗
DOE OSTI2025

High-Fidelity CFD Assessments of Flow Resistance in a 61-Pin Wire-Wrapped Assembly with Partially Blocked Channels

The examination of thermal-hydraulic behaviors in wire-wrapped rod bundles continues to be an active area of research. The sodium fast reactor, a prominent candidate in next-generation nuclear designs, utilizes a hexagonal configuration of wire-wrapped fuel pins. Here, the potential for channel blockage within this compact arrangement poses a significant safety challenge, spurring a number of recent experimental and computational investigations to evaluate its impact on flow and heat transfer. The present work aims to benchmark the high-fidelity NekRS computational fluid dynamics (CFD) solver in predicting the pressure drops associated with substantial blockages, using available experimental data as a reference. A 61-pin wire-wrapped fuel assembly with two flow blockage configurations has been simulated and investigated at a range of low to moderate Reynolds numbers (487 ≤ Re ≤ 14 600). The NekRS solver demonstrates an exponential reduction of spatial discretization error with increasing polynomial order. The high level of agreement between the numerical results and measured data confirms the accuracy and consistency of the present numerical approach. This benchmark study establishes the capability of NekRS to perform reliable hydrodynamic simulations for sodium fast reactor applications and supports its use in design, licensing, and safety analyses.

CFD Benchmarking↗
DOE OSTI2023

Overview of Liquid-Metal PFC R&D at the University of Illinois Urbana-Champaign

The design and implementation of future flowing liquid-lithium plasma-facing components (LLPFCs) will be dependent on several factors. Of course, one of the most important is the need to be able to deal with high heat fluxes incident on the surface of the LLPFCs, but there are also several other important liquid-metal behaviors that have been identified for their critical impact on the feasibility of a LLPFC. One of these is the ability to constantly wet 100% of the plasma-facing component area and the best way to achieve that. Another key point is knowing and understanding the erosion and corrosion of the surfaces subject to a flowing liquid-lithium system and the ability for hydrogen and helium uptake by the system. The Center for Plasma Material Interactions (CPMI) has been tasked with looking at these various issues. The Mock-up Entry module for EAST device was used to investigate wetting and erosion effects and to design a suitable distribution and collection system with a liquid-lithium loop. The vapor shielding effects of lithium on the surface were also modeled and studied. A model coupling CRANE, an open-source global reaction network solver, and Zapdos, a plasma transport solver, is being developed to better understand the dynamics of the vapor cloud. Experiments on the Magnum-PSI at the Dutch Institute for Fundamental Energy Research have been carried out to study the vapor shielding effect and obtain experimental benchmarks to verify the model. Also, initial experiments using the Hybrid Illinois Device for Research and Applications have been performed to understand the pumping effects of lithium on helium. Experiments with a drop of liquid lithium (~100 mg) into a helium plasma have shown the ability of lithium to take out the cold recycling helium gas as well as hydrogen and oxygen impurity gases. The improvement in plasma performance was significant, and further understanding of this effect will have impacts on how future LLPFCs will be designed. Further investigation into the exact mechanism for helium pumping by lithium needs to be performed in the future. Finally, this paper presents a summary of the results obtained at the CPMI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature (T e ) and ion temperature (T i ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts T e by 46%, with much larger variance than MMM, and underpredicts T i by 25%. As the ratio of kinetic to magnetic field pressure (β) is increased across NSTX discharges, TGLF predicts lower T e and significant flattening of the T i profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both T e and T i are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Furthermore, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.

Lestz, J. B. [General Atomics, San Diego, CA (Unit↗