Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Partitioning”

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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Tempering of low-temperature bainite

Electron microscopy, X-ray diffraction, and atom probe tomography have been used to identify the changes which occur during the tempering of a carbide-free bainitic steel transformed at 473 K (200 °C). Partitioning of solute between ferrite and thin-films of retained austenite was observed on tempering at 673 K (400 °C) for 30 minutes. After tempering at 673 K (400 °C) and 773 K (500 °C) for 30 minutes, cementite was observed in the form of nanometre scale precipitates. Here, proximity histograms showed that the partitioning of solutes other than silicon from the cementite was slight at 673 K (400 °C) and more obvious at 773 K (500 °C). In both cases, the nanometre scale carbides are greatly depleted in silicon.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Perspectives on Quenching and Tempering 4340 Steel

Steels are ubiquitous due to their affordability and the landscape of useful properties that can be generated for engineering applications. But to further expand the performance envelope, one must be able to understand and control microstructure development by alloying and processing. In this paper we use multiscale, advanced characterization to better understand the structural and chemical evolution of AISI 4340 steel after quenching and tempering (Q&T), including the role of quench rate and short-time, isothermal tempering below 573 K (300 °C), with an emphasis on carbide formation. We compare the microstructure and/or property changes produced by conventional tempering to those produced by higher temperature, short-time “rapid” tempering. We underscore that no single characterization technique can fully capture the subtle microstructure changes like carbon redistribution, transition carbide and/or cementite formation, and retained austenite decomposition that occur during Q&T. Only the use of multiple techniques begins to unravel these complexities. After controlled fast or slow quenching, η transition carbides clearly exist in the microstructure, likely associated with autotempering of this high martensite start temperature (Ms) steel. Isothermal tempering below 598 K (325 °C) results in the relief of carbon supersaturation in the martensite, primarily by the formation of η transition carbides that exhibit a range of carbon levels, seemingly without substitutional element partitioning between the carbide and matrix phases. Hägg transition carbide is present between 300 °C and 325 °C. After conventional tempering at or above 598 K (325 °C) for 2 h, cementite is predominant, but small amounts of cementite are also present in other conditions, even after quenching. Previous work has indicated that silicon (Si) and substitutional elements partition between the cementite, which initially forms under paraequilibrium conditions, and the matrix. Phosphorous (P) may also be preferentially located at cementite/matrix interfaces after high temperature tempering. Slower quench rates result in greater amounts of retained austenite compared to those after fast quenching, which we attribute to increased austenite stability resulting from “autopartitioning”. Rapid, high temperature tempering is also found to diminish tempered martensite embrittlement (TME) believed to be associated with the extent of austenite decomposition, resulting in mechanical properties not attainable by conventional tempering, which may have important implications with respect to industrial heat treatment processes like induction tempering. Controlling the amount and stability of retained austenite is not only relevant to the properties of Q&T steels, but also next-generation advanced high strength steels (AHSS) with austenite/martensite mixtures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Mechanism for Carbon Depletion at Bondline of High-Frequency Electric-Resistance-Welded X70 Pipeline Steel

The bondline of electric-resistance-welded (ERW) linepipe steel, often etched white (i.e., ferrite) in optical microscopy, is generally believed to be carbon depleted. The mechanism for the carbon depletion, however, is not fully understood by researchers. To this end, here we used atom probe tomography (APT) to measure elemental segregation of the as-welded and post-weld heat-treated bondline regions of X70 linepipe welds. The thin vertical features at the bondline in the as-welded condition were identified as carbon-rich martensite-austenite (M-A) constituents, and the majority ferrite phase in the bondline was identified as carbon-depleted ferrite. Following the post-weld normalization, all alloying elements, except Nb and Mo, are homogenized across the bondline and heat-affected zone. The carbon depletion in the ERW bondline was accurately measured. A new mechanism for carbon depletion has been proposed using Scheil calculations of elemental partitioning during weld formation. Segregation of elements in the heat-affected zone was shown to follow the negligible partitioning local equilibrium (NPLE) kinetics for bainite transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Influence of Molybdenum on Rapid Solidification Microstructure and Microsegregation in Primary Ferrite Solidified Stainless Steel

Due to the widespread use and familiarity of austenitic stainless steels in many industries these alloys are often the first materials used when adapting new processing techniques. The rapid cooling and solidification rates that occur during fusion based additive manufacturing processes like powder bed fusion cause a fundamental shift in the solidification behavior of the material. While the effects that Cr and Ni have on the microstructure have been fairly well documented across a range of solidification conditions, the impact of alloying elements like Mo are not as well understood at rapid solidification rates. For this study, four custom stainless steel alloy feedstocks were made with targeted modification of the Cr and Mo concentrations to maintain a constant Cr/Ni eq ratio of ~1.7. Two-piston splat quenching was used to produce rapidly solidified samples and the solidification mode, phase, and cell size were investigated. Solidification rates were estimated to be between 0.4 and 1.5m/s. A comparison of the microsegregation and partitioning behavior of Cr, Ni and Mo in ferrite solidified material was also performed for the different alloy compositions. Variations of the Mo concentration between each feedstock showed no measurable influence on the solidification/cooling rates, solidification morphology, or the partitioning of Mo. Furthermore, higher concentrations of Mo were found to significantly increase the amount of ferrite phase in the microstructure by suppressing the ferrite to austenite massive transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Probing the Dissociation of Protein Complexes by Means of Gas-Phase H/D Exchange Mass Spectrometry

Gas-phase hydrogen/deuterium exchange measured by mass spectrometry (gas-phase HDX-MS) is a fast method to probe the conformation of protein ions. The use of gas-phase HDX-MS to investigate the structure and interactions of protein complexes is however mostly unharnessed. Ionizing proteins under conditions that maximize preservation of their native structure (native MS) enables the study of solution-like conformation for milliseconds after electrospray ionization (ESI), which enables the use of ND{sub 3}-gas inside the mass spectrometer to rapidly deuterate heteroatom-bound non-amide hydrogens. Here, we explored the utility of gas-phase HDX-MS to examine protein-protein complexes and inform on their binding surface and the structural consequences of gas-phase dissociation. Protein complexes ranging from 24 kDa dimers to 395 kDa 24mers were analyzed by gas-phase HDX-MS with subsequent collision-induced dissociation (CID). The number of exchangeable sites involved in complex formation could, therefore, be estimated. For instance, dimers of cytochrome c or α-lactalbumin incorporated less deuterium/subunit than their unbound monomer counterparts, providing a measure of the number of heteroatom-bound side-chain hydrogens involved in complex formation. We furthermore studied if asymmetric charge-partitioning upon dissociation of protein complexes caused intermolecular H/D migration. In larger multimeric protein complexes, the dissociated monomer showed a significant increase in deuterium. This indicates that intermolecular H/D migration occurs as part of the asymmetric partitioning of charge during CID. We discuss several models that may explain this increase deuterium content and find that a model where only deuterium involved in migrating charge can account for most of the deuterium enrichment observed on the ejected monomer. In summary, the deuterium content of the ejected subunit can be used to estimate that of the intact complex with deviations observed for large complexes accounted for by charge migration. .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Generalized trigonometric Fourier-series method with automatic time step control for non-linear point kinetics equations

The trigonometric Fourier-series method (TFS) is generalized to provide approximate solutions for non-linear point kinetics equations with feedback using varying step sizes. This method can provide a very stable solution against the size of the discrete time step allowing much larger step sizes to be used. Systems of the point kinetics equations are solved using Fourier-series expansion over a partition of the total time interval. The approximate solution requires determining the series coefficients over each time step in that partition. These coefficients are determined using the high-order derivatives of the solution vector at the beginning of the time step introducing a system of linear algebraic equations to be solved at each step. This system is similar to the Vandermonde system. Two successive orders of the partial sums are used to estimate the local truncation error. This error and some other constrains are used to produce the largest step size allowable at each step while keeping the error within a specific tolerance. The process of calculating suitable step sizes should be automatic and inexpensive. Convergence and stability of the proposed method are discussed and a new formula is introduced to maintain stability. The proposed method solves the general linear and non-linear kinetics problems. The method has been applied to five different types of reactivities including step/ramp insertions with temperature feedback. The method is seemed to be valid for larger time intervals than those used in the conventional numerical integration, and is thus useful in reducing computing time. Computational results are found to be consistent with the analysis, they demonstrate that the convergence of the iteration scheme can be accelerated and the resulting computing times can be greatly reduced while maintain computational accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Deformation mechanisms in a precipitation-strengthened ferritic super alloy revealed by in situ neutron dffraction studies at elevated temperatures

The ferritic superalloy Fe–10Ni–6.5Al–10Cr–3.4Mo strengthened by ordered (Ni,Fe)AlB2-type precipitates is a candidate material for ultra-supercritical steam turbine applications above 923 K. Despite earlier success in improving its room-temperature ductility, the creep resistance of this material at high temperatures needs to be further improved, which requires a fundamental understanding of the high-temperature deformation mechanisms at the scales of individual phases and grains. In situ neutron diffraction has been utilized to investigate the lattice strain evolution and the microscopic load-sharing mechanisms during tensile deformation of this ferritic superalloy at elevated temperatures. Finite-element simulations based on the crystal plasticity theory are employed and compared with the experimental results, both qualitatively and quantitatively. Based on these interphase and intergranular load-partitioning studies, it is found that the deformation mechanisms change from dislocation slip to those related to dislocation climb, diffusional flow and possibly grain boundary sliding, below and above 873 K, respectively. Insights into microstructural design for enhancing creep resistance are also discussed.

nuclear (including radiation effects), defects, me↗
DOE OSTI2021

Enhanced work hardening from oxygen-stabilized ω precipitates in an aged metastable β Ti-Nb alloy

High levels of oxygen in solid solution in Ti alloys are considered detrimental to mechanical properties because of embrittlement concerns. In metastable β titanium alloys, the formation of isothermal ω precipitates is also known to cause severe embrittlement and ductility reduction. However, oxygen has been shown to partition to the ω phase during ageing, and this partitioning behavior may potentially impact ω’s mechanical contribution. Using micropillar compression, we compared the deformation behavior of Ti-20Nb (at. %) with oxygen-stabilized ω precipitates to the behavior of oxygen-free specimens. The oxygen-stabilized microstructures showed increased compressive yield strength and enhanced work hardening behavior compared to oxygen-free specimens. In the absence of oxygen, the compressed pillars showed slip band formation and catastrophic failure, and transmission electron microscopy imaging revealed that ω precipitates were sheared within the continuous deformation channels resulting in slip localization. In contrast, oxygen-stabilized ω precipitates were harder to shear and the formation of continuous deformation channels was suppressed during compression, leading to improved work hardening behavior up to 15% strain. Furthermore, this counter-intuitive role of oxygen may offer design strategies to address the significant embrittlement and loss of ductility observed for ω-strengthened β Ti alloys without oxygen and avenues to expand the use of β Ti alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Local quanta, unitary inequivalence, and vacuum entanglement

In this work we develop a formalism for describing localised quanta for a real-valued Klein–Gordon field in a one-dimensional box [0,R]. We quantise the field using non-stationary local modes which, at some arbitrarily chosen initial time, are completely localised within the left or the right side of the box. In this concrete set-up we directly face the problems inherent to a notion of local field excitations, usually thought of as elementary particles. Specifically, by computing the Bogoliubov coefficients relating local and standard (global) quantisations, we show that the local quantisation yields a Fock representation of the Canonical Commutation Relations (CCR) which is unitarily inequivalent to the standard one. In spite of this, we find that the local creators and annihilators remain well defined in the global Fock space F{sup G}, and so do the local number operators associated to the left and right partitions of the box. We end up with a useful mathematical toolbox to analyse and characterise local features of quantum states in F{sup G}. Specifically, an analysis of the global vacuum state |0{sub G}〉∈F{sup G} in terms of local number operators shows, as expected, the existence of entanglement between the left and right regions of the box. The local vacuum |0{sub L}〉∈F{sup L}, on the contrary, has a very different character. It is neither cyclic (with respect to any local algebra of operators) nor separating and displays no entanglement between left and right partitions. Further analysis shows that the global vacuum also exhibits a distribution of local excitations reminiscent, in some respects, of a thermal bath. We discuss how the mathematical tools developed herein may open new ways for the analysis of fundamental problems in local quantum field theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Braiding statistics and link invariants of bosonic/fermionic topological quantum matter in 2+1 and 3+1 dimensions

Topological Quantum Field Theories (TQFTs) pertinent to some emergent low energy phenomena of condensed matter lattice models in 2+1 and 3+1 dimensions are explored. Many of our TQFTs are highly-interacting without free quadratic analogs. Some of our bosonic TQFTs can be regarded as the continuum field theory formulation of Dijkgraaf–Witten twisted discrete gauge theories. Other bosonic TQFTs beyond the Dijkgraaf–Witten description and all fermionic TQFTs (namely the spin TQFTs) are either higher-form gauge theories where particles must have strings attached, or fermionic discrete gauge theories obtained by gauging the fermionic Symmetry-Protected Topological states (SPTs). We analytically calculate both the Abelian and non-Abelian braiding statistics data of anyonic particle and string excitations in these theories, where the statistics data can one-to-one characterize the underlying topological orders of TQFTs. Namely, we derive path integral expectation values of links formed by line and surface operators in these TQFTs. The acquired link invariants include not only the familiar Aharonov–Bohm linking number, but also Milnor triple linking number in 3 dimensions, triple and quadruple linking numbers of surfaces, and intersection number of surfaces in 4 dimensions. We also construct new spin TQFTs with the corresponding knot/link invariants of Arf(–Brown–Kervaire), Sato–Levine and others. We propose a new relation between the fermionic SPT partition function and the Rokhlin invariant. As an example, we can use these invariants and other physical observables, including ground state degeneracy, reduced modular $S^{xy}$ and $T^{xy}$ matrices, and the partition function on $RP^3$ manifold, to identify all $\nu \in Z_8$ classes of 2+1 dimensional gauged $Z_2$-Ising-symmetric $Z^f_2$-fermionic Topological Superconductors (realized by stacking v layers of a pair of chiral and anti-chiral -wave superconductors [ p + ip and p - ip ], where boundary supports non-chiral Majorana–Weyl modes) with continuum spin-TQFTs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Superintegrable systems from block separation of variables and unified derivation of their quadratic algebras

Highlights: • A new method for constructing D-dimensional minimally superintegrable systems. • Superintegrable systems from block coordinate separation of variables. • Unified derivation of quadratic symmetry algebras via gauge transformations. - Abstract: We present a new method for constructing D-dimensional minimally superintegrable systems based on block coordinate separation of variables. We give two new families of superintegrable systems with N (N≤D) singular terms of the partitioned coordinates and involving arbitrary functions. These Hamiltonians generalize the singular oscillator and Coulomb systems. We derive their exact energy spectra via separation of variables. We also obtain the quadratic algebras satisfied by the integrals of motion of these models. We show how the quadratic symmetry algebras can be constructed by novel application of the gauge transformations from those of the non-partitioned cases. We demonstrate that these quadratic algebraic structures display a universal nature to the extent that their forms are independent of the functions in the singular potentials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Sectoral and growth rate control on elemental uptake by individual calcite crystals

Heterogeneous distribution of trace elements (impurities) within individual calcite crystals is a phenomenon commonly observed in natural and laboratory systems. Changes in thermodynamic intensive parameters (mostly chemical potential and temperature) cannot always explain the inhomogeneous impurity patterns in calcite crystals and it has been suggested that growth rate and crystallographic orientation may exert strong effects on elemental incorporation into calcite. In addition, there are a number of experimental studies on micro-scale element (E) distribution between non-equivalent pairs of calcite vicinal faces (known as sectoral zoning); however, the variability of partition coefficients (e.g., K E = (E/Ca) calcite /(E/Ca) fluid ) within individual crystals remains undetermined. In this study, we have extended the work on elemental distribution between crystal sectors to evaluation of partition coefficients of trace and minor elements (Li, B, Mg, and Sr) in calcite crystal faces (10–14) and (01-12), whose growth rates were assessed. Growth entrapment model (GEM) and lattice strain theory were applied to explain K E heterogeneity by varying near-surface diffusivity of Mg and Sr and by varying surface enrichment factor for Li, B, Mg, and Sr. Decoupling of sectoral and growth rate effects reveals that sectoral zoning plays a key role in elemental distribution. More specifically, K Li and K B vary by more than one order of magnitude and K Mg varies by a factor of two within individual crystal faces. Strontium sectoral distribution is different from those of Li, B, and Mg and K Sr varies by up to a factor of two. Finally, these behaviors likely reflect different mechanisms of incorporation of Li, B, Mg, and Sr into calcite.

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

General field evaluation in high-order meshes on GPUs

Robust and scalable function evaluation at any arbitrary point in the finite/spectral element mesh is required for querying the partial differential equation solution at points of interest, comparison of solution between different meshes, and Lagrangian particle tracking. This is a challenging problem, particularly for high-order unstructured meshes partitioned in parallel with MPI, as it requires identifying the element that overlaps a given point and computing the corresponding reference space coordinates. Here, we present a robust and efficient technique for general field evaluation in large-scale high-order meshes with quadrilaterals and hexahedra. In the proposed method, a combination of globally partitioned and processor-local maps are used to first determine a list of candidate MPI ranks, and then locally candidate elements that could contain a given point. Next, element-wise bounding boxes further reduce the list of candidate elements. Finally, Newton’s method with trust region is used to determine the overlapping element and corresponding reference space coordinates. Since GPU-based architectures have become popular for accelerating computational analyses using meshes with tensor-product elements, specialized kernels have been developed to utilize the proposed methodology on GPUs. The method is also extended to enable general field evaluation on surface meshes. The paper concludes by demonstrating the use of the proposed method in various applications ranging from mesh-to-mesh transfer during r-adaptivity to Lagrangian particle tracking.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Controlling the numerical Cerenkov instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction

In this study we present a customized finite-difference-time-domain (FDTD) Maxwell solver for the particle-in-cell (PIC) algorithm. The solver is customized to effectively eliminate the numerical Cerenkov instability (NCI) which arises when a plasma (neutral or non-neutral) relativistically drifts on a grid when using the PIC algorithm. We control the EM dispersion curve in the direction of the plasma drift of a FDTD Maxwell solver by using a customized higher order finite difference operator for the spatial derivative along the direction of the drift (1ˆ direction). We show that this eliminates the main NCI modes with moderate |k 1 |, while keeps additional main NCI modes well outside the range of physical interest with higher |k 1 |. These main NCI modes can be easily filtered out along with first spatial aliasing NCI modes which are also at the edge of the fundamental Brillouin zone. The customized solver has the possible advantage of improved parallel scalability because it can be easily partitioned along 1ˆ which typically has many more cells than other directions for the problems of interest. We show that FFTs can be performed locally to current on each partition to filter out the main and first spatial aliasing NCI modes, and to correct the current so that it satisfies the continuity equation for the customized spatial derivative. This ensures that Gauss’ Law is satisfied. Lastly, we present simulation examples of one relativistically drifting plasma, of two colliding relativistically drifting plasmas, and of nonlinear laser wakefield acceleration (LWFA) in a Lorentz boosted frame that show no evidence of the NCI can be observed when using this customized Maxwell solver together with its NCI elimination scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Scalable implementation of polynomial filtering for density functional theory calculation in PARSEC

In this work, we present an efficient implementation of polynomial filtering methods in PARSEC, a real-space pseudopotential based Kohn–Sham density functional theory solver. The implementation described here improves upon a Chebyshev-filtered subspace iteration algorithm used in the previous version of PARSEC. We present a hybrid polynomial filtering scheme that combines Chebyshev-filtered subspace iteration and a spectrum slicing method that partitions the spectrum into several spectral slices and uses bandpass-filtered subspace iteration to compute approximate eigenpairs within each interior slice simultaneously. We describe a procedure to partition a spectrum and construct polynomial filters. We also discuss a number of practical issues such as the use of appropriate data layouts for carrying out the computation on a two-dimensional process grid and how to achieve good load balance by allocating an appropriate number of process groups to each spectral slice. Numerical examples are presented to demonstrate the effectiveness of the hybrid polynomial filtering method as well as the superior parallel scalability of spectrum slicing in comparison to that of Chebyshev-filtered subspace iteration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A multi-center quadrature scheme for the molecular continuum

We report a common way to evaluate electronic integrals for polyatomic molecules is to use Becke’s partitioning scheme (Becke and Chem, 1988) in conjunction with overlapping grids centered at each atomic site. The Becke scheme was designed for integrands that fall off rapidly at large distances, such as those approximating bound electronic states. When applied to states in the electronic continuum, however, Becke scheme exhibits slow convergence and it is highly redundant. Here, we present a modified version of Becke scheme that is applicable to functions of the electronic continuum, such as those involved in molecular photoionization and electron–molecule scattering, and which ensures convergence and efficiency comparable to those realized in the calculation of bound states. In this modified scheme, the atomic weights already present in Becke’s partition are smoothly switched off within a range of few bond lengths from their respective nuclei, and complemented by an asymptotically unitary weight. The atomic integrals are evaluated on small spherical grids, centered on each atom, with size commensurate to the support of the corresponding atomic weight. The residual integral of the interstitial and long-range region is evaluated with a central master grid. The accuracy of the method is demonstrated by evaluating integrals involving integrands containing Gaussian Type Orbitals and Yukawa potentials, on the atomic sites, as well as spherical Bessel functions centered on the master grid. These functions are representative of those encountered in realistic electron-scattering and photoionization calculations in polyatomic molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Accelerated impurity solver for DMFT and its diagrammatic extensions

Here, we present ComCTQMC, a GPU accelerated quantum impurity solver. It uses the continuous-time quantum Monte Carlo (CTQMC) algorithm wherein the partition function is expanded in terms of the hybridisation function (CT-HYB). ComCTQMC supports both partition and worm-space measurements, and it uses improved estimators and the reduced density matrix to improve observable measurements whenever possible. ComCTQMC efficiently measures all one and two-particle Green's functions, all static observables which commute with the local Hamiltonian, and the occupation of each impurity orbital. ComCTQMC can solve complex-valued impurities with crystal fields that are hybridized to both fermionic and bosonic baths. Most importantly, ComCTQMC utilizes graphical processing units (GPUs), if available, to dramatically accelerate the CTQMC algorithm when the Hilbert space is sufficiently large. We demonstrate acceleration by a factor of over 600 (100) in a simulation of δ-Pu at 600 K with (without) crystal fields. In easier problems, the GPU offers less impressive acceleration or even decelerates the CTQMC. Here we describe the theory, algorithms, and structure used by ComCTQMC in order to achieve this set of features and level of acceleration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

INGRID: An interactive grid generator for 2D edge plasma modeling

A fusion boundary-plasma domain is defined by axisymmetric magnetic surfaces where the geometry is often complicated by the presence of one or more X-points; and modeling boundary plasmas usually relies on computational grids that account for the magnetic field geometry. The new grid generator INGRID (Interactive Grid Generator) presented in this work is a Python-based code for calculating grids for fusion boundary plasma modeling, for a variety of configurations with one or two X-points in the domain. INGRID first performs partitioning over the domain consisting of a small number of patches conforming to the magnetic field and wall geometry; then it generates a subgrid on each of the patches and joins them into a global grid. This domain partitioning strategy makes possible a uniform treatment of various configurations with one or two X-points in the domain. This includes single-null, double-null, and other configurations with two X-points in the domain. The INGRID design allows generating grids either interactively, via a parameter-file driven GUI, or using a non-interactive script-controlled workflow. Results of testing demonstrate that INGRID is a flexible, robust, and user-friendly grid-generation tool for fusion boundary-plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗