DEVELOPMENT OF PTO-SIM: A POWER PERFORMANCE MODULE FOR THE OPEN-SOURCE WAVE ENERGY CONVERTER CODE WEC-SIM.
Abstract not provided.
Ruehl, Kelley Michelle↗
SEARCH · Engineering Papers
Results for “source 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.
HELIOS: AN OPEN-SOURCE, GPU-ACCELERATED RADIATIVE TRANSFER CODE FOR SELF-CONSISTENT EXOPLANETARY ATMOSPHERES
We present the open-source radiative transfer code named HELIOS, which is constructed for studying exoplanetary atmospheres. In its initial version, the model atmospheres of HELIOS are one-dimensional and plane-parallel, and the equation of radiative transfer is solved in the two-stream approximation with nonisotropic scattering. A small set of the main infrared absorbers is employed, computed with the opacity calculator HELIOS-K and combined using a correlated-k approximation. The molecular abundances originate from validated analytical formulae for equilibrium chemistry. We compare HELIOS with the work of Miller-Ricci and Fortney using a model of GJ 1214b, and perform several tests, where we find: model atmospheres with single-temperature layers struggle to converge to radiative equilibrium; k-distribution tables constructed with ≳0.01 cm{sup −1} resolution in the opacity function (≲10{sup 3} points per wavenumber bin) may result in errors ≳1%–10% in the synthetic spectra; and a diffusivity factor of 2 approximates well the exact radiative transfer solution in the limit of pure absorption. We construct “null-hypothesis” models (chemical equilibrium, radiative equilibrium, and solar elemental abundances) for six hot Jupiters. We find that the dayside emission spectra of HD 189733b and WASP-43b are consistent with the null hypothesis, while the latter consistently underpredicts the observed fluxes of WASP-8b, WASP-12b, WASP-14b, and WASP-33b. We demonstrate that our results are somewhat insensitive to the choice of stellar models (blackbody, Kurucz, or PHOENIX) and metallicity, but are strongly affected by higher carbon-to-oxygen ratios. The code is publicly available as part of the Exoclimes Simulation Platform (exoclime.net).
79 ASTRONOMY AND ASTROPHYSICS↗
Digital Controller For Emergency Beacon
Prototype digital controller intended for use in 406-MHz emergency beacon. Undergoing development according to international specifications, 406-MHz emergency beacon system includes satellites providing worldwide monitoring of beacons, with Doppler tracking to locate each beacon within 5 km. Controller turns beacon on and off and generates binary codes identifying source (e.g., ship, aircraft, person, or vehicle on land). Codes transmitted by phase modulation. Knowing code, monitor attempts to communicate with user, monitor uses code information to dispatch rescue team appropriate to type and locations of carrier.
Ivancic, William D.↗
Evaluation of the COHORT 2 Monte Carlo shielding code - Gamma ray transport using source direction biasing
Gamma ray transport portion of COHORT 2 Monte Carlo shielding code
Bloomfield, H. S.↗
Automation of FARM from Alpha Phase to Beta Phase
Integrated energy systems (IES) combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework and the RAVEN plugin for grid and capacity optimization (HERON) are used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES, by respecting the limits on the production variables and the corresponding rates of variation (explicit constraints). Besides, there are other process variables whose evolution needs to be bounded to avoid damaging the components (e.g., condensers, heat exchangers, steam generators, etc.) or degrading the process efficiency (e.g., electrolysis in the hydrogen production process). To avoid violating these latter limits (implicit constraints), a proof-of-concept HERON validator based on Feasible Actuator Range Modifier (FARM-Alpha) was developed by Argonne National Laboratory in January 2021. This FARM-Alpha validator calculates the evolution of process variables on whom the implicit constraints are placed, and then provides feedback to HERON dispatcher to adjust the power setpoints of three IES components, i.e., Balance of Plant, Secondary Energy Source, and Thermal Energy Storage, so as to meet both the explicit and implicit constraints. FARM-Alpha was designed to assess the performance of FARM as a HERON validator only, i.e., the list of components and implicit operational constraints were hard-coded within the source code. The lack of flexibility of the corresponding software structure does not allow the deployment in production environment. This report describes the development and the implementation of an enhanced version of the FARM-based validator (FARM-Beta), which ensures more flexibility for the end user in modeling multiple IES configurations and scenarios. Several test cases of the power dispatch problem were then selected to demonstrate the capabilities offered by FARM-beta. The test cases illustrate the efficiency of the closed-loop optimization scheme and the capability to calculate set-point trajectories satisfying both explicit and implicit constraints.
24 POWER TRANSMISSION AND DISTRIBUTION↗
MCNP6.3: A Year in Review [Slides]
This presentation discusses the past year and the many accomplishments of the MCNP6.3 team. It states that the MCNP6.3 release is imminent and that the approved, final documents have been making their way to the website. The code executables and source are already packaged up for distribution and the new installer is being finalized and tested now, for all platforms. The package will be sent to RSICC before the end of October 2022. Additionally, the presentation discusses some things to think about as MCNP6.3 is requested and/or used. For example, many of the recent efforts have been focused on making development, updates, and distribution of the code and documents more robust and streamlined. They will be revising/updating documents more frequently than ever before and they will be exploring avenues to distribute official patches to MCNP6.3. As they explain, this allows them to be more responsive to bugs and issues that are identified. In conclusion, the application of a patch to MCNP6.3 will require having the source code.
96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
HARD: A performance portable radiation hydrodynamics code based on FleCSI framework
Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.
97 MATHEMATICS AND COMPUTING↗
Numerical simulation of involute-plate research reactor flow behavior using RANS, LES and DNS
This paper investigates the flow behavior of involute-plate research reactors by performing Reynolds-Averaged Navier Stokes simulation (RANS), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) of the channel flow between fuel plates. By modeling turbulence with different numerical approaches, this study provides data with three levels of fidelity. For the RANS simulation, three widely used turbulence models, i.e., k-ε, k-ω, Reynolds Stress Turbulence model (RST) are applied by using the commercial CFD code STAR-CCM +. For LES and DNS, the open-source CFD code, Nek5000, is used given its outstanding scalability on High Performance Computer (HPC) and high-order technique. The results from RANS simulations are compared with that from LES and DNS for benchmarking. Both macroscale parameters and turbulence statistics, such as velocity magnitude, lateral velocity and turbulence kinetic energy, are presented and analyzed. The results from RANS simulation achieve good agreement with LES and DNS on velocity and turbulence kinetic energy prediction. The RST turbulence model predicts the most similar flow pattern of lateral velocity as compared to LES and DNS. The Lambda-2 (λ2) criterion with a reasonable threshold is used to demonstrate the instantaneous vortices distribution in the involute channel from both LES and DNS calculation. The DNS simulation captures more detailed turbulence especially near the corner, which explains the discrepancy between LES and DNS results near the corner. The normalized RMS error are defined and calculated to assess the performance of those turbulence models. The RST model captures the anisotropic feature of turbulence, which enable it to outperform other turbulence models for predicting the flow behavior in an involute channel. Although some discrepancies are found between LES and DNS results in the corner, the overall deviations between LES and DNS are found to be small. In conclusion, given that the computational cost of DNS calculation is an order of magnitude higher, using LES data for benchmarking RANS model is a cost-effective approach.
DNS↗
Fourier-Bessel Particle-In-Cell (FBPIC) v0.1.0
The Fourier-Bessel Particle-In-Cell code is a scientific simulation software for relativistic plasma physics. It is a Particle-In-Cell code whose distinctive feature is to use a spectral decomposition in cylindrical geometry. This decomposition allows to combine the advantages of spectral 3D Cartesian PIC codes (high accuracy and stability) and those of finite-difference cylindrical PIC codes with azimuthal decomposition (orders-of-magnitude speedup when compared to 3D simulations). The code is built on Python and can run both on CPU and GPU (the GPU runs being typically 1 or 2 orders of magnitude faster than the corresponding CPU runs.) The code has the exact same output format as the open-source PIC codes Warp and PIConGPU (openPMD format: openpmd.org) and has a very similar input format as Warp (Python script with many similarities). There is therefore tight interoperability between Warp and FBPIC, and this interoperability will increase even more in the future.
Lehe, Remi↗
An Analytical Tool for Fracture Control
The structural integrity of critical NASA space flight hardware is ensured through a damage tolerance program which is dictated in the fracture control requirements. The goal of a fracture control program is to ensure proper material selection and control, control of working stress levels, use of fracture resistant design concepts, manufacturing process control, and the use of qualified inspection procedures. The damage tolerance philosophy assumes that all components contain initial defects and the damage accumulation process is thus entirely crack propagation. The service capability of fracture critical components is determined by analysis and/or test. A large percentage of NASA fracture critical hardware is certified by analysis only and it is therefore imperative that the analysis code utilized contain the latest state of the art technology. Not all space contractors and vendors are able to develop and maintain an analysis code which meets this requirement. A response to this need has been fulfilled by the development of the NASA/FLAGRO (NASGRO) fracture mechanics code here at the Johnson Space Center. The NASGRO code is based primarily in Linear Elastic Fracture Mechanics (LEFM) and is continually updated to meet to needs of the user community. Development assistance comes from other NASA centers as well as the European Space Agency (ESA). The code currently has a vast materials data file (over 350 different entries) which defines the material strength, fracture toughness, crack growth rate, and crack growth threshold. The code also allows for user defined inputs of data from other sources. The code also provides a large library of stress intensity factor solutions for numerous geometries and loading conditions.
Patin, Raymond M.↗
V and V of the HTR Code Package (HCP) as a comprehensive HTR steady state and transient safety analysis framework - 18664
The HTR code package (HCP) allows for the simulation of several safety-related aspects of a High Temperature Reactor core in a highly integrated manner. HCP currently couples the thermo-fluid dynamics and time dependent neutronics code MGT, the spectrum code TRISHA, the burn up code TNT, the fuel management code SHUFLE and the fission product release code STACY. During its development, state-of-the-art programming techniques and standards are applied. The cross sections in HCP are generated by a 0D- or 1D solver featuring an innovative approach to treat the double heterogeneity of pebble fuel. Also major improvements have been made to optimize the nuclear data library based on ENDF/B-VII. Two advanced fuel management models are provided by the code module SHUFLE that go far beyond the capabilities of existing system codes like VSOP. The source term analysis code module STACY is coupled to HCP providing release rate calculations with a high spatial resolution making use of the nuclide densities provided by TNT. An outstanding new feature of HCP is the possibility to simulate long-term operation scenarios based on OTTO or MEDUL fuel shuffling schemes as well as selected transients in one integrated code package. Even alternating steady state/transient/steady state simulations are possible. So with HCP different fuel strategies and their influence on various kinds of accidents can be examined with one consistent reactor model. This paper provides an overview of the development status of the HCP and reports about selected benchmark results. It is demonstrated that the new system code HCP is capable to replace existing stand-alone codes like VSOP, TINTE/MGT, FRESCO or PANAMA while introducing new features, which so far to our knowledge were not available in the field of HTR safety research. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
(U) Release of SOURCES4D Featuring Sensitivity Capabilities
SOURCES is a code for computing neutron source rates and spectra from spontaneous fission (including delayed neutrons) and (α,n) reactions. The first release since 2002, SOURCES4D, is now available. The main updates include new output that allows the user to calculate, in post-processing, first and second derivatives of the (α,n) source rate density and spectrum with respect to nuclide densities in a homogeneous material and first derivatives of the (α,n) source rate density and spectrum with respect to nuclide stopping powers and (α,n) cross sections (nuclear data). The value of π was made more accurate in order to eliminate negative spontaneous-fission sources. Otherwise, SOURCES4D has made no changes to the physics or data of neutron source calculations used in SOURCES4C. First and second derivatives are calculated in an example problem.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Proton irradiation on materials
A computer code is developed by utilizing a radiation transport code developed at NASA Langley Research Center to study the proton radiation effects on materials which have potential application in NASA's future space missions. The code covers the proton energy from 0.01 Mev to 100 Gev and is sufficient for energetic protons encountered in both low earth and geosynchronous orbits. With some modification, the code can be extended for particles heavier than proton as the radiation source. The code is capable of calculating the range, stopping power, exit energy, energy deposition coefficients, dose, and cumulative dose along the path of the proton in a target material. The target material can be any combination of the elements with atomic number ranging from 1 to 92, or any compound with known chemical composition. The generated cross section for a material is stored and is reused in future to save computer time. This information can be utilized to calculate the proton dose a material would receive in an orbit when the radiation environment is known. It can also be used to determine, in the laboratory, the parameters such as beam current of proton and irradiation time to attain the desired dosage for accelerated ground testing of any material. It is hoped that the present work be extended to include polymeric and composite materials which are prime candidates for use as coating, electronic components, and structure building. It is also desirable to determine, for ground testing these materials, the laboratory parameters in order to simulate the dose they would receive in space environments. A sample print-out for water subject to 1.5 Mev proton is included as a reference.
Chang, C. Ken↗
YNOGK: A NEW PUBLIC CODE FOR CALCULATING NULL GEODESICS IN THE KERR SPACETIME
Following the work of Dexter and Agol, we present a new public code for the fast calculation of null geodesics in the Kerr spacetime. Using Weierstrass's and Jacobi's elliptic functions, we express all coordinates and affine parameters as analytical and numerical functions of a parameter p, which is an integral value along the geodesic. This is the main difference between our code and previous similar ones. The advantage of this treatment is that the information about the turning points does not need to be specified in advance by the user, and many applications such as imaging, the calculation of line profiles, and the observer-emitter problem, become root-finding problems. All elliptic integrations are computed by Carlson's elliptic integral method as in Dexter and Agol, which guarantees the fast computational speed of our code. The formulae to compute the constants of motion given by Cunningham and Bardeen have been extended, which allow one to readily handle situations in which the emitter or the observer has an arbitrary distance from, and motion state with respect to, the central compact object. The validation of the code has been extensively tested through applications to toy problems from the literature. The source FORTRAN code is freely available for download on our Web site http://www1.ynao.ac.cn/{approx}yangxl/yxl.html.
79 ASTRONOMY AND ASTROPHYSICS↗
Foundations of plasma standards
Abstract The field of low-temperature plasmas (LTPs) excels by virtue of its broad intellectual diversity, interdisciplinarity and range of applications. This great diversity also challenges researchers in communicating the outcomes of their investigations, as common practices and expectations for reporting vary widely in the many disciplines that either fall under the LTP umbrella or interact closely with LTP topics. These challenges encompass comparing measurements made in different laboratories, exchanging and sharing computer models, enabling reproducibility in experiments and computations using traceable and transparent methods and data, establishing metrics for reliability, and in translating fundamental findings to practice. In this paper, we address these challenges from the perspective of LTP standards for measurements, diagnostics, computations, reporting and plasma sources. This discussion on standards, or recommended best practices, and in some cases suggestions for standards or best practices, has the goal of improving communication, reproducibility and transparency within the LTP field and fields allied with LTPs. This discussion also acknowledges that standards and best practices, either recommended or at some point enforced, are ultimately a matter of judgment. These standards and recommended practices should not limit innovation nor prevent research breakthroughs from having real-time impact. Ultimately, the goal of our research community is to advance the entire LTP field and the many applications it touches through a shared set of expectations.
data and reaction mechanisms↗
CINDER08 and Scripting Environment for Accelerator Activation Problems - Paper 101
CINDER is a transmutation code, developed originally at the Bettis Atomic Power Laboratory, which predicts induced radioactivity due to either reactor or accelerator-driven radiation sources. The code and associated data libraries have been expanded and improved over the years at Los Alamos National Laboratory. CINDER08 is the latest update in the series and represents a complete rewrite of CINDER, improving it in many areas. Beginning with CINDER90, several Perl scripts were developed to treat multi-cell problems in combination with the radiation transport code MCNPX. These scripts simplified the preparation of CINDER input and made the code much more user-friendly. In concert with the development of CINDER08, these scripting tools have also been updated. The script developers are intensive users of CINDER and had significant input into the development of CINDER08. CINDER08, through use of these scripts, also works with output from MCNP5, which is useful for problems that do not require the extended particle capabilities of MCNPX and do not require particle energies outside of the tabular physics regime, and with MCNP6. A new version of CINDER, CINDER08, incorporates many improvements to the coding and the data libraries. Updated versions of scripting tools improve the code's usability and simplify the use of CINDER output as the basis of gamma source terms for subsequent calculations. (authors)
61 RADIATION PROTECTION AND DOSIMETRY↗
RAM-SCB Runs on Request at CCMC
The Ring current–Atmosphere interactions Model (RAM) with Self-Consistent magnetic (B) field (SCB) combines a large-scale kinetic model of ring current plasma with a three-dimensional (3-D) force-balanced model of the terrestrial magnetic field. RAM-SCB simulates the evolution of major ion species (H + , O + , and He + by default) and electrons as a function of azimuth, radial distance, energy, and pitch angle. Simulation outputs include the Dst index, and pressure and differential flux of the modeled species, thus providing benefit as a science code and to inform surface charging hazard within the model domain. Version 2.2 of this open-source simulation code has now—as of 8 August 2023—been made available for Runs-On-Request via the Community Coordinated Modeling Center.
79 ASTRONOMY AND ASTROPHYSICS↗
Estimating an airborne dipole source using 3-D wavefield simulations and seismic receivers on the ground
Acoustic signals can couple to the ground, providing an opportunity to use seismic stations to investigate airborne sources. The study of Bishop et al. used wavefield simulations in a fluid-solid medium to quantify the role of topography on the seismic (ground) recordings of a monopole source in the air. We build upon this study by linking wavefield forward modelling with the source estimation code MTUQ, which can accommodate point forces or moment tensors in a solid medium, as well as sources in the air (or water) if they are enabled by the forward-modelling solver. We perform a series of synthetic numerical experiments to demonstrate that a dipole airborne source can be estimated using ground-based receivers, even within the presence of realistic topography. We investigate the influence of receiver coverage, topography, and assumed source location on the estimated results. The established capabilities raise the prospects for future efforts to estimate dipole sources in 3-D models that include heterogeneity in the air and the earth in addition to topography.
58 GEOSCIENCES↗