Development of MPACT for Full-Core Simulations of MAGNOX Gas-Cooled Nuclear Reactors
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Multiphysics simulations for nuclear reactor analysis are usually performed by resorting to operator splitting and fixed point iterations between single-physics solvers. This enables the separate solution of each physics, such as neutronics, fuel performance, and thermal hydraulics, on meshes tailored to the requirements of the respective numerical discretizations of the equations. As the equations are coupled, several fields must be transferred between single-physics solves. Projecting fields between meshes while preserving order of accuracy, conservation properties, and mapping non-overlapping geometries is a complex endeavor. This conference paper will present the transfers as implemented in MOOSE, which can handle arbitrary meshes, arbitrary mappings, conservation of integral quantities, and are made to scale with distributed simulations on both ends of the transfers. Their adequacy for advanced nuclear reactor multiphysics coupling is shown through examples and numerical studies.
Dynamical response of water exposed to x-rays is of utmost importance in a wealth of science areas. We exposed isolated water isotopologues to short x-ray pulses from a free-electron laser and detected momenta of all produced ions in coincidence. By combining experimental results and theoretical modeling, we identify significant structural dynamics with characteristic isotope effects in H 2 O 2+ , D 2 O 2+ , and HDO 2+ , such as asymmetric bond elongation and bond-angle opening, leading to two-body or three-body fragmentation on a timescale of a few femtoseconds. A method to disentangle the sequences of events taking place upon the consecutive absorption of two x-ray photons is described. The obtained deep look into structural properties and dynamics of dissociating water isotopologues provides essential insights into the underlying mechanisms.
We have studied two regioisomeric terrylenediimide (TDI) dimers in which the 1-positions of two TDIs are linked via 1,3- or 1,4-phenylene spacers, mTDI 2 and pTDI, respectively. The nature and the dynamics of the multiexciton state are tuned by altering the through-bond electronic couplings in the ground and excited states and by changing the solvent environment. Our results show that controlling the electronic coupling between the two chromophores by an appropriate choice of linker can result in independent triplet state formation, even though the initial correlated triplet pair state is confined to a dimer. Moreover, even in polar solvents, if the electronic coupling is strong, the correlated triplet pair state is observed prior to symmetry-breaking charge separation. These results point out the close relationship between the singlet, correlated triplet pair, and charge transfer states in molecular dimers.
Two schemes for coupling gyrokinetic simulations of microturbulence in tokamaks are proposed. The first scheme is based on an additive Schwarz domain decomposition. In this work, we show that, because the goal of turbulence is long-time averages of the dynamics rather than temporal accuracy, the iteration to self-consistency across domains, which is typically required by Schwarz schemes, can be avoided, thereby accelerating the computation. Second, we propose a coupling scheme that relies entirely on the addition of source terms, leaving the boundary conditions arbitrary. The practical motivations for such a scheme are discussed, and forms of the source terms that ensure consistency and stability are derived. The schemes are tested on a nonlinear, one-dimensional model problem, and the first scheme is further tested on the Hasegawa–Wakatani model.
This paper provides modeling evidence for a new variant of microturbulence instability that regulates the electron thermal fluxes in the transport barrier region of DIII-D "high β p " discharges. This scenario is being investigated as a possible scenario for future steady-state reactors, and thus it is essential to understand the physics processes which govern its performance in order to confidently extrapolate to future devices. Additional analysis examining the extrapolation of this scenario to ITER shows that while the new microtearing mode identified here is unlikely to be unstable in that regime (due to increased temperatures and reduced collisionality), the large shift in magnetic axis (Shafranov shift) that arises in these plasmas may still favor stabilization of common ballooning instabilities which peak on the low field side, but destabilize electrostatic slab modes on the high field side.
Femtosecond X-ray pump X-ray probe experiments are currently possible at free electron lasers such as the linac coherent light source (LCLS), which opens new opportunities for studying solvated transition metal-ligand complexes. In order to make the most effective use of these kinds of experiments, it is necessary to determine which chemical properties an X-ray pump pulse will measure. In this paper, we have combined electron cascade calculations and excited-state time-dependent density functional theory calculations to predict the initial state prepared by an X-ray pump and the subsequent X-ray probe spectra at the Fe K-edge in the solvated model transition metal complex, K 4 Fe II (CN) 6 . We find several key spectral features that report on the ligand-field splitting, solvent-solute interactions, and the 3p and 3d electron interactions. We then show how these features could be measured in an experiment.
We have measured, analyzed and simulated the ground state valence photoelectron (PES), X-ray absorption (XAS), X-ray photoelectron (XPS) as well as normal and resonant Auger-Meitner electron (AES) spectra of oxazole at the carbon, oxygen and nitrogen K-edge in order to understand its electronic structure. Experimental data is compared to theoretical calculations performed at the coupled cluster, restricted active space perturbation theory to second order and time-dependent density functional levels of theory. We demonstrate that both N and O K-edge XAS are sensitive to the amount of dynamical electron correlation included in the theoretical description, and that for a complete description of XPS, additional orbital correlation and orbital relaxation effects need to be considered. The normal AES are dominated by a singlet excitation channel and well described by theory. Furthermore, the resonant AES, however, are more complicated. While the participator decay channels, dominating at higher kinetic energies, are well described by coupled cluster theory, spectator channels can only be described satisfactorily using a method that combines restricted active space perturbation theory to second order for the bound part and a one-center approximation for the continuum.
A model of an all-plasmonic switch based on a graphene waveguide integrated with a stub nanoresonator loaded with a quantum dot is considered. Using the full-wave numerical simulation, it is shown that successive on/off switching of the control surface plasmon – polariton wave at the input of the waveguide leads to a change in the phase of the signal plasmon – polariton wave in the nanoresonator by π. This causes reversible switching of the device operation regime – from almost complete blocking to stable transmission of the signal plasmon – polariton wave through the waveguide. The effect is implemented at rates of ∼0.5 THz for infrared electromagnetic waves localised in a device with a switching region of 40 × 20 nm. (paper)
Photosynthetic organisms rely on sophisticated photoprotective mechanisms to prevent oxidative damage under high or fluctuating solar illumination. Cyanobacteria, which have evolved a unique, water-soluble light-harvesting complex—the phycobilisome—achieve photoprotection through a photoactivatable quencher called the Orange Carotenoid Protein (OCP). Phycobilisomes are highly symmetric and modular, formed by hierarchical assembly of conserved subunits into diverse geometries ranging from simple bundles to elaborate fan- or bouquet-like macromolecular architectures. Although OCP is known to provide photoprotection across species of cyanobacteria with different phycobilisome structures, it is not known whether or how these structural variations relate to changes in the photoprotective function of OCP. For example, OCP was recently discovered to bind as a dimer at two specific instances of an abundant structural motif on the tricylindrical phycobilisome of Synechocystis sp. PCC 6803, yet these sites are sterically inaccessible on a more common pentacylindrical phycobilisome ( Anabaena sp. PCC 7120). To understand how structural modularity and binding specificity contribute to conservation of OCP binding sites and function across different phycobilisome architectures, here we compare experimentally measured photophysical states accessible to these prototypical tricylindrical and pentacylindrical phycobilisomes, with and without OCP, at the single-molecule level. Together with Monte Carlo simulations of exciton transfer in OCP-quenched phycobilisomes, our results suggest that OCP binds at distinct and specific sites in each type of phycobilisome, yet provides nearly identical quenching strength to both phycobilisomes. Our findings highlight the utility of modular phycobilisome structures in balancing robust conservation of photoprotective function with adaptability of site-specific binding across species.
The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.
The CASL reactor simulation package VERA has been adapted to provide high-fidelity simulation capabilities for modeling source range detector response during subcritical reactor configurations. New features include the activation and shuffling of secondary-source assemblies, use of burned fuel neutron emission data from the ORIGEN depletion solver to the MPACT deterministic neutron transport solver, allowance of user-defined sources in MPACT based on material composition, ability to solve the subcritical source-driven system with neutron multiplication using the MPACT diffusion solver, and transfer of the calculated fission source from MPACT to the continuous-energy Monte Carlo solver Shift for final detector response evaluation using the CADIS methodology for variance reduction. These new capabilities were benchmarked against Watts Bar Unit 1 plant operating data for the first few fuel loading steps and were found to demonstrate excellent agreement with the measured data.
Not Available
Additive manufacturing (AM) methods are currently being explored for applications in nuclear reactors to make advanced reactors more efficient, safe, and reliable. The Transformational Challenge Reactor (TCR) program has explored AM for nuclear by designing a high-temperature gas reactor (HTGR) using an AM silicon carbide fuel form with uranium nitride–tristructural isotropic fuel. In this study, we detail the design process for the TCR fuel form’s coolant channels using computational fluid dynamics models with conjugate heat transfer. Additionally, this work discusses how these models were interfaced with other design teams, project milestones, and the agile design method used to mature the reactor design. The methodology deployed was able to create a channel design with lower maximum fuel temperatures and thermal stresses in the fuel form over traditional channel designs that can be manufactured subtractively. These results were achieved with only small manufacturing penalties. Results are discussed and presented on lessons learned for designing AM components for nuclear reactors. Finally, areas of opportunity are discussed for advanced design tools to further automate design activities and optimize reactors with fewer built-in assumptions.
Not Available
Here we investigate using deep learning, a type of machine-learning algorithm employing multiple layers of artificial neurons, for the mathematical representation of multigroup cross sections for use in the Griffin reactor multiphysics code for two-step deterministic neutronics calculations. A three-dimensional fuel element typical of a high-temperature gas reactor as well as a two-dimensional sodium-cooled fast reactor lattice are modeled using the Serpent Monte Carlo code, and multigroup macroscopic cross sections are generated for various state parameters to produce a training data set and a separate validation data set. A fully connected, feedforward neural network is trained using the open-source PyTorch machine-learning framework, and its accuracy is compared against the standard piecewise linear interpolation model. Additionally, we provide in this work a generic technique for propagating the cross-section model errors up to the k eff using sensitivity coefficients with the first-order uncertainty propagation rule. Quantifying the eigenvalue error due to the cross-section regression errors is especially practical for appropriately selecting the mathematical representation of the cross sections. We demonstrate that the artificial neural network model produces lower errors and therefore enables better accuracy relative to the piecewise linear model when the cross sections exhibit nonlinear dependencies; especially when a coarse grid is employed, where the errors can be halved by the artificial neural network. However, for linearly dependent multigroup cross sections as found for the sodium-cooled fast reactor case, a simpler linear regression outperforms deeper networks.