Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear matrix elements”

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 289 records · Page 16

Deformation Mechanism Transition in Additively Manufactured Compositionally Graded Fe-Base Alloys

Microstructure-dependent deformation and fracture behavior was investigated for an additively manufactured compositionally graded alloy (CGA) printed using the laser-directed energy deposition (L-DED) method to explore an alternative approach for dissimilar metal joints in nuclear energy systems. The electron backscatter diffraction (EBSD) maps from scanning electron microscopy (SEM) display a clear microstructural transition with decreasing austenite-forming elements (Ni and Mn), from an austenite (γ) dominant structure, to a complex composite structure containing ferrite (α), martensite (α') and retained austenite, and then to a fully ferritic structure. EBSD data were recorded in situ during tensile testing in SEM, and the evolution of the deformation mechanism and microstructure was characterized using Kikuchi diffraction pattern analysis. Complementary analysis for high-resolution features was also performed using scanning transmission electron microscopy (STEM). The Ni/Mn-rich austenite-dominant microstructures showed a complex deformation mechanism of two-step martensitic transformation (γ→ε→α'), whereas the minor austenite phase retained in the ferrite and/or martensite matrix showed a single transformation route (γ→α'). Ordinary dislocation glide and twinning via partial dislocation glide were observed in the austenite deformation. Meanwhile, the ferrite and martensite grains deformed mainly by ordinary dislocation slips and grain rotation. Furthermore, static tensile fracture was also highly dependent on local composition and phase constituents.

36 MATERIALS SCIENCE↗

Exponential time differencing scheme for mass transport and depletion in molten salt reactors

This work extends the capability previously shown for addressing the problem of computing depletion and mass transport calculations in molten salt reactors (MSRs) by calculating matrix exponentials. Additional algorithms are implemented to compute the matrix exponential and the action of the matrix exponential on a matrix. These algorithms include two methods based on the Pade approximation, a Taylor series method, and three methods based on Cauchy's integral formula. In addition to the added matrix exponential solvers, a variable-order total variation diminishing scheme is applied to the convective flux approximation to provide enhanced accuracy. Finally, a simplified MSR problem is shown for each of the exponential time differencing solvers along with classical backwards differencing integrators. The results show excellent convergence for exponential time differencing methods. Computation time is a key element for selecting the optimal solver in these problems, and this work shows that Pade and Cauchy-based solvers may provided the fastest and most accurate solutions. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

147 Nd Quantification Using HSCCC-Purified Samples

Quantifying the fission product 147 Nd in nuclear debris samples is an important component of post-detonation nuclear forensics. The most accurate quantifications are obtained when Nd is purified from all other fission products, actinides, activation products, and environmental matrix contained within the debris. In this study, a recently developed method for Nd purification was tested, purifying 147 Nd from solutions of mixed fission products using high-speed counter-current chromatography (HSCCC). Importantly, the new method allowed for faster elution of Nd from the column as compared with established high performance liquid chromatography (HPLC) methods, and resulted in accurate/precise 147 Nd quantification by gamma-ray spectrometry. While the up-front equipment costs associated with HSCCC may be higher, its operational costs are on par with those of HPLC (solvents, extractants, power). Gas-flow proportional beta decay counting revealed contamination from the nearest neighbor lanthanide 143 Pr (a gamma-silent radioisotope) in the HSCCC-purified samples, but the activity contribution from 147 Nd could still be quantified. Remarkably consistent elution profiles were observed for the HSCCC method, spanning rare earth element (REE) loadings of more than 10 orders of magnitude (tracer to mmol quantities). In conclusion, the reliability and speed of the new method suggest utility for the rapid separation and quantification of 147 Nd in unknown samples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Converged ab initio calculations of heavy nuclei

We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of E 3max =28, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of 132 Sn converged to within a few MeV with respect to the E 3max truncation. In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum J rel > 9/2 to the ground-state energy is dominated by a basis-truncation artifact, which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of 132 Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Galactic oscillator symmetry

Riemann ellipsoids model rotating galaxies when the galactic velocity field is a linear function of the Cartesian coordinates of the galactic masses. In nuclear physics, the kinetic energy in the linear velocity field approximation is known as the collective kinetic energy. But, the linear approximation neglects intrinsic degrees of freedom associated with nonlinear velocity fields. To remove this limitation, the theory of symplectic dynamical symmetry is developed for classical systems. A classical phase space for a self-gravitating symplectic system is a co-adjoint orbit of the noncompact group SP(3,R). The degenerate co-adjoint orbit is the 12 dimensional homogeneous space Sp(3,R)/U(3), where the maximal compact subgroup U(3) is the symmetry group of the harmonic oscillator. The Hamiltonian equations of motion on each orbit form a Lax system X = (X,F), where X and F are elements of the symplectic Lie algebra. The elements of the matrix X are the generators of the symplectic Lie algebra, viz., the one-body collective quadratic functions of the positions and momenta of the galactic masses. The matrix F is composed from the self-gravitating potential energy, the angular velocity, and the hydostatic pressure. Solutions to the hamiltonian dynamical system on Sp(3,R)/U(3) are given by symplectic isospectral deformations. The Casimirs of Sp(3,R), equal to the traces of powers of X, are conserved quantities.

Rosensteel, George↗

Efficient and Flexible Sensitivity Matrix Computation for Adaptive Electrical Capacitance Volume Tomography

Electrical capacitance tomography is a widely used sensor modality for flow imaging in many industrial settings. Adaptive Electrical Capacitance Volume Tomography (AECVT) extends the capabilities of traditional ECT by enabling direct volumetric imaging and an improved resolution. Construction of the sensitivity matrix is a necessary step to obtain flow images. This step requires computation of the electric field inside the sensing domain, which is done via a typical field solver such as the finite element method. In this work, we present an efficient and flexible method to construct the sensitivity matrix for Adaptive Electrical Capacitance Volume Tomography (AECVT) based on individual electrode segment excitations and their judicious combination to form desired matrix elements. We illustrate how the proposed method yields the same sensitivity matrix as the traditional method but at a much lower computational cost. Once all segment contributions are obtained, we also indicate how the proposed method, unlike the traditional approach, can generate the sensitivity matrix on demand for an arbitrary combination of synthetic electrodes and obviating the need for any additional field computations. Finally, we present image reconstruction results for two different experimental scenarios where the mutual capacitance data and the corresponding sensitivity vectors are obtained through the proposed measurement combination scheme.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Construction of a Code Verification Matrix for Heat Conduction With Finite Element Code Applications

When establishing the pedigree of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model. During this process, numerical results on various meshes are systematically compared to a reference analytic solution. The selection of analytic solutions can be a laborious process, as it is difficult to establish adequate code confidence without performing redundant work. In this report we address this issue by applying a physics-based process that establishes a set of reference problems. In this process, code simulation options are categorized and systematically tested, which ensures that gaps in testing are easily identified and addressed. The resulting problems are primarily intended for code verification analysis but may also be useful for comparison to other simulation codes, troubleshooting activities, or training exercises. The process is used to select fifteen code verification problems relevant for the one-dimensional steady-state heat conduction equation. These problems are applicable to a wide variety of simulation tools, but, in this work, a demonstration is performed using the finite element-based nuclear fuel performance code BISON. Convergence to the analytic solution at the theoretical rate is quantified for a selection of the problems, which establishes a baseline pedigree for the code. Not only can this standard set of conduction solutions be used for verification of other codes, but also the physics-based process for selecting problems can be utilized to quantify and expand testing for any simulation tool.

42 ENGINEERING↗

Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}$ 2 from lattice QCD at physical quark masses

We present a lattice quantum chromodynamics calculation of the 𝑥-dependent isovector quark helicity parton distribution function (PDF) of the proton in the large momentum effective theory (LaMET) framework. Through operator product expansion (OPE) we also extract the $\tilde{d}$ 2 moment of the twist-3 PDF 𝑔 𝑇 ⁡(𝑥) for the first time in the $\overline{MS}$ scheme, which is proportional to the average color Lorentz force experienced by the quark in the proton. This calculation is performed on a lattice of spacing 𝑎 =0.076 fm at physical quark masses. The quasi-PDF matrix elements are measured in proton states boosted to momenta 𝑃 𝑧 ={0,0.25,1.02,1.53} GeV. We first extract the lowest few helicity PDF moments from the renormalization-group (RG) invariant ratios of the matrix elements with OPE. Combined with the matrix elements relevant for 𝑔 𝑇 ⁡(𝑥), we obtain $\tilde{d}$$^{u-d}_2$⁡(2 GeV) =0.0024⁢(46) at next-to-leading order in $\overline{MS}$. Then, the helicity quasi-PDF matrix elements are renormalized in the hybrid scheme with linear renormalon resummation and Fourier transformed to the 𝑥-space after an asymptotic extrapolation. The quasi-PDF is perturbatively matched to the $\overline{MS}$ PDF with RG and threshold resummations at next-to-leading power and next-to-next-to-leading logarithmic accuracies. After resummations, we determine the PDF in the region 𝑥 ∈[0.25,0.75]. The end-point regions are then parametrized, combined with the LaMET prediction at moderate 𝑥, and fitted to the short-distance matrix elements in coordinate space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Characteristics Approach to the Finite Element Method

Herein, we present a new method for solving the linear Boltzmann transport equation. Two commonly used and well-understood methods for solving partial differential equations are the method of characteristics (MOC) and the finite element method (FEM). We propose a new method that combines the fundamental concept of the FEM with the analytic solution from the MOC to obtain coefficients for the FEM basis function expansion. Traditionally, coefficients for the FEM basis function expansion are obtained via matrix inversion. Instead, we solve for the coefficients with the MOC and represent the underlying fields with the basis function expansion using these coefficients. We provide a convergence study for our method with results from two sets of FEM basis functions: Gauss-Legendre and Gauss-Lobatto sets. We also compare two different variations of our method categorized as short characteristics and intermediate characteristics.

42 ENGINEERING↗

Comparison of Fission Product Release Predictions using PARFUME and BISON with Results from the AGR-3/4 Irradiation Experiment

The PARFUME (PARticle Fuel ModEl) fuel performance modeling code and the BISON nuclear fuel performance application built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library were used to predict the fission product release from tristructural isotropic (TRISO) coated fuel particles and compacts during the third and fourth irradiation experiment of the Advanced Gas Reactor (AGR-3/4) Fuel Development and Qualification Program. The fuel performance modeling codes PARFUME and BISON modeled the AGR-3/4 irradiation experiment using the fuel compact time-averaged volume averaged (TAVA) daily temperatures for a total irradiation duration of 369.1 effective full power days (EFPD) to predict the release fraction of the fission product silver (Ag-110m) from a representative TRISO-coated fuel particle from AGR-3/4 compacts. Post-irradiation examination (PIE) measurements provided data on the release of these fission products in the compacts outside of the silicon carbide (SIC) layer. The PARFUME and BISON results were then compared to the silver release measured from compact gamma scanning. The results showed good agreement between PARFUME and BISON but both codes under-predicted the silver release fraction for all the compacts. In addition, BISON was used to model and predict the fission product concentration radial profile outside of the compacts in capsules’ inner and outer rings. These rings were either comprised of matrix and/or structural graphite. To obtain the concentration profiles of silver, cesium, and strontium, a sorption isotherm model was developed in BISON to capture the effects of fission product transport across the gaps between the concentric rings. The general shape of the concentration radial profiles as calculated by BISON were similar in the inner ring (IR) but varied in the outer ring (OR) depending on the fission product of interest or capsule temperature. Using this methodology and model, BISON now has the capability to aid in developing new fission product diffusion coefficients for matrix or structural graphite materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Carbide-based fuel assembly for thermal propulsion applications

Carbide-based fuel assembly includes outer structural member of ceramic matrix composite material, the interior surface of which is lined in higher temperature regions with an insulation layer of porous refractory ceramic material. Continuous insulation layer extends the length of the fuel assembly or separate insulation layer sections have a thickness increasing step-wise along the length of the fuel assembly from upper (inlet) section towards bottom (outlet) section. A fuel element positioned inward of the insulation layer and between support meshes has a fuel composition including HALEU and the form of a plurality of individual elongated fuel bodies or one or more fuel monolith bodies containing coolant flow channels. Fuel assemblies are distributively arranged in a moderator block, with upper end of the outer structural member attached to an inlet for propellant and lower end of the outer structural member operatively interfaced with a nozzle forming a nuclear thermal propulsion reactor.

Barringer, Eric A.↗

Carbide-based fuel assembly for thermal propulsion applications

Carbide-based fuel assembly includes outer structural member of ceramic matrix composite material (e.g., SiC—SiC composite), insulation layer of porous refractory ceramic material (e.g., zirconium carbide with open-cell foam structure or fibrous zirconium carbide), and interior structural member of refractory ceramic-graphite composite material (e.g., zirconium carbide-graphite or niobium carbide-graphite). Spacer structures between various layers provide a defined and controlled spacing relationship. A fuel element bundle positioned between support meshes includes a plurality of distributively arranged fuel elements or a solid, unitary fuel element with coolant channels, each having a fuel composition including high assay, low enriched uranium (HALEU). Fuel assemblies are distributively arranged in a moderator block and the upper end of the outer structural member is attached to a metallic inlet tube for hydrogen propellant and the lower end of the outer structural member is interfaced with a support plate, forming a nuclear thermal propulsion reactor.

Barringer, Eric A.↗

A thermodynamics-based damage model for the non-linear mechanical behavior of SiC/SiC ceramic matrix composites in irradiation and thermal environments

A damage model is developed and validated with experimental data for the non-linear mechanical behavior of SiC/SiC composite materials in nuclear applications. Cyclic thermal and mechanical loading associated with neutron irradiation effects of these composites leads to wide-spread and progressive micro-cracking that leads to loss of thermal conductivity and further enhancement of thermo-mechanical damage. A physics-based model of wide-spread micro-cracking is developed within the thermodynamic framework of continuum damage mechanics. Evolution equations for damage parameters that describe the growth of continuum damage are developed, where the material variables are obtained from experiments. The model novelty is in coupling mechanical, thermal, and irradiation damage through a consistent thermodynamic framework, including loss of thermal conductivity due to the evolution of mechanically induced micro-cracks. A number of thermo-mechanical experiments were conducted to confirm model assumptions. The model is shown to be validated with out-of-pile experiments, and then implemented using commercial finite element code COMSOL to the fuel cladding problem with normal and high radiation dose cases.

Materials Science↗

A Review of Advanced Test Reactor Fuel and Assessment of Its Compatibility with the ZIRCEX Chlorination Process

Advanced Test Reactor (ATR) fuel has been identified as a resource for high-assay low-enriched uranium (HALEU) production. A survey was performed on the published literature describing ATR fuel. The geometry of the fuel is complex; different parts of the fuel compact experience differing neutron flux and burnup. The literature is sparse, and access is controlled. Therefore, fundamental studies of fuel reprocessing must use a model fuel that represents the main chemical and structural features. Advanced chlorination, or chlorination with sulfur-chlorine bearing reagents is being investigated as way to separate the fuel from metal matrix alloys. A UAl x alloy will be fabricated with x = 3, 4, and 5. The potential chlorination of individual UAl x intermetallics will be assessed in the advanced chlorination process of Al-8001 and Al-6061 as well as a representative mixture. Initial studies will track the alloying elements of the Al, which are Si, Fe, Cu, Mn, Mg, Cr, Zn, and Ti, in addition to the U itself. Further studies will include fission product simulants. Because advanced chlorination solvents include sulfur, the chemistry of sulfur with major and minor constituents will also be investigated. The experimental work accompanied by neutronic calculations will allow the assessment of the feasibility of advanced chlorination to separate aluminum from uranium. If bench-scale testing appears promising, then small-scale tests in shielded facilities with irradiated cladding, lightly irradiated fuel, and spent nuclear fuel are recommended to track the complete inventory of fissile actinides, fission product impurities, and reagent solids and liquids.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Low and moderate x gluon contribution to exclusive Compton scattering processes

We revisit the high energy semi-classical description of the exclusive processes DVCS, TCS, and Double DVCS by explicitly keeping track of the Feynman x dependence in both the hard and the hadronic matrix elements. This is achieved by a modification of the standard shock wave approximation to derive the effective Feynman rules, which leads to a generic expression on which we then perform a partial twist expansion to get rid of quantities suppressed by the proper physical scales. We obtain a compact factorized master formula that can be used to investigate the Bjorken limit at leading twist. In particular, we recover the full one-loop result in the collinear limit for pure gluon exchange with the target. Finally, we discuss the subtleties in taking the simultaneous collinear and small x limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Candidate Core Designs for the Transformational Challenge Reactor

Early cycle activities under the Transformational Challenge Reactor (TCR) program focused on analyzing and maturing four reactor core design concepts: two fast-spectrum systems and two thermal-spectrum systems. A rapid, iterative approach has been implemented through which designs can be modified and analyzed and subcomponents can be manufactured in parallel over time frames of weeks rather than months or years. To meet key program initiatives (e.g., timeline, material use), several constraints—including fissile material availability (less than 250 kg of HALEU), component availabilities, materials compatibility, and additive manufacturing capabilities—were factored into the design effort, yielding small (less than one cubic meter in volume) cores with near-term viability. The fast-spectrum designs did not meet the fissile material constraint, so the thermal-spectrum systems became the primary design focus. Since significant progress has been made on advanced moderator materials (YH x ) under the TCR program, gas-cooled thermal-spectrum systems using less than 250 kg of HALEU that occupy less than 1 m 3 are now feasible. The designs for two of these systems have been evolved and matured. In both thermal-spectrum design concepts, bidirectional coolant flow is used. Coolant flows down through YH x moderator elements and is reversed in a bottom manifold and core support structure, and then flows up though or around the fuel elements. The main difference between the two thermal-spectrum design concepts is the fuel elements—one uses traditional UO 2 ceramic fuel, and the other uses UN-bearing TRISO fuel particles embedded inside a SiC matrix. Finally, core neutronics and thermal performance for these systems are assessed and summarized herein.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Isospin mixing and the cubic isobaric multiplet mass equation in the lowest T=2, A=32 quintet

The isobaric multiplet mass equation (IMME) is known to break down in the first T=2,A=32 isospin quintet. In this work we combine high-resolution experimental data with state-of-the-art shell-model calculations to investigate isospin mixing as a possible cause for this violation. The experimental data are used to validate isospin-mixing matrix elements calculated with newly developed shell-model Hamiltonians. Our analysis shows that isospin mixing with nonanalog T=1 states contributes to the IMME breakdown, making the requirement of an anomalous cubic term inevitable for the multiplet.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Postirradiation characterization of palladium as an additive for fuel cladding chemical interaction mitigation in metallic fuel

This work describes the microstructural and elemental characterization of irradiated metallic fuels containing palladium as an additive. The use of additives has been proposed to control Fuel-Cladding Chemical Interaction (FCCI) and thus to promote higher fuel utilization (i.e., higher burnup). In this work, Pd has been investigated as a potential additive to metallic fuel to bind lanthanides, impeding their migration and attack on the cladding. The influence of Pd on the microstructure, chemistry and performance of metallic fuel has been characterized via scanning electron microscopy for two metallic fuel designs—namely, annular and solid fuel. Pd was observed to play an important role in the chemistry of the fuel. Indeed, the addition of Pd leads to the formation of new phases. Pd was detected to combine not only with the lanthanides, as intended, but also with Zr, a main element of the fuel matrix. While Pd proved to be effective in preventing lanthanide migration and their attack on the cladding, the Pd-Zr compound may potentially lead to other unexpected fuel-performance issues, such as the formation of low-melting point phases and increased unalloyed U available for FCCI interaction with Fe in the cladding. Even the increase of Zr to 13wt%. did not completely mitigate this adverse phenomenon generated by the Pd-Zr interaction. Furthermore, the efficacy of using this additive needs further investigation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗