Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random number generators”

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 325 records · Page 18

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↗

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.↗

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↗

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↗

Silhouette-Informed Trajectory Generation Through a Wire Maze for Small UAS

Current rapidly-exploring random tree (RRT) algorithms rely on proximity query packages that often include collision checkers, tolerance verification, and distance computation algorithms for the generation of safe paths. In this paper, we broaden the information available to the path-planning algorithm by incorporating silhouette information of nearby obstacles in conflict. A silhouette-informed tree (SIT) is generated through the flight-safe region of a wire maze for a single unmanned aerial system (UAS). The silhouette is used to extract local geometric information of nearby obstacles and provide path alternatives around these obstacles. Thus, focusing the search for the generation of new tree branches near these obstacles, and decreasing the number of samples required to explore the narrow corridors within the wire maze. The SIT is then processed to extract a path that connects the initial location of the UAS with the goal, reduce the number of line segments in this path if possible, and smooth the resulting path using Pythagorean Hodograph Bezier curves. To ensure that the smoothed path remains in the flight-safe region of the configuration space, a tolerance verification algorithm for Bezier curves and convex polytopes in three dimensions is proposed. Lastly, temporal specifications are imposed on the smoothed path in the shape of an arbitrary speed profile.

Puig-Navarro, Javier↗

GLAD-M35: a joint P and S global tomographic model with uncertainty quantification

We present our third and final generation joint P and S global adjoint tomography (GLAD) model, GLAD-M35, and quantify its uncertainty based on a low-rank approximation of the inverse Hessian. Starting from our second-generation model, GLAD-M25, we added 680 new earthquakes to the database for a total of 2160 events. New P-wave categories are included to compensate for the imbalance between P- and S-wave measurements, and we enhanced the window selection algorithm to include more major-arc phases, providing better constraints on the structure of the deep mantle and more than doubling the number of measurement windows to 40 million. Two stages of a Broyden–Fletcher–Goldfarb–Shanno (BFGS) quasi-Newton inversion were performed, each comprising five iterations. With this BFGS update history, we determine the model’s standard deviation and resolution length through randomized singular value decomposition.

58 GEOSCIENCES↗

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population’s mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn .

59 BASIC BIOLOGICAL SCIENCES↗

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population’s mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn .

59 BASIC BIOLOGICAL SCIENCES↗

Effects of Dose Error and Sample Size on Sonic Boom Dose-response Curves

NASA will soon be collecting noise-annoyance community survey data as the X-59 aircraft flies supersonically over several communities in the USA. Sparse measurements of the X-59 sonic thumps will be used together with physics-based simulations to estimate noise doses at survey participant locations. These dose estimates have associated error that affects the accuracy of modeled dose-response curves, which can result in misestimation of annoyance. The precision in dose-response curves is also a consideration in selecting the number of survey participants. To enable pretest studies of dose error and precision, simulated dose-response data were generated based on NASA’s Quiet Supersonic Flights 2018 test. The data included various degrees of dose error and sample size. Frequentist multilevel logistic regression models were fit to the true and perturbed dose-response data. Simple proportional relationships were identified between the model parameters and the perturbation standard deviation. The summary dose-response curves illustrate the impact on accuracy if dose error is not accounted for in the model. The precision in the dose-response curves is also shown as the number of participants and degree of participation is varied. Finally, sampling variability is illustrated by showing the dose-response curves for several replicates with random draws of participants and errors.

X-59↗

Effects of Dose Error and Sample Size on Sonic Boom Dose-Response Curves

NASA will soon be collecting noise-annoyance community survey data as the X-59 aircraft flies supersonically over several communities in the USA. Sparse measurements of the X-59 sonic thumps will be used together with physics-based simulations to estimate noise doses at survey participant locations. These dose estimates have associated error that affects the accuracy of modeled dose-response curves, which can result in misestimation of annoyance. The precision in dose-response curves is also a consideration in selecting the number of survey participants. To enable pretest studies of dose error and precision, simulated dose-response data were generated based on NASA’s Quiet Supersonic Flights 2018 test. The data included various degrees of dose error and sample size. Frequentist multilevel logistic regression models were fit to the true and perturbed dose-response data. Simple proportional relationships were identified between the model parameters and the perturbation standard deviation. The summary dose-response curves illustrate the impact on accuracy if dose error is not accounted for in the model. The precision in the dose-response curves is also shown as the number of participants and degree of participation is varied. Finally, sampling variability is illustrated by showing the dose-response curves for several replicates with random draws of participants and errors.

X-59↗

Some aspects of the aeroacoustics of high-speed jets

Some of the background to contemporary jet aeroacoustics is addressed. Then scaling laws for noise generation by low-Mach-number airflows and by turbulence convected at 'not so low' Mach number is reviewed. These laws take into account the influence of Doppler effects associated with the convection of aeroacoustic sources. Next, a uniformly valid Doppler-effect approximation exhibits the transition, with increasing Mach number of convection, from compact-source radiation at low Mach numbers to a statistical assemblage of conical shock waves radiated by eddies convected at supersonic speed. In jets, for example, supersonic eddy convection is typically found for jet exit speeds exceeding twice the atmospheric speed of sound. The Lecture continues by describing a new dynamical theory of the nonlinear propagation of such statistically random assemblages of conical shock waves. It is shown, both by a general theoretical analysis and by an illustrative computational study, how their propagation is dominated by a characteristic 'bunching' process. That process associated with a tendency for shock waves that have already formed unions with other shock waves to acquire an increased proneness to form further unions - acts so as to enhance the high-frequency part of the spectrum of noise emission from jets at these high exit speeds.

Lighthill, James↗

An alternative resource sharing scheme for land mobile satellite services

A preliminary comparison between the two competing channelization concepts for the Land Mobile Satellite Services (LMSS), namely frequency division (FD) and code division (CD), is presented. Both random access and demand-assigned approaches are considered under these concepts. The CD concept is compared with the traditional FD concept based on the system consideration and a projected traffic model. It is shown that CD is not particularly attractive for the first generation Mobile Satellite Services because of the spectral occupancy of the network bandwidth. However, the CD concept is a viable alternative for future systems such as the personal access satellite system (PASS) in the Ka-band spectrum where spectral efficiency is not of prime concern. The effects of power robbing and voice activity factor are incorporated. It was shown that the traditional rule of thumb of dividing the number of raw channels by the voice activity factor to obtain the effective number of channels is only valid asymptotically as the aggregated traffic approaches infinity.

Yan, Tsun-Yee↗

Benchmarking the performance of a high-Q cavity qudit using random unitaries

High-coherence cavity resonators are excellent resources for encoding quantum information in higher-dimensional Hilbert spaces, moving beyond traditional qubit-based platforms. A natural strategy is to use the Fock basis to encode information in qudits. One can perform quantum operations on the cavity mode qudit by coupling the system to a non-linear ancillary transmon qubit. However, the performance of the cavity-transmon device is limited by the noisy transmons. It is, therefore, important to develop practical benchmarking tools for these qudit systems in an algorithm-agnostic manner. We gauge the performance of these qudit platforms using sampling tests such as the heavy output generation test as well as the linear cross-entropy benchmark, by way of simulations of such a system subject to realistic dominant noise channels. We use selective number-dependent arbitrary phase and unconditional displacement gates as our universal gateset. Our results show that contemporary transmons comfortably enable controlling a few tens of Fock levels of a cavity mode. This framework allows benchmarking even higher dimensional qudits as those become accessible with improved transmons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computational Power of Random Quantum Circuits in Arbitrary Geometries

Empirical evidence for a gap between the computational powers of classical and quantum computers has been provided by experiments that sample the output distributions of two-dimensional quantum circuits. Many attempts to close this gap have utilized classical simulations based on tensor network techniques, and their limitations shed light on the improvements to quantum hardware required to frustrate classical simulability. In particular, quantum computers having in excess of approximately 50 qubits are primarily vulnerable to classical simulation due to restrictions on their gate fidelity and their connectivity, the latter determining how many gates are required (and, therefore, how much infidelity is suffered) in generating highly entangled states. Here, we describe recent hardware upgrades to Quantinuum’s H2 quantum computer, enabling it to operate on up to 56 qubits with arbitrary connectivity and 99.843(5)% two-qubit gate fidelity. We define a class of circuits with random geometries that become hard to classically simulate in very low depth and implement them utilizing the flexible connectivity of H2. A careful analysis demonstrating the fast saturation of classical simulation complexity with depth indicates that H2 can yield data well beyond the reach of state-of-the art classical simulation methods at unprecedented fidelities. We find that the considerable difficulty of classically simulating H2 is likely limited only by qubit number, demonstrating the promise and scalability of the quantum charge-coupled device architecture as continued progress is made toward building larger machines. Published by the American Physical Society 2025

DeCross, M.↗

Complexity-calibrated benchmarks for machine learning reveal when prediction algorithms succeed and mislead

Abstract Recurrent neural networks are used to forecast time series in finance, climate, language, and from many other domains. Reservoir computers are a particularly easily trainable form of recurrent neural network. Recently, a “next-generation” reservoir computer was introduced in which the memory trace involves only a finite number of previous symbols. We explore the inherent limitations of finite-past memory traces in this intriguing proposal. A lower bound from Fano’s inequality shows that, on highly non-Markovian processes generated by large probabilistic state machines, next-generation reservoir computers with reasonably long memory traces have an error probability that is at least $$\sim 60\%$$ ∼ 60 % higher than the minimal attainable error probability in predicting the next observation. More generally, it appears that popular recurrent neural networks fall far short of optimally predicting such complex processes. These results highlight the need for a new generation of optimized recurrent neural network architectures. Alongside this finding, we present concentration-of-measure results for randomly-generated but complex processes. One conclusion is that large probabilistic state machines—specifically, large $$\epsilon$$ ϵ -machines—are key to generating challenging and structurally-unbiased stimuli for ground-truthing recurrent neural network architectures.

97 MATHEMATICS AND COMPUTING↗

Application of Machine Learning Algorithms to the Study of Noise Artifacts in Gravitational-Wave Data

The sensitivity of searches for astrophysical transients in data from the Laser Interferometer Gravitationalwave Observatory (LIGO) is generally limited by the presence of transient, non-Gaussian noise artifacts, which occur at a high-enough rate such that accidental coincidence across multiple detectors is non-negligible. Furthermore, non-Gaussian noise artifacts typically dominate over the background contributed from stationary noise. These "glitches" can easily be confused for transient gravitational-wave signals, and their robust identification and removal will help any search for astrophysical gravitational-waves. We apply Machine Learning Algorithms (MLAs) to the problem, using data from auxiliary channels within the LIGO detectors that monitor degrees of freedom unaffected by astrophysical signals. Terrestrial noise sources may manifest characteristic disturbances in these auxiliary channels, inducing non-trivial correlations with glitches in the gravitational-wave data. The number of auxiliary-channel parameters describing these disturbances may also be extremely large; high dimensionality is an area where MLAs are particularly well-suited. We demonstrate the feasibility and applicability of three very different MLAs: Artificial Neural Networks, Support Vector Machines, and Random Forests. These classifiers identify and remove a substantial fraction of the glitches present in two very different data sets: four weeks of LIGO's fourth science run and one week of LIGO's sixth science run. We observe that all three algorithms agree on which events are glitches to within 10% for the sixth science run data, and support this by showing that the different optimization criteria used by each classifier generate the same decision surface, based on a likelihood-ratio statistic. Furthermore, we find that all classifiers obtain similar limiting performance, suggesting that most of the useful information currently contained in the auxiliary channel parameters we extract is already being used. Future performance gains are thus likely to involve additional sources of information, rather than improvements in the MLAs themselves.

gravitational-wave data↗

Step Bunching: Influence of Impurities and Solution Flow

Step bunching results in striations even at relatively early stages of its development and in inclusions of mother liquor at the later stages. Therefore, eliminating step bunching is crucial for high crystal perfection. At least 5 major effects causing and influencing step bunching are known: (1) Basic morphological instability of stepped interfaces. It is caused by concentration gradient in the solution normal to the face and by the redistribution of solute tangentially to the interface which redistribution enhances occasional perturbations in step density due to various types of noise; (2) Aggravation of the above basic instability by solution flowing tangentially to the face in the same directions as the steps or stabilization of equidistant step train if these flows are antiparallel; (3) Enhanced bunching at supersaturation where step velocity v increases with relative supersaturation s much faster than linear. This v(s) dependence is believed to be associated with impurities. The impurities of which adsorption time is comparable with the time needed to deposit one lattice layer may also be responsible for bunching; (4) Very intensive solution flow stabilizes growing interface even at parallel solution and step flows; (5) Macrosteps were observed to nucleate at crystal corners and edges. Numerical simulation, assuming step-step interactions via surface diffusion also show that step bunching may be induced by random step nucleation at the facet edge and by discontinuity in the step density (a ridge) somewhere in the middle of a face. The corresponding bunching patterns produce the ones observed in experiment. The nature of step bunching generated at the corners and edges and by dislocation step sources, as well as the also relative importance and interrelations between mechanisms 1-5 is not clear, both from experimental and theoretical standpoints. Furthermore, several laws controlling the evolution of existing step bunches have been suggested, though unambiguous conclusions are still missing. Addressing these issues is the major goal of the present project. The theory addressing the above problem, experimental methods, several figures which include: (1) the spatial wave numbers at which the system is neutrally stable as a function of growth velocity for linear kinetics and supersaturation for nonlinear kinetics; (2) a schematic of the experiment of lysozyme crystal growing under conditions of natural convection; (3) fluctuations in time, t, of the normal growth rate, R(t), vicinal slope, p(t) and Fourier Spectra of R(t), discussions and conclusions are presented.

Chernov, A. A.↗

Two Solvers for Tractable Temporal Constraints with Preferences

A number of reasoning problems involving the manipulation of temporal information can naturally be viewed as implicitly inducing an ordering of potential local decisions involving time on the basis of preferences. Soft temporal constraints problems allow to describe in a natural way scenarios where events happen over time and preferences are associated to event distances and durations. In general, solving soft temporal problems require exponential time in the worst case, but there are interesting subclasses of problems which are polynomially solvable. We describe two solvers based on two different approaches for solving the same tractable subclass. For each solver we present the theoretical results it stands on, a description of the algorithm and some experimental results. The random generator used to build the problems on which tests are performed is also described. Finally, we compare the two solvers highlighting the tradeoff between performance and representational power.

Rossi, F.↗