Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron transport equation”

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 55 records · Page 3

Analytical benchmark solution for 1-D neutron transport coupled with thermal conduction and material expansion

In this work we present fully analytical solutions for a class of finite, homogeneous, 1-D slab benchmark problems with nonlinear temperature feedback effects. The proposed class of benchmarks include multiplicative 1-D neutron transport (limited to quasistatic S{sub 2} with μ = ±1) coupled with thermal conduction, convection, Doppler broadening, and expansion effects along the length of the slab. This class of benchmark models, along with the corresponding analytical solutions, are valuable for validating multiphysics analysis frameworks that support coupled neutronics/thermal/structural calculations. Analytical solutions for the benchmarks are obtained by introducing an ansatz that the equilibrium flux and temperature distributions in the slab have the same shape. Specific values for the total microscopic cross section, σ{sub t,0}, and conductive heat transfer coefficient, h, that satisfy the assumed ansatz are then determined. A discussion of the procedure for generating benchmark models and analytical solutions is provided, along with numerical results for an example set of model parameters and thoughts on practical applications of the benchmarks for multiphysics code validation. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions

More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.

42 ENGINEERING↗

Consistent Pl Analysis of Aqueous Uranium-235 Critical Assemblies

The lethargy-dependent equations of the consistent Pl approximation to the Boltzmann transport equation for slowing down neutrons have been used as the basis of an IBM 704 computer program. Some of the effects included are (1) linearly anisotropic center of mass elastic scattering, (2) heavy element inelastic scattering based on the evaporation model of the nucleus, and (3) optional variation of the buckling with lethargy. The microscopic cross-section data developed for this program covered 473 lethargy points from lethargy u = 0 (10 Mev) to u = 19.8 (0.025 ev). The value of the fission neutron age in water calculated here is 26.5 square centimeters; this value is to be compared with the recent experimental value given as 27.86 square centimeters. The Fourier transform of the slowing-down kernel for water to indium resonance energy calculated here compared well with the Fourier transform of the kernel for water as measured by Hill, Roberts, and Fitch. This method of calculation has been applied to uranyl fluoride - water solution critical assemblies. Theoretical results established for both unreflected and fully reflected critical assemblies have been compared with available experimental data. The theoretical buckling curve derived as a function of the hydrogen to uranium-235 atom concentration for an energy-independent extrapolation distance was successful in predicting the critical heights of various unreflected cylindrical assemblies. The critical dimensions of fully water-reflected cylindrical assemblies were reasonably well predicted using the theoretical buckling curve and reflector savings for equivalent spherical assemblies.

Fieno, Daniel↗

(U) Transport Corrections Implemented in SENSMG

Truncating the spherical harmonics expansion of the neutron scattering source in the Boltzmann transport equation leads to truncation error. Methods for approximately correcting for this truncation error are called transport corrections. The PARTISN multigroup discrete ordinates neutron transport code has three transport correction options. These options can be specified in the SENSMG multigroup neutron sensitivity code and passed to PARTISN for the neutron transport.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of deterministic transport methods for low energy neutrons for shielding in space

Transport of low energy neutrons associated with the galactic cosmic ray cascade is analyzed in this dissertation. A benchmark quality analytical algorithm is demonstrated for use with BRYNTRN, a computer program written by the High Energy Physics Division of NASA Langley Research Center, which is used to design and analyze shielding against the radiation created by the cascade. BRYNTRN uses numerical methods to solve the integral transport equations for baryons with the straight-ahead approximation, and numerical and empirical methods to generate the interaction probabilities. The straight-ahead approximation is adequate for charged particles, but not for neutrons. As NASA Langley improves BRYNTRN to include low energy neutrons, a benchmark quality solution is needed for comparison. The neutron transport algorithm demonstrated in this dissertation uses the closed-form Green's function solution to the galactic cosmic ray cascade transport equations to generate a source of neutrons. A basis function expansion for finite heterogeneous and semi-infinite homogeneous slabs with multiple energy groups and isotropic scattering is used to generate neutron fluxes resulting from the cascade. This method, called the FN method, is used to solve the neutral particle linear Boltzmann transport equation. As a demonstration of the algorithm coded in the programs MGSLAB and MGSEMI, neutron and ion fluxes are shown for a beam of fluorine ions at 1000 MeV per nucleon incident on semi-infinite and finite aluminum slabs. Also, to demonstrate that the shielding effectiveness against the radiation from the galactic cosmic ray cascade is not directly proportional to shield thickness, a graph of transmitted total neutron scalar flux versus slab thickness is shown. A simple model based on the nuclear liquid drop assumption is used to generate cross sections for the galactic cosmic ray cascade. The ENDF/B V database is used to generate the total and scattering cross sections for neutrons in aluminum. As an external verification, the results from MGSLAB and MGSEMI were compared to ANISN/PC, a routinely used neutron transport code, showing excellent agreement. In an application to an aluminum shield, the FN method seems to generate reasonable results.

Ganapol, Barry↗

An open-source hybrid unstructured mesh - CAD fusion multiphysics analysis workflow in SALAMANDER

Plasma facing components in fusion devices will endure extreme neutron and heat fluxes. To facilitate their design using simulation tools, the open-source Fusion Module, Fusion ENergy Integrated multiphys-X (FENIX) framework is being developed to model these components with a high-fidelity multi-physics multi-dimensional approach. It can iteratively resolve couplings between all the physics at play, from neutron radiation, to thermomechanics, to near-wall plasma dynamics. This framework is based on the Multiphysics Object Oriented Simulation Environment (MOOSE), which is developed by a collaboration of US National Laboratories since 2008, for advanced nuclear, geomechanics simulations and other applications. FENIX couples numerous simulation tools, including OpenMC, the Tritium Migration Analysis Program v8, the NekRS CFD software, and most MOOSE modules. For the coupling of radiation transport and other physics, FENIX supports a hybrid workflow between Computer Assisted Design (CAD) and unstructured mesh geometries. The CAD can be generated from skinning the unstructured mesh, to enable a coarse geometry for efficient particle transport, but still resolving the local material compositions and temperature gradients. Neutron transport is performed using DAGMC on the CAD, and Cardinal, integrated in FENIX, maps tallied quantities, such as the heat deposition or the tritium generation rates, from a tally volumetric mesh to the other physics’ unstructured mesh. This coupling was exercised on a simplified tokamak geometry, coupling neutron transport with the heat conduction equation, and on a monoblock divertor problem, coupling additionally with tritium migration. Mesh convergence studies highlight the importance of the mapping conservativeness. Coupling with thermo-mechanics is further enabled by the generalization of the approach to moving meshes. The presentation will include these coupled analysis as well as an update on status of the FENIX framework.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational attributes of the integral form of the equation of transfer

Difficulties can arise in radiative and neutron transport calculations when a highly anisotropic scattering phase function is present. In the presence of anisotropy, currently used numerical solutions are based on the integro-differential form of the linearized Boltzmann transport equation. This paper, departs from classical thought and presents an alternative numerical approach based on application of the integral form of the transport equation. Use of the integral formalism facilitates the following steps: a reduction in dimensionality of the system prior to discretization, the use of symbolic manipulation to augment the computational procedure, and the direct determination of key physical quantities which are derivable through the various Legendre moments of the intensity. The approach is developed in the context of radiative heat transfer in a plane-parallel geometry, and results are presented and compared with existing benchmark solutions. Encouraging results are presented to illustrate the potential of the integral formalism for computation. The integral formalism appears to possess several computational attributes which are well-suited to radiative and neutron transport calculations.

Frankel, J. I.↗

Transport @ LANL: An Overview [Slides]

Deterministic neutron transport requires discretization of the Boltzmann transport equation in space, angle, energy, and time. Discretization in energy is typically performed with the multigroup approximation, which tries to preserve continuous energy cross section information by using a weighted average over discrete energy groups. There are two important requirements for accurate multigroup cross sections: appropriate energy group boundaries and an accurate weight function. The former can be selected with appropriate optimizers; the latter either requires flux solutions from comparable simulations or expert knowledge.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Space to Air High-Altitude Region Adjoint Neutron Transport

Neutrons from an atmospheric nuclear explosion can be detected by sensors in orbit. Current tools for characterizing the neutron energy spectrum assume a known source and use forward transport to recreate the detector response. In realistic scenarios the true source is unknown, making this an inefficient, iterative approach. In contrast, the adjoint approach directly solves for the source spectrum, enabling source reconstruction. The time–energy fluence at the satellite and adjoint transport equation allow a Monte Carlo method to characterize the neutron source’s energy spectrum directly in a new model: the Space to High-Altitude Region Adjoint (SAHARA) model. A new adjoint source event estimator was developed in SAHARA to find feasible solutions to the neutron transport problem given the constraints of the adjoint environment. This work explores SAHARA’s development and performance for mono-energetic and continuous neutron energy sources. In general, the identified spectra were shifted towards energies approximately 5% lower than the true source spectra, but SAHARA was able to capture the correct spectral shapes. Continuous energy sources, including real-world sources Fat Man and Little Boy, resulted in identifiable spectra that could have been produced by the same distribution as the true sources as demonstrated by two-dimensional (2D) Kolmogorov–Smirnov tests.

LaMere, Zachary W.↗

Inferring small-sample-reactivity worth from oscillation experiment with auto-corrected local flux

To estimate reactivity changes as a function of time from the signal of detectors count rates, one usually uses point kinetic equations. However, when detectors are located inside the core, it is often necessary to compute a correction factor to account for the flux-shape deformation (i.e. the change of detector efficiency that may have been induced by neutron absorbers movements). Various simulation methods have been developed depending on transients observed. They range from the resolution of static equations up to full-time-dependent neutron transport in 3D. However, those methods require high computational-cost calculations of local reaction rates that rely on nuclear data evaluations. Therefore, a 'measured' reactivity is never determined through pure measurements. Taking advantage of the known oscillation movement of a small sample at the center of the MINERVE reactor, the innovative algorithm presented in this paper derives the correction factor needed to get no time drift on reactivity steps. Derived correction factors are in good agreement with those obtained independently with TRIPOLI-4 through sensitivity analysis. Low sensitivity to the sample composition and moderate sensitivity to the detectors position reinforce the reliability of those results. This application paves the way to almost pure nuclear data reactivity measurements. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Transient Multi-Level method for Monte Carlo reactor statics calculations

The Transient Multi-Level (TML) method is applied to a time-dependent Monte Carlo transport solver to offload some of the computational burden of the expensive Monte Carlo solve to lower-order Coarse Mesh Finite Difference (CMFD) and Exact Point Kinetics Equations (EPKE) solvers via factorization of the neutron flux at the transport and CMFD levels using the Predictor Corrector Quasi-Static Method (PCQM). The Monte Carlo transient is solved by a modified fission source iteration scheme that introduces a single transient source bank. The method is implemented in the production-level Monte Carlo code, Shift, and verified with prescribed reactivity ramps from the two-dimensional version of the C5G7-TD reactor benchmark. The results show that, as compared to other quasi-static methods, the TML reduces the stochastic noise inherent to the transient Monte Carlo solver by factors of ~2 to 6 for various norm comparisons of the reactor power amplitude. Finally, the TML additionally reduces the number of Monte Carlo evaluations needed to simulate the transient, leading to roughly an order of magnitude improvement in CPU time relative to the standard PCQM for the problems tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An Improved Neutron Transport Algorithm for Space Radiation

A low-energy neutron transport algorithm for use in space radiation protection is developed. The algorithm is based upon a multigroup analysis of the straight-ahead Boltzmann equation by using a mean value theorem for integrals. This analysis is accomplished by solving a realistic but simplified neutron transport test problem. The test problem is analyzed by using numerical and analytical procedures to obtain an accurate solution within specified error bounds. Results from the test problem are then used for determining mean values associated with rescattering terms that are associated with a multigroup solution of the straight-ahead Boltzmann equation. The algorithm is then coupled to the Langley HZETRN code through the evaporation source term. Evaluation of the neutron fluence generated by the solar particle event of February 23, 1956, for a water and an aluminum-water shield-target configuration is then compared with LAHET and MCNPX Monte Carlo code calculations for the same shield-target configuration. The algorithm developed showed a great improvement in results over the unmodified HZETRN solution. In addition, a two-directional solution of the evaporation source showed even further improvement of the fluence near the front of the water target where diffusion from the front surface is important.

Heinbockel, John H.↗

The cosmic ray electron spectrum and its modulation from 1968 through 1972

Over the past five years we have measured the energy spectrum of primary cosmic ray electrons with both a balloon-borne and a satellite absorption spectrometer. All of the balloon flights used identical equipment that was launched each summer from Fort Churchill, Manitoba, Canada. The satellite, OGO-5, has been in an eccentric orbit since March 1968. Together these instruments provide the electron spectrum over a range of energy from 20 MeV to 20 GeV. This wide range and the substantial span of time covered by the measurements permit a detailed study of the solar modulation of electrons. These results are compared with the modulation of the nuclear components as observed by a neutron monitor and interpreted using the cosmic ray transport equation.

Fulks, G.↗

Investigating the AGN-201M research reactor's unique dominance ratio

The AGN-201M reactor is a small, thermal spectrum reactor located at the University of New Mexico (UNM). It is moderated by polyethylene, reflected by graphite, with fuel comprised of uranium microspheres embedded in polyethylene plates that are separated by an aluminum baffle. The dominance ratio of a multiplying system is the ratio of the first higher-mode eigenvalue to the fundamental, k{sub 1}/k{sub 0}. It is a key parameter in characterizing the rate of convergence, analyzing system stability, and determining how tightly coupled the system's neutron fission is. The purpose of this paper is to investigate the system's unusually low dominance ratio of 0.632. This investigation varied the reactor's physical parameters, such as the fuel geometry, fuel density, and reflector thickness, and examined their impact on the reactor's dominance ratio. Additionally, this study examined neutronically similar systems in hopes of identifying common causes for systems with low dominance ratios. The reason why the dominance ratio of the reactor is small compared to larger thermal reactors was determined to be because of its significantly smaller geometry and the composition of the fuel plates. The reflector's effect on the dominance ratio is small in comparison to the other factors, but was found to have a non-zero effect. Furthermore, the AGN-201M was found to have a significantly lower dominance ratio than systems in which it shares a very high (c{sub k} >95%) degree of neutronic similarity.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development for Integrated System-Level Analysis Capabilities in SAM for Molten Salt Reactors

In recent years, there has been renewed interest in Molten Salt Reactors (MSRs) for their potential advantages compared to reactors that rely on solid fuel. In response to such interest, many methods and codes have been developed to capture the unique features of MSRs. Among them, the System Analysis Module (SAM) is a modern system analysis tool that provides fast-running, modest-fidelity, whole-plant transient analysis capabilities, essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. For liquid-fuel MSRs, the complex physics and chemistry involved in MSR operation—such as reactor kinetics, fluid flow, heat transfer, and salt composition dynamics—pose significant challenges for system-level modeling. Specific modeling capabilities are needed for system-level transient simulation. This paper presents recent advancements in SAM capability enhancements for system-level modeling of MSRs, focusing on improved simulation fidelity, computational efficiency, and multi-physics integration. Key enhancements include the development of species transport, Delayed Neutron Precursor (DNP) drift, modified Point Kinetics Equations (PKE), decay heat modeling, key fission product behavior, salt corrosion, and thermal-hydraulic coupling, as well as code robustness and performance enhancements for MSR applications. The code enhancement allows for better predictive accuracy in safety analysis, transient behavior, and operational optimization, thus supporting the design and licensing of next-generation MSRs. Results from case studies are presented to demonstrate the benefits of these enhancements in accurately capturing key reactor transient behaviors.

Hu, Rui (ORCID:0000000237712920)↗

Thermal Effects in Binary Neutron Star Mergers

Abstract We study the impact of finite-temperature effects in numerical-relativity simulations of binary neutron star mergers with microphysical equations of state and neutrino transport in which we vary the effective nucleon masses in a controlled way. We find that, as the specific heat is increased, the merger remnants become colder and more compact due to the reduced thermal pressure support. Using a full Bayesian analysis, we demonstrate that this effect will be measurable in the postmerger gravitational wave signal with next-generation observatories at signal-to-noise ratios of 15.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MPACT Software Management Plan (V.4.3)

The MPACT code solves a discretized form of the Boltzmann transport equation on a wide variety of geometries and is distributed with a multigroup neutron cross section library. MPACT provides an advanced geometrically resolved neutral-particle transport capability to solve the flux distribution throughout the entire problem geometry, and it can model the isotopic depletion, decay, and activation of materials. The flux solution in MPACT is provided using a 2D/1D synthesis method within the framework of the 3D coarse mesh finite difference (CMFD) method for which axial and radial correction factors are obtained from 2D method of characteristics (MOC) and 1D nodal expansion method (NEM), PN, or SN. Other key characteristics of the MPACT code include the subgroup method and the embedded self-shielding method (ESSM) for resonance treatment, depletion capability based on the ORIGEN exponential matrix method, and a simplified thermal-hydraulics method for temperature/fluid feedback. The sole purpose of the simplified feedback model is to provide a mechanism for testing during code development and to provide a limited capability for educational applications. Work performed at the code level supports the VERA-QA-001, quality assurance program plan (QAPP) and VERA-QA-002, VERA Software Quality Assurance Plan.

97 MATHEMATICS AND COMPUTING↗

Linear source method of characteristics in HELIOS2

The Studsvik Scandpower code HELIOS2 is a fast lattice transport code that has gained wide use for reactor fuel cycle analysis in commercial, regulatory, and research environments. It contains two transport solvers - Collision Probabilities (CP) and Method of Characteristics (MoC) - with a general geometry ability and generalized depletion for a variety of fuel types, experimental facilities, and isotope production applications. The HELIOS2 lattice code has been extended with the development of a linear source MoC option. The linear source equations are derived within the scope of the HELIOS2 general geometry and non-cyclic ray tracing. The implementation was tested for fine and coarse spatial mesh models of the 2D C5G7 MOX computational benchmark, which provides a pure and extreme test of a transport solver. The linear source yields at least comparable accuracy with reduced computational resource requirements relative to the flat-source approximation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗