Search NASA⌕ Search

SEARCH · Search NASA

Results for “reach code”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Improved Safety Injection Flow Map Associated with Target RCS Depressurization to Maintain Core Coolability of OPR1000

In Korea, Severe Accident Management Guidance (SAMG) for Optimized Power Reactor (OPR1000) has been developed through the research results related to the WOG SAMG, in which the operation phase shifts from the Emergency Operation Procedure (EOP) to the SAMG when the Core Exit Temperature (CET) reaches 923 K. Among 7 mitigation strategies envisioned, the third management strategy, injecting coolant into Reactor Coolant System (RCS) by depressurizing the RCS is one of the key strategies to fortify the in-vessel retention capability. As an aid to such mitigation strategy, a safety injection (SI) flow map has been developed by Lee et al. for predicting required flow rate associated with target RCS depressurization to recover in-vessel coolability. The mechanistic SI flow map used the CET increase rate and core water level decrease rate obtained from MELCOR code simulation. It should be noted that during the severe accident, the total heat generated is sum of the decay heat and the oxidation heat. Thus, it is necessary to consider additional heat sources by oxidation reaction between cladding, structure, control rod, and hot steam prior to the injection of coolant. It was believed that the SI model estimated required flow rate rather conservatively by considering the expected amount decay heat and the accumulated heat in the core based on the RCS pressure. However, it was confirmed that the SI model developed by Lee et al. could not represent the total heat removal especially since start of oxidation. Under transient conditions, the heat sources accumulated in the heat structures (fuel, cladding, supporting structures) should be embarked in the model for accurate contribution of the total heat source to the required flow rate. Therefore, the objective of this study is to improve the previous SI flow map predicting the injection flow rate associated with target RCS depressurization to satisfy the in-vessel coolability during severe accidents. Additional terms were developed to reflect accurate heat accumulating and included in the existing model. Resulting SI flow map includes the target RCS depressurization, which enables discharge of the required flow rate by utilizing the flow rate-pressure curve with operation of two high pressure safety injections (HPSIs). Also recent MELCOR simulation result of a postulated steam generator tube rupture (SGTR) accident verified the accuracy of the improved map.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

Improved Moderator Performance Calculations at SNS

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory reached a nominal operating power of 1 MW in September, 2009. Even at ~2/3 of its designed capacity, SNS already becomes the brightest pulsed neutron source in the world, and a premier facility for material and life science research. To provide accurate target performance estimation and to meet future development needs, we had developed an as-built model for the operating target system of SNS. Based on this asbuilt model, the moderator performance at SNS was recalculated taking advantage of the improved physics and tallying method in the Monte Carlo simulation code, MCNPX, and most recent scattering kernels. The neutron flux intensity spectra, emission time distributions of neutron pulses and moderator brightness profiles regarding the SNS moderator performance are presented and discussed. Compared to the previous results, the current calculations show ~10-20% lower in flux with improved emission time distributions of neutron pulses. The difference between the two data sets is analyzed, as well as other factors impacting the moderator performance at SNS. The current calculation results will serve as source terms for the instrument performance simulations at SNS.

Lu, Wei↗
DOE OSTI2015

Evaluation of noise pollution in urban traffic hubs—Noise maps and measurements

A study was made of some of the main traffic hubs in a Latin American metropolis, in order to determine the presence or absence of noise by means of noise measurements and acoustic mapping. To characterize noise in the evaluated road stretches, 232 measurements were taken at different points. The Predictor software package was used for the noise mapping calculations. Noise sensitive areas, e.g., hospitals, were identified in the evaluated road stretches. Noise maps were calculated for two hospitals, showing the current levels of noise that reach their facades. Hypothetical scenarios were simulated by making changes in the composition of traffic and total number of vehicles, and an assessment was made of the potential influence of these modifications in reducing the noise levels reaching the facades of the buildings in question. The simulations indicated that a 50% reduction in total traffic flow, or a 50% reduction in heavy vehicle traffic flow, would reduce the noise levels by about 3 dB(A). - Highlights: • Evaluation of noise pollution in urban traffic hubs • Street systems • Environmental noise impacts • Noise mapping.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

SN 2021adxl: A luminous nearby interacting supernova in an extremely low-metallicity environment

SN 2021adxl is a slowly evolving, luminous, Type IIn supernova with asymmetric emission line profiles, similar to the well-studied SN 2010jl. We present extensive optical, near-ultraviolet, and near-infrared photometry and spectroscopy covering ∼1.5 years post discovery. SN 2021adxl occurred in an unusual environment, atop a vigorously star-forming region that is offset from its host galaxy core. The appearance of Ly α and O II , as well as the compact core, would classify the host of SN 2021adxl as a “Blueberry” galaxy, analogous to higher redshift, low-metallicity, star-forming dwarf “Green Pea” galaxies. Using several abundance indicators, we find a metallicity of the explosion environment of only ∼0.1 Z ⊙ , the lowest reported metallicity for a Type IIn SN environment. SN 2021adxl reaches a peak magnitude of M r ≈ −20.2 mag and since discovery, SN 2021adxl has faded by only ∼4 magnitudes in the r band with a cumulative radiated energy of ∼1.5 × 10 50 erg over 18 months. SN 2021adxl shows strong signs of interaction with a complex circumstellar medium, seen by the detection of X-rays, revealed by the detection of coronal emission lines, and through multi-component hydrogen and helium profiles. In order to further understand this interaction, we model the H α profile using a Monte Carlo electron scattering code. The blueshifted high-velocity component is consistent with emission from a radially thin spherical shell resulting in the broad emission components due to electron scattering. Using the velocity evolution of this emitting shell, we find that the SN ejecta collide with circumstellar material of at least ∼5 M ⊙ assuming a steady-state mass-loss rate of ∼4 − 6 × 10 −3 M ⊙ yr −1 for the first ∼200 days of evolution. SN 2021adxl was last observed to be slowly declining at ∼0.01 mag d −1 , and if this trend continues, SN 2021adxl will remain observable after its current solar conjunction. Continuing the observations of SN 2021adxl may reveal signatures of dust formation or an infrared excess, similar to that seen for SN 2010jl.

Brennan, S. J. (ORCID:0000000313256235)↗
DOE OSTI2015

An Approach to Enhance pnetCDF Performance in Environmental Modeling Applications

I/O has been considered as a bottleneck in parallel applications. The software package, pnetCDF which works with parallel file systems, was developed to address this issue and provide parallel I/O capability. This study examines the performance of a novel approach which performs data aggregation along either row or column dimension of the spatial domain, and then applies the pnetCDF parallel I/O paradigm. The test was done with three different domain sizes which represents small, moderately large and large data domains, using a small scale Community Multi-scale Air Quality model (CMAQ) mocked up code. The examination includes comparing I/O performance with traditional serial I/O technique, straight application of pnetCDF, and the data aggregation along row and column dimension before applying pnetCDF. After the comparison, “optimal” I/O configurations for this new novel approach were quantified. Data aggregation along the row dimension (pnetCDFcr) works better than along the column dimension (pnetCDFcc) although it may perform slightly worse than straight the pnetCDF method with a small number of processors. When the number of processors becomes larger, pnetCDFcr out performs pnetCDF significantly. If the number of processors keeps increasing, pnetCDF reaches a point that the performance is even worse than the serial I/O technique. This new approach has also been tested on a real application where it performs two times better than the straight pnetCDF paradigm.

I/O↗
DOE OSTI2015

The Acoustic Emission signal acquired by the microphones placed in the CABRI test device along the fourteen last R.I.A. experiments: an example of reproducible research in nuclear science

The CABRI facility is an experimental nuclear reactor of the French Atomic Energy Commission. It is located at the Cadarache Research Centre in southern France and it is designed to act as a support to the French nuclear infrastructure. The purpose of the new testing programme termed, 'CABRI International Programme' (CIP) is to study the behaviour of PWR-type fuel rods at high burnup, equipped with an 'advanced' cladding, under Reactivity Initiated Accident (RIA) conditions (such as the scenario of a control rod ejection). Within the framework of this programme, piloted and funded by the French Institute of Nuclear Radioprotection and Safety (IRSN), ten tests are to be conducted with a frequency of two tests per year. The LPRE laboratory of the CEA which is in charge of the Preparation, realisation and breakdown of the test results studies the possibility to set up a new test analysis based on the processing of signals coming from sensors placed within the test equipment. During the experimental phase, the behaviour of the fuel element generates acoustic waves which can be detected by two microphones placed upstream and downstream of the test device. Studies showed the interest to realize temporal and spectral analyses on these signals by showing the existence of signatures which can be correlated with physical phenomena as the rod failure or the test shutdown (i.e. scram). The work presented in this article results from the will to consolidate these studies. Since the main phenomenon to be tracked is the fuel rod failure, the aim would be to highlight specific events which would have been precursors of the rod failure in order to use in the future these signals for further interpretation. Such an antecedent works resumption leads to a better understanding of the experimental needs and constitutes a good initial state to prepare the new very fast digital data acquisition systems. Although all the raw data are accessible in the form of text files, analyses and graphics representations were not straightforward to reproduce from the ancient studies since that, on one hand, people who were in charge of the original work left the laboratory and on the other hand because it is not easy when the time passes, even with our own work, to be able to remember the steps of data manipulations and the exact setup: - During the ancient experiments the use of analog data acquisition systems required to digitize tapes to be able to realize computer treatments. That had had for consequence to lose the initial dating. This one must be correctly edited to do temporal comparisons. - Analyses require functions for calculations whose parameters has to be well-known to reach the same results. We thus wished to manage our workflow in the idea that it can be easily reproducible on all the experiments. The object of the work presented in this article was to put in practice this strong bind between the data, treatments and generation of the document in order not to hesitate to do the iteration principle in action. We do not have to be afraid by the data driven analyses. According to the philosophy of the literate programming, the text of the technical document is woven with the computer code that produces all the printed output as tables, graphs for the study eliminating hence the unrealistic cut and paste. This difficulty is not specific to the nuclear domain. For many years, researchers have been worked out solutions to this mundane issue. And, presently, new technologies and high-level programming languages offer us actual answers. We will firstly present the tools applied in our laboratory to implement this workflow, then we will describe the global perception carried out to continue the study of the Acoustic Emission signals recorded by the two microphones during the fourteen last CABRI R.I.A. test.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Electron-optical systems for planar gyrotrons

The methodology of designing an electron-optical system (EOS) that forms sheet helical electron beams (HEBs) for high-power gyrotrons is developed. As an example, we consider the EOS for a 140-GHz gyrotron operated at the first harmonic of the cyclotron frequency with an accelerating voltage of 50 kV, a beam current of 30 A, and a magnetic field compression of 36. A planar geometry of the magnetron-injection gun (MIG) is suggested. The adiabatic theory of MIGs modified for the planar geometry of EOS is used for preliminary estimations of MIG parameters. Numerical simulation of the HEB properties based on the CST STUDIO SUITE 3D code is performed to find the optimal configuration of a planar MIG. The accuracy of the calculated data is discussed. The main factors that affect the HEB quality are considered. It is shown that a sheet HEB with a pitch-factor of 1.3 and velocity spread not exceeding 25%–30% can be formed; this is quite acceptable for high-efficiency operation of modern gyrotrons. Calculation of the beam-wave interaction with the obtained HEB parameters proved that a high output power with a sufficiently good efficiency of about 20% can be reached. Simulations show the feasibility of the experimental implementation of a novel planar EOS and its use in short-wave planar gyrotrons. The developed technique can be used for the study and optimization of planar gyrotrons of different frequency bands and power levels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

UBVR{sub c} I{sub c} ANALYSIS OF THE RECENTLY DISCOVERED TOTALLY ECLIPSING EXTREME MASS RATIO BINARY V1853 ORIONIS, AND A STATISTICAL LOOK AT 25 OTHER EXTREME MASS RATIO SOLAR-TYPE CONTACT BINARIES

We present precision CCD light curves, a period study, photometrically derived standard magnitudes, and a five-color simultaneous Wilson code solution of the totally eclipsing, yet shallow amplitude (A{sub v} {approx} 0.4 mag) eclipsing, binary V1853 Orionis. It is determined to be an extreme mass ratio, q = 0.20, W-type W UMa overcontact binary. From our standard star observations, we find that the variable is a late-type F spectral-type dwarf, with a secondary component of about 0.24 solar masses (stellar type M5V). Its long eclipse duration (41 minutes) as compared to its period, 0.383 days, attests to the small relative size of the secondary. Furthermore, it has reached a Roche lobe fill-out of {approx}50% of its outer critical lobe as it approaches its final stages of binary star evolution, that of a fast spinning single star. Finally, a summary of about 25 extreme mass ratio solar-type binaries is given.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Modeling jet and outflow feedback during star cluster formation

Powerful jets and outflows are launched from the protostellar disks around newborn stars. These outflows carry enough mass and momentum to transform the structure of their parent molecular cloud and to potentially control star formation itself. Despite their importance, we have not been able to fully quantify the impact of jets and outflows during the formation of a star cluster. The main problem lies in limited computing power. We would have to resolve the magnetic jet-launching mechanism close to the protostar and at the same time follow the evolution of a parsec-size cloud for a million years. Current computer power and codes fall orders of magnitude short of achieving this. In order to overcome this problem, we implement a subgrid-scale (SGS) model for launching jets and outflows, which demonstrably converges and reproduces the mass, linear and angular momentum transfer, and the speed of real jets, with ∼1000 times lower resolution than would be required without the SGS model. We apply the new SGS model to turbulent, magnetized star cluster formation and show that jets and outflows (1) eject about one-fourth of their parent molecular clump in high-speed jets, quickly reaching distances of more than a parsec, (2) reduce the star formation rate by about a factor of two, and (3) lead to the formation of ∼1.5 times as many stars compared to the no-outflow case. Most importantly, we find that jets and outflows reduce the average star mass by a factor of ∼ three and may thus be essential for understanding the characteristic mass of the stellar initial mass function.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Pre-Equilibrium De-Excitations for Neutrino-Nucleus Interactions

The Deep Underground Neutrino Experiment (DUNE) is sensitive to MeV-scale energy depositions from low-energy astrophysical neutrinos, including those from core-collapse supernovae. Interpreting these detector signals requires accurate modeling of the nuclear de-excitation that follows from the neutrino-nucleus interaction. The MARLEY (Model of Argon Reaction Low-Energy Yields) event generator specializes in the low-energy regime. MARLEY currently assumes the residual nucleus equilibrates immediately after the primary interaction. This omits the intermediate pre-equilibrium stage in which energy redistributes among nucleons until statistical equilibrium is reached. While pre-equilibrium effects are well established for nucleon-induced reactions, they have not previously been studied for neutrino-nucleus interactions. This work addresses that gap by implementing a two-component exciton model, which is the first dedicated treatment of pre-equilibrium de-excitation for neutrino-nucleus interactions, restructured around MARLEY's existing class hierarchy to prepare for direct integration, including particle-hole state densities, internal transition rates, and pre-equilibrium particle emission. The calculations show encouraging agreement with the TALYS-2.2 nuclear reaction code for neutron-nucleus interactions. We further propose a concrete integration path into the full MARLEY event generator, including derived class structure and an extended event record for pre-equilibrium vertices in support of future reweighting. Remaining work focuses on refining the emission width calculation, adding $\gamma$-ray emission, and completing this integration to quantify the impact of pre-equilibrium effects on the expected low-energy neutrino signals in DUNE and similar experiments.

Visser, Erin [Michigan State U., East Lansing (mai↗
DOE OSTI2024

Integrated modeling of boron powder injection for real-time plasma-facing component conditioning

An integrated modeling framework for investigating the application of solid boron (B) powder injection for real-time surface conditioning of plasma-facing components (PFCs) in tokamak environments is presented. Utilizing the DIII-D impurity powder dropper (IPD) setup, this study simulates B powder injection scenarios ranging from milligrams to tens of milligrams per second, corresponding to boron flux rates of 10 20 –10 21 B/s in standard L-mode conditions. The comprehensive modeling approach combines EMC3-EIRENE for simulating the deuterium plasma background and the Dust Injection Simulator (DIS) for the ablation and transport of the boron powder particles. EMC3 trace impurity fluid modeling results show substantial boron transport to the inboard lower divertor, predominantly influenced by the main ion plasma flow. The dependency on powder particle size (5-250 µm) was found to be insignificant for the scenario considered. The effects of erosion and redeposition were considered to reconcile the discrepancies with experimental observations, which saw substantial deposition on the outer divertor plasma-facing components. For this purpose, the WallDYN3D code was updated to include boron sources within the plasma domain and integrated into the modeling framework. The mixed-material migration modeling shows evolving boron deposition patterns, suggesting the formation of mixed B-C layers or predominantly B coverage depending on the powder mass flow rate. While the modeling outcomes at lower B injection rates tend to align with DIII-D experimental observations, the prediction of near-pure boron layers at higher rates has yet to be experimentally verified in the carbon environment of the DIII-D tokamak. The extensive reach of boron layers found in the modeling suggests the need for modeling that encompasses the entire wall geometry for more accurate experimental correlations. This integrated approach sets a precedent for analyzing and applying real-time in-situ boron coating techniques in advanced tokamak scenarios, potentially extendable to ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The structural, electronic and phonon behavior of CsPbI{sub 3}: A first principles study

Metal halide perovskites are optoelectronic materials that have attracted enormous attention as solar cells with power conversion efficiencies reaching 20%. The benefit of using hybrid compounds resides in their ability to combine the advantage of these two classes of compounds: the high mobility of inorganic materials and the ease of processing of organic materials. In spite of the growing attention of this new material, very little is known about the electronic and phonon properties of the inorganic part of this compounds. A theoretical study of structural, electronic and phonon properties of metal-halide cubic perovskite, CsPbI{sub 3} is presented, using first-principles calculations with planewave pseudopotential method as personified in PWSCF code. In this approach local density approximation (LDA) is used for exchange-correlation potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗
DOE OSTI2017

Calculations vs. Measurements for Remnant Dose Rates from SNS Spent Structures

The Spallation Neutron Source (SNS) in Oak Ridge, Tennessee, is an accelerator driven neutron scattering facility for materials research. Presently SNS is capable to operate at 1.4 MW proton beam power incident on a mercury target with a proton beam energy of 1 GeV and 60 Hz repetition rate. SNS target system components are periodically replaced because they reach their end-of-life due to radiation induced material damage. Target vessel, which houses mercury target, is exchanged about two-three times per year and the proton beam window (PBW) is exchanged every two – three years.Each spent structure that leaves the SNS site requires supporting documentation with radionuclide inventory and dose rate prediction for the time of the transportation. Neutronics analyses are performed, assuming realistic irradiation history and decay case to ensure that the container/package, housing the structure, is compliant with the waste management regulations. Analyses are complex due to geometry, multi-code usage and following data treatment.To validate analyses, measurements of dose rates from the spent target vessel # 13 and PBW module #5 were performed. Neutronics analyses were performed to calculate residual dose rates from both structures for the time of measurements.

Popova, Irina I.↗
DOE OSTI2014

Transport simulations of linear plasma generators with the B2.5-Eirene and EMC3-Eirene codes

Linear plasma generators are cost effective facilities to simulate divertor plasma conditions of present and future fusion reactors. For this research, the codes B2.5-Eirene and EMC3-Eirene were extensively used for design studies of the planned Material Plasma Exposure eXperiment (MPEX). Effects on the target plasma of the gas fueling and pumping locations, heating power, device length, magnetic configuration and transport model were studied with B2.5-Eirene. Effects of tilted or vertical targets were calculated with EMC3-Eirene and showed that spreading the incident flux over a larger area leads to lower density, higher temperature and off-axis profile peaking in front of the target. In conclusion, the simulations indicate that with sufficient heating power MPEX can reach target plasma conditions that are similar to those expected in the ITER divertor. B2.5-Eirene simulations of the MAGPIE experiment have been carried out in order to establish an additional benchmark with experimental data from a linear device with helicon wave heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A structural model of the long-term degradation of the concrete biological shield

The concrete biological shield (CBS) of light water reactors is exposed to high neutron radiation dose in the long term, which may lead to the degradation of the concrete’s mechanical properties. Given the important shielding role of the CBS, it is necessary to investigate the irradiation effects at the structural scale and provide estimates of the damage extent from the wall’s inner surface to study potential license renewals. For this purpose, we developed a mechanical model accounting for radiation-induced expansion, creep, and damage in concrete using the Grizzly finite element code, informed by ex-core neutron flux calculations using the VERA tool. The model was applied to a 3D CBS structure represented by the CBS wall, a steel liner, reinforcement bars, and a concrete base mat and evaluated damage at 40, 60, and 80 years of operation. The VERA model predicted a maximum fluence of approximately 2 x 10 19 ncm -2 at 80 years of operation. The results showed that damage is highest at the inner surface of the CBS wall and gradually decreases with depth. It extends beyond the rebar after 60 years and reaches a depth of approximately 12 cm at 80 years.

42 ENGINEERING↗
DOE OSTI2014

The coyote universe extended: Precision emulation of the matter power spectrum

Modern sky surveys are returning precision measurements of cosmological statistics such as weak lensing shear correlations, the distribution of galaxies, and cluster abundance. To fully exploit these observations, theorists must provide predictions that are at least as accurate as the measurements, as well as robust estimates of systematic errors that are inherent to the modeling process. In the nonlinear regime of structure formation, this challenge can only be overcome by developing a large-scale, multi-physics simulation capability covering a range of cosmological models and astrophysical processes. As a first step to achieving this goal, we have recently developed a prediction scheme for the matter power spectrum (a so-called emulator), accurate at the 1% level out to k ∼ 1 Mpc{sup –1} and z = 1 for wCDM cosmologies based on a set of high-accuracy N-body simulations. It is highly desirable to increase the range in both redshift and wavenumber and to extend the reach in cosmological parameter space. To make progress in this direction, while minimizing computational cost, we present a strategy that maximally reuses the original simulations. We demonstrate improvement over the original spatial dynamic range by an order of magnitude, reaching k ∼ 10 h Mpc{sup –1}, a four-fold increase in redshift coverage, to z = 4, and now include the Hubble parameter as a new independent variable. To further the range in k and z, a new set of nested simulations run at modest cost is added to the original set. The extension in h is performed by including perturbation theory results within a multi-scale procedure for building the emulator. This economical methodology still gives excellent error control, ∼5% near the edges of the domain of applicability of the emulator. A public domain code for the new emulator is released as part of the work presented in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗