Three-dimensional multigroup diffusion code ANDEX based on nodal method for cartesian geometry
An analytic polynomial nodal method using partial currents has been derived for the solution of multigroup neutron diffusion equations in three-dimensional (3-D) cartesian geometry. This method is characterized by expressing the source and leakage terms in an auxiliary 1-D diffusion equation by quadratic polynomials and solving it analytically. Based on this method, we have developed a 3-D multigroup diffusion code ANDEX, and applied to 2-D LWR and 3-D FBR models. The results of keff, power distributions and computing time have been compared with those of finite difference method calculations. (author)
Development and validation of Apros multigroup nodal diffusion model
Rintala, Antti
2015-01-01
The development of a steady state and transient multigroup nodal diffusion model for process simulation software Apros was continued and the models were validated. The initial implementation of the model was performed in 2009 and it has not been under continuous development afterwards. Some errors in the steady state model were corrected. The transient model was found to be incorrect. The solution method of the transient model was derived, and the program code not common with the steady s...
FINELM: a multigroup finite element diffusion code
FINELM is a FORTRAN IV program to solve the Neutron Diffusion Equation in X-Y, R-Z, R-theta, X-Y-Z and R-theta-Z geometries using the method of Finite Elements. Lagrangian elements of linear or higher degree to approximate the spacial flux distribution have been provided. The method of dissections, coarse mesh rebalancing and Chebyshev acceleration techniques are available. Simple user defined input is achieved through extensive input subroutines. The input preparation is described followed by a program structure description. Sample test cases are provided. (Auth.)
1-D EQUILIBRIUM DISCRETE DIFFUSION MONTE CARLO
T. EVANS; ET AL
2000-08-01
We present a new hybrid Monte Carlo method for 1-D equilibrium diffusion problems in which the radiation field coexists with matter in local thermodynamic equilibrium. This method, the Equilibrium Discrete Diffusion Monte Carlo (EqDDMC) method, combines Monte Carlo particles with spatially discrete diffusion solutions. We verify the EqDDMC method with computational results from three slab problems. The EqDDMC method represents an incremental step toward applying this hybrid methodology to non-equilibrium diffusion, where it could be simultaneously coupled to Monte Carlo transport.
Converged accelerated finite difference scheme for the multigroup neutron diffusion equation
Computer codes involving neutron transport theory for nuclear engineering applications always require verification to assess improvement. Generally, analytical and semi-analytical benchmarks are desirable, since they are capable of high precision solutions to provide accurate standards of comparison. However, these benchmarks often involve relatively simple problems, usually assuming a certain degree of abstract modeling. In the present work, we show how semi-analytical equivalent benchmarks can be numerically generated using convergence acceleration. Specifically, we investigate the error behavior of a 1D spatial finite difference scheme for the multigroup (MG) steady-state neutron diffusion equation in plane geometry. Since solutions depending on subsequent discretization can be envisioned as terms of an infinite sequence converging to the true solution, extrapolation methods can accelerate an iterative process to obtain the limit before numerical instability sets in. The obtained results have been compared to the analytical solution to the 1D multigroup diffusion equation when available, using FORTRAN as the computational language. Finally, a slowing down problem has been solved using a cascading source update, showing how a finite difference scheme performs for ultra-fine groups (104 groups) in a reasonable computational time using convergence acceleration. (authors)
Multigroup finite element-boundary element method for neutron diffusion
Full text: The finite element method (FEM) is an efficient method used for the solution of partial differential equations (PDE's) of engineering physics due to its symmetric, sparse and positive-definite coefficient matrix. FEM has been successfully applied for the solution of multigroup neutron transport and diffusion equations since 1970's. The boundary element method (BEM), on the other hand, is a newer method and is unique among the numerical methods used for the solution of PDE's with its property of confining the unknowns only to the boundaries of homogeneous regions, thus, greatly reducing matrix dimensions. The first application of BEM to the neutron diffusion equation (NDE) dates back to 1985 and many researchers are currently working in this area. Although BEM is known to have the desirable property of being an internal-mesh free method, this advantage is lost in some of its application to the NDE due to the existence of fission source volume integrals in fissionable regions unless domain-decomposition methods are used. To exploit the favorable properties of both FEM and BEM, a hybrid FE/BE method has been recently proposed for reflected systems treated by one or two-group diffusion theories in a recent paper co-authored by the first author. In this work, the hybrid FE/BE method for reflected systems is generalized to multigroup diffusion theory. The core is treated by FEM to preserve the high accuracy of FEM in such neutron-producing regions. Using a boundary integral equation formerly proposed by the second author, BEM, is utilized for the discretization of the reflector, thus, eliminating the internal mesh completely for this nonfissionable region. The multigroup FE/BE method has been implemented in our recently developed FORTRAN program. The program is validated by comparison of the calculated effective multiplication factor and the group fluxes with their analytical counterparts for a two-group reflected system. Comparison of these results and
A code called COMESH based on corner mesh finite difference scheme has been developed to solve multigroup diffusion theory equations. One can solve 1-D, 2-D or 3-D problems in Cartesian geometry and 1-D (r) or 2-D (r-z) problem in cylindrical geometry. On external boundary one can use either homogeneous Dirichlet (θ-specified) or Neumann (∇θ specified) type boundary conditions or a linear combination of the two. Internal boundaries for control absorber simulations are also tackled by COMESH. Many an acceleration schemes like successive line over-relaxation, two parameter Chebyschev acceleration for fission source, generalised coarse mesh rebalancing etc., render the code COMESH a very fast one for estimating eigenvalue and flux/power profiles in any type of reactor core configuration. 6 refs. (author)
The subject is divided in two parts: In the first part a nodal method has been worked out to solve the steady state multigroup diffusion equation. This method belongs to the same set of nodal methods currently used to calculate the exact fission powers and neutron fluxes in a very short computing time. It has been tested on a two dimensional idealized reactors. The effective multiplication factor and the fission powers for each fuel element have been calculated. The second part consists in studying and mastering the multigroup diffusion code DAHRA - a reduced version of DIANE - a two dimensional code using finite difference method
FINELM: a multigroup finite element diffusion code. Part I
The author presents a two dimensional code for multigroup diffusion using the finite element method. It was realized that the extensive connectivity which contributes significantly to the accuracy, results in a matrix which, although symmetric and positive definite, is wide band and possesses an irregular profile. Hence, it was decided to introduce sparsity techniques into the code. The introduction of the R-Z geometry lead to a great deal of changes in the code since the rotational invariance of the removal matrices in X-Y geometry did not carry over in R-Z geometry. Rectangular elements were introduced to remedy the inability of the triangles to model essentially one dimensional problems such as slab geometry. The matter is discussed briefly in the text in the section on benchmark problems. This report is restricted to the general theory of the triangular elements and to the sparsity techniques viz. incomplete disections. The latter makes the size of the problem that can be handled independent of core memory and dependent only on disc storage capacity which is virtually unlimited. (Auth.)
1 - Description of program or function: PARTISN (Parallel, Time-Dependent SN) is the evolutionary successor to CCC-0547/DANTSYS. User input and cross section formats are very similar to that of DANTSYS. The linear Boltzmann transport equation is solved for neutral particles using the deterministic (SN) method. Both the static (fixed source or eigenvalue) and time-dependent forms of the transport equation are solved in forward or adjoint mode. Vacuum, reflective, periodic, white, or inhomogeneous boundary conditions are solved. General anisotropic scattering and inhomogeneous sources are permitted. PARTISN solves the transport equation on orthogonal (single level or block-structured AMR) grids in 1-D (slab, two-angle slab, cylindrical, or spherical), 2-D (X-Y, R-Z, or R-T) and 3-D (X-Y-Z or R-Z-T) geometries. 2 - Methods:PARTISN numerically solves the multigroup form of the neutral-particle Boltzmann transport equation. The discrete-ordinates form of approximation is used for treating the angular variation of the particle distribution. For curvilinear geometries, diamond differencing is used for angular discretization. The spatial discretizations may be either low-order (diamond difference or Adaptive Weighted Diamond Difference (AWDD)) or higher-order (linear discontinuous or exponential discontinuous). Negative fluxes are eliminated by a local set-to-zero-and-correct algorithm for the diamond case (DD/STZ). Time differencing is Crank-Nicholson (diamond), also with a set-to-zero fix-up scheme. Both inner and outer iterations can be accelerated using the diffusion synthetic acceleration method, or transport synthetic acceleration can be used to accelerate the inner iterations. The diffusion solver uses either the conjugate gradient or multigrid method. Chebyshev acceleration of the fission source is used. The angular source terms may be treated either via standard PN expansions or Galerkin scattering. An option is provided for strictly positive scattering sources
Highlights: → Coupled neutron and gamma transport is considered in the multigroup diffusion approximation. → The model accommodates fission, up- and down-scattering and common neutron-gamma interactions. → The exact solution to the diffusion equation in a heterogeneous media of any number of regions is found. → The solution is shown to parallel the one-group case in a homogeneous medium. → The discussion concludes with a heterogeneous, 2 fuel-plate 93.2% enriched reactor fuel benchmark demonstration. - Abstract: The angular flux for the 'rod model' describing coupled neutron/gamma (n, γ) diffusion has a particularly straightforward analytical representation when viewed from the perspective of a one-group homogeneous medium. Cast in the form of matrix functions of a diagonalizable matrix, the solution to the multigroup equations in heterogeneous media is greatly simplified. We shall show exactly how the one-group homogeneous medium solution leads to the multigroup solution.
Multigroup neutron transport equation in the diffusion and P1 approximation
Investigations of the properties of the multigroup transport operator, width and without delayed neutrons in the diffusion and P1 approximation, is performed using Keldis's theory of operator families as well as a technique . recently used for investigations into the properties of the general linearized Boltzmann operator. It is shown that in the case without delayed neutrons, multigroup transport operator in the diffusion and P1 approximation possesses a complete set of generalized eigenvectors. A formal solution to the initial value problem is also given. (author)
APPLE, Plot of 1-D Multigroup Neutron Flux and Gamma Flux and Reaction Rates from ANISN
A - Description of problem or function: The APPLE-2 code has the following functions: (1) It plots multi-group energy spectra of neutron and/or gamma ray fluxes calculated by ANISN, DOT-3.5, and MORSE. (2) It gives an overview plot of multi-group neutron fluxes calculated by ANISN and DOT-3.5. The scalar neutron flux phi(r,E) is plotted with the spatial parameter r linear along the Y-axis, logE along the X-axis and log phi(r,E) in the Z direction. (3) It calculates the spatial distribution and region volume integrated values of reaction rates using the scalar flux calculated with ANISN and DOT-3.5. (4) Reaction rate distribution along the R or Z direction may be plotted. (5) An overview plot of reaction rates or scalar fluxes summed over specified groups may be plotted. R(ri,zi) or phi(ri,zi) is plotted with spatial parameters r and z along the X- and Y-axes in an orthogonal coordinate system. (6) Angular flux calculated by ANISN is rearranged and a shell source at any specified spatial mesh point may be punched out in FIDO format. The shell source obtained may be employed in solving deep penetration problems with ANISN, when the entire reactor system is divided into two or more parts and the neutron fluxes in two adjoining parts are connected by using the shell source. B - Method of solution: (a) The input data specification is made as simple as possible by making use of the input data required in the radiation transport code. For example, geometry related data in ANISN and DOT are transmitted to APPLE-2 along with scalar flux data so as to reduce duplicity and errors in reproducing these data. (b) Most the input data follow the free form FIDO format developed at Oak Ridge National Laboratory and used in the ANISN code. Furthermore, the mixture specifying method used in ANISN is also employed by APPLE-2. (c) Libraries for some standard response functions required in fusion reactor design have been prepared and are made available to users of the 42-group neutron
SIXTUS-2. A two dimensional multigroup diffusion theory code in hexagonal geometry. Pt. 1
A new algorithm for solving the 2-dimensional multigroup diffusion equations in hexagonal geometry is described. It is based on three novel ideas: analytic intranodal solutions, use of the group irreducible representations and an explicit scheme for solving the response matrix equations. The resulting computer code SIXTUS-2 has been found to be very accurate and effective. (Auth.)
Second order time evolution of the multigroup diffusion and P1 equations for radiation transport
Highlights: → An existing multigroup transport algorithm is extended to be second-order in time. → A new algorithm is presented that does not require a grey acceleration solution. → The two algorithms are tested with 2D, multi-material problems. → The two algorithms have comparable computational requirements. - Abstract: An existing solution method for solving the multigroup radiation equations, linear multifrequency-grey acceleration, is here extended to be second order in time. This method works for simple diffusion and for flux-limited diffusion, with or without material conduction. A new method is developed that does not require the solution of an averaged grey transport equation. It is effective solving both the diffusion and P1 forms of the transport equation. Two dimensional, multi-material test problems are used to compare the solution methods.
Tredit A 3-D multigroup diffusion theory simulator for hexagonal fuel assembly cores
A multigroup 3-D reactor core simulator based on neutron diffusion theory, called TREDIT has been developed for Light Water Reactors (LWRs). It considers triangle shaped meshes in X-Y plane and variable mesh spacing in Z-direction. Thus it is especially suited for designing and analysing LWR cores with hexagonal fuel assemblied like the Russian WWER reactors. When fuel assembly cross-sections in multigroup form are input as fitted constants, the computer code TREDIT can build up core burnup distribution with power distribution computed for initial reactor conditions. The results of this code have been compared with another diffusion theory based code and found satisfactory. Xenon feedback effects on core power distribution are demonstrated. (author)
Cassandre : a two-dimensional multigroup diffusion code for reactor transient analysis
CASSANDRE is a two-dimensional (x-y or r-z) finite element neutronics code with thermohydraulics feedback for reactor dynamics prior to the disassembly phase. It uses the multigroup neutron diffusion theory. Its main characteristics are the use of a generalized quasistatic model, the use of a flexible multigroup point-kinetics algorithm allowing for spectral matching and the use of a finite element description. The code was conceived in order to be coupled with any thermohydraulics module, although thermohydraulics feedback is only considered in r-z geometry. In steady state criticality search is possible either by control rod insertion or by homogeneous poisoning of the coolant. This report describes the main characterstics of the code structure and provides all the information needed to use the code. (Author)
FINELM: a multigroup finite element diffusion code. III
The authors describe the formalism of the finite element for r-theta geometry. The treatment of a four sided element with two circular arcs and two radial sides is presented. At the centre of the disc this element is collapsed into a three sided element by reducing the inner arc length to zero. The singularity of the leakage term of such an element at r=0 is described. The modelling of control and safety rods by introducing boundary conditions interior to the reactor is presented. This has the effect of converting the diffusion region from a simply connected to a multiply connected region. Dirichlet conditions with prescribed flux values zero are imposed on nodes in the interior of the non diffusion zones in order to decouple these nodes from the rest. The method is tested by calculating the effectiveness of control and safety rods for a 500 MW Th high temperature reactor. The results are compared with those obtained with DIFGEN using x-y geometry. The DIFGEN results and the problem data are by courtsey of Hochtemperatur Reaktor Bau AG, Mannheim. (Auth.)
GIS-BASED 1-D DIFFUSIVE WAVE OVERLAND FLOW MODEL
KALYANAPU, ALFRED [Los Alamos National Laboratory; MCPHERSON, TIMOTHY N. [Los Alamos National Laboratory; BURIAN, STEVEN J. [NON LANL
2007-01-17
This paper presents a GIS-based 1-d distributed overland flow model and summarizes an application to simulate a flood event. The model estimates infiltration using the Green-Ampt approach and routes excess rainfall using the 1-d diffusive wave approximation. The model was designed to use readily available topographic, soils, and land use/land cover data and rainfall predictions from a meteorological model. An assessment of model performance was performed for a small catchment and a large watershed, both in urban environments. Simulated runoff hydrographs were compared to observations for a selected set of validation events. Results confirmed the model provides reasonable predictions in a short period of time.
Development of a three-dimensional multigroup nodal diffusion code for the LMR
STEP is a three-dimensional multigroup nodal diffusion code for the neutronics analysis of the LMR core and accepts microscopic cross section data. Material cross sections are obtained by summing the product of atom densities and microscopic cross sections over all isotopes comprising the material. STEP contains a thermal-hydraulics module which enables feedback effects from both fuel temperature and coolant temperature changes. Numerical results of the STEP code over the KALIMER core (392 MWt) agree well with those of DIF-3D. And it has been observed that the thermal-hydraulics module is working properly
An effective method of solving the multigroup diffusion problem in hexagonal geometry. Part I
An effective method of solving two-dimensional multigroup diffusion equations in hexagonal geometry is described. The method is based on the following two ideas: nodal approach, and expansion of one-dimensional neutron fluxes inside the node into polynomials up to the third order. The resulting relations for the interface-averaged partial currents, node-averaged fluxes and flux moments are used in computer code NEHEX. The code was found to be an accurate and effective computational tool. Its description and validation against reference benchmark problems will be published as Part II of this report. (author) 1 fig., 1 tab., 9 refs
Simulate-HEX - The multi-group diffusion equation in hexagonal-z geometry
The multigroup diffusion equation is solved for the hexagonal-z geometry by dividing each hexagon into 6 triangles. In each triangle, the Fourier solution of the wave equation is approximated by 8 plane waves to describe the intra-nodal flux accurately. In the end an efficient Finite Difference like equation is obtained. The coefficients of this equation depend on the flux solution itself and they are updated once per power/void iteration. A numerical example demonstrates the high accuracy of the method. (authors)
FEM-BABEL, 3-D Multigroup Neutron Diffusion by Galerkin Method
1 - Nature of the physical problem solved: This program computes the three-dimensional multigroup neutron diffusion equation using the finite element method. 2 - Method of solution: The equation is solved using a solution algorithm based on a Galerkin-type scheme. Prism and box-shaped finite elements are used. The resulting equation system is solved using the successive over-relaxation method and the inner iterations are accelerated by a coarse mesh re-balancing technique. 3 - Restrictions on the complexity of the problem: Any down-scattering of neutrons is allowed but up-scattering and region-dependent fission spectra are not permitted
The good features of the Analytic Function Expansion Nodal (AFEN) method are utilized to develop a practical scheme for the multigroup diffusion problems, in combination with the polynomial expansion nodal (PEN) method. The thermal group fluxes exhibiting strong gradients are solved by the AFEN method, while the fast group fluxes that are smoother than the thermal group fluxes are solved by the PEN method. The scheme is developed for cores of rectangular and hexagonal geometries. In particular, to model the fast group fluxes in the hexagonal geometry by the PEN method, a polynomial function set which shows good performance in accuracy and numerical stability is derived, in premiere. (author)
Nonlinear diffusion acceleration (NDA) can improve the performance of a neutron transport solver significantly especially for the multigroup eigenvalue problems. The high-order transport equation and the transport-corrected low-order diffusion equation form a nonlinear system in NDA, which can be solved via a Picard iteration. The consistency of the correction of the low-order equation is important to ensure the stabilization and effectiveness of the iteration. It also makes the low-order equation preserve the scalar flux of the high-order equation. In this paper, the consistent correction for a particular discretization scheme, self-adjoint angular flux (SAAF) formulation with discrete ordinates method (SN) and continuous finite element method (CFEM) is proposed for the multigroup neutron transport equation. Equations with the anisotropic scatterings and a void treatment are included. The Picard iteration with this scheme has been implemented and tested with RattleSNake, a MOOSE-based application at INL. Convergence results are presented. (authors)
Most of the neutron diffusion codes use numerical methods giving accurate results in structured meshes. However, the application of these methods in unstructured meshes to deal with complex geometries is not straightforward and it may cause problems of stability and convergence of the solution. By contrast, the Finite Volume Method (FVM) is easily applied to unstructured meshes and is typically used in the transport equations due to the conservation of the transported quantity within the volume. In this paper, the FVM algorithm implemented in the ARB Partial Differential Equations Solver has been used to discretize the multigroup neutron diffusion equation to obtain the matrices of the generalized eigenvalue problem, which has been solved by means of the SLEPc library. Nevertheless, these matrices could be large for fine meshes and the eigenvalue problem resolution could require a high calculation time. Therefore, a transformation of the generalized eigenvalue problem into a standard one is performed in order to reduce the calculation time. (author)
Solution of the 1D kinetic diffusion equations using a reduced nodal cubic scheme
In this work it is described a novel method to solve the multi-group time-dependent diffusion equations based on a nodal cubic space interpolation in addition to the application of quadrature rules simplifying the stiffness and mass matrices arising in a finite element procedure. Numerical results for a well known benchmark problem are also provided. (authors)
In this work we developed a software to model and generate results in tables and graphs of one-dimensional neutron transport problems in multi-group formulation of energy. The numerical method we use to solve the problem of neutron diffusion is analytic, thus eliminating the truncation errors that appear in classical numerical methods, e.g., the method of finite differences. This numerical analytical method increases the computational efficiency, since they are not refined spatial discretization necessary because for any spatial discretization grids used, the numerical result generated for the same point of the domain remains unchanged unless the rounding errors of computational finite arithmetic. We chose to develop a computational application in MatLab platform for numerical computation and program interface is simple and easy with knobs. We consider important to model this neutron transport problem with a fixed source in the context of shielding calculations of radiation that protects the biosphere, and could be sensitive to ionizing radiation
Lozano Montero, Juan Andrés; García Herranz, Nuria; Ahnert Iglesias, Carolina; Aragonés Beltrán, José María
2008-01-01
In this work we address the development and implementation of the analytic coarse-mesh finite-difference (ACMFD) method in a nodal neutron diffusion solver called ANDES. The first version of the solver is implemented in any number of neutron energy groups, and in 3D Cartesian geometries; thus it mainly addresses PWR and BWR core simulations. The details about the generalization to multigroups and 3D, as well as the implementation of the method are given. The transverse integration procedure i...
An extended version of Hassitt's one-dimension Multi-group diffusion programme has been prepared which allows for a maximum of twenty-two energy groups rather than sixteen. It also permits the use of drum storage by programme and data up to the full machine capacity of 1024 sectors, rather than 478 sectors. Some minor corrections have been made, A binary tape of an 830-sector version has been prepared (Winfrith P.5272). (author)
The fine mesh diffusion formulation is extended to deal with multigroup 3-D problems in rectangular geometries. The formulation includes interface discontinuity factors per cell type, pre-calculated from transport solutions. The iterative scheme, aiming to an efficient parallel implementation in memory distributed multi-processors, is based on domain decomposition in the 4 possible sets of 4 neighbor quarters of assemblies. The alternate dissections achieve convergence to the exact boundary conditions, while attenuating high frequency noise. Whole core convergence is accelerated in the long wavelength effects by a consistent high-order analytical nodal solution performed by the ANDES solver. A neutronics - thermal-hydraulics iterative scheme is also developed to compute best estimate results, by coupling at the detailed cell-subchannel scale the COBAYA3 code with several TH subchannel codes. The numerical performance and convergence rates are verified by computing pin-cell scale solutions for the OECD/NEA/USNRC PWR MOX/UO2 Core Transient Benchmark in 8 energy groups and heterogeneous assemblies. The cell-subchannel scale neutronics and thermal-hydraulics coupling, allows the verification of the effects of the detailed TH feedbacks on cross-sections and, thus, on fuel pin powers, calculated here for a 3D color-set of two different fuel types of the previous benchmark, using COBAYA3 and COBRA-3C. (authors)
Three-dimensional h-adaptivity for the multigroup neutron diffusion equations
Wang, Yaqi
2009-04-01
Adaptive mesh refinement (AMR) has been shown to allow solving partial differential equations to significantly higher accuracy at reduced numerical cost. This paper presents a state-of-the-art AMR algorithm applied to the multigroup neutron diffusion equation for reactor applications. In order to follow the physics closely, energy group-dependent meshes are employed. We present a novel algorithm for assembling the terms coupling shape functions from different meshes and show how it can be made efficient by deriving all meshes from a common coarse mesh by hierarchic refinement. Our methods are formulated using conforming finite elements of any order, for any number of energy groups. The spatial error distribution is assessed with a generalization of an error estimator originally derived for the Poisson equation. Our implementation of this algorithm is based on the widely used Open Source adaptive finite element library deal.II and is made available as part of this library\\'s extensively documented tutorial. We illustrate our methods with results for 2-D and 3-D reactor simulations using 2 and 7 energy groups, and using conforming finite elements of polynomial degree up to 6. © 2008 Elsevier Ltd. All rights reserved.
FEM-2D, 2-D MultiGroup Diffusion in X-Y Geometry
1 - Nature of physical problem solved: FEM-2D solves the two-dimensional diffusion equation in x-y geometry. This is done by the finite elements method. 2 - Method of solution: FEM-2D uses triangular elements with first and second order Lagrange approximations. The systems equations are formulated in multigroup form and solved by Cholesky procedure which operates only on nonzero elements. Various acceleration techniques are available for the outer iteration. Fluxes along various lines and rates in arbitrary zones may be output. 3 - Restrictions on the complexity of the problem: The code uses variable dimensioning. Thus, the problem size is restricted by the largest array which usually is the systems matrix. Fluxes of all groups are kept in memory. This might become another restrictive data set for a large number of groups. The validity of the results is restricted by the approximations used. FEM-2D requires a finite element net which allows the approximation of fluxes by at most parabolas. The node distribution should be more dense in areas of heavy flux changes (near absorbers or the reflector)
Multigroup radiation hydrodynamics with flux-limited diffusion and adaptive mesh refinement
González, Matthias; Commerçon, Benoît; Masson, Jacques
2015-01-01
Radiative transfer plays a key role in the star formation process. Due to a high computational cost, radiation-hydrodynamics simulations performed up to now have mainly been carried out in the grey approximation. In recent years, multi-frequency radiation-hydrodynamics models have started to emerge, in an attempt to better account for the large variations of opacities as a function of frequency. We wish to develop an efficient multigroup algorithm for the adaptive mesh refinement code RAMSES which is suited to heavy proto-stellar collapse calculations. Due to prohibitive timestep constraints of an explicit radiative transfer method, we constructed a time-implicit solver based on a stabilised bi-conjugate gradient algorithm, and implemented it in RAMSES under the flux-limited diffusion approximation. We present a series of tests which demonstrate the high performance of our scheme in dealing with frequency-dependent radiation-hydrodynamic flows. We also present a preliminary simulation of a three-dimensional p...
The Analytic Coarse-Mesh Finite-Difference method is developed in detail for multi-group and multi-dimensional diffusion calculations, including the general and particular modal solutions in the complex space for any number of groups. For rectangular multidimensional geometries, the Chao's generalized relations with transverse integration provide a high-order approximation of the ACMFD method, where all energy groups are coupled by matrix-vector FD relations and the errors are limited to the ones incurred by the interpolation of the transverse interface currents, in a non-linear iterative scheme. The implementation of the method in a multigroup 3D rectangular geometry nodal solver called ANDES is discussed, pointing out the encapsulation achieved for integration of the solver as an optional module within larger code systems. The performance of the ANDES solver in 3D rectangular (X-Y-Z) geometry and multi-groups is verified by its application to several 2D-3D model and international benchmarks (NEA-OECD), with given diffusion cross section sets in few-groups (2 to 8). The extensive verification, always required for new methods and codes, shows a quite fast convergence of ANDES in both the eigenvalue and transverse leakage iteration loops and with the nodal coarse-mesh size, allowing to reach the conclusion that quite high accuracy is achieved with rather large nodes, one node or four nodes per PWR fuel assembly, as compared with reference solutions obtained with fine-mesh finite-difference diffusion calculations using mesh sizes 64 to 128 times smaller than the ANDES nodes. (authors)
Highlights: ► We develop a 2-D, multigroup neutron/adjoint diffusion computer code based on GFEM. ► The spatial discretization is performed using unstructured triangle elements. ► Multiplication factor, flux/adjoint and power distribution are outputs of the code. ► Sensitivity analysis to the number, arrangement and size of elements is performed. ► We proved that the developed code is a reliable tool to solve diffusion equation. -- Abstract: Various methods for solving the forward/adjoint equation in hexagonal and rectangular geometries are known in the literatures. In this paper, the solution of multigroup forward/adjoint equation using Finite Element Method (FEM) for hexagonal and rectangular reactor cores is reported. The spatial discretization of equations is based on Galerkin FEM (GFEM) using unstructured triangle elements. Calculations are performed for both linear and quadratic approximations of the shape function; based on which results are compared. Using power iteration method for the forward and adjoint calculations, the forward and adjoint fluxes with the corresponding eigenvalues are obtained. The results are then benchmarked against the valid results for IAEA-2D, BIBLIS-2D and IAEA-PWR benchmark problems. Convergence rate of GFEM in linear and quadratic approximations of the shape function are calculated and results are quantitatively compared. A sensitivity analysis of the calculations to the number and arrangement of elements has been performed.
Shestakov, A I; Offner, S R
2007-03-02
We present a scheme to solve the nonlinear multigroup radiation diffusion (MGD) equations. The method is incorporated into a massively parallel, multidimensional, Eulerian radiation-hydrodynamic code with adaptive mesh refinement (AMR). The patch-based AMR algorithm refines in both space and time creating a hierarchy of levels, coarsest to finest. The physics modules are time-advanced using operator splitting. On each level, separate 'level-solve' packages advance the modules. Our multigroup level-solve adapts an implicit procedure which leads to a two-step iterative scheme that alternates between elliptic solves for each group with intra-cell group coupling. For robustness, we introduce pseudo transient continuation ({Psi}tc). We analyze the magnitude of the {Psi}tc parameter to ensure positivity of the resulting linear system, diagonal dominance and convergence of the two-step scheme. For AMR, a level defines a subdomain for refinement. For diffusive processes such as MGD, the refined level uses Dirichet boundary data at the coarse-fine interface and the data is derived from the coarse level solution. After advancing on the fine level, an additional procedure, the sync-solve (SS), is required in order to enforce conservation. The MGD SS reduces to an elliptic solve on a combined grid for a system of G equations, where G is the number of groups. We adapt the 'partial temperature' scheme for the SS; hence, we reuse the infrastructure developed for scalar equations. Results are presented. We consider a multigroup test problem with a known analytic solution. We demonstrate utility of {Psi}tc by running with increasingly larger timesteps. Lastly, we simulate the sudden release of energy Y inside an Al sphere (r = 15 cm) suspended in air at STP. For Y = 11 kT, we find that gray radiation diffusion and MGD produce similar results. However, if Y = 1 MT, the two packages yield different results. Our large Y simulation contradicts a long-standing theory
VARI-QUIR-3, 2-D Multigroup Steady-State Neutron Diffusion in X-Y R-Z or R-Theta Geometry
1 - Nature of physical problem solved: The steady-state, multigroup, two-dimensional neutron diffusion equations are solved in x-y, r-z, and r-theta geometry. 2 - Method of solution: A Gauss-Seidel type of solution with inner and outer iterations is used. The source is held constant during the inner iterations
A multigroup diffusion theory code, TRIHEX-3D, has been developed for hexagonal lattice core analyses. For 2-D problems one can use hexagonal or triangular centre-mesh finite difference (FD) schemes. The geometrical description of the problem is for hexagonal geometry only. Subdivision of each hexagon into uniform triangles is facilitated by a built-in auto-triangularisati on procedure. One can analyse any symmetric part of the core or the whole core as well. Reflective (30deg, 60deg, 90deg, 120deg and 180deg) and rotational (60deg, 120deg and 180deg) symmetry boundary conditions are allowed. For 3-D problems one can use a direct 3-D FDM or an axial flux synthesis method. TRIHEX-3D can be used for the core design problems of VVER type of hexagonal lattice cores. The code has been validated against a LMFBR SNR-300 benchmark problem. (author). 8 tabs., 9 figs., 9 refs., 5 appendixes
EQUIVA, Few-Group Diffusion Parameter for PWR Reflector Region by 1-D Transport Calculation
1 - Description of program or function: EQUIVA-1 generates few-group equivalent diffusion theory parameters for pressurized water (PWR) reflector regions from the results of simple one-dimensional (slab) multigroup transport calculations. The one-dimensional 'normalized generalized equivalence theory' (NGET) method is used for this purpose. Equivalent parameters can be generated for any number of condensed energy groups and for different material regions, such as for the explicit baffle and water reflector of the radial reflector of a PWR, or for the homogenized baffle/reflector. EQUIVA-2 generates environment-insensitive equivalent diffusion theory parameters for pressurized water (PWR) reflector regions from the results of simple one-dimensional (slab) multigroup transport calculations. The one-dimensional reflector models have been implemented, namely the few-group NGET-RM method and the two-group KOEBKE-RM method. Equivalent parameters can be generated for one homogenized region only, although this region itself may be constituted by any number of component regions. 2 - Method of solution: EQUIVA-1: Analytic functions of non-symmetric real matrices are computed by means of a spectral analysis method. EQUIVA-2: Analytic functions of non-symmetric real matrices are computed by means of a spectral analysis method. Component region slab response matrices are combined using rigorous addition rules to obtain a global response matrix for a single region/node. 3 - Restrictions on the complexity of the problem: EQUIVA-1.1 is a variably dimensioned code and there is no restriction on the number of energy groups, etc. The size of the problem is restricted only by the computer core storage available. EQUIVA-2.0 is a variably dimensioned code and there is no restriction on the number of energy groups, etc. The size of the problem is restricted only by the computer core storage available
The objective of this work is to obtain an analytical solution of the neutron diffusion kinetic equation in one-dimensional cartesian geometry, to monoenergetic and multigroup problems. These equations are of the type stiff, due to large differences in the orders of magnitude of the time scales of the physical phenomena involved, which make them difficult to solve. The basic idea of the proposed method is applying the spectral expansion in the scalar flux and in the precursor concentration, taking moments and solving the resulting matrix problem by the Laplace transform technique. Bearing in mind that the equation for the precursor concentration is a first order linear differential equation in the time variable, to enable the application of the spectral method we introduce a fictitious diffusion term multiplied by a positive value which tends to zero. This procedure opened the possibility to find an analytical solution to the problem studied. We report numerical simulations and analysis of the results obtained with the precision controlled by the truncation order of the series. (author)
Prinsloo, Rian Hendrik
2006-01-01
Nodal diffusion methods have been used extensively in nuclear reactor calculations specifically for their performance advantage, but also their superior accuracy. In this work a nodal diffusion method is developed for three-dimensional cylindrical geometry. Recent developments in the Pebble Bed Modular Reactor (PBMR) project have sparked renewed interest in the application of different modelling methods to its design, and naturally included in this effort is a nodal method for ...
Today, we can use a computer cluster consist of a few hundreds CPUs with reasonable budget. Such computer system enables us to do detailed modeling of reactor core. The detailed modeling will improve the safety and the economics of a nuclear reactor by eliminating un-necessary conservatism or missing consideration. To take advantage of such a cluster computer, efficient parallel algorithms must be developed. Mechanical structure analysis community has studied the domain decomposition method to solve the stress-strain equation using the finite element methods. One of the most successful domain decomposition method in terms of robustness is FETI-DP. We have modified the original FETI-DP to solve the eigenvalue problem for the multi-group diffusion problem in previous study. In this study, we report the result of recent modification to handle the three-dimensional subdomain partitioning, and the sub-domain multi-group problem. Modified FETI-DP algorithm has been successfully applied for the eigenvalue problem of multi-group neutron diffusion equation. The overall CPU time is decreasing as number of sub-domains (partitions) is increasing. However, there may be a limit in decrement due to increment of the number of primal points will increase the CPU time spent by the solution of the global equation. Even distribution of computational load (criterion a) is important to achieve fast computation. The subdomain partition can be effectively performed using suitable graph theory partition package such as MeTIS
Highlights: ► The multi-group IDE-NDK was solved numerically in 2D-Cartesian geometry. ► The progressive basic polynomial (BPn) methods showed no numerical oscillations. ► The BP2 algorithm showed good accuracy and efficiency. -- Abstract: The multi-group time-integro-differential equations of the neutron diffusion kinetics (IDE-NDK) was solved numerically in 2D Cartesian geometry with the use of the basic-progressive polynomial approximation (BPn). Two applications were computed: a ramp, and an instantaneous change of the thermal removal macroscopic cross sections of the driver material of the 2D-TWGL benchmark problems. The BP2 algorithm showed good accuracy when compared with the results of other codes. BPn did not show the undesirable oscillations that appeared in other codes.
Benchmarking with the multigroup diffusion high-order response matrix method
The benchmarking capabilities of the high-order response matrix eigenvalue method, which was developed more than a decade ago, are demonstrated by means of the numerical analysis of a variety of two-dimensional Cartesian geometry light-water reactor test problems. These problems are typical of those generally used for the benchmarking of coarse-mesh (nodal) diffusion methods and the numerical results show that the high-order response matrix eigenvalue method is well suited to be used as an alternative to fine-mesh finite-difference and refined mesh nodal methods for the purpose of generating reference solutions to such problems. (author)
In this work we address the development and implementation of the analytic coarse-mesh finite-difference (ACMFD) method in a nodal neutron diffusion solver called ANDES. The first version of the solver is implemented in any number of neutron energy groups, and in 3D Cartesian geometries; thus it mainly addresses PWR and BWR core simulations. The details about the generalization to multigroups and 3D, as well as the implementation of the method are given. The transverse integration procedure is the scheme chosen to extend the ACMFD formulation to multidimensional problems. The role of the transverse leakage treatment in the accuracy of the nodal solutions is analyzed in detail: the involved assumptions, the limitations of the method in terms of nodal width, the alternative approaches to implement the transverse leakage terms in nodal methods - implicit or explicit -, and the error assessment due to transverse integration. A new approach for solving the control rod 'cusping' problem, based on the direct application of the ACMFD method, is also developed and implemented in ANDES. The solver architecture turns ANDES into an user-friendly, modular and easily linkable tool, as required to be integrated into common software platforms for multi-scale and multi-physics simulations. ANDES can be used either as a stand-alone nodal code or as a solver to accelerate the convergence of whole core pin-by-pin code systems. The verification and performance of the solver are demonstrated using both proof-of-principle test cases and well-referenced international benchmarks
Development of 3D multi-group neutron diffusion code for hexagonal geometry
Based on the theory of new flux expansion nodal method to solve the neutron diffusion equations, the intra-nodal fluence rate distribution was expanded in a series of analytic basic functions for each group. In order to improve the accuracy of calculation result, continuities of neutron fluence rate and current were utilized across the nodal surfaces. According to the boundary conditions, the iteration method was adopted to solve the diffusion equation, where inner iteration speedup method is Gauss-Seidel method and outer is Lyusternik-Wagner. A new speedup method (one-outer-iteration and multi-inner-iteration method) was proposed according to the characteristic that the convergence speed of multiplication factor is faster than that of neutron fluence rate and the update of inner iteration matrix is slow. Based on the proposed model, the code HANDF-D was developed and tested by 3D two-group vver440 benchmark, experiment 2 of HFETR, 3D four-group thermal reactor benchmark, and 3D seven-group fast reactor benchmark. The numerical results show that HANDF-D can predict accurately the multiplication factor and nodal powers. (authors)
Improvement of the axial diffusion solver of DeCART employing 1-D transport solution
Instead of the 3-D transport equation, DeCART solves a transverse leakage coupled radial transport and axial diffusion equations to obtain an approximate 3-D transport solution. In this paper, two of the approximations contained in DeCART related with diffusion constants and cell homogenization are exmained, and practical improvement schemes are suggested. To overcome the diffusion approximation used in the axial direction, a current conservation scheme based on the axial 1-D transport solution is introduced. To overcome the cell homogenization effect, a plane height refinement scheme is employed near the axial core/reflector boundary where homogenization constants vary significantly in the axial direction within the plane. These schemes are evaluated by solving the 3-D VENUS-2 MOX core benchmark. The current conservation and plane height refinement schemes bring about 280 pcm and 100 pcm improvement in k-eff, respectively, and about 390 pcm in total, but trivial effects in the power distribution
A variational nodal expansion method for the solution of multigroup neutron diffusion equations
An accurate neutronics analysis method is needed for light water reactor core monitoring systems to efficiently operate the core with a smaller margin to limiting parameters. It is also required in in-core fuel management systems to optimize the core loading patterns, and the fuel designs with a higher reliability. When mixed oxide fuel or much higher burnup fuel is used, a new higher order nodal method seems necessary to introduce. Based on these considerations, a new nodal diffusion method for the neutronics analysis of light water reactor cores has been developed. The method is based on an approximation of neutron fluxes by expanding them with a set of functions defined within a node. The expansion coefficients are determined in such a way that the solution becomes the most accurate approximation to the exact solution by utilizing the variational principle. The expansion functions are obtained only from single assembly diffusion calculations. The present method includes no homogenization procedure, and the assembly heterogeneity effect on neutron fluxes is taken into account in a consistent way. The intra-nodal pin-power distribution can also be determined in a consistent way with high accuracy. The present method was implemented in a two-dimensional nodal code, and tested for benchmark cases. The results proved that the accuracy of the present method was excellent. The root mean square errors of both nodal powers and nodal maximum pin powers were observed to be less than 1%. The computing time of the code was measured to be about 3% of the reference, fine-mesh calculation. A three-dimensional version is currently being developed, and since the heterogeneity effect is of less importance in axial direction, a more efficient calculation method can be adopted for the axial solution of the neutron flux. The new method can be used as a ''plug-in'' module to existing core simulators to increase the accuracy of the neutronics part of existing core models, including the
HEXPEDITE: A net current multigroup nodal diffusion method for hexagonal-z geometry
The feasibility of a nodal diffusion algorithm for hexagonal cores was first demonstrated by Duracz and by Lawrence. They implemented a polynomial method with partial currents for internode coupling. Following them, several authors introduced variants of the expansion technique. Wagner developed an analytical method; however, like all previous authors, he still used partial currents for internode coupling and a response matrix solution approach. Very recently, another polynomial model with net currents expressed in terms of transverse-integrated fluxes and a nodal integral method based on coordinate transformations were presented. A transformation-group method was also introduced. In this paper, a hexagonal-z method similar in approach to that of the Cartesian geometry ILLICO is presented. The new method uses an analytical solution of the transverse-integrated equations, net currents for internode coupling, and a global coupling solution scheme different from that of the methods discussed earlier. An extension that treats explicitly the in-node spatial dependence of cross sections is also introduced
Al-Chalabi, Rifat M. Khalil
1997-09-01
Development of an improvement to the computational efficiency of the existing nested iterative solution strategy of the Nodal Exapansion Method (NEM) nodal based neutron diffusion code NESTLE is presented. The improvement in the solution strategy is the result of developing a multilevel acceleration scheme that does not suffer from the numerical stalling associated with a number of iterative solution methods. The acceleration scheme is based on the multigrid method, which is specifically adapted for incorporation into the NEM nonlinear iterative strategy. This scheme optimizes the computational interplay between the spatial discretization and the NEM nonlinear iterative solution process through the use of the multigrid method. The combination of the NEM nodal method, calculation of the homogenized, neutron nodal balance coefficients (i.e. restriction operator), efficient underlying smoothing algorithm (power method of NESTLE), and the finer mesh reconstruction algorithm (i.e. prolongation operator), all operating on a sequence of coarser spatial nodes, constitutes the multilevel acceleration scheme employed in this research. Two implementations of the multigrid method into the NESTLE code were examined; the Imbedded NEM Strategy and the Imbedded CMFD Strategy. The main difference in implementation between the two methods is that in the Imbedded NEM Strategy, the NEM solution is required at every MG level. Numerical tests have shown that the Imbedded NEM Strategy suffers from divergence at coarse- grid levels, hence all the results for the different benchmarks presented here were obtained using the Imbedded CMFD Strategy. The novelties in the developed MG method are as follows: the formulation of the restriction and prolongation operators, and the selection of the relaxation method. The restriction operator utilizes a variation of the reactor physics, consistent homogenization technique. The prolongation operator is based upon a variant of the pin power
The report documents the computer code block VENTURE designed to solve multigroup neutronics problems with application of the finite-difference diffusion-theory approximation to neutron transport (or alternatively simple P1) in up to three-dimensional geometry. It uses and generates interface data files adopted in the cooperative effort sponsored by the Reactor Physics Branch of the Division of Reactor Research and Development of the Energy Research and Development Administration. Several different data handling procedures have been incorporated to provide considerable flexibility; it is possible to solve a wide variety of problems on a variety of computer configurations relatively efficiently
In the present contribution we discuss the solution of the two-dimensional multi-group neutron kinetic equation in cylindrical geometry. The solution is obtained in analytical representation. To this end the scalar flux is extended in terms of the eigenfunctions associated to the respective problem in Cartesian geometry. Taking moments and using orthogonality properties of the eigenfunctions we get a matrix differential equation for the expansion coefficients which has a known solution. We apply this methodology for the neutron kinetic diffusion equation and show numerical results for two-energy groups. (author)
This report presents the HEXAGA-III-programme solving multi-group time-independent real and/or adjoint neutron diffusion equations for three-dimensional-triangular-z-geometry. The method of solution is based on the AGA two-sweep iterative method belonging to the family of factorization techniques. An arbitrary neutron scattering model is permitted. The report written for users provides the description of the programme input and output and the use of HEXAGA-III is illustrated by a sample reactor problem. (orig.)
The computer code block VENTURE, designed to solve multigroup neutronics problems with application of the finite-difference diffusion-theory approximation to neutron transport (or alternatively simple P1) in up to three-dimensional geometry is described. A variety of types of problems may be solved: the usual eigenvalue problem, a direct criticality search on the buckling, on a reciprocal velocity absorber (prompt mode), or on nuclide concentrations, or an indirect criticality search on nuclide concentrations, or on dimensions. First-order perturbation analysis capability is available at the macroscopic cross section level
Vondy, D.R.; Fowler, T.B.; Cunningham, G.W.
1977-11-01
The report documents the computer code block VENTURE designed to solve multigroup neutronics problems with application of the finite-difference diffusion-theory approximation to neutron transport (or alternatively simple P/sub 1/) in up to three-dimensional geometry. It uses and generates interface data files adopted in the cooperative effort sponsored by the Reactor Physics Branch of the Division of Reactor Research and Development of the Energy Research and Development Administration. Several different data handling procedures have been incorporated to provide considerable flexibility; it is possible to solve a wide variety of problems on a variety of computer configurations relatively efficiently.
Oliveira, F.R.; Vilhena, Marco T.; Bodmann, B.E.J., E-mail: fernando.rodrigues@ufrgs.br, E-mail: bardo.bodmann@ufrgs.br [Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, RS (Brazil); Carvalho, F., E-mail: fernando@nuclear.ufrj.br [Coordenacao dos Cursos de Pos-Graduacao em Engenharia (COPPE/UFRJ), Rio de Janeiro, RJ (Brazil). Institute Alberto Luiz Coimbra
2015-07-01
In the present contribution we discuss the solution of the two-dimensional multi-group neutron kinetic equation in cylindrical geometry. The solution is obtained in analytical representation. To this end the scalar flux is extended in terms of the eigenfunctions associated to the respective problem in Cartesian geometry. Taking moments and using orthogonality properties of the eigenfunctions we get a matrix differential equation for the expansion coefficients which has a known solution. We apply this methodology for the neutron kinetic diffusion equation and show numerical results for two-energy groups. (author)
The multi-group integro-differential equations of the neutron diffusion kinetics (IDE-NDK) was presented and solved numerically in multi-slab geometry with the use of the progressive polynomial approximation. Four applications were computed: a positive ramp, a negative ramp, a sinusoidal and an instantaneous change of thermal macroscopic cross-sections in an 120 slab-nuclear reactor for a 2 prompt-group model. The results showed good accuracy for the developed non-iterative algorithms. It was shown the advantage of using the IDE-NDK over the traditional partial differential equations of the neutron diffusion kinetics from an accuracy point of view. Finite difference algorithms were also developed to obtain initial conditions and to make desired comparisons.
Quantum Diffusion on Molecular Tubes: Universal Scaling of the 1D to 2D Transition
Chuang, Chern; Lee, Chee Kong; Moix, Jeremy M.; Knoester, Jasper; Cao, Jianshu
2016-05-01
The transport properties of disordered systems are known to depend critically on dimensionality. We study the diffusion coefficient of a quantum particle confined to a lattice on the surface of a tube, where it scales between the 1D and 2D limits. It is found that the scaling relation is universal and independent of the temperature, disorder, and noise parameters, and the essential order parameter is the ratio between the localization length in 2D and the circumference of the tube. Phenomenological and quantitative expressions for transport properties as functions of disorder and noise are obtained and applied to real systems: In the natural chlorosomes found in light-harvesting bacteria the exciton transfer dynamics is predicted to be in the 2D limit, whereas a family of synthetic molecular aggregates is found to be in the homogeneous limit and is independent of dimensionality.
The energy spectra of fast and thermal neutrons from fission reactions in the FZJ code TINTE are modelled by two broad energy groups. Present demands for increased numerical accuracy led to the question of how precise the 2-group approximation is compared to a multi-group model. Therefore a new simulation program called MGT (Multi Group TINTE) has recently been developed which is able to handle up to 43 energy groups. Furthermore, an internal spectrum calculation for the determination of cross-sections can be performed for each time step and location within the reactor. In this study the multi-group energy models are compared to former calculations with only two energy groups. Different scenarios (normal operation and design-basis accidents) have been defined for a high temperature pebble bed reactor design with annular core. The effect of an increasing number of energy groups on safety-related parameters like the fuel and coolant temperature, the nuclear heat source or the xenon concentration is studied. It has been found that for the studied scenarios the use of up to 8 energy groups is a good trade-off between precision and a tolerable amount of computing time. (orig.)
Methodology for 3-D calculation analysis of nuclear reactor cell with axial symmetry and finite mesh step is described. This methodology is based on the axial leakage calculation analysis method that has been developed for nuclear reactor with closed lattice like VVER-type. The trial functions that are used at full core level of nuclear reactor calculation analysis are defined. Connection between core reactor equation and the definition of trial functions is given. Importance of different trial functions from the point of view the full reactor core calculation is analyzed. If we deal with the case when reactor has strong neutron flux gradients caused with regularization rods it is important to take into account the influence of neutron spectrum into axial leakage. So this paper focuses upon just multi-group approach to obtain matrixes that are defined with trial functions values and with boundary conditions. Previous numerical results of comparison of the matrixes elements analytically obtained and matrix elements obtained with described methodology are given. Analytical expressions for two-group matrix elements are considered as verification results for multi-group numerical scheme. (authors)
A modular computer code system called FEMSYN has been developed to solve the multigroup diffusion theory equations. The various methods that are incorporated in FEMSYN are (i) finite difference method (FDM) (ii) finite element method (FEM) and (iii) single channel flux synthesis method (SCFS). These methods are described in detail in parts II, III and IV of the present report. In this report, a comparison of the accuracy and the speed of different methods of solution for some benchmark problems are reported. The input preparation and listing of sample input and output are included in the Appendices. The code FEMSYN has been used to solve a wide variety of reactor core problems. It can be used for both LWR and PHWR applications. (author)
Ceolin, C.; Schramm, M.; Vilhena, M.T.; Bodmann, B.E.J., E-mail: celina.ceolin@gmail.com, E-mail: marceloschramm@hotmail.com, E-mail: vilhena@pq.cnpq.br, E-mail: bardo.bodmann@ufrgs.br [Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, RS (Brazil). Programa de Pos-Graduacao em Engenharia Mecanica
2013-07-01
In this work the authors solved the steady state neutron diffusion equation for a multi-layer slab assuming the multi-group energy model. The method to solve the equation system is based on a expansion in Taylor Series, which was proven to be useful in [1] [2] [3]. The results obtained can be used as initial condition for neutron space kinetics problems. The neutron scalar flux was expanded in a power series, and the coefficients were found by using the ordinary differential equation and the boundary and interface conditions. The effective multiplication factor k was evaluated using the power method [4]. We divided the domain into several slabs to guarantee the convergence with a low truncation order. We present the formalism together with some numerical simulations. (author)
Ceolin, Celina; Schramm, Marcelo; Bodmann, Bardo Ernst Josef; Vilhena, Marco Tullio Mena Barreto de [Universidade Federal do Rio Grande do Sul, Porto Alegre (Brazil). Programa de Pos-Graduacao em Engenharia Mecanica; Bogado Leite, Sergio de Queiroz [Comissao Nacional de Energia Nuclear, Rio de Janeiro (Brazil)
2014-11-15
In this work the authors solved the steady state neutron diffusion equation for a multi-layer slab assuming the multi-group energy model. The method to solve the equation system is based on an expansion in Taylor Series resulting in an analytical expression. The results obtained can be used as initial condition for neutron space kinetics problems. The neutron scalar flux was expanded in a power series, and the coefficients were found by using the ordinary differential equation and the boundary and interface conditions. The effective multiplication factor k was evaluated using the power method. We divided the domain into several slabs to guarantee the convergence with a low truncation order. We present the formalism together with some numerical simulations.
1 - Description of program or function: DIF3D solves multigroup diffusion theory eigenvalue, adjoint, fixed source and criticality (concentration search) problems in 1-, 2- and 3-space dimensions for orthogonal (rectangular or cylindrical), triangular and hexagonal geometries. Anisotropic diffusion coefficients are permitted. Flux and power density maps by mesh cell and region-wise balance integrals are provided. Although primarily designed for fast reactor problems, up-scattering and internal black boundary conditions are also treated. The DIF3D8.0/VARIANT8.0 release differs from the previous DIF3D7.0 release in that it includes a significantly expanded set of solution techniques using variational nodal methods. DIF3D's nodal option solves the multigroup steady state neutron diffusion equation in two- and three-dimensional hexagonal and cartesian geometries and solves the transport equation in two-and three-dimensional cartesian geometries. Eigenvalue, adjoint, fixed source and criticality (concentration) search problems are permitted as are anisotropic diffusion coefficients. Flux and power density maps by mesh cell and region-wise balance integrals are provided. Although primarily designed for fast reactor problems, up-scattering and for finite difference option only internal black boundary conditions are also treated. VARIANT solves the multigroup steady-state neutron diffusion and transport equations in two- and three-dimensional Cartesian and hexagonal geometries using variational nodal methods. The transport approximations involve complete spherical harmonic expansions up to order P5. Eigenvalue, adjoint, fixed source, gamma heating, and criticality (concentration) search problems are permitted. Anisotropic scattering is treated, and although primarily designed for fast reactor problems, up-scattering options are also included. Related and Auxiliary Programs: DIF3D reads and writes the standard interface files specified by the Committee on Computer Code
Retroviral intasomes search for a target DNA by 1D diffusion which rarely results in integration.
Jones, Nathan D; Lopez, Miguel A; Hanne, Jeungphill; Peake, Mitchell B; Lee, Jong-Bong; Fishel, Richard; Yoder, Kristine E
2016-01-01
Retroviruses must integrate their linear viral cDNA into the host genome for a productive infection. Integration is catalysed by the retrovirus-encoded integrase (IN), which forms a tetramer or octamer complex with the viral cDNA long terminal repeat (LTR) ends termed an intasome. IN removes two 3'-nucleotides from both LTR ends and catalyses strand transfer of the recessed 3'-hydroxyls into the target DNA separated by 4-6 bp. Host DNA repair restores the resulting 5'-Flap and single-stranded DNA (ssDNA) gap. Here we have used multiple single molecule imaging tools to determine that the prototype foamy virus (PFV) retroviral intasome searches for an integration site by one-dimensional (1D) rotation-coupled diffusion along DNA. Once a target site is identified, the time between PFV strand transfer events is 470 ms. The majority of PFV intasome search events were non-productive. These observations identify new dynamic IN functions and suggest that target site-selection limits retroviral integration. PMID:27108531
More accurate modelling of physical phenomena involved in present and future nuclear reactors requires a multi-scale and multi-physics approach. This challenge can be accomplished by the coupling of best-estimate core-physics, thermal-hydraulics and multi-physics solvers. In order to make viable that coupling, the current trends in reactor simulations are along the development of a new generation of tools based on user-friendly, modular, easily linkable, faster and more accurate codes to be integrated in common platforms. These premises are in the origin of the NURESIM Integrated Project within the 6. European Framework Program, which is envisaged to provide the initial step towards a Common European Standard Software Platform for nuclear reactors simulations. In the frame of this project and to reach the above-mentioned goals, a 3-D multigroup nodal solver for neutron diffusion calculations called ANDES (Analytic Nodal Diffusion Equation Solver) has been developed and tested in-depth in this Thesis. ANDES solves the steady-state and time-dependent neutron diffusion equation in three-dimensions and any number of energy groups, utilizing the Analytic Coarse-Mesh Finite-Difference (ACMFD) scheme to yield the nodal coupling equations. It can be applied to both Cartesian and triangular-Z geometries, so that simulations of LWR as well as VVER, HTR and fast reactors can be performed. The solver has been implemented in a fully encapsulated way, enabling it as a module to be readily integrated in other codes and platforms. In fact, it can be used either as a stand-alone nodal code or as a solver to accelerate the convergence of whole core pin-by-pin code systems. Verification of performance has shown that ANDES is a code with high order definition for whole core realistic nodal simulations. In this paper, the methodology developed and involved in ANDES is presented. (authors)
Hernandez-Charpak, J.; Hoogeboom-Pot, K.; Anderson, E.; Murnane, M.; Kapteyn, H.; Nardi, D.
2014-03-01
How is thermal transport affected by spatial confinement in nanoscale systems? In past work we and others demonstrated that the Fourier Law of heat diffusion fails for length scales smaller than the mean free path of the energy carriers in a material. Here we probe how the transition from macroscopic diffusive behavior of phonons through the quasi-ballistic regime is different for 1D and 2D nano-confined hot spots. We study a series of periodic nickel lines (1D) and dots (2D) with linewidths varying from 750 to 30 nm deposited on both sapphire and silicon substrates. The thermal relaxation of these femtosecond-laser-excited nanostructures is monitored by the diffraction of extreme ultraviolet (EUV) light obtained from tabletop high harmonic generation. The short wavelength of EUV light, combined with the coherence and ultrashort pulses of high harmonic sources, provides a unique and powerful probe for nanostructured materials on their intrinsic length and time scales. The relaxation dynamics are linked to an effective thermal boundary resistivity with the assistance of multi-physics finite element analysis to quantify the stronger deviation from macroscopic diffusive behavior as a function of nanostructure linewidth in 2D hot spots compared to 1D. This work was supported by SRC Contract 2012-OJ-2304, by NSF Award No.: DGE 1144083, and used facilities provided by the NSF Engineering Research Center in EUV Science and Technology.
Iterative 2-D/1-D methods for the 3-D neutron diffusion calculation
To remedy the problems arising from assembly homogenization and de-homogenization, several efforts have been made to solve directly the heterogeneous problem with a fine mesh and to reduce the computational burden by coupling 2-D planar with 1-D axial solutions using a Transverse Leakage (TL) coupling. However, the potential for a numerical instability at a small axial mesh size has been observed. Lee et al. showed that one of the two existing methods, method A, is mathematically unstable at a small mesh size while the other, method B, is always stable. They also proposed a new method for a 2-D/1-D coupling, method C, and they showed that it is always stable and it provides the best performance in terms of the convergence rate. In this paper another algorithm, method D, is proposed and its stability is also investigated
Very high-order finite volume scheme for 1D convection diffusion problem: the implicit case
Machado, Gaspar; Clain, Stephane; Pereira, Rui
2012-01-01
We present a fourth-order in space, second-order in time finite volume scheme for transient convection diffusion problem based on the Polynomial Reconstruction Operator and a Crank-Nicholson method. A detailed description of the scheme is provided and we perform numerical tests to highlight the performance of the method in comparison with the classical Patankar method.
Ferizi, Uran; Schneider, Torben; Alipoor, Mohammad; Eufracio, Odin; Fick, Rutger H J; Deriche, Rachid; Nilsson, Markus; Loya-Olivas, Ana K; Rivera, Mariano; Poot, Dirk H J; Ramirez-Manzanares, Alonso; Marroquin, Jose L; Rokem, Ariel; Pötter, Christian; Dougherty, Robert F; Sakaie, Ken; Wheeler-Kingshott, Claudia; Warfield, Simon K; Witzel, Thomas; Wald, Lawrence L; Raya, José G; Alexander, Daniel C
2016-01-01
A large number of mathematical models have been proposed to describe the measured signal in diffusion-weighted (DW) magnetic resonance imaging (MRI) and infer properties about the white matter microstructure. However, a head-to-head comparison of DW-MRI models is critically missing in the field. To address this deficiency, we organized the "White Matter Modeling Challenge" during the International Symposium on Biomedical Imaging (ISBI) 2015 conference. This competition aimed at identifying the DW-MRI models that best predict unseen DW data. in vivo DW-MRI data was acquired on the Connectom scanner at the A.A.Martinos Center (Massachusetts General Hospital) using gradients strength of up to 300 mT/m and a broad set of diffusion times. We focused on assessing the DW signal prediction in two regions: the genu in the corpus callosum, where the fibres are relatively straight and parallel, and the fornix, where the configuration of fibres is more complex. The challenge participants had access to three-quarters of t...
A New 2D-Transport, 1D-Diffusion Approximation of the Boltzmann Transport equation
Larsen, Edward
2013-06-17
The work performed in this project consisted of the derivation, implementation, and testing of a new, computationally advantageous approximation to the 3D Boltz- mann transport equation. The solution of the Boltzmann equation is the neutron flux in nuclear reactor cores and shields, but solving this equation is difficult and costly. The new “2D/1D” approximation takes advantage of a special geometric feature of typical 3D reactors to approximate the neutron transport physics in a specific (ax- ial) direction, but not in the other two (radial) directions. The resulting equation is much less expensive to solve computationally, and its solutions are expected to be sufficiently accurate for many practical problems. In this project we formulated the new equation, discretized it using standard methods, developed a stable itera- tion scheme for solving the equation, implemented the new numerical scheme in the MPACT code, and tested the method on several realistic problems. All the hoped- for features of this new approximation were seen. For large, difficult problems, the resulting 2D/1D solution is highly accurate, and is calculated about 100 times faster than a 3D discrete ordinates simulation.
Selected neutron reaction nuclear data libraries and photon-atomic interaction cross section libraries for elements of interest to the IAEA's program on Fusion Evaluated Nuclear Data Library (FENDL) have been processed into MATXSR format using the NJOY system on the VAX4000 computer of the IAEA. This document lists the resulting multigroup data libraries. All the multigroup data generated are available cost-free upon request from the IAEA Nuclear Data Section. (author). 9 refs
Highlights: • New type of multi-level rebalancing approach for nodal transport. • Generally improved and more mesh-independent convergence behavior. • Importance for intended regime of 3D pin-by-pin core computations. - Abstract: A new multi-level surface rebalancing (MLSR) approach has been developed, aimed at enabling an improved non-linear acceleration of nodal flux iteration convergence in 3D steady-state and transient reactor simulation. This development is meant specifically for anticipating computational needs for solving envisaged multi-group diffusion-like SPN calculations with enhanced mesh resolution (i.e. 3D multi-box up to 3D pin-by-pin grid). For the latter grid refinement regime, the previously available multi-level coarse mesh rebalancing (MLCMR) strategy has been observed to become increasingly inefficient with increasing 3D mesh resolution. Furthermore, for very fine 3D grids that feature a very fine axial mesh as well, non-convergence phenomena have been observed to emerge. In the verifications pursued up to now, these problems have been resolved by the new approach. The novelty arises from taking the interface current balance equations defined over all Cartesian box edges, instead of the nodal volume-integrated process-rate balance equation, as an appropriate restriction basis for setting up multi-level acceleration of fine grid interface current iterations. The new restriction strategy calls for the use of a newly derived set of adjoint spectral equations that are needed for computing a limited set of spectral response vectors per node. This enables a straightforward determination of group-condensed interface current spectral coupling operators that are of crucial relevance in the new rebalancing setup. Another novelty in the approach is a new variational method for computing the neutronic eigenvalue. Within this context, the latter is treated as a control parameter for driving another, newly defined and numerically more fundamental
D3D and D3E, branches of a computer program, solve two- and three-dimensional real and ajoint stationary multigroup neutron diffusion equations by approximating the differential equations by finite difference equations. The discrete grid is a mesh edged one, so that the neutron fluxes are calculated on surfaces separating zones to which different physical conditions apply. Different options allow to treat homogeneous, i.e. eigenvalue problems as well as inhomogeneous, i.e. external source driven problems. The linear algebraic system of the difference equations is solved by the outer and inner iterations method. An outer iteration of the homogeneous problem is the power iteration with the fission source, whereas the outer iteration of the inhomogeneous problem is an iteration with the fission source. Within the process of an outer iteration the group fluxes are determined by inner iterations, either via block overrelaxation or a method of conjugate gradients. (orig./HP)
The Suppression of Energy Discretization Errors in Multigroup Transport Calculations
Larsen, Edward
2013-06-17
The Objective of this project is to develop, implement, and test new deterministric methods to solve, as efficiently as possible, multigroup neutron transport problems having an extremely large number of groups. Our approach was to (i) use the standard CMFD method to "coarsen" the space-angle grid, yielding a multigroup diffusion equation, and (ii) use a new multigrid-in-space-and-energy technique to efficiently solve the multigroup diffusion problem. The overall strategy of (i) how to coarsen the spatial and energy grids, and (ii) how to navigate through the various grids, has the goal of minimizing the overall computational effort. This approach yields not only the fine-grid solution, but also coarse-group flux-weighted cross sections that can be used for other related problems.
The Suppression of Energy Discretization Errors in Multigroup Transport Calculations
The Objective of this project is to develop, implement, and test new deterministric methods to solve, as efficiently as possible, multigroup neutron transport problems having an extremely large number of groups. Our approach was to (i) use the standard CMFD method to 'coarsen' the space-angle grid, yielding a multigroup diffusion equation, and (ii) use a new multigrid-in-space-and-energy technique to efficiently solve the multigroup diffusion problem. The overall strategy of (i) how to coarsen the spatial an energy grids, and (ii) how to navigate through the various grids, has the goal of minimizing the overall computational effort. This approach yields not only the fine-grid solution, but also coarse-group flux-weighted cross sections that can be used for other related problems.
AMPX-77, Modular System for Coupled Neutron-Gamma Multigroup Cross-Sections from ENDF/B-5
1 - Description of program or function: The AMPX system is a system of computer programs (modules) capable of producing coupled multigroup neutron-gamma-ray cross section sets. The system is one of the standards for producing multigroup neutron, gamma-ray production, gamma-ray interaction, and coupled neutron-gamma cross-section sets from ENDF data. AMPX-produced cross sections can be used directly with a variety of diffusion theory, discrete ordinates, and Monte Carlo radiation transport computer codes. A one-dimensional Sn calculation capability is provided for general use and for cross section collapsing. Treatments are included for resonance self-shielding effects. 2 - Method of solution: The system includes a full range of features needed to: (1) produce multigroup neutron, gamma-ray production, and/or gamma-ray interaction cross-section data, (2) resonance self-shield, (3) spectrally collapse, (4) convert cross-section libraries from one format to another format, (5) execute a one- dimensional (1-D) discrete-ordinates calculation, and (6) perform miscellaneous cross section-operations. 3 - Restrictions on the complexity of the problem: The principal restriction is the availability of adequate core storage. All large modules are variably dimensioned. Certain modules will automatically use external storage (disk,tape), if in-core storage is inadequate. While these procedures are of little consequence on today's large computers with 'virtual memory' capabilities, they can be important when small-core PC's or workstations are used
Coupling of Nod1D and HOTCHANNEL: static case
In this work the joining of the programs Nod1D and HOTCHANNEL, developed in the National Polytechnic Institute (IPN) and in the Electrical Research Institute (IIE) respectively is described. The first one allows to study the neutronic of a nuclear reactor and the second one allows to carry out the analysis of hot channel of a Boiling Water Reactor (BWR). Nod1 D is a program that it solves by nodal methods type finite element those diffusion equations in multigroup, and it is the static part of Nod Kin that it solves the diffusion equation in their time dependent part. For another side HOTCHANNEL is based on a mathematical model constituted by four conservation equations (two of mass conservation, one of motion quantity and one of energy), which are solved applying one discretization in implicit finite differences. Both programs have been verified in independent form using diverse test problems. In this work the modifications that were necessary to carry out to both for obtaining a coupled program that it provides the axial distribution of the neutron flux, the power, the burnup and the void fraction, among others parameters as much as neutronic as thermal hydraulics are described. Those are also mentioned limitations, advantages and disadvantages of the final product to which has been designated Nod1 D-HotChn. Diverse results for the Cycle 1 of the Laguna Verde Unit 1 reactor of the Nucleo electric central comparing them with those obtained directly with the CoreMasterPresto code are provided. (Author)
Coupling of Nod1D and HOTCHANNEL: static case; Acoplamiento de Nod1D y HOTCHANNEL: caso estatico
Gomez T, A.M. [IPN-ESFM, 07738 Mexico D.F. (Mexico); Ovando C, R. [IIE-Gcia. de Energia Nuclear, Cuernavaca, Morelos (Mexico)]. e-mail: rovando@iie.org.mx
2003-07-01
In this work the joining of the programs Nod1D and HOTCHANNEL, developed in the National Polytechnic Institute (IPN) and in the Electrical Research Institute (IIE) respectively is described. The first one allows to study the neutronic of a nuclear reactor and the second one allows to carry out the analysis of hot channel of a Boiling Water Reactor (BWR). Nod1 D is a program that it solves by nodal methods type finite element those diffusion equations in multigroup, and it is the static part of Nod Kin that it solves the diffusion equation in their time dependent part. For another side HOTCHANNEL is based on a mathematical model constituted by four conservation equations (two of mass conservation, one of motion quantity and one of energy), which are solved applying one discretization in implicit finite differences. Both programs have been verified in independent form using diverse test problems. In this work the modifications that were necessary to carry out to both for obtaining a coupled program that it provides the axial distribution of the neutron flux, the power, the burnup and the void fraction, among others parameters as much as neutronic as thermal hydraulics are described. Those are also mentioned limitations, advantages and disadvantages of the final product to which has been designated Nod1 D-HotChn. Diverse results for the Cycle 1 of the Laguna Verde Unit 1 reactor of the Nucleo electric central comparing them with those obtained directly with the CoreMasterPresto code are provided. (Author)
We study the physical-vapor-deposition of 1D bismuth nanostructures. Bi nanowire elongating along [012] and/or [110] direction as well as anisotropic Bi nano-columns are physical-vapor-deposited successfully. The coexistence and competition of surface diffusion and geometric shielding are critical to their formation as well as growth mode transition among them. Since physical-vapor-deposition is a vacuum process, we make use of it to fabricate the ohmic contact to prevent the damage to the bismuth nanostructures brought by the etching to their thick surface oxide layer. (paper)
The multigroup neutronics model of NuStar's 3D core code EGRET
As a key component of NuStar's core analysis system for PWR application, EGRET is designed to perform steady-state coupled neutronic/hydraulic analysis of PWRs. This paper presents EGRET's unique 3D nodal diffusion model and 2D pin power reconstruction (PPR) model. Unlike the practice in most of today's production codes that iteratively solves the global 3D coarse-mesh problem and the local axially 1D fine-mesh problem to handle the axial heterogeneity within a node caused by fuel grid and partially-inserted control rod, EGRET resolves the issue by inventing a new nodal technology and introducing the adaptive meshing technique to follow the movement of control rod tip. The new nodal method employs fine-mesh heterogeneous calculation with coarse-mesh transverse coupling such that the axial heterogeneous nodes can be explicitly modeled in exact geometry and directly incorporated into the scheme of transversely coupled coarse-mesh nodal methods. Each axial channel can have its own fine-mesh division without the need of dividing the whole core into radially coupled fine-meshes. There is no need to do 1D fine-mesh and 3D coarse-mesh iteration either. While for the PPR model, EGRET adopts a group-decoupled direct fitting method, which avoids both the complication of constructing 2D analytic multigroup flux solution and any group-coupled iteration. Another unique feature of the PPR model is that it fully utilizes all the information available from 3D core calculation into the downstream PPR process. Particularly, for the first time, the 1D profiles of transversely-integrated fluxes are utilized as the additional conditions to reconstruct pin power. Numerical results of series of benchmark problems verify the good performance of EGRET's unique multi-group neutronics model. (author)
Petersen, Claudio Z. [Universidade Federal de Pelotas, Capao do Leao (Brazil). Programa de Pos Graduacao em Modelagem Matematica; Bodmann, Bardo E.J.; Vilhena, Marco T. [Universidade Federal do Rio Grande do Sul, Porto Alegre, RS (Brazil). Programa de Pos-graduacao em Engenharia Mecanica; Barros, Ricardo C. [Universidade do Estado do Rio de Janeiro, Nova Friburgo, RJ (Brazil). Inst. Politecnico
2014-12-15
In the present work we solve in analytical representation the three dimensional neutron kinetic diffusion problem in rectangular Cartesian geometry for homogeneous and bounded domains for any number of energy groups and precursor concentrations. The solution in analytical representation is constructed using a hierarchical procedure, i.e. the original problem is reduced to a problem previously solved by the authors making use of a combination of the spectral method and a recursive decomposition approach. Time dependent absorption cross sections of the thermal energy group are considered with step, ramp and Chebyshev polynomial variations. For these three cases, we present numerical results and discuss convergence properties and compare our results to those available in the literature.
Hinata, Hirofumi; Kataoka, Tomoya
2016-08-15
We propose a belt transect setting strategy for mark-recapture experiments (MREs) to evaluate the time-independent 1D diffusion coefficient (〈Dp0〉) of marine litter in the cross-shore direction that determines the backwashing flux of the litter, based on two-year MREs for plastic floats (PFs) on Wadahama Beach, Nii-jima Island, Japan. When the alongshore width of the belt transect (Lt) was of the order of, or longer than, the length scale of wave-induced nearshore current circulation (Lc), the PFs were rarely transported alongshore across the selected transects prior to being backwashed offshore. Thus, the transect residence time became longer and showed a much weaker dependence on the transect position, in contrast to when Lt was even shorter than Lc. We therefore obtained the diffusion coefficients close to the value of (〈Dp0〉) when we set Lt to the order of, or longer than, Lc. PMID:27263978
A multigroup treatment of radiation transport
A multi-group radiation package is outlined which will accurately handle radiation transfer problems in laser-produced plasmas. Bremsstrahlung, recombination and line radiation are included as well as fast electron Bremsstrahlung radiation. The entire radiation field is divided into a large number of groups (typically 20), which diffuse radiation energy in real space as well as in energy space, the latter occurring via electron-radiation interaction. Using this model a radiation transport code will be developed to be incorporated into MEDUSA. This modified version of MEDUSA will be used to study radiative preheat effects in laser-compression experiments at the Central Laser Facility, Rutherford Laboratory. The model is also relevant to heavy ion fusion studies. (author)
Procedure to Generate the MPACT Multigroup Library
The CASL neutronics simulator MPACT is under development for the neutronics and T-H coupled simulation for the light water reactor. The objective of this document is focused on reviewing the current procedure to generate the MPACT multigroup library. Detailed methodologies and procedures are included in this document for further discussion to improve the MPACT multigroup library.
Procedure to Generate the MPACT Multigroup Library
Kim, Kang Seog [Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
2015-12-17
The CASL neutronics simulator MPACT is under development for the neutronics and T-H coupled simulation for the light water reactor. The objective of this document is focused on reviewing the current procedure to generate the MPACT multigroup library. Detailed methodologies and procedures are included in this document for further discussion to improve the MPACT multigroup library.
Laginha Silva, Patricia; Martins, Flávio A.; Boski, Tomász; Sampath, Dissanayake M. R.
2010-05-01
processes. In this viewpoint the system is broken down into its fundamental components and processes and the model is build up by selecting the important processes regardless of its time and space scale. This viewpoint was only possible to pursue in the recent years due to improvement in system knowledge and computer power (Paola, 2000). The primary aim of this paper is to demonstrate that it is possible to simulate the evolution of the sediment river bed, traditionally studied with synthetic models, with a process-based hydrodynamic, sediment transport and morphodynamic model, solving explicitly the mass and momentum conservation equations. With this objective, a comparison between two mathematical models for alluvial rivers is made to simulate the evolution of the sediment river bed of a conceptual 1D embayment for periods in the order of a thousand years: the traditional synthetic basin infilling aggregate diffusive type model based on the diffusion equation (Paola, 2000), used in the "synthesist" viewpoint and the process-based model MOHID (Miranda et al., 2000). The simulation of the sediment river bed evolution achieved by the process-based model MOHID is very similar to those obtained by the diffusive type model, but more complete due to the complexity of the process-based model. In the MOHID results it is possible to observe a more comprehensive and realistic results because this type of model include processes that is impossible to a synthetic model to describe. At last the combined effect of tide, sea level rise and river discharges was investigated in the process based model. These effects cannot be simulated using the diffusive type model. The results demonstrate the feasibility of using process based models to perform studies in scales of 10000 years. This is an advance relative to the use of synthetic models, enabling the use of variable forcing. REFERENCES • Briggs, L.I. and Pollack, H.N., 1967. Digital model of evaporate sedimentation. Science, 155, 453
The spatial eigenfunction expansion method is used to solve the multigroup time-dependent diffusion equation when the absorption cross-section in the thermal group is a function of time. An expression for the multi region reactor transfer function is obtained. Some numerical results for two energy groups are also presented. (author)
A discretization of the multigroup PN radiative transfer equation on general meshes
Hermeline, F.
2016-05-01
We propose and study a finite volume method of discrete duality type for discretizing the multigroup PN approximation of radiative transfer equation on general meshes. This method is second order-accurate on a very large variety of meshes, stable under a Courant-Friedrichs-Lewy condition and it preserves naturally the diffusion asymptotic limit.
A 3D multigroup transport kinetics code in hexagonal geometry for fast reactor transient analysis
A description of the 3D multigroup diffusion/transport kinetics code HEXNODYN is given and numerical results are reported. HEXNODYN couples time integration by the quasi-static method with space integration by HEXNOD's analytic (diffusion option) or discrete ordinates (transport option) nodal method. An equivalent hexagonal version of the KfK rod ejection problem has been set up to validate the diffusion option by comparison with available 2D diffusion codes. The transport option has been validated by comparison with the diffusion option. Numerical results indicate that the diffusion option may be considered as fully validated while the transport version is at least internally consistent
A new method for the calculation of diffusion coefficients with Monte Carlo
This paper presents a new Monte Carlo-based method for the calculation of diffusion coefficients. One distinctive feature of this method is that it does not resort to the computation of transport cross sections directly, although their functional form is retained. Instead, a special type of tally derived from a deterministic estimate of Fick's Law is used for tallying the total cross section, which is then combined with a set of other standard Monte Carlo tallies. Some properties of this method are presented by means of numerical examples for a multi-group 1-D implementation. Calculated diffusion coefficients are in general good agreement with values obtained by other methods. (author)
A New Method for the Calculation of Diffusion Coefficients with Monte Carlo
Dorval, Eric
2014-06-01
This paper presents a new Monte Carlo-based method for the calculation of diffusion coefficients. One distinctive feature of this method is that it does not resort to the computation of transport cross sections directly, although their functional form is retained. Instead, a special type of tally derived from a deterministic estimate of Fick's Law is used for tallying the total cross section, which is then combined with a set of other standard Monte Carlo tallies. Some properties of this method are presented by means of numerical examples for a multi-group 1-D implementation. Calculated diffusion coefficients are in general good agreement with values obtained by other methods.
A numerical model for multigroup radiation hydrodynamics
We present in this paper a multigroup model for radiation hydrodynamics to account for variations of the gas opacity as a function of frequency. The entropy closure model (M1) is applied to multigroup radiation transfer in a radiation hydrodynamics code. In difference from the previous grey model, we are able to reproduce the crucial effects of frequency-variable gas opacities, a situation omnipresent in physics and astrophysics. We also account for the energy exchange between neighbouring groups which is important in flows with strong velocity divergence. These terms were computed using a finite volume method in the frequency domain. The radiative transfer aspect of the method was first tested separately for global consistency (reversion to grey model) and against a well-established kinetic model through Marshak wave tests with frequency-dependent opacities. Very good agreement between the multigroup M1 and kinetic models was observed in all tests. The successful coupling of the multigroup radiative transfer to the hydrodynamics was then confirmed through a second series of tests. Finally, the model was linked to a database of opacities for a Xe gas in order to simulate realistic multigroup radiative shocks in Xe. The differences with the previous grey models are discussed.
In reactor physics, calculation schemes with deterministic codes are validated with respect to a reference Monte Carlo code. The remaining biases are attributed to the approximations and models induced by the multigroup theory (self-shielding models and expansion of the scattering law using Legendre polynomials) to represent physical phenomena (resonant absorption and scattering anisotropy respectively). This work focuses on the relevance of a polynomial expansion to model the scattering law. Since the outset of reactor physics, the latter has been expanded on a truncated Legendre polynomial basis. However, the transfer cross sections are highly anisotropic, with non-zero values for a very small range of the cosine of the scattering angle. Besides, the finer the energy mesh and the lighter the scattering nucleus, the more exacerbated is the peaked shape of this cross section. As such, the Legendre expansion is less suited to represent the scattering law. Furthermore, this model induces negative values which are non-physical. In this work, various scattering laws are briefly described and the limitations of the existing model are pointed out. Hence, piecewise-constant functions have been used to represent the multigroup scattering cross section. This representation requires a different model for the diffusion source. The discrete ordinates method which is widely employed to solve the transport equation has been adapted. Thus, the finite volume method for angular discretization has been developed and implemented in Paris environment which hosts the Sn solver, Snatch. The angular finite volume method has been compared to the collocation method with Legendre moments to ensure its proper performance. Moreover, unlike the latter, this method is adapted for both the Legendre moments and the piecewise-constant functions representations of the scattering cross section. This hybrid-source method has been validated for different cases: fuel cell in infinite lattice
Chacón Rebollo, Tomás
2015-03-01
This paper introduces a variational multi-scale method where the sub-grid scales are computed by spectral approximations. It is based upon an extension of the spectral theorem to non necessarily self-adjoint elliptic operators that have an associated base of eigenfunctions which are orthonormal in weighted L2 spaces. This allows to element-wise calculate the sub-grid scales by means of the associated spectral expansion. We propose a feasible VMS-spectral method by truncation of this spectral expansion to a finite number of modes. We apply this general framework to the convection-diffusion equation, by analytically computing the family of eigenfunctions. We perform a convergence and error analysis. We also present some numerical tests that show the stability of the method for an odd number of spectral modes, and an improvement of accuracy in the large resolved scales, due to the adding of the sub-grid spectral scales.
1 - Description of problem or function: KENO is a multigroup, Monte Carlo criticality code containing a special geometry package which allows easy description of systems composed of cylinders, spheres, and cuboids (rectangular parallelepipeds) arranged in any order with only one restriction. They cannot be rotated or translated. Each geometrical region must be described as completely enclosing all regions interior to it. For systems not describable using this special geometry package, the program can use the generalized geometry package (GEOM) developed for the O5R Monte Carlo code. It allows any system that can be described by a collection of planes and/or quadratic surfaces, arbitrarily oriented and intersecting in arbitrary fashion. The entire problem can be mocked up in generalized geometry, or one generalized geometry unit or box type can be used alone or in combination with standard KENO units or box types. Rectangular arrays of fissile units are allowed with or without external reflector regions. Output from KENO consists of keff for the system plus an estimate of its standard deviation and the leakage, absorption, and fissions for each energy group plus the totals for all groups. Flux as a function of energy group and region and fission densities as a function of region are optional output. KENO-4: Added features include a neutron balance edit, PICTURE routines to check the input geometry, and a random number sequencing subroutine written in FORTRAN-4. 2 - Method of solution: The scattering treatment used in KENO assumes that the differential neutron scattering cross section can be represented by a P1 Legendre polynomial. Absorption of neutrons in KENO is not allowed. Instead, at each collision point of a neutron tracking history the weight of the neutron is reduced by the absorption probability. When the neutron weight has been reduced below a specified point for the region in which the collision occurs, Russian roulette is played to determine if the
Multigroup albedo method applied to gamma radiation shielding
The Albedo method, when applied to shielding calculations, is characterized by following the radiation through the materials, determining the reflected, absorbed and transmitted fractions of the incident current, independently of flux calculations. The excellent results obtained to neutron shielding cases in which the diffusion approximation could be applied motivated this work, where the method was applied in order to develop a multigroup and multilayered algorithm. A gamma radiation shielding simulation was carried out to a system constituted by three infinite slabs of varied materials and six energy groups. The results obtained by Albedo Method were the same generated by ANISN, a consecrated deterministic nuclear code. Concludingly, this work demonstrates the validity of Albedo Method to gamma radiation shielding analysis through its agreement with the full Transport Equation. (author)
Parallel computation of multigroup reactivity coefficient using iterative method
One of the research activities to support the commercial radioisotope production program is a safety research target irradiation FPM (Fission Product Molybdenum). FPM targets form a tube made of stainless steel in which the nuclear degrees of superimposed high-enriched uranium. FPM irradiation tube is intended to obtain fission. The fission material widely used in the form of kits in the world of nuclear medicine. Irradiation FPM tube reactor core would interfere with performance. One of the disorders comes from changes in flux or reactivity. It is necessary to study a method for calculating safety terrace ongoing configuration changes during the life of the reactor, making the code faster became an absolute necessity. Neutron safety margin for the research reactor can be reused without modification to the calculation of the reactivity of the reactor, so that is an advantage of using perturbation method. The criticality and flux in multigroup diffusion model was calculate at various irradiation positions in some uranium content. This model has a complex computation. Several parallel algorithms with iterative method have been developed for the sparse and big matrix solution. The Black-Red Gauss Seidel Iteration and the power iteration parallel method can be used to solve multigroup diffusion equation system and calculated the criticality and reactivity coeficient. This research was developed code for reactivity calculation which used one of safety analysis with parallel processing. It can be done more quickly and efficiently by utilizing the parallel processing in the multicore computer. This code was applied for the safety limits calculation of irradiated targets FPM with increment Uranium
Establishment of multi-groups atomic parametric database
A method is given to establish multi-groups atomic parametric database for multi-groups radiation transport equation. The equation can be used in calculating the X-ray radiation from plasma. Several methods to check the calculation of the multi-groups database is also given. A 20 groups atomic parametric database of Au element with grid of 20 (plasma density) x 20 (electron temperature) x 20 (photon temperature) is given too
A Note on Multigroup Comparisons Using SAS PROC CALIS
Jones-Farmer, L. Allison; Pitts, Jennifer P.; Rainer, R. Kelly
2008-01-01
Although SAS PROC CALIS is not designed to perform multigroup comparisons, it is believed that SAS can be "tricked" into doing so for groups of equal size. At present, there are no comprehensive examples of the steps involved in performing a multigroup comparison in SAS. The purpose of this article is to illustrate these steps. We demonstrate…
Cross section probability tables in multi-group transport calculations
The use of cross section probability tables in multigroup transport calculations is presented. Emphasis is placed on how probability table parameters are generated in a multigroup cross section processor and how existing transport codes must be modifed to use them. In order to illustrate the accuracy obtained by using probability tables, results are presented for a variety of neutron and photon transport problems
In multigroup calculations of reactivity and sensitivity coefficients, methodical errors can appear if the interdependence of multigroup constants is not taken into account. For this effect to be taken into account, so-called implicit components of the aforementioned values are introduced. A simple technique for computing these values is proposed. It is based on the use of subgroup parameters.
Multigroup neutron dose calculations for proton therapy
We have developed tools for the preparation of coupled multigroup proton/neutron cross section libraries. Our method is to use NJOY to process evaluated nuclear data files for incident particles below 150 MeV and MCNPX to produce data for higher energies. We modified the XSEX3 program of the MCNPX code system to produce Legendre expansions of scattering matrices generated by sampling the physics models that are comparable to the output of the GROUPR routine of NJOY. Our code combines the low and high energy scattering data with user input stopping powers and energy deposition cross sections that we also calculated using MCNPX. Our code also calculates momentum transfer coefficients for the library and optionally applies an energy straggling model to the scattering cross sections and stopping powers. The motivation was initially for deterministic solution of space radiation shielding calculations using Attila, but noting that proton therapy treatment planning may neglect secondary neutron dose assessments because of difficulty and expense, we have also investigated the feasibility of multi group methods for this application. We have shown that multigroup MCNPX solutions for secondary neutron dose compare well with continuous energy solutions and are obtainable with less than half computational cost. This efficiency comparison neglects the cost of preparing the library data, but this becomes negligible when distributed over many multi group calculations. Our deterministic calculations illustrate recognized obstacles that may have to be overcome before discrete ordinates methods can be efficient alternatives for proton therapy neutron dose calculations
Multigroup neutron dose calculations for proton therapy
Kelsey Iv, Charles T [Los Alamos National Laboratory; Prinja, Anil K [Los Alamos National Laboratory
2009-01-01
We have developed tools for the preparation of coupled multigroup proton/neutron cross section libraries. Our method is to use NJOY to process evaluated nuclear data files for incident particles below 150 MeV and MCNPX to produce data for higher energies. We modified the XSEX3 program of the MCNPX code system to produce Legendre expansions of scattering matrices generated by sampling the physics models that are comparable to the output of the GROUPR routine of NJOY. Our code combines the low and high energy scattering data with user input stopping powers and energy deposition cross sections that we also calculated using MCNPX. Our code also calculates momentum transfer coefficients for the library and optionally applies an energy straggling model to the scattering cross sections and stopping powers. The motivation was initially for deterministic solution of space radiation shielding calculations using Attila, but noting that proton therapy treatment planning may neglect secondary neutron dose assessments because of difficulty and expense, we have also investigated the feasibility of multi group methods for this application. We have shown that multigroup MCNPX solutions for secondary neutron dose compare well with continuous energy solutions and are obtainable with less than half computational cost. This efficiency comparison neglects the cost of preparing the library data, but this becomes negligible when distributed over many multi group calculations. Our deterministic calculations illustrate recognized obstacles that may have to be overcome before discrete ordinates methods can be efficient alternatives for proton therapy neutron dose calculations.
Multigroup Free-atom Doppler-broadening Approximation. Theory
Gray, Mark Girard [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2015-11-06
Multigroup cross sections at a one target temperature can be Doppler-broadened to multigroup cross sections at a higher target temperature by matrix multiplication if the group structure suf- ficiently resolves the original temperature continuous energy cross section. Matrix elements are the higher temperature group weighted averages of the integral over the lower temperature group boundaries of the free-atom Doppler-broadening kernel. The results match theory for constant and 1/v multigroup cross sections at 618 lanl group structure resolution.
Modelling and simulations of macroscopic multi-group pedestrian flow
Mahato, Naveen K; Tiwari, Sudarshan
2016-01-01
We consider a multi-group microscopic model for pedestrian flow describing the behaviour of large groups. It is based on an interacting particle system coupled to an eikonal equation. Hydrodynamic multi-group models are derived from the underlying particle system as well as scalar multi-group models. The eikonal equation is used to compute optimal paths for the pedestrians. Particle methods are used to solve the macroscopic equations. Numerical test cases are investigated and the models and, in particular, the resulting evacuation times are compared for a wide range of different parameters.
Multigroup fast fission factor treatment in a thermal reactor lattice
A multigroup procedure for the studies of the fast fission effects in the thermal reactor lattice and the calculation of the fast fission factor was developed. The Monte Carlo method and the multigroup procedure were combined to calculate the fast neutron interaction and backscattering effects in a reactor lattice. A set of probabilities calculated by the Monte Carlo method gives a multigroup spectrum of neutrons coming from the moderator and entering the fuel element. Thus, the assumptions adopted so far in defining and calculating the fast fission factor has been avoided, and a new definition including the backscattering and interaction effects in a reactor lattice have been given. (author)
Multigroup cross section library; WIMS library
The WIMS library has been extensively used in thermal reactor calculations. This multigroup constants library was originally developed from the UKNDL in the late 60's and has been updated in 1986. This library has been distributed with the WIMS-D code by NEA data bank. The references to WIMS library in literature are the 'old' which is the original as developed by the AEA Winfrith and the 'new' which is the current 1986 WIMS library. IAEA has organised a CRP where a new and fully updated WIMS library will soon be available. This paper gives an overview of the definitions of the group constants that go into any basic nuclear data library used for reactor calculations. This paper also outlines the contents of the WIMS library and some of its shortcomings
Multi-group neutron transport theory
Multi-group neutron transport theory. In the paper the general theory of the application of the K. M. Case method to N-group neutron transport theory in plane geometry is given. The eigenfunctions (distributions) for the system of Boltzmann equations have been derived and the completeness theorem has been proved. By means of general solution two examples important for reactor and shielding calculations are given: the solution of a critical and albedo problem for a slab. In both cases the system of singular integral equations for expansion coefficients into a full set of eigenfunction distributions has been reduced to the system of Fredholm-type integral equations. Some results can be applied also to some spherical problems. (author)
Application of equivalence methods on Monte Carlo method based homogenization multi-group constants
The multi-group constants generated via continuous energy Monte Carlo method do not satisfy the equivalence between reference calculation and diffusion calculation applied in reactor core analysis. To the satisfaction of the equivalence theory, general equivalence theory (GET) and super homogenization method (SPH) were applied to the Monte Carlo method based group constants, and a simplified reactor core and C5G7 benchmark were examined with the Monte Carlo constants. The results show that the calculating precision of group constants is improved, and GET and SPH are good candidates for the equivalence treatment of Monte Carlo homogenization. (authors)
A code to calculate multigroup constants for fast neutron reactor
KQCS-2 code is a new improved version of KQCS code, which was designed to calculate multigroup constants for fast neutron reactor. The changes and improvements on KQCS are described in this paper. (author)
Consistent Multigroup Theory Enabling Accurate Course-Group Simulation of Gen IV Reactors
Rahnema, Farzad; Haghighat, Alireza; Ougouag, Abderrafi
2013-11-29
The objective of this proposal is the development of a consistent multi-group theory that accurately accounts for the energy-angle coupling associated with collapsed-group cross sections. This will allow for coarse-group transport and diffusion theory calculations that exhibit continuous energy accuracy and implicitly treat cross- section resonances. This is of particular importance when considering the highly heterogeneous and optically thin reactor designs within the Next Generation Nuclear Plant (NGNP) framework. In such reactors, ignoring the influence of anisotropy in the angular flux on the collapsed cross section, especially at the interface between core and reflector near which control rods are located, results in inaccurate estimates of the rod worth, a serious safety concern. The scope of this project will include the development and verification of a new multi-group theory enabling high-fidelity transport and diffusion calculations in coarse groups, as well as a methodology for the implementation of this method in existing codes. This will allow for a higher accuracy solution of reactor problems while using fewer groups and will reduce the computational expense. The proposed research represents a fundamental advancement in the understanding and improvement of multi- group theory for reactor analysis.
Radiation Transport for Explosive Outflows: A Multigroup Hybrid Monte Carlo Method
Wollaeger, Ryan T; Graziani, Carlo; Couch, Sean M; Jordan, George C; Lamb, Donald Q; Moses, Gregory A
2013-01-01
We explore the application of Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) to radiation transport in strong fluid outflows with structured opacity. The IMC method of Fleck & Cummings is a stochastic computational technique for nonlinear radiation transport. IMC is partially implicit in time and may suffer in efficiency when tracking Monte Carlo particles through optically thick materials. The DDMC method of Densmore accelerates an IMC computation where the domain is diffusive. Recently, Abdikamalov extended IMC and DDMC to multigroup, velocity-dependent neutrino transport with the intent of modeling neutrino dynamics in core-collapse supernovae. Densmore has also formulated a multifrequency extension to the originally grey DDMC method. In this article we rigorously formulate IMC and DDMC over a high-velocity Lagrangian grid for possible application to photon transport in the post-explosion phase of Type Ia supernovae. The method described is suitable for a large variety of non-mono...
A general multigroup formulation of the analytic nodal method
In this paper the theoretical description of an alternative approach to the Analytic Nodal Method is given, in which a full multigroup formulations is developed. This approach differs from the well known QUANDRY approach in three aspects. Firstly, a notation which is more widely used in Quantum Mechanics has been adopted to enable a clear and concise presentation of this multigroup approach. A basis transformation is then used to reduce the directional equations to a scalar form and finally, Green's secondary identity is used to rewrite each of the resulting scalar equations in a form which eventually leads to a response matrix, as opposed to using classical methods to actually solve the coupled multigroup directional equations
Scattering approach to classical quasi-1D transport
Kogan, Eugene
1996-01-01
General dynamical transport of classical particles in disordered quasi-1D samples is viewed in the framework of scattering approach. Simple equation for the transfer-matrix is obtained within this unified picture. In the case of diffusive transport the solution of this equation exactly coincides with the solution of diffusion equation.
Smith, L.A.; Gallmeier, F.X. [Oak Ridge Institute for Science and Energy, TN (United States); Gehin, J.C. [Oak Ridge National Lab., TN (United States)] [and others
1995-05-01
The FOEHN critical experiment was analyzed to validate the use of multigroup cross sections and Oak Ridge National Laboratory neutronics computer codes in the design of the Advanced Neutron Source. The ANSL-V 99-group master cross section library was used for all the calculations. Three different critical configurations were evaluated using the multigroup KENO Monte Carlo transport code, the multigroup DORT discrete ordinates transport code, and the multigroup diffusion theory code VENTURE. The simple configuration consists of only the fuel and control elements with the heavy water reflector. The intermediate configuration includes boron endplates at the upper and lower edges of the fuel element. The complex configuration includes both the boron endplates and components in the reflector. Cross sections were processed using modules from the AMPX system. Both 99-group and 20-group cross sections were created and used in two-dimensional models of the FOEHN experiment. KENO calculations were performed using both 99-group and 20-group cross sections. The DORT and VENTURE calculations were performed using 20-group cross sections. Because the simple and intermediate configurations are azimuthally symmetric, these configurations can be explicitly modeled in R-Z geometry. Since the reflector components cannot be modeled explicitly using the current versions of these codes, three reflector component homogenization schemes were developed and evaluated for the complex configuration. Power density distributions were calculated with KENO using 99-group cross sections and with DORT and VENTURE using 20-group cross sections. The average differences between the measured values and the values calculated with the different computer codes range from 2.45 to 5.74%. The maximum differences between the measured and calculated thermal flux values for the simple and intermediate configurations are {approx} 13%, while the average differences are < 8%.
Multigroup Confirmatory Factor Analysis: Locating the Invariant Referent Sets
French, Brian F.; Finch, W. Holmes
2008-01-01
Multigroup confirmatory factor analysis (MCFA) is a popular method for the examination of measurement invariance and specifically, factor invariance. Recent research has begun to focus on using MCFA to detect invariance for test items. MCFA requires certain parameters (e.g., factor loadings) to be constrained for model identification, which are…
Application de la methode des sous-groupes au calcul Monte-Carlo multigroupe
Martin, Nicolas
effects of the scattering reaction consistent with the subgroup method. In this study, we generalize the Discrete Angle Technique, already proposed for homogeneous, multigroup cross sections, to isotopic cross sections on the form of probability tables. In this technique, the angular density is discretized into probability tables. Similarly to the cross-section case, a moment approach is used to compute the probability tables for the scattering cosine. (4) The introduction of a leakage model based on the B1 fundamental mode approximation. Unlike deterministic lattice packages, most Monte Carlo-based lattice physics codes do not include leakage models. However the generation of homogenized and condensed group constants (cross sections, diffusion coefficients) require the critical flux. This project has involved the development of a program into the DRAGON framework, written in Fortran 2003 and wrapped with a driver in C, the GANLIB 5. Choosing Fortran 2003 has permitted the use of some modern features, such as the definition of objects and methods, data encapsulation and polymorphism. The validation of the proposed code has been performed by comparison with other numerical methods: (1) The continuous-energy Monte Carlo method of the SERPENT code. (2) The Collision Probability (CP) method and the discrete ordinates (SN) method of the DRAGON lattice code. (3) The multigroup Monte Carlo code MORET, coupled with the DRAGON code. Benchmarks used in this work are representative of some industrial configurations encountered in reactor and criticality-safety calculations: (1)Pressurized Water Reactors (PWR) cells and assemblies. (2) Canada-Deuterium Uranium Reactors (CANDU-6) clusters. (3) Critical experiments from the ICSBEP handbook (International Criticality Safety Benchmark Evaluation Program).
We witnessed an initial hyped period and enthusiasm on carbon nano tubes in the 1990s later went through a significant expansion into nano tubes of other materials (metal di chalcogenides, boron nitride, etc.) as well as various nano wires and nano rods. While much of the hype might have gone, the research on one-dimensional (1D) nano materials has matured as one of the most active research areas within the nano science and nano technology community, flourishing with ample, exciting, and new research opportunities. Just like any other research frontier, researchers working in the 1D nano materials field are constantly striving to develop new fundamental science as well as potential applications. It remains a common belief that versatility and tunability of 1D nano materials would challenge many new rising tasks coming from our resource and energy demanding modern society. The traditional semiconductor industry has produced so many devices and systems from transistors, sensors, lasers, and LEDs to more sophisticated solar panels, which are now part of our daily lives. By down sizing the core components or parts to 1D form, one might wonder how fundamentally the dimensionality and morphology would impact the device performance, this is, as always, requiring us to fully understand the structure-property relationship in 1D nano materials. It may be equally crucial in connecting discovery-driven fundamental science to market-driven technology industry concerning potentially relevant findings derived from these novel materials. The importance of a platform that allows active researchers in this field to present their new development in a timely and efficient manner is therefore self-evident. Following the success of two early special issues devoted to 1D nano materials, this is the third one in a row organized by the same group of guest editors, attesting that such a platform has been well received by the readers
Multi-group calculations for fast reactors
The paper deals with various causes of error in calculations. The first part sets out the mathematical approximations (diffusion approximation, Sn method, etc.), the numerical resolution methods (effect of integration step), the models used, and the implications of these various factors in the determination of the principal characteristics of a fast neutron reactor. The second part studies the effect on reactivity of variations of element cross-sections, using various fuels, in a reactor of rather hard spectrum. (author)
SCALE 6 computes problem-dependent multigroup (MG) cross sections through a combination of the conventional Bondarenko shielding-factor method and a deterministic pointwise (PW) transport calculation of the fine-structure spectra in the resolved resonance and thermal energy ranges. The PW calculation is performed by the CENTRM code using a 1-D cylindrical Wigner-Seitz model with the white boundary condition instead of the real rectangular cell shape to represent a lattice unit cell. The pointwise fluxes computed by CENTRM are not exact because a 1-D model is used for the transport calculation, which introduces discrepancies in the MG self-shielded cross sections, resulting in some deviation in the eigenvalue. In order to solve this problem, the method of characteristics (MOC) has been applied to enable the CENTRM PW transport calculation for a 2-D square pin cell. The computation results show that the new BONAMI/CENTRM-MOC procedure produces very precise self-shielded cross sections compared to MCNP reaction rates.
Multigroup cross sections of resonant nuclei considering moderator mass differences
The multigroup constants library MGCL in the nuclear criticality safety evaluation code system JACS has been produced by the Bondarenko method to treat self-shielding effects. For estimating errors of this treatment, the multigroup cross sections of MGCL are compared with those obtained by precise treatment, i.e. with the weighted cross sections by ultra-fine spectra of neutron. The precise calculations are made for homogeneous mixtures of a resonant nucleus (235U, 238U, 239Pu, 240Pu, 242Pu or 56Fe) and a fictitious moderator nucleus with mass number 1, 12 or 200. The ultra-fine spectrum is calculated by the RABBLE code. Distinct differences are found in the self-shielding factors by comparisons between both treatments. Moreover, as the mass number increases, depressions of the self-shielding factor at the resonance peaks and its enhancements at the window of resonances are observed. (author)
Cyclotron radiation by a multi-group method
A multi-energy group technique is developed to study conditions under which cyclotron radiation emission can shift a Maxwellian electron distribution into a non-Maxwellian; and if the electron distribution is non-Maxwellian, to study the rate of cyclotron radiation emission as compared to that emitted by a Maxwellian having the same mean electron density and energy. The assumptions in this study are: the electrons should be in an isotropic medium and the magnetic field should be uniform. The multi-group technique is coupled into a multi-group Fokker-Planck computer code to study electron behavior under the influence of cyclotron radiation emission in a self-consistent fashion. Several non-Maxwellian distributions were simulated to compare their cyclotron emissions with the corresponding energy and number density equivalent Maxwellian distribtions
Multigroup-multiwaves Lisrel modeling in tourist satisfaction analysis
Cristina Bernini; Silvia Cagnone
2013-01-01
The paper analyzes the influence of tourist heterogeneity on the Tourist Local System Overall Satisfaction and its changes over time. We investigate two aspects: if different tourists segmented according to their trip motivation (seaside, conference and sport) show the same pattern of evaluation toward some relevant features of the TLS and if the evaluation scheme is dynamic. At this aim, a Multigroup-Multiwaves Lisrel model is estimated on a data set from the Tourist Satisfaction Survey, con...
Optimal calculational schemes for solving multigroup photon transport problem
A scheme of complex algorithm for solving multigroup equation of radiation transport is suggested. The algorithm is based on using the method of successive collisions, the method of forward scattering and the spherical harmonics method, and is realized in the FORAP program (FORTRAN, BESM-6 computer). As an example the results of calculating reactor photon transport in water are presented. The considered algorithm being modified may be used for solving neutron transport problems
Nuclear data and multigroup methods in fast reactor calculations
The work deals with fast reactor multigroup calculations, and the efficient treatment of basic nuclear data, which serves as raw material for the calculations. Its purpose is twofold: to build a computer code system that handles a large, detailed library of basic neutron cross section data, (such as ENDF/B-III) and yields a compact set of multigroup cross sections for reactor calculations; to use the code system for comparative analysis of different libraries, in order to discover basic uncertainties that still exist in the measurement of neutron cross sections, and to determine their influence upon uncertainties in nuclear calculations. A program named NANICK which was written in two versions is presented. The first handles the American basic data library, ENDF/B-III, while the second handles the German basic data library, KEDAK. The mathematical algorithm is identical in both versions, and only the file management is different. This program calculates infinitely diluted multigroup cross sections and scattering matrices. It is complemented by the program NASIF that calculates shielding factors from resonance parameters. Different versions of NASIF were written to handle ENDF/B-III or KEDAK. New methods for evaluating in reactor calculations the long term behavior of the neutron flux as well as its fine structure are described and an efficient calculation of the shielding factors from resonance parameters is offered. (B.G.)
A nodal expansion method for solving the multigroup SP3 equations in the reactor code DYN3D
The core model DYN3D which has been developed for three-dimensional analyses of steady states and transients in thermal reactors with quadratic or hexagonal fuel assemblies is based on nodal methods for the solution of the two-group neutron diffusion equation. Loading cores with higher content of MOX fuel, the increase of the fuel cycle length and new types of reactors are challenging for these standard methods. A nodal expansion method for solving the equations of the simplified P3 approximation (SP3) of the multigroup transport equation was developed to improve the accuracy of the DYN3D code. In this paper, the method used in DYN3D-SP3 is described. It is applied for the pin-wise calculation of a steady state of the OECD/NEA and U.S. NRC PWR MOX/UO2 Core Transient Benchmark. The eigenvalue keff, assembly powers and the pin powers are computed. The results calculated with different approaches including diffusion theory are compared with the reference solution obtained from a heterogeneous transport calculation with the code DeCART. Different approaches of the diffusion coefficient used in the SP3 equations are investigated. The SP3 results obtained with the transport cross section of multigroup diffusion theory show the smallest deviations from the reference solution. These deviations are in the same order as the results of the code DORT, whereas the DORT and DYN3D calculations were carried out with the same library of group constants for homogenized pin cells. (authors)
MUXS: a code to generate multigroup cross sections for sputtering calculations
This report documents MUXS, a computer code to generate multigroup cross sections for charged particle transport problems. Cross sections generated by MUXS can be used in many multigroup transport codes, with minor modifications to these codes, to calculate sputtering yields, reflection coefficients, penetration distances, etc
SERKON program for compiling a multigroup library to be used in BETTY calculation
A SERKON-type program was written to compile data sets generated by FEDGROUP-3 into a multigroup library for BETTY calculation. A multigroup library was generated from the ENDF/B-IV data file and tested against the TRX-1 and TRX-2 lattices with good results. (author)
Variational nodal solution algorithms for multigroup criticality problems
Variational nodal transport methods are generalized for the treatment of multigroup criticality problems. The generation of variational response matrices is streamlined and automated through the use of symbolic manipulation. A new red-black partitioned matrix algorithm for the solution of the within-group equations is formulated and shown to be at once both a regular matrix splitting and a synthetic acceleration method. The methods are implemented in X- Y geometry as a module of the Argonne National Laboratory code DIF3D. For few group problems highly accurate P3 eigenvalues are obtained with computing times comparable to those of an existing interface-current nodal transport method
Status of multigroup cross-section data for shielding applications
Multigroup cross-section libraries for shielding applications in formats for direct use in discrete ordinates or Monte Carlo codes have long been a part of the Data Library Collection (DLC) of the Radiation Shielding Information Center (RSIC). In recent years libraries in more flexible and comprehensive formats, which allow the user to derive his own problem-dependent sets, have been added to the collection. The current status of both types is described, as well as projections for adding data libraries based on ENDF/B-V
Multigroup Free-atom Doppler-broadening Approximation. Experiment
Gray, Mark Girard [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2015-11-06
The multigroup energy Doppler-broadening approximation agrees with continuous energy Dopplerbroadening generally to within ten percent for the total cross sections of ^{1}H,^{ 56}Fe, and ^{235}U at 250 lanl. Although this is probably not good enough for broadening from room temperature through the entire temperature range in production use, it is better than any interpolation scheme between temperatures proposed to date, and may be good enough for extrapolation from high temperatures. The method deserves further study since additional improvements are possible.
Nonparametric Multi-group Membership Model for Dynamic Networks
Kim, Myunghwan; Leskovec, Jure
2013-01-01
Relational data-like graphs, networks, and matrices-is often dynamic, where the relational structure evolves over time. A fundamental problem in the analysis of time-varying network data is to extract a summary of the common structure and the dynamics of the underlying relations between the entities. Here we build on the intuition that changes in the network structure are driven by the dynamics at the level of groups of nodes. We propose a nonparametric multi-group membership model for dynami...
SNAP - a three dimensional neutron diffusion code
This report describes a one- two- three-dimensional multi-group diffusion code, SNAP, which is primarily intended for neutron diffusion calculations but can also carry out gamma calculations if the diffusion approximation is accurate enough. It is suitable for fast and thermal reactor core calculations and for shield calculations. SNAP can solve the multi-group neutron diffusion equations using finite difference methods. The one-dimensional slab, cylindrical and spherical geometries and the two-dimensional case are all treated as simple special cases of three-dimensional geometries. Numerous reflective and periodic symmetry options are available and may be used to reduce the number of mesh points necessary to represent the system. Extrapolation lengths can be specified at internal and external boundaries. (Author)
Highlights: • A new AFEN code, MGANSP3, is developed for simplified P3 (SP3) calculations. • Surface averaged partial currents are used for coupling the nodes. • Coarse group rebalancing method is applied to increase the speed of calculations. • Four benchmark problems are used to examine the accuracy of the MGANSP3 code. - Abstract: In this study, a new analytic function expansion nodal (AFEN) method was developed to solve multi-group and three dimensional neutron simplified P3 equations (SP3) in reactor cores with rectangular fuel assemblies. In this method, the intranodal fluxes are expanded into a set of analytic basis functions for each group and moment. The nodes are coupled through the surface averaged partial currents at each nodal interface. Thus, six boundary conditions at each group and Legendre moments have been considered. Coarse group rebalancing (CGR) method was applied to increase the speed of code calculations. The code takes few-groups cross sections produced by a lattice code such as WIMS and calculates the effective multiplication factor, zeroth and second moments of the flux in multi-group energy, reactivity, and the relative power density at each fuel assembly. The numerical results for different benchmark problems demonstrate that solution of SP3 equations by our AFEN method improves both effective multiplication factor (keff) and power distribution compared to our AFEN diffusion method, especially in heterogeneous geometry and mixed-oxide (MOX) fuel problems
BETA-S, Multi-Group Beta-Ray Spectra
1 - Description of program or function: BETA-S calculates beta-decay source terms and energy spectra in multigroup format for time-dependent radionuclide inventories of actinides, fission products, and activation products. Multigroup spectra may be calculated in any arbitrary energy-group structure. The code also calculates the total beta energy release rate from the sum of the average beta-ray energies as determined from the spectral distributions. BETA-S also provides users with an option to determine principal beta-decaying radionuclides contributing to each energy group. The CCC-545/SCALE 4.3 (or SCALE4.2) code system must be installed on the computer before installing BETA-S, which requires the SCALE subroutine library and nuclide-inventory generation from the ORIGEN-S code. 2 - Methods:Well-established models for beta-energy distributions are used to explicitly represent allowed, and 1., 2. - and 3. -forbidden transition types. Forbidden non-unique transitions are assumed to have a spectral shape of allowed transitions. The multigroup energy spectra are calculated by numerically integrating the energy distribution functions using an adaptive Simpson's Rule algorithm. Nuclide inventories are obtained from a binary interface produced by the ORIGEN-S code. BETA-S calculates the spectra for all isotopes on the binary interface that have associated beta-decay transition data in the ENSDF-95 library, developed for the BETA-S code. This library was generated from ENSDF data and contains 715 materials, representing approximately 8500 individual beta transition branches. 3 - Restrictions on the complexity of the problem: The algorithms do not treat positron decay transitions or internal conversion electrons. The neglect of positron transitions in inconsequential for most applications involving aggregate fission products, since most of the decay modes are via electrons. The neglect of internal conversion electrons may impact on the accuracy of the spectrum in the low
Generation of subgroup parameters from JENDL-2 based multigroup data set for FBR core materials
Subgroup method gives a more accurate treatment to the resonance absorption in nuclear reactors, especially when it is heterogeneous, than the usual multigroup method. An algorithm has been developed based on a modified form of Roth's procedure, to calculate subgroup parameters, from the multigroup table of self-shielding factors given against a set of temperatures and dilution cross sections. A program SPART has been written with this algorithm, and it has been used to generate subgroup parameters for some important fast reactor core materials from the JENDL-2 based multigroup set, recently created and validated at IGCAR. In this report, the algorithm is discussed, and the subgroup parameters generated are presented. (author)
Reflector modelling of small high leakage cores making use of multi-group nodal equivalence theory
This research focuses on modelling reflectors in typical material testing reactors (MTRs). Equivalence theory is used to homogenise and collapse detailed transport solutions to generate equivalent nodal parameters and albedo boundary conditions for reflectors, for subsequent use in full core nodal diffusion codes. This approach to reflector modelling has been shown to be accurate for two-group large commercial light water reactor (LWR) analysis, but has not been investigated for MTRs. MTRs are smaller, with much larger leakage, environment sensitivity and multi-group spectrum dependencies than LWRs. This study aims to determine if this approach to reflector modelling is an accurate and plausible homogenisation technique for the modelling of small MTR cores. The successful implementation will result in simplified core models, better accuracy and improved efficiency of computer simulations. Codes used in this study include SCALE 6.1, OSCAR-4 and EQUIVA (the last two codes are developed and used at Necsa). The results show a five times reduction in calculational time for the proposed reduced reactor model compared to the traditional explicit model. The calculated equivalent parameters however show some sensitivity to the environment used to generate them. Differences in the results compared to the current explicit model, require more careful investigation including comparisons with a reference result, before its implementation can be recommended. (authors)
Adaptive matrix formation (AMF) method has been developed for the numerical solution of the transient multigroup neutron diffusion and delayed precursor equations in two- and three-dimensional geometry. The method is applied to a general class of two- and three- dimensional problems. The results of numerical experiments, as well as comparison with space-time experimental results indicate that the method is accurate and that the two- and three-dimensional calculations can be performed at 'reasonable' computer costs. Moreover, the AMF method offers the flexibility of using smaller time steps between flux shape calculations to achieve a specified accuracy and capability, without encountering numerical problems that occur in the other conventional methods. There is a large considerable saving in computer time and costs due to the partitioning of the matrix adopted in the presented AMF method. The two- and three-dimensional problems were analyzed with the present calculations model to illustrate the accuracy and stability of the method. Furthermore, the stability of the investigated method has been tested for sinusoidal, ramp, and step-change reactivity insertions. The results are in a good agreement with those of the other less approximate methods, including the problems in which the reflector zone is perturbed
Development and verification of a nodal approach for solving the multigroup SP{sub 3} equations
Beckert, C. [Forschungszentrum Dresden-Rossendorf, Institute of Safety Research, P.O. Box 51 01 19, D-01314 Dresden (Germany); Grundmann, U. [Forschungszentrum Dresden-Rossendorf, Institute of Safety Research, P.O. Box 51 01 19, D-01314 Dresden (Germany)], E-mail: U.Grundmann@fzd.de
2008-01-15
The core model DYN3D which has been developed for three-dimensional analyses of steady states and transients in thermal reactors with quadratic or hexagonal fuel assemblies is based on nodal methods for the solution of the two-group neutron diffusion equation. Loading cores with higher content of MOX fuel, the increase of the fuel cycle length, and the consideration of new reactor types are challenging for these standard methods. A nodal expansion method for solving the equations of the simplified P{sub 3} (SP{sub 3}) approximation of the multigroup transport equation was developed to improve the accuracy of the DYN3D code. The method described in the paper is verified with pinwise calculations of a steady state of the OECD/NEA and US NRC PWR MOX/UO{sub 2} Core Transient Benchmark. The used 16-group cross section library was generated for DORT calculations with homogenized pin cells. Two different approximations of the diffusion coefficient which occurs in the within-group form of the SP{sub 3} equations are investigated. Using the transport cross section for the calculation of the diffusion coefficient gives much better results than those obtained with the removal cross section. The improvement of the results in comparison to a pinwise diffusion calculation is shown. The results are compared with the DORT and the heterogeneous reference solution of the code DeCART. Concerning the SP{sub 3} calculation using the diffusion coefficient based on the transport cross section (DYN3D-SP3-TR) the deviations of the eigenvalue k{sub eff} and the assembly powers from the transport solutions of DORT and DeCART are in the same order as those between the two transport solutions themselves. The improvement of the DYN3D-SP3-TR results in comparison to the diffusion calculation is presented. As the DYN3D-SP3-TR and DORT calculations are performed with homogenized pin cells, the pin powers of the two calculations are closer to each other than to the pin powers of the DeCART solution
Development and verification of a nodal approach for solving the multigroup SP3 equations
The core model DYN3D which has been developed for three-dimensional analyses of steady states and transients in thermal reactors with quadratic or hexagonal fuel assemblies is based on nodal methods for the solution of the two-group neutron diffusion equation. Loading cores with higher content of MOX fuel, the increase of the fuel cycle length, and the consideration of new reactor types are challenging for these standard methods. A nodal expansion method for solving the equations of the simplified P3 (SP3) approximation of the multigroup transport equation was developed to improve the accuracy of the DYN3D code. The method described in the paper is verified with pinwise calculations of a steady state of the OECD/NEA and US NRC PWR MOX/UO2 Core Transient Benchmark. The used 16-group cross section library was generated for DORT calculations with homogenized pin cells. Two different approximations of the diffusion coefficient which occurs in the within-group form of the SP3 equations are investigated. Using the transport cross section for the calculation of the diffusion coefficient gives much better results than those obtained with the removal cross section. The improvement of the results in comparison to a pinwise diffusion calculation is shown. The results are compared with the DORT and the heterogeneous reference solution of the code DeCART. Concerning the SP3 calculation using the diffusion coefficient based on the transport cross section (DYN3D-SP3-TR) the deviations of the eigenvalue keff and the assembly powers from the transport solutions of DORT and DeCART are in the same order as those between the two transport solutions themselves. The improvement of the DYN3D-SP3-TR results in comparison to the diffusion calculation is presented. As the DYN3D-SP3-TR and DORT calculations are performed with homogenized pin cells, the pin powers of the two calculations are closer to each other than to the pin powers of the DeCART solution. To estimate the contribution of
Multigroup covariance matrices for fast-reactor studies
This report presents the multigroup covariance matrices based on the ENDF/B-V nuclear data evaluations. The materials and reactions have been chosen according to the specifications of ORNL-5517. Several cross section covariances, other than those specified by that report, are included due to the derived nature of the uncertainty files in ENDF/B-V. The materials represented are Ni, Cr, 16O, 12C, Fe, Na, 235U, 238U, 239Pu, 240Pu, 241Pu, and 10B (present due to its correlation to 238U). The data have been originally processed into a 52-group energy structure by PUFF-II and subsequently collapsed to smaller subgroup strutures. The results are illustrated in 52-group correlation matrix plots and tabulated into thirteen groups for convenience
Multigroup-multiwaves Lisrel modeling in tourist satisfaction analysis
Cristina Bernini
2013-05-01
Full Text Available The paper analyzes the influence of tourist heterogeneity on the Tourist Local System Overall Satisfaction and its changes over time. We investigate two aspects: if different tourists segmented according to their trip motivation (seaside, conference and sport show the same pattern of evaluation toward some relevant features of the TLS and if the evaluation scheme is dynamic. At this aim, a Multigroup-Multiwaves Lisrel model is estimated on a data set from the Tourist Satisfaction Survey, conducted in Rimini from 2004 to 2006 by the Faculty of Statistics – University of Bologna. The analysis shows that tourist evaluation scheme toward Rimini is quite similar among groups and over time, suggesting that differences among tourists do not affect TLS satisfaction.
Intragroup Socialization for Adult Korean Adoptees: A Multigroup Analysis
Kimberly J. Langrehr
2014-06-01
Full Text Available The purpose of the current study was to test a model of socialization among a sample of adult Korean adoptees. Based on the tenants of homophily and social identity theory, it was hypothesized that participants’ early racial and ethnic socialization experiences would account for their current intragroup friendships as adults, and that this relationship would be mediated by early intragroup contact and moderated by early ethnic identity status. The two ethnic and racial socialization variables (i.e., ethnic heritage activities and racial in-exposure significantly accounted for participants’ relationships with other Korean adoptees and nonadopted Koreans, and the effects were partially explained by early intragroup contact. Results of multigroup testing indicated the proposed socialization model was non-invariant across groups, such that the effects of ethnic heritage activities on intragroup contact and the effect of racial in-exposure on friendships with Korean adoptees were significantly different based on early ethnic identity status.
MORET: Version 4.B. A multigroup Monte Carlo criticality code
MORET 4 is a three dimensional multigroup Monte Carlo code which calculates the effective multiplication factor (keff) of any configurations more or less complex as well as reaction rates in the different volumes of the geometry and the leakage out of the system. MORET 4 is the Monte Carlo code of the APOLLO2-MORET 4 standard route of CRISTAL, the French criticality package. It is the most commonly used Monte Carlo code for French criticality calculations. During the last four years, the MORET 4 team has developed or improved the following major points: modernization of the geometry, implementation of perturbation algorithms, source distribution convergence, statistical detection of stationarity, unbiased variance estimation and creation of pre-processing and post-processing tools. The purpose of this paper is not only to present the new features of MORET but also to detail clearly the physical models and the mathematical methods used in the code. (author)
Multigroup representation of fusion product orbits in a plasma column
A method is derived for describing the time-depending behavior of α particles produced in a radially nonuniform slender plasma column as a distribution function among the possible orbits. A multigroup numerical approximation is introduced to analyze the development of the distribution function and its moments. Results are presented of calculations of the time-dependent α-particle energy spectrum and radial density, energy, and electron heating profiles in plasma columns with radii comparable to the α Larmor radius. This technique allows calculation of the α particle history at much more rapid rates than allowed by Monte Carlo technuques: The characteristic time scale is the α-electron slowing-down time rather than the cyclotron period
Starting from the radiation transport equation for homogeneous, refractive lossy media, we derive the corresponding time-dependent multifrequency diffusion equations. Zeroth and first moments of the transport equation couple the energy density, flux and pressure tensor. The system is closed by neglecting the temporal derivative of the flux and replacing the pressure tensor by its diagonal analogue. The radiation equations are coupled to a diffusion equation for the matter temperature. We are interested in modeling heating and cooling of silica (SiO2), at possibly rapid rates. Hence, in contrast to related work, we retain the temporal derivative of the radiation field. We derive boundary conditions at a planar air-silica interface taking account of reflectivities obtained from the Fresnel relations that include absorption. The spectral dimension is discretized into a finite number of intervals leading to a system of multigroup diffusion equations. Three simulations are presented. One models cooling of a silica slab, initially at 2500 K, for 10 s. The other two are 1D and 2D simulations of irradiating silica with a CO2 laser, λ = 10.59 μm. In 2D, a laser beam (Gaussian profile, r0 = 0.5 mm for 1/e decay) shines on a disk (radius = 0.4, thickness = 0.4 cm).
Kalpakkam multigroup cross section set for fast reactor applications - status and performance
This report documents the status of the presently created set of multigroup constants at Kalpakkam. The list of nuclides processed and the details of multigroup structure are given. Also included are the particulars of dilutions and temperatures for each nuclide in the multigroup cross section set for which self shielding factors have been calculated. Using this new multigroup cross section set, measured integral quantities such as K-eff, central reaction rate ratios, central reactivity worths etc. were calculated for a few fast critical benchmark assemblies and the calculated values of neutronic parameters obtained were compared with those obtained using the available Cadarache cross section library and those published in literature for ENDF/B-IV based set and Japanese evaluated nuclear data library (JENDL). The details of analyses are documented along with the conclusions. (author). 17 refs., 12 tabs
MCMG: a 3-D multigroup P3 Monte Carlo code and its benchmarks
In this paper a 3-D Monte Carlo multigroup neutron transport code MCMG has been developed from a coupled neutron and photon transport Monte Carlo code MCNP. The continuous-energy cross section library of the MCNP code is replaced by the multigroup cross section data generated by the transport lattice code, such as the WIMS code. It maintains the strong abilities of MCNP for geometry treatment, counting, variance reduction techniques and plotting. The multigroup neutron scattering cross sections adopt the Pn (n ≤ 3) approximation. The test results are in good agreement with the results of other methods and experiments. The number of energy groups can be varied from few groups to multigroup, and either macroscopic or microscopic cross section can be used. (author)
One-Dimensional (1-D) Nanoscale Heterostructures
Guozhen SHEN; Di CHEN; Yoshio BANDO; Dmitri GOLBERG
2008-01-01
One-dimensional (1-D) nanostructures have been attracted much attention as a result of their exceptional properties, which are different from bulk materials. Among 1-D nanostructures, 1-D heterostructures with modulated compositions and interfaces have recently become of particular interest with respect to potential applications in nanoscale building blocks of future optoelectronic devices and systems. Many kinds of methods have been developed for the synthesis of 1-D nanoscale heterostructures. This article reviews the most recent development, with an emphasize on our own recent efforts, on 1-D nanoscale heterostructures, especially those synthesized from the vapor deposition methods, in which all the reactive precursors are mixed together in the reaction chamber. Three types of 1-D nanoscale heterostructures, defined from their morphologies characteristics, are discussed in detail, which include 1-D co-axial core-shell heterostructures, 1-D segmented heterostructures and hierarchical heterostructures. This article begins with a brief survey of various methods that have been developed for synthesizing 1-D nanoscale heterostructures and then mainly focuses on the synthesis, structures and properties of the above three types of nanoscale heterostructures. Finally, this review concludes with personal views towards the topic of 1-D nanoscale heterostructures.
Coste-Delclaux, M
2006-03-15
This document describes the improvements carried out for modelling the self-shielding phenomenon in the multigroup transport code APOLLO2. They concern the space and energy treatment of the slowing-down equation, the setting up of quadrature formulas to calculate reaction rates, the setting-up of a method that treats directly a resonant mixture and the development of a sub-group method. We validate these improvements either in an elementary or in a global way. Now, we obtain, more accurate multigroup reaction rates and we are able to carry out a reference self-shielding calculation on a very fine multigroup mesh. To end, we draw a conclusion and give some prospects on the remaining work. (author)
A variational finite element-spherical harmonics method is presented for the solution of the even-parity multigroup equations with anisotropic scattering and sources. It is shown that by using a simple and natural formulation the numerical implementation of the method for any desired geometry is greatly eased and the anisotropy of scatter treated without any difficulty. Numerical examples demonstrate the ability of the resulting code to solve geometrically complex multigroup problems. (Author)
Radiation Transport for Explosive Outflows: A Multigroup Hybrid Monte Carlo Method
Wollaeger, Ryan T.; van Rossum, Daniel R.; Graziani, Carlo; Couch, Sean M.; Jordan, George C., IV; Lamb, Donald Q.; Moses, Gregory A.
2013-12-01
We explore Implicit Monte Carlo (IMC) and discrete diffusion Monte Carlo (DDMC) for radiation transport in high-velocity outflows with structured opacity. The IMC method is a stochastic computational technique for nonlinear radiation transport. IMC is partially implicit in time and may suffer in efficiency when tracking MC particles through optically thick materials. DDMC accelerates IMC in diffusive domains. Abdikamalov extended IMC and DDMC to multigroup, velocity-dependent transport with the intent of modeling neutrino dynamics in core-collapse supernovae. Densmore has also formulated a multifrequency extension to the originally gray DDMC method. We rigorously formulate IMC and DDMC over a high-velocity Lagrangian grid for possible application to photon transport in the post-explosion phase of Type Ia supernovae. This formulation includes an analysis that yields an additional factor in the standard IMC-to-DDMC spatial interface condition. To our knowledge the new boundary condition is distinct from others presented in prior DDMC literature. The method is suitable for a variety of opacity distributions and may be applied to semi-relativistic radiation transport in simple fluids and geometries. Additionally, we test the code, called SuperNu, using an analytic solution having static material, as well as with a manufactured solution for moving material with structured opacities. Finally, we demonstrate with a simple source and 10 group logarithmic wavelength grid that IMC-DDMC performs better than pure IMC in terms of accuracy and speed when there are large disparities between the magnitudes of opacities in adjacent groups. We also present and test our implementation of the new boundary condition.
Variational methods in steady state diffusion problems
Classical variational techniques are used to obtain accurate solutions to the multigroup multiregion one dimensional steady state neutron diffusion equation. Analytic solutions are constructed for benchmark verification. Functionals with cubic trial functions and conservational lagrangian constraints are exhibited and compared with nonconservational functionals with respect to neutron balance and to relative flux and current at interfaces. Excellent agreement of the conservational functionals using cubic trial functions is obtained in comparison with analytic solutions
MCNP - transport calculations in ducts using multigroup albedo coefficients
In this work, the use of multigroup albedo coefficients in Monte Carlo calculations of particle reflection and transmission by ducts is investigated. The procedure consists in modifying the MCNP code so that an albedo matrix computed previously by deterministic methods or Monte Carlo is introduced into the program to describe particle reflection by a surface. This way it becomes possible to avoid the need of considering particle transport in the duct wall explicitly, changing the problem to a problem of transport in the duct interior only and reducing significantly the difficulty of the real problem. The probability of particle reflection at the duct wall is given, for each group, as the sum of the albedo coefficients over the final groups. The calculation is started by sampling a source particle and simulating its reflection on the duct wall by sampling a group for the emerging particle. The particle weight is then reduced by the reflection probability. Next, a new direction and trajectory for the particle is selected. Numerical results obtained for the model are compared with results from a discrete ordinates code and results from Monte Carlo simulations that take particle transport in the wall into account. (author)
Macroscopic multigroup constants for accelerator driven system core calculation
The high-level wastes stored in facilities above ground or shallow repositories, in close connection with its nuclear power plant, can take almost 106 years before the radiotoxicity became of the order of the background. While the disposal issue is not urgent from a technical viewpoint, it is recognized that extended storage in the facilities is not acceptable since these ones cannot provide sufficient isolation in the long term and neither is it ethical to leave the waste problem to future generations. A technique to diminish this time is to transmute these long-lived elements into short-lived elements. The approach is to use an Accelerator Driven System (ADS), a sub-critical arrangement which uses a Spallation Neutron Source (SNS), after separation the minor actinides and the long-lived fission products (LLFP), to convert them to short-lived isotopes. As an advanced reactor fuel, still today, there is a few data around these type of core systems. In this paper we generate macroscopic multigroup constants for use in calculations of a typical ADS fuel, take into consideration, the ENDF/BVI data file. Four energy groups are chosen to collapse the data from ENDF/B-VI data file by PREPRO code. A typical MOX fuel cell is used to validate the methodology. The results are used to calculate one typical subcritical ADS core. (author)
Dorval, Eric
2016-01-01
Neutron transport calculations by Monte Carlo methods are finding increased application in nuclear reactor simulations. In particular, a versatile approach entails the use of a 2-step pro-cedure, with Monte Carlo as a few-group cross section data generator at lattice level, followed by deterministic multi-group diffusion calculations at core level. In this thesis, the Serpent 2 Monte Carlo reactor physics burnup calculation code is used in order to test a set of diffusion coefficient model...
Multi-group SP3 approximation for simulation of a three-dimensional PWR rod ejection accident
Highlights: • The multi-group SP3 method developed and implemented in PARCS for the MOX analysis. • The verifications were performed in 2D and 3D, 2G and MG, diffusion and transport, with and without feedback. • All results show consistency with the reference results obtained from the ANL PN transport code VARIANT for steady-state and transport calculations. • It was found that the SP3 angular approximation captures sufficient transport effects for both steady-state and transient, and provides essentially the same results as the VARIANT P5 method. • From the transient results of the full-core problem, it was noted that MG is more conservative than 2G, and P1 is more conservative than SP3. - Abstract: Previous researchers have shown that the simplified P3 (SP3) approximation is capable of providing sufficiently high accuracy for both static and transient simulations for reactor core analysis with considerably less computational expense than higher order transport methods such as the discrete ordinate or the full spherical harmonics methods. The objective of this paper is to provide a consistent comparison of two-group (2G) and multi-group (MG) diffusion and SP3 transport for rod ejection accident (REA) in a practical light water reactor (LWR) problem. The analysis is performed on two numerical benchmarks, a 3 × 3 assembly mini-core and a full pressurized water reactor (PWR) core. The calculations were performed using pin homogenized and assembly homogenized cross sections for a series of benchmarks of increasing difficulty, in two-dimensional (2D) and three-dimensional (3D), 2G and MG, diffusion and transport, as well as with and without feedback. All results show consistency with the reference results obtained from higher-order methods. It is demonstrated that the analyzed problems show small group-homogenization effects, but relatively significant transport effects which are satisfactorily addressed by the SP3 transport method. The sensitivity tests
Development of a multi-group SN transport calculation code with unstructured tetrahedral meshes
This paper reviews the computational methods used in the MUST (Multi-group Unstructured geometry SN Transport) code for solving the multi-group Sn transport equation in general geometries and describes the status of development of MUST. MUST solves the multi-group transport equation with unstructured tetrahedral meshes for modeling complicated geometrical problems. For tetrahedral mesh generation, input generation, and output visualization, we developed a management program where the mesh generation is based on Gmsh and TetGen that are open softwares. The geometrical modeling is done with the commercial CAD softwares such as CATIA. MUST uses the discontinuous finite element method (DFEM) and two-sub cell balance methods with linear discontinuous expansion (LDEM-SCB) to spatially discretize the transport equation. We applied MUST to three neutron and gamma coupled test problems for testing MUST. (author)
Current theories for approximating the effects of stochastic media on radiation transport assume very limited physics such as one dimension, constant grey opacities, and no material energy balance equation. When applied to more complex physical problems, the standard theory fails to match the results from direct numerical simulations. This work presents the first direct numerical simulations of multigroup radiation transport coupled to a material temperature equation in a 2D stochastic medium that are compared to closures proposed by various authors. After extending it from grey to multigroup physics, one closure that is not commonly used successfully models the results in dilute systems where one material comprises less than 5% of the total. This closure is more accurate for related grey transport problems than it is for the multigroup problem. When the specific heats are material- and temperature-dependent, it is much more difficult to fit the direct numerical solutions with an approximate closure.
Optimal control in multi-group coupled within-host and between-host models
Eric Numfor
2016-03-01
Full Text Available We formulate and then analyze a multi-group coupled within-host model of ODEs and between-host model of ODE and first-order PDEs, using the Human Immunodeficiency Virus (HIV for illustration. The basic reproduction number of the multi-group coupled epidemiological model is derived, steady states solutions are calculated and stability analysis of equilbria is investigated. An optimal control problem for our model with drug treatment on the multi-group within-host system is formulated and analyzed. Ekeland's principle is used in proving existence and uniqueness of an optimal control pair. Numerical simulations based on the semi-implicit finite difference schemes and the forward-backward sweep iterative method are obtained.
Development of a Multi-Group Neutron Cross Section Library Generation System for PWR
Kim, Kang Seog; Hong, Ser Gi; Song, Jae Seung; Lee, Kyung Hoon; Cho, Jin Young; Kim, Ha Yong; Koo, Bon Seung; Shim, Hyung Jin; Park, Sang Yoon
2008-10-15
This report describes a generation system of multi-group cross section library which is used in the KARMA lattice calculation code. In particular, the theoretical methodologies, program structures, and input preparations for the constituent programs of the system are described in detail. The library generation system consists of the following five programs : ANJOY, GREDIT, MERIT, SUBDATA, and LIBGEN. ANJOY generates automatically the NJOY input files and two batch files for automatic NJOY run for all the nuclides considered. The automatic NJOY run gives TAPE 23 (PENDF output file of BROADR module of NJOY) and TAPE24 (GENDF output file of GROUPR module of NJOY) files for each nuclide. GREDIT prepares a formatted multi-group cross section file in which the cross sections are tabulated versus temperature and background cross section after reading the TAPE24 file. MERIT generates the hydrogen equivalence factors and the resonance integral tables by solving the slowing down equation with ultra-fine group cross sections which are prepared with the TAPE 23 file. SUBDATA generates the subgroup data including subgroup levels and weights after reading the MERIT output file. Finally, LIBGEN generates the final multi-group library file by assembling the data prepared in the previous steps and by reading the other data such as fission product yield data and decay data.The multi-group cross section library includes general multi-group cross sections, resonance data, subgroup data, fission product yield data, kappa-values (energy release per fission), and all the data which are required in the depletion calculation. The addition or elimination of the cross sections for some nuclides can be easily done by changing the LIBGEN input file if the general multi-group cross section and the subgroup data files are prepared.
Development of a Multi-Group Neutron Cross Section Library Generation System for PWR
This report describes a generation system of multi-group cross section library which is used in the KARMA lattice calculation code. In particular, the theoretical methodologies, program structures, and input preparations for the constituent programs of the system are described in detail. The library generation system consists of the following five programs : ANJOY, GREDIT, MERIT, SUBDATA, and LIBGEN. ANJOY generates automatically the NJOY input files and two batch files for automatic NJOY run for all the nuclides considered. The automatic NJOY run gives TAPE 23 (PENDF output file of BROADR module of NJOY) and TAPE24 (GENDF output file of GROUPR module of NJOY) files for each nuclide. GREDIT prepares a formatted multi-group cross section file in which the cross sections are tabulated versus temperature and background cross section after reading the TAPE24 file. MERIT generates the hydrogen equivalence factors and the resonance integral tables by solving the slowing down equation with ultra-fine group cross sections which are prepared with the TAPE 23 file. SUBDATA generates the subgroup data including subgroup levels and weights after reading the MERIT output file. Finally, LIBGEN generates the final multi-group library file by assembling the data prepared in the previous steps and by reading the other data such as fission product yield data and decay data.The multi-group cross section library includes general multi-group cross sections, resonance data, subgroup data, fission product yield data, kappa-values (energy release per fission), and all the data which are required in the depletion calculation. The addition or elimination of the cross sections for some nuclides can be easily done by changing the LIBGEN input file if the general multi-group cross section and the subgroup data files are prepared
Application of diffusion theory to the transport of neutral particles in fusion plasmas
It is shown that the widely held view that diffusion theory can not provide good accuracy for the transport of neutral particles in fusion plasmas is misplaced. In fact, it is shown that multigroup diffusion theory gives quite good accuracy as compared to the transport theory. The reasons for this are elaborated and some of the physical and theoretical reasons which make the multigroup diffusion theory provide good accuracy are explained. Energy dependence must be taken into consideration to obtain a realistic neutral atom distribution in fusion plasmas. There are two reasons for this; presence of either is enough to necessitate an energy dependent treatment. First, the plasma temperature varies spatially, and second, the ratio of charge-exchange to total plasma-neutral interaction cross section (c) is not close to one. A computer code to solve the one-dimensional multigroup diffusion theory in general geometry (slab, cylindrical and spherical) has been written for use on Cray computers, and its results are compared with those from the one-dimensional transport code ANISN to support the above finding. A fast, compact and versatile two-dimensional finite element multigroup diffusion theory code, FINAT, in X-Y and R-Z cylindrical/toroidal geometries has been written for use on CRAY computers. This code has been compared with the two dimensional transport code DOT-4.3. The accuracy is very good, and FENAT runs much faster compared even to DOT-4.3 which is a finite difference code
Social exploration of 1D games
Valente, Andrea; Marchetti, Emanuela
2013-01-01
In this paper the apparently meaningless concept of a 1 dimensional computer game is explored, via netnography. A small number of games was designed and implemented, in close contact with online communities of players and developers, providing evidence that 1 dimension is enough to produce intere...... interesting gameplay, to allow for level design and even to leave room for artistic considerations on 1D rendering. General techniques to re-design classic 2D games into 1D are also emerging from this exploration....
Two-dimensional multigroup finite element calculation of fast reactor in diffusion approximation
When a linear element of triangular shape is used the actual finite element calculation is relatively simple. Extensive programs for mesh generation were written for easy inputting the configuration of reactors. A number of other programs were written for plotting neutron flux fields in individual groups, the power distribution, spatial plotting of fields, etc. The operation of selected programs, data preparation and operating instructions are described and examples given of data and results. All programs are written in GIER ALGOL. The used method and the developed programs have demonstrated that they are a useful instrument for the calculation of criticality and the distribution of neutron flux and power of both fast and thermal reactors. (J.B.)
Ceolin, Celina; Vilhena, Marco T.; Bodmann, Bardo E.J., E-mail: vilhena@pq.cnpq.b, E-mail: bardo.bodmann@ufrgs.b [Universidade Federal do Rio Grande do Sul (PROMEC/UFRGS), Porto Alegre, RS (Brazil). Programa de Pos-Graduacao em Engenharia Mecanica; Alvim, Antonio Carlos Marques, E-mail: alvim@nuclear.ufrj.b [Universidade Federal do Rio de Janeiro (PEN/COPPE/UFRJ), RJ (Brazil). Coordenacao dos Programas de Pos-Graduacao de Engenharia. Programa de Energia Nuclear
2011-07-01
The authors solved analytically the neutron kinetic equations in a homogeneous slab, assuming the multi group energy model and six delayed neutron precursor groups by the Generalized Integral Laplace Transform Technique (GILTT) for a multi-layered slab. To this end, averaged values for the nuclear parameters in the multi-layered slab are used and the solution is constructed following the idea of Adomian's decomposition method upon reducing the heterogeneous problem to a set of recursive problems with constant parameters in the multi-layered slab. More specifically, the corrections that render the initially homogeneous problem into a heterogeneous one are plugged into the equation as successive source terms. To the best of our knowledge this sort of solution is novel and not found in literature. We further present some numerical simulations. (author)
Adoption as a Social Marker: The Diffusion of Products in a Multigroup Environment
Smaldino, Paul E; Hillis, Vicken; Bednar, Jenna
2015-01-01
Social identities are among the key factors driving social behavior in complex societies. Recent attention to social identity in consumer behavior indicates that individuals avoid products that might signal membership in an outgroup. Yet the population-level consequences of adoption as identity signaling are largely unknown. Whereas previous work has focused on asymmetric attraction and repulsion among groups with different social identities, here we consider the spread of innovations in a structured population in which there are multiple groups who don't want to be mistaken for one another, using both analytical and agent-based modeling. This formal analysis, the first to take the spatial structure of a population into account, allows us to consider empirically important factors, including demographic skew and communication scale, that likely influence overall patterns of adoption. We find that as products become emergent social markers, aversion to outgroup-associated products can decrease overall patterns ...
Intracellular facilitated diffusion: searchers, crowders and blockers
Brackley, C A; Marenduzzo, D
2013-01-01
In bacteria, regulatory proteins search for a specific DNA binding target via "facilitated diffusion": a series of rounds of 3D diffusion in the cytoplasm, and 1D linear diffusion along the DNA contour. Using large scale Brownian dynamics simulations we find that each of these steps is affected differently by crowding proteins, which can either be bound to the DNA acting as a road block to the 1D diffusion, or freely diffusing in the cytoplasm. Macromolecular crowding can strongly affect mechanistic features such as the balance between 3D and 1D diffusion, but leads to surprising robustness of the total search time.
Jones, K.; Johnston, R.; Manley, D.J.; Owen, D.; Charlton, C.
2015-01-01
We develop and apply a multilevel modeling approach that is simultaneously capable of assessing multigroup and multiscale segregation in the presence of substantial stochastic variation that accompanies ethnicity rates based on small absolute counts. Bayesian MCMC estimation of a log-normal Poisson
DIAMANT2 - A multigroup neutron transport program for triangular and hexagonal geometry
DIAMANT2 evolved out of the DIAMANT-code. DIAMANT2 solves the multigroup neutron transport equation in planar geometry using the Ssub(N) method. Spatial discretization is accomplished by taking finite differences on a meshgrid composed of equilateral triangles. This report contains a detailed documentation of the program and the input description. (orig./HJ)
A finite element option for the MARC transport/ diffusion theory computer code
The MARC multigroup transport/diffusion theory computer code has been extended to include a finite element option. The facility is available in two-dimensional geometry and has a novel feature in allowing high order polynomial approximations to the flux using an automated computer procedure. (U.K.)
YORP torques with 1D thermal model
Breiter, Slawomir; Czekaj, Maria
2010-01-01
A numerical model of the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect for objects defined in terms of a triangular mesh is described. The algorithm requires that each surface triangle can be handled independently, which implies the use of a 1D thermal model. Insolation of each triangle is determined by an optimized ray-triangle intersection search. Surface temperature is modeled with a spectral approach; imposing a quasi-periodic solution we replace heat conduction equation by the Helmholtz equation. Nonlinear boundary conditions are handled by an iterative, FFT based solver. The results resolve the question of the YORP effect in rotation rate independence on conductivity within the nonlinear 1D thermal model regardless of the accuracy issues and homogeneity assumptions. A seasonal YORP effect in attitude is revealed for objects moving on elliptic orbits when a nonlinear thermal model is used.
Mezzacappa, A; Bruenn, S W; Blondin, J M; Guidry, M W; Strayer, M R; Umar, A S
1996-01-01
We investigate neutrino-driven convection in core collapse supernovae and its ramifications for the explosion mechanism. We begin with an ``optimistic'' 15 solar mass precollapse model, which is representative of the class of stars with compact iron cores. This model is evolved through core collapse and bounce in one dimension using multigroup (neutrino-energy--dependent) flux-limited diffusion (MGFLD) neutrino transport and Lagrangian hydrodynamics, providing realistic initial conditions for the postbounce convection and evolution. Our two-dimensional simulation begins at 106 ms after bounce at a time when there is a well-developed gain region, and proceeds for 400 ms. We couple two-dimensional (PPM) hydrodynamics to one-dimensional MGFLD neutrino transport. At 225 ms after bounce we see large-scale convection behind the shock, characterized by high-entropy, mushroom-like, expanding upflows and dense, low-entropy, finger-like downflows. The upflows reach the shock and distort it from sphericity. The radial c...
1D ferrimagnetism in homometallic chains
Coronado Miralles, Eugenio; Gómez García, Carlos José; Borrás Almenar, Juan José
1990-01-01
The magnetic properties of the cobalt zigzag chain Co(bpy)(NCS)2 (bpy=2,2′‐bipyridine) are discussed on the basis of an Ising‐chain model that takes into account alternating Landé factors. It is emphasized, for the first time, that a homometallic chain containing only one type of site can give rise to a 1D ferrimagneticlike behavior. ,
A multigroup library HENDL2.1/SS (Hybrid Evaluated Nuclear Data Library/Self-Shielding) based on ENDF/B-VII.0 evaluate data has been generated using Bondarenko and flux calculator method for the correction of self-shielding effect of neutronics analyses. To validate the reliability of the multigroup library HENDL2.1/SS, transport calculations for fusion-fission hybrid system FDS-I were performed in this paper. It was verified that the calculations with the HENDL2.1/SS gave almost the same results with MCNP calculations and were better than calculations with the HENDL2.0/MG which is another multigroup library without self-shielding correction. The test results also showed that neglecting resonance self-shielding caused underestimation of the Keff, neutron fluxes and waste transmutation ratios in the multigroup calculations of FDS-I.
Zou Jun, E-mail: jzou@ipp.ac.cn [Institute of Plasma Physics, Chinese Academy of Sciences, Hefei, Anhui 230031 (China); School of Nuclear Science and Technology, University of Science and Technology of China, Hefei, Anhui 230027 (China); He Zhaozhong; Zeng Qin; Qiu Yuefeng; Wang Minghuang [Institute of Plasma Physics, Chinese Academy of Sciences, Hefei, Anhui 230031 (China); School of Nuclear Science and Technology, University of Science and Technology of China, Hefei, Anhui 230027 (China)
2010-12-15
A multigroup library HENDL2.1/SS (Hybrid Evaluated Nuclear Data Library/Self-Shielding) based on ENDF/B-VII.0 evaluate data has been generated using Bondarenko and flux calculator method for the correction of self-shielding effect of neutronics analyses. To validate the reliability of the multigroup library HENDL2.1/SS, transport calculations for fusion-fission hybrid system FDS-I were performed in this paper. It was verified that the calculations with the HENDL2.1/SS gave almost the same results with MCNP calculations and were better than calculations with the HENDL2.0/MG which is another multigroup library without self-shielding correction. The test results also showed that neglecting resonance self-shielding caused underestimation of the K{sub eff}, neutron fluxes and waste transmutation ratios in the multigroup calculations of FDS-I.
A hybrid multigroup/continuous-energy Monte Carlo algorithm is developed for solving the Boltzmann-Fokker-Planck equation. This algorithm differs significantly from previous charged-particle Monte Carlo algorithms. Most importantly, it can be used to perform both forward and adjoint transport calculations, using the same basic multigroup cross-section data. The new algorithm is fully described, computationally tested, and compared with a standard condensed history algorithm for coupled electron-photon transport calculations
FAYEZ MOUSTAFA MOAWAD, RAGAB
2016-01-01
[EN] The neutron diffusion equation is an approximation of the neutron transport equation that describes the neutron population in a nuclear reactor core. In particular, we will consider here VVER-type reactors which use the neutron diffusion equation discretized on hexagonal meshes. Most of the simulation codes of a nuclear power reactor use the multigroup neutron diffusion equation to describe the neutron distribution inside the reactor core.To study the stationary state of a reactor, the r...
Correction of multigroup cross sections for resolved resonance interference in mixed absorbers
The effect that interference between resolved resonances has on averaging multigroup cross sections is examined for thermal reactor-type problems. A simple and efficient numerical scheme is presented to correct a preprocessed multigroup library for interference effects. The procedure is implemented in a design oriented lattice physics computer code and compared with rigorous numerical calculations. The approximate method for computing resonance interference correction factors is applied to obtaining fine-group cross sections for a homogeneous uranium-plutonium mixture and a uranium oxide lattice. It was found that some fine group cross sections are changed by more than 40% due to resonance interference. The change in resonance interference correction factors due to burnup of a PWR fuel pin is examined and found to be small. The effect of resolved resonance interference on collapsed broad-group cross sections for thermal reactor calculations is discussed
Processing ENDF/B-V uncertainty data into multigroup covariance matrices
The purpose of this work is to develop and demonstrate the capability of processing Evaluated Nuclear Data File, system B, version five (ENDF/B-V) uncertainty data into multigroup covariance matrices. These covariances may then be folded with sensitivity coefficients to obtain uncertainties in selected integral parameters such as K-effective and breeding ratio. The project consisted of separating the previous uncertainty processor (PUFF) from the basic nuclear data cross section processor (MINX), updating the uncertanty processor to theENDF/B-V format, programming the processor for new uncertainty data, and demonstrating the processor capabilities by producing a multigroup covariance library. These capabilities were verified in various ways including hand calculations and comparisons with other known results. A computer code named PUFF-II was written to perform the task described above
MC2-2: a code to calculate fast neutron spectra and multigroup cross sections
MC2-2 is a program to solve the neutron slowing down problem using basic neutron data derived from the ENDF/B data files. The spectrum calculated by MC2-2 is used to collapse the basic data to multigroup cross sections for use in standard reactor neutronics codes. Four different slowing down formulations are used by MC2-2: multigroup, continuous slowing down using the Goertzel-Greuling or Improved Goertzel-Greuling moderating parameters, and a hyper-fine-group integral transport calculation. Resolved and unresolved resonance cross sections are calculated accounting for self-shielding, broadening and overlap effects. This document provides a description of the MC2-2 program. The physics and mathematics of the neutron slowing down problem are derived and detailed information is provided to aid the MC2-2 user in preparing input for the program and implementation of the program on IBM 370 or CDC 7600 computers
A Multigroup Method for the Calculation of Neutron Fluence with a Source Term
Heinbockel, J. H.; Clowdsley, M. S.
1998-01-01
Current research on the Grant involves the development of a multigroup method for the calculation of low energy evaporation neutron fluences associated with the Boltzmann equation. This research will enable one to predict radiation exposure under a variety of circumstances. Knowledge of radiation exposure in a free-space environment is a necessity for space travel, high altitude space planes and satellite design. This is because certain radiation environments can cause damage to biological and electronic systems involving both short term and long term effects. By having apriori knowledge of the environment one can use prediction techniques to estimate radiation damage to such systems. Appropriate shielding can be designed to protect both humans and electronic systems that are exposed to a known radiation environment. This is the goal of the current research efforts involving the multi-group method and the Green's function approach.
Updated multi-group cross sections of minor actinides with improved resonance treatment
The study of minor actinide in transmutation reactors and other future applications makes resonance self-shielding treatment a significant issue for criticality and isotope depletion. Resonance treatment for minor actinides has been carried out by subgroup method with improved interference effect through interference correction. Subgroup data was generated using RMET21 and GENP codes along with multi-group cross section data by NJOY nuclear data processing system. Updated multi-group cross section data library for a neutron transport code nTRACER was compared with solutions from MCNPX. The resonance interaction of uranium with minor actinides has been included by modified interference treatment of interference correction in subgroup methodology. The comparison of cross sections and multiplication factor in pin and assembly problems showed significant improvement from systematic resonance treatment especially for 237Np and 243Am. (author)
A conservative multi-group approach to the Boltzmann equations for reactive gas mixtures
Bisi, M.; Rossani, A.; Spiga, G.
2015-11-01
Starting from a simple kinetic model for a quaternary mixture of gases undergoing a bimolecular chemical reaction, multi-group integro-differential equations are derived for the particle distribution functions of all species. The procedure takes advantage of a suitable probabilistic formulation, based on the underlying collision frequencies and transition probabilities, of the relevant reactive kinetic equations of Boltzmann type. Owing to an appropriate choice of a sufficiently large number of weight functions, it is shown that the proposed multi-group equations are able to fulfil exactly, at any order of approximation, the correct conservation laws that must be inherited from the original kinetic equations, where speed was a continuous variable. Future developments are also discussed.
PHISICS multi-group transport neutronic capabilities for RELAP5
Epiney, A.; Rabiti, C.; Alfonsi, A.; Wang, Y.; Cogliati, J.; Strydom, G. [Idaho National Laboratory (INL), 2525 N. Fremont Ave., Idaho Falls, ID 83402 (United States)
2012-07-01
PHISICS is a neutronic code system currently under development at INL. Its goal is to provide state of the art simulation capability to reactor designers. This paper reports on the effort of coupling this package to the thermal hydraulic system code RELAP5. This will enable full prismatic core and system modeling and the possibility to model coupled (thermal-hydraulics and neutronics) problems with more options for 3D neutron kinetics, compared to the existing diffusion theory neutron kinetics module in RELAP5 (NESTLE). The paper describes the capabilities of the coupling and illustrates them with a set of sample problems. (authors)
TRIMARAN: a three dimensional multigroup P1 Monte Carlo code for criticallity studies
TRIMARAN is developed for safety analysis of nuclar components containing fissionnable materials: shipping casks, storage and cooling pools, manufacture and reprocessing plants. It solves the transport equation by Monte Carlo method in general three dimensional geometry with multigroup P1 approximation. A special representation of cross sections and numbers has been developed in order to reduce considerably the computing cost and allow this three dimensional code to compete with standard numerical program used in parametric studies
Ford, W.E. III; Roussin, R.W.; Petrie, L.M.; Diggs, B.R.; Comolander, H.E.
1979-01-01
Contents of the IBM version of the APMX system distributed by the Radiation Shielding Information Center (APMX-II) are described. Sample problems which demonstrate the procedure for implementing AMPX-II modules to generate point cross sections; generate multigroup neutron, photon production, and photon interaction cross sections for various transport codes; collapse multigroup cross sections; check, edit, and punch multigroup cross sections; and execute a one-dimensional discrete ordinates transport calculation are detailed. 25 figures, 9 tables.
2D/1D approximations to the 3D neutron transport equation. I: Theory
A new class of '2D/1D' approximations is proposed for the 3D linear Boltzmann equation. These approximate equations preserve the exact transport physics in the radial directions x and y and diffusion physics in the axial direction z. Thus, the 2D/1D equations are more accurate approximations of the 3D Boltzmann equation than the conventional 3D diffusion equation. The 2D/1D equations can be systematically discretized, to yield accurate simulation methods for 3D reactor core problems. The resulting solutions will be more accurate than 3D diffusion solutions, and less expensive to generate than standard 3D transport solutions. In this paper, we (i) show that the simplest 2D/1D equation has certain desirable properties, (ii) systematically discretize this equation, and (iii) derive a stable iteration scheme for solving the discrete system of equations. In a companion paper [1], we give numerical results that confirm the theoretical predictions of accuracy and iterative stability. (authors)
For satisfaction of future global customer needs, dedicated efforts are being coordinated internationally and pursued continuously at AREVA NP. The currently ongoing CONVERGENCE project is committed to the development of the ARCADIAR next generation core simulation software package. ARCADIAR will be put to global use by all AREVA NP business regions, for the entire spectrum of core design processes, licensing computations and safety studies. As part of the currently ongoing trend towards more sophisticated neutronics methodologies, an SP3 nodal transport concept has been developed for ARTEMIS which is the steady-state and transient core simulation part of ARCADIAR. For enabling a high computational performance, the SPN calculations are accelerated by applying multi-level coarse mesh re-balancing. In the current implementation, SP3 is about 1.4 times as expensive computationally as SP1 (diffusion). The developed SP3 solution concept is foreseen as the future computational workhorse for many-group 3D pin-by-pin full core computations by ARCADIAR. With the entire numerical workload being highly parallelizable through domain decomposition techniques, associated CPU-time requirements that adhere to the efficiency needs in the nuclear industry can be expected to become feasible in the near future. The accuracy enhancement obtainable by using SP3 instead of SP1 has been verified by a detailed comparison of ARTEMIS 16-group pin-by-pin SPN results with KAERI's DeCart reference results for the 2D pin-by-pin Purdue UO2/MOX benchmark. This article presents the accuracy enhancement verification and quantifies the achieved ARTEMIS-SP3 computational performance for a number of 2D and 3D multi-group and multi-box (up to pin-by-pin) core computations. (authors)
Reference calculations on critical assemblies with Apollo2 code working with a fine multigroup mesh
The objective of this thesis is to add to the multigroup transport code APOLLO2 the capability to perform deterministic reference calculations, for any type of reactor, using a very fine energy mesh of several thousand groups. This new reference tool allows us to validate the self-shielding model used in industrial applications, to perform depletion calculations, differential effects calculations, critical buckling calculations or to evaluate precisely data required by the self shielding model. At its origin, APOLLO2 was designed to perform routine calculations with energy meshes around one hundred groups. That is why, in the current format of cross sections libraries, almost each value of the multigroup energy transfer matrix is stored. As this format is not convenient for a high number of groups (concerning memory size), we had to search out a new format for removal matrices and consequently to modify the code. In the new format we found, only some values of removal matrices are kept (these values depend on a reconstruction precision choice), the other ones being reconstructed by a linear interpolation, what reduces the size of these matrices. Then we had to show that APOLLO2 working with a fine multigroup mesh had the capability to perform reference calculations on any assembly geometry. For that, we successfully carried out the validation with several calculations for which we compared APOLLO2 results (obtained with the universal mesh of 11276 groups) to results obtained with Monte Carlo codes (MCNP, TRIPOLI4). Physical analysis led with this new tool have been very fruitful and show a great potential for such an R and D tool. (author)
Hybrid method of deterministic and probabilistic approaches for multigroup neutron transport problem
A hybrid method of deterministic and probabilistic methods is proposed to solve Boltzmann transport equation. The new method uses a deterministic method, Method of Characteristics (MOC), for the fast and thermal neutron energy ranges and a probabilistic method, Monte Carlo (MC), for the intermediate resonance energy range. The hybrid method, in case of continuous energy problem, will be able to take advantage of fast MOC calculation and accurate resonance self shielding treatment of MC method. As a proof of principle, this paper presents the hybrid methodology applied to a multigroup form of Boltzmann transport equation and confirms that the hybrid method can produce consistent results with MC and MOC methods. (authors)
REX1-87, Multigroup Neutron Cross-Sections from ENDF/B
1 - Description of program or function: The program calculates self- shielding factors for reactor applications from a pre-processed (linearized) evaluated nuclear data file in the ENDF/B format. 2 - Method of solution: Bondarenko definition of multigroup self- shielding factors invoking narrow resonance treatment is used. 3 - Restrictions on the complexity of the problem: a) Maximum no. of energy group is 620. b) Only the built-in forms of the weighting functions can be chosen. c) The program is strictly limited to resolved resonance region from physical considerations
Multi-group pin power reconstruction method based on colorset form functions
A multi-group pin power reconstruction method that fully exploits nodal information obtained from global coarse mesh solution has been developed. It expands the intra-nodal flux distributions into nonseparable semi-analytic basis functions, and a colorset based form function generating method is proposed, which can accurately model the spectral interaction occurring at assembly interface. To demonstrate its accuracy and applicability to realistic problems, the new method is tested against two benchmark problems, including a mixed-oxide fuel problem. The results show that the new methods is comparable in accuracy to fine-mesh methods. (authors)
Specifications for a two-dimensional multi-group scattering code: ALCI
This report describes the specifications of the ALCI programme. This programme resolves the system of difference equations similar to the homogeneous problem of multigroup neutron scattering, with two dimensions in space, in the three geometries XY, RZ, RΘ. It is possible with this method to calculate geometric and composition criticalities and also to calculate the accessory problem on demand. The maximum number of points dealt with is 6000. The maximum permissible number of groups is 12. The internal iterations are treated by the method of alternating directions. The external iterations are accelerated using the extrapolation method due to Tchebychev. (authors)
Global analysis on a class of multi-group SEIR model with latency and relapse.
Wang, Jinliang; Shu, Hongying
2016-02-01
In this paper, we investigate the global dynamics of a multi-group SEIR epidemic model, allowing heterogeneity of the host population, delay in latency and delay due to relapse distribution for the human population. Our results indicate that when certain restrictions on nonlinear growth rate and incidence are fulfilled, the basic reproduction number R0 plays the key role of a global threshold parameter in the sense that the long-time behaviors of the model depend only on R0. The proofs of the main results utilize the persistence theory in dynamical systems, Lyapunov functionals guided by graph-theoretical approach. PMID:26776266
Modification of the resonance treatment in multigroup neutron slowing-down codes
The previously reported computer codes GRACE and BETTY for resonance treatment in the multigroup neutron slowing-down processes have been improved, employing the new results of resonance absorption calculations. The total resonance integral formulae were changed, 239Pu resonance integral data were included in the library of group constants and the selection of partial resonance integral distribution functions was automatized. The users of the GRACE and BETTY codes are provided with a more credible and more comfortable resonance treatment. Explicit description of modification of user's manuals is given. (D.P.)
On the completeness of the multigroup eigenfunctions set of a reactor system Boltzmann operator
An example is given, which illustrates how the set of the eigenfunctions shifts from incompleteness to completeness when a coupling relationship is established between the spectrum of the neutrons produced by fission and the energy of the neutrons which generate the fissions. The proposed method allows one to complete the set of eigenfunctions of the Boltzmann operator in the multigroup case. That, in principle, enlarges the possibility to apply the SM, Standard Method, and the GSM, Generalized Standard Method, to any problem in reactor physics, regardless of the number of energy groups. (author)
Geospatial Data Fusion and Multigroup Decision Support for Surface Water Quality Management
Sun, A. Y.; Osidele, O.; Green, R. T.; Xie, H.
2010-12-01
Social networking and social media have gained significant popularity and brought fundamental changes to many facets of our everyday life. With the ever-increasing adoption of GPS-enabled gadgets and technology, location-based content is likely to play a central role in social networking sites. While location-based content is not new to the geoscience community, where geographic information systems (GIS) are extensively used, the delivery of useful geospatial data to targeted user groups for decision support is new. Decision makers and modelers ought to make more effective use of the new web-based tools to expand the scope of environmental awareness education, public outreach, and stakeholder interaction. Environmental decision processes are often rife with uncertainty and controversy, requiring integration of multiple sources of information and compromises between diverse interests. Fusing of multisource, multiscale environmental data for multigroup decision support is a challenging task. Toward this goal, a multigroup decision support platform should strive to achieve transparency, impartiality, and timely synthesis of information. The latter criterion often constitutes a major technical bottleneck to traditional GIS-based media, featuring large file or image sizes and requiring special processing before web deployment. Many tools and design patterns have appeared in recent years to ease the situation somewhat. In this project, we explore the use of Web 2.0 technologies for “pushing” location-based content to multigroups involved in surface water quality management and decision making. In particular, our granular bottom-up approach facilitates effective delivery of information to most relevant user groups. Our location-based content includes in-situ and remotely sensed data disseminated by NASA and other national and local agencies. Our project is demonstrated for managing the total maximum daily load (TMDL) program in the Arroyo Colorado coastal river basin
Xiaoming Fan
2014-01-01
Full Text Available We discuss multigroup SIRS (susceptible, infectious, and recovered epidemic models with random perturbations. We carry out a detailed analysis on the asymptotic behavior of the stochastic model; when reproduction number ℛ0>1, we deduce the globally asymptotic stability of the endemic equilibrium by measuring the difference between the solution and the endemic equilibrium of the deterministic model in time average. Numerical methods are employed to illustrate the dynamic behavior of the model and simulate the system of equations developed. The effect of the rate of immunity loss on susceptible and recovered individuals is also analyzed in the deterministic model.
On the dynamics of a class of multi-group models for vector-borne diseases
Iggidr, Abderrahman; Sallet, Gauthier; Souza, Max O.
2016-01-01
The resurgence of vector-borne diseases is an increasing public health concern, and there is a need for a better understanding of their dynamics. For a number of diseases, e.g. dengue and chikungunya, this resurgence occurs mostly in urban environments, which are naturally very heterogeneous, particularly due to population circulation. In this scenario, there is an increasing interest in both multi-patch and multi-group models for such diseases. In this work, we study the dynamics of a vector...
Hursin, Mathieu [School of Nuclear Engineering, Purdue University, 400 Central Drive, IN 47907 (United States); Xiao Shanjie [School of Nuclear Engineering, Purdue University, 400 Central Drive, IN 47907 (United States); Jevremovic, Tatjana [School of Nuclear Engineering, Purdue University, 400 Central Drive, IN 47907 (United States)]. E-mail: tatjanaj@purdue.edu
2006-09-15
This paper summarizes the theoretical and numerical aspects of the AGENT code methodology accurately applied for detailed three-dimensional (3D) multigroup steady-state modeling of neutron interactions in complex heterogeneous reactor domains. For the first time we show the fine-mesh neutron scalar flux distribution in Purdue research reactor (that was built over forty years ago). The AGENT methodology is based on the unique combination of the three theories: the method of characteristics (MOC) used to simulate the neutron transport in two-dimensional (2D) whole core heterogeneous calculation, the theory of R-functions used as a mathematical tool to describe the true geometry and fuse with the MOC equations, and one-dimensional (1D) higher-order diffusion correction of 2D transport model to account for full 3D heterogeneous whole core representation. The synergism between the radial 2D transport and the 1D axial transport (to take into account the axial neutron interactions and leakage), called the 2D/1D method (used in DeCART and CHAPLET codes), provides a 3D computational solution. The unique synergism between the AGENT geometrical algorithm capable of modeling any current or future reactor core geometry and 3D neutron transport methodology is described in details. The 3D AGENT accuracy and its efficiency are demonstrated showing the eigenvalues, point-wise flux and reaction rate distributions in representative reactor geometries. The AGENT code, comprising this synergism, represents a building block of the computational system, called the virtual reactor. Its main purpose is to perform 'virtual' experiments and demonstrations of various mainly university research reactor experiments.
This paper summarizes the theoretical and numerical aspects of the AGENT code methodology accurately applied for detailed three-dimensional (3D) multigroup steady-state modeling of neutron interactions in complex heterogeneous reactor domains. For the first time we show the fine-mesh neutron scalar flux distribution in Purdue research reactor (that was built over forty years ago). The AGENT methodology is based on the unique combination of the three theories: the method of characteristics (MOC) used to simulate the neutron transport in two-dimensional (2D) whole core heterogeneous calculation, the theory of R-functions used as a mathematical tool to describe the true geometry and fuse with the MOC equations, and one-dimensional (1D) higher-order diffusion correction of 2D transport model to account for full 3D heterogeneous whole core representation. The synergism between the radial 2D transport and the 1D axial transport (to take into account the axial neutron interactions and leakage), called the 2D/1D method (used in DeCART and CHAPLET codes), provides a 3D computational solution. The unique synergism between the AGENT geometrical algorithm capable of modeling any current or future reactor core geometry and 3D neutron transport methodology is described in details. The 3D AGENT accuracy and its efficiency are demonstrated showing the eigenvalues, point-wise flux and reaction rate distributions in representative reactor geometries. The AGENT code, comprising this synergism, represents a building block of the computational system, called the virtual reactor. Its main purpose is to perform 'virtual' experiments and demonstrations of various mainly university research reactor experiments
Analysis list: Nr1d2 [Chip-atlas[Archive
Full Text Available Nr1d2 Liver + mm9 http://dbarchive.biosciencedbc.jp/kyushu-u/mm9/target/Nr1d2.1.tsv... http://dbarchive.biosciencedbc.jp/kyushu-u/mm9/target/Nr1d2.5.tsv http://dbarchive.biosciencedbc.jp/kyushu-u/mm9/target/Nr1d2....10.tsv http://dbarchive.biosciencedbc.jp/kyushu-u/mm9/colo/Nr1d2.Liver.tsv http://dbarchive.biosciencedbc.jp/kyushu-u/mm9/colo/Liver.gml ...
A new derivation of Akcasu's 'MLP' equations for 1-D particle transport in stochastic media
This paper presents a new derivation of Akcasu's modified Levermore-Pomraning (MLP) model, which estimates the ensemble-averaged angular flux for particle transport problems in 1-D geometrically random media. The significant new feature of the MLP equations is that, unlike the earlier Levermore-Pomraning (LP) model, the MLP equations are exact for certain classes of problems with scattering. We also show, via asymptotic analyses, that the MLP equations are accurate in the atomic mix and diffusion limits
The discrete ordinates method is the most powerful and generally used deterministic method to obtain approximate solutions of the Boltzmann transport equation. A finite element formulation, utilizing a canonical form of the transport equation, is here developed to obtain both integral and pointwise solutions to neutron transport problems. The formulation is based on the use of linear triangles. A general treatment of anisotropic scattering is included by employing discrete ordinates-like approximations. In addition, multigroup source outer iteration techniques are employed to perform group-dependent calculations. The ability of the formulation to reduce substantially ray effects and its ability to perform streaming calculations are demonstrated by analyzing a series of test problems. The anisotropic scattering and multigroup treatments used in the development of the formulation are verified by a number of one-dimensional comparisons. These comparisons also demonstrate the relative accuracy of the formulation in predicting integral parameters. The applicability of the formulation to nonorthogonal planar geometries is demonstrated by analyzing a hexagonal-type lattice. A small, high-leakage reactor model is analyzed to investigate the effects of varying both the spatial mesh and order of angular quadrature. This analysis reveals that these effects are more pronounced in the present formulation than in other conventional formulations. However, the insignificance of these effects is demonstrated by analyzing a realistic reactor configuration. In addition, this final analysis illustrates the importance of incorporating anisotropic scattering into the finite element formulation. 8 tables, 29 figures
Comparative calculations of the experimental benchmark of iron sphere with Cf source have been performed in order to assess the sensibility of the calculations of neutron transmission through iron media to different multigroup libraries generated on the base of ENDF/B-6 and ENDF/B-4. Similar calculations and comparison of the neutron flux passed through media typical as geometry and material compositions for the WWER-1000 and WWER-440 vessels have been carried out. Except the already well-known problem dependent libraries, the new libraries BGL-440 and BGL-1000 generated on the base of ENDF/B-6 for the WWER-440 and WWER-1000 RPV neutron fluence calculations have been applied. The solving of neutron transport through iron media using ENDF/B-6 data gives better consistency with the experiment than using ENDF/B-4. The latter underestimate the experimental fluxes more substantially in the energy range above 2 MeV and the evaluations of the neutron flux responses for the WWER vessel surveillance is preferably to be carried out by the appropriate BGL library. Key words: neutron transport, multigroup neutron cross section libraries
Verification of a Multi-group Cross Section Library for Burnup Calculation
Daing, Aung Tharn; Kim, Myung Hyun [Kyung Hee Univ., Yongin (Korea, Republic of); Joo, Hang Yu [Seoul National Univ., Seoul (Korea, Republic of)
2013-05-15
Despite satisfying the estimation of the neutronic parameters without depletion to some extent, it still requires detailed investigation of the behavior of a fuel with strong neutron absorber over its operating life time by nTRACER, the direct whole core calculation code with the conventional semi Predictor-Corrector method. This study is mainly focused on the verification of the newly generated multi-group library for burnup calculation by nTRACER through the analysis of its performance of depletion calculation of UO{sub 2} fuel with strong neutron absorbers such as Gadolinium. Firstly, the depletion calculation results of nTRACER are presented by comparing the evolution of k-inf and the inventories of commonly found important isotopes as a function of burnup in the cases of gadolinia(GAD)-bearing fuel pin and fuel assembly (FA) with those of MCNPX-version.2.6.0. The newly generated multi-group library for burnup calculation by nTRACER was verified through GAD-bearing fuel after the new approach of resonance treatment had been employed. Though very good agreement in the overall effect reflected on the multiplication factor of FA at BOC, the evolution of k-inf along fuel irradiation history was systematically well underestimated by nTRACER when compared to Monte Carlo results.
Multi-group unified nodal method with two-group coarse-mesh finite difference formulation
The one-node kernels of the unified nodal method (UNM) which were originally developed for two-group (2G) problems are extended to solve multi-group (MG) problems within the framework of the 2G coarse-mesh finite difference (CMFD) formulation. The analytic nodal method (ANM) kernel of UNM is reformulated for the MG application by adopting the Pade approximation to avoid the similarity transform required to diagonalize the G x G buckling matrix. In addition, a one-node semi-analytic nodal method (SANM) kernel which is considered adequate for multi-group calculations is also integrated into the UNM formulation by expressing it in the form consistent with the other UNM kernels. As an efficient global solution framework, the 2G CMFD formulation with dynamic group condensation and prolongation is established and the performance of the various MG kernels is examined using various static and transient benchmark problems. It turns out that the SANM kernel is the best one for MG problems not only because it retains accuracy comparable to MGANM with a shorter computing time but also because its accuracy or its convergence does not depend on the eigenvalue range of the buckling matrix of the system. The 2G CMFD formulation with MG one-node UNM kernels turns out to be very effective in that it conveniently accelerates the MG source iteration
The group-level consequences of sexual conflict in multigroup populations.
Omar Tonsi Eldakar
Full Text Available In typical sexual conflict scenarios, males best equipped to exploit females are favored locally over more prudent males, despite reducing female fitness. However, local advantage is not the only relevant form of selection. In multigroup populations, groups with less sexual conflict will contribute more offspring to the next generation than higher conflict groups, countering the local advantage of harmful males. Here, we varied male aggression within- and between-groups in a laboratory population of water striders and measured resulting differences in local population growth over a period of three weeks. The overall pool fitness (i.e., adults produced of less aggressive pools exceeded that of high aggression pools by a factor of three, with the high aggression pools essentially experiencing no population growth over the course of the study. When comparing the fitness of individuals across groups, aggression appeared to be under stabilizing selection in the multigroup population. The use of contextual analysis revealed that overall stabilizing selection was a product of selection favoring aggression within groups, but selected against it at the group-level. Therefore, this report provides further evidence to show that what evolves in the total population is not merely an extension of within-group dynamics.
Historical review of group diffusion computation
Diffusion theory neutron flux computations are the backbone of reactor physics design studies. Early one-space-dimension computations of simple configurations with Marchand and Friden desk computers were swept with the computer revolution into detailed three-space-dimension calculations of complex reactor configurations on today's parallel and vector machines. Methods for modeling reactors have evolved during the past 45 yr. The multigroup diffusion theory first described by Ehrlich and Hurwitz has been embellished with transport theory corrections and accounting for mesh effects in discrete modeling. Spectral variations within a few-group framework were accounted for by space-dependent group constants based on estimated hardness of spectrum. A more sophisticated procedure used overlapping energy groups
P1 adaptation of TRIPOLI-4 code for the use of 3D realistic core multigroup cross section generation
In this paper, we discuss some improvements we recently implemented in the Monte-Carlo code TRIPOLI-4 associated with the homogenization and collapsing of subassemblies cross sections. The improvement offered us another approach to get critical multigroup cross sections with Monte-Carlo method. The new calculation method in TRIPOLI-4 tries to ensure the neutronic balances, the multiplicative factors and the critical flux spectra for some realistic geometries. We make it by at first improving the treatment of the energy transfer probability, the neutron excess weight and the neutron fission spectrum. This step is necessary for infinite geometries. The second step which will be enlarged in this paper is aimed at better dealing with the multigroup anisotropy distribution law for finite geometries. Usually, Monte-Carlo homogenized multi-group cross sections are validated within a core calculation by a deterministic code. Here, the validation of multigroup constants will also be carried out by Monte-Carlo core calculation code. Different subassemblies are tested with the new collapsing method, especially for the fast neutron reactors subassemblies. (authors)
Liquid Salt Cooled Reactors (LSCRs) are high temperature reactors, cooled by liquid salt, with a TRISO-particle based fuel in a solid form (stationary fuel elements or circulating fuel pebbles); this paper is focusing on the former. In either case, due to the double heterogeneity, core physics analysis require different considerations with more complex approaches than LWRs core physics calculations. Additional challenges appear when using the multi-group approach. In this paper we examine the use of SCALE6.1.1. Double heterogeneity may be accounted for through the Dancoff factor, however, SCALE6.1.1 does not provide an automated method to calculate Dancoff Factors for fuel planks with TRISO fuel particles. Therefore, depletion with continuous energy Monte Carlo Transport (CE depletion) in SCALE6.2 beta was used to generate MC Dancoff factors for multi-group calculations. MCDancoff corrected multi-group depletion agrees with the results for CE depletion within ±100 pcm, and within ±2σ. Producing MCDancoff factors for multi-group (MG) depletion calculations is necessary to LSCR analysis because CE depletion runtime and memory requirements are prohibitive for routine use. MG depletion with MCDancoff provides significantly shorter runtime and lower memory requirements while providing results of acceptable accuracy. (author)
Cai, Li; Pénéliau, Yannick; Diop, Cheikh M.; Malvagi, Fausto
2014-06-01
In this paper, we discuss some improvements we recently implemented in the Monte-Carlo code TRIPOLI-4® associated with the homogenization and collapsing of subassemblies cross sections. The improvement offered us another approach to get critical multigroup cross sections with Monte-Carlo method. The new calculation method in TRIPOLI-4® tries to ensure the neutronic balances, the multiplicative factors and the critical flux spectra for some realistic geometries. We make it by at first improving the treatment of the energy transfer probability, the neutron excess weight and the neutron fission spectrum. This step is necessary for infinite geometries. The second step which will be enlarged in this paper is aimed at better dealing with the multigroup anisotropy distribution law for finite geometries. Usually, Monte-Carlo homogenized multi-group cross sections are validated within a core calculation by a deterministic code. Here, the validation of multigroup constants will also be carried out by Monte-Carlo core calculation code. Different subassemblies are tested with the new collapsing method, especially for the fast neutron reactors subassemblies.
The FOEHN critical experiments were analyzed to validate the use of multigroup cross sections in the design of the Advanced Neutron Source. Eleven critical configurations were evaluated using the KENO, DORT, and VENTURE neutronics codes. Eigenvalue and power density profiles were computed and show very good agreement with measured values
A - Nature of physical problem solved: The ANISN system treats neutron and gamma transport in one- dimensional plane, spherical and cylinder geometry. The multigroup cross sections prepared by the programs LIANE and SUPERTOG are processed by the program RETTOG, which produces a binary library with Legendre expansions. The binary library can be updated and edited with the program LGR/B. The photon multigroup cross sections are created with the program GAMLEG/A. If the bulk of the data is too large, the program TAPEMA produces a special group-by-group library. The volume sources are calculated from a reduced set of input data and punched in a format suitable for input to ANISN, using the program PRESOU. ANISN calculates fluxes by groups, space intervals, angle and any number of reaction rates. The energy and space dependent fluxes are stored on tape and can be reprocessed, edited and plotted with the program ANISEX, which also permits to calculate supplementary reaction rates. The program ANISN can condense cross sections into a reduced number of groups. The ANISN system is used as a reference system for the evaluation of approximation methods (space-diffusion or point kernel) or for the preparation of multigroup libraries for two-dimensional transport codes (DOT). In particular it is used for shielding problems with high attenuation in water reactors and fast reactors. ANISN-E solves the same problems as the original ANISN code. Some modifications concern weighted cross sections output and fixed distributed sources input/output. ANISN-E (CCC-0082/09): The CYBER 175 version of ANISN-E also contains the free-format input capability. ANISN-JR extends the applicability of the original ANISN code for shielding analyses by adding options of calculating the reaction rates distributions from detector response, generating the volume- flux weighted cross sections in arbitrary regions or zones and plotting the neutron or gamma-ray spectra and the reaction rates distributions
Development of advanced nodal diffusion methods for modern computer architectures
A family of highly efficient multidimensional multigroup advanced neutron-diffusion nodal methods, ILLICO, were implemented on sequential, vector, and vector-concurrent computers. Three-dimensional realistic benchmark problems can be solved in vectorized mode in less than 0.73 s (33.86 Mflops) on a Cray X-MP/48. Vector-concurrent implementations yield speedups as high as 9.19 on an Alliant FX/8. These results show that the ILLICO method preserves essentially all of its speed advantage over finite-difference methods. A self-consistent higher-order nodal diffusion method was developed and implemented. Nodal methods for global nuclear reactor multigroup diffusion calculations which account explicitly for heterogeneities in the assembly nuclear properties were developed and evaluated. A systematic analysis of the zero-order variable cross section nodal method was conducted. Analyzing the KWU PWR depletion benchmark problem, it is shown that when burnup heterogeneities arise, ordinary nodal methods, which do not explicitly treat the heterogeneities, suffer a significant systematic error that accumulates. A nodal method that treats explicitly the space dependence of diffusion coefficients was developed and implemented. A consistent burnup-correction method for nodal microscopic depletion analysis was developed
MPI version of NJOY and its application to multigroup cross-section generation
Multigroup cross-section libraries are needed in performing neutronics calculations. These libraries are referred to as broad-group libraries. The number of energy groups and group structure are highly dependent on the application and/or user's objectives. For example, for shielding calculations, broad-group libraries such as SAILOR and BUGLE with 47-neutron and 20-gamma energy groups are used. The common procedure to obtain a broad-group library is a three-step process: (1) processing pointwise ENDF (PENDF) format cross sections; (2) generating fine-group cross sections; and (3) collapsing fine-group cross sections to broad-group. The NJOY code is used to prepare fine-group cross sections by processing pointwise ENDF data. The code has several modules, each one performing a specific task. For instance, the module RECONR performs linearization and reconstruction of the cross sections, and the module GROUPR generates multigroup self-shielded cross sections. After fine-group, i.e., groupwise ENDF (GENDF), cross sections are produced, cross sections are self-shielded, and a one-dimensional transport calculation is performed to obtain flux spectra at specific regions in the model. These fluxes are then used as weighting functions to collapse the fine-group cross sections to obtain a broad-group cross-section library. The third step described is commonly performed by the AMPX code system. SMILER converts NJOY GENDF filed to AMPX master libraries, AJAX collects the master libraries. BONAMI performs self-shielding calculations, NITAWL converts the AMPX master library to a working library, XSDRNPM performs one-dimensional transport calculations, and MALOCS collapses fine-group cross sections to broad-group. Finally, ALPO is used to generate ANISN format libraries. In this three-step procedure, generally NJOY requires the largest amount of CPU time. This time varies depending on the user's specified parameters for each module, such as reconstruction tolerances, temperatures
Nelson, Adam
Multi-group scattering moment matrices are critical to the solution of the multi-group form of the neutron transport equation, as they are responsible for describing the change in direction and energy of neutrons. These matrices, however, are difficult to correctly calculate from the measured nuclear data with both deterministic and stochastic methods. Calculating these parameters when using deterministic methods requires a set of assumptions which do not hold true in all conditions. These quantities can be calculated accurately with stochastic methods, however doing so is computationally expensive due to the poor efficiency of tallying scattering moment matrices. This work presents an improved method of obtaining multi-group scattering moment matrices from a Monte Carlo neutron transport code. This improved method of tallying the scattering moment matrices is based on recognizing that all of the outgoing particle information is known a priori and can be taken advantage of to increase the tallying efficiency (therefore reducing the uncertainty) of the stochastically integrated tallies. In this scheme, the complete outgoing probability distribution is tallied, supplying every one of the scattering moment matrices elements with its share of data. In addition to reducing the uncertainty, this method allows for the use of a track-length estimation process potentially offering even further improvement to the tallying efficiency. Unfortunately, to produce the needed distributions, the probability functions themselves must undergo an integration over the outgoing energy and scattering angle dimensions. This integration is too costly to perform during the Monte Carlo simulation itself and therefore must be performed in advance by way of a pre-processing code. The new method increases the information obtained from tally events and therefore has a significantly higher efficiency than the currently used techniques. The improved method has been implemented in a code system
Highlights: • Multigroup nuclear data are sampled based on multivariate normal distributions. • Multigroup perturbation factors are applied to pointwise-ACE nuclear data. • Samples of perturbed pointwise-ACE nuclear data are generated by NUSS for MCNPX. • Variances in MCNPX outputs due to perturbed samples of ACE data are quantified. • NUSS is verified with TSUNAMI and MCNPX PERT CARD sensitivity/uncertainty methods. - Abstract: Stochastic sampling (SS) method for quantifying nuclear data uncertainties is accomplished by using perturbed nuclear data in routine neutronics calculations and determining the variance of output parameters due to the input nuclear data uncertainties. Existing SS-based methods have demonstrated the feasibility and efficiency of propagating uncertainties in multigroup nuclear data. However, in fields such as criticality safety assessment, pointwise representation of nuclear data is more appropriate in order to corroborate the increasing safety demand and best-estimate modeling capabilities. In this work, an SS-based tool, called NUSS is implemented which perturbs pointwise ACE-formatted nuclear data using multigroup nuclear data covariance. The use of pointwise ACE-formatted nuclear data in NUSS can accommodate flexible multigroup covariance structures and allows for nuclear data uncertainty propagation through the continuous/pointwise-energy transport code MCNPX. As a first step of the NUSS development and verification, uncertainty contributions from 239Pu and 235U nuclear data were assessed for Jezebel (Pu-fueled) and Godiva (U-fueled) fast-spectrum criticality benchmarks. NUSS results are compared to those by other uncertainty quantification methods such as TSUNAMI and MCNPX PERT CARD. Next, Light Water Reactor (LWR) pin cell models from the OECD/NEA UAM Phase-1 benchmark were analyzed. Results of cross section and kinf uncertainties in consideration of different nuclear data covariance libraries are presented
Facilitated diffusion buffers noise in gene expression
Schoech, Armin; Zabet, Nicolae Radu
2014-01-01
Transcription factors perform facilitated diffusion (3D diffusion in the cytosol and 1D diffusion on the DNA) when binding to their target sites to regulate gene expression. Here, we investigated the influence of this binding mechanism on the noise in gene expression. Our results showed that, for biologically relevant parameters, the binding process can be represented by a two-state Markov model and that the accelerated target finding due to facilitated diffusion leads to a reduction in both ...
AIRDIF, Neutron and Gamma Doses from Nuclear Explosion by 2-D Air Diffusion
1 - Description of problem or function: AIRDIF is a two-dimensional atmospheric radiation diffusion code designed to calculate neutron and gamma doses in the environment of a nuclear explosion. It calculates radiation fluxes in one-dimensional homogeneous air, or two-dimensional variable density air. The results are limited by the assumptions inherent in diffusion theory: the region of interest must be large compared to the radiation mean free path, the spatial flux gradients must not be steep, flux varies linearly with the cosine of the direction angle. The code requires as input data neutron and gamma source spectra, coupled neutron-gamma multigroup cross sections, and, for two- dimensional problems, a set of mass integral scaling (MIS) coefficients. These latter are calculated from an AIRDIF output flux file for a one-dimensional problem by the auxiliary program MISFIT, using a least squares fitting technique to Murphy's radiation transmission equation. MISFIT can also be used to calculate one- dimensional MIS doses. The MIS coefficients and doses can be input to AIRDIF, in two- dimensional mode to calculate 2-D fluxes, doses and K-factors (the ratio of 2-D to 1-D dose). Alternatively the 2-D doses and K-factors may be computed using the output 2-D flux file of a previous AIRDIF run using the auxiliary program DOSCOMP. 2 - Method of solution: Un-collided particle flux is determined from an analytic expression describing exponential attenuation with distance. Diffusion theory is used for the flux, using un-collided flux as a source term. A central collided differencing technique is used to reduce the diffusion equation to a matrix equation, which is solved by the Successive Line Over-relaxation (SLOR) method. Total flux is calculated as the sum of collided and un-collided components. To maintain a mesh interval which has the same relationship to mean free path at all heights, an expanding non-orthogonal coordinate system is used. In homogeneous air this system
Greene, N.M.; Ford, W.E. III; Petrie, L.M.; Arwood, J.W.
1992-10-01
AMPX-77 is a modular system of computer programs that pertain to nuclear analyses, with a primary emphasis on tasks associated with the production and use of multigroup cross sections. AH basic cross-section data are to be input in the formats used by the Evaluated Nuclear Data Files (ENDF/B), and output can be obtained in a variety of formats, including its own internal and very general formats, along with a variety of other useful formats used by major transport, diffusion theory, and Monte Carlo codes. Processing is provided for both neutron and gamma-my data. The present release contains codes all written in the FORTRAN-77 dialect of FORTRAN and wig process ENDF/B-V and earlier evaluations, though major modules are being upgraded in order to process ENDF/B-VI and will be released when a complete collection of usable routines is available.
AMPX-77 is a modular system of computer programs that pertain to nuclear analyses, with a primary emphasis on tasks associated with the production and use of multigroup cross sections. AH basic cross-section data are to be input in the formats used by the Evaluated Nuclear Data Files (ENDF/B), and output can be obtained in a variety of formats, including its own internal and very general formats, along with a variety of other useful formats used by major transport, diffusion theory, and Monte Carlo codes. Processing is provided for both neutron and gamma-my data. The present release contains codes all written in the FORTRAN-77 dialect of FORTRAN and wig process ENDF/B-V and earlier evaluations, though major modules are being upgraded in order to process ENDF/B-VI and will be released when a complete collection of usable routines is available
The MGPRAKTINETs computer code for the BESM-6 computer intended for calculation of zone average trmal neutron group fluxes and functionals is described. The neutron spatial-energy distribution in a multizone cyllindrically-symmetric reactor cell is calculated by the operator splitting method. For the solution of the spatial part of the problem the method of surface pseudosources (Gsub(N)-approximation) in approximation of plane derivatives from the energy neutron current is employed. The energy part of the problem is solved in a multigroup approximation. Computer code efficiency has been demonstrated by calculation of two-zone cells with internal and external sources of the cell with on additional absorber and RBMK cell with reduction of the latter to cylindrical geometry. It is shown that the approximation of plane derivatives of neutron energy current allows calculating reactor cell characteristics with a sufficient for design calculations accuracy
Analyzing Average and Conditional Effects with Multigroup Multilevel Structural EquationModels
Axel Mayer
2014-04-01
Full Text Available Conventionally, multilevel analysis of covariance (ML-ANCOVA has been therecommended approach for analyzing treatment effects in quasi-experimental multilevel designswith treatment application at the cluster-level. In this paper, we introduce the generalizedML-ANCOVA with linear effect functions that identifies average and conditional treatment effectsin the presence of treatment-covariate interactions. We show how the generalized ML-ANCOVAmodel can be estimated with multigroup multilevel structural equation models that offerconsiderable advances compared to traditional ML-ANCOVA. The proposed model takes intoaccount measurement error in the covariates, sampling error in contextual covariates,treatment-covariate interactions, and stochastic predictors. We illustrate the implementation ofML-ANCOVA with an example from educational effectiveness research where we estimateaverage and conditional effects of early transition to secondary schooling on readingcomprehension.
ETOA, ABBN Multigroup Constants from ENDF/B for Fast Reactors
1 - Nature of physical problem solved: Production of ABBN type group constants up to 70 groups for fast reactor calculations, reading ENDF/B library as input. 2 - Method of solution: The multigroup method of Bondarenko et al. is used for processing basic nuclear data. Calculational algorithms for an unresolved resonance region are the same as those in the MC2 code. For a resolved resonance region, an ultrafine energy structure dependent on a level scheme is adopted. 3 - Restrictions on the complexity of the problem: Maximum number of: energy groups: 70; sigma0 values: 6; temperatures: 5. Self-shielding factors for an unrealistically low value of sigma0 are not guaranteed because of the approximations used in the unresolved resonance region
Recent validation experience with multigroup cross-section libraries and scale
This paper will discuss the results obtained and lessons learned from an extensive validation of new ENDF/B-V and ENDF/B-VI multigroup cross-section libraries using analyses of critical experiments. The KENO V. a Monte Carlo code in version 4.3 of the SCALE computer code system was used to perform the critical benchmark calculations via the automated SCALE sequence CSAS25. The cross-section data were processed by the SCALE automated problem-dependent resonance-processing procedure included in this sequence. Prior to calling KENO V.a, CSAS25 accesses BONAMI to perform resonance self-shielding for nuclides with Bondarenko factors and NITAWL-II to process nuclides with resonance parameter data via the Nordheim Integral Treatment
Approximate analytical solution of two-dimensional multigroup P-3 equations
Iterative solution of multigroup spherical harmonics equations reduces, in the P-3 approximation and in two-dimensional geometry, to a problem of solving an inhomogeneous system of eight ordinary first order differential equations. With appropriate boundary conditions, these equations have to be solved for each energy group and in each iteration step. The general solution of the corresponding homogeneous system of equations is known in analytical form. The present paper shows how the right-hand side of the system can be approximated in order to derive a particular solution and thus an approximate analytical expression for the general solution of the inhomogeneous system. This combined analytical-numerical approach was shown to have certain advantages compared to the finite-difference method or the Lie-series expansion method, which have been used to solve similar problems. (orig./RW)
A Method to Solve Multigroup P3 Equations in Cylindrical Geometry
To determine the space-energy distribution of thermal neutrons in a reactor cell a combination of the spherical harmonics method and multigroup procedure has been chosen. In P-3 approximation and cylindrical geometry such a scheme implies the solution of an inhomogeneous system of six ordinary first order differential equations. The general solution of the corresponding homogeneous system is known in analytical form. The present work shows how the free term of the system can be approximated in order to find a particular solution, and thus the general solution, of the inhomogeneous system. The procedure has been applied to calculate thermal spectra in a number of different reactor cells. Some results are presented and discussed. (author)
Approximate analytical solution of two-dimensional multigroup P-3 equations
Iterative solution of multigroup spherical harmonics equations reduces, in the P-3 approximation and in two-dimensional geometry, to a problem of solving an inhomogeneous system of eight ordinary first order differential equations. With appropriate boundary conditions, these equations have to be solved for each energy group and in each iteration step. The general solution of the corresponding homogeneous system of equations is known in analytical form. The present paper shows how the right-hand side of the system can be approximated in order to derive a particular solution and thus an approximate analytical expression for the general solution of the inhomogeneous system. This combined analytical-numerical approach was shown to have certain advantages compared to the finite-difference method or the Lie-series expansion method, which have been used to solve similar problems. (author)
Jones, Kelvyn; Johnston, Ron; Manley, David; Owen, Dewi; Charlton, Chris
2015-12-01
We develop and apply a multilevel modeling approach that is simultaneously capable of assessing multigroup and multiscale segregation in the presence of substantial stochastic variation that accompanies ethnicity rates based on small absolute counts. Bayesian MCMC estimation of a log-normal Poisson model allows the calculation of the variance estimates of the degree of segregation in a single overall model, and credible intervals are obtained to provide a measure of uncertainty around those estimates. The procedure partitions the variance at different levels and implicitly models the dependency (or autocorrelation) at each spatial scale below the topmost one. Substantively, we apply the model to 2011 census data for London, one of the world's most ethnically diverse cities. We find that the degree of segregation depends both on scale and group. PMID:26487190
Global dynamics of a novel multi-group model for computer worms
In this paper, we study worm dynamics in computer networks composed of many autonomous systems. A novel multi-group SIQR (susceptible-infected-quarantined-removed) model is proposed for computer worms by explicitly considering anti-virus measures and the network infrastructure. Then, the basic reproduction number of worm R0 is derived and the global dynamics of the model are established. It is shown that if R0 is less than or equal to 1, the disease-free equilibrium is globally asymptotically stable and the worm dies out eventually, whereas, if R0 is greater than 1, one unique endemic equilibrium exists and it is globally asymptotically stable, thus the worm persists in the network. Finally, numerical simulations are given to illustrate the theoretical results. (general)
On the feasibility of a homogenised multi-group Monte Carlo method in reactor analysis
The use of homogenised multi-group cross sections to speed up Monte Carlo calculation has been studied to some extent, but the method is not widely implemented in modern calculation codes. This paper presents a calculation scheme in which homogenised material parameters are generated using the PSG continuous-energy Monte Carlo reactor physics code and used by MORA, a new full-core Monte Carlo code entirely based on homogenisation. The theory of homogenisation and its implementation in the Monte Carlo method are briefly introduced. The PSG-MORA calculation scheme is put to practice in two fundamentally different test cases: a small sodium-cooled fast reactor (JOYO) and a large PWR core. It is shown that the homogenisation results in a dramatic increase in efficiency. The results are in a reasonably good agreement with reference PSG and MCNP5 calculations, although fission source convergence becomes a problem in the PWR test case. (authors)
The French 'CEA 86' multigroup cross-section library and its integral qualification
This paper describe the up-dated 99 groups library of the APOLLO French neutron computer code, the denominated 'CEA 86' library. The multigroup cross-section sets are based on the more recent nuclear data evaluations. The THEMIS code was generally used for the JEF-1 processing. In order to account for recent differential measurements and to improve the consistency between calculation and integral experiments, we produced our own CEA evaluations for the actinide nuclides: 235U, 238U, 239Pu, 240Pu, 241Am. This new APOLLO library was checked against critical experiments and PWR measurements: computed Conversion Factor, Reactivity Coefficients, Multiplication Factor, and Pu build-up are now in good agreement with LWR experimental results. PWR Pu recycling calculations, as does as HCLWR design studies, are also significantly improved. (author)
Release of the mtmg01ex NDI Neutron Multigroup Data Library
We have released the multi-temperature neutron multigroup transport library mtmg01ex, consisting of 181 isotope tables from mtmg01 and 18 element tables calculated from the isotope tables, all at 15 temperatures. These data, based primarily on the evaluations that produced the lanl2006 library, include gamma production and americium branching data. They were subjected to our standard production library testing. Because there are still known problems with and unanswered questions about multi-temperature data, including data size and load time issues, we do not recommend this data for general use; however, its quality is good enough for production release, and we request user help in addressing the remaining problems.
The solution of the multigroup neutron transport equation using spherical harmonics
A solution of the multi-group neutron transport equation in up to three space dimensions is presented. The flux is expanded in a series of unnormalised spherical harmonics. Using the various recurrence formulae a linked set of first order differential equations is obtained for the moments psisup(g)sub(lm)(r), γsup(g)sub(lm)(r). Terms with odd l are eliminated resulting in a second order system which is solved by two methods. The first is a finite difference formulation using an iterative procedure, secondly, in XYZ and XY geometry a finite element solution is given. Results for a test problem using both methods are exhibited and compared. (orig./RW)
The importance of accounting for resonance self-screening effects in multigroup cross sections when calculating fast reactors and neutron shields is considered. Formulae for averaging cross sections over resonance features with the account of anisotropy for scattering with large energy losses are derived. The model calculations of neutron fluxes have been performed for a U-H mixture (rhosub(H)/rhosub(U)=0.1), a U-Fe-H mixture and for the latter with rhosub(5)/rhosub(Fe)=0.01-0.5. It is concluded that in hydrogen-containing reactors the effect may be significant if the core contains iron in large quantities. The cross section averaging is considered for 3 systems: the KBR-2 critical assembly, spherical model of a large breeder, critical sphere of UO2 with 30% enrichment. The scattering anisotropy changes the multiplication factors of the first two systems by about 0.3%
Zhang, Xi; Showman, Adam P.
2015-11-01
Most of the current atmospheric chemistry models for planets (e.g., Krasnopolsky & Parshev 1981; Yung & Demore 1982; Yung, Allen & Pinto 1984; Lavvas et al. 2008; Zhang et al. 2012) and exoplanets (e.g., Line, Liang & Yung 2010; Moses et al. 2011; Hu & Seager 2014) adopt a one-dimensional (1D) chemical-diffusion approach in the vertical coordinate. Although only a crude approximation, these 1D models have succeeded in explaining the global-averaged vertical profiles of many chemical species in observations. One of the important assumptions of these models is that all chemical species are transported via the same eddy diffusion profile--that is, the assumption is made that the eddy diffusivity is a fundamental property of the dynamics alone, and does not depend on the chemistry. Here we show that, as also noticed in the Earth community (e.g., Holton 1986), this “homogenous eddy diffusion” assumption generally breaks down. We first show analytically why the 1D eddy diffusivity must generally depend both on the horizontal eddy mixing and the chemical lifetime of the species. This implies that the long-lived species and short-lived chemical species will generally exhibit different eddy diffusion profiles, even in a given atmosphere with identical dynamics. Next, we present tracer-transport simulations in a 2D chemical-diffusion-advection model (Shia et al. 1989; Zhang, Shia & Yung 2013) and a 3D general circulation model (MITgcm, e.g., Liu & Showman 2013), for both rapid-rotating planets and tidally-locked exoplanets, to further explore the effect of chemical timescales on the eddy diffusivity. From the 2D and 3D simulation outputs, we derive effective 1D eddy diffusivity profiles for chemical tracers exhibiting a range of chemical timescales. We show that the derived eddy diffusivity can depend strongly on the horizontal eddy mixing and chemistry, although the dependences are more complex than the analytic model predicts. Overall, these results suggest that
Multigroup Albedo Method applied to coupled neutron-gamma radiations shielding
Shielding calculations for neutron-gamma radiation are usually done by using the full Theory of Transport or the Monte Carlo Techniques. After some works based on the Albedo Method, the shielding calculations for neutron-gamma radiation have a reliable tool with great didactical value which shows its clarity and simplicity for the resolution of cases that involve neutrons and photon shielding in nonmultiplying media. The excellent results of these works have motivated the elaboration and the development of this study that will be presented in this dissertation. The balance of a neutronic current entering a shield of two layers considering the coupling neutron-gamma will be determined by the Albedo Method. The shield will be composed of a layer of iron and another one of manganese with 10 cm of thickness each. The arrays of the materials coefficients will be obtained from the ANISN code. ANISN is a one dimensional deterministic code that is based on transport equation. The final results obtained by the Albedo Method will be compared with the ANISN results for an order of angular quadrature S2. The angular quadrature S2 admits that the radiation has two routes in the same direction what better describes the Albedo Method behavior. The results obtained by using the Albedo Method show an excellent agreement with the values predicted by the adopted deterministic code ANISN. Due to the excellent results, the multigroup Albedo Method should be applied to the shielding calculations with multiple layers. In conclusion the multigroup Albedo Method has the great ability in solving shielding problems concerning to the Nuclear Engineering. (author)
XNWLUP, Graphical user interface to plot WIMS-D library multigroup cross sections
1 - Description of program or function: XnWlup is a computer program with user-friendly graphical interface to help the users of WIMS-D library to enable quick visualisation of the plots of the energy dependence of the multigroup cross sections of any nuclide of interest. This software enables the user to generate and view the histogram of 69 multi-group cross sections as a function of neutron energy under Microsoft Windows environment. This software is designed using Microsoft Visual C++ and Microsoft Foundation Classes Library. IAEA1395/05: New features of version 3.0: - Plotting absorption and fission cross sections of resonant nuclide after applying the self-shielding cross section. - Plotting the data of Resonant Integral table, as a function of dilution cross section for a selected temperature and for a given energy group. - Plotting the data of Resonant Integral table, as a function of temperature for a selected background dilution cross section and for a given energy group. - Clearing all the graphs except one graph from the display screen is easily done by using a tool bar button. - Displaying the coordinate of the cursor point with appropriate units. 2 - Methods: XnWlup helps to obtain histogram plots of the values of cross section data of an element/isotope available as 69-group WIMS-D library as a function of energy bins. The software XnWlup is developed with this graphical user interface in order to help those users who frequently refer to the WIMS-D library cross section data of neutron-nuclear reactions. The software also helps to produce handbook of WIMS-D cross sections
The discrete ordinates method is the most powerful and generally used deterministic method to obtain approximate solutions of the Boltzmann transport equation. However, as presently formulated, it is both restricted to orthogonal geometries and susceptible to producing ray effects. In this work, a finite element formulation, utilizing a canonical form of the transport equation, is developed to obtain both integral and pointwise solutions to neutron transport problems. To facilitate its application to nonorthogonal planar geometries, the formulation is based on the use of linear triangles. A general treatment of anisotropic scattering is included in the formulation by employing discrete ordinates like approximations. In addition, multigroup source outer iteration techniques are employed to perform group dependent calculations. The ability of the formulation to substantially reduce ray effects and its ability to perform streaming calculations are demonstrated by analyzing a series of test problems. The anisotropic scattering and multigroup treatments used in the development of the formulation are verified by a number of one-dimensional comparisons. These comparisons also demonstrate the relative accuracy of the formulation in predicting integral parameters. The applicability of the formulation to nonorthogonal planar geometries is demonstrated by analyzing a hexagonal type lattice. A small high leakage reactor model is analyzed to investigate the effects of varying both the spatial mesh and order of angular quadrature. This analysis reveals that these effects are more pronounced in the present formulation than in other conventional formulations. However, the insignificance of these effects is demonstrated by analyzing a realistic reactor configuration. In addition, this final analysis illustrates the importance of incorporating anisotropic scattering into the finite element formulation
ZZ AMZ, 70-Group 40 Isotope Multigroup Library for Fast Reactor Calculation
1 - Description of program or function: format: EXPANDA; number of groups: 70-group library of multigroup constants; nuclides: H-1, Be-9, B-10, B-11, C-12, O-16, N-23, Mg, Al-27, Si, Ti, V, Cr, Mn-55, Fe, Ni, Cu, Ga, Zr, Nb-93, Mo, In-115, Sn, Pb, Th-232, Pa-233, U-233, U-234, U-235, U-236, U-238, Pu-238, Pu-239, Pu-240, Pu-241, Pu-242, Am-241, and lumped fission products of U-233, U-235, Pu-239. origin: ENDF/B-IV and ENDF/B-V; weighting spectrum: Fission products inventories for BBR reactor at 360 and 600 days of irradiation were calculated and used as weighting function. AMZ is a 70-group library of multigroup constants for the fast reactor nuclear design code EXPANDA. Data is stored for three temperatures (300 K, 900 K, 2100 K) and for seven background cross sections. The following isotopes are available: H1, Be9, B10, B11, C12, O16, N23, Mg, Al27, Si, Ti, V, Cr, Mn55, Fe, Ni, Cu, Ga, Zr, Nb93, Mo, In115, Sn, Pb, Th232, Pa233, U233, U234, U235, U236, U238, Pu238, Pu239, Pu240, Pu241, Pu242, Am241, and lumped fission products of U233, U235, Pu239. 2 - Method of solution: Nuclear cross sections, transfer matrices, and self-shielding factors were generated from ENDF/B-IV data using the codes NJOY (PSR-0171) and RGENDF
Polyspecific transporters from the organic anion transporting polypeptide (OATP/Oatp) superfamily mediate the uptake of a wide range of compounds. In zebrafish, Oatp1d1 transports conjugated steroid hormones and cortisol. It is predominantly expressed in the liver, brain and testes. In this study we have characterized the transport of xenobiotics by the zebrafish Oatp1d1 transporter. We developed a novel assay for assessing Oatp1d1 interactors using the fluorescent probe Lucifer yellow and transient transfection in HEK293 cells. Our data showed that numerous environmental contaminants interact with zebrafish Oatp1d1. Oatp1d1 mediated the transport of diclofenac with very high affinity, followed by high affinity towards perfluorooctanesulfonic acid (PFOS), nonylphenol, gemfibrozil and 17α-ethinylestradiol; moderate affinity towards carbaryl, diazinon and caffeine; and low affinity towards metolachlor. Importantly, many environmental chemicals acted as strong inhibitors of Oatp1d1. A strong inhibition of Oatp1d1 transport activity was found by perfluorooctanoic acid (PFOA), chlorpyrifos-methyl, estrone (E1) and 17β-estradiol (E2), followed by moderate to low inhibition by diethyl phthalate, bisphenol A, 7-acetyl-1,1,3,4,4,6-hexamethyl-1,2,3,4 tetrahydronapthalene and clofibrate. In this study we identified Oatp1d1 as a first Solute Carrier (SLC) transporter involved in the transport of a wide range of xenobiotics in fish. Considering that Oatps in zebrafish have not been characterized before, our work on zebrafish Oatp1d1 offers important new insights on the understanding of uptake processes of environmental contaminants, and contributes to the better characterization of zebrafish as a model species. - Highlights: • We optimized a novel assay for determination of Oatp1d1 interactors • Oatp1d1 is the first SLC characterized fish xenobiotic transporter • PFOS, nonylphenol, diclofenac, EE2, caffeine are high affinity Oatp1d1substrates • PFOA, chlorpyrifos
Popovic, Marta; Zaja, Roko [Laboratory for Molecular Ecotoxicology, Division for Marine and Environmental Research, Rudjer Boskovic Institute, Bijenicka 54, 10 000 Zagreb (Croatia); Fent, Karl [University of Applied Sciences Northwestern Switzerland, School of Life Sciences, Gründenstrasse 40, CH-4132 Muttenz (Switzerland); Swiss Federal Institute of Technology (ETH Zürich), Department of Environmental System Sciences, Institute of Biogeochemistry and Pollution Dynamics, CH-8092 Zürich (Switzerland); Smital, Tvrtko, E-mail: smital@irb.hr [Laboratory for Molecular Ecotoxicology, Division for Marine and Environmental Research, Rudjer Boskovic Institute, Bijenicka 54, 10 000 Zagreb (Croatia)
2014-10-01
Polyspecific transporters from the organic anion transporting polypeptide (OATP/Oatp) superfamily mediate the uptake of a wide range of compounds. In zebrafish, Oatp1d1 transports conjugated steroid hormones and cortisol. It is predominantly expressed in the liver, brain and testes. In this study we have characterized the transport of xenobiotics by the zebrafish Oatp1d1 transporter. We developed a novel assay for assessing Oatp1d1 interactors using the fluorescent probe Lucifer yellow and transient transfection in HEK293 cells. Our data showed that numerous environmental contaminants interact with zebrafish Oatp1d1. Oatp1d1 mediated the transport of diclofenac with very high affinity, followed by high affinity towards perfluorooctanesulfonic acid (PFOS), nonylphenol, gemfibrozil and 17α-ethinylestradiol; moderate affinity towards carbaryl, diazinon and caffeine; and low affinity towards metolachlor. Importantly, many environmental chemicals acted as strong inhibitors of Oatp1d1. A strong inhibition of Oatp1d1 transport activity was found by perfluorooctanoic acid (PFOA), chlorpyrifos-methyl, estrone (E1) and 17β-estradiol (E2), followed by moderate to low inhibition by diethyl phthalate, bisphenol A, 7-acetyl-1,1,3,4,4,6-hexamethyl-1,2,3,4 tetrahydronapthalene and clofibrate. In this study we identified Oatp1d1 as a first Solute Carrier (SLC) transporter involved in the transport of a wide range of xenobiotics in fish. Considering that Oatps in zebrafish have not been characterized before, our work on zebrafish Oatp1d1 offers important new insights on the understanding of uptake processes of environmental contaminants, and contributes to the better characterization of zebrafish as a model species. - Highlights: • We optimized a novel assay for determination of Oatp1d1 interactors • Oatp1d1 is the first SLC characterized fish xenobiotic transporter • PFOS, nonylphenol, diclofenac, EE2, caffeine are high affinity Oatp1d1substrates • PFOA, chlorpyrifos
NMR 1D-imaging of water infiltration into meso-porous matrices
It is shown that coupling nuclear magnetic resonance (NMR) 1D-imaging with the measure of NMR relaxation times and self-diffusion coefficients can be a very powerful approach to investigate fluid infiltration into porous media. Such an experimental design was used to study the very slow seeping of pure water into hydrophobic materials. We consider here three model samples of nuclear waste conditioning matrices which consist in a dispersion of NaNO3 (highly soluble) and/or BaSO4 (poorly soluble) salt grains embedded in a bitumen matrix. Beyond studying the moisture progression according to the sample depth, we analyze the water NMR relaxation times and self-diffusion coefficients along its 1D-concentration profile to obtain spatially resolved information on the solution properties and on the porous structure at different scales. It is also shown that, when the relaxation or self-diffusion properties are multimodal, the 1D-profile of each water population is recovered. Three main levels of information were disclosed along the depth-profiles. They concern (i) the water uptake kinetics, (ii) the salinity and the molecular dynamics of the infiltrated solutions and (iii) the microstructure of the water-filled porosities: open networks coexisting with closed pores. All these findings were fully validated and enriched by NMR cryo-poro-metry experiments and by performing environmental scanning electronic microscopy observations. Surprisingly, results clearly show that insoluble salts enhance the water progression and thereby increase the capability of the material to uptake water. (authors)
Lozano Montero, Juan Andrés; Aragonés Beltrán, José María; García Herranz, Nuria
2009-01-01
More accurate modelling of physical phenomena involved in present and future nuclear reactors requires a multi-scale and multi-physics approach. This challenge can be accomplished by the coupling of best-estimate core-physics, thermal-hydraulics and multi-physics solvers. In order to make viable that coupling, the current trends in reactor simulations are along the development of a new generation of tools based on user-friendly, modular, easily linkable, faster and more accurate codes to be i...
The Gain Properties of 1-D Active Photonic Crystal
无
2003-01-01
The terminology 'ID frequency'(w ID) is proposed after analyzing the 1D active photonic crystal based on the transfer matrix method. The relationship between wID and the structure parameters of the photonic crystal is investigated.
Popovic, Marta; Zaja, Roko; Fent, Karl; Smital, Tvrtko
2014-10-01
Polyspecific transporters from the organic anion transporting polypeptide (OATP/Oatp) superfamily mediate the uptake of a wide range of compounds. In zebrafish, Oatp1d1 transports conjugated steroid hormones and cortisol. It is predominantly expressed in the liver, brain and testes. In this study we have characterized the transport of xenobiotics by the zebrafish Oatp1d1 transporter. We developed a novel assay for assessing Oatp1d1 interactors using the fluorescent probe Lucifer yellow and transient transfection in HEK293 cells. Our data showed that numerous environmental contaminants interact with zebrafish Oatp1d1. Oatp1d1 mediated the transport of diclofenac with very high affinity, followed by high affinity towards perfluorooctanesulfonic acid (PFOS), nonylphenol, gemfibrozil and 17α-ethinylestradiol; moderate affinity towards carbaryl, diazinon and caffeine; and low affinity towards metolachlor. Importantly, many environmental chemicals acted as strong inhibitors of Oatp1d1. A strong inhibition of Oatp1d1 transport activity was found by perfluorooctanoic acid (PFOA), chlorpyrifos-methyl, estrone (E1) and 17β-estradiol (E2), followed by moderate to low inhibition by diethyl phthalate, bisphenol A, 7-acetyl-1,1,3,4,4,6-hexamethyl-1,2,3,4 tetrahydronapthalene and clofibrate. In this study we identified Oatp1d1 as a first Solute Carrier (SLC) transporter involved in the transport of a wide range of xenobiotics in fish. Considering that Oatps in zebrafish have not been characterized before, our work on zebrafish Oatp1d1 offers important new insights on the understanding of uptake processes of environmental contaminants, and contributes to the better characterization of zebrafish as a model species. PMID:25088042