A domain-specific analysis system for examining nuclear reactor simulation data for light-water and sodium-cooled fast reactors

Building a new generation of fission reactors in the United States presents many technical and regulatory challenges. One important challenge is the need to share and present results from new high-fidelity, high-performance simulations in an easily usable way. Since modern multiscale, multi-physics simulations can generate petabytes of data, they will require the development of new techniques and methods to reduce the data to familiar quantities of interest (e.g., pin powers, temperatures) with a more reasonable resolution and size. Furthermore, some of the results from these simulations may be new quantities for which visualization and analysis techniques are not immediately available in the community and need to be developed. This paper describes a new system for managing high-performance simulation results in a domain-specific way that naturally exposes quantities of interest for light water and sodium-cooled fast reactors. It describes requirements to build such a system and the technical challenges faced in its development at all levels (simulation, user interface, etc.). An example comparing results from two different simulation suites for a single assembly in a light-water reactor is presented, along with a detailed discussion of the system's requirements and design.Comment: Article on NiCE's Reactor Analyzer. 23 pages. Keywords: modeling, simulation, analysis, visualization, input-outpu

Structure of solar coronal loops: from miniature to large-scale

We will use new data from the High-resolution Coronal Imager (Hi-C) with unprecedented spatial resolution of the solar corona to investigate the structure of coronal loops down to 0.2 arcsec. During a rocket flight Hi-C provided images of the solar corona in a wavelength band around 193 A that is dominated by emission from Fe XII showing plasma at temperatures around 1.5 MK. We analyze part of the Hi-C field-of-view to study the smallest coronal loops observed so far and search for the a possible sub-structuring of larger loops. We find tiny 1.5 MK loop-like structures that we interpret as miniature coronal loops. These have length of the coronal segment above the chromosphere of only about 1 Mm and a thickness of less than 200 km. They could be interpreted as the coronal signature of small flux tubes breaking through the photosphere with a footpoint distance corresponding to the diameter of a cell of granulation. We find loops that are longer than 50 Mm to have a diameter of about 2 arcsec or 1.5 Mm, consistent with previous observations. However, Hi-C really resolves these loops with some 20 pixels across the loop. Even at this greatly improved spatial resolution the large loops seem to have no visible sub-structure. Instead they show a smooth variation in cross-section. The fact that the large coronal loops do not show a sub-structure at the spatial scale of 0.1 arcsec per pixel implies that either the densities and temperatures are smoothly varying across these loops or poses an upper limit on the diameter of strands the loops might be composed of. We estimate that strands that compose the 2 arcsec thick loop would have to be thinner than 15 km. The miniature loops we find for the first time pose a challenge to be properly understood in terms of modeling.Comment: Accepted for publication in A&A (Jun 19, 2013), 11 pages, 10 figure

Emergent singular solutions of non-local density-magnetization equations in one dimension

We investigate the emergence of singular solutions in a non-local model for a magnetic system. We study a modified Gilbert-type equation for the magnetization vector and find that the evolution depends strongly on the length scales of the non-local effects. We pass to a coupled density-magnetization model and perform a linear stability analysis, noting the effect of the length scales of non-locality on the system's stability properties. We carry out numerical simulations of the coupled system and find that singular solutions emerge from smooth initial data. The singular solutions represent a collection of interacting particles (clumpons). By restricting ourselves to the two-clumpon case, we are reduced to a two-dimensional dynamical system that is readily analyzed, and thus we classify the different clumpon interactions possible.Comment: 19 pages, 13 figures. Submitted to Phys. Rev.

Global-in-time solutions for the isothermal Matovich-Pearson equations

In this paper we study the Matovich-Pearson equations describing the process of glass fiber drawing. These equations may be viewed as a 1D-reduction of the incompressible Navier-Stokes equations including free boundary, valid for the drawing of a long and thin glass fiber. We concentrate on the isothermal case without surface tension. Then the Matovich-Pearson equations represent a nonlinearly coupled system of an elliptic equation for the axial velocity and a hyperbolic transport equation for the fluid cross-sectional area. We first prove existence of a local solution, and, after constructing appropriate barrier functions, we deduce that the fluid radius is always strictly positive and that the local solution remains in the same regularity class. To the best of our knowledge, this is the first global existence and uniqueness result for this important system of equations

Quantum Monte Carlo Studies of Relativistic Effects in Light Nuclei

Relativistic Hamiltonians are defined as the sum of relativistic one-body kinetic energy, two- and three-body potentials and their boost corrections. In this work we use the variational Monte Carlo method to study two kinds of relativistic effects in the binding energy of 3H and 4He. The first is due to the nonlocalities in the relativistic kinetic energy and relativistic one-pion exchange potential (OPEP), and the second is from boost interaction. The OPEP contribution is reduced by about 15% by the relativistic nonlocality, which may also have significant effects on pion exchange currents. However, almost all of this reduction is canceled by changes in the kinetic energy and other interaction terms, and the total effect of the nonlocalities on the binding energy is very small. The boost interactions, on the other hand, give repulsive contributions of 0.4 (1.9) MeV in 3H (4He) and account for 37% of the phenomenological part of the three-nucleon interaction needed in the nonrelativistic Hamiltonians.Comment: 33 pages, RevTeX, 11 PostScript figures, submitted to Physical Review

Detecting chaos in particle accelerators through the frequency map analysis method

The motion of beams in particle accelerators is dominated by a plethora of non-linear effects which can enhance chaotic motion and limit their performance. The application of advanced non-linear dynamics methods for detecting and correcting these effects and thereby increasing the region of beam stability plays an essential role during the accelerator design phase but also their operation. After describing the nature of non-linear effects and their impact on performance parameters of different particle accelerator categories, the theory of non-linear particle motion is outlined. The recent developments on the methods employed for the analysis of chaotic beam motion are detailed. In particular, the ability of the frequency map analysis method to detect chaotic motion and guide the correction of non-linear effects is demonstrated in particle tracking simulations but also experimental data.Comment: Submitted for publication in Chaos, Focus Issue: Chaos Detection Methods and Predictabilit

On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems

A consequent approach is proposed to construct symplectic force-gradient algorithms of arbitrarily high orders in the time step for precise integration of motion in classical and quantum mechanics simulations. Within this approach the basic algorithms are first derived up to the eighth order by direct decompositions of exponential propagators and further collected using an advanced composition scheme to obtain the algorithms of higher orders. Contrary to the scheme by Chin and Kidwell [Phys. Rev. E 62, 8746 (2000)], where high-order algorithms are introduced by standard iterations of a force-gradient integrator of order four, the present method allows to reduce the total number of expensive force and its gradient evaluations to a minimum. At the same time, the precision of the integration increases significantly, especially with increasing the order of the generated schemes. The algorithms are tested in molecular dynamics and celestial mechanics simulations. It is shown, in particular, that the efficiency of the new fourth-order-based algorithms is better approximately in factors 5 to 1000 for orders 4 to 12, respectively. The results corresponding to sixth- and eighth-order-based composition schemes are also presented up to the sixteenth order. For orders 14 and 16, such highly precise schemes, at considerably smaller computational costs, allow to reduce unphysical deviations in the total energy up in 100 000 times with respect to those of the standard fourth-order-based iteration approach.Comment: 23 pages, 2 figures; submitted to Phys. Rev.

Fourth Order Algorithms for Solving the Multivariable Langevin Equation and the Kramers Equation

We develop a fourth order simulation algorithm for solving the stochastic Langevin equation. The method consists of identifying solvable operators in the Fokker-Planck equation, factorizing the evolution operator for small time steps to fourth order and implementing the factorization process numerically. A key contribution of this work is to show how certain double commutators in the factorization process can be simulated in practice. The method is general, applicable to the multivariable case, and systematic, with known procedures for doing fourth order factorizations. The fourth order convergence of the resulting algorithm allowed very large time steps to be used. In simulating the Brownian dynamics of 121 Yukawa particles in two dimensions, the converged result of a first order algorithm can be obtained by using time steps 50 times as large. To further demostrate the versatility of our method, we derive two new classes of fourth order algorithms for solving the simpler Kramers equation without requiring the derivative of the force. The convergence of many fourth order algorithms for solving this equation are compared.Comment: 19 pages, 2 figure

Towards an experimental von Karman dynamo: numerical studies for an optimized design

Numerical studies of a kinematic dynamo based on von Karman type flows between two counterrotating disks in a finite cylinder are reported. The flow has been optimized using a water model experiment, varying the driving impellers configuration. A solution leading to dynamo action for the mean flow has been found. This solution may be achieved in VKS2, the new sodium experiment to be performed in Cadarache, France. The optimization process is described and discussed, then the effects of adding a stationary conducting layer around the flow on the threshold, on the shape of the neutral mode and on the magnetic energy balance are studied. Finally, the possible processes involved into kinematic dynamo action in a von Karman flow are reviewed and discussed. Among the possible processes we highlight the joint effect of the boundary-layer radial velocity shear and of the Ohmic dissipation localized at the flow/outer-shell boundary