Search NASA⌕ Search

SEARCH · Search NASA

Results for “stochastic control”

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 847 records · Page 47

Improving Balance Function Using Low Levels of Electrical Stimulation of the Balance Organs

Crewmembers returning from long-duration space flight face significant challenges due to the microgravity-induced inappropriate adaptations in balance/ sensorimotor function. The Neuroscience Laboratory at JSC is developing a method based on stochastic resonance to enhance the brain s ability to detect signals from the balance organs of the inner ear and use them for rapid improvement in balance skill, especially when combined with balance training exercises. This method involves a stimulus delivery system that is wearable/portable providing imperceptible electrical stimulation to the balance organs of the human body. Stochastic resonance (SR) is a phenomenon whereby the response of a nonlinear system to a weak periodic input signal is optimized by the presence of a particular non-zero level of noise. This phenomenon of SR is based on the concept of maximizing the flow of information through a system by a non-zero level of noise. Application of imperceptible SR noise coupled with sensory input in humans has been shown to improve motor, cardiovascular, visual, hearing, and balance functions. SR increases contrast sensitivity and luminance detection; lowers the absolute threshold for tone detection in normal hearing individuals; improves homeostatic function in the human blood pressure regulatory system; improves noise-enhanced muscle spindle function; and improves detection of weak tactile stimuli using mechanical or electrical stimulation. SR noise has been shown to improve postural control when applied as mechanical noise to the soles of the feet, or when applied as electrical noise at the knee and to the back muscles.

Bloomberg, Jacob↗

MHD-induced SOL filaments and divertor heat flux striations in NSTX

Magnetohydrodynamic (MHD) activity in NSTX has been found to induce striated heat flux patterns at the divertor. These striated heat fluxes were previously interpreted as strike-point splitting caused by stochastic magnetic lobes near the X-point. Recent analysis indicates that MHD can directly generate filaments in the scrape-off layer (SOL), and the MHD-related filaments are consistent with the rotation and the frequency of the MHD. Magnetic field line tracing from the mid-plane of the MHD-related filaments to the divertor is consistent with the location of striated heat fluxes on the divertor. Additionally, gas puff imaging (GPI) diagnostic data shows the generation of regular blobs following each edge harmonic oscillation (EHO) toroidal mode number n = 1 perturbation, which indicates that the MHD perturbation directly drives the blobs. These blobs are responsible for the divertor heat flux striations. Here, this finding provides a new framework for understanding MHD-induced edge transport and has important implications for heat-load control in NSTX-U and future spherical tokamaks.

MHD↗

Cosmogenic 3 He exposure dating in mafic rocks by “Virtual mineral separation” of pyroxene

We describe a “virtual mineral separation” method for measuring the cosmogenic 3 He concentration in pyroxene in mafic rocks that consist mainly of plagioclase and pyroxene, without physically separating the minerals. This approach is significantly faster and more cost-effective than the conventional method, which requires physical separation and purification of pyroxene grains by time-consuming and labor-intensive crushing, acid cleaning, magnetic separation, HF etching, and handpicking under a microscope. The premise of the method is that helium diffusivity is much higher in plagioclase than in pyroxene, so controlled preheating of a mixed whole-rock sample can degas 3 He from plagioclase while retaining all 3 He in pyroxene. A second heating step releases all 3 He from pyroxene for measurement. To then obtain a 3 He concentration in pyroxene rather than the whole rock, we determine the pyroxene weight fraction in the sample using X-ray computed tomography (CT). A comparison of 3 He concentrations in pyroxene measured using virtual mineral separation with those measured in the same samples by physical mineral separation in previous work shows no evidence of systematic bias between the methods. Virtual mineral separation greatly simplifies the workflow for 3 He exposure-dating of mafic rocks, reduces time, effort, and cost, and permits measurements on very small samples. This enables new emerging applications of exposure dating, such as quantifying stochastic surface processes, ecosystem studies, and potential subglacial bedrock exposure dating.

Geosciences↗

Tailoring resonant magnetic perturbation to optimize fast-ion confinement during ELM control in KSTAR

Abstract 3D resonant magnetic perturbation (RMP) is one promising way to control edge localized modes that can cause excessive material erosion of tokamak first walls. However, RMP can lead to undesired degradation of plasma confinement, including fast-particle losses, which can impact the performance and safety of the reactor. This work investigates the optimization of the poloidal spectrum of the 3D field to optimize fast ion confinement during edge localized mode (ELM) suppression. In the initial step, the validity of the modeling framework is tested against experimental data. Simulations successfully replicate an increase in poloidal limiter temperature with different poloidal spectra. Then, the simulation shows improvement of fast ion confinement with a reduction of core resonant response, while edge resonant magnetic fields are maintained above the threshold to sustain the ELM suppression. Reduction of the core resonant fields keeps the Kolmogorov–Arnold–Moser surface and reduces the fast particle losses due to the stochastic magnetic field lines. The results highlight the potential of edge localization of the resonant fields to enhance the performance of fusion reactors, but further investigation is needed to improve the validation of this approach.

3D field↗

From Coherence to Function: Exploring the Connection in Chemical Systems

The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories, such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate, the pathway from reactants to products, is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic–vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photoinduced electron transfer, singlet–triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into the rapid decoherence, or loss of phase and amplitude correlations, of vibrational coherences along promoter vibrations during subpicosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from the Franck–Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics rather than by the laser pulse in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine–molybdenum complex hosting a nonreactive dinitrogen substrate. In conclusion, a slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck–Condon active terpyridine vibration to a dinitrogen stretching vibration, energizing the dinitrogen bond.

Electrical energy↗

Functional Traits Resolve Mechanisms Governing the Assembly and Distribution of Nitrogen-Cycling Microbial Communities in the Global Ocean

Microorganisms drive much of the marine nitrogen (N) cycle, which jointly controls the primary production in the global ocean. However, our understanding of the microbial communities driving the global ocean N cycle remains fragmented. Focusing on “who is doing what, where, and how?”, this study draws a clear picture describing the global biogeography of marine N-cycling microbial communities by utilizing the Tara Oceans shotgun metagenomes. The marine N-cycling communities are highly variable taxonomically but relatively even at the functional trait level, showing clear functional redundancy properties. The functional traits and taxonomic groups are shaped by the same set of geo-environmental factors, among which, depth is the major factor impacting marine N-cycling communities, differentiating mesopelagic from epipelagic communities. Latitudinal diversity gradients and distance-decay relationships are observed for taxonomic groups, but rarely or weakly for functional traits. The composition of functional traits is strongly deterministic as revealed by null model analysis, while a higher degree of stochasticity is observed for taxonomic composition. Integrating multiple lines of evidence, in addition to drawing a biogeographic picture of marine N-cycling communities, this study also demonstrated an essential microbial ecological theory—determinism governs the assembly of microbial communities performing essential biogeochemical processes; the environment selects functional traits rather than taxonomic groups; functional redundancy underlies stochastic taxonomic community assembly.

59 BASIC BIOLOGICAL SCIENCES↗

Core Design of the Holos-Quad Microreactor

The Holos-Quad micro-reactor concept, developed by HolosGen LLC, is equipped with a 22 MWt (Mega-Watt thermal) core and an integral power conversion system converting the core thermal energy into approximately 10 MWe (Mega-Watt electric). This design can be configured to support a wide range of applications. It is a very innovative high-temperature gas-cooled reactor concept using TRI-structural ISOtropic particle fuel (TRISO) distributed in graphite hexagonal blocks, cooled with helium in a direct Brayton cycle independently executed by four Subcritical Power Modules (SPMs) fitted into a hardened 40-foot container whose dimensions are in compliance with ISO shipping containers requirements. In FY2019 HolosGen LLC was awarded by the Department of Energy Advanced Research Project Agency-Energy (DOE ARPA-E) under the MEITNER funding program. As part of the MEITNER award, the Argonne National Laboratory (ANL) contributed expertise through two specialized teams: The “Design Team” and the “Resource Team”. The Design Team was dedicated to validate feasibility of the Holos-Quad core and to optimize its core design through neutronics analyses. The Resource Team was dedicated to feasibility verification via high-fidelity codes of Holos-Quad thermal-hydraulic, heat transfer, shielding, and structural aspects. This report summarizes the activities conducted by ANL Design Team. A rigorous design approach based on multi-criteria optimization and code-to-code comparison involving stochastic and high-fidelity deterministic solutions was developed and employed at several evolutionary stages of the Holos-Quad design. Several generations of the Holos-Quad core were designed within this project before converging to the current full-scale Gen 2+ design that is detailed in this report. Figure EA-1 illustrates a cross-sectional view of Gen 2+ Holos-Quad core configuration, and Figure EA-2 provides a simplified perspective view of 1-of-4 SPMs. The Holos-Quad uses four thermal-hydraulically independent SPMs locked into stationary positions during power operation, surrounded by BeO reflector and structural component fully comprised within the dimensional constraints represented by traditional ISO containers. One of the benefits of this approach is to enable transportation of each SPM promptly after irradiation in shielded containers. The core is designed to operate for approximately 8 full-power years while the reactivity controls and power conversion system enable load-following operations. The reactivity controls are represented by independent, diversified, and redundant reactivity control systems based on control drums and redundant sets of shutdown rods. The high-fidelity simulation tools were used to assess detailed power and flux distributions of the three-dimensional full-core or quarter-core of the Gen 2+ configuration. Single-physics and multi-physics simulations of the neutronics code PROTEUS and the thermal-hydraulic code System Analysis Module (SAM) were performed to analyze the Holos design configurations with detailed high-fidelity solutions. The design work performed confirmed feasibility of the Holos-Quad concept, provided realistic design description for detailed design of the operational system, and identified several core design improvements to be further considered for future reactor development activities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Turbulence and the Stabilization Principle

Further results of research, reported in several previous NASA Tech Briefs articles, were obtained on a mathematical formalism for postinstability motions of a dynamical system characterized by exponential divergences of trajectories leading to chaos (including turbulence). To recapitulate: Fictitious control forces are introduced to couple the dynamical equations with a Liouville equation that describes the evolution of the probability density of errors in initial conditions. These forces create a powerful terminal attractor in probability space that corresponds to occurrence of a target trajectory with probability one. The effect in ordinary perceived three-dimensional space is to suppress exponential divergences of neighboring trajectories without affecting the target trajectory. Con sequently, the postinstability motion is represented by a set of functions describing the evolution of such statistical quantities as expectations and higher moments, and this representation is stable. The previously reported findings are analyzed from the perspective of the authors Stabilization Principle, according to which (1) stability is recognized as an attribute of mathematical formalism rather than of underlying physics and (2) a dynamical system that appears unstable when modeled by differentiable functions only can be rendered stable by modifying the dynamical equations to incorporate intrinsic stochasticity.

Zak, Michail↗

Lower hybrid waves excited through linear mode coupling and the heating of ions in the auroral and subauroral magnetosphere

Two related aspects of the linear mode conversion mechanism are considered: (1) the controlled heating of suprathermal ions in the ionosphere and magnetosphere over a powerful VLF/ELF transmitter using lower hybrid (LH) waves excited by the electromagnetic (EM) transmitter signals through linear mode conversion; (2) the excitation of intense LH waves in the auroral regions through the linear conversion of VLF/ELF EM auroral hiss, and the subsequent heating of ions by the excited LH waves. Experimental results are used to demonstrate that strong LH waves can be excited by EM waves down to frequencies as low as 2 kHz in the subauroral low-altitude magnetosphere. It is also shown that an approximately 250 kW VLF/ELF transmitter operating in the subauroral region at 2 kHz could excite sufficiently intense LH waves to heat suprathermal 1 eV H(+) ions and 16 eV O(+) ions to roughly 50 eV in a region extending in altitude from 1000 to 5000 km above the transmitter with a horizontal scale of about 500 km. If coherent wave stochastic heating occurs, the O(+) energy gain could be as large as 200 eV.

Bell, T. F.↗

Efficient Distance-based Global Sensitivity Analysis for Terrestrial Ecosystem Modeling

Sensitivity analysis in terrestrial ecosystem modeling is important for understanding controlling processes, guiding model development, and targeting new observations to reduce parameter and prediction uncertainty. Complex and computationally expensive terrestrial ecosystem models (TEM) limit the number of ensemble simulations, requiring sophisticated and efficient methods to analyze sensitivities of multiple model responses to different types of parameter uncertainties. In this study, we propose a distance-based global sensitivity analysis (DGSA) method. DGSA first classifies model response samples into a small set of discrete classes and then calculates the distance between parameter frequency distributions in different classes to measure the parameter sensitivity. The principle is that, if the parameter distribution is the same in each class, then the model response is insensitive to the parameter, while a large difference in the distributions indicates the parameter is influential to the response. Built on this idea, DGSA can be applied to analyze sensitivity of a single and a group of responses to different kinds of parameter uncertainties including continuous, discrete and even stochastic. Besides the main-effect sensitivity from a single parameter, DGSA can also quantify the sensitivity from parameter interactions. Additionally, DGSA is computationally efficient which can use a small number of model evaluations to obtain an accurate and statistically significant result. We applied DGSA to two TEMs, one having eight parameters and three kinds of model responses, and the other having 47 parameters and a long-period response. We demonstrated that DGSA can be used for sensitivity problems with multiple responses and high-dimensional parameters efficiently.

Lu, Dan↗

Stabilization of the 81-channel coherent beam combination using machine learning

We develop a rapidly converging algorithm for stabilizing a large channel-count diffractive optical coherent beam combination. An 81-beam combiner is controlled by a novel, machine-learning based, iterative method to correct the optical phases, operating on an experimentally calibrated numerical model. A neural-network is trained to detect phase errors based on interference pattern recognition of uncombined beams adjacent to the combined one. Due to the non-uniqueness of solutions in the full space of possible phases, the network is trained within a limited phase perturbation/error range. This also reduces the number of samples needed for training. Simulations have proven that the network can converge in one step for small phase perturbations. When the trained neural-network is applied to a realistic case of 360 degree full range, an iterative scheme exploits random walking at the beginning, with the accuracy of prediction on phase feedback direction, to allow the neural-network to step into the training range for fast convergence. This neural-network-based iterative method of phase detection works tens of times faster than the commonly used stochastic parallel gradient descent approach (SPGD) using a single-detector and random dither when both are tested with random phase perturbations.

Wang, Dan↗

Nonequilibrium steady state and heat transport in nonlinear open quantum systems: Stochastic influence action and functional perturbative analysis

In this paper, we show that a nonequilibrium steady state (NESS) exists at late times in open quantum systems with weak nonlinearity by following its nonequilibrium dynamics with a perturbative analysis. We consider an oscillator chain containing three-types of anharmonicity: cubic α- and quartic β-type Fermi–Pasta–Ulam–Tsingou (FPUT) nearest-oscillator interactions and the on-site (pinned) Klein–Gordon (KG) quartic self-interaction. Assuming weak nonlinearity, we introduce a stochastic influence action approach to the problem and obtain the energy flows in different junctures across the chain. The formal results obtained here can be used for quantum transport problems in weakly nonlinear quantum systems. For α-type anharmonicity, we observe that the first-order corrections do not play any role in the thermal transport in the NESS of the configuration we considered. For KG and β-types anharmonicity, we work out explicitly the case of two weakly nonlinearly coupled oscillators, with results scalable to any number of oscillators. We examine the late-time energy flows from one thermal bath to the other via the coupled oscillators, and show that both the zeroth- and the first-order contributions of the energy flows become constant in time at late times, signaling the existence of a late-time NESS to first order in nonlinearity. Our perturbative calculations provide a measure of the strength of nonlinearity for nonlinear open quantum systems, which may help control the mesoscopic heat transport distinct from or close to linear transport. Furthermore, our results also give a benchmark for the numerical challenge of simulating heat transport. Our setup and predictions can be implemented and verified by investigating heat flow in an array of Josephson junctions in the limit of large Josephson energy with the platform of circuit QED.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tracking oxygen vacancy migration in memristor devices using operando hard X-ray photoelectron spectroscopy

Fabrication of ultrathin (sub-2 nm) oxide semiconductor memristors poses the fundamental challenge of achieving oxide growth with atomic precision in terms of electronic structure and defect formation. Recently, ultrathin memristors consisting of bilayers of mixed Al 2 O 3 and MgO atomic layers were fabricated using an in vacuo atomic layer deposition process. This approach offers a unique platform for precise atomic control of oxygen vacancies in the device in which the vacancies are introduced by placing MgO atomic layers between pristine Al 2 O 3 layers. In this work, we present a systematic operando Hard X-ray PhotoElectron Spectroscopy (HAXPES) study of the switching of such memristors, combined with complementary current-voltage and capacitance-voltage (C-V) measurements. We used a memristor stack of Al 2 O 3 /MgO atomic layers, with the MgO-containing oxide deposited on the bottom Al metal electrode and a pure Al 2 O 3 layer below the top Pd electrode. HAXPES analysis shows a substantial change in the chemical shift of the Aluminum oxide when switching between the ”OFF” and ”ON” states indicative of a redistribution of oxygen vacancies in the device active layer. Interestingly, subsequent switching to the OFF state shows hysteretic behavior indicating the retention of some oxygen vacancies in the top Al 2 O 3 layer. This vacancy retention can be correlated with the stochastic behavior of the switching voltage observed in these devices. C–V measurements show a clear frequency-dependent response in the OFF state, consistent with enhanced polarization and vacancy trapping at low frequencies.

47 OTHER INSTRUMENTATION↗

A Q-Learning-Based Approach for Distributed Beam Scheduling in mmWave Networks

We consider the problem of distributed downlink beam scheduling and power allocation for millimeter-Wave (mmWave) cellular networks where multiple base stations (BSs) belonging to different service operators share the same unlicensed spectrum with no centralized control or explicit coordination among them. Our goal is to design efficient distributed beam scheduling and power allocation schemes such that the network-level payoff, defined as the weighted sum of the total throughput and a power penalization term, can be maximized. To this end, we propose a distributed scheduling approach to power allocation and adaptation for efficient interference management over the shared spectrum by modeling each BS as an independent Q-learning agent. As a baseline, we compare the proposed approach to the state-of-the-art non-cooperative game-based approach which was previously developed for the same problem. We conduct extensive experiments under various scenarios to verify the effect of multiple factors on the performance of both approaches. Experiment results show that the proposed approach adapts well to different interference conditions by learning from experience and thus can achieve higher payoff than the game-based approach. The proposed approach can also be integrated into our previously developed Lyapunov stochastic optimization framework for the purpose of network utility maximization with optimality guarantee. As a result, the weights in the payoff function can be automatically and optimally determined by the virtual queue values from the sub-problems derived from the Lyapunov optimization framework.

99 GENERAL AND MISCELLANEOUS↗

Tectonic-Climate Interactions in Action Orogenic Belts: Quantification of Dynamic Topography with SRTM data

This project was undertaken to examine the approach to steady state in collisional mountain belts. Although the primary thrust of this grant was to look at larger collisional mountain belts, such as the Himalaya, the Tien Shan, and Southern Alps, we began by looking at smaller structures represented by growing and propagating folds. Like ranges that are evolving toward a topographic steady state, these folds undergo a series of morphologic changes as they are progressively uplifted and eroded. We wanted to document the nature of these changes and to try to discern some of the underlying controls on them. We initially focused on the Wheeler Ridge anticline in southern California. Subsequently, we progressed to looking at the topographic development and the effects of differential uplift and glaciation on the Kyrgyz Range in the northern Tien Shan. This range is unusual inasmuch as it is transformed along its length from a simple uplift with a largely preserved Mesozoic erosion surface arching across it to a highly dissected and heavily glaciated uplift in the region where uplift has been sustained at higher rates over longer intervals. In efforts to understand the distribution of erosion rates at 10(exp 3) - 10(exp 5) year time scales, cosmogenic radionuclide (CRN) concentrations have been gaining increasingly widespread usage (Brown et al., 1995; Riebe et al., 2004; Riebe et al., 2001; Vance et al., 2003). Most studies to date, however, have been conducted in slowly eroding ranges. In rapidly eroding mountains where landslides deliver most of the sediments to the rivers, we hypothesized that CRN concentrations could be highly perturbed by the stochastic processes of landsliding. Therefore, we undertook the development of a numerical model that simulated the effects of both landsliding and grain-by-grain attrition within fluvial catchments. This modeling effort has shown the effects of catchment size and erosion rate on CRN concentrations and allows a prediction of where to sample to obtain the optimal erosion rate estimates using CRN techniques. Finally, we developed computational techniques to operate on DEMs to extract useful information that would enable quantification of climate-erosion interactions. In particular, we worked on rapid techniques to define catchments of any given range of sizes, to extract channel gradients, to combine precipitation information to calculate discharge, and to utilize various stream-power models to determine the erosional energy within any given catchment within a transect. We briefly describe results from Wheeler Ridge, the Kyrgyz Range, the Nepal Himalaya, and our numerical modeling.

Burbank, Douglas W.↗

Relationships Among Intrinsic Properties of Ordinary Chondrites: Oxidation State, Bulk Chemistry, Oxygen-isotopic Composition, Petrologic Type, and Chondrule Size

The properties of ordinary chondrites (OC) reflect both nebular and asteroidal processes. OC are modeled here as having acquired nebular water, probably contained within phyllosilicates, during agglomeration. This component had high Ai70 and acted like an oxidizing agent during thermal metamorphism. The nebular origin of this component is consistent with negative correlations in H, L, and LL chondrites between oxidation state (represented by olivine Fa) and bulk concentration ratios of elements involved in the metal-silicate fractionation (e.g., NdSi, Ir/Si, Ir/Mn, Ir/Cr, Ir/Mg, Ni/Mg, As/Mg, Ga/Mg). LL chondrites acquired the greatest abundance of phyllosilicates with high (delta)O-17 among OC (and thus became the most oxidized group and the one with the heaviest O isotopes); H chondrites acquired the lowest abundance, becoming the most reduced OC group with the lightest O isotopes. Chondrule precursors may have grown larger and more ferroan with time in each OC agglomeration zone. Nebular turbulence may have controlled the sizes of chondrule precursors. H-chondrite chondrules (which are the smallest among OC) formed from the smallest precursors. In each OC region, low-FeO chondrules formed before high-FeO chondrules during repeated episodes of chondrule formation. During thermal metamorphism, phyllosilicates were dehydrated; the liberated water oxidized metallic Fe-Ni. This caused correlated changes with petrologic type including decreases in the modal abundance of metal, increases in olivine Fa and low-Ca pyroxene Fs, increases in the olivine/pyroxene ratio, and increases in the kamacite Co and Ni contents. As water (with its heavy 0 isotopes) was lost during metamorphism, inverse correlations between bulk (delta)O-18 and bulk (delta)O-17 with petrologic type were produced. The H5 chondrites that were ejected from their parent body approx.7.5 Ma ago during a major impact event probably had been within a few kilometers of each other since they accreted approx.4.5 Ga ago. There are significant differences in the olivine compositional distributions among these rocks; these reflect stochastic nebular sampling of the oxidant (Le., phyllosilicates with high (delta)O-17) on a 0.1-1 km scale during agglomeration.

Rubin, Alan E.↗

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids—A case study in diamond

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater–Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a three-dimensional APOLLO3 neutrons deterministic scheme for the CABRI reactor

CABRI is an experimental reactor to study the fuel behavior during reactivity injection transients. These transients being highly multiphysics, the development of suitable modeling and simulation tools to simulate them is important for the optimization of the tests and the control of the experimental conditions. This paper focuses on the development of an APOLLO3 deterministic core calculation tool dedicated to the CABRI transient analysis. It represents the first stage of the incremental process for the implementation of a multiphysics time-dependent modeling of the CABRI transient. The neutron calculation scheme is based on a classical two-step approach. The first step consists of a 281-energy group calculation flux with the TDT-MOC (Method Of Characteristics) solver for cross-section space and energy (23 groups) collapsing for the CABRI different assembly clusters. The bias on a 2D core neutron calculation due to the self-shielding calculation and collapsing on a restricted pattern are investigated thanks to a comparison with a direct full 2D calculation on a quarter of core. The second step relies on a pin-resolved transport 3D transport core calculation with the SN solver MINARET. A progressive numerical validation process is followed to quantify the calculation biases on reactivity and reaction rates at each step using reference calculations with the stochastic code TRIPOLI4. The next development stage toward a multiphysics scheme will be the implementation of the 3D-kinetics equation resolution and the coupling with a core thermal-hydraulics model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗