Simulating the neutrino flux from the Spallation Neutron Source for the COHERENT experiment
Not Available
SEARCH · Engineering Papers
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.
Not Available
In this work we present fully analytical solutions for a class of finite, homogeneous, 1-D slab benchmark problems with nonlinear temperature feedback effects. The proposed class of benchmarks include multiplicative 1-D neutron transport (limited to quasistatic S{sub 2} with μ = ±1) coupled with thermal conduction, convection, Doppler broadening, and expansion effects along the length of the slab. This class of benchmark models, along with the corresponding analytical solutions, are valuable for validating multiphysics analysis frameworks that support coupled neutronics/thermal/structural calculations. Analytical solutions for the benchmarks are obtained by introducing an ansatz that the equilibrium flux and temperature distributions in the slab have the same shape. Specific values for the total microscopic cross section, σ{sub t,0}, and conductive heat transfer coefficient, h, that satisfy the assumed ansatz are then determined. A discussion of the procedure for generating benchmark models and analytical solutions is provided, along with numerical results for an example set of model parameters and thoughts on practical applications of the benchmarks for multiphysics code validation. (authors)
The present work follows a recent survey of airframe noise prediction methodologies. In that survey, Lighthill s acoustic analogy was identified as the most prominent analytical basis for current approaches to airframe noise research. Within this approach, a problem is typically modeled with the Ffowcs Williams and Hawkings (FW-H) equation, for which a geometry-independent solution is obtained by means of the use of the free-space Green function (FSGF). Nonetheless, the aeroacoustic literature would suggest some interest in the use of tailored or exact Green s function (EGF) for aerodynamic noise problems involving solid boundaries, in particular, for trailing edge (TE) noise. A study of possible applications of EGF for prediction of broadband noise from turbulent flow over an airfoil surface and the TE is, therefore, the primary topic of the present work. Typically, the applications of EGF in the literature have been limited to TE noise prediction at low Mach numbers assuming that the normal derivative of the pressure vanishes on the airfoil surface. To extend the application of EGF to higher Mach numbers, the uniqueness of the solution of the wave equation when either the Dirichlet or the Neumann boundary condition (BC) is specified on a deformable surface in motion. The solution of Lighthill s equation with either the Dirichlet or the Neumann BC is given for such a surface using EGFs. These solutions involve both surface and volume integrals just like the solution of FW-H equation using FSGF. Insight drawn from this analysis is evoked to discuss the potential application of EGF to broadband noise prediction. It appears that the use of a EGF offers distinct advantages for predicting TE noise of an airfoil when the normal pressure gradient vanishes on the airfoil surface. It is argued that such an approach may also apply to an airfoil in motion. However, for the prediction of broadband noise not directly associated with a trailing edge, the use of EGF does not appear to offer any advantages over the use of FSGF at the present stage of development. It is suggested here that the applications of EGF for airframe noise analysis be continued. As an example pertinent to airframe noise prediction, the Fast Scattering Code of NASA Langley is utilized to obtain the EGF numerically on the surface of a three dimensional wing with a flap and leading edge slat in uniform rectilinear motion. The interpretation and use of these numerical Green functions are then discussed.
Neural manifolds summarize the intrinsic structure of the information encoded by a population of neurons. Advances in experimental techniques have made simultaneous recordings from multiple brain regions increasingly commonplace, raising the possibility of studying how these manifolds relate across populations. However, when the manifolds are nonlinear and possibly code for multiple unknown variables, it is challenging to extract robust and falsifiable information about their relationships. We introduce a framework, called the method of analogous cycles, for matching topological features of neural manifolds using only observed dissimilarity matrices within and between neural populations. We demonstrate via analysis of simulations and in vivo experimental data that this method can be used to correctly identify multiple shared circular coordinate systems across both stimuli and inferred neural manifolds. Conversely, the method rejects matching features that are not intrinsic to one of the systems. Further, as this method is deterministic and does not rely on dimensionality reduction or optimization methods, it is amenable to direct mathematical investigation and interpretation in terms of the underlying neural activity. We thus propose the method of analogous cycles as a suitable foundation for a theory of cross-population analysis via neural manifolds.
The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-kinetic-ballooning mode (KBM) model scaling typically by at least 25%. Ion-scale ($k_yp_s$ < 1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δ$f$ gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β$_e$) is ∼60% below the KBM onset threshold and the electron temperature gradient mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β$_e$) and normalized electron density and temperature scale lengths,${a}$/${L_n}_e$, ${a}$/${L_T}_n$using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode, ion-temperature gradient mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local E × B shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.
To mitigate adverse effects from molten corium following a reactor pressure vessel failure (RPVF), some new reactor designs employ a core catcher and a sacrificial material (SM), such as ceramic or concrete, to stabilize the molten corium and avoid containment breach. Existing reactors cannot easily be modified to include these SMs but could be modified to allow injectable cooling materials. Current reactor designs are limited to using water to stabilize the corium, but this can create other issues such as reaction of water with the concrete forming hydrogen gas. Here the novel SM proposed is a granular carbonate mineral that can be used in existing light water reactor plants. The granular carbonate will decompose when exposed to heat, inducing an endothermic reaction to quickly solidify the corium in place and producing a mineral oxide and carbon dioxide. Corium spreading is a complex process strongly influenced by coupled chemical reactions, including decay heat from the corium, phase change, and reactions between the concrete containment and available water. A recently completed Sandia National Laboratories laboratory directed research and development (LDRD) project focused on two research areas: experiments to demonstrate the feasibility of the novel SM concept, and modeling activities to determine the potential applications of the concept to actual nuclear plants. Small-scale experiments using lead oxide (PbO) as a surrogate for molten corium demonstrate that the reaction of the SM with molten PbO results in a fast solidification of the melt due to the endothermic carbonate decomposition reaction and the formation of open pore structures in the solidified PbO from CO 2 released during the decomposition. A simplified carbonate decomposition model was developed to predict thermal decomposition of carbonate mineral in contact with corium. This model was incorporated into MELCOR, a severe accident nuclear reactor code. A full-plant MELCOR simulation suggests that by the introduction of SM to the reactor cavity prior to RPVF ex-vessel accident progression, e.g., core-concrete interaction and core spreading on the containment floor, could be delayed by at least 15 h; this may be enough for additional accident management to be implemented to alleviate the situation.
Appendix D of Title 40 Part 61 of the US Code of Federal Regulations (CFR) provides a procedure that US Department of Energy (US DOE) facility owners and operators can use to estimate radionuclide emissions to the atmosphere for dose calculations instead of measuring emissions for minor sources under the 40 CFR Part 61, Subpart H, National Emission Standards for Emissions of Radionuclides Other Than Radon From Department of Energy Facilities, regulation. The procedure assumes that any radioactive material heated above 100 °C is completely vaporized and emitted to the atmosphere. In 1991, the US DOE Oak Ridge Reservation (ORR) requested approval to use different release fractions (RFs) for uranium because of its high melting and boiling points. In response to the request, the US Environmental Protection Agency (US EPA) Region IV approved the use of modified RFs for elemental uranium provided no reaction had taken place to alter its chemical form. In 2015, the ORR requested approval to use different RFs for tungsten, again because of its high melting and boiling points. EPA Region IV approved the use of modified RFs for heated radioactive tungsten metal. In accordance with the two precedents set for heating uranium and radioactive tungsten metals, in 2016, the ORR requested approval to use modified RFs in a similar fashion for other radioactive solid metals and compounds with melting and boiling points above 500 °C that might be heated above 100 °C in future research projects and experiments. EPA Region IV again granted approval to use modified RFs for the list of compounds. This note discusses the proposed modified RFs and their development.
From the late 1960s to 1985, Lawrence Livermore National Laboratory (LLNL) undertook a series of ‘pulsed sphere’ experiments for 32 materials involving 148 different experiments using 75 different spheres to measure the neutron leakage spectra for different materials from 14-MeV neutrons generated by 3 H(d, n) 4 He reactions induced by an incident D+ beam from the Insulated Core Transformer (ICT) accelerator at LLNL. This deuteron beam was focused to impinge on tritium loaded onto a titanium substrate within a low mass target assembly placed within the center of a spherical shell constructed from the material of interest. The neutron time of flight was measured at a detector outside a collimator. The program was intended to be comprehensive with numerous targets of various thicknesses spanning the periodic table from H to 239 Pu. Five different neutron detectors were utilized to avoid possible systematic errors due to differences in detector responses. Further, these measurements were designed with the aim that the simple geometry of the measurements could be easily simulated to validate Monte Carlo transport codes and nuclear data libraries. As a result, these benchmarks have been extensively used in the validation of ENDL, ENDF and other libraries. Two experiments selected as benchmark cases in this evaluation are: 1) target assembly without material sphere (blank run), and 2) target assembly with a hollow polyethylene sphere, 1.8 mean free path (mfp) or 16.56 cm thick.
Ionizing radiation is a major environmental variable for cells on Earth, and so organisms have adapted to either prevent or to repair damages caused by it, primarily from the appearance and accumulation of reactive oxygen species (ROS). In this study, we measured the differential gene expression in Deinococcus radiodurans UWO298 cultures deprived of background ionizing radiation (IR) while growing 605 m underground at the Waste Isolation Pilot Plant (WIPP), reducing the dose rate from 72.1 to 0.9 nGy h –1 from control to treatment, respectively. This reduction in IR dose rate delayed the entry into the exponential phase of the IR-shielded cultures, resulting in a lower biomass accumulation for the duration of the experiment. The RNASeq-based transcriptome analysis showed the differential expression of 0.2 and 2.7% of the D. radiodurans genome after 24 and 34 h of growth in liquid culture, respectively. Gene expression regulation after 34 h was characterized by the downregulation of genes involved in folding newly synthesized and denatured/misfolded proteins, in the assimilation of nitrogen for amino acid synthesis and in the control of copper transport and homeostasis to prevent oxidative stress. We also observed the upregulation of genes coding for proteins with transport and cell wall assembly roles. These results show that D. radiodurans is sensitive to the absence of background levels of ionizing radiation and suggest that its transcriptional response is insufficient to maintain optimal growth.
Previously obtained results are extended to helium ambient gas to help in the assessment of performance in the atmospheres of the major planets. Measurements confirm benefits derived form detonating propellant over conventional chemical rocket propulsion. Benefits however, are reduced at high pressures in low molecular weight gas. Numerical calculations by means of a monodimensional hydrodynamic code follow the different trends obtained experimentally for high and low molecular weight gas and also offer new insights on the time behavior of the process.
Recent developments in the use of numerical methods for fluid flow simulation show an increasing tendency to use numerical operators that can, in various ways, be factored. Use of methods having this property often increases the accuracy and efficiency of computer codes. Tracing the factorization property provides a unity to the basic concepts involved in the development of cyclic reduction, predictor-corrector, splitting, fast Fourier transform, and pseudospectral methods.-
The author has identified the following significant results. The S190B color photography is as useful as LANDSAT data for the mapping of color differences in the rocks and soils of the terrain. An S192 ratio of 0.79 - 0.89 and 0.93 - 1.05 micron bands produced an apparently successful delineation of ferrous, ferric, and other materials, in agreement with theory and ratio code studies. From an analysis of S191 data, basalt and dacite were separated on the basis of differences in spectral emissivity in the 8.3 - 12 micron region.
The finite element method is used to investigate the static behavior of laminated composite flat plates and cylindrical shells. The analysis incorporates the effects of transverse shear deformation in each layer through the assumption that the normals to the undeformed layer midsurface remain straight but need not be normal to the mid-surface after deformation. A digital computer program was developed to perform the required computations. The program includes a very efficient equation solution code which permits the analysis of large size problems. The method is applied to the problem of stretching and bending of a perforated curved plate.
It is suggested that asymptotic results for the behavior of thin shells can be incorporated in a general computer code for the analysis of a complex shell structure. The advantage when compared to existing finite difference or finite element codes is a substantial reduction in computational labor with the capability of working to a specified level of accuracy. A reduction in user preparation time and dependance on user judgment is also gained, since mesh spacing can be internally generated. The general theory is described in this paper, as well as the implementation in the computer code FAST 1 (Functional Algorithm for Shell Theory) for the analysis of the general axisymmetric shell structure with axisymmetric loading.
A substantially modified version of the NASA-Lewis Chemical Equilibrium Code was recently developed. The modifications were designed to extend the power and convenience of the Code as a tool for performing combustor analysis for MHD systems studies. The effect of the programming details is described from a user point of view.
Any process which may be modeled statistically for a simple plasma particle may be modeled by the particle in cell technique. The success of the calculation is then dependent on having a large enough number of particles that the statistical treatment is meaningful. Thus it is possible to include the effects of secondary emission, backscattering, charge sticking and possibly dielectric breakdown, photoemission, and spallation. The first plasma dielectric interaction included in the computer code for simulating plasma insulator interactions is secondary electron emission. A calculated current density vs. voltage curve is presented and compared to an experimental curve.
The results of an experimental study of aerodynamic (surface velocity) and heat transfer distributions over the surfaces of two different, highly loaded, low-solidity contemporary turbine vane designs are presented. The aerodynamic configurations of the two vanes were carefully selected to emphasize fundamental differences in the character of the solution surface pressure distributions and the consequent effect on surface heat transfer distributions. The experimental measurements were made in moderate-temperature, three-vane cascades under steady-state conditions. The principal independent parameters (Mach number, Reynolds number, turbulence intensity, and wall-to-gas temperature ratio) were varied over ranges consistent with actual engine operation, and the heat matrix was structured to provide an assessment of the independent influence of each parameter. These measurements are intended to serve as verification data for a parallel analytical code development effort. The results of this parallel effort are briefly reviewed, and the principal conclusions to date are summarized.
The results of a numerical simulation of the time-averaged inviscid flow field through the blade rows of a multiblade row turboprop configuration are presented. The governing equations are outlined along with a discussion of the solution procedure and coding strategy. Numerical results obtained from a simulation of the flow field through a modern high-speed turboprop will be shown.