Search NASA⌕ Search

SEARCH · Search NASA

Results for “probability of initiation”

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 73 records · Page 4

Time-dependent relaxation of observables in complex quantum systems

Abstract We consider time-dependent relaxation of observables in quantum systems of chaotic and regular type. Using statistical arguments and exact numerical solutions we show that the spread of the initial wave function in the Hilbert space and the main characteristics of evolution of observables have certain generic features. The study compares examples of regular dynamics, a completely chaotic case of the Gaussian orthogonal ensemble, a bosonic system with random interactions, and a fully realistic case of the time evolution of various initial non-stationary states in the nuclear shell model. In the case of the Gaussian orthogonal ensemble we show that the survival probability obtained analytically also fully defines the relaxation timescale of observables. This is not the case in general. Using the realistic nuclear shell model and the quadrupole moment as an observable we demonstrate that the relaxation time is significantly longer than defined by the survival probability of the initial state. The full analysis does not show the presence of an analog of the Lyapunov exponent characteristic for examples of classical chaos.

Volya, Alexander (ORCID:0000000217656466)↗

Reactive Resonances in N+N2 Exchange Reaction

Rich reactive resonances are found in a 3D quantum dynamics study of the N + N2 exchange reaction using a recently developed ab initio potential energy surface. This surface is characterized by a feature in the interaction region called Lake Eyring , that is, two symmetric transition states with a shallow minimum between them. An L2 analysis of the quasibound states associated with the shallow minimum confirms that the quasibound states associated with oscillations in all three degrees of freedom in Lake Eyring are responsible for the reactive resonances in the state-to-state reaction probabilities. The quasibound states, mostly the bending motions, give rise to strong reasonance peaks, whereas other motions contribute to the bumps and shoulders in the resonance structure. The initial state reaction probability further proves that the bending motions are the dominating factors of the reaction probability and have longer life times than the stretching motions. This is the first observation of reactive resonances from a "Lake Eyring" feature in a potential energy surface.

Wang, Dunyou↗

3D probabilistic fracture mechanics / computational fluid dynamics simulation of a reactor pressure vessel under transient conditions

Reactor pressure vessels (RPVs) are safety-critical light-water-reactor components that, under irradiation, experience long-term material degradation in the form of embrittlement. This can increase their susceptibility to fracture under thermal-shock conditions, which could occur during off-normal transients such as loss-of-coolant accidents (LOCAs). During a LOCA, the most severe conditions for the RPV occur when emergency core cooling water is injected through the cold legs into the water-and-steam-filled RPV. The rapid cooling of the downcomer and internal RPV surface causes decreased temperature and elevated thermally driven tensile stresses in the RPV wall. This, combined with long-term material embrittlement, may cause fracture initiation at pre-existing flaws, challenging the integrity of the RPV. Assessing RPV integrity during transients with a large spatial variation in the coolant temperature requires a modeling approach that considers the effects of spatially varying coolant temperature on the fracture probability of a population of flaws distributed throughout the RPV, accounting for spatially varying embrittlement. Here, the present study addresses this need by demonstrating first-of-its-kind coupling of a high-fidelity 3D computational fluid dynamics code with 3D probabilistic fracture mechanics This was accomplished using representative models of a pressurized-water reactor subjected to small- and medium-break LOCA conditions, both of which can result in large spatial temperature variations. While the observed impact of accounting for 3D effects was minimal under the small-break LOCA this study indicates a significant increase in the probability of fracture initiation under the medium-break LOCA when 3D effects are considered, relative to a spatially uniform cooling scenario.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Max-independent set and the quantum alternating operator ansatz

he maximum-independent set (MIS) problem of graph theory using the quantum alternating operator ansatz is studied. We perform simulations on the Rigetti Forest simulator for the square ring, K 2,3 , and K3,3 graphs and analyze the dependence of the algorithm on the depth of the circuit and initial states. The probability distribution of observation of the feasible states representing maximum-independent sets is observed to be asymmetric for the MIS problem, which is unlike the Max-Cut problem where the probability distribution of feasible states is symmetric. For asymmetric graphs, it is shown that the algorithm clearly favors the independent set with the larger number of elements even for finite circuit depth. Finally, we also compare the approximation ratios for the algorithm when we choose different initial states for the square ring graph and show that it is dependent on the choice of the initial state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Activated recombinative desorption: A potential component in mechanisms of spacecraft glow

The concept of activated recombination of atomic species on surfaces can explain the production of vibrationally and translationally excited desorbed molecular species. Equilibrium statistical mechanics predicts that the molecular quantum state distributions of desorbing molecules is a function of surface temperature only when the adsorption probability is unity and independent of initial collision conditions. In most cases, the adsorption probability is dependent upon initial conditions such as collision energy or internal quantum state distribution of impinging molecules. From detailed balance, such dynamical behavior is reflected in the internal quantum state distribution of the desorbing molecule. This concept, activated recombinative desorption, may offer a common thread in proposed mechanisms of spacecraft glow. Using molecular beam techniques and equipment available at Los Alamos, which includes a high translational energy 0-atom beam source, mass spectrometric detection of desorbed species, chemiluminescence/laser induced fluorescence detection of electronic and vibrationally excited reaction products, and Auger detection of surface adsorbed reaction products, a fundamental study of the gas surface chemistry underlying the glow process is proposed.

Cross, J. B.↗

Predicting Lightning Initiation using Deep Learning

Lightning occurrence presents safety challenges to people and property. The main challenge with lightning safety is that the majority of guidance is reactive. In other words, lightning has to have already occurred nearby before a person will respond and take shelter. Further, most injuries or fatalities occur as the storm approaches, or as it's moving away, when rainfall may not be present at the time of the flash. Thus, this project develops a physically-based deep learning model to produce lightning probabilities out to 15 minutes. The deep learning model combines a Convolutional Neural Network (CNN) with a Long Short-Term Memory (LSTM) network to capture both the spatial and temporal evolution of storms to predict the probability that lightning initiation will occur in the next 15 minutes. The model combines radar reflectivity, correlation coefficient and differential reflectivity to inferred storm hydrometer type and precipitation phase, which aids in the identification of electrification processes. The model is trained with data from the Geostationary Lightning Mapper (GLM), which is a near infrared sensor onboard the GOES-R series of satellites that measures optical brightness from lightning. This presentation will provide an overview of the project.

Andrew T White↗

The early ultraviolet spectral evolution of Nova Cygni 1992

We report on the high resolution ultraviolet line profile evolution of Nova Cygni 1992 using the International Ultraviolet Explorer (IUE) satellite and the Goddard High Resolution Spectrograph (GHRS) on the Hubble Space Telescope (HST). Our observations cover the period from discovery in 1992 February through 1993 April. The initial expansion velocity, derived from P Cyg profiles on Mg II and other resonance lines, was about 4500 km/s. All resonance transitions displayed P Cyg profiles around the time of UV maximum. We first discuss the spectral development as the initially optically thick atomic absorption curtain became optically thin. We then present the interpretation of the nebular spectral stage. The high resolution line profiles show that the shell rapidly developed nearly symmetric knots, or filaments, of emission. The optically thin transitions show an emission weighted expansion velocity of about 1800 km/s, consistent with the velocity inferred from the radio and infrared data. The emission lines display a symmetric set of filaments that appeared as soon as the line profiles became optically thin. The GHRS observations demonstrate that the filaments were located in the inner, slower moving parts of the ejecta. These probably record an initial low wave number instability imposed on the ejecta at the time of outburst. We suggest that the likely cause is a Rayleigh-Taylor instability. The emission at later stages is dominated by the filaments and, since this phenomenon is seen in virtually all novae, nebular line formation codes will have to account for the density inhomgeneities in order to adequately determine abundances for nove ejecta.

Shore, Steven N.↗

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

Initialization and Restart in Stochastic Local Search: Computing a Most Probable Explanation in Bayesian Networks

For hard computational problems, stochastic local search has proven to be a competitive approach to finding optimal or approximately optimal problem solutions. Two key research questions for stochastic local search algorithms are: Which algorithms are effective for initialization? When should the search process be restarted? In the present work we investigate these research questions in the context of approximate computation of most probable explanations (MPEs) in Bayesian networks (BNs). We introduce a novel approach, based on the Viterbi algorithm, to explanation initialization in BNs. While the Viterbi algorithm works on sequences and trees, our approach works on BNs with arbitrary topologies. We also give a novel formalization of stochastic local search, with focus on initialization and restart, using probability theory and mixture models. Experimentally, we apply our methods to the problem of MPE computation, using a stochastic local search algorithm known as Stochastic Greedy Search. By carefully optimizing both initialization and restart, we reduce the MPE search time for application BNs by several orders of magnitude compared to using uniform at random initialization without restart. On several BNs from applications, the performance of Stochastic Greedy Search is competitive with clique tree clustering, a state-of-the-art exact algorithm used for MPE computation in BNs.

Mengshoel, Ole J.↗

Structural and stratigraphic features and ERS 1 synthetic aperture radar backscatter characteristics of ice growing on shallow lakes in NW Alaska, winter 1991-1992

Changes in Earth Remote-Sensing Satellite (ERS) 1 C band synthetic aperture radar (SAR) backscatter intensity (sigma(exp 0)) from ice growing on shallow tundra lakes at three locations in NW Alaska are described. Ice core analysis shows that all lakes on the coast at Barrow the ice, whether floating or frozen to the bottom, includes an inclusion-free layer overlying a layer of ice with tubular bubbles oriented parallel to the direction of growth. The clear ice may also be overlain by a discontinuous layer of bubbly snow ice. Backscatter is low (-16 to -22 dB) at the time of initial ice formation, probably due to the specular nature of the upper and lower ice surfaces causing the radar pulse to be reflected away from the radar. As the ice thickens during the autumn, backscatter rises steadily. Once the ice freezes to the lake bottom, regardless of the presence of foward scattering tubular bubbles, low backscatter values of -17 to -18 dB are caused by absorption of the radar signal in the lake bed. For ice that remains afloat all winter the ice-water interface and the tubular bubbles combine, presumably via an incoherent double-bounce mechanism, to cause maximum backscatter values of the order of -6 to -7 dB. The sigma(exp 0) saturates at -6 to -7 dB before maximum ice thickness and tubular bubble content are attained. A simple ice growth model suggests that the layer of ice with tubular bubbles need be only a few centimeters thick midway through the growth season to cause maximum backscatter from floating ice. During the spring thaw a previously unreported backscatter reversal is observed on the floating and grounded portions of the coastal lakes but not on the lakes farther inland. This reversal may be related to the ice surface topography and wetness plus the effects of a longer, cooler melt period by the coast. Time series of backscatter variations from shallow tundra lakes are a record of (1) the development of tubular bubbles in the ice and, by association, changes in the gas content of the underlying water and (2) the freezing of ice to the bottoms of the lakes and therefore lake bathymetry and water availability. SAR is also able to detect the onset of lake ice growth in autumn and the initiation of the spring thaw and thus has potential for monitoring high-altitude lake ice growth and decay processes in relation to climate variability.

Jeffries, M. O.↗

Caught in headlights: Captive white-tailed deer responses to variations in vehicle lighting during imminent collision scenarios

Vehicle collisions with deer (Odocoileus spp.) cause billions of dollars in damages and injure thousands of drivers every year in the United States, and few mitigation methods have proven effective. However, recent research suggests that vehicle lighting might influence white-tailed deer (Odocoileus virginianus; hereafter, deer) responses to oncoming vehicles. Most new vehicles are manufactured with light emitting diode (LED) headlights which differ in total radiance and wavelength of light emitted compared to the previous industry standard of tungstenhalogen (halogen) headlights. Also, frontal vehicle illumination through rear-facing lighting has shown promise in enhancing deer responses to vehicles, but its effectiveness has not been tested under various headlight conditions (headlight type or intensity). As such, it remains unclear how these aspects of vehicle lighting affect deer responses to an approaching vehicle. We used 23 captive, wild-type deer to investigate how variations in vehicle lighting affect deer responses to an approaching vehicle at night, when most collisions occur. We released deer into a 95 m long, 3 m wide chute and approached them from the opposite end with an electric golf cart outfitted with two versions of stock 2017–2020 Ford Fusion headlights (LED and halogen) and a 51 cm rear-facing lightbar to test how vehicle lighting affected deer avoidance behaviors in an imminent, head-on collision scenario. Each deer received eight lighting treatments consisting of unique combinations of headlight type (LED vs. halogen), light intensity (low vs. high beam), and rear-facing lighting (lightbar off vs. on). We measured deer alert and flight behavior using infrared videography. We found that the halogen, high beam, lightbar off treatment had the greatest probability of evoking an alert response. Furthermore, when the lightbar was off, high beams appeared to increase alert probability for halogen headlights. Also, we found evidence that high beam, halogen headlights tend to increase alert probability over high beam, LED headlights, when the lighbar was off. We found no effect of our lighting treatments on deer alert distance, flight probability, or flight initiation distance. Across all behavioral responses, the random effect deer ID explained 0.86–9.19 × more variation than our lighting treatments, reflecting large differences in responses among deer. Overall, we found that vehicle lighting can impact deer behavior during an imminent, head-on collision scenario, although lighting was ineffective at increasing favorable flight behaviors. Future research should investigate how vehicle lighting treatments affect free-ranging, wild deer in a variety of real-world scenarios and at longer approach distances.

White-tailed deer (Odocoileus virginianus) Deer be↗

Control-Affine Schrödinger Bridge and Generalized Bohm Potential

From a stochastic control perspective, the Schrödinger bridge is a density-valued continuous curve parameterized by time that connects a given pair of initial and terminal probability densities via minimum effort controlled Brownian motion. The control-affine Schrödinger bridge extends this idea to a generic control-affine Itô diffusion, possibly with an additive state cost. Here, in this letter, we recast the necessary conditions of optimality for the control-affine Schrödinger bridge problem as a two point boundary value problem for a quantum mechanical Schrödinger PDE with complex potential. This complex-valued potential is a generalization of the real-valued Bohm potential in quantum mechanics. Our derived potential is akin to the optical potential in nuclear physics where the real part of the potential encodes elastic scattering (transmission of wave function), and the imaginary part encodes inelastic scattering (absorption of wave function). The key takeaway is that the process noise that drives the evolution of probability densities induces an absorbing medium in the evolution of wave function. These results make new connections between control theory and non-equilibrium statistical mechanics through the lens of quantum mechanics.

Markov processes↗

Reactor Pressure Vessel Fracture Mechanics Development and Concrete Application Testing for Grizzly

The Grizzly code is being developed to address degradation issues in nuclear reactor structures and components. For light-water reactors, Grizzly currently has capabilities to simulate degradation processes and their effects on structural integrity in two key areas: reactor pressure vessels (RPVs) and reinforced concrete structures. This report documents improvements made to Grizzly’s ability to address both of these structural systems. For RPVs, the reduced-order models (ROMs) used in fracture mechanics calculations have been expanded to allow their application over a broader range of the parameter space than was permitted by the previous models. The ROMs currently used in Grizzly for the evaluation of flaws that are fully embedded within the RPV (as opposed to surface-breaking flaws) are based on a model that is known to be conservative, indicating higher stress intensity factors than would be obtained from direct simulations. A more accurate model that eliminates these excess conservatisms has been recently included in the American Society of Mechanical Engineers Boiler and Pressure Vessel Code but was not applicable for flaws near the RPV surface, which is where the most critical flaws are usually located. That model has recently been extended for increased applicability in this near-surface region. The ROMs for embedded flaws in the Grizzly code have been expanded to include these recent extensions, which permit their use in a much broader set of cases than previously possible. Direct 3D simulations have been used to check these ROMs and have shown good agreement in most cases, although there are still some cases that need further investigation. There are considerable benefits to using these these more accurate and less conservative ROMs for embedded flaws. On a benchmark probabilistic fracture mechanics problem tested here, the conditional probability of fracture initiation computed for a population of flaws in a single plate in an RPV decreased by over a factor of 3. To address aging in reinforced concrete structures, a capability to simulate multiple degradation mechanisms, including alkali-silica reaction and radiation-induced volumetric expansion has been developed in Grizzly over the past several years. This had previously been demonstrated on laboratory-scale specimens but not on full-scale nuclear concrete structures with reinforcement. To demonstrate the applicability of Grizzly to the analysis of large-scale structures of interest, a full 3D model of a representative reinforced concrete structure, including a complex arrangement of reinforcing bars, was developed and demonstrated in Grizzly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials and Components. Third Annual Progress Report

Additive manufacturing (AM) of high-strength corrosion resistance alloys for nuclear energy applications, such as stainless steel and Inconel, is currently based on laser powder bed fusion (LPBF) process. Some of the challenges with using LPBF method for nuclear manufacturing include the possibility of introducing pores into metallic structures. Probability of crack initiation at the pore depends on size, shape, and orientation of the defect. Pulsed Infrared Thermography Imaging (PIT) provides a capability for non-destructive evaluation (NDE) of sub-surface defects in arbitrary size structures. The PIT method is based on recording material surface temperature transients with infrared (IR) camera following thermal pulse delivered on material surface with flash light. The PIT method has advantages for NDE of actual AM structures because the method involves one-sided non-contact measurements and fast processing of large sample areas captured in one image. Following initial qualification of an AM component for deployment in a nuclear reactor, a PIT system can also be used for in-service nondestructive evaluation (NDE) applications. In this report, we describe recent progress in enhancing PIT capabilities in detecting microscopic subsurface defects in metals, and classifying shapes and orientation of pores in thermal images. For detection of microscopic defects in PIT imaging data, we have developed Spatial Temporal Denoised Thermal Source Separation (STDTSS) unsupervised machine learning (ML) image processing algorithm. We show that flat bottom hole (FBH) defects as small as 200µm in SS316 and IN718 specimens, can be detected with STDTSS algorithm. To the best of our knowledge, these are the smallest detected defects which are reported in literature. For classification of defects shapes, we have previously developed thermal tomography (TT) algorithm to obtain depth reconstructions of material defects from data cube of sequentially recorded surface temperatures. However, interpretation of TT images is non-trivial because of blurring with increasing depth. To address this challenge, we have developed a deep learning convolutional neural network (CNN) to classify size and orientation subsurface defects in simulated TT images.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Probabilistic Fracture Mechanics and Optimum Fracture Control Analytical Procedures for a Reusable Solid Rocket Motor Case

A methodology for the reliability analysis of a reusable solid rocket motor case is discussed. The analysis is based on probabilistic fracture mechanics and probability distribution for initial flaw sizes. The developed reliability analysis is used to select the structural design variables of the solid rocket motor case on the basis of minimum expected cost and specified reliability bounds during the projected design life of the case. Effects of failure prevention plans such as nondestructive inspection and the material erosion between missions are also considered in the developed procedure for selection of design variables. The reliability-based procedure can be modified to consider other similar structures of reusable space vehicle systems with different failure prevention plans.

Hanagud, S.↗

Solution to a gene divergence problem under arbitrary stable nucleotide transition probabilities

A nucleic acid chain, L nucleotides in length, with the specific base sequence B(1)B(2) ... B(L) is defined by the L-dimensional vector B = (B(1), B(2), ..., B(L)). For twelve given constant non-negative transition probabilities that, in a specified position, the base B is replaced by the base B' in a single step, an exact analytical expression is derived for the probability that the position goes from base B to B' in X steps. Assuming that each base mutates independently of the others, an exact expression is derived for the probability that the initial gene sequence B goes to a sequence B' = (B'(1), B'(2), ..., B'(L)) after X = (X(1), X(2), ..., X(L)) base replacements. The resulting equations allow a more precise accounting for the effects of Darwinian natural selection in molecular evolution than does the idealized (biologically less accurate) assumption that each of the four nucleotides is equally likely to mutate to and be fixed as one of the other three. Illustrative applications of the theory to some problems of biological evolution are given.

Holmquist, R.↗

Reusable solid rocket motor case - Optimum probabilistic fracture control

A methodology for the reliability analysis of a reusable solid rocket motor case is discussed in this paper. The analysis is based on probabilistic fracture mechanics and probability distribution for initial flaw sizes. The developed reliability analysis can be used to select the structural design variables of the solid rocket motor case on the basis of minimum expected cost and specified reliability bounds during the projected design life of the case. Effects on failure prevention plans such as nondestructive inspection and the material erosion between missions can also be considered in the developed procedure for selection of design variables. The reliability-based procedure that has been discussed in this paper can easily be modified to consider other similar structures of reusable space vehicle systems with different fracture control plans.

Hanagud, S.↗

Schiaparelli Basin, Mars - Morphology, Tectonics and infilling history

Schiaparelli Basin (an elliptical basin 460 x 400 km in diameter) is the best preserved large impact structure on Mars. Multiple periods of basin infilling by lava and eolian material can be recognized. Numerous ridges within the basin show a strong northwest-southeast orientation, consistent with a regional stress field induced by the loading of the lithosphere by Tharsis. Ridge location is, however, apparently dominated by a buried peak ring. Earth-based radar measurements of Schiaparelli indicate that the ejecta deposits are best preserved on topographic highs southeast of the basin, while eolian materials have preferentially accumulated in the lower (western) side of the basin. Lyot (215 km diameter) is believed to be a good analogue to the pre-infilled Schiaparelli, and indicates that Schiaparelli probably had an initially hummocky floor topography and both ballistic and ground-flow ejecta deposits.

Mouginis-Mark, P. J.↗