Search NASA⌕ Search

SEARCH · Search NASA

Results for “triangularity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,333 records · Page 74

Critical Configuration and Physics Measurements for Assemblies of U(93.15)O 2 Fuel Rods (1.506-cm Pitch)

A series of critical experiments were completed from 1962–1965 at Oak Ridge National Laboratory’s (ORNL’s) Critical Experiments Facility (CEF) in support of the Medium-Power Reactor Experiments (MPRE) program. In the late 1950s, efforts were made to study “power plants for the production of electrical power in space vehicles.”(a) The MPRE program was a part of those efforts and studied the feasibility of a stainless-steel system, boiling potassium 1 MW(t), or about 140 kW(e), reactor. The program was carried out in [fiscal years] 1964, 1965, and 1966. A summary of the program’s effort was compiled in 1967.a The delayed critical experiments were a mockup of a small, potassium-cooled space power reactor for validation of reactor calculations and reactor physics methods. Initial experiments, performed in November and December of 1962, consisted of a core of unmoderated stainless-steel tubes, each containing 26 UO 2 fuel pellets, surrounded by a graphite reflector. Measurements were performed to determine critical reflector arrangements, relative fission-rate distributions, and cadmium ratio distributions. Subsequent experiments used beryllium reflectors and also measured the reactivity for various materials placed in the core. “The [assemblies were built] on [a] vertical assembly machine so that the movable part was the core and bottom reflector” (see Reference 1). The experiment studied in this evaluation was the second of the series and had the fuel rods in a 1.506-cm-triangular pitch. One critical configuration was found (see Reference 3). Once the critical configuration had been achieved, various measurements of reactivity, relative axial and radial activation rates of 235U,bc and cadmium ratios were performed. The cadmium ratio, reactivity, and activation rate measurements performed on the critical configuration are described in Sections 1.3, 1.4, and 1.7, respectively.

1.506-cm pitch↗

Measure Guideline: Optimizing the Configuration of Flexible Duct Junction Boxes

This measure guideline offers additional recommendations to heating, ventilation, and air conditioning (HVAC) system designers for optimizing flexible duct, constant-volume HVAC systems using junction boxes within Air Conditioning Contractors of America (ACCA) Manual D guidance (Rutkowski, H. Manual D -- Residential Duct Systems, 3rd edition, Version 1.00. Arlington, VA: Air Conditioning Contractors of America, 2009.). IBACOS used computational fluid dynamics software to explore and develop guidance to better control the airflow effects of factors that may impact pressure losses within junction boxes among various design configurations (Beach, R., Prahl, D., and Lange, R. CFD Analysis of Flexible Duct Junction Box Design. Golden, CO: National Renewable Energy Laboratory, submitted for publication 2013). These recommendations can help to ensure that a system aligns more closely with the design and the occupants' comfort expectations. Specifically, the recommendations described herein show how to configure a rectangular box with four outlets, a triangular box with three outlets, metal wyes with two outlets, and multiple configurations for more than four outlets. Designers of HVAC systems, contractors who are fabricating junction boxes on site, and anyone using the ACCA Manual D process for sizing duct runs will find this measure guideline invaluable for more accurately minimizing pressure losses when using junction boxes with flexible ducts.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The Physics of Basis For A Conservative Physics And Conservative Technology Tokamak Power Plant, ARIES-ACT2

The conservative physics and conservative technology tokamak power plant ARIES-ACT2 has a major radius of 9.75 m at aspect ratio of 4.0, strong shaping with elongation of 2.2 and triangularity of 0.63. The no wall βN reaches ~ 2.4, limited by n=1 external kink mode, and can be extended to 3.2 with a stabilizing shell behind the ring structure shield. The bootstrap current fraction is 77% with a q95 of 8.0, requiring about ~ 4.0 MA of external current drive. This current is supplied with 30 MW of ICRF/FW and 80 MW of negative ion NB. Up to 1.0 MA can be driven with LH with no wall, and 1.5 or more MA can be driven with a stabilizing shell. EC was examined and is most effective for safety factor control over ρ ~ 0.2-0.6 with 20 MW. The pedestal density is ~ 0.65x10 20 /m 3 and the temperature is ~ 9.0 keV. The H98 factor is 1.25, n/n Gr = 1.3, and the net power to LH threshold power is 1.3-1.4 in the flattop. Due to the high toroidal field and high central temperature the cyclotron radiation loss was found to be high depending on the first wall reflectivity.

(MHD) RF Heating↗

The Physics Basis For An Advanced Physics And Advanced Technology Tokamak Power Plant Configuration, ARIES-ACT1

The advanced physics and advanced technology tokamak power plant ARIES-ACT1 has a major radius of 6.25 m at aspect ratio of 4.0, toroidal field of 6.0 T, strong shaping with elongation of 2.2 and triangularity of 0.63. The broadest pressure cases reached wall stabilized β N ~ 5.75, limited by n=3 external kink mode requiring a conducting shell at b/a = 0.3, and requiring plasma rotation, feedback, and or kinetic stabilization. The medium pressure peaking case reached βN = 5.28 with BT = 6.75, while the peaked pressure case reaches β N < 5.15. Fast particle MHD stability shows that the alpha particles are unstable, but this leads to redistribution to larger minor radius rather than loss from the plasma. Edge and divertor plasma modeling show that about 75% of the power to the divertor can be radiated with an ITER-like divertor geometry, while over 95% can be radiated in a stable detached mode with an orthogonal target and wide slot geometry. The bootstrap current fraction is 91% with a q 95 of 4.5, requiring about ~ 1.1 MA of external current drive. This current is supplied with 5 MW of ICRF/FW and 40 MW of LHCD. EC was examined and is most effective for safety factor control over ρ ~ 0.2-0.6 with 20 MW. The pedestal density is ~ 0.9x10 20 /m 3 and the temperature is ~ 4.4 keV. The H 98 factor is 1.65, n/n Gr = 1.0, and the net power to LH threshold power is 2.8-3.0 in the flattop.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High Fidelity Computational Analysis of CO2 Trapping at Pore Scales (Final Technical Report)

With an alarming rise in carbon dioxide (CO2) emission from anthropogenic sources, CO2 sequestration has become an attractive choice to mitigate the emission. Some popular storage media for CO 2 are oil reservoirs, deep coal-bed, and deep oceanic-beds. These have been used for the long term CO 2 storage. Due to special lowering viscosity and surface tension property of CO 2 , it has been widely used for enhanced oil recovery. The sites for CO 2 sequestration or enhanced oil recovery mostly consist of porous rocks. Lack of knowledge of molecular mobility under confinement and molecule-surface interactions between CO2 and natural porous media results in generally governed by unpredictable absorption kinetics and total absorption capacity for injected fluids, and therefore, constitutes barriers to the deployment of this technology. Therefore, it is important to understand the flow dynamics of CO 2 through the porous microstructures at the finest scale (pore-scale) to accurately predict the storage potential and long-term dynamics of the sequestered CO 2 . This report discusses about pore-network flow modeling approach using variational method and analyzes simulated results this method simulations at pore-scales for idealized network and using Berea Sandstone CT scanned images. Variational method provides a promising way to study the kinetic behavior and storage potential at the pore scale in the presence of other phases. The current study validates variational solutions for single and two-phase Newtonian and single phase non-Newtonian flow through angular pores for special geometries whose analytical and/or empirical solutions are known. The hydraulic conductance for single phase flow through a triangular duct was also validated against empirical results derived from lubricant theory.

54 ENVIRONMENTAL SCIENCES↗

Investigations of surface structural, dynamical, and magnetic properties of systems exhibiting multiferroicity, and topological phases by helium scattering spectroscopies (Final Report)

We propose to investigate the surface structural, dynamics and magnetic properties of the novel class of topological insulator crystals, as well as crystals that exhibit multiferroicity, magnetoelectricity and thermoelectricity. Topological insulators (TIs) are a new class of insulators in which a bulk gap for electronic excitations is generated because of the strong spin-orbit coupling inherent to these systems. These materials are distinguished from ordinary insulators by the presence of gapless metallic surface states, resembling chiral edge modes in quantum Hall systems, but with unconventional spin textures. These exotic metallic states are formed by topological conditions that also render the electrons travelling on such surfaces insensitive to scattering by impurities. The electronic quasi-particles populating the topological surface state are Dirac fermions; they have a linear dispersion and thus are massless just like photons. We propose to investigate the interaction of these massless Dirac fermions with the massive lattice in the newly discovered crystals, Bi 2 Se 3 , Bi 2 Te 3 and Sb 2 Te 3 . We shall use inelastic helium beam scattering from surfaces to search for related signatures in surface phonon dispersions mappings that cover the entire surface Brillouin zone of these materials. Our recent investigations of the (001) surface of the multiferroic crystals (Li/Na)Cu 2 O 2 revealed an anomalous surface structural behavior where surface Cu$^{2+}$ row rise above the surface plane as the crystal was cooled. Subsequent worming revealed the onset of a thermally activated incommensurate surface phase, driven by the elevated rows. We are currently investigating the structure of the magnetic phases in these quasi-one-dimensional magnetic rows. Multiferroics are excellent candidates for large magnetoelectric response. We propose to extend this investigation to the class of delafossites which are also multiferroics and have been investigated as good candidates for thermoelectric power devices. They are also typical triangular lattice antiferromagnets with geometric magnetic frustration that leads to helimagnetic structures.

36 MATERIALS SCIENCE↗

Measure Guideline: Optimizing the Configuration of Flexible Duct Junction Boxes

This measure guideline offers additional recommendations to heating, ventilation, and air conditioning (HVAC) system designers for optimizing flexible duct, constant-volume HVAC systems using junction boxes within Air Conditioning Contractors of America (ACCA) Manual D guidance. IBACOS used computational fluid dynamics software to explore and develop guidance to better control the airflow effects of factors that may impact pressure losses within junction boxes among various design configurations. These recommendations can help to ensure that a system aligns more closely with the design and the occupants' comfort expectations. Specifically, the recommendations described herein show how to configure a rectangular box with four outlets, a triangular box with three outlets, metal wyes with two outlets, and multiple configurations for more than four outlets. Designers of HVAC systems, contractors who are fabricating junction boxes on site, and anyone using the ACCA Manual D process for sizing duct runs will find this measure guideline invaluable for more accurately minimizing pressure losses when using junction boxes with flexible ducts.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Stationary Liquid Fuel Fast Reactor

For effective burning of hazardous transuranic (TRU) elements of used nuclear fuel, a transformational advanced reactor concept named SLFFR (Stationary Liquid Fuel Fast Reactor) was proposed based on stationary molten metallic fuel. The fuel enters the reactor vessel in a solid form, and then it is heated to molten temperature in a small melting heater. The fuel is contained within a closed, thick container with penetrating coolant channels, and thus it is not mixed with coolant nor flow through the primary heat transfer circuit. The makeup fuel is semi- continuously added to the system, and thus a very small excess reactivity is required. Gaseous fission products are also removed continuously, and a fraction of the fuel is periodically drawn off from the fuel container to a processing facility where non-gaseous mixed fission products and other impurities are removed and then the cleaned fuel is recycled into the fuel container. A reference core design and a preliminary plant system design of a 1000 MWt TRU- burning SLFFR concept were developed using TRU-Ce-Co fuel, Ta-10W fuel container, and sodium coolant. Conservative design approaches were adopted to stay within the current material performance database. Detailed neutronics and thermal-fluidic analyses were performed to develop a reference core design. Region-dependent 33-group cross sections were generated based on the ENDF/B-VII.0 data using the MC2-3 code. Core and fuel cycle analyses were performed in theta-r-z geometries using the DIF3D and REBUS-3 codes. Reactivity coefficients and kinetics parameters were calculated using the VARI3D perturbation theory code. Thermo-fluidic analyses were performed using the ANSYS FLUENT computational fluid dynamics (CFD) code. Figure 0.1 shows a schematic radial layout of the reference 1000 MWt SLFFR core, and Table 0.1 summarizes the main design parameters of SLFFR-1000 loop plant. The fuel container is a 2.5 cm thick cylinder with an inner radius of 87.5 cm. The fuel container is penetrated by twelve hexagonal control assembly (CA) guide tubes, each of which has 3.0 mm thickness and 69.4 mm flat-to-flat outer distance. The distance between two neighboring CA guide tube is selected to be 26 cm to provide an adequate space for CA driving systems. The fuel container has 18181 penetrating coolant tubes of 6.0 mm inner diameter and 2.0 mm thickness. The coolant tubes are arranged in a triangular lattice with a lattice pitch of 1.21 cm. The fuel, structure, and coolant volume fractions inside the fuel container are 0.386, 0.383, and 0.231, respectively. Separate steel reflectors and B4C shields are used outside of the fuel container. Six gas expansion modules (GEMs) of 5.0 cm thickness are introduced in the radial reflector region. Between the radial reflector and the fuel container is a 2.5 cm sodium gap. The TRU inventory at the beginning of equilibrium cycle (BOEC) is 5081 kg, whereas the TRU inventory at the beginning of life (BOL) was 3541 kg. This is because the equilibrium cycle fuel contains a significantly smaller fissile fraction than the LWR TRU feed. The fuel inventory at BOEC is composed of 34.0 a/o TRU, 41.4 a/o Ce, 23.6 a/o Co, and 1.03 a/o solid fission products. Since uranium-free fuel is used, a theoretical maximum TRU consumption rate of 1.011 kg/day is achieved. The semi-continuous fuel cycle based on the 300-batch, 1- day cycle approximation yields a burnup reactivity loss of 26 pcm/day, and requires a daily reprocessing of 32.5 kg of SLFFR fuel. This yields a daily TRU charge rate of 17.45 kg, including a makeup TRU feed of 1.011 kg recovered from the LWR used fuel. The charged TRU-Ce-Co fuel is composed of 34.4 a/o TRU, 40.6 a/o Ce, and 25.0 a/o Co.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Focused Research Group in Correlated Electron and Complex Materials

While the remarkable physical properties of correlated and complex electronic materials hold great promise for technological applications, one of the key values of the research in this field is its profound impact on fundamental physics. The transition metal oxides, pnictides, and chalcogenides play a key role and occupy an especially important place in this field. The basic reason is that the outer shell of transition metals contains the atomic d-orbitals that have small spatial extent, but not too small to behave as localized orbtials. These d-electrons therefore have a small wave function overlap in a solid, e.g. in an octahedral environment, and form energy bands that are relatively narrow and on the scale of the short-range intra-atomic Coulomb repulsion (Hubbard U). In this intermediate correlation regime lies the challenge of the many-body physics responsible for new and unconventional physical properties. The study of correlated electron and complex materials represents both the challenge and the vitality of condensed matter and materials physics and often demands close collaborations among theoretical and experimental groups with complementary techniques. Our team has a track record and a long-term research goal of studying the unusual complexities and emergent behaviors in the charge, spin, and orbital sectors of the transition metal compounds in order to gain basic knowledge of the quantum electronic states of matter. During the funding period of this grant, the team continued their close collaborations between theory, angle-resolved photoemission spectroscopy, and scanning tunneling microscopy and made significant progress and contributions to the field of iron-based superconductors, copper-oxide high-temperature superconductors, triangular lattice transition metal oxide cobaltates, strontium ruthenates, spin orbital coupled iridates, as well as topological insulators and other topological quantum states of matter. These results include both new discoveries and the resolution to outstanding and unresolved issues. It should be emphasized that the DOE funding provided the crucial support for the close and meaningful collaborations of the focused research group that go far beyond simply putting the research papers from each group together. Indeed, the majority of the publications involved multiple PIs and collaborations between theory and experiments.

36 MATERIALS SCIENCE↗

Pulsed-field magnetization of frustrated S = 1/2 Cu(pyrimidine) 1.5 (H 2 O)(BF 4 ) 2

Cu(pym) 1.5 (H 2 O)(BF 4 ) 2 (pym = pyrimidine) was synthesized and its structure determined by synchrotron single crystal X-ray diffraction. The compound contains S = 1/2 Cu(II) ions arranged in a distorted triangular array (Fig. 1). Each Cu(II) ion is coordinated to three pym ligands, two weakly held BF 4 - anions and one H 2 O. To get a sense to the extent (i.e., strength) of possible frustrated exchange interactions in this new compound we measured the magnetization of Cu(pym) 1.5 (H 2 O)(BF 4 ) 2 in pulsed magnetic fields up to 60 T.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

P13631-E002PF: Pulsed field magnetostriction of Ba 2 CoTeO 6

Ba 2 CoTeO 6 is an insulating material consisting of two magnetic subsystems referred as A and B build of S=1/2 spins. Subsystem A is considered to be a Heisenberg-like antiferromagnet (AFM) on a triangular lattice and subsystem B is a J 1 -J 2 Ising-like AFM on a honey-comb lattice. The magnetic phase transitions are observed at T N1 = 12 K and T N2 =3K for A and B respectively. The application of magnetic fields unveils a rich phase diagram that varies depending on the direction of the applied field for either H ll c or H perp c. To date the phase diagram has been investigated by means of specific heat measurement up to 9T and susceptibility measurements with a SQUID magnetometer up to 7T. Magnetization measured in pulsed magnetic fields up to 60T at 1.3K and 4.2K reveal several steps and plateaux occurring at varying critical fields depending on the crystallographic direction. Common to the magnetization parallel c and perpendicular c is the saturation above 40T.

36 MATERIALS SCIENCE↗

Carbon Dioxide Shuttling Thermochemical Storage Using Strontium Carbonate

Phase I concludes with significant progress made towards the SunShot ELEMENTS goals of high energy density, high power density, and high temperature by virtue of a SrO/SrCO 3 based material. A detailed exploration of sintering inhibitors has been conducted and relatively stable materials supported by YSZ or SrZO 3 have been identified as the leading candidates. In 15 cycle runs using a 3 hour carbonation duration, several materials demonstrated energy densities of roughly 1500 MJ/m 3 or greater. The peak power density for the most productive materials consistently exceeded 40 MW/m 3 —an order of magnitude greater than the SOPO milestone. The team currently has a material demonstrating nearly 1000 MJ/m 3 after 100 abbreviated (1 hour carbonation) cycles. A subsequent 8 hour carbonation after the 100 cycle test exhibited over 1500 MJ/m 3 , which is evidence that the material still has capacity for high storage albeit with slower kinetics. Kinetic carbonation experiments have shown three distinct periods: induction, kinetically-controlled, and finally a diffusion-controlled period. In contrast to thermodynamic equilibrium prediction, higher carbonation temperatures lead to greater conversions over a 1 hour periods, as diffusion of CO 2 is more rapid at higher temperatures. A polynomial expression was fit to describe the temperature dependence of the linear kinetically-controlled regime, which does not obey a traditional Arrhenius relationship. Temperature and CO 2 partial pressure effects on the induction period were also investigated. The CO 2 partial pressure has a strong effect on the reaction progress at high temperatures but is insignificant at temperatures under 900°C. Tomography data for porous SrO/SrCO 3 structures at initial stage and after multiple carbonation/decomposition cycles have been obtained. Both 2D slices and 3D reconstructed representations have been obtained from the raw data. Porosity and surface area have been computed based on tomography. The two parameters, permeability and Forchheimer constants, in the Darcy’s law for linear and non-linear range have been obtained from the fluid dynamics simulation. The structure is nearly isotropic from the hydrodynamics point of view. Thermal conductivity at elevated temperature has been obtained from the 3-D pore-scale simulations. The computed thermal conductivities of the structure compare well with an existing correlation. Computational models have been developed for solving the energy and mass transport equations with chemical reaction on the surface and inside the particles. The chemical kinetics or rate law is needed on the particle scale; however it is generally measured using a lumped model based on the average of reacted as well as unreacted mass. Many simulations have been conducted on a 2-D model to determine the effects of those pore scale parameters on the heat release based on the measured data. The simulation results using 2-D model show that average volumetric heat flow rate of 4 MW/m 3 in the first 5 minutes can be obtained. The insight gained from the 2-D simulation will be useful for the next stage simulation based on realistic 3-D structures. A probabilistic analysis has been conducted on the baseline system of recovering heat from the hot CO 2 evolved from the endothermic reaction and compressing the cooled CO 2 into a storage tank at 20-21 bar. Using a Monte Carlo simulation using 2,000,000 samples and both uniform and triangular probability distributions of four variables (heat exchanger pinch points, compressor efficiency, and power block efficiency), exergetic efficiency of 90% appears to be very likely.

14 SOLAR ENERGY↗

User Manual for the PROTEUS Mesh Tools

PROTEUS is built around a finite element representation of the geometry for visualization. In addition, the PROTEUS-SN solver was built to solve the even-parity transport equation on a finite element mesh provided as input. Similarly, PROTEUS-MOC and PROTEUS-NEMO were built to apply the method of characteristics on unstructured finite element meshes. Given the complexity of real world problems, experience has shown that using commercial mesh generator to create rather simple input geometries is overly complex and slow. As a consequence, significant effort has been put into place to create multiple codes that help assist in the mesh generation and manipulation. There are three input means to create a mesh in PROTEUS: UFMESH, GRID, and NEMESH. At present, the UFMESH is a simple way to generate two-dimensional Cartesian and hexagonal fuel assembly geometries. The UFmesh input allows for simple assembly mesh generation while the GRID input allows the generation of Cartesian, hexagonal, and regular triangular structured grid geometry options. The NEMESH is a way for the user to create their own mesh or convert another mesh file format into a PROTEUS input format. Given that one has an input mesh format acceptable for PROTEUS, we have constructed several tools which allow further mesh and geometry construction (i.e. mesh extrusion and merging). This report describes the various mesh tools that are provided with the PROTEUS code giving both descriptions of the input and output. In many cases the examples are provided with a regression test of the mesh tools. The most important mesh tools for any user to consider using are the MT_MeshToMesh.x and the MT_RadialLattice.x codes. The former allows the conversion between most mesh types handled by PROTEUS while the second allows the merging of multiple (assembly) meshes into a radial structured grid. Note that the mesh generation process is recursive in nature and that each input specific for a given mesh tool (such as .axial or .merge) can be used as “mesh” input for any of the mesh tools discussed in this manual.

97 MATHEMATICS AND COMPUTING↗

CUBIT: Geometry and Mesh Generation Toolkit (15.1 User Documentation)

CUBIT is a full-featured software toolkit for robust generation of two- and three-dimensional finite element meshes (grids) and geometry preparation. Its main goal is to reduce the time to generate meshes, particularly large hex meshes of complicated, interlocking assemblies. It is a solid-modeler based preprocessor that meshes volumes and surfaces for finite element analysis. Mesh generation algorithms include quadrilateral and triangular paving, 2D and 3D mapping, hex sweeping and multi-sweeping, tetrahedral meshing, and various special purpose primitives. CUBIT contains many algorithms for controlling and automating much of the meshing process, such as automatic scheme selection, interval matching, sweep grouping, and also includes state-of-the-art smoothing algorithms.

97 MATHEMATICS AND COMPUTING↗

Particle Mass Flow Control for High-Temperature Concentrating Solar Receivers

This report summarizes the results of a two-year project funded by the U.S. Department of Energy's Solar Energy Technologies Office (SuNLaMP 1506) to evaluate the performance of high-temperature (>700 °C) particle receivers for concentrating solar power (see Appendix A for project information). In the first year, novel particle release patterns were designed and tested to increase the effective solar absorptance of the particle curtain. Modeling results showed that increasing the magnitude and frequency of different wave-like patterns increased the effective absorptance and thermal efficiency by several percentage points, depending on the mass flow rate. Tests showed that triangular-wave, square-wave, and parallel-curtain particle release patterns could be implemented and maintained at flow rates of ~10 kg/s/m. The second year of the project focused on the development and testing of particle mass-flow control and measurement methods. An automated slide gate controlled by the outlet temperature of the particles was designed and tested. Testing demonstrated that the resolution accuracy of the slide-gate positioning was less than ~1 mm, and the speed of the slide gate enabled rapid adjustments to accommodate changes in the irradiance to maintain a desired outlet temperature range. Different in-situ particle mass-flow measurement techniques were investigated, and two were tested. The in-situ microwave sensor was found to be unreliable and sensitive to variations in particle flow patterns. However, the in-situ weigh hopper using load cells was found to provide reliable and repeatable measurements of real-time in-situ particle mass flow. On-sun tests were performed to determine the thermal efficiency of the receiver as a function of mass flow rate, particle temperature, and irradiance. Models of the tests were also developed and compared to the tests.

14 SOLAR ENERGY↗

A Low Order Mimetic Finite Difference Method for Resistive Magnetohydrodynamics in 2D

We developed a low order Mimetic Finite Difference (MFD) discretization for the equations of magnetohydrodynamics (MHD) in two space dimensions. These equations describe the evolution of the electric and magnetic fields in the presence of prescribed velocity field.The method is designed to work on general polygonal meshes and preserves the divergence-free condition on the magnetic field. The electric field is discretized at the vertices/nodes and the magnetic field uses edge-based discretization. The method reconstructs the magnetic field to extract nodal values necessary to approximate some terms present in Ohm’s law. We test the robustness of our numerical scheme on three different types of meshes: with triangular elements, quadrilateral and unstructured polyhedrons obtain from a Voronoi tesselation. Analysis of the convergence for each of the aforementioned mesh types is presented. We finish with a test problem that shows the method is capable of modelling magnetic reconnection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

PROTEUS-NODAL User Manual

The PROTEUS-NODAL code is being developed as a three-dimensional (3D) neutron transport code based on homogenized geometry models of nuclear reactors. The code at present includes two solution options: diffusion theory and a nodal simplified P 3 (SP 3 ). The diffusion solution option is based upon the DIF3D-VARIANT PN transport methodology where the transport options have not been implemented due to the desire to create a SN based solution methodology. The diffusion solver works for Cartesian, hexagonal, and triangular-z geometries, while the SP 3 solution option is currently available for hexagonal-z geometry only. Both diffusion and SP 3 capabilities provide comparable solutions to DIF3D-VARIANT when using similar approximations. This document provides brief descriptions of inputs and outputs of the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Steady-State Thermal-Hydraulic Analysis and Bowing Reactivity Evaluation Methods Based on Neutron and Gamma Transport Calculations

With the advancement of computer technology, the variational nodal transport code VARIANT is now routinely used for sodium-cooled fast reactor (SFR) design analysis. In particular, VARIANT transport calculations are usually performed for whole core neutronics analyses and for fuel cycle analyses with REBUS-3. However, the steady-state thermal hydraulics analysis with the SE2-ANL coed and the assembly bowing reactivity calculation with the structural analysis with the NUBOW-3D code are still performed with the DIF3D finite difference diffusion theory code. Furthermore, the DIF3D diffusion calculation for these applications is performed with only six triangular meshes per hexagonal assembly, and the SE2-ANL code determines the pin power distributions approximately by assuming that the power distribution within an assembly is separable in the radial and axial directions. In order to overcome the limitations and to improve the accuracy of the existing methods, it was proposed to develop the relevant computational methods based on the transport calculations with the VARIANT and PROTEUS-SN codes. The objectives of this study are to develop a new heating calculation method based on the coupled neutron and gamma transport calculations with VARIANT, to implement the new heating calculation method into the SE2-ANL steady-state thermal-hydraulics analysis code with improved numerical algorithms and additional capabilities for automatic flow allocation calculations, and to develop a method to calculate bowing reactivity coefficients based on the VARIANT core solutions and the PROTEUS-SN assembly solutions, which can be used in NUBOW-3D calculations. A new, coupled neutron and gamma heating calculation procedure has been developed based on VARIANT transport calculations. In the new heating calculation procedure, the neutron flux distribution is first determined by solving an eigenvalue transport problem using the VARIANT transport code. Using the calculated neutron flux, the intra-nodal gamma source distributions are calculated to be consistent with the neutron flux distribution with the GAMSOR code. For this, the GAMSOR code has been updated to generate the intra-nodal gamma source distributions in the form of the VARSRC dataset of VARIANT. With this gamma source distribution, the gamma flux distribution is determined by solving a fixed source gamma transport problem using VARIANT again. With the neutron and gamma flux distributions, the heat generation rates are evaluated at individual fuel pins and assembly duct walls, without assuming separable axial and radial profiles for the intra-assembly power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗