Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Randomized methods”

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 1,135 records · Page 63

Monte Carlo Dropout Uncertainty Quantification of Long Short-Term Memory Autoencoder Anomaly Detection in a Liquid Sodium Cold Trap

Advanced high-temperature fluid reactors, such as sodium-cooled fast reactors (SFRs) and molten salt–cooled reactors (MSCRs), require coolant purification systems to prevent fluid contamination and local freezing that can lead to plugging. Liquid sodium purification can be achieved with a cold trap, where the sodium temperature is reduced to a near-freezing point to precipitate out impurities. Automation of monitoring of the cold trap performance with machine learning algorithms can aid in early detection of incipient anomalies. An efficient approach to loss-of-coolant–type anomaly detection in a cold trap monitored with more than two dozen thermal-hydraulic sensors consists of a long short-term memory (LSTM) autoencoder. This work develops the uncertainty quantification of the LSTM autoencoder performance for cold trap anomaly detection using the Monte Carlo (MC) dropout method. The MC dropout methodology creates a distribution of sister distributions that all slightly differ from each other because of random neurons being turned off for testing. The variances of the sister network distributions are used to make an uncertainty interval. Our analysis shows that the uncertainty in the autoencoder performance is largest near the peak of the anomaly signal. Using the MC dropout method, we investigate the uncertainty in the anomaly detection with missing sensor inputs. This capability allows the reactor operator to evaluate resilience of the anomaly detection system and to make informed decisions about continuity of operation in the event of sensor failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computational Study of Low-Speed Boundary-Layer Laminar-Turbulent Transition over a 30º Swept Plate in a Channel

Reliable boundary-layer laminar-turbulent transition predictions are critical for the design of aircraft that use natural or controlled laminar flow to achieve low drag and better fuel efficiency. During cruise flight through low turbulence, the three-dimensional boundary layer over swept surfaces can transition prematurely due to stationary crossflow instabilities induced by surface roughness. Parabolized stability equations analyses have been successfully used to explain the different phases of the transition process under well controlled perturbations such as steps or discrete roughness elements. However, methodologies are still lacking to model the boundary layer receptivity to distributed random roughness, and the following non-linear instability stages leading to transition. For engineering applications, several semi-empirical methods based on linear amplification factors or local helicity have found growing adoption, but their reliability is limited by their level of simplification and the sparsity of data they are based on. To address this gap, a new study of the effect of distributed random surface roughness on crossflow transition was started, using a 30◦ swept plate as a test vehicle. The plate, which is mounted in the 2 Foot × 3 Foot Low-Speed Boundary-Layer Channel at the NASA Langley Research Center, has been tested with micro-scale roughness of varied geometry near the leading edge. In order to analyze the flow in the test section and prepare Direct Numerical Simulations with resolved roughness, Reynolds-Averaged Navier-Stokes (RANS) simulations of the flow over the plate installed in the wind tunnel test section were performed using the Langtry-Menter transition model with stationary crossflow extension. This paper presents the results of these simulations, with comparison to the experiment.

ARMD↗

Probabilistic Forecasting of Generators Startups and Shutdowns in the MISO System Based on Random Forest

Solving security constrained unit commitment (SCUC) problems to plan an economical generation schedule for day-head electricity market has been an important research topic in recent years. Mixed integer programming method (MIP), the-state-of-art approach for solving SCUC problem, is known computationally hard when the number of binary status variables is large. In this paper, a machine learning-based algorithm - random forest (RF), was applied to forecast the startups (SU) and shutdowns (SD) hours of generators, based on historical hourly system condition observations in the Midcontinent Independent System Operator (MISO) system. The main purpose is to reduce the number of binary status variables, by fixing the SU/SD hours to a narrow range of high confidence. This would significantly reduce the size of the decision space, and therefore speed up SCUC solutions with reduced uncertainty.

Lin, Xinming↗

Selective norepinephrine reuptake inhibition as a human model of orthostatic intolerance

BACKGROUND: Observations in patients with functional mutations of the norepinephrine transporter (NET) gene suggest that impaired norepinephrine uptake may contribute to idiopathic orthostatic intolerance. METHODS AND RESULTS: We studied the effect of the selective NET blocker reboxetine and placebo in a randomized, double-blind, crossover fashion on cardiovascular responses to cold pressor testing, handgrip testing, and a graded head-up tilt test (HUT) in 18 healthy subjects. In a subset, we determined isoproterenol and phenylephrine sensitivities. Subjects ingested 8 mg reboxetine or placebo 12 hours and 1 hour before testing. In the supine position, heart rate was 65+/-2 bpm with placebo and 71+/-3 bpm with reboxetine. At 75 degrees HUT, heart rate was 84+/-3 and 119+/-4 bpm with placebo and with reboxetine (P<0.0001). Mean arterial pressure was 85+/-2 with placebo and 91+/-2 mm Hg with reboxetine while supine (P<0.01) and 88+/-2 mm Hg and 90+/-3 mm Hg at 75 degrees HUT. Blood pressure responses to cold pressor and handgrip testing were attenuated with reboxetine. Reboxetine increased the sensitivity to the chronotropic effect of isoproterenol and the pressor effect of phenylephrine. Vasovagal reactions occurred in 9 subjects on placebo and in 1 subject on reboxetine. CONCLUSIONS: Selective NET blockade creates a phenotype that resembles idiopathic orthostatic intolerance. This observation supports the hypothesis that disordered norepinephrine uptake mechanisms can contribute to human cardiovascular disease. Our study also suggests that NET inhibition might be useful in preventing vasovagal reactions.

Non-NASA Center↗

Rover mast calibration, exact camera pointing, and camara handoff for visual target tracking

This paper presents three technical elements that we have developed to improve the accuracy of the visual target tracking for single-sol approach-and-instrument placement in future Mars rover missions. An accurate, straightforward method of rover mast calibration is achieved by using a total station, a camera calibration target, and four prism targets mounted on the rover. The method was applied to Rocky8 rover mast calibration and yielded a 1.1-pixel rms residual error. Camera pointing requires inverse kinematic solutions for mast pan and tilt angles such that the target image appears right at the center of the camera image. Two issues were raised. Mast camera frames are in general not parallel to the masthead base frame. Further, the optical axis of the camera model in general does not pass through the center of the image. Despite these issues, we managed to derive non-iterative closed-form exact solutions, which were verified with Matlab routines. Actual camera pointing experiments aver 50 random target image paints yielded less than 1.3-pixel rms pointing error. Finally, a purely geometric method for camera handoff using stereo views of the target has been developed. Experimental test runs show less than 2.5 pixels error on high-resolution Navcam for Pancam-to-Navcam handoff, and less than 4 pixels error on lower-resolution Hazcam for Navcam-to-Hazcam handoff.

camera handoff↗

A survey of unsupervised learning methods for high-dimensional uncertainty quantification in black-box-type problems

Constructing surrogate models for uncertainty quantification (UQ) on complex partial differential equations (PDEs) having inherently high-dimensional O(10 n ), n ≥ 2, stochastic inputs (e.g., forcing terms, boundary conditions, initial conditions) poses tremendous challenges. The “curse of dimensionality” can be addressed with suitable unsupervised learning techniques used as a pre-processing tool to encode inputs onto lower-dimensional subspaces while retaining its structural information and meaningful properties. In this work, we review and investigate thirteen dimension reduction methods including linear and nonlinear, spectral, blind source separation, convex and non-convex methods and utilize the resulting embeddings to construct a mapping to quantities of interest via polynomial chaos expansions (PCE). Here, we refer to the general proposed approach as manifold PCE (m-PCE), where manifold corresponds to the latent space resulting from any of the studied dimension reduction methods. To investigate the capabilities and limitations of these methods we conduct numerical tests for three physics-based systems (treated as black-boxes) having high-dimensional stochastic inputs of varying complexity modeled as both Gaussian and non-Gaussian random fields to investigate the effect of the intrinsic dimensionality of input data. We demonstrate both the advantages and limitations of the unsupervised learning methods and we conclude that a suitable m-PCE model provides a cost-effective approach compared to alternative algorithms proposed in the literature, including recently proposed expensive deep neural network-based surrogates and can be readily applied for high-dimensional UQ in stochastic PDEs.

42 ENGINEERING↗

Fluctuation analysis of relativistic nucleus-nucleus collisions in emulsion chambers

An analytical technique was developed for identifying enhanced fluctuations in the angular distributions of secondary particles produced from relativistic nucleus-nucleus collisions. The method is applied under the assumption that the masses of the produced particles are small compared to their linear momenta. The importance of particles rests in the fact that enhanced fluctuations in the rapidity distributions is considered to be an experimental signal for the creation of the quark-gluon-plasma (QGP), a state of nuclear matter predicted from the quantum chromodynamics theory (QCD). In the approach, Monte Carlo simulations are employed that make use of a portable random member generator that allow the calculations to be performed on a desk-top computer. The method is illustrated with data taken from high altitude emulsion exposures and is immediately applicable to similar data from accelerator-based emulsion exposures.

Mcguire, Stephen C.↗

Dimension-adaptive machine learning-based quantum state reconstruction

Here, we introduce an approach for performing quantum state reconstruction on systems of n qubits using a machine learning-based reconstruction system trained exclusively on m qubits, where m ≥ n. This approach removes the necessity of exactly matching the dimensionality of a system under consideration with the dimension of a model used for training. We demonstrate our technique by performing quantum state reconstruction on randomly sampled systems of one, two, and three qubits using machine learning-based methods trained exclusively on systems containing at least one additional qubit. The reconstruction time required for machine learning-based methods scales significantly more favorably than the training time; hence this technique can offer an overall saving of resources by leveraging a single neural network for dimension-variable state reconstruction, obviating the need to train dedicated machine learning systems for each Hilbert space.

42 ENGINEERING↗

AI and extreme scale computing to learn and infer the physics of higher order gravitational wave modes of quasi-circular, spinning, non-precessing black hole mergers

We use artificial intelligence (AI) to learn and infer the physics of higher order gravitational wave modes of quasi-circular, spinning, non precessing binary black hole mergers. We trained AI models using 14 million waveforms, produced with the surrogate model NRHybSur3dq8, that include modes up to $\ell$ ≤ 4 and (5,5), except for (4,0) and (4,1), that describe binaries with mass-ratios $\textit{q}$ ≤ 8, individual spins $s^z_{\{1,2\}} \in$[–0.8,0.8], and inclination angle $θ \in$ [0,π]. Our probabilistic AI surrogates can accurately constrain the mass-ratio, individual spins, effective spin, and inclination angle of numerical relativity waveforms that describe such signal manifold. We compared the predictions of our AI models with Gaussian process regression, random forest, k-nearest neighbors, and linear regression, and with traditional Bayesian inference methods through the PyCBC Inference toolkit, finding that AI outperforms all these approaches in terms of accuracy, and are between three to four orders of magnitude faster than traditional Bayesian inference methods. Our AI surrogates were trained within 3.4 hours using distributed training on 1,536 NVIDIA V100 GPUs in the Summit supercomputer.

79 ASTRONOMY AND ASTROPHYSICS↗

Improvised Field Expedient Method for Renal Replacement Therapy in a Porcine Model of Acute Kidney Injury

ABSTRACT Objective: Dialysis patients may not have access to conventional renal replacement therapy (RRT) following disasters. We hypothesized that improvised renal replacement therapy (ImpRRT) would be comparable to continuous renal replacement therapy (CRRT) in a porcine acute kidney injury model. Methods: Following bilateral nephrectomies and 2 hours of caudal aortic occlusion, 12 pigs were randomized to 4 hours of ImpRRT or CRRT. In the ImpRRT group, blood was circulated through a dialysis filter using a rapid infuser to collect the ultrafiltrate. Improvised replacement fluid, made with stock solutions, was infused pre-pump. In the CRRT group, commercial replacement fluid was used. During RRT, animals received isotonic crystalloids and norepinephrine. Results: There were no differences in serum creatinine, calcium, magnesium, or phosphorus concentrations. While there was a difference between groups in serum potassium concentration over time ( P < 0.001), significance was lost in pairwise comparison at specific time points. Replacement fluids or ultrafiltrate flows did not differ between groups. There were no differences in lactate concentration, isotonic crystalloid requirement, or norepinephrine doses. No difference was found in electrolyte concentrations between the commercial and improvised replacement solutions. Conclusion: The ImpRRT system achieved similar performance to CRRT and may represent a potential option for temporary RRT following disasters.

Public, Environmental & Occupational Health↗

Cross-correlation spectra in interacting quantum dot systems

Two-color spin-noise spectroscopy of interacting electron spins in singly charged semiconductor quantum dots provides information on the interquantum dot interactions. We investigate the spin cross-correlation function in a quantum dot ensemble employing a modified semiclassical approach. Spin-correlation functions are calculated using a Hamilton quaternion approach that maintains local quantum mechanical properties of the spins. This method takes into account the effects of the nuclear-electric quadrupolar interactions, the randomness of the coupling constants, and the variation of the electron g factor on the spin-noise power spectra. We demonstrate that the quantum dot ensemble can be mapped on an effective two-quantum dot problem and discuss how the characteristic length scale of the interdot interaction modifies the low-frequency cross-correlation spectrum. Here, we argue that details on the interaction strength distribution can be extracted from the cross-correlation spectrum when applying a longitudinal or a transversal external magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Microscopic predictions for the production of neutron-rich nuclei in the reaction Yb 176 + Yb 176

Background: Production of neutron-rich nuclei is of vital importance to both understanding nuclear structure far from stability and to informing astrophysical models of the rapid neutron capture process (r-process). Multinucleon transfer (MNT) in heavy-ion collisions offers a possibility to produce neutron-rich nuclei far from stability. Purpose: The 176 Yb + 176 Yb reaction has been suggested as a potential candidate to explore the neutron-rich region surrounding the principal fragments. The current study has been conducted with the goal of providing guidance for future experiments wishing to study this (or similar) system. Methods: Time-dependent Hartree-Fock (TDHF) and its time-dependent random-phase approximation (TDRPA) extension are used to examine both scattering and MNT characteristics in 176 Yb + 176 Yb. TDRPA calculations are performed to compute fluctuations and correlations of the neutron and proton numbers, allowing for estimates of primary fragment production probabilities. Results: Both scattering results from TDHF and transfer results from the TDRPA are presented for different energies, orvientations, and impact parameters. In addition to fragment composition, scattering angles and total kinetic energies, as well as correlations between these observables are presented. Conclusions: 176 Yb + 176 Yb appears to be an interesting probe for the midmass neutron-rich region of the chart of nuclides. Lastly, the predictions of both TDHF and TDRPA are speculative, and will benefit from future experimental results to test the validity of this approach to studying MNT in heavy, symmetric collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effects of quasiparticle-vibration coupling on Gamow-Teller strength and 𝛽 decay with the Skyrme proton-neutron finite-amplitude method

We adapt the proton-neutron finite-amplitude method, which in its original form is an efficient implementation of the Skyrme quasiparticle random phase approximation, to include the coupling of quasiparticles to like-particle phonons. The approach allows us to add beyond-quasiparticle random-phase approximation correlations to computations of Gamow-Teller strength and 𝛽-decay rates in deformed nuclei for the first time. We test the approach in several deformed isotopes for which measured strength distributions are available. Furthermore, the additional correlations dramatically improve agreement with the data, and will lead to improved global 𝛽-decay rates.

Beta decay↗

Scattering of Magnons at Graphene Quantum-Hall-Magnet Junctions

Motivated by recent nonlocal transport studies of quantum-Hall-magnet (QHM) states formed in monolayer graphene’s N = 0 Landau level, here we study the scattering of QHM magnons by gate-controlled junctions between states with different integer filling factors ν . For the ν = 1 | – 1 | 1 geometry we find that magnons are weakly scattered by electric potential variation in the junction region, and that the scattering is chiral when the junction lacks a mirror symmetry. For the ν = 1 | 0 | 1 geometry, we find that kinematic constraints completely block magnon transmission if the incident angle exceeds a critical value. Our results explain the suppressed nonlocal–voltage signals observed in the ν = 1 | 0 | 1 case. We use our theory to propose that valley waves generated at ν = – 1 | 1 junctions and magnons can be used in combination to probe the spin or valley flavor structure of QHM states at integer and fractional filling factors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Toward Accelerating Discovery via Physics-Driven and Interactive Multifidelity Bayesian Optimization

Both computational and experimental material discovery bring forth the challenge of exploring multidimensional and often nondifferentiable parameter spaces, such as phase diagrams of Hamiltonians with multiple interactions, composition spaces of combinatorial libraries, processing spaces, and molecular embedding spaces. Often these systems are expensive or time consuming to evaluate a single instance, and hence classical approaches based on exhaustive grid or random search are too data intensive. This resulted in strong interest toward active learning methods such as Bayesian optimization (BO) where the adaptive exploration occurs based on human learning (discovery) objective. However, classical BO is based on a predefined optimization target, and policies balancing exploration and exploitation are purely data driven. In practical settings, the domain expert can pose prior knowledge of the system in the form of partially known physics laws and exploration policies often vary during the experiment. Here, we propose an interactive workflow building on multifidelity BO (MFBO), starting with classical (data-driven) MFBO, then expand to a proposed structured (physics-driven) structured MFBO (sMFBO), and finally extend it to allow human-in-the-loop interactive interactive MFBO (iMFBO) workflows for adaptive and domain expert aligned exploration. These approaches are demonstrated over highly nonsmooth multifidelity simulation data generated from an Ising model, considering spin–spin interaction as parameter space, lattice sizes as fidelity spaces, and the objective as maximizing heat capacity. Detailed analysis and comparison show the impact of physics knowledge injection and real-time human decisions for improved exploration with increased alignment to ground truth. Here, the associated notebooks allow to reproduce the reported analyses and apply them to other systems.

97 MATHEMATICS AND COMPUTING↗

A novel statistical methodology for quantifying the spatial arrangements of axons in peripheral nerves

A thorough understanding of the neuroanatomy of peripheral nerves is required for a better insight into their function and the development of neuromodulation tools and strategies. In biophysical modeling, it is commonly assumed that the complex spatial arrangement of myelinated and unmyelinated axons in peripheral nerves is random, however, in reality the axonal organization is inhomogeneous and anisotropic. Present quantitative neuroanatomy methods analyze peripheral nerves in terms of the number of axons and the morphometric characteristics of the axons, such as area and diameter. In this study, we employed spatial statistics and point process models to describe the spatial arrangement of axons and Sinkhorn distances to compute the similarities between these arrangements (in terms of first- and second-order statistics) in various vagus and pelvic nerve cross-sections. We utilized high-resolution transmission electron microscopy (TEM) images that have been segmented using a custom-built high-throughput deep learning system based on a highly modified U-Net architecture. Our findings show a novel and innovative approach to quantifying similarities between spatial point patterns using metrics derived from the solution to the optimal transport problem. We also present a generalizable pipeline for quantitative analysis of peripheral nerve architecture. Our data demonstrate differences between male- and female-originating samples and similarities between the pelvic and abdominal vagus nerves.

59 BASIC BIOLOGICAL SCIENCES↗

Complex Disorder in Type-I Clathrates: Synthesis and Structural Characterization of A8GaxSn46−x (A = Rb, Cs; 6.9 < x < 7.5)

Exploratory studies in the systems Rb–Ga–Sn and Cs–Ga–Sn yielded the cubic type-I clathrates with refined compositions Rb8GaxSn46−x and Cs8GaxSn46−x (6.9 < x < 7.5). Nearly single-phase materials with good crystallinity were obtained from stoichiometric reactions of the elements. The structures were characterized by means of single-crystal X-ray diffraction methods. Both Rb8GaxSn46−x and Cs8GaxSn46−x represents cases, where a Group 13 element randomly substitutes a Group 14 element in the structure. The extent of Ga/Sn mixing is apparently governed by the drive of the system to achieve an optimal valence electron count, and hence, Rb8GaxSn46−x and Cs8GaxSn46−x (x ≈ 8) can be regarded as Zintl phases. This notion is supported by structure refinements on a multitude of single-crystal X-ray diffraction data, which also confirm that both types of cages in the cubic type-I structure are fully occupied by Rb and Cs atoms. The open-framework, comprised of 46 nodes per formula unit, adapts to the incorporation of nearly eight Ga atoms within the matrix of Sn, whereby small, short-range distortions result. The exact nature of these effects is still unclear, as so far, the structural variations could only be modeled as both positional and occupational disorder at one of three framework sites. Since vacancies in the structures of the binary type-I clathrates A8Sn46−x☐x (A = Rb, Cs; ☐ = missing Sn atom) are also known to cause local distortions, the latter were also synthesized with the same protocols used for the synthesis of A8GaxSn46−x and structurally re-analyzed. The results from the latter studies confirm that homogeneity issues abound, and that the final structures/compositions are an intricate function of the experimental conditions.

36 MATERIALS SCIENCE↗