Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core simulation”

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,351 records · Page 75

Online randomized interpolative decomposition with a posteriori error estimator for temporal PDE data reduction

Traditional low-rank approximation is a powerful tool for compressing large data matrices that arise in simulations of partial differential equations (PDEs), but suffers from high computational cost and requires several passes over the PDE data. The compressed data may also lack interpretability thus making it difficult to identify feature patterns from the original data. Here, to address these issues, we present an online randomized algorithm to compute the interpolative decomposition (ID) of large-scale data matrices in situ. Compared to previous randomized IDs that used the QR decomposition to determine the column basis, we adopt a streaming ridge leverage score-based column subset selection algorithm that dynamically selects proper basis columns from the data and thus avoids an extra pass over the data to compute the coefficient matrix of the ID. In particular, we adopt a single-pass error estimator based on the non-adaptive Hutch++ algorithm to provide real-time error approximation for determining the best coefficients. As a result, our approach only needs a single pass over the original data and thus is suitable for large and high-dimensional matrices stored outside of core memory or generated in PDE simulations. A strategy to improve the accuracy of the reconstructed data gradient, when desired, within the ID framework is also presented. We provide numerical experiments on turbulent channel flow and ignition simulations, and on the NSTX Gas Puff Image dataset, comparing our algorithm with the offline ID algorithm to demonstrate its utility in real-world applications.

Column subset selection↗

Comparing XRISM Cluster Velocity Dispersions with Predictions from Cosmological Simulations: Are Feedback Models Too Ejective?

The dynamics of the intracluster medium (ICM), the hot plasma that fills galaxy clusters, are shaped by gravity-driven cluster mergers and feedback from supermassive black holes (SMBHs) in the cluster cores. XRISM measurements of ICM velocities in several clusters offer insights into these processes. We compare XRISM measurements for nine galaxy clusters (Virgo, Perseus, Centaurus, Hydra-A, PKS 0745–19, A2029, Coma, A2319, and Ophiuchus) with predictions from three state-of-the-art cosmological simulation suites, TNG-Cluster, the Three Hundred Project GADGET-X, and GIZMO-SIMBA, that employ different models of feedback. In cool cores, XRISM reveals systematically lower velocity dispersions than the simulations predict, with all 10 measurements below the median simulated values by a factor of 1.5–1.7 on average and all falling within the bottom 10% of the predicted distributions. The observed kinetic-to-total pressure ratio is also lower, with a median value of 2.2%, compared to the predicted 5.0%–6.5% for the three simulations. Outside the cool cores and in non-cool-core (NCC) clusters, simulations show better agreement with XRISM measurements, except for the outskirts of the relaxed, cool-core cluster A2029, which exhibits an exceptionally low kinetic pressure support (<1%), with none of the simulated systems in either of the three suites reaching such low levels. The NCC Coma and A2319 exhibit dispersions at the lower end but within the simulated spread. Our comparison suggests that the three numerical models may overestimate the kinetic effects of SMBH feedback in cluster cores. Additional XRISM observations of NCC clusters will clarify if there is a systematic tension in the gravity-dominated regime as well.

galaxy clusters↗

Systematic and simulation-free coarse graining of homopolymer melts: A structure-based study

We propose a systematic and simulation-free strategy for coarse graining of homopolymer melts, where each chain of N{sub m} monomers is uniformly divided into N segments, with the spatial position of each segment corresponding to the center-of-mass of its monomers. We use integral-equation theories suitable for the study of equilibrium properties of polymers, instead of many-chain molecular simulations, to obtain the structural and thermodynamic properties of both original and coarse-grained (CG) systems, and quantitatively examine how the effective pair potentials between CG segments and the thermodynamic properties of CG systems vary with N. Our systematic and simulation-free strategy is much faster than those using many-chain simulations, thus effectively solving the transferability problem in coarse graining, and provides the quantitative basis for choosing the appropriate N-values. It also avoids the problems caused by finite-size effects and statistical uncertainties in many-chain simulations. Taking the simple hard-core Gaussian thread model [K. S. Schweizer and J. G. Curro, Chem. Phys. 149, 105 (1990)] as the original system, we demonstrate our strategy applied to structure-based coarse graining, which is quite general and versatile, and compare in detail the various integral-equation theories and closures for coarse graining. Our numerical results show that the effective CG potentials for various N and closures can be collapsed approximately onto the same curve, and that structure-based coarse graining cannot give thermodynamic consistency between original and CG systems at any N < N{sub m}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-linear dynamics of compound sawteeth in tokamaks

Compound sawteeth is studied with the XTOR-2F code. Non-linear full 3D magnetohydrodynamic simulations show that the plasma hot core is radially displaced and rotates during the partial crash, but is not fully expelled out of the q = 1 surface. Partial crashes occur when the radius of the q = 1 surface exceeds a critical value, at fixed poloidal beta. This critical value depends on the plasma elongation. The partial crash time is larger than the collapse time of an ordinary sawtooth, likely due to a weaker diamagnetic stabilization. This suggests that partial crashes result from a competition between destabilizing effects such as the q = 1 radius and diamagnetic stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improvement of the edge rotation diagnostic spectrum analysis via simulation

The edge rotation diagnostic (ERD) system has been developed on the Joint Texas Experimental Tokamak to measure the edge toroidal rotation velocity by observing the shifted wavelength of carbon V (C V 227.09 nm). Since the measured spectrum is an integrated result along the viewing line from the plasma core to the edge, a method via simulation has been developed to analyze the ERD spectrum. With the necessary parameters such as C V radiation profile and the ion temperature profile, a local rotation profile at the normalized minor radius of 0.5-1 is obtained.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Status of NASA's Space Launch System

NASA's Space Launch System (SLS) continued to make significant progress in 2015 and 2016, completing hardware and testing that brings NASA closer to a new era of deep space exploration. Programmatically, SLS completed Critical Design Review (CDR) in 2015. A team of independent reviewers concluded that the vehicle design is technically and programmatically ready to move to Design Certification Review (DCR) and launch readiness in 2018. Just five years after program start, every major element has amassed development and flight hardware and completed key tests that will lead to an accelerated pace of manufacturing and testing in 2016 and 2017. Key to SLS' rapid progress has been the use of existing technologies adapted to the new launch vehicle. The existing fleet of RS-25 engines is undergoing adaptation tests to prove it can meet SLS requirements and environments with minimal change. The four-segment shuttle-era booster has been modified and updated with a fifth propellant segment, new insulation, and new avionics. The Interim Cryogenic Upper Stage is a modified version of an existing upper stage. The first Block I SLS configuration will launch a minimum of 70 metric tons (t) of payload to low Earth orbit (LEO). The vehicle architecture has a clear evolutionary path to more than 100t and, ultimately, to 130t. Among the program's major 2015-2016 accomplishments were two booster qualification hotfire tests, a series of RS-25 adaptation hotfire tests, manufacturing of most of the major components for both core stage test articles and first flight tank, delivery of the Pegasus core stage barge, and the upper stage simulator. Renovations to the B-2 test stand for stage green run testing was completed at NASA Stennis Space Center. This year will see the completion of welding for all qualification and flight EM-1 core stage components and testing of flight avionics, completion of core stage structural test stands, casting of the EM-1 solid rocket motors, additional testing of RS-25 engines and flight engine controllers This paper will discuss these and other technical and programmatic successes and challenges over the past year and provide a preview of work ahead before the first flight of this new capability.

Honeycutt, John↗

Verification of High-Fidelity Neutronics Code PROTEUS for C5G7 Benchmark Problems

Under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the 3D whole core transport code PROTEUS has been developed at Argonne National Laboratory to perform the detailed 3D transport calculation for various reactor types including irregular or complex geometric reactors. PROTEUS includes two distinct transport solvers: SN and MOC. Among them, the MOC solver based on the extruded geometry, thus named MOCEX, can produce the faithful 3-D transport solutions by employing the rigorous and robust 2D/3D MOC formulation. This formulation overcomes the drawbacks of a full 3D MOC method requiring enormous memory and computational resources and resolves the instability issues of the 2D/1D method as well. The solution of MOCEX converges to that of a Monte Carlo code as the discretization of angle, space and energy is refined, unlike the 2D/1D method. As a verification and validation (V and V) effort of PROTEUS-MOCEX, this paper presents verification test results using the well-known 3D C5G7 benchmark problems by solving and comparing eigenvalues and pin powers with MCNP solutions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thorium Fuel Options for Sustained Transuranic Burning in Pressurized Water Reactors - 12381

As described in companion papers, Westinghouse is proposing the adoption of a thorium-based fuel cycle to burn the transuranics (TRU) contained in the current Used Nuclear Fuel (UNF) and transition towards a less radio-toxic high level waste. A combination of both light water reactors (LWR) and fast reactors (FR) is envisaged for the task, with the emphasis initially posed on their TRU burning capability and eventually to their self-sufficiency. Given the many technical challenges and development times related to the deployment of TRU burners fast reactors, an interim solution making best use of the current resources to initiate burning the legacy TRU inventory while developing and testing some technologies of later use is desirable. In this perspective, a portion of the LWR fleet can be used to start burning the legacy TRUs using Th-based fuels compatible with the current plants and operational features. This analysis focuses on a typical 4-loop PWR, with 17x17 fuel assembly design and TRUs (or Pu) admixed with Th (similar to U-MOX fuel, but with Th instead of U). Global calculations of the core were represented with unit assembly simulations using the Linear Reactivity Model (LRM). Several assembly configurations have been developed to offer two options that can be attractive during the TRU transmutation campaign: maximization of the TRU transmutation rate and capability for TRU multi-recycling, to extend the option of TRU recycling in LWR until the FR is available. Homogeneous as well as heterogeneous assembly configurations have been developed with various recycling schemes (Pu recycle, TRU recycle, TRU and in-bred U recycle etc.). Oxide as well as nitride fuels have been examined. This enabled an assessment of the potential for burning and multi-recycling TRU in a Th-based fuel PWR to compare against other more typical alternatives (U-MOX and variations thereof). Results will be shown indicating that Th-based PWR fuel is a promising option to multi-recycle and burn TRU in a thermal spectrum, while satisfying top-level operational and safety constraints. Various assembly designs have been proposed to assess the TRU burning potential of Th-based fuel in PWRs. In addition to typical homogeneous loading patterns, heterogeneous configurations exploiting the breeding potential of thorium to enable multiple cycles of TRU irradiation and burning have been devised. The homogeneous assembly design, with all pins featuring TRU in Th, has the benefit of a simple loading pattern and the highest rate of TRU transmutation, but it can be used only for a few cycles due to the rapid rise in the TRU content of the recycled fuel, which challenges reactivity control, safety coefficients and fuel handling. Due to its simple loading pattern, such assembly design can be used as the first step of Th implementation, achieving up to 3 times larger TRU transmutation rate than conventional U-MOX, assuming same fraction of MOX assemblies in the core. As the next step in thorium implementation, heterogeneous assemblies featuring a mixed array of Th-U and Th-U-TRU pins, where the U is in-bred from Th, have been proposed. These designs have the potential to enable burning an external supply of TRU through multiple cycles of irradiation, recovery (via reprocessing) and recycling of the residual actinides at the end of each irradiation cycle. This is achieved thanks to a larger breeding of U from Th in the heterogeneous assemblies, which reduces the TRU supply and thus mitigates the increase in the TRU core inventory for the multi-recycled fuel. While on an individual cycle basis the amount of TRU burned in the heterogeneous assembly is reduced with respect to the homogeneous design, TRU burning rates higher than single-pass U-MOX fuel can still be achieved, with the additional benefits of a multi-cycle transmutation campaign recycling all TRU isotopes. Nitride fuel, due its higher density and U breeding potential, together with its better thermal properties, ideally suits the objectives and constraints of the heterogeneous assemblies. However, significant technological advancements must be made before nitride fuels can be employed in an LWR: its water resistance needs to be improved and a viable technology to enrich N in N-15 must be devised. Moreover, for the nitride heterogeneous configurations examined in this study, the enhancement in TRU burning performance is achieved not only by replacing oxide with nitride fuel, but also by increasing the fuel rod size. This latter modification, allowed by the high thermal conductivity of nitride fuel, leads however to a very tight lattice, which may challenge reactor coolant pumps and assembly hold-down mechanisms, the former through an increase in core pressure drop and the latter through an increase in assembly lift-off forces. To alleviate these issues, while still achieving the large fuel-to-moderator ratios resulting from using tight lattices, wire wraps could be used in place of grid spacers. For tight lattices, typical grid spacers are hard to manufacture and their replacement with wire wraps is known to allow for a pressure drop reduction by at least 2 times. The studies, while certainly very preliminary, provide a starting point to devise an optimum strategy for TRU transmutation in Th-based PWR fuel. The viability of the scheme proposed depends on the timely phasing in of the associated technologies, with proper lead time and to solve the many challenges. These challenges are certainly substantial, and make the current once-through U-based scheme pursued in the US by far a more practical (and cheaper) option. However, when compared to other transmutation schemes, the proposed one has arguably similar challenges and unknowns with potentially bigger rewards. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Gripper Test Assembly Status Report (FY2020)

Refueling systems for fast reactors are designed to handle fresh and used core assemblies (fuel, reflector, and shield core assemblies) within the reactor vessel in an opaque coolant environment without visual reference. These refueling machines are designed to work in a sodium (or other fast reactor coolant) and argon vapor space environment and are engineered with the rotatable plug system to allow for the movement of fresh and spent fuel into and out of the reactor core. The refueling machines are a critical component in any reactor and thus need to undergo extensive testing in a prototypic environment to ensure that they will meet all of the system functions and requirements. Argonne has developed an innovative compact refueling system design for the Advanced Fast Reactor-100 that is based upon some mechanisms used in previous reactor designs, such as the U.K.’s Prototype Fast Reactor (PFR) and some mechanisms that have not been used in sodium. This compact refueling machine supports the reduction in size of the AFR-100’s reactor vessel, and if fully developed, would support and inform the development of the in-vessel refueling machines for such commercial reactors as the GEH PRISM reactor plant, the ARC Clean Energy’s ARC-100 reactor, and the Natrium reactor, among others. This refueling system is a vital component of a fast reactor that supports reducing the cost of the reactor and increasing its reliability. During the development of the compact fuel handling machine conceptual design for the AFR-100, a lack of testing data for many mechanical components in sodium under typical in-reactor loads and conditions was discovered. The reduction in lifetime of the various mechanical components in the liquid sodium environment needs to be quantified versus the calculated component lifetimes under normal conditions in the testing while they are subjected to typical loading profiles experienced in the past. The Gripper Test Assembly discussed here includes a full-size gripper device with appropriate mechanical features that will be tested in sodium to provide this testing data. The Gripper Test Assembly is used to test various mechanical fuel handling components submerged in high temperature liquid sodium. These mechanisms are gears, bearings, gripper jaws and head, universal joints and shafts, ball screws, among others. These components will be tested under the typical sodium environmental conditions experienced during refueling operations with appropriate loading conditions that simulate the removal and insertion of core assemblies into a fast reactor grid plate structure. This Gripper Test Assembly is the second in a series of refueling system mechanisms developed for testing in sodium. The first test assembly is the Gear Test Assembly (GTA) which was used to test the performance of gears and bearings operating in sodium. Because of the successful testing conducted with the GTA, it was decided to continue with the development of the Gripper Test Assembly which uses the same gears and bearings tested in GTA. Using the data collected during operation of this gripper test assembly, lifetime reduction factors of the various mechanical components can be calculated for the material combinations selected. These lifetime reduction factors can be used in the design of future mechanical systems which operate in these environments to accurately predict component end of life. In addition, an understanding of the ability of these components and mechanisms to operate under-sodium with the chosen materials will be accomplished. Again, this Gripper Test Assembly is the follow-on test article to the Gear Test Assembly which was testing the ability of gears and bearings operating in a sodium environment. Once fabricated and qualified, it will be tested in the Mechanism Engineering Test Loop (METL) facility located in Building 308 at Argonne National Laboratory.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The MELCOR Analysis of Chinshan Nuclear Power Plant Spent Fuel Pool for Fukushima-like Accident

The safety analysis of the nuclear power plant (NPP) is very important work in the NPP safety. Especially after the Fukushima NPP accident occurred, the importance of NPP safety analysis has been raised and there is more concern for the safety of NPPs in Taiwan. Because the earthquake and tsunami occurred, the cooling system of spent fuel pool failed and the safety issue of spent fuel pool generated in Japan's Fukushima NPP. Chinshan NPP is the first NPP in Taiwan which is BWR/4 plant and OLTP (Original Licensed Thermal Power) for each unit is 1775 MWt. Chinshan NPP finished SPU (stretch power uprate) and the operating power is 103.66% of the OLTP, which is 1840 MWt now. After Fukushima NPP accident occurred, in order to concern the safety of Chinshan NPP spent fuel pool, we performed the safety analysis of spent fuel pool by using MELCOR/SNAP which also assumed the cooling system of spent fuel pool failed. The geometry of Chinshan NPP spent fuel pool is 12.17 m x 7.87 m x 11.61 m and the initial condition is 60 deg. (water temperature) / 1.013 x 10{sup 5} Pa. And, the total power of the fuels is roughly 8.9 MWt initially. MELCOR is developed by Sandia National Laboratories. MELCOR is a fully integrated, engineering-level computer code that models the progression of severe accidents in light water reactor NPPs, including containment response and leakage of activity to the environment. A broad spectrum of severe accident phenomena in both BWRs and PWRs is treated in MELCOR in a unified framework. These include thermalhydraulic response in the reactor coolant system, reactor cavity, containment, and confinement buildings; core heatup, degradation, and relocation; core-concrete attack; hydrogen production, transport, and combustion; fission product release and transport behavior. Current uses of MELCOR include estimation of severe accident source terms and their sensitivities and uncertainties in a variety of applications. According to Dr. Carbajo's paper, MELCOR was used to perform the study of spent fuel pool of Fukushima Daiichi Unit 4. Niina E. Koenoenen used MELCOR to establish the spent fuel pool model of a Nordic BWR. The loss-of-pool-cooling accidents was simulated and analyzed by this model. The above studies indicate that MELCOR is capable of handling the simulation of the spent fuel pool. An increasing number of researchers are using MELCOR code to analyze test facilities and nuclear power plants. Jun Wang et al. used MELCOR to perform core degradation simulation and MELCOR results were compared with CORA experimental data. Through this work, they reviewed the performance of MELCOR COR package in detail. Tae Woon Kim et al. developed MELCOR model of APR1400 (an advanced pressurized water reactor). They used this model to perform the LBLOCA transient simulation and sensitivity study. Longze Li et al. established the MELCOR model of CPR1000 (a Chinese pressurized reactor 1000-MW power plant). The station blackout (SBO) with failure of the steam generator (SG) safety relief valve (SRV) transient was simulated and analyzed by this model. SNAP is a graphic user interface program that processes the inputs and outputs of MELCOR. In addition, there is an animation function in SNAP which can present the animation of analysis results. There are two main steps in this research. The first step was the establishment of Chinshan NPP spent fuel pool MELCOR/SNAP model. And the transient analysis of Chinshan NPP spent fuel pool MELCOR/SNAP model under the cooling system failure condition was performed. The next step was the comparison among the MELCOR, TRACE and CFD data. TRACE and CFD data were from INER reports. In addition, the animation model of Chinshan NPP spent fuel pool was presented by using the animation function of SNAP with MELCOR/SNAP analysis results. This study has developed MELCOR/SNAP model of Chinshan NPP spent fuel pool successfully. By using the above model, the safety analysis of the spent fuel pool was performed under the cooling system of spent fuel pool failed condition. The analysis results of MELCOR, TRACE and CFD were similar in this case. It indicated that there was a respectable accuracy in MELCOR/SNAP model. The analysis results depicted that the uncovered of the fuels occurred at 2.6 day and the metal-water reaction of fuels occurred roughly at 3.7 day after the cooling system failed. The above results indicated that the failure of cladding occurred after 3.7 day. This study's results can help to evaluate the safety issue of Chinshan NPP spent fuel pool. In addition, the cooling system and some mitigation measures of spent fuel pool is not simulated in the above model. However, we will do it in the future. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Polyamorphic Transformations in Fe-Ni-C Liquids: Implications for Chemical Evolution of Terrestrial Planets

During the formation of the Earth's core, the segregation of metallic liquids from silicate mantle should have left behind evident geochemical imprints on both the mantle and the core. Some distinctive geochemical signatures of the mantle-derived rocks likely own their origin to the metal-silicate differentiation of the primitive Earth, setting our planet apart from undifferentiated meteorites as well as terrestrial planets or moons isotopically and compositionally. Understanding the chemical evolution of terrestrial planetary bodies requires knowledge on properties of both liquid iron alloys and silicates equilibrating under physicochemical conditions pertinent to the deep magma ocean. Here we report experimental and computational results on the pressure-induced structural evolution of iron-nickel liquids alloyed with carbon. Our X-ray diffraction experiments up to 7.3 gigapascals (GPa) demonstrate that Fe-Ni (Fe 90 Ni 10 ) liquids alloyed with 3 and 5 wt % carbon undergo a polyamorphic liquid structure transition at approximately 5 GPa. Corroborating the experimental observations, our first-principles molecular dynamic calculations reveal that the structural transitions result from the marked prevalence of three-atom face-sharing polyhedral connections in the liquids at >5 GPa. The structure and polyamorphic transitions of liquid iron-nickel-carbon alloys govern their physical and chemical properties and may thus cast fresh light on the chemical evolution of terrestrial planets and moons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Updates to the Generation of Physics Data Inputs for MAMMOTH Simulations of the Transient Reactor Test Facility - FY2016

The INL is currently evolving the modeling and simulation (M&S) capability that will enable improved core operation as well as design and analysis of TREAT experiments. This M&S capability primarily uses MAMMOTH, a reactor physics application being developed under Multi-physics Object Oriented Simulation Environment (MOOSE) framework. MAMMOTH allows the coupling of a number of other MOOSE-based applications. This second year of work has been devoted to the generation of a deterministic reference solution for the full core, the preparation of anisotropic diffusion coefficients, the testing of the SPH equivalence method, and the improvement of the control rod modeling. In addition, this report includes the progress made in the modeling of the M8 core configuration and experiment vehicle since January of this year.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Simulation of reflooding on two parallel heated channel by TRACE

In case of Loss-Of-Coolant accident (LOCA) in a Boiling Water Reactor (BWR), heat generated in the nuclear fuel is not adequately removed because of the decrease of the coolant mass flow rate in the reactor core. This fact leads to an increase of the fuel temperature that can cause damage to the core and leakage of the radioactive fission products. In order to reflood the core and to discontinue the increase of temperature, an Emergency Core Cooling System (ECCS) delivers water under this kind of conditions. This study is an investigation of how the power distribution between two channels can affect the process of reflooding when the emergency water is injected from the top of the channels. The peak cladding temperature (PCT) on LOCA transient for different axial level is determined as well. A thermal-hydraulic system code TRACE has been used. A TRACE model of the two heated channels has been developed, and three hypothetical cases with different power distributions have been studied. Later, a comparison between a simulated and experimental data has been shown as well.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Suppression of fine-structure splitting and oscillator strength of sodium D-line in a Debye plasma

We investigate theoretically the influence of static plasma screening on relativistic spin-orbit interaction-induced fine-structure splitting of the D-line doublet arising from the transitions 3p{sub 1/2}–3s{sub 1/2} and 3p{sub 3/2}–3s{sub 1/2} of the valence electron of a sodium atom embedded in a model plasma environment. The many-electron atomic problem is formulated first as an effective one-electron problem in which the interaction between the optically active valence electron and the atomic ion core is represented by an accurate parametric model potential including core-polarization correction, and then the plasma effect on the atomic system is simulated by the Debye-screening model for the valence-core interaction. It is observed that the magnitude of spin-orbit energy shift reduces for both the upper component 3p{sub 3/2} and the lower component 3p{sub 1/2} with increasing plasma screening strength, thereby reducing the spin-orbit energy separation between these two components as the screening becomes stronger. As a consequence, the magnitude of fine-structure splitting between the D{sub 1} and D{sub 2} line energies of sodium drops significantly with stronger plasma screening. The optical (absorption) oscillator strength for 3s → 3p transition is seen to reduce with stronger screening and this leads to a screening-induced gradual suppression of the 3p → 3s spontaneous decay rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic flux pumping in 3D nonlinear magnetohydrodynamic simulations

A self-regulating magnetic flux pumping mechanism in tokamaks that maintains the core safety factor at q≈1, thus preventing sawteeth, is analyzed in nonlinear 3D magnetohydrodynamic simulations using the M3D-C1 code. In these simulations, the most important mechanism responsible for the flux pumping is that a saturated (m=1,n=1) quasi-interchange instability generates an effective negative loop voltage in the plasma center via a dynamo effect. It is shown that sawtoothing is prevented in the simulations if β is sufficiently high to provide the necessary drive for the (m=1,n=1) instability that generates the dynamo loop voltage. In conclusion, the necessary amount of dynamo loop voltage is determined by the tendency of the current density profile to centrally peak which, in our simulations, is controlled by the peakedness of the applied heat source profile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cosmic-Ray Feedback Heating of the Intracluster Medium

Active galactic nuclei (AGNs) play a central role in solving the decades-old cooling-flow problem. Although there is consensus that AGNs provide the energy to prevent catastrophically large star formation, one major problem remains: How is the AGN energy thermalized in the intracluster medium (ICM)? We perform a suite of three-dimensional magnetohydrodynamical adaptive mesh refinement simulations of AGN feedback in a cool core cluster including cosmic rays (CRs). CRs are supplied to the ICM via collimated AGN jets and subsequently disperse in the magnetized ICM via streaming, and interact with the ICM via hadronic, Coulomb, and streaming instability heating. We find that CR transport is an essential model ingredient at least within the context of the physical model considered here. When streaming is included, (i) CRs come into contact with the ambient ICM and efficiently heat it, (ii) streaming instability heating dominates over Coulomb and hadronic heating, (iii) the AGN is variable and the atmosphere goes through low-/high-velocity dispersion cycles, and, importantly, (iv) CR pressure support in the cool core is very low and does not demonstrably violate observational constraints. However, when streaming is ignored, CR energy is not efficiently spent on the ICM heating and CR pressure builds up to a significant level, creating tension with the observations. Overall, we demonstrate that CR heating is a viable channel for the AGN energy thermalization in clusters and likely also in ellipticals, and that CRs play an important role in determining AGN intermittency and the dynamical state of cool cores.

79 ASTRONOMY AND ASTROPHYSICS↗

BISON Simulation of MARVEL Fuel Performance

BISON simulation of MARVEL fuel rod performance was conducted for normal operation. Fuel temperatures, plenum pressures and cumulative damage factors were calculated for five core zones. This report presents the results. BISON simulation for transient cases are planned. When the results are available, this report will be updated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Data Collection Methods for Validation of Advanced Multi-Resolution Fast Reactor Simulations

In pool-type Sodium Fast Reactors (SFR) the regions most susceptible to thermal striping are the upper instrumentation structure (UIS) and the intermediate heat exchanger (IHX). This project experimentally and computationally (CFD) investigated the thermal mixing in the region exiting the reactor core to the UIS. The thermal mixing phenomenon was simulated using two vertical jets at different velocities and temperatures as prototypic of two adjacent channels out of the core. Thermal jet mixing of anticipated flows at different temperatures and velocities were investigated. Velocity profiles are measured throughout the flow region using Ultrasonic Doppler Velocimetry (UDV), and temperatures along the geometric centerline between the jets were recorded using a thermocouple array. CFD simulations, using COMSOL, were used to initially understand the flow, then to design the experimental apparatus and finally to compare simulation results and measurements characterizing the flows. The experimental results and CFD simulations show that the flow field is characterized into three regions with respective transitions, namely, convective mixing, (flow direction) transitional, and post-mixing. Both experiments and CFD simulations support this observation. For the anticipated SFR conditions the flow is momentum dominated and thus thermal mixing is limited due to the short flow length associated from the exit of the core to the bottom of the UIS. This means that there will be thermal striping at any surface where poorly mixed streams impinge; rather unless lateral mixing is ‘actively promoted out of the core, thermal striping will prevail. Furthermore we note that CFD can be considered a ‘separate effects (computational) test’ and is recommended as part of any integral analysis. To this effect, poorly mixed streams then have potential impact on the rest of the SFR design and scaling, especially placement of internal components, such as the IHX that may see poorly mixed streams. Finally, due to lack or infrastructural support for carrying out sodium experiments, only water experiments and CFD studies were realized in, an otherwise sodium approved facility.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗