Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “C CODES”

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,063 records · Page 59

Sierra/SolidMechanics 5.0 Example Problems Manual

Presented in this document are tests that exist in the Sierra/SolidMechanics example problem suite, which is a subset of the Sierra / SM regression and performance test suite. These examples showcase common and advanced code capabilities. A wide variety of other regression and verification tests exist in the Sierra / SM test suite that are not included in this manual.

97 MATHEMATICS AND COMPUTING↗

Verification of RESRAD-OFFSITE Code (V.4)

This report documents the verification of RESRAD-OFFSITE Version 4.0 and describes, where necessary, the verification of the following: • The data comprising the standard dose and risk coefficient libraries in the RESRAD database files Master_dcf_ICRP07.mdb and Master_dcf_2k.mdb. • The extraction and transfer of the data from the selected database file to the computational code by the RESRAD-OFFSITE 4.0 interface, ResOWin.exe. • The different processes that are modeled by the main computational code in RESRAD OFFSITE 4.0, ResOMain.exe. • The data displayed in the graphical and text reports. Many verifications were performed as part of the quality assurance quality control program associated with the development and release of RESRAD-OFFSITE 4.0, namely: • developer testing, • internal independent testing, and • release testing. Some were also performed in response to questions from users regarding the performance of the code. The main text of the report focuses on summarizing a subset of those tests, both independent and developer tests that verified the computations performed by the code. The verifications included in this report served as the basis for the development of the release tests of the computational executables and provided the quantitative results to be compared with the code output. The input and output interfaces and the data transfers between the various executables of the code were tested while performing the verification testing. They were tested intentionally during release testing. This report also provides some basic information to help in understanding the activities that were verified. The report: • outlines the components of RESRAD-OFFSITE 4.0 and the interconnections between these components, • outlines the processes modeled by the computational code, • provides summary figures and tables to offer confirmation of the verification of the computational components of the code, • reproduces the verifiers’ reports, if available, in individual appendices, • refers to the previous verification report (Yu et al. 2011) for more details about some of the verifications, and • reproduces the test cases and the testers’ reports from the release testing in individual appendices, when possible.

54 ENVIRONMENTAL SCIENCES↗

New Grids of Pure-hydrogen White Dwarf NLTE Model Atmospheres and the HST/STIS Flux Calibration

Nonlocal thermodynamic equilibrium (NLTE) calculations of hot white dwarf (WD) model atmospheres are the cornerstone of modern flux calibrations for the Hubble Space Telescope and for the CALSPEC database. These theoretical spectral energy distributions (SEDs) provide the relative flux versus wavelength, and only the absolute flux level remains to be set by reconciling the measured absolute flux of Vega in the visible with the Midcourse Space Experiment values for Sirius in the mid-IR. The most recent SEDs calculated by the tlusty and tmap NLTE model atmosphere codes for the primary WDs G191-B2B, GD 153, and GD 71 show improved agreement to 1% from 1500 Å to 30 μm, in comparison to the previous 1% consistency only from 2000 Å to 5 μm. These new NLTE models of hot WDs now provide consistent flux standards from the far UV to the mid-IR.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining Physical Models at Gigabar Pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative information. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. The goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predictions of core plasma performance for the SPARC tokamak

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with Ion Cyclotron Range-of-Frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach Q~11, well above its mission objective of Q>2. To build confidence that SPARC will be successful, physics-based integrated modeling has also been performed. The TRANSP code coupled with the theory-based TGLF turbulence model and EPED predictions for pedestal stability find that Q~9 is attainable in standard H-mode operation and confirms Q>2 operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modeled with the Monte-Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Raw seedling and sapling census data collected from June 2015 - April 2019 within a tropical wet forest in Puerto Rico

This data package contains data used to evaluate the effects of experimental warming (4˚C above ambient), drought, and hurricane disturbance on species richness, diversity, and composition of understory plant communities across ontogeny at the Tropical Responses to Altered Climate Experiment (TRACE) in Luquillo, Puerto Rico. The raw data in this package include a file for the species composition matrix of woody seedlings (>10 cm and <20 cm height; TRACE_composition_seedling.csv) and saplings (>20 cm height; TRACE_composition_sapling.csv) within the six TRACE plots, collected during six census dates from June 2015 to April 2019. Also included is also a file with the species names associated with the species codes found in the composition files (TRACE_species_codes.csv), and a data dictionary file describing the columns in the composition files (data_dictionary_composition_files.csv). All data files are in csv format, with the exception of the metadata which is in xml format.

54 ENVIRONMENTAL SCIENCES↗

Heating effects on pyrogenic organic matter mass loss and temperature in a sand matrix in a lab setting from 2022

This dataset contains data associated with the preprint “A laboratory method for simulating the effects of subsequent fires on pyrogenic organic matter at different exposure depths in a sand matrix” (https://doi.org/10.1101/2024.07.30.605814), which has been accepted for publication in the International Journal of Wildland Fire. The dataset was generated to investigate the effects of subsequent fires on preexisting pyrogenic organic matter (PyOM) at different soil depths. Specifically, we designed a laboratory method to simulate soil heating from above, applying realistic heat flux profiles (High, Low, Control) derived from jack pine (Pinus banksiana Lamb.) log burns. Using a cone calorimeter, these profiles were applied to buried jack pine PyOM (produced at 350°C) in a sand matrix pyrocosm to examine PyOM mass loss at different depths (surface, 1 cm, 5 cm). Each treatment (e.g., High heat flux at 1 cm depth) was applied once to five replicates in a pyrocosm.The dataset includes mass loss fraction of PyOM and time-series temperature profiles (measured with type K thermocouples). Additionally, calibration datasets for high and low heat flux profiles are provided, including programmed heat flux profile data and three sets of calibration results. All datasets are in .csv format and were analyzed in R using code found at https://github.com/MengmengLuo/Reburning-pyrogenic-organic-matter-A-laboratory-method-for-dosing-dynamic-heat-fluxes-from-above. This dataset supports research on PyOM transformations during subsequent fires, and methodological advancements in simulating soil heating. It can be used to explore fire-induced PyOM alterations across different fire intensities and soil depths.

54 ENVIRONMENTAL SCIENCES↗

Hot Spot Evolution Measured by High-Resolution X-Ray Spectroscopy at the National Ignition Facility

Evolution of the hot spot plasma conditions was measured using high-resolution x-ray spectroscopy at the National Ignition Facility. The capsules were filled with DD gas with trace levels of Kr and had either a high-density-carbon (HDC) ablator or a tungsten (W)-doped HDC ablator. Time-resolved measurement of the Kr He β spectra, absolutely calibrated by a simultaneous time-integrated measurement, allows inference of the electron density and temperature through observing Stark broadening and the relative intensities of dielectronic satellites. By matching the calculated hot spot emission using a collisional-radiative code to experimental observations, the hot spot size and areal density are determined. Here these advanced spectroscopy techniques further reveal the effect of W dopant in the ablator on the hot spot parameters for their improved implosion performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design Guidance for High Temperature Concentrating Solar Power Components

This report provides guidance for the design of components for concentrating power facilities operating at high temperatures and undergoing high, secondary thermal stresses relative to the applied primary pressure stress. The design rules were developed for the design of Generation 3 CSP tubular receivers manufactured from Alloy 740H, but are generally applicable to a wide range of component types undergoing similar loads. Part 1 of the report provides procedural design rules for components to be used in conjunction with the 2019 edition of the ASME Boiler & Pressure Vessel Code. Part 2 provides corresponding design data for Alloy 740H. Part 3 of the report is a commentary describing the rationale behind the design rules and the data underlying the design material properties. Finally, Part 4 provides an extensive set of worked sample problems detailing the application of the rules to CSP components.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Sierra/SolidMechanics 4.56 Example Problems Manual

Presented in this document are tests that exist in the Sierra/SolidMechanics example problem suite, which is a subset of the Sierra/SM regression and performance test suite. These examples showcase common and advanced code capabilities. A wide variety of other regression and verification tests exist in the Sierra/SM test suite that are not included in this manual.

97 MATHEMATICS AND COMPUTING↗

Sierra/SolidMechanics 4.58 Example Problems Manual

Presented in this document are tests that exist in the Sierra/SolidMechanics example problem suite, which is a subset of the Sierra/SM regression and performance test suite. These examples showcase common and advanced code capabilities. A wide variety of other regression and verification tests exist in the Sierra/SM test suite that are not included in this manual.

97 MATHEMATICS AND COMPUTING↗

Sierra/SolidMechanics 5.2 Example Problems Manual

Presented in this document are tests that exist in the Sierra/SolidMechanics example problem suite, which is a subset of the Sierra/SM regression and performance test suite. These examples showcase common and advanced code capabilities. A wide variety of other regression and verification tests exist in the Sierra/SM test suite that are not included in this manual.

97 MATHEMATICS AND COMPUTING↗

Sierra/SolidMechanics 5.18: Example Problems Manual

Presented in this document are tests that exist in the Sierra/SolidMechanics example problem suite, which is a subset of the Sierra/SM regression and performance test suite. These examples showcase common and advanced code capabilities. A wide variety of other regression and verification tests exist in the Sierra/SM test suite that are not included in this manual.

97 MATHEMATICS AND COMPUTING↗

The Santa Barbara Binary-disk Code Comparison

We have performed numerical calculations of a binary interacting with a gas disk, using 11 different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, 0.2(a 2 Ω B /v) is the kinematic viscosity. The largest discrepancies between codes resulted from the dimensionality of the setup (3D vs. 2D disks). We find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.

79 ASTRONOMY AND ASTROPHYSICS↗

A Novel Two Stage Collimation Unit for Fermilab Booster

A new two-stage collimation unit (2SC) for Fermilab Booster will be installed during 2024 summer shutdown. It is a supplementary collimator for existing single stage Booster collimators. Design details of this 2SC adapted to Booster conditions are described. Results of beam dynamics simulations on collimation efficiency and evaluation of collimator shielding with MARS code on this new system are presented. The analysis on beam dynamics shows that the system will help to collimate up to about 1% of the beam in a confined way and, prompt and residual activation meet Fermilab Radiological Control limits in current as well as during PIP-II era high intensity operation of the Booster.

43 PARTICLE ACCELERATORS↗

Speeding-up fuzzing through directional seeds

Abstract Fuzzing is an automated process for discovering inputs in a program that may trigger unexpected behavior. Today, fuzzing has become a standard practice for the discovery of bugs and security vulnerabilities. However, the main issue with such practices is that the exploration of the input space of programs can often be prohibitively expensive. Therefore, several alternative fuzzing strategies have been introduced during the last few years. Some fuzzing techniques rely on human expertise to provide a plausible set of initial input examples, namely, seeds. However, the process of handcrafting seeds for fuzzing purposes often becomes strenuous for humans as it requires a deeper understanding of the Program-Under-Test (PUT). Also, the use of known inputs to programs often does not trigger vulnerable program behavior or may not reach potentially vulnerable code locations. To address those issues, we propose a seed generation framework that enables Human-In-The-Loop (HITL) directed fuzzing where the human assumes a more active role in the creation of seeds that can penetrate and assess desired locations of the PUT. Our proposed framework uses Symbolic Execution (SE) to generate seeds that exercise paths to target program locations. Moreover, our framework enables the visualization of the explored execution paths in the binary of the PUT for the generated seeds. We evaluated our approach on a set of 12 carefully designed C programs with diverse characteristics that mimic real-world programs. The experimental results show the effectiveness of the proposed approach in improving the performance of standard fuzzing tools such as the American Fuzzy Lop ("Image missing" <#comment/> ). Specifically, our solution can generate seeds that substantially enhance the performance of the fuzzer, achieving speedups ranging from $$1.46\times $$ 1.46 × to $$68.53\times $$ 68.53 × for branch conditions, $$1.39\times $$ 1.39 × to $$254.62\times $$ 254.62 × for branch depths, $$14,879.59\times $$ 14 , 879.59 × to $$30,295.88\times $$ 30 , 295.88 × for branch widths over traditional seeds. Additionally, the speedup increases with the number of target function ranging from $$12,260\times $$ 12 , 260 × to $$22,856.07\times $$ 22 , 856.07 × over traditional seeds while only requiring less than 15 seconds on average for the seed generation step.

97 MATHEMATICS AND COMPUTING↗

tapir: A tool for topologies, amplitudes, partial fraction decomposition and input for reductions

The demand for precision predictions in the field of high energy physics has dramatically increased over recent years. Experiments conducted at the LHC, as well as precision measurements at the intensity frontier such as Belle II require equally precise theoretical predictions to make full use of the acquired data. To match the experimental precision, second-, third- and, for certain quantities, even higher-order calculations in perturbative quantum field theory are required. To facilitate such calculations, computer software automating as many steps as possible is required. Yet, each calculation poses different challenges and thus, a high level of configurability is required. In this context we present tapir: a tool for identification, manipulation and minimization of Feynman integral families. It is designed to integrate in toolchains based on the computer algebra system FORM, the use of which is common practice in the field. tapir can be used to reduce the complexity of multi-loop problems with cut-filters, topology mapping, partial fraction decomposition and alike. Program Title:tapir CPC Library link to program files:https://doi.org/10.17632/ptc9t46xyn.1 Developer's repository link:https://gitlab.com/tapir-devs/tapir Licensing provisions: GPLv3 Programming language:python 3, C++ Nature of problem: Multi-loop computations require the automatization of a large number of different tasks related to Feynman integral topologies. Among them are the identification and minimization of integral topologies, partial fraction decomposition of topologies in the case of linearly dependent propagators as well as mapping scalar products of loop momenta to scalar functions. Solution method: The minimization of topologies is performed by comparison of their respective Nickel indices [1], even further minimization utilizes Pak's algorithm [2]. To efficiently map scalar products of loop momenta to scalar functions FORM [3] code is generated. Additional comments including restrictions and unusual features: Minimization based on Pak's algorithm slows down for many lines and scales. A coarser minimization using the Nickel indices, however, is still possible. [1]B. Nickel, D. Meiron, G.A.J. Baker, Compilation of 2-pt and 4-pt graphs for continuous spin model, Report, University of Guelph, 1977.[2]A. Pak, J. Phys. Conf. Ser. 368 (2012) 012049, https://doi.org/10.1088/1742-6596/368/1/012049, arXiv:1111.0868.[3]B. Ruijl, T. Ueda, J. Vermaseren, FORM version 4.2, arXiv:1707.06453, 7 2017.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Novel SOLPS-ITER simulations of X-point target and snowflake divertors

Abstract The design and understanding of alternative divertor configurations may be crucial for achieving acceptable steady-state heat and particle material loads for magnetic confinement fusion reactors. Multiple X-point alternative divertor geometries such as snowflakes and X-point targets have great potential in reducing power loads, but have not yet been simulated widely in codes with kinetic neutrals. This paper discusses recent changes made to the SOLPS-ITER code to allow for the simulation of X-point target and low-field side snowflake divertor geometries. Snowflake simulations using this method are presented, in addition to the first SOLPS-ITER simulation of the X-point target. Analysis of these results show reasonable consistency with the simple modelling and theoretical predictions, supporting the validity of the methodology implemented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗