Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pointing Performance Analysis”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Virtual Dynamic Grid Impedance and Its Impacts on Harmonics and Stability of Inverter Based Resources Plant

Nowadays, a large number of inverter-based resources (IBRs) are integrated into the grid at a single connection point as an IBR plant. In this article, a virtual dynamic grid impedance concept is introduced to evaluate the harmonics and stability for grid integration of an IBR plant containing multiple IBRs. First, a detailed theoretical study is conducted to build a foundation of the virtual grid impedance concept, which is a dynamic impedance that changes with the number of IBRs added into an IBR plant. Then, based on the new virtual dynamic grid impedance concept, frequency spectrum analysis is performed to explore harmonic impact at the interconnection point for an IBR plant with IBRs having L, LC, and LCL filters, respectively. An electromagnetic transient (EMT) simulation model of a grid-connected IBR plant is developed to explore the harmonics and stability of the IBR plant connected to the grid as well as the reliable operation of IBRs within the plant from the novel virtual equivalent dynamic grid impedance point of view. Hardware experiments are conducted to validate the EMT simulation evaluation. The results show that the number of IBRs added into an IBR plant influences the grid impedance, i.e., grid strength, and the grid impedance variation has multiple impacts on the IBR plant and IBRs within the plant depending on the grid-connected filters of the IBRs.

grid impedance↗
DOE OSTI2023

The catalog-to-cosmology framework for weak lensing and galaxy clustering for LSST

We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC), and designed for cosmology analyses using LSST data. In this paper, we present the pipeline for the so-called 3x2pt analysis -- a combination of three two-point functions that measure the auto- and cross-correlation between galaxy density and shapes. We perform the analysis both in real and harmonic space using TXPipe and other LSST-DESC tools. We validate the pipeline using Gaussian simulations and show that it accurately measures data vectors and recovers the input cosmology to the accuracy level required for the first year of LSST data under this simplified scenario. We also apply the pipeline to a realistic mock galaxy sample extracted from the CosmoDC2 simulation suite (Korytov et al. 2019). TXPipe establishes a baseline framework that can be built upon as the LSST survey proceeds. Furthermore, the pipeline is designed to be easily extended to science probes beyond the 3x2pt analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Smart Charging for Electric Ride-Hailing Vehicles using Renewables: A San Francisco Case Study

Charging large fleets of electric ride-hailing vehicles (ERVs) is a complex matter that could serve different objectives: lower carbon dioxide emissions, lower monetary expenditures, or maximize solar photovoltaics (PV) energy consumption. Currently, it is unclear how each of those objectives could impact the business and performance of a ride-hailing fleet. In order to fill this gap, this article employs a dynamic transportation model: a smart charging simulation that combines agent-based, discrete-event, and system dynamic modelling by comparing the above-mentioned objectives in separate scenarios. The results show that each scenario successfully manages to shift between 34% and 87% of all load to hours of the day when the objectives of those scenarios are met. Therefore, in comparison to the baseline, smart charging can save between 5% and 26% of monthly emissions and between 4% and 57% of monthly expenditures. The solar PV scenario, however, results in the highest savings, while ensuring profitable economics via net metering in the short- as well as long term. Finally, the sensitivity analysis points to important trade-offs between several fleet performance metrics. The article concludes by giving business and policy recommendations for maximising the economic, energy and environmental efficiency of large ERV fleets.

ADVANCED PROPULSION SYSTEMS,SOLAR ENERGY↗
DOE OSTI2022

Dark Energy Survey Year 3 results: Cosmological constraints from galaxy clustering and weak lensing

We present the first cosmology results from large-scale structure using the full 5000 deg 2 of imaging data from the Dark Energy Survey (DES) Data Release 1. We perform an analysis of large-scale structure combining three two-point correlation functions ( 3 × 2 pt ): (i) cosmic shear using 100 million source galaxies, (ii) galaxy clustering, and (iii) the cross-correlation of source galaxy shear with lens galaxy positions, galaxy–galaxy lensing. To achieve the cosmological precision enabled by these measurements has required updates to nearly every part of the analysis from DES Year 1, including the use of two independent galaxy clustering samples, modeling advances, and several novel improvements in the calibration of gravitational shear and photometric redshift inference. The analysis was performed under strict conditions to mitigate confirmation or observer bias; we describe specific changes made to the lens galaxy sample following unblinding of the results and tests of the robustness of our results to this decision. We model the data within the flat Λ CDM and w CDM cosmological models, marginalizing over 25 nuisance parameters. We find consistent cosmological results between the three two-point correlation functions; their combination yields clustering amplitude S 8 = 0.77 6 - 0.017 + 0.017 and matter density Ω m = 0.33 9 - 0.031 + 0.032 in Λ CDM , mean with 68% confidence limits; S 8 = 0.77 5 - 0.024 + 0.026 , Ω m = 0.35 2 - 0.041 + 0.035 , and dark energy equation-of-state parameter w = - 0.9 8 - 0.20 + 0.32 in w CDM . These constraints correspond to an improvement in signal-to-noise of the DES Year 3 3 × 2 pt data relative to DES Year 1 by a factor of 2.1, about 20% more than expected from the increase in observing area alone. This combination of DES data is consistent with the prediction of the model favored by the Planck 2018 cosmic microwave background (CMB) primary anisotropy data, which is quantified with a probability-to-exceed p = 0.13 –0.48. We find better agreement between DES 3 × 2 pt and Planck than in DES Y1, despite the significantly improved precision of both. When combining DES 3 × 2 pt data with available baryon acoustic oscillation, redshift-space distortion, and type Ia supernovae data, we find p = 0.34 . Combining all of these datasets with Planck CMB lensing yields joint parameter constraints of S 8 = 0.81 2 - 0.008 + 0.008 , Ω m = 0.30 6 - 0.005 + 0.004 , h = 0.68 0 - 0.003 + 0.004 , and ∑ m ν < 0.13 eV (95% C.L.) in Λ CDM ; S 8 = 0.81 2 - 0.008 + 0.008 , Ω m = 0.30 2 - 0.006 + 0.006 , h = 0.68 7 - 0.007 + 0.006 , and w = - 1.03 1 - 0.027 + 0.030 in w CDM .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Automated segmentation and analysis of point clouds of pier foundations using Pier Inspection and Evaluation Report (PIER)

Pier foundations are commonly used in locations with unstable soil or where other types of foundations are unsuitable or cost prohibitive. A pier foundation consists of vertical columns to support the structure and elevate it above the ground. Common materials for pier foundations include masonry, concrete, timber, and steel. The methods for accurate placement of pier foundations have remained relatively unchanged for decades. For simple installations, construction chalk lines are used to layout the locations of piers to ensure accurate placement and elevation. For more complex installations, surveying instruments operated by trained professionals are employed to accurately locate piers and assess correct elevation before construction. After installation, another survey may need to be performed to assess the quality of the as-built foundation. However, with the advent of terrestrial laser scanners (TLS), the means now exist for contractors to conduct their own assessments of as-built foundations. The major barrier preventing contractors from performing their own assessments of as-built foundation quality is the segmentation and analysis of point cloud data, a skill that often requires a trained user. The objective of this research is to develop a software tool (PIER: Pier Inspection and Evaluation Report) to enable automated segmentation and analysis of point clouds of pier foundations. In this paper, the automated segmentation and analysis algorithms are detailed. A mockup lay out of pier foundations was built using concrete masonry units, and the algorithms were tested to evaluate performance. Limitations of the current algorithms and future research direction are discussed.

Turki, Amine [ORNL]↗
DOE OSTI2021

Experimental Measurements of the Wake of a Sphere at Subcritical Reynolds Numbers

This work experimentally investigated the flow phenomena and vortex structures in the wake of a sphere located in a water loop at Reynolds numbers of Re = 850, 1,250, and 1,700. Velocity fields in the wake region were obtained by applying the time-resolved stereoscopic particle image velocimetry (TR-SPIV) technique. From the acquired TR-SPIV velocity vector fields, the statistical values of mean and fluctuating velocities were computed. Spectral analysis, two-point velocity–velocity cross-correlation, proper orthogonal decomposition (POD) and vortex identification analyses were also performed. The velocity fields show a recirculation region that decreases in length with an increase of Reynolds numbers. The power spectra from the spectral analysis had peaks corresponding to a Strouhal number of St = 0.2, which is a value commonly found in the literature studies of flow over a sphere. The two-point cross-correlation analysis revealed elliptical structures in the wake, with estimated integral length scales ranging between 12% and 63% of the sphere diameter. The POD analysis revealed the statistically dominant flow structures that captured the most flow kinetic energy. It is seen that the flow kinetic energy captured in the smaller scale flow structures increased as Reynolds number increased. The POD modes contained smaller structure as the Reynolds number increased and as mode order increased. In addition, spectral analysis performed on the POD temporal coefficients revealed peaks corresponding to St = 0.2, similar to the spectral analysis on the fluctuating velocity. The ability of POD to produce low-order reconstructions of the flow was also utilized to facilitate vortex identification analysis, which identified average vortex sizes of 0.41D for Re1, 0.33D for Re2, and 0.32D for Re3.

Engineering↗
DOE OSTI2022

Accelerating cosmological inference with Gaussian processes and neural networks – an application to LSST Y1 weak lensing and galaxy clustering

ABSTRACT Studying the impact of systematic effects, optimizing survey strategies, assessing tensions between different probes and exploring synergies of different data sets require a large number of simulated likelihood analyses, each of which cost thousands of CPU hours. In this paper, we present a method to accelerate cosmological inference using emulators based on Gaussian process regression and neural networks. We iteratively acquire training samples in regions of high posterior probability which enables accurate emulation of data vectors even in high dimensional parameter spaces. We showcase the performance of our emulator with a simulated 3×2 point analysis of LSST-Y1 with realistic theoretical and systematics modelling. We show that our emulator leads to high-fidelity posterior contours, with an order of magnitude speed-up. Most importantly, the trained emulator can be re-used for extremely fast impact and optimization studies. We demonstrate this feature by studying baryonic physics effects in LSST-Y1 3×2 point analyses where each one of our MCMC runs takes approximately 5 min. This technique enables future cosmological analyses to map out the science return as a function of analysis choices and survey strategy.

Astronomy & Astrophysics↗
DOE OSTI2022

Suppression of ferromagnetism governed by a critical lattice parameter in CeTiGe 3 with hydrostatic pressure or V substitution

We combine structural and magnetic measurements to compare the different magnetic phase diagrams between the pressure and substitution studies in CeTiGe 3 . We report on the structural, magnetic, and electrical transport properties of single crystals of CeTi 1–x V x Ge 3 (x = 0, 0.1, 0.2, 0.3, 0.4, 0.9, and 1), and of polycrystalline samples (x = 0.5, 0.6, 0.7, 0.8), as well as structural properties of CeTiGe 3 under pressure up to 9 GPa. The ferromagnetic ordering in CeTiGe 3 is suppressed with V doping in CeTi 1–x V x Ge 3 , and suggests a possible ferromagnetic quantum critical point near x = 0.45. We perform a detailed crystalline electric field (CEF) analysis, and the magnetic susceptibility data in pure CeTiGe 3 and CeVGe 3 can be well explained by the CEF model. The proposed CEF energy levels suggest that there is a gradual change of the ground state from |±5/2 in CeTiGe 3 to |±1/2 in CeVGe 3 , and a suppression of CEF splitting energies near the quantum critical region. When hydrostatic pressure is used instead of chemical substitution, the quantum critical point is avoided by the appearance of magnetic phases above around 4.1 GPa. In the substitution study, the ferromagnetic and antiferromagnetic regions are well separated, whereas they touch in the pressure study. We observe a different trend in the temperature dependence of the resistivity maximum in both studies, suggesting that the CEF splitting energy is suppressed by V substitution but enhanced by pressure. Here, we also observe different responses in lattice constants between the two studies, highlighting the fact that substitution effects cannot be reduced to chemical pressure effects only. Nevertheless, when the magnetic phase diagrams of both hydrostatic pressure and substitution are compared, we find a common critical lattice constant c = 5.78 Å where the ferromagnetic ordering is suppressed in both studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The Gross-Neveu-Yukawa archipelago

We perform a bootstrap analysis of a mixed system of four-point functions of bosonic and fermionic operators in parity-preserving 3d CFTs with O(N) global symmetry. Our results provide rigorous bounds on the scaling dimensions of the O(N)-symmetric Gross-Neveu-Yukawa (GNY) fixed points, constraining these theories to live in isolated islands in the space of CFT data. We focus on the cases N = 1, 2, 4, 8, which have applications to phase transitions in condensed matter systems, and compare our bounds to previous analytical and numerical results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Data transfers for full core heterogeneous reactor high- fidelity multiphysics studies

Multiphysics simulations for nuclear reactor analysis are usually performed by resorting to operator splitting and fixed point iterations between single-physics solvers. This enables the separate solution of each physics, such as neutronics, fuel performance, and thermal hydraulics, on meshes tailored to the requirements of the respective numerical discretizations of the equations. As the equations are coupled, several fields must be transferred between single-physics solves. Projecting fields between meshes while preserving order of accuracy, conservation properties, and mapping non-overlapping geometries is a complex endeavor. This conference paper will present the transfers as implemented in MOOSE, which can handle arbitrary meshes, arbitrary mappings, conservation of integral quantities, and are made to scale with distributed simulations on both ends of the transfers. Their adequacy for advanced nuclear reactor multiphysics coupling is shown through examples and numerical studies.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Dark Energy Survey Year 6 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing

We present cosmology results combining galaxy clustering and weak gravitational lensing measured in the full six years (Y6) of observations by the Dark Energy Survey (DES) covering $\sim$5000 deg$^2$. We perform a large-scale structure analysis using three two-point correlation functions (3$\times$2pt): (i) cosmic shear from 140 million source galaxy shapes, (ii) galaxy clustering of 9 million lens galaxy positions, and (iii) galaxy-galaxy lensing from their cross-correlation. We model the data in flat $Λ$CDM and $w$CDM cosmologies. The combined analysis yields $S_8\equiv σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.789^{+0.012}_{-0.012}$ and matter density $Ω_{\rm m} = 0.333^{+0.023}_{-0.028}$ in $Λ$CDM (68% CL), where $σ_8$ is the clustering amplitude. These constraints show a (full-space) parameter difference of 1.8$σ$ from a combination of cosmic microwave background (CMB) primary anisotropy datasets from Planck 2018, ACT-DR6, and SPT-3G DR1. Projected only into $S_8$ the difference is $2.6σ$. In $w$CDM the Y6 3$\times$2pt results yield $S_8 = 0.782^{+0.021}_{-0.020}$, $Ω_{\rm m} = 0.325^{+0.032}_{-0.035}$, and dark energy equation-of-state parameter $w = -1.12^{+0.26}_{-0.20}$. For the first time, we combine all DES dark-energy probes: 3$\times$2pt, SNe Ia, BAO and Clusters. In $Λ$CDM this combination yields a $2.8σ$ parameter difference from the CMB. When combining DES 3$\times$2pt with other low-redshift datasets (DESI DR2 BAO, DES SNe Ia, SPT clusters), we find a 2.3$σ$ parameter difference with CMB. A joint fit of Y6 3$\times$2pt, CMB, and those low-redshift datasets produces the tightest $Λ$CDM constraints to date: $S_8 = 0.806^{+0.006}_{-0.007}$, $Ω_{\rm m} = 0.302^{+0.003}_{-0.003}$, $h = 0.683^{+0.003}_{-0.002}$, and $\sum m_ν< 0.14$ eV (95% CL). In $w$CDM, this combination yields $w = -0.981^{+0.021}_{-0.022}$.

Abbott, T. M.C. [Cerro-Tololo InterAmerican Obs.]↗
DOE OSTI2025

Implications of point defect accumulation on UO 2 thermal conductivity and fission gas release under accelerated fuel irradiation

Evaluation of thermal properties is a crucial factor for nuclear fuel performance. During reactor operation, the accumulation of fission products and irradiation-induced lattice defects are responsible for degradation in thermal conductivity. Consequently, it affects fuel temperature and fission gas release (FGR) among other Multiphysics processes important for economics and safety analysis. We analyze the implications of point defects (PD) accumulation described using a rate theory (RT) Model on lattice thermal conductivity of UO 2 . Here, we demonstrate that fission rate-dependent point defect concentrations have the largest impact on in-pile thermal conductivity in the periphery of light water reactor fuels below a temperature threshold governed by the migration barrier of defects. Our analysis provides a mechanistic description of this phenomena which current fuel performance codes treat empirically. The reduction of thermal conductivity in the low -temperature rim region acts as additional thermal resistance and leads to a temperature notably larger than suggested by Lucuta thermal conductivity correlation. These effects are anticipated to have notable impacts when fuels are exposed to accelerated radiation. The impact of such point defect-informed treatment of thermal conductivity on fuel performance is evaluated by a detailed analysis of fission gas behavior and its release. We consider several models capturing different stages of fission gas bubble evolution and fission gas release (FGR). Finally, a new fission rate-dependent correction to the Lucuta correlation is proposed. The results show a significant reduction in thermal conductivity at the fuels’ periphery and an increase in fuel centerline temperature specifically at low burnups. Ultimately a modified LC shows a higher FGR compared to the original LC, while the acceleration process results in a reduction in overall FGR.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Charged point defect benchmark of Hybrid and GGA-PBE

Data used by the publication "High-throughput calculations of charged point defect properties with semi-local density functional theory - performance benchmarks for materials screening applications." This work presented an in-depth benchmark analysis of automated, semi-local point defect calculations with a-posteriori corrections, compared to 245 “gold standard” hybrid calculations previously published. We considered three different a-posteriori correction sets for semi-local calculations, implemented in a fully automated workflow, and consider the qualitative and quantitative differences for four different categories of defect information: thermodynamic transition levels, formation energies, fermi levels, and dopability limits. We highlighted the type of qualitative information about point defect properties that can be extracted from high-throughput calculations based on semi-local DFT methods, while also demonstrating the limits of quantitative accuracy that can be achieved by these approaches.

Broberg, Danny↗
DOE OSTI2022

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Whole program Adaptive Error Detection and Mitigation

Our project investigated methods to detect soft errors in computing systems. Such detection schemes must have low overheads and also generate no false positives. Our solutions focused on control-space error detection based on changing addressing schemes so that faults tend to cascade – rather than being masked. Our solutions also addressed data-space error detection by synthesizing error detectors based on rigorous floating-point error analysis.

42 ENGINEERING↗
DOE OSTI2020

Stability Analysis of Inline ZFP Compression for Floating-Point Data in Iterative Methods

Currently, the dominating constraint in many high performance computing applications is data capacity and bandwidth, in both internode communications and even moreso in intranode data motion. A new approach to address this limitation is to make use of data compression in the form of a compressed data array. Storing data in a compressed data array and converting to standard IEEE-754 types as needed during a computation can reduce the pressure on bandwidth and storage. However, repeated conversions (lossy compression and decompression) introduce additional approximation errors, which need to be shown to not significantly affect the simulation results. Here, we extend recent work that analyzed the error of a single use of compression and decompression of the ZFP compressed data array representation to the case of time-stepping and iterative schemes, where an advancement operator is repeatedly applied in addition to the conversions. We show that the accumulated error for iterative methods involving fixed-point and time evolving iterations is bounded under standard constraints. An upper bound is established on the number of additional iterations required for the convergence of stationary fixed-point iterations. An additional analysis of traditional forward and backward error of stationary iterative methods using ZFP compressed arrays is also presented. The results of several 1D, 2D, and 3D test problems are provided to demonstrate the correctness of the theoretical bounds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

RELAP5-3D validation studies based on the High Temperature Test facility

In the spring and summer of 2019, experiments were conducted at the High Temperature Test Facility (HTTF) that form the basis of an upcoming high-temperature gas-cooled reactor (HTGR) thermal hydraulics (T/H) benchmark. HTTF is an integral effects test facility for HTGR T/H modeling validation. This paper presents RELAP5-3D models of two of those experiments: PG-27, a pressurized conduction cooldown (PCC); and PG-29, a depressurized conduction cooldown (DCC). These models used the RELAP5-3D model of HTTF originally developed by Paul Bayless as a starting point. The sensitivity analysis and uncertainty quantification code, RAVEN was used to perform calibration studies for the steady-state portion of PG-27. Here we developed four PG-27 calibrations based on steady-state conditions. These calibrations all used an effective thermal conductivity equal to 36 % of the measured thermal conductivity, but they differed with respect to the frictional pressure drops and radial conduction models. These models all captured the trends in steady-state temperature distributions and transient temperature behavior well. All four calibrations show room for improvement in predicting the transient temperature rise. The smallest error in temperature rise during the transient was a 21 % underprediction, and the largest was a 48 % underprediction. The errors in transient temperature rise are largely a result of a mismatch in power density between the RELAP5-3D model and the experiment due to the location of active heater rods along the boundary between heat structures in the model. The best of these calibrations was applied to PG-29 to model the DCC. Once again, temperatures during the transient were underpredicted but trends in temperature were captured. The RELAP5-3D model captured trends in the data but could not reproduce measured temperatures exactly. This result is not attributed to deficiencies in the experimental data or to RELAP5–3D itself. Rather, this result likely arises due to the some of the assumptions and decisions made when the RELAP5-3D model was first developed, prior to the execution of HTTF experiments. An agreement in prediction of temperature trends but challenges reproducing HTTF temperatures within measurement uncertainty is consistent with previous analyses of HTTF in the literature. Future RELAP5-3D validation activities centered around HTTF may be able to provide greater insight into the code’s capabilities for HTGR modeling with a more finely nodalized model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

PyCrystalField : software for calculation, analysis and fitting of crystal electric field Hamiltonians

PyCrystalField is a Python software package for calculating single-ion crystal electric field (CEF) Hamiltonians. This software can calculate a CEF Hamiltonian ab initio from a point charge model for any transition or rare earth ion in either the J basis or the LS basis, perform symmetry analysis to identify nonzero CEF parameters, calculate the energy spectrum and observables such as neutron spectrum and magnetization, and fit CEF Hamiltonians to any experimental data. Here, the theory, implementation and examples of its use are discussed.

36 MATERIALS SCIENCE↗