Benchmark Analysis for Space Radiation Transport Codes in Thick Shielding for GCR Environments
No abstract available
Slaba, Tony C.↗
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.
MM_Mix: An Equation of State Generation Code for Dense Mixtures
Explore the source record for details and available documents.
Physics - Condensed matter physics↗
Site-level Foliar C, N, delta13C data from samples collected during field survey associated with NEON AOP survey, East River, CO 2018
Carbon (C) and nitrogen (N) weight percent concentrations were obtained from the bulk foliar samples collected across East River, Washington Gulch, Slate River, and Coal Creek watersheds in Gunnison, Colorado during the summer of 2018. These samples were collected from sun-lit leaves of meadow, shrub, and tree vegetation types within 72 hours of the National Ecological Observatory Network Airborne Observation Platform (NEON AOP) hyperspectral survey. These samples were dried and then run on an elemental analyzer coupled to a mass spectrometer for C and N concentrations and delta13C. The intention of these data were to provide ground truth data to develop foliar nutrient models across the imaging spectroscopy domain. Additional associated site metadata and vegetation fractional cover information can be found in the associated data releases (DOI: 10.15485/1618130). Here we provide site level C, N, and delta 13C data in the CN_Results_Foliar.csv. All of these samples were also registered with IGSN codes, and the details from that registration with metadata for the samples is included in foliar_sample_igsn_metadata.csv.
54 ENVIRONMENTAL SCIENCES↗
Benchmarking Monte Carlo codes for the modelling of low-energy neutron production target reactions
The increasing adoption of accelerator-based neutron sources (ABNS) for applications including neutron capture therapy (NCT) research has highlighted the need for accurate simulation tools. Precise modelling of the neutron production target is crucial to ensure that simulated predictions of neutron beam characteristics used for subsequent beam shaping assembly design are reliable. This work presents a comprehensive benchmarking of four widely-used Monte Carlo codes - Geant4, PHITS, FLUKA (CERN), and MCNP - for modelling low-energy neutron production target reactions. Using their recommended physics models and cross-section libraries, we evaluate each code’s performance in simulating four beam-target reactions: 7 Li(p,n) 7 Be, 9 Be(p,n) 9 B, 9 Be(d,n) 10 B, and C(d,n)N. Predictions of neutron yield, angular distributions, and energy spectra are compared against available thick target experimental data. Results show varying levels of agreement between the codes depending on the reaction type, energy range, and beam characteristics. Geant4, MCNP and PHITS are the overall best performing codes for the simulation of total neutron yield and yield in the forward direction across most reactions. Across energies where experimental benchmarks exist, inter-code discrepancies in total and forward-directed yield are typically 10 to 30%, with larger deviations at near-threshold incident ion energies. PHITS provides the best overall reproduction of experimental spectra, particularly for the 9 Be(p,n) 9 B reaction. Additionally, PHITS demonstrates superior computational performance for most reactions. These findings provide valuable guidance for ABNS design, highlighting the strengths and limitations of each code for the simulation of low-energy neutron production reactions.
43 PARTICLE ACCELERATORS↗
Venus Global Reference Atmospheric Model (Venus-GRAM) Upgrades
Introduction: The Venus Global Reference Atmospheric Model (Venus-GRAM) is one of the most widely used engineering models of Venus’ atmosphere. The Venus-GRAM upgrades are being developed by NASA Marshall Space Flight Center (MSFC) and NASA Langley Research Center (LaRC). This presentation will provide details regarding the upgrades that have been made to Venus-GRAM and the ongoing objectives, tasks, and milestones related to the GRAM upgrades funded by the NASA Science Mission Directorate (SMD). Venus-GRAM: Venus-GRAM is an engineering-oriented atmospheric model that estimates mean values and statistical variations of the atmospheric properties of Venus. Venus-GRAM provides mean values and variability for any point in the atmosphere as well as seasonal, geographic, and altitude variations. Venus-GRAM outputs include atmospheric density, temperature, pressure, winds, and chemical composition along a user-defined path. It is extensively used by the engineering community because of its ability to create realistic dispersions. GRAMs have been integrated into high fidelity flight dynamic simulations of launch, entry, descent and landing (EDL), aerobraking and aerocapture. GRAMs are currently available for Earth, Mars, Venus, Neptune, Titan, and Uranus. The lower atmosphere model in Venus-GRAM (up to 250 km) is based on the Venus International Reference Atmosphere (VIRA) [1]. The Venus-GRAM thermosphere (250 to 1000 km) is based on a MSFC-developed model [2] which assumes an isothermal temperature profile initialized using VIRA conditions at 250 km [3]. The VIRA version included in Venus-GRAM includes Pioneer Venus Orbiter and Probe data as well as Venera probe data, but it does not include a solid planet model or a high-resolution gravity model [4]. Venus-GRAM Upgrade Status: Code Modernization. Venus-GRAM has been rearchitected from Fortran to a common object-oriented C++ framework called the GRAM Suite. This new architecture creates a common GRAM library of data models and utilities. The first C++ release of the rearchitected Venus-GRAM is a straight conversion from the latest Fortran version. Model Upgrades. The focus of the model upgrade task is to improve the atmosphere models in the existing GRAMs and to establish a foundation for developing GRAMs for additional destinations. The GRAM ephemeris has been upgraded to the NASA Navigation and Ancillary Information Facility (NAIF) SPICE toolkit (version N0066). The calculation of the speed of sound has also been improved in the GRAMs. In FY20, the GRAM project established contracts to improve the model data within Venus-GRAM. Hampton University is developing an empirical global model for Venus. The University of Wisconsin is reanalyzing the Venus Express radio occultation observations and analyzing the Akatsuki thermal imaging data. Upgraded Venus-GRAM Release. GRAM Suite Version 1.3 will be released in September 2021 and will contain the rearchitected Venus-GRAM, including the common GRAM framework and planet–specific code. A User Guide and Programmer’s Manual are released with all GRAMs. Conclusions: GRAMs are frequently used toolsets and vital in assessing effects of atmospheres on interplanetary spacecraft during the program life cycle process. Releases of the GRAM Suite, upgrades of the existing planetary GRAMs, and development of new planetary GRAMs are ongoing. Venus-GRAM atmosphere model upgrades will be included in the next phase of GRAM tasks. References: [1] Kliore, A. J. et al. (1985) ASR, 5, 11, 1-304. [2] Justh, H. L. et al. (2006) AIAA/AAS Astrodynamics Specialist Conference & Exhibit, Abstract AIAA-2006-6394. [3] Guide to Reference and Standard Atmosphere Models, BSR/AIAA G-003-2010. [4] Limaye, S. S. (2012), LPSC VEXAG Townhall Meeting. Acknowledgments: The authors gratefully acknowledge support from the NASA SMD.
atmospheric models↗
A Navier-Stokes solver for turbomachinery applications
A computer code for solving the Reynolds-averaged full Navier-Stokes equations has been developed and applied using H- and C-type grids. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. The integration in time is based on an explicit four-stage Runge-Kutta scheme. Local time stepping, variable coefficient implicit residual smoothing, and a full multigrid method have been implemented to accelerate steady-state calculations. A grid independence analysis is presented for a transonic rotor blade. Comparisons with experimental data show that the code is an accurate viscous solver and can give very good blade-to-blade predictions for engineering applications.
Arnone, A.↗
Implicit LES of Turbulent, Separated Flow: Wall-Mounted Hump Configuration
Direct simulations (ILES) of turbulent, separated flow over the wall-mounted hump configuration is conducted to investigate the physics of separated flows. A chord-based Reynolds number of Re(sub c) = 47,500 is set up, with a turbulent in flow of Re(sub theta) = 1,400 (theta/c = 3%). FDL3DI, a code that solves the compressible Navier-Stokes equations using high- order compact-difference scheme and filter, with the standard recycling/rescaling method of turbulence generation, is used. Two different configurations of the upper-wall are analyzed, and results are compared with both a higher Re(sub c) (= 936,000, Re(sub theta) = 7,200, theta/c = 0.77%) experiment for major flow features, and RANS (k-omega SST) results. A lower Rec allows for DNS-like mesh resolution, and an adequately wide span. Both ILES and RANS show delayed reattachment compared to experiment, and significantly higher skin friction in the forebody of the hump, as expected. The upper-wall shape influences the C(sub p) distribution only. Results from this study are being used to setup higher Rec (lower theta/c) ILES.
Implicit LES↗
The CubeSat Communication Platform (CCP) – Mission Overview and ConOps
This paper presents the CubeSat Communications Platform (CCP) mission concept, architecture, and development. The CCP is a CubeSat demonstration mission, being developed at the University of Alaska Fairbanks to improve satellite communication capabilities. The CCP payloads include a Software Defined Radio (SDR) with Variable Coded Modulation (VCM) protocols, and an S-band phased array antenna. The mission will test the performance of the VCM protocols versus conventional fixed modulation and coding schemes, relative to the percentage of the Shannon Channel Capacity that each scheme achieves over entire satellite passes. The mission will also test the performance of the phased array antenna, with regard to beamforming and interrogator scanning accuracy. The CCP Mission is collaborating with NASA's Near Space Network (NSN) to demonstrate DVB-S2 VCM and achievable maximum data rate in NASA S-band 5 MHz channel. NSN currently supports missions that communicate with fixed channel codes, modulations, and symbol rates, resulting in a constant data rate that does not adapt to the dynamic link margin. VCM adapts to the dynamics of the link to increase information throughput by changing modulation and coding when the signal-to-noise ratio (SNR) is high. The CCP will be the first mission to demonstrate VCM with NSN ground stations.
phased antenna array↗
FlowPM: Distributed TensorFlow implementation of the FastPM cosmological N-body solver
Here, we present FlowPM, a Particle-Mesh (PM) cosmological N-body code implemented in Mesh-TensorFlow for GPU-accelerated, distributed, and differentiable simulations. We implement and validate the accuracy of a novel multi-grid scheme based on multiresolution pyramids to compute large-scale forces efficiently on distributed platforms. We explore the scaling of the simulation on large-scale supercomputers and compare it with corresponding Python based PM code, finding on an average 10x speed-up in terms of wallclock time. We also demonstrate how this novel tool can be used for efficiently solving large scale cosmological inference problems, in particular reconstruction of cosmological fields in a forward model Bayesian framework with hybrid PM and neural network forward model. We provide skeleton code for these examples and the entire code is publicly available at https://github.com/modichirag/flowpm[Formula presented].
79 ASTRONOMY AND ASTROPHYSICS↗
Upgrade of hardware controls for the STAR experiment at RHIC
The STAR experiment has been delivering significant physics results for more than 20 years. Stable operation of the experiment was achieved by using a robust controls system based on the Experimental Physics and Industrial Control System (EPICS). Now an object-oriented approach with Python libraries, adapted for EPICS software, is going to replace the procedural-based EPICS C libraries previously used at STAR. Advantages of the new approach include stability of operation, code reduction and straightforward project documentation. The first two sections of this paper introduce the STAR experiment, give an overview of the EPICS architecture, and present the use of Python for controls software. Therefore, specific examples, as well as upgrades of user interfaces, are outlined in the following sections.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Achieving performance portability in Gaussian basis set density functional theory on accelerator based architectures in NWChemEx
The numerical integration of the exchange–correlation (XC) potential is one of the primary computational bottlenecks in Gaussian basis set Kohn–Sham density functional theory (KS-DFT). To achieve optimal performance and accuracy, care must be taken in this numerical integration to preserve local sparsity as to allow for near linear weak scaling with system size. This leads to an integration scheme with several performance critical kernels which must be hand optimized for each architecture of interest. As the set of available accelerator hardware goes more diverse, a key challenge for developers of KS-DFT software is to maintain performance portability across a wide range of computational architectures. In this article, we examine a modular software design pattern which decouples the implementation details of performance critical kernels from the expression of high-level algorithmic workflows in a device-agnostic language such as C++; thus allowing for developers to target existing and emerging accelerator hardware within a single code base. We consider the efficacy of such a design pattern in the numerical integration of the XC potential by demonstrating its ability to achieve performance portability across a set of accelerator architectures which are representative of those on current and future U.S. Department of Energy Leadership Computing Facilities.
97 MATHEMATICS AND COMPUTING↗
Electron temperature gradient instability and transport analysis in NSTX and NSTX-U plasmas
Extensive linear and nonlinear simulations to study electron temperature gradient (ETG) stability and thermal transport in National Spherical Torus Experiment (NSTX) and NSTX-U plasmas were performed using the fully electromagnetic gyrokinetic code CGYRO. Linear simulations were performed to determine ETG thresholds in different discharges, showing that ETG modes in spherical tokamaks can present different scalings compared to conventional aspect-ratio tokamaks. Nonlinear gyrokinetic simulations were conducted for selected cases to calculate electron thermal transport and compare to experimental values. Results are also compared with those of ETG modes in the multi-mode model and the Trapped-Gyro-Landau-Fluid reduced model codes, to better understand their applicability in spherical tokamaks.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Extended Physics-Informed Neural Networks (XPINNs): A Generalized Space-Time Domain Decomposition Based Deep Learning Framework for Nonlinear Partial Differential Equations
Here we propose a generalized space-time domain decomposition approach for the physics-informed neural networks (PINNs) to solve nonlinear partial differential equations (PDEs) on arbitrary complex-geometry domains. The proposed framework, named eXtended PINNs ( X P I N N s ), further pushes the boundaries of both PINNs as well as conservative PINNs (cPINNs), which is a recently proposed domain decomposition approach in the PINN framework tailored to conservation laws. Compared to PINN, the XPINN method has large representation and parallelization capacity due to the inherent property of deployment of multiple neural networks in the smaller subdomains. Unlike cPINN, XPINN can be extended to any type of PDEs. Moreover, the domain can be decomposed in any arbitrary way (in space and time), which is not possible in cPINN. Thus, XPINN offers both space and time parallelization, thereby reducing the training cost more effectively. In each subdomain, a separate neural network is employed with optimally selected hyperparameters, e.g., depth/width of the network, number and location of residual points, activation function, optimization method, etc. A deep network can be employed in a subdomain with complex solution, whereas a shallow neural network can be used in a subdomain with relatively simple and smooth solutions. We demonstrate the versatility of XPINN by solving both forward and inverse PDE problems, ranging from one-dimensional to three-dimensional problems, from time-dependent to time-independent problems, and from continuous to discontinuous problems, which clearly shows that the XPINN method is promising in many practical problems. The proposed XPINN method is the generalization of PINN and cPINN methods, both in terms of applicability as well as domain decomposition approach, which efficiently lends itself to parallelized computation. The XPINN code is available on h t t p s : / / g i t h u b . c o m / A m e y a J a g t a p / X P I N N s .
97 MATHEMATICS AND COMPUTING↗
Noise variance estimator for high-rate coded telemetry
Noise variance extreme value estimator for biorthogonally block coded high rate telemetry system for Mariner missions
Posner, E. C.↗
An experimental documentation of a hypersonic shock-wave turbulent boundary layer interaction flow: With and without separation
An experiment, thoroughly documenting the flow field resulting from the interaction of a shock wave with a nonadiabatic hypersonic turbulent boundary layer, is described. Detailed mean flow and surface data are presented for two shock strengths resulting in attached and separated flows, respectively. The surface measurements include continuous pressure, shear and heat-flux distributions upstream, in, and downstream of the interaction regions. At closely spaced intervals along the surface, boundary-layer profiles of static and pitot pressure and total temperature were obtained from which velocity, density and static temperature profiles were derived. The data are presented in both graphical and tabular form. These data are of sufficient detail to validate advanced computer codes and their associated turbulence models.
Kussoy, M. I.↗
Study of the modifications needed for effective operation NASTRAN on IBM virtual storage computers
The necessary modifications were determined to make NASTRAN operational under virtual storage operating systems (VS1 and VS2). Suggested changes are presented which will make NASTRAN operate more efficiently under these systems. Estimates of the cost and time involved in design, coding, and implementation of all suggested modifications are included.
Mccormick, C. W.↗
Optimal source codes for geometrically distributed integer alphabets
An approach is shown for using the Huffman algorithm indirectly to prove the optimality of a code for an infinite alphabet if an estimate concerning the nature of the code can be made. Attention is given to nonnegative integers with a geometric probability assignment. The particular distribution considered arises in run-length coding and in encoding protocol information in data networks. Questions of redundancy of the optimal code are also investigated.
Gallager, R. G.↗
Syndrome-source-coding and its universal generalization
A method of using error-correcting codes to obtain data compression, called syndrome-source-coding, is described in which the source sequence is treated as an error pattern whose syndrome forms the compressed data. It is shown that syndrome-source-coding can achieve arbitrarily small distortion with the number of compressed digits per source digit arbitrarily close to the entropy of a binary memoryless source. A 'universal' generalization of syndrome-source-coding is formulated which provides robustly effective distortionless coding of source ensembles. Two examples are given, comparing the performance of noiseless universal syndrome-source-coding to (1) run-length coding and (2) Lynch-Davisson-Schalkwijk-Cover universal coding for an ensemble of binary memoryless sources.
Ancheta, T. C., Jr.↗