Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D equilibrium”

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 91 records · Page 5

On the frequency bifurcations of the MHD startup modes in NSTX

The observed bifurcations of the low frequency (<50 kHz) and low toroidal periodicity (n < 5) magnetohydrodynamic (MHD) activity often present in the initial part of the National Spherical Tokamak Experiment (NSTX) discharges can be explained by the evolution of the radial profile of the safety factor (q=rB φ /RB θ ) crossing multiple rational surfaces in the core. Important performance limiting instability mechanisms in the NSTX spherical tokamak are often linked to low frequency and low-n MHD activity. They are quite common in long-pulse NSTX plasmas. They can be present at the beginning of the plasma current flat-top, at the end of the discharge or during the whole duration, and they have been observed to deleteriously impact performance over a wide range of q 95 . An interesting feature observed in some NSTX discharges is the presence of a bifurcation in the frequency of the low n modes, as low as n = 1, that have frequencies comparable to the plasma core rotation divided by n. Equilibrium reconstructions constrained by magnetic diagnostics data and motional stark effect pitch angle radial profiles suggest that the observed bifurcations are linked to a fast evolving minimum value of q. Finally, 3D non-linear resistive MHD simulations show that these modes are ideal and exist as non-resonant before the correspondent rational surface enters the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A versatile approach to stabilize liquid-liquid interfaces using surfactant self-assembly

Stabilizing liquid–liquid interfaces, whether between miscible or immiscible liquids, is crucial for a wide range of applications, including energy storage, microreactors, and biomimetic structures. In this study, a versatile approach for stabilizing the water-oil interface is presented using the morphological transitions that occur during the self-assembly of anionic, cationic, and nonionic surfactants mixed with fatty acid oils. The morphological transitions underlying this approach are characterized and extensively studied through small-angle X-ray scattering (SAXS), rheometry, and microscopy techniques. Dissipative particle dynamics (DPD) as a simulation tool is adopted to investigate these morphological transitions both in the equilibrium ternary system as well as in the dynamic condition of the water-oil interface. Such a versatile strategy holds promise for enhancing applications such as liquid-in-liquid 3D printing. Moreover, it has the potential to revolutionize a wide range of fields where stabilizing liquid–liquid interfaces not only offers unprecedented opportunities for fine-tuning nanostructural morphologies but also imparts interesting practical features to the resulting liquid shapes. These features include perfusion capabilities, self-healing, and porosity, which could have significant implications for various industries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Systematic analysis of magnetic equilibrium reconstruction with eddy currents on LTX- β

Magnetic equilibrium reconstruction in the Lithium Tokamak Experiment-Beta (LTX-β) is complicated by strong eddy currents and toroidal asymmetries arising from its segmented, close-fitting conducting shell. In the earlier experiment LTX, these three-dimensional (3D) conductor effects significantly distorted diagnostic signals and challenged conventional axisymmetric reconstruction methods. In this work, we demonstrate that accurate plasma equilibria can be recovered on LTX-β by incorporating a small number of dominant eddy current modes derived from realistic conductor models. Using the open-source $\tt{TokaMaker}$ Grad–Shafranov solver, we reconstruct equilibria across a systematically selected set of LTX-β discharges and validate them against the legacy $\tt{PSI-Tri}$ hybrid 2D-3D code. Our results show that fully 2D $\tt{TokaMaker}$ reconstructions achieve significantly improved agreement with flux loop and Mirnov probe measurements, reducing total chi-squared fitting errors, especially during startup. Finally, we develop new hybrid 2D-3D $\tt{TokaMaker}$ reconstructions by integrating the $\tt{ThinCurr}$ 3D eddy current model, which we optimize for significant further reduction of chi-squared in most scenarios. These findings underscore the importance of realistic wall-current modeling in short-pulse tokamaks, and establish a physics-based reconstruction framework that is extensible to devices with complex passive structures and 3D wall interactions.

Grad–Shafranov solver↗

Non-symmetric pinning of topological defects in living liquid crystals

Topological defects, such as vortices and disclinations, play a crucial role in spatiotemporal organization of equilibrium and non-equilibrium systems. The defect immobilization or pinning is a formidable challenge in the context of the out-of-equilibrium system, like a living liquid crystal, a suspension of swimming bacteria in lyotropic liquid crystal. Here we control the emerged topological defects in a living liquid crystal by arrays of 3D-printed microscopic obstacles (pillars). Our studies show that while –1/2 defects may be easily immobilized by the pillars, +1/2 defects remain motile. Due to attraction between oppositely charged defects, positive defects remain in the vicinity of pinned negative defects, and the diffusivity of positive defects is significantly reduced. Experimental findings are rationalized by computational modeling of living liquid crystals. Our results provide insight into the engineering of active systems via targeted immobilization of topological defects.

42 ENGINEERING↗

Measuring Multicomponent Adsorption of Tracer Gases on Natural Zeolites

A natural clinoptilolite sample near the Nevada National Security Site was obtained to study adsorption and retardation on gas transport. Of interest is understanding the competition for adsorption sites that may reduce tracer gas adsorption relative to single-component measurements, which may be affected by the multi-scale pore structure of clinoptilolite. Clinoptilolite has three distinct domains of pore size distributions ranging from nanometers to micrometers: micropores with 0.4–0.7 nm diameters, measured on powders by CO2 adsorption at 273 K, representing the zeolite cages; mesopores with 4–200 nm diameters, observed using liquid nitrogen adsorption at 77 K; and macropores with 300–1000 nm diameters, measured by mercury injection on rock chips (~ 100 mesh), likely representing the microfractures. These pore size distributions are consistent with X-ray computed tomography (CT) and focused ion beam scanning electron microscope (FIB-SEM) images, which are used to construct the three-dimensional (3D) pore network to be used in future gas transport modeling. To quantify tracer gas adsorption in this multi-scale pore structure and multicomponent gas species environment, natural zeolite samples initially in equilibrium in air were exposed to a mixture of tracer gases. As the tracer gases diffuse and adsorb in the sample, the remaining tracer gases outside the sample fractionate. Using a quadrupole mass spectrometer to quantify this fractionation, the degree of adsorption of tracer gases in the multicomponent gas environment and multi-scale pore structure is assessed. The major finding is that Kr reaches equilibrium much faster than Xe in the presence of ambient air, which leads to more Kr uptake than Xe over limited exposure periods. When the clinoptilolite chips were exposed to humid air, the adsorption capability decreases significantly for both Xe and Kr with relative humidity (RH) as low as 3%. Both Xe and Kr reaches equilibrium faster at higher RH. The different, unexpected, adsorption behavior for Xe and Kr is due to their kinetic diameters similar to the micropores in clinoptilolite which makes it harder for Xe to access compared to Kr.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Resolving ECRH deposition broadening due to edge turbulence in DIII-D

Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Vidyut3d: A Non-Equilibrium Plasma Modeling Tool [SWR-24-101]

Vidyut3d is a massively-parallel plasma-fluid solver for low-temperature plasmas (LTPs) that supports both local field (LFA) and local mean energy (LMEA) approximations, as well as complex gas and surface-phase chemistry. The solver supports 2D and 3D domains, and uses AMReX's adaptive mesh refinement capabilities to increase the grid resolution around complex structures (e.g. streamer heads and sheaths) while maintaining a tractable problem size. Vidyut specializes in simulating various types of gas-phase discharges, as well as plasma/surface interactions and surface chemistry (e.g. for plasma-mediated catalysis applications). The solver also supports hybrid CPU/GPU parallelization strategies, and has demonstrated excellent scaling on various HPC architectures for problem sizes consisting of O(100 M) control volumes.

Sitaraman, Hariswaran↗

A new moment-based general-relativistic neutrino-radiation transport code: Methods and first applications to neutron star mergers

ABSTRACT We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in υ/c, retaining all non-linear neutrino–matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino–matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim }70\,$ ms. The latter is the longest dynamical space-time, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the post-merger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries, and equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of temperature and density evolution in MHD simulations of HIT-SI

The helicity injected torus-steady inductive (HIT-SI) experiment uses steady inductive helicity injection to form a spheromak equilibrium and sustain the structure against resistive decay. Helicity injection is performed using two half-tori “injectors” connected to the main plasma volume, whose fields are oscillated in an AC manner. The properties of the sustained spheromak equilibrium have been experimentally observed to vary with the frequency of the injector oscillation, producing higher current gains and more-symmetric and outwardly shifted current centroids with higher frequency. A computational scan of injector frequency using the 3D MHD code PSI-Tet, which models the entire HIT-SI plasma volume including the injectors, has been performed, including a comparison of the results using the full Hall MHD model to those obtained using a simplified “zero-beta” (constant temperature and density) model. The results of both PSI-Tet models are also compared with experimental data and with simulations using the NIMROD code, which does not model the injector regions. The results of the PSI-Tet simulations show that the average temperature and current gain increase with injector frequency, in agreement with experimental trends. The simulations also show qualitative changes in the dynamics of several quantities with increasing injector frequency, such as density oscillations and current evolution. However, the outward shift and symmetrizing of the current centroid, observed experimentally, are not observed in these MHD simulations, indicating that unresolved or excluded dynamics may be important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates

We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high- and low-density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Sparse measurement medical CT reconstruction using multi-fused block matching denoising priors

A major challenge for medical X-ray CT imaging is reducing the number of X-ray projections to lower radiation dosage and reduce scan times without compromising image quality. However these under-determined inverse imaging problems rely on the formulation of an expressive prior model to constrain the solution space while remaining computationally tractable. Traditional analytical reconstruction methods like Filtered Back Projection (FBP) often fail with sparse measurements, producing artifacts due to their reliance on the Shannon-Nyquist Sampling Theorem. Consensus Equilibrium, which is a generalization of Plug and Play, is a recent advancement in Model-Based Iterative Reconstruction (MBIR), has facilitated the use of multiple denoisers are prior models in an optimization free framework to capture complex, non-linear prior information. However, 3D prior modelling in a Plug and Play approach for volumetric image reconstruction requires long processing time due to high computing requirement. Instead of directly using a 3D prior, this work proposes a BM3D Multi Slice Fusion (BM3D-MSF) prior that uses multiple 2D image denoisers fused to act as a fully 3D prior model in Plug and Play reconstruction approach. Our approach does not require training and are thus able to circumvent ethical issues related with patient training data and are readily deployable in varying noise and measurement sparsity levels. In addition, reconstruction with the BM3D-MSF prior achieves similar reconstruction image quality as fully 3D image priors, but with significantly reduced computational complexity. We test our method on clinical CT data and demonstrate that our approach improves reconstructed image quality.

Hossain, Maliha [ORNL]↗

Thermochemical Nonequilibrium Modeling in a Continuous-Galerkin, Finite-Element Framework

The presented work discusses the implementation, verification, and validation of Park's two-temperature model in a scalable, computational fluid dynamics (CFD) code developed at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The implementation of Park's two-temperature model was verified through 0D test cases involving an adiabatic reactor and a nitrogen thermal bath. The implementation was then validated through comparisons with other validated CFD codes and experimental data on a hypersonic cylinder and double cones. These are standard benchmark test cases for thermochemical non-equilibrium (TCNE) modeling, and all data are shared publicly. The verification and validation results showed that ORNL's in-house CFD code could model complex, high-speed flow problems with and without TCNE modeling. This work is essential for future research involving 3D shock wave/boundary layer interactions (SBLIs).

Nutter, Nicole↗

The 3D Kinematics of the Orion Nebula Cluster: NIRSPEC-AO Radial Velocities of the Core Population

The kinematics and dynamics of stellar and substellar populations within young, still-forming clusters provide valuable information for constraining theories of formation mechanisms. Using Keck II NIRSPEC+AO data, we have measured radial velocities for 56 low-mass sources within 4' of the core of the Orion Nebula Cluster (ONC). We also remeasure radial velocities for 172 sources observed with SDSS/APOGEE. These data are combined with proper motions measured using HST ACS/WFPC2/WFC3IR and Keck II NIRC2, creating a sample of 135 sources with all three velocity components. The velocities measured are consistent with a normal distribution in all three components. We measure intrinsic velocity dispersions of (${\sigma }_{{v}_{\alpha }}$, ${\sigma }_{{v}_{\delta }}$, ${\sigma }_{{v}_{r}}$) = (1.64 ± 0.12, 2.03 ± 0.13, ${2.56}_{-0.17}^{+0.16}$) km s –1 . Our computed intrinsic velocity dispersion profiles are consistent with the dynamical equilibrium models from Da Rio et al. (2014) in the tangential direction but not in the line-of-sight direction, possibly indicating that the core of the ONC is not yet virialized, and may require a nonspherical potential to explain the observed velocity dispersion profiles. We also observe a slight elongation along the north–south direction following the filament, which has been well studied in previous literature, and an elongation in the line-of-sight to tangential velocity direction. These 3D kinematics will help in the development of realistic models of the formation and early evolution of massive clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Localizing resonant magnetic perturbations for edge localized mode control in KSTAR

An external 3D magnetic perturbation typically drives a resonant response at the rational surfaces from the core to the edge of tokamak plasmas, due to strong mode coupling and amplification. In this paper, we present a method to isolate the edge from core resonant fields using the ideal perturbed equilibrium code and to design an edge-localized resonant magnetic perturbation (RMP) for effective edge localized mode (ELM) control. A robust feature of the edge-localized RMP is the curtailed response to the field at the low-field-side (LFS) midplane, as opposed to typical RMPs which strongly resonate with the LFS fields. This emphasizes the importance of off-midplane coils to improve ELM control without provoking a large core response that could lead to devastating instabilities. The conceptual design of new ELM control coils based on the edge-localized RMP in KSTAR shows how this new insight can be utilized to enhance the efficiency of our ELM suppression capabilities. Simple window-pane coils matching the edge-localized resonant mode structure substantially expand in the ELM suppression window beyond the existing coil. Further optimization using the flexible optimized coils using space-curves code leads to additional enhancement in the edge-localized control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Drift orbit islands of energetic particles due to 3D fields in ITER

In this study, we report a comparative study of the orbit islands formed by the guiding center drift motion of test energetic particles (EPs), including both high-energy deuterium ions and fusion-born α-particles, with respect to the magnetic islands formed by tracing the field lines, for an ITER plasma representing the 15 MA baseline scenario. The particle drift orbit is modified by the presence of static resonant magnetic perturbations (RMPs) of different toroidal mode numbers n (n = 1 and 3 in this study), forming orbit islands in the Poincaré plane by passing EPs. The key findings are as follows. (i) With an n = 1 RMP field, the size of the orbit islands combined with the total RMP field, including the plasma response, is about three times smaller than that obtained by assuming the corresponding vacuum field. (ii) The orbit island size is not sensitive to the EP energy, and is comparable to those of the magnetic islands. (iii) Passing EPs with outward radial orbital drift form orbit islands that shift inward with respect to the corresponding magnetic islands, and vice versa. This radial shift, ΔΨ, measured in the normalized equilibrium poloidal flux, is quantified as ΔΨ ≃ 0.03X for the ITER baseline plasma and assuming n = 1 RMPs, with $X\equiv {(M/{M}_{\mathrm{p}})}^{1/2}{(E[\mathrm{M}\mathrm{e}\mathrm{V}])}^{1/2}{Z}^{-1}$, where M is the particle mass, Mp is the proton mass, E is the particle energy, and Z is the particle charge number. (iv) Trapped EPs do not form drift orbit islands, even in the presence of the 3D fields produced by the 90 kAt ELM control coil current in ITER, and thus possess good confinement properties. Nevertheless, the deformation of the trapped particle drift orbits in the Poincaré plane shows that the canonical toroidal angular momentum of EPs is no longer conserved in the presence of the RMP fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SCALE Modeling of the Sodium Cooled Fast-Spectrum Advanced Burner Test Reactor

This report documents the modeling and simulation of a sodium-cooled fast reactor (SFR) as part of a U.S. Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with the Oak Ridge National Laboratory code SCALE and the Sandia National Laboratories (SNL) code MELCOR. Based on publicly available benchmark specifications, a fully heterogeneous 3D SCALE model of the 250 MWth Advanced Burner Test Reactor (ABTR) was developed to demonstrate SCALE’s capabilities for full-core reactivity analysis, fuel inventory prediction, and decay heat analysis of an SFR. The benchmark specifications contain modeling details for the ABTR core at the beginning of equilibrium cycle (BOEC) at operating conditions; they were derived from a 2006 preconceptual design report produced by Argonne National Laboratory. The ABTR was designed to demonstrate reactor-based transmutation of transuranics, that is, to “burn” transuranics recovered from light-water reactor (LWR) spent fuel. The ABTR’s fuel is designed to operate in 4 month cycles using uranium/transuranic (U/TRU) metallic fuel, with a TRU content of approximately 20%, at a conversion ratio of approximately 0.6. Various reactivity calculations were performed with SCALE for the ABTR and, where possible, compared with results available in the open literature. Additionally, SCALE was used to perform a full-core depletion calculation over the 4 month cycle to obtain the nuclide inventory at the end of equilibrium cycle (EOEC). These nuclide inventories, decay heat, power profiles, and reactivity feedback coefficients at EOEC represent the initial conditions for analyzing severe accident scenarios with MELCOR.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Tungsten–carbon surface evolution and erosion modeling for a small angle slot divertor in DIII-D

Abstract We modeled tungsten–carbon mixed surface evolution, sputtering erosion, and transport for the tungsten coated region of a small angle slot (SAS) divertor design for the DIII-D tokamak. This divertor concept aims to achieve a closed slot dissipative plasma to minimize heat load and surface erosion, and to study high- Z material performance. Our advanced simulations use coupled ITMC-DYN material mixing/response and 3D full kinetic REDEP/WBC erosion/redeposition code packages, with divertor plasma solution from the SOLPS-ITER package with 4 MW power input. The SAS design geometry and resulting in-slot plasma parameters cause significant differences in predicted sputter/transport from a conventional divertor. For 2% C/D incident plasma ratio, an equilibrium mixed C/W surface is attained at ∼30 s of discharge, from wall sputtered carbon transported to the 10 cm long tungsten divertor region. Tungsten remains exposed to the plasma, but the evolved surface composition varies with different C/D ratios. Tungsten is primarily sputtered from the mixed surface by impinging carbon ions in the +1 to +4 charge states, with some self-sputtering. Redeposition of sputtered tungsten to the divertor is significant, ∼80% along the higher plasma temperature attached plasma SAS entrance region, but this is less than the typically near-unity values for a conventional divertor. Plasma-incident carbon is highly backscattered (∼50%) from the mixed surface, with little redeposition (<10%); this helps maintain tungsten in the surface sputter zone. Carbon is mainly sputtered from the mixed surface by D + ions, also with low redeposition (∼10%–30%). Finally, the modeling shows non-zero but low sputtered tungsten current from the divertor to the core plasma direction. These results appear favorable for effective testing of a tungsten-containing SAS divertor in DIII-D, and extrapolation of mixed-material evolution/response findings to the analogous low- Z /high- Z , Be/W, ITER plasma facing system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MEUMAPPS (C++ Version)

Many materials, metal alloys in particular, have features on the on micrometer or nanometer scale that have a large impact on the properties of the material. These features are known as the microstructure of the material. Understanding why and how the microstructure forms in a material is of fundamental scientific interest as well as of significant technological interest. The capability to predict microstructure evolution in a material allows the intentional design of microstructures and hence the intentional design of material properties. The phase-field method is one of the leading methods for predicting microstructure evolution. One of the most significant problems for phase-field models is their computational expense. Even limited phase-field simulations can easily require thousands of CPU core-hours to complete, which significantly limits their use. This code provides both a general framework for creating scalable, GPU-accelerated phase-field model applications as well as several applications themselves. The code is capable of using hundreds of GPUs efficiently, which greatly reduces the time required to perform simulations. The code is written with an emphasis on performance portability, that is the ability for the code to run efficiently on a number of different computing architectures without modification of the source code. The performance portability of this code is primarily enabled through the use of two libraries, Kokkos (performance portable data structures and execution patterns) and heFFTe (performance portable distributed 3D fast Fourier transforms). The code consists of a core library, applications, and tests. The core library includes shared functionality between applications. This includes interfaces with fast Fourier transform (FFT) libraries such as heFFTe, data structures based on Kokkos, file input and output capabilities, and a solver for infinitesimal strain mechanical equilibrium problems. Five applications are included in the code. The flagship application is the MEUMAPPS-SS application, which implements the Kim-Kim-Suzuki phase-field model for precipitation for an arbitrary number of phases and components in a metal alloy. Five simpler applications are also included that solve the Eshelby inclusion problem, Allen-Cahn equation, the coupled Allen-Cahn and diffusion equations, and the Cahn-Hilliard equation. The code includes two applications to solve the Cahn-Hilliard equation, one with constant-step-size first-order time integration and the second with adaptive high-order time integration.

DeWitt, Stephen [Oak Ridge National Lab. (ORNL), O↗