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,459 records · Page 81

Spin-lattice dynamics of surface vs core magnetization in Fe nanoparticles

Magnetization of clusters is often simulated using atomistic spin dynamics for a fixed lattice. Coupled spinlattice dynamics simulations of the magnetization of nanoparticles have up to now neglected the change in the size of the atomic magnetic moments near surfaces. We show that the introduction of variable magnetic moments leads to a better description of experimental data for the magnetization of small Fe nanoparticles. To this end we divide atoms into a surface-near shell and a core with bulk properties. Here, it is demonstrated that both the magnitude of the shell magnetic moment and the exchange interactions need to be modified to obtain a fair representation of the experimental data. This allows for a reasonable description of the average magnetic moment versus cluster size, and also the cluster magnetization versus temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

CEFR simulation using diffusion code system RAST-F

In this study, the RAST-F code system, which is based on the two-step approach that couples a multi-group cross-section generation Monte-Carlo (MC) code and a multi-group nodal diffusion solver, was used for the neutronic simulation of the CEFR start-up experiments. The numerical solution of the RAST-F system was verified against the full core MC solution MCS at all control rods fully inserted and withdrawn states. The RAST-F solution of the selected experimental simulations was compared against the measurement data. A good agreement between RAST-F and MCS solutions was observed with less than 120 pcm discrepancies and 1.2% root-mean-square error in terms of k{sub eff} and power distribution, respectively. Meanwhile, the RAST-F result agreed well with the experimental data within two-sigma of experimental uncertainty. The good agreement of these results indicates that RAST-F can be used to neutronic steady-state simulations for small core-size SFR, which was challenged to deterministic code systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Coupling TRACE with nodal neutronics code Ants using the Exterior Communications Interface and VTT's multi-physics driver Cerberus

In this paper a recently developed coupling between system code TRACE and nodal neutronics code Ants is presented. The TRAC/RELAP Advanced Computational Engine (TRACE) is a reactor system code developed by U.S. Nuclear Regulatory Commission (NRC). The code has been designed to analyse loss-of-coolant accidents (LOCAs), operational transients and other accident scenarios in light water reactors. Ants is the reduced order nodal neutronics code of the Kraken framework. The diffusion solution method in Ants is based on the Analytic Function Expansion Nodal method (AFEN) and Flux Expansion Nodal Method (FENM). Rectangular, hexagonal and triangular geometries are supported. Both steady state and burnup problems including micro-depletion can be simulated. The coupling is implemented using VTT's multi-physics driver Cerberus and the Exterior Communications Interface (ECI) which is supplied with TRACE and enables the coupling of TRACE without source code modifications. The main focus in the paper is on describing the TRACE-Ants coupling in detail in order to share to other code developers what we have learned from the ECI based coupling process. The coupled code system is briefly tested by simulating a simple transient in a SMR core initiated by control rod movement. The results of the test calculation look reasonable and more complicated test cases will be modelled in the near future.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

REMNANTS OF BINARY WHITE DWARF MERGERS

We carry out a comprehensive smooth particle hydrodynamics simulation survey of double-degenerate white dwarf binary mergers of varying mass combinations in order to establish correspondence between initial conditions and remnant configurations. We find that all but one of our simulation remnants share general properties such as a cold, degenerate core surrounded by a hot disk, while our least massive pair of stars forms only a hot disk. We characterize our remnant configurations by the core mass, the rotational velocity of the core, and the half-mass radius of the disk. We also find that some of our simulations with very massive constituent stars exhibit helium detonations on the surface of the primary star before complete disruption of the secondary. However, these helium detonations are insufficiently energetic to ignite carbon, and so do not lead to prompt carbon detonations.

79 ASTRONOMY AND ASTROPHYSICS↗

Core compressor exit stage study. 1: Aerodynamic and mechanical design

The effect of aspect ratio on the performance of core compressor exit stages was demonstrated using two three stage, highly loaded, core compressors. Aspect ratio was identified as having a strong influence on compressors endwall loss. Both compressors simulated the last three stages of an advanced eight stage core compressor and were designed with the same 0.915 hub/tip ratio, 4.30 kg/sec (9.47 1bm/sec) inlet corrected flow, and 167 m/sec (547 ft/sec) corrected mean wheel speed. The first compressor had an aspect ratio of 0.81 and an overall pressure ratio of 1.357 at a design adiabatic efficiency of 88.3% with an average diffusion factor or 0.529. The aspect ratio of the second compressor was 1.22 with an overall pressure ratio of 1.324 at a design adiabatic efficiency of 88.7% with an average diffusion factor of 0.491.

Burdsall, E. A.↗

A Numerical Study on the Effect of Facesheet-Core Disbonds on the Buckling Load of Curved Honeycomb Sandwich Panels

A numerical study on the effect of facesheet-core disbonds on the post-buckling response of curved honeycomb sandwich panels is presented herein. This work was conducted as part of the development of a damage tolerance approach for the next-generation Space Launch System heavy lift vehicle payload fairing. As such, the study utilized full-scale fairing barrel segments as the structure of interest. The panels were composed of carbon fiber reinforced polymer facesheets and aluminum honeycomb core. The panels were analyzed numerically using the finite element method. Facesheet and core nodes in a predetermined circular region were detached to simulate a disbond induced via low-speed impact between the outer mold line facesheet and honeycomb core. Surface-to-surface contact in the disbonded region was invoked to prevent interpenetration of the facesheet and core elements. The diameter of this disbonded region was varied and the effect of the size of the disbond on the post-buckling response was observed. A significant change in the slope of the edge load-deflection response was used to determine the onset of global buckling and corresponding buckling load.

structural failure↗

Core spectroscopy of oxazole

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.

74 ATOMIC AND MOLECULAR PHYSICS↗

Interplay between hydrogen, temperature, and character angle on dissociated dislocation energies in Fe–Ni–Cr austenitic stainless steels

Dislocation energy has an important role in the mechanical performance of structural metals. While dislocation energies cannot be fully obtained from continuum theories due to the contribution of the dislocation core, they have been calculated via atomistic simulations in elemental metals. However, constraints on the local atomic environments have prevented the use of such approaches in systems that incorporate alloying or interstitial solutes. In this work, we develop robust molecular dynamics methods to resolve these issues through a geometric construction of dislocation dipoles and the calculation of time-averaged energies. Furthermore, we apply these methods to calculate dislocation energies (including core energies) in an Fe 70 Ni 11 Cr 19 austenitic steel at a variety of character angles, hydrogen concentrations, temperatures, and dipole spacings. The resulting highly converged energies show an excellent agreement with continuum expressions. Overall, hydrogen concentrations up to 1.0 % do not have a significant effect on the elastic parameters and dislocation energy. The methods and insights derived in this work have the potential to facilitate the calculation of dislocation energies in a wide range of systems, and to guide our understanding of hydrogen embrittlement.

Alloyed systems↗

Deprojection of X-ray data in galaxy clusters: confronting simulations with observations

ABSTRACT Numerical simulations with varying realism indicate an emergent principle−multiphase condensation and large cavity power occur when the ratio of the cooling time to the free-fall time (tcool/tff) falls below a threshold value close to 10. Observations indeed show cool-core signatures when this ratio falls below 20–30, but the prevalence of cores with tcool/tff ratio below 10 is rare as compared to simulations. In X-ray observations, we obtain projected spectra from which we have to infer radial gas density and temperature profiles. Using idealized models of X-ray cavities and multiphase gas in the core and 3D hydro jet-ICM simulations, we quantify the biases introduced by deprojection based on the assumption of spherical symmetry in determining tcool/tff. We show that while the used methods are able to recover the tcool/tff ratio for relaxed clusters, they have an uncertainty of a factor of 2−3 in systems containing large cavities (≳ 20 kpc). We also show that the mass estimates from these methods, in the absence of X-ray spectra close to the virial radius, suffer from a degeneracy between the virial mass (M200) and the concentration parameter (c) in the form of M200c2 ≈ constant. Additionally, the lack of soft-X-ray (≲ 0.5 keV) coverage and poor spatial resolution makes us overestimate min(tcool/tff) by a factor of few in clusters with min(tcool/tff) ≲ 5. This bias can largely explain the lack of cool-core clusters with min(tcool/tff) ≲ 5.

Sarkar, Kartick C. (ORCID:0000000277678472)↗

Verification of bubble tracking method and DNS examinations of single- and two-phase turbulent channel flows

Direct numerical simulation (DNS) has been regarded as a reliable data source for the development and validation of turbulence models along with experiments. The realization of DNS usually involves a very fine mesh that should be able to resolve all relevant turbulence scales down to Kolmogorov scale [1]. As the most computationally expensive approach compared to other CFD techniques, DNS applications used to be limited to flow studies at very low Reynolds numbers. Thanks to the tremendous growth of computing power over the past decades, the simulation capability of DNS has now started overlapping with some of the most challenging engineering problems. One of those examples in nuclear engineering is the turbulent coolant flow inside reactor cores. Coupled with interface tracking methods (ITM), the simulation capability of DNS can be extended to more complicated two-phase flow regimes. Departure from nucleate boiling (DNB) is the limiting critical heat flux phenomena for the majority of accidents that are postulated to occur in pressurized water reactors (PWR) [2]. As one of the major modeling and simulation (M&S) challenges pursued by CASL, the prediction capability is being developed for the onset of DNB utilizing multiphase-CFD (M-CFD) approach. DNS (coupled with ITM) can be employed to provide closure law information for the multiphase flow modeling at CFD scale. In the presented work, research groups at NCSU and UND will focus on applying different ITM to different geometries. Higher void fraction flow analysis at reactor prototypical conditions will be performed, and novel analysis methods will be developed, implemented and verified for the challenging flow conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

PELICAN Status Report for Year 2021

The Pressure Drop Experimental Loop for Investigations of Core Assemblies in Advanced Nuclear Reactors, PELICAN, is a novel experimental facility designed and constructed to recreate the prototypic full scale hydraulic flow conditions for the Versatile Test Reactor (VTR) currently under development by the U.S. Department of Energy. To support VTR design efforts, PELICAN is used to replicate VTR’s primary heat transport system (PHTS) in the core assembly region. The measurement objectives for data generated from this test facility are to mimic the prototypic VTR assembly and to meet the validation needs for code calculations and simulations of the reference design for the VTR core. These objectives focus on the need to validate pressure drop results across the various segments of the fuel assembly, including the axial reflectors, fuel, and transition regions, as they relate directly to the pumping power and safety behavior of the reactor. In the following sections, the major milestones for the experimental work performed on PELICAN over the calendar year (CY) 2021 are presented. In the next section, we discuss the assembly and shakedown testing of the facility. In Section 3, an analysis of the errors in the major measurements is performed. In Section 4 and 5 we present the results from the first set of test articles placed inside the loop. The design of the remaining test articles as part of the full assembly are shown in Section 6 and the conclusion and path forward for the following year, CY22 is given in the final section.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fast ion confinement in the presence of core magnetic islands in Wendelstein 7-X

The effect of magnetic islands in the core region of Wendelstein 7-X (W7-X) on fast ion confinement is explored through simulations with the BEAMS3D code. A magnetic configuration where the n/m = 5/5 island chain is shifted to r/a ~ 0.7 allows the exploration of core island physics in W7-X. The control coil system on W7-X allows the tuning of the island size either increasing the island width or decreasing it. A coupling of the BEAMS3D code to the FIELDLINES code provides a versatile mechanism for incorporating magnetic islands and stochastic regions into the BEAMS3D code. Collisionless simulations suggest that the presence of core islands degrade the confinement of passing particles in the region of the island chain. Full neutral beam simulations of W7-X show a similar behavior with confinement decreasing as the island width is increased. Comparisons between a vacuum magnetic field and low beta HINT2 simulation are made showing similar fast ion behavior. Measurements of lost fast ions in W7-X confirm this trend with the control coil suppressed island configuration showing lower losses than that with no control coils applied. Simulations of fast ion wall loads are performed suggesting no drastic change in loss pattern and a slight reduction in losses with minimized islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-Physics Demonstration Problem with the SHARP Reactor Simulation Toolkit

This report describes to employ SHARP to perform a first-of-a-kind analysis of the core radial expansion phenomenon in an SFR. This effort required significant advances in the framework Multi-Physics Demonstration Problem with the SHARP Reactor Simulation Toolkit used to drive the coupled simulations, manipulate the mesh in response to the deformation of the geometry, and generate the necessary modified mesh files. Furthermore, the model geometry is fairly complex, and consistent mesh generation for the three physics modules required significant effort. Fully-integrated simulations of a 7-assembly mini-core test problem have been performed, and the results are presented here. Physics models of a full-core model of the Advanced Burner Test Reactor have also been developed for each of the three physics modules. Standalone results of each of the three physics modules for the ABTR are presented here, which provides a demonstration of the feasibility of the fully-integrated simulation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Towards Interactive, Reproducible Analytics at Scale on HPC Systems

The growth in scientific data volumes has resulted in a need to scale up processing and analysis pipelines using High Performance Computing (HPC) systems. These workflows need interactive, reproducible analytics at scale. The Jupyter platform provides core capabilities for interactivity but was not designed for HPC systems. In this paper, we outline our efforts that bring together core technologies based on the Jupyter Platform to create interactive, reproducible analytics at scale on HPC systems. Our work is grounded in a real world science use case-applying geophysical simulations and inversions for imaging the subsurface. Our core platform addresses three key areas of the scientific analysis workflow-reproducibility, scalability, and interactivity. We describe our implemention of a system, using Binder, Science Capsule, and Dask software. We demonstrate the use of this software to run our use case and interactively visualize real-Time streams of HDF5 data.

containers↗

University Research in Support of TREAT Modeling and Simulation, FY 2016

Idaho National Laboratory 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 the Multi-physics Object Oriented Simulation Environment (MOOSE) framework. MAMMOTH allows the coupling of a number of other MOOSE-based applications. In support of this research, INL is working with four universities to explore advanced solution methods that will complement or augment capabilities in MAMMOTH. This report consists of a collection of year end summaries of research from the universities performed in support of TREAT modeling and simulation. This research was led by Prof. Sedat Goluoglu at the University of Florida, Profs. Jim Morel and Jean Ragusa at Texas A&M University, Profs. Benoit Forget and Kord Smith at Massachusetts Institute of Technology, Prof. Leslie Kerby of Idaho State University and Prof. Barry Ganapol of University of Arizona. A significant number of students were supported at various levels though the projects and, for some, also as interns at INL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hardware-Based Demonstration of Temperature Control Functions for Reactor Systems

Establishing autonomy in reactor control systems has become essential for the expansion of nuclear technologies. Thermal regulation in particular remains crucial for maintaining stable operation and ensuring the integrity of fuel. To alleviate public skepticism of the safety of nuclear reactors, demonstrating control over this key factor is pivotal. Utilizing electric heat pads to simulate the heat released in a reactor core, thermocouples for temperature monitoring, and an Arduino micro programmable logic controller (PLC) for control, a hardware-based demonstration of a reactor heating system validates the efficacy of reactor control over this key parameter. To improve precision, a proportional-integral-derivative (PID) algorithm was implemented in the heating control loop to ensure meticulous control of reactor functions. In addition, the integration of this physical system with a digital simulator tool such as RELAP5-3D establishes a foundation for a comprehensive testing environment. This allows for a refinement of temperature control under various simulated reactor conditions, bringing another layer of reliability to the operation of the system. By facilitating a physical demonstration of reactor thermal management and control strategies, this project provides a foundation for expanded testing and educational outreach. Ultimately, this system advances the broader goal of demonstrating the safety and viability of autonomous reactor operations, contributing to public trust and future reactor deployment.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A lightweight method for evaluating in situ workflow efficiency

Performance evaluation is crucial to understanding the behavior of scientific workflows. In this study, we target an emerging type of workflow, called in situ workflows. These workflows tightly couple components such as simulation and analysis to improve overall workflow performance. To understand the tradeoffs of various configurable parameters for coupling these heterogeneous tasks, namely simulation stride, and component placement, separately monitoring each component is insufficient to gain insights into the entire workflow behavior. Through an analysis of the state-of-the-art research, we propose a lightweight metric, derived from a defined in situ step, for assessing resource usage efficiency of an in situ workflow execution. By applying this metric to a synthetic workflow, which is parameterized to emulate behaviors of a molecular dynamics simulation, we explore two possible scenarios (Idle Simulation and Idle Analyzer) for the characterization of in situ workflow execution. In addition to preliminary results from a recently published study [11], we further exploit the proposed metric to evaluate a practical in situ workflow with a real molecular dynamics application, i.e., GROMACS. Here, experimental results show that the in transit placement (analytics on dedicated nodes) sustains a higher frequency for performing in situ analysis compared to the helper-core configuration (analytics co-allocated with simulation).

97 MATHEMATICS AND COMPUTING↗

Global particle simulation of lower hybrid wave propagation and mode conversion in tokamaks

Particle-in-cell simulation of lower hybrid (LH) waves in core plasmas is presented with a realistic electron-to-ion mass ratio in toroidal geometry. Due to the fact that LH waves mainly interact with electrons to drive the current, ion dynamic is described by cold fluid equations for simplicity, while electron dynamic is described by drift kinetic equations. This model could be considered as a new method to study LH waves in tokamak plasmas, which has advantages in nonlinear simulations. The mode conversion between slow and fast waves is observed in the simulation when the accessibility condition is not satisfied, which is consistent with the theory. The poloidal spectrum upshift and broadening effects are observed during LH wave propagation in the toroidal geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗