Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coding standards”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Understanding DERMS

With the rapid deployment of distributed energy resources (DER), there is a high level of interest in how these devices can be integrated with utility operations at all levels for management and monitoring purposes. This integration is challenging in that the number of devices is high and that ownership is often that of a customer or third party. Industry stakeholders first began to address DER integration by identifying and standardizing the functions that individual DER can perform autonomously, a distributed manner. Device-level functions like “voltage ride-through”, “volt-var”, “frequency-watt”, and “dynamic reactive current” were designed and documented and are now supported by communication standards and grid codes worldwide, making grid-supportive capabilities mandatory for new interconnections. This was a first step. A necessary step, but not sufficient to achieve end-to-end integration of DER with the grid. A substantial gap was recognized between the granular controls of individual DER and the type of organized services needed for grid support, integration with distribution management systems (DMS) and grid operations. In 2012, stakeholders began working to define a common set of grid-supportive services and means to integrate large quantities of DERs into utility systems in a way that is practical, sustainable and extensible. Because multiple parties could be involved; utilities, DMS providers, DER aggregators and facility/microgrid controller providers are working together.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2012

Disposition of Radioisotope Thermoelectric Generators Currently Located at the Oak Ridge National Laboratory - 12232

Under the American Recovery and Reinvestment Act (ARRA), the U.S. Department of Energy (DOE) awarded SEC Federal Services Corporation (SEC) a 34-building demolition and disposal (D and D) project at the Oak Ridge National Laboratory (ORNL) that included the disposition of six Strontium (Sr-90) powered Radioisotope Thermoelectric Generators (RTGs) stored outside of ORNL Building 3517. Disposition of the RTGs is very complex both in terms of complying with disposal facility waste acceptance criteria (WAC) and U.S. Department of Transportation (DOT) requirements for packaging and transportation in commerce. Two of the RTGs contain elemental mercury which requires them to be Land Disposal Restrictions (LDR) compliant prior to disposal. In addition, all of the RTGs exceed the Class C waste concentration limits under Nuclear Regulatory Commission (NRC) Waste Classification Guidelines. In order to meet the LDR requirements and Nevada National Security Site (NNSS) WAC, a site specific treatability variance for mercury was submitted to the U.S. Environmental Protection Agency (EPA) to allow macro-encapsulation to be an acceptable treatment standard for elemental mercury. By identifying and confirming the design configuration of the mercury containing RTGs, the SEC team proved that the current configuration met the macro-encapsulation standard of 40 Code of Federal Regulations (CFR) 268.45. The SEC Team also worked with NNSS to demonstrate that all radioisotope considerations are compliant with the NNSS low-level waste (LLW) disposal facility performance assessment and WAC. Lastly, the SEC team determined that the GE2000 Type B cask met the necessary size, weight, and thermal loading requirements for five of the six RTGs. The sixth RTG (BUP-500) required a one-time DOT shipment exemption request due to the RTG's large size. The DOT exemption justification for the BUP-500 relies on the inherent robust construction and material make-up of the BUP- 500 RTG. DOE-ORO, SEC, and the entire SEC RTG team are nearing the conclusion of the Sr-90 RTG disposition challenge - a legacy now 50 years in the making. Over 600,000 Ci of Sr-90 waste await disposal and its removal from ORNL will mark an historical moment in the clean-up of the cold-war legacy in the ORNL central industrial area. Elimination (i.e., removal) of the RTGs will reduce security risks at ORNL and disposal will permanently eliminate security risks. The RTGs will eventually decay to benign levels within a reasonable timeframe relative to radiological risks posed by long-lived isotopes. The safety authorization basis at ORNL Building 3517 will be reduced enabling greater operational flexibility in future clean-out and D and D campaigns. Upon disposition the Department of Energy will realize reduced direct and indirect surveillance and maintenance costs that can be reapplied to accelerated and enhanced clean-up of the Oak Ridge Reservation. At present, waste profiles for the RTGs are developed and under review by NNSS RWAP staff and approval authorities. Disposition schedule is driven by the availability of compliant shipping casks necessary to safely transport the RTGs from ORNL to NNSS. The first disposal of the RCA RTG is expected in April 2012 and the remaining RTGs disposed in 2012 and 2013. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

Secondary neutron yields from thick-target GCR accelerator experiments - 25434

This work overviews recent experiments conducted at the NASA Space Radiation Laboratory which sought to experimentally verify model predictions of an optimal spacecraft shielding thickness and provide data which can be used to quantify uncertainties in the physics models within the standard particle transport codes which presently, in some places, poorly predict data and whose results can vary significantly from code to code. Now that the multi-year experiment campaign is complete, some of the overarching trends in the data can be studied; though this work only provides a cursory look at these overarching trends. A small set of the neutron yields measured in these experiments are showcased here. The experiment employed two thick targets along a beamline to emulate the walls of a spacecraft being bombarded by galactic cosmic rays. Time-of-flight techniques were used to characterize neutrons produced in the upstream target while a novel pulse height deconvolution methodology was developed and implemented to characterize neutrons produced in the downstream target. (authors)

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

CMakePPLang: An object-oriented extension to CMake

CMakePPLang is an object-oriented extension to the CMake language written entirely using the original CMake language, with the goal of making projects built on CMake easier to create and maintain. That said, CMakePPLang has different coding practices, paradigms, and standards than the original CMake language, much in the same way that C++ coding differs from C coding despite some level of interoperability. Currently, CMakePPLang is used within the CMakePP organization as the foundation for two in-progress projects: CMakeTest and CMaize. CMakeTest provides a solution for unit testing CMake and CMakePPLang code. CMaize is a CMake tool to simplify interoperability between projects and writing their build systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Electric Vehicle Charging for Residential and Commercial Energy Codes: Technical Brief

Numerous studies show that sales of electric vehicles (EVs) have grown consistently over recent years in the U.S. The U.S. Energy Information Administration (EIA) estimated 3 million EVs were on the road in 2022, and the Edison Electric Institute (EEI) forecasts a total of 26.4 million EVs on the road by 2030. Based on this forecast, EEI projects the need for an additional 12.9 million EV charge ports by 2030. If EV charging infrastructure fails to keep pace with sales of EVs it could result in consumers stranded without options to power their vehicles. EVs are capable of providing substantial benefits to the consumers. EVs are less expensive to operate than conventional internal combustion engine vehicles, have lower maintenance costs, and have the convenience of fueling (charging) at home or work. Studies conducted in California show that costs associated with installing EV charging infrastructure can be substantially more expensive for retrofit scenarios compared to new construction, making inclusion of EV infrastructure in new construction codes a cost-effective policy option to increase infrastructure to meet growing demands. PNNL tracks adoption of mandatory EV provisions across the U.S. As of December 20, 2024, 12 states (California, Oregon, Washington, Colorado, New Mexico, Illinois, Maryland, Delaware, New Jersey, Rhode Island, Massachusetts and Vermont) and 53 local governments have added EV provisions to their building codes, local ordinances and zoning requirements. Originally published in 2022, this tech brief has been revised to align with recent model energy code committee discussions and published EV infrastructure code language. This technical brief summarizes market trends, costs and benefits, and provides sample code language for EV charging infrastructure for consideration to be included in model codes, such as the International Energy Conservation Code (IECC) and ANSI/ASHRAE/IES Standard 90.1, as well as directly by states and local governments in their building codes. The technical brief summarizes related efforts undertaken by states and local governments, and builds upon language considered during the 2021 and 2024 IECC development cycles.

2021 IECC↗
DOE OSTI2018

Energy Savings Analysis of the Proposed NYStretch-Energy Code 2018

This study was conducted by the Pacific Northwest National Laboratory (PNNL) in support of the stretch energy code development led by the New York State Energy Research and Development Authority (NYSERDA). In 2017 NYSERDA developed its 2016 Stretch Code Supplement to the 2016 New York State Energy Conservation Construction Code (hereinafter referred to as “NYStretch-Energy”). NYStretch-Energy is intended as a model energy code for statewide voluntary adoption that anticipates other code advancements culminating in the goal of a statewide Net Zero Energy Code by 2028. Since then, NYSERDA continues to develop the NYStretch-Energy Code 2018 edition. To support the effort, PNNL conducted energy simulation analysis to quantify the energy savings of proposed commercial provisions of the NYStretch-Energy Code (2018) in New York. The focus of this project is the 20% improvement over existing commercial model energy codes. A key requirement of the proposed stretch code is that it be ‘adoptable’ as an energy code, meaning that it must align with current code scope and limitations, and primarily impact building components that are currently regulated by local building departments. It is largely limited to prescriptive measures, which are what most building departments and design projects are most familiar with. This report describes a set of energy-efficiency measures (EEMs) that demonstrate 20% energy savings over ANSI/ASHRAE/IES Standard 90.1-2013 (ASHRAE 2013) across a broad range of commercial building types and all three climate zones in New York. In collaboration with New Building Institute, the EEMs were developed from national model codes and standards, high-performance building codes and standards, regional energy codes, and measures being proposed as part of the on-going code development process. PNNL analyzed these measures using whole building energy models for selected prototype commercial buildings and multifamily buildings representing buildings in New York. Section 2 of this report describes the analysis methodology, including the building types and construction area weights update for this analysis, the baseline, and the method to conduct the energy saving analysis. Section 3 provides detailed specifications of the EEMs and bundles. Section 4 summarizes the results of individual EEMs and EEM bundles by building type, energy end-use and climate zone. Appendix A documents detailed descriptions of the selected prototype buildings. Appendix B provides energy end-use breakdown results by building type for both the baseline code and stretch code in all climate zones.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2018

Final Energy Savings Analysis of the Proposed NYStretch-Energy Code 2018

This study was conducted by the Pacific Northwest National Laboratory (PNNL) in support of the stretch energy code development led by the New York State Energy Research and Development Authority (NYSERDA). In 2017 NYSERDA developed its 2016 Stretch Code Supplement to the 2016 New York State Energy Conservation Construction Code (hereinafter referred to as “NYStretch-Energy”). NYStretch-Energy is intended as a model energy code for statewide voluntary adoption that anticipates other code advancements culminating in the goal of a statewide Net Zero Energy Code by 2028. Since then, NYSERDA continues to develop the NYStretch-Energy Code 2018 edition (NYSERDA 2018). To support the effort, PNNL conducted energy simulation analysis to quantify the energy savings of proposed commercial provisions of the NYStretch-Energy Code 2018 in New York. The focus of this project is to significantly improve the energy efficiency over existing commercial model energy codes. A key requirement of the proposed stretch code is that it be ‘adoptable’ as an energy code, meaning that it must align with current code scope and limitations, and primarily impact building components that are currently regulated by local building departments. It is largely limited to prescriptive measures, which are what most building departments and design projects are most familiar with. This report describes a set of energy-efficiency measures (EEMs) that demonstrate significant energy savings over ANSI/ASHRAE/IES Standard 90.1-2013 (ASHRAE 2013) across a broad range of commercial building types and all three climate zones in New York. In collaboration with New Building Institute, the EEMs were developed from national model codes and standards, high-performance building codes and standards, regional energy codes, and measures being proposed as part of the on-going code development process. PNNL analyzed these measures using whole building energy models for selected prototype commercial buildings and multifamily buildings representing buildings in New York.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Commercial Zero Code Plug-In: Zero Energy and Operational Emissions Overlay for Model Energy Codes (Technical Brief)

Model energy codes (MEC) describe requirements needed to demonstrate building energy performance compliance. They can be readily adopted by states and local jurisdictions to support desirable energy efficiency investment in new buildings and major renovations. This report provides commercial building energy code language as an overlay to current MEC, which is recognized by the U.S. DOE as ASHRAE Standard 90.1 2022. The code language provides a performance-based compliance path for achieving net zero energy or net zero operational energy emissions buildings with MEC. The path requires meeting two compliance metric target values: 1) a required minimum level of efficiency and 2) a measure of zero energy or emissions. The supporting documentation illustrates the magnitude of needed efficiency improvement and the additional offsets to be achieved from renewable energy sources.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗
DOE OSTI2013

GRID-BASED EXPLORATION OF COSMOLOGICAL PARAMETER SPACE WITH SNAKE

We present a fully parallelized grid-based parameter estimation algorithm for investigating multidimensional likelihoods called Snake, and apply it to cosmological parameter estimation. The basic idea is to map out the likelihood grid-cell by grid-cell according to decreasing likelihood, and stop when a certain threshold has been reached. This approach improves vastly on the 'curse of dimensionality' problem plaguing standard grid-based parameter estimation simply by disregarding grid cells with negligible likelihood. The main advantages of this method compared to standard Metropolis-Hastings Markov Chain Monte Carlo methods include (1) trivial extraction of arbitrary conditional distributions; (2) direct access to Bayesian evidences; (3) better sampling of the tails of the distribution; and (4) nearly perfect parallelization scaling. The main disadvantage is, as in the case of brute-force grid-based evaluation, a dependency on the number of parameters, N{sub par}. One of the main goals of the present paper is to determine how large N{sub par} can be, while still maintaining reasonable computational efficiency; we find that N{sub par} = 12 is well within the capabilities of the method. The performance of the code is tested by comparing cosmological parameters estimated using Snake and the WMAP-7 data with those obtained using CosmoMC, the current standard code in the field. We find fully consistent results, with similar computational expenses, but shorter wall time due to the perfect parallelization scheme.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

WIMS/PANTHER analysis of UO{sub 2}/MOX cores using embedded super-cells

This paper describes a method of analysing PWR UO{sub 2}MOX cores with WIMS/PANTHER. Embedded super-cells, run within the reactor code, are used to correct the standard methodology of using 2-group smeared data from single assembly lattice calculations. In many other codes the weakness of this standard approach has been improved for MOX by imposing a more realistic environment in the lattice code, or by improving the sophistication of the reactor code. In this approach an intermediate set of calculations is introduced, leaving both lattice and reactor calculations broadly unchanged. The essence of the approach is that the whole core is broken down into a set of 'embedded' super-cells, each extending over just four quarter assemblies, with zero leakage imposed at the assembly mid-lines. Each supercell is solved twice, first with a detailed multi-group pin-by-pin solution, and then with the standard single assembly approach. Correction factors are defined by comparing the two solutions, and these can be applied in whole core calculations. The restriction that all such calculations are modelled with zero leakage means that they are independent of each other and of the core-wide flux shape. This allows parallel pre-calculation for the entire cycle once the loading pattern has been determined, in much the same way that single assembly lattice calculations can be pre-calculated once the range of fuel types is known. Comparisons against a whole core pin-by-pin reference demonstrates that the embedding process does not introduce a significant error, even after burnup and refuelling. Comparisons against a WIMS reference demonstrate that a pin-by-pin multi-group diffusion solution is capable of capturing the main interface effects. This therefore defines a practical approach for achieving results close to lattice code accuracy, but broadly at the cost of a standard reactor calculation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Impacts of Model Building Energy Codes

The Department of Energy (DOE) Building Energy Codes Program (BECP) periodically evaluates national and state-level impacts associated with energy codes in residential and commercial buildings. Pacific Northwest National Laboratory (PNNL), funded by DOE, conducted an assessment of the prospective impacts of national model building energy codes from 2010 through 2040. A previous PNNL study evaluated the impact of the Building Energy Codes Program. A 2016 study looked more broadly at overall code impacts and this report describes the methodology used for the assessment and presents the impacts in terms of energy savings, consumer cost savings, and reduced emissions at the state level and at aggregated levels. In 2021, DOE conducted an interim and limited update to its 2016 study to evaluate potential building code updates using the 2016 methodology. That interim update includes estimated savings resulting from updates to the model energy codes, including the ANSI/ASHRAE/IES Standard 90.1-2016 (ASHRAE 90.1-2016) and 2019 editions, as well as the 2018 and 2021 International Energy Conservation Code (IECC). In 2023, DOE developed a fully updated report that includes code updates (ASHRAE 90.1-2019 and 2021 IECC), as well as additional enhancements and updates, including updated energy prices, annual floorspace additions, state code adoption dates, and emission factors, among others. This current version is another fully updated report that includes code updates (ASHRAE 90.1-2022 and 2024 IECC), as well as additional enhancements and updates, including updated energy prices, state code adoption dates, emission factors, and renewable energy contribution among others. Energy codes follow a three-phase cycle that starts with the development of a new model code, proceeds with the adoption of the new code by states and local jurisdictions, and finishes when the new code is implemented and builders, architects, and engineers are required to comply with the new provisions. The development of new model code editions creates the potential for increased energy savings. After a new model code is adopted, potential savings are realized in the field when new buildings (or additions and alterations) are constructed to comply with the new code. The contributions of all three phases are crucial to the overall impact of codes and are considered in this assessment. Figure ES.1 schematically describes the analysis framework. Energy savings are expressed in terms of energy use intensity (EUI) in the figure.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2016

Building Energy Codes: Policy Overview and Good Practices

Globally, 32% of total final energy consumption is attributed to the building sector. To reduce energy consumption, energy codes set minimum energy efficiency standards for the building sector. With effective implementation, building energy codes can support energy cost savings and complementary benefits associated with electricity reliability, air quality improvement, greenhouse gas emission reduction, increased comfort, and economic and social development. This policy brief seeks to support building code policymakers and implementers in designing effective building code programs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2012

Understanding Energy Code Acceptance within the Alaska Building Community

This document presents the technical assistance provided to the Alaska Home Financing Corporation on behalf of PNNL regarding the assessment of attitudes toward energy codes within the building community in Alaska. It includes a summary of the existing situation and specific assistance requested by AHFC, the results of a questionnaire designed for builders surveyed in a suburban area of Anchorage, interviews with a lender, a building official, and a research specialist, and recommendations for future action by AHFC.

2009 International Energy Conservation Code↗
DOE OSTI2023

Cross-section measurement of the 82 Kr ($p,γ$) 83 Rb reaction in inverse kinematics

The total cross section of the 82 Kr(p,γ) 83 Rb reaction was measured for the first time at effective center-of-mass energies between 2.4 and 3.0 MeV, within the relevant Gamow window for the astrophysical γ process. The experiment took place at the National Superconducting Cyclotron Laboratory at Michigan State University using the ReA facility. A 82 Kr beam was directed onto a hydrogen gas cell located at the center of the Summing NaI(Tl) (SuN) detector. Here, the obtained spectra were analyzed using the γ-summing technique and the extracted cross section was compared to standard statistical model calculations using the non-smoker and talys codes. The comparison indicates that standard statistical model calculations tend to overproduce the cross section of the 82 Kr(p,γ) 83 Rb reaction relative to the experimentally measured values. Furthermore, the experimental data were used to provide additional constraints on the nuclear level density and the γ-ray strength function used in the statistical model calculations.

59 ≤ A ≤ 89↗
DOE OSTI2020

Construction of a Code Verification Matrix for Heat Conduction With Finite Element Code Applications

When establishing the pedigree of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model. During this process, numerical results on various meshes are systematically compared to a reference analytic solution. The selection of analytic solutions can be a laborious process, as it is difficult to establish adequate code confidence without performing redundant work. In this report we address this issue by applying a physics-based process that establishes a set of reference problems. In this process, code simulation options are categorized and systematically tested, which ensures that gaps in testing are easily identified and addressed. The resulting problems are primarily intended for code verification analysis but may also be useful for comparison to other simulation codes, troubleshooting activities, or training exercises. The process is used to select fifteen code verification problems relevant for the one-dimensional steady-state heat conduction equation. These problems are applicable to a wide variety of simulation tools, but, in this work, a demonstration is performed using the finite element-based nuclear fuel performance code BISON. Convergence to the analytic solution at the theoretical rate is quantified for a selection of the problems, which establishes a baseline pedigree for the code. Not only can this standard set of conduction solutions be used for verification of other codes, but also the physics-based process for selecting problems can be utilized to quantify and expand testing for any simulation tool.

42 ENGINEERING↗
DOE OSTI2016

Matter power spectrum and the challenge of percent accuracy

Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day N -body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used N -body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at k ≤1 h Mpc{sup −1} and to within three percent at k ≤10 h Mpc{sup −1}. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for k ≤ 2 h Mpc{sup −1}. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than L =0.5 h {sup −1}Gpc to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of M {sub p}=10{sup 9} h {sup −1}M{sub ⊙} is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗