SEARCH · Search NASA

Results for “Random number generation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Optimal calibration of gates in trapped-ion quantum computers

Abstract To harness the power of quantum computing, it is essential that a quantum computer provide maximal possible fidelity for a quantum circuit. To this end, much work has been done in the context of qubit routing or embedding, i.e., mapping circuit qubits to physical qubits based on gate performance metrics to optimize the fidelity of execution. Here, we take an alternative approach that leverages a unique capability of a trapped-ion quantum computer, i.e., the all-to-all qubit connectivity. We develop a method to determine a fixed number (budget) of quantum gates that, when calibrated, will maximize the fidelity of a batch of input quantum programs. This dynamic allocation of calibration resources on randomly accessible gates, determined using our heuristics, increases, for a wide range of calibration budget, the average fidelity from 70% or lower to 90% or higher for a typical batch of jobs on an 11-qubit device, in which the fidelity of calibrated and uncalibrated gates are taken to be 99% and 90%, respectively. Our heuristics are scalable, more than 2.5 orders of magnitude faster than a randomized method for synthetic benchmark circuits generated based on real-world use cases.

Maksymov, Andrii (ORCID:0000000282442851)↗

DOE OSTI2019

Effect of Mach number and volume fraction in air-shock interacting with a bed of randomly distributed spherical particles

Shock-particle interaction is a fundamental problem in many engineering applications, with the dynamics being heavily influenced by the incident-shock Mach number and the particle volume fraction. In this paper we present fully resolved inviscid simulations of an incident-shock wave traveling through a bed of randomly distributed spherical particles. We vary the strength of the incident shock along with the particle volume fraction in order to study the complex wave interaction during shock-particle interaction. In this study we are interested in the early-time behavior during which the particles do not move and hence in our simulations all the particles are fixed in space. We compute the streamwise average of flow field quantities to generate the x - t contour plots to study the unsteady oscillations inside the particle bed. We observe that the transmitted shock slows down under certain conditions and it is partly due to tortuosity and partly due to weakening caused by energy dissipation. We also present the force histories of the streamwise drag and lift forces for all the particles. The random distribution of particles leads to high variability in the drag force experienced by the particles. We compute the mean peak drag force as a function of the streamwise location to study the mean behavior of the transmitted shock. Based on our findings, we propose simple modifications to improve the current point-particle models used in Euler-Lagrange simulations of shock interacting with a bed of particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1979

Equilibrium velocities of a planetesimal population

The random velocities of planetesimal populations specified by maximum and minimum masses and a power-law number density are calculated iteratively. The calculation is based on two different physical models involving ratios of rates: (1) excitation of kinetic energy by gravitational perturbation and elastic collision equal to damping of kinetic energy by inelastic collisions; and (2) excitation of kinetic energy a ratio b (about 3 usually) to doubling of mass. Both physical models are developed approximately and then more precisely. The approximate model (2) agrees with results of Safronov (1972). Both precise models are applied to populations approximating those generated by the calculations of Greenberg et al. (1978), Cox and Lewis (1979), and Wetherill (1979).

Kaula, W. M.↗

NASA NTRS2002

Measurements of Turbulence Attenuation by a Dilute Dispersion of Solid Particles in Homogeneous Isotropic Turbulence

This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact the addition of gravity as a variable parameter may help us to better understand the physics of turbulence attenuation. The experiments are conducted in a turbulence chamber capable of producing stationary or decaying isotropic turbulence with nearly zero mean flow and Taylor microscale Reynolds numbers up to nearly 500. The chamber is a 410 mm cubic box with the corners cut off to make it approximately spherical. Synthetic jet turbulence generators are mounted in each of the eight corners of the box. Each generator consists of a loudspeaker forcing a plenum and producing a pulsed jet through a 20 mm diameter orifice. These synthetic jets are directed into ejector tubes pointing towards the chamber center. The ejector tubes increase the jet mass flow and decrease the velocity. The jets then pass through a turbulence grid. Each of the eight loudspeakers is forced with a random phase and frequency. The resulting turbulence is highly Isotropic and matches typical behavior of grid turbulence. Measurements of both phases are acquired using particle image velocimetry (PIV). The gas is seeded with approximately 1 micron diameter seeding particles while the solid phase is typically 150 micron diameter spherical glass particles. A double-pulsed YAG laser and a Kodak ES-1.0 10-bit PIV camera provide the PIV images. Custom software is used to separate the images into individual images containing either gas-phase tracers or large particles. Modern high-resolution PIV algorithms are then used to calculate the velocity field. A large set of image pairs are acquired for each case, then the results are averaged both spatially and over the ensemble of acquired images. The entire apparatus is mounted in two racks which are carried aboard NASA's KC-135 Flying Microgravity Laboratory. The rack containing the turbulence chamber, the laser head, and the camera floats freely in the airplane cabin (constrained by competent NASA personnel) to minimize g-jitter.

Eaton, John↗

DOE OSTI2011

Using Psychokinesis to Explore the Nature of Quantum Randomness

In retrocausation different causal events can produce different successor events, yet a successor event reflecting a particular cause occurs before the causal event does. It is sometimes proposed that the successor event is determined by propagation of the causal effect backwards in time via the dynamical equations governing the events. However, because dynamical equations are time reversible, the evolution of the system is not subject to change. Therefore, the backward propagation hypothesis implies that what may have seemed to be an arbitrary selection of a causal factor was in reality predetermined.Yet quantum randomness can be used to determine the causal factor, and a quantum random event is ordinarily thought of as being arbitrarily generated. So we must ask, when quantum random events occur, are they arbitrary (subject to their probabilistic constraints) or are they predetermined?Because psychokinesis (PK) can act on quantum random events, it can be used as a probe to explore questions such as the above. It is found that if quantum random events are predetermined (aside from the action of PK), certain types of experimental design can show enhanced PK through the use of precognition. Actual experiments are examined and compared, and most of those for which the design is especially suitable for showing this effect had unusually low p values for the number of trials. It is concluded that either the experimenter produced a remarkably strong experimenter effect or quantum random events are predetermined, thereby enabling enhanced PK in suitable experimental designs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Monomer-dimer problem on random planar honeycomb lattice

We consider the monomer-dimer (MD) problem on a random planar honeycomb lattice model, namely, the random multiple chain. This is a lattice system with non-periodic boundary condition, whose generating process is inspired by the growth of single walled zigzag carbon nanotubes. By applying algebraic and combinatorial techniques we establish a calculating expression of the MD partition function for bipartite graphs, which corresponds to the permanent of a matrix. Further, by using the transfer matrix argument we show that the computing problem of the permanent of high order matrix can be converted into some lower order matrices for this family of lattices, based on which we derive an explicit recurrence formula for evaluating the MD partition function of multiple chains and random multiple chains. Finally, we analyze the expectation of the number of monomer-dimer arrangements on a random multiple chain and the asymptotic behavior of the annealed MD entropy when the multiple chain becomes infinite in width and length, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Helicon modes in uniform plasmas. III. Angular momentum

Helicons are electromagnetic waves with helical phase fronts propagating in the whistler mode in magnetized plasmas and solids. They have similar properties to electromagnetic waves with angular momentum in free space. Helicons are circularly polarized waves carrying spin angular momentum and orbital angular momentum due to their propagation around the ambient magnetic field B{sub 0}. These properties have not been considered in the community of researchers working on helicon plasma sources, but are the topic of the present work. The present work focuses on the field topology of helicons in unbounded plasmas, not on helicon source physics. Helicons are excited in a large uniform laboratory plasma with a magnetic loop antenna whose dipole axis is aligned along or across B{sub 0}. The wave fields are measured in orthogonal planes and extended to three dimensions (3D) by interpolation. Since density and B{sub 0} are uniform, small amplitude waves from loops at different locations can be superimposed to generate complex antenna patterns. With a circular array of phase shifted loops, whistler modes with angular and axial wave propagation, i.e., helicons, are generated. Without boundaries radial propagation also arises. The azimuthal mode number m can be positive or negative while the field polarization remains right-hand circular. The conservation of energy and momentum implies that these field quantities are transferred to matter which causes damping or reflection. Wave-particle interactions with fast electrons are possible by Doppler shifted resonances. The transverse Doppler shift is demonstrated. Wave-wave interactions are also shown by showing collisions between different helicons. Whistler turbulence does not always have to be created by nonlinear wave-interactions but can also be a linear superposition of waves from random sources. In helicon collisions, the linear and/or orbital angular momenta can be canceled, which results in a great variety of field topologies. The work will be contrasted to the research on helicon plasma sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Particle Track Classification Using Quantum Associative Memory (Final Technical Report)

This project explored the use of quantum-assisted algorithms for pattern matching in sub-atomic physics experiments. Pattern matching algorithms are commonly employed to prune data of random noise and to help discriminate between signals generated by particle tracks of interest and signals generated by background events. The quantum-assisted algorithms explored in this project were based on an Ising formulation of quantum associative model (QAMM) recall and quantum content-addressable memory (QCAM) recall. The recall is performed by comparing a probe pattern with those stored in a library of patterns encoded in the QAMM/QCAM model. The classification accuracy of QAMM and QCAM recall was determined as a function of detector resolution, noise, and efficiency and pattern density, where pattern density is defined as the ratio of the number of reference signal patterns encoded in the library to each pattern’s length. We found that QAMM achieved high classification accuracy when applied to datasets with low pattern density. QCAM achieved high classification accuracy for datasets with high pattern density and was found to be more robust to detector noise. The project methodology and results are described in detail in our arXiv preprint (arXiv:2011.11848) . This project was conducted by scientists at the Johns Hopkins University Applied Physics Laboratory and Oak Ridge National Laboratory from August 2018 to August 2020 and was supported by DOE grant DE-SC0019497.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1995

Probabilistic Analysis of Large-Scale Composite Structures Using the IPACS Code

An investigation was performed to ascertain the feasibility of using IPACS (Integrated Probabilistic Assessment of Composite Structures) for probabilistic analysis of a composite fan blade, the development of which is being pursued by various industries for the next generation of aircraft engines. A model representative of the class of fan blades used in the GE90 engine has been chosen as the structural component to be analyzed with IPACS. In this study, typical uncertainties are assumed in the level, and structural responses for ply stresses and frequencies are evaluated in the form of cumulative probability density functions. Because of the geometric complexity of the blade, the number of plies varies from several hundred at the root to about a hundred at the tip. This represents a extremely complex composites application for the IPACS code. A sensitivity study with respect to various random variables is also performed.

Lemonds, Jeffrey↗

DOE OSTI2015

Improved Sampling Algorithms in the Risk-Informed Safety Margin Characterization Toolkit

The RISMC approach is developing advanced set of methodologies and algorithms in order to perform Probabilistic Risk Analyses (PRAs). In contrast to classical PRA methods, which are based on Event-Tree and Fault-Tree methods, the RISMC approach largely employs system simulator codes applied to stochastic analysis tools. The basic idea is to randomly perturb (by employing sampling algorithms) timing and sequencing of events and internal parameters of the system codes (i.e., uncertain parameters) in order to estimate stochastic parameters such as core damage probability. This approach applied to complex systems such as nuclear power plants requires to perform a series of computationally expensive simulation runs given a large set of uncertain parameters. These types of analysis are affected by two issues. Firstly, the space of the possible solutions (a.k.a., the issue space or the response surface) can be sampled only very sparsely, and this precludes the ability to fully analyze the impact of uncertainties on the system dynamics. Secondly, large amounts of data are generated and tools to generate knowledge from such data sets are not yet available. This report focuses on the first issue and in particular employs novel methods that optimize the information generated by the sampling process by sampling unexplored and risk-significant regions of the issue space: adaptive (smart) sampling algorithms. They infer system response from surrogate models constructed from existing samples and predict the most relevant location of the next sample. It is therefore possible to understand features of the issue space with a small number of carefully selected samples. In this report, we will present how it is possible to perform adaptive sampling using the RISMC toolkit and highlight the advantages compared to more classical sampling approaches such Monte-Carlo. We will employ RAVEN to perform such statistical analyses using both analytical cases but also another RISMC code: RELAP-7.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Mapping Vegetation at Species Level with High-Resolution Multispectral and Lidar Data Over a Large Spatial Area: A Case Study with Kudzu

Mapping vegetation species is critical to facilitate related quantitative assessment, and mapping invasive plants is important to enhance monitoring and management activities. Integrating high-resolution multispectral remote-sensing (RS) images and lidar (light detection and ranging) point clouds can provide robust features for vegetation mapping. However, using multiple sources of high-resolution RS data for vegetation mapping on a large spatial scale can be both computationally and sampling intensive. Here, we designed a two-step classification workflow to potentially decrease computational cost and sampling effort and to increase classification accuracy by integrating multispectral and lidar data in order to derive spectral, textural, and structural features for mapping target vegetation species. We used this workflow to classify kudzu, an aggressive invasive vine, in the entire Knox County (1362 km2) of Tennessee (U.S.). Object-based image analysis was conducted in the workflow. The first-step classification used 320 kudzu samples and extensive, coarsely labeled samples (based on national land cover) to generate an overprediction map of kudzu using random forest (RF). For the second step, 350 samples were randomly extracted from the overpredicted kudzu and labeled manually for the final prediction using RF and support vector machine (SVM). Computationally intensive features were only used for the second-step classification. SVM had constantly better accuracy than RF, and the producer’s accuracy, user’s accuracy, and Kappa for the SVM model on kudzu were 0.94, 0.96, and 0.90, respectively. SVM predicted 1010 kudzu patches covering 1.29 km2 in Knox County. We found the sample size of kudzu used for algorithm training impacted the accuracy and number of kudzu predicted. The proposed workflow could also improve sampling efficiency and specificity. Our workflow had much higher accuracy than the traditional method conducted in this research, and could be easily implemented to map kudzu in other regions as well as map other vegetation species.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2017

CME Velocity and Acceleration Error Estimates Using the Bootstrap Method

The bootstrap method is used to determine errors of basic attributes of coronal mass ejections (CMEs) visually identified in images obtained by the Solar and Heliospheric Observatory (SOHO) mission's Large Angle and Spectrometric Coronagraph (LASCO) instruments. The basic parameters of CMEs are stored, among others, in a database known as the SOHO/LASCO CME catalog and are widely employed for many research studies. The basic attributes of CMEs (e.g. velocity and acceleration) are obtained from manually generated height-time plots. The subjective nature of manual measurements introduces random errors that are difficult to quantify. In many studies the impact of such measurement errors is overlooked. In this study we present a new possibility to estimate measurements errors in the basic attributes of CMEs. This approach is a computer-intensive method because it requires repeating the original data analysis procedure several times using replicate datasets. This is also commonly called the bootstrap method in the literature. We show that the bootstrap approach can be used to estimate the errors of the basic attributes of CMEs having moderately large numbers of height-time measurements. The velocity errors are in the vast majority small and depend mostly on the number of height-time points measured for a particular event. In the case of acceleration, the errors are significant, and for more than half of all CMEs, they are larger than the acceleration itself.

Michalek, Grzegorz↗

DOE OSTI2020

Data-driven Alignment of the Tracker for the Mu2e Experiment

An alignment of the Mu2e Tracker is investigated using the Mu2e Offline Software and Millepede-II. Simulated misalignments of Tracker planes were inserted between Monte-Carlo event generation and before reconstruction using randomly generated $(x,y,z)$-translations and rotations of the rigid-body planes in the detector model. Cosmic ray muons, with no magnetic field present in the Detector Solenoid, will follow a straight-line trajectory. A maximum-likelihood fit is used to reconstruct a track and determine its parameters. Tracks traversing at least four planes, with ten registered straw hits, and a maximum time residual of $20$ ns are selected to be used in an alignment sample. The cuts were chosen to discard tracks with low information content, such as those with high levels of multiple scattering or those with low numbers of measurements contributing to a single track fit. Millepede-II produces the resulting alignment constants according to a simultaneous least-squares fit of the alignment parameters and the track parameters for all input tracks. The improved alignment constants are used to repeat the process and recover an increased number of suitable tracks, further improving the estimate of the alignment constants in subsequent iterations. A performance study of the process was conducted on a misaligned Tracker in the case of translation-only, and translation plus rotational misalignments. In both cases, the correct alignment constants were recovered to within statistical error after three global fit iterations. Additionally, a continuous integration system was developed for the \texttt{Mu2e/Offline} GitHub repository, introducing a convenient automated testing process for all proposed changes (``pull requests"). The tests developed included a full build of the software, followed by tests of the validation routines, geometry overlap checking, and multi-threading mode. Additionally, code quality checks and low-statistics physics comparison checks a re also available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2012

NSE abundance data

A novel method of calculating nuclear statistical equilibrium (NSE) is presented. Basic equations are carefully solved using arbitrary precision arithmetic. A special interpolation procedure is then used to retrieve all abundances using tabulated results for neutrons and protons, together with basic nuclear data. Proton and neutron abundance tables, basic nuclear data, and partition functions for nuclides used in the calculations are provided. A simple interpolation algorithm using pre-calculated p and n abundances tabulated as functions of kT, {rho} and Y{sub e} is outlined. Unique properties of this method are: (1) ability to pick up out of NSE selected nuclei only, (2) computational time scaling linearly with number of re-calculated abundances, (3) relatively small amount of stored data: only two large tables, (4) slightly faster than solving the NSE equations using traditional Newton-Raphson methods for small networks (few tens of species); superior for huge (800-3000) networks, (5) does not require initial guess; works well on random input, (6) can be tailored to specific application, (7) ability to use third-party NSE solvers to obtain fully compatible tables, and (8) encapsulation of the NSE code for bug-free calculations. A range of applications for this approach is possible: covering tests of traditional NSE Newton-Raphson codes, generating starting values, code-to-code verification, and possible replacement of the old legacy procedures in supernova simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2014

Contrasting Metabolism in Perenniating Structures of Upland and Lowland Switchgrass Plants Late in the Growing Season

Background: Switchgrass (Panicum virgatum L.) is being developed as a bioenergy crop for many temperate regions of the world. One way to increase biomass yields is to move southern adapted lowland cultivars to more northern latitudes. However, many southerly adapted switchgrass germplasm can suffer significant winter kill in northerly climes. Materials and Methods: Here, we have applied next-generation sequencing in combination with biochemical analyses to query the metabolism of crowns and rhizomes obtained from two contrasting switchgrass cultivars. Crowns and rhizomes from field-grown lowland (cv Kanlow) and upland (cv Summer) switchgrass cultivars were collected from three randomly selected post-flowering plants. Summer plants were senescing, whereas Kanlow plants were not at this harvest date. Results: Principal component analysis (PCA) differentiated between both the Summer and Kanlow transcriptomes and metabolomes. Significant differences in transcript abundances were detected for 8,050 genes, including transcription factors such as WRKYs and those associated with phenylpropanoid biosynthesis. Gene-set enrichment analyses showed that a number of pathways were differentially up-regulated in the two populations. For both populations, protein levels and enzyme activities agreed well with transcript abundances for genes involved in the phenylpropanoid pathway that were upregulated in Kanlow crowns and rhizomes. The combination of these datasets suggests that dormancy-related mechanisms had been triggered in the crowns and rhizomes of the Summer plants, whereas the crowns and rhizomes of Kanlow plants had yet to enter dormancy. Conclusions: Delayed establishment of dormancy at more northerly latitudes could be one factor that reduces wintersurvival in the high-yielding Kanlow plants. Understanding the cellular signatures that accompany the transition to dormancy can be used in the future to select plants with improved winter hardiness.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2001

Determining the Mass Density Along Magnetic Field Lines from Toroida Eigenfrequencies

Toroidal eigenfrequencies can be used to remotely sense the equatorial mass density rho(sub eq) and the density dependence along a magnetic field line. Here we present improvements to the method of Schulz [1996], which allows rho(sub eq) and the power law index alpha (for mass density along a field line proportional to R(sup -alpha), where R is the radial distance from the center of the Earth) to be determined from the y intercept and slope of a plot of toroidal frequency versus toroidal harmonic number n. Our modifications include a model form for eigenfrequencies with a fractional precision of 0.0005 for -6 less than or = alpha less than or = 6 and 2 less than or = L less than or = 8 (accuracy is doubtful beyond L = 5) and an iterative procedure for getting more accurate results than those found using Schulz's method. In addition, we do an analysis of the effect of random measurement errors. Observed frequencies need to be accurate to approx. 6% (3%) of the fundamental frequency in order to determine rho(sub eq) (alpha) to a precision of 30% (unity). We then apply our method to data generated using the Global Core Plasma Model for plasmaspheric mass density; our analysis demonstrates clearly bow the alpha index represents the mass density dependence on the outer part of the field line (R/(LR(sub E)) greater than or approx. 2/3).

Denton, R.↗

NASA NTRS2014

An Evaluation of a Passively Cooled Cylindrical Spectrometer Array in Lunar Orbit

This thesis will evaluate a passively cooled cylindrical spectrometer array in lunar orbit characterizing the thermal response in order to provide context for decision-making to scientists and engineers. To provide perspective on thermal issues and controls of space science instruments, a background search of historical lunar missions is provided. Next, a trial science mission is designed and analyzed which brings together the elements of the background search, lunar orbit environment and passive cooling. Finally, the trial science mission analysis results are provided along with the conclusions drawn. Scintillators are materials that when struck by particle radiation, absorb the particle energy which is then reemitted as light in or near the visible range. Nuclear astrophysics utilizes scintillating materials for observation of high-energy photons which are generated by sources such as solar flares, supernovae and neutron stars. SPMs are paired with inorganic scintillators to detect the light emitted which is converted into electronic signals. The signals are captured and analyzed in order to map the number and location of the high-energy sources. The SPM is utilized as it has single photon sensitivity, low voltage requirements and a fast response. SPMs are also compact, relatively inexpensive and allow the usage of lower-cost scintillating materials within the spectrometer. These characteristics permit large-area arrays while lowering cost and power requirements. The ability of a spectrometer to record and identify the interaction of high-energy photons for scientific return is not a trivial matter. Background noise is generated when particles that have not originated from the desired distant source impact the spectrometer. Additionally, thermally induced electrical signals are randomly generated within the SPM even in the absence of light which is referred to as dark current. Overcoming these obstacles requires greater light emittance and energy resolution with reduced dark current. Strong scintillation photon emittance ensures that low energy impacts will produce enough visible photons to be detected by the SPM. Higher energy resolution will ensure that single photon impacts can be distinguished from others of similar wavelength and energy; reduced dark current decreases the generation of random signals not associated with a photon impact. Increasing efficiency in each of these properties in a spectrometer comprised of inorganic scintillators and SPMs requires low temperatures. Low temperature maintenance in a lunar environment presents many unique challenges of its own. Even with the accumulated successes of past missions, the lunar environment remains a thermal challenge for engineers. The lunar orbit thermal environment is driven by radiation from three sources, direct solar radiation, reflected solar radiation from the lunar surface (albedo) and lunar radiation (Clawson 2002). Direct solar radiation values are consistent with those seen in Earth orbit (1325 W/m2) (Clawson 2002). The percentage of solar radiation reflected from the moon is consistently very low with the moon's dark regolith covered surface absorbing nearly 90% of the incident light (Clawson 2002). Yet, it is this absorption that gives the lunar orbit environment one of its most difficult thermal attributes as the absorbed solar radiation is released from the lunar surface as infrared radiation (IR). IR is of a wavelength that is readily absorbed by surfaces designed to function as radiation emitters. It is practical to therefore "choose radiator locations and spacecraft attitude to minimize radiator views to the lunar surface, when possible...pointing the radiator towards the sun to some extent, to minimize its view to the lunar surface, is frequently preferable. (Clawson 2002)" Additionally, the amount of direct solar radiation, lunar IR and albedo an orbiting satellite receives varies from one side of the moon to the other as the moon blocks the sun from view. This environment produces large temperature variations in a satellite's instrumentation, control electronics and propulsion systems which must be understood to characterize operating temperature envelopes.

Waggoner, Jason↗

DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗