Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “topological data analysis”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantum textures of the many-body wavefunctions in magic-angle graphene

Interactions among electrons create novel many-body quantum phases of matter with wavefunctions that often reflect electronic correlation effects, broken symmetries, and novel collective excitations. A wide range of quantum phases has been discovered in MATBG, including correlated insulating, unconventional superconducting, and magnetic topological phases. The lack of microscopic information, including precise knowledge of possible broken symmetries, has thus far hampered our understanding of MATBG's correlated phases. Here we use high-resolution scanning tunneling microscopy to directly probe the wavefunctions of the correlated phases in MATBG. The squares of the wavefunctions of gapped phases, including those of the correlated insulators, pseudogap, and superconducting phases, show distinct patterns of broken symmetry with a √3 x √3 super-periodicity on the graphene atomic lattice that has a complex spatial dependence on the moiré superlattice scale. We introduce a symmetry-based analysis to describe our measurements of the wavefunctions of MATBG's correlated phases with a set of complex-valued local order parameters. For the correlated insulators in MATBG, at fillings of ±2 electrons per moiré unit cell relative to charge neutrality, we compare the observed quantum textures to those expected for proposed theoretical ground states. In typical MATBG devices, the textures of correlated insulators' wavefunctions closely match those of the theoretically proposed IKS order, while in ultra-low-strain samples our data has local symmetries like those of a T-IVC phase. Here, we also study the wavefunction of MATBG's superconducting state, revealing strong signatures of intervalley coherence that can only be distinguished from those of the insulator with our phase-sensitive measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

p +p - electroproduction off protons in deuterium: measurements of quasi-free cross sections and kinematic probing of final state interactions

In this study, the process of pi+pi- electroproduction off protons bound in deuterium nuclei is explored. The exploration is performed through the analysis of experimental data on electron scattering off the deuteron target, collected in Hall B at Jefferson Lab with the CLAS detector. As a main result, the set of integrated and single-differential cross sections of the reaction yvp(n) -> p'(n')pi+pi- is obtained. The cross sections are extracted in the quasi-free regime, which implies that only events not affected by final state interactions are subject to selection. The measurements are performed in the kinematic region of the invariant mass W from 1.3 GeV to 1.825 GeV and photon virtuality Q2 from 0.4 GeV2 to 1 GeV2. Sufficient experimental statistics allows narrow binning, i.e. 25 MeV in W and 0.05 GeV2 in Q2, while maintaining an adequate statistical uncertainty. The extraction of quasi-free cross sections is accompanied by the kinematic probing of FSI between the reaction final hadrons and the spectator neutron in the aforementioned exclusive channel. In this probing the distributions of missing quantities are used in order to investigate the relative spread of events with FSI along the reaction phase space, trace the difference of FSI manifestations in different reaction topologies, reveal details on alterations of the hadron momentum in FSI, and isolate FSI contributions of various final hadrons. The performed examination is also capable of retrieving information on some underlying FSI mechanisms, among which the process of resonance formation in the intermediate state of pion-neutron interactions is particularly remarkable.

Skorodumina, Iuliia↗
DOE OSTI2017

Four coordination polymers based on 5-tert-butyl isophthalic acid and rigid bis(imidazol-1yl)benzene linkers: Synthesis, luminescence detection of acetone and optical properties

Four coordination polymers including, [Co(µ-Htbip){sub 2}(µ-dib)]{sub n} (1), [Co(µ-tbip)(µ-dmib){sub 0.5}]{sub n} (2), [Zn{sub 2}(µ-tbip)(µ{sub 3}-tbip)(µ-dmib){sub 1.5}]{sub n} (3) and [Cd(µ{sub 3}-tbip)(µ-dib){sub 0.5} (H{sub 2}O)]{sub n} (4) (tbip: 5-tert-butylisophthalate, dib: 1,4-bis(imidazol-1yl)benzene, dmib: 1,4-bis(imidazol-1yl)-2,5-dimethylbenzene), were hydrothermally synthesized and characterized by elemental analysis, IR spectra, single crystal and powder X-ray diffraction and thermal analysis (TG/DTA). The structural diversity is observed depending on ligands and coordination number of metal centers in the synthesized complexes. The tbip ligand displayed five different coordination modes in its complexes. In 1 and 2, complex 1 is 3D framework with the dia topology while complex 2 has 2D structure with the sql topology depending on coordination geometries of Co ions. Complex 3 is 3D framework with the fsh 4,6-conn topology and complex 4 has 2D 4-connected sql topology. Photoluminescent properties of complex 3 dispersed in various organic solvents were investigated and the results showed that 3 dispersed in methanol could be used as a fluorescent sensor for the detection of acetone. Moreover, thermal and optical properties of the complexes were also studied. - Graphical abstract: Four coordination polymers were hydrothermally synthesized and characterized by various techniques. The complexes showed the structural diversity depending on ligands and coordination number of metal centers. The tbip ligand displayed four different coordination modes in its complexes. In 1 and 2, complexes 1 and 2 are 3D and 2D structures with the dia and sql topologies depending on coordination geometries of Co ions, respectively. Complexes 3 and 4 are 3D and 2D structures with the fsh 4,6-conn and sql topology, respectively. Photoluminescent properties of complex 3 dispersed in various organic solvents were investigated and the results showed that 3 dispersed in methanol could be used as a fluorescent sensor for the detection of acetone. Moreover, thermal and optical properties of the complexes were also studied. - Highlights: • Four new 2D and 3D coordination polymers with 5-tert-butyl isophthalic acid and rigid bis(imidazol-1yl)benzene linkers. • The structural diversity depending on ligands and coordination number of metal centers. • Fluorescent sensor for the detection of acetone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2019

HERMeS Thruster Magnetic Field Topology Optimization Study: Performance, Stability, and Wear Results

NASA's Hall Effect Rocket with Magnetic Shielding (HERMeS) 12.5kW Technology Demonstration Unit-1 (TDU-1) has been the subject of extensive technology maturation in preparation for flight system development. The TDU-1 thruster implements a magnetically shielded field topology and has demonstrated the elimination of the discharge channel erosion as a life limiting mechanism. Extensive wear testing the TDU Hall thrusters has identified the thruster front pole covers as the next life limiting component. This effort aims to explore and investigate alternate magnetic field topologies to assess whether reductions in the front pole cover erosion can be attained while still maintaining very low erosion rates on the discharge channel walls. Four candidate magnetic field topologies that reduce the effectiveness of the shielding along the discharge channel walls with the intent to also reduce the erosion rates along the front pole covers were designed. Three of the four candidate magnetic field topologies (B1, B2, and B4) have been manufactured and were subjected to an extensive test campaign that includes laser induced fluorescence (LIF), performance, stability, wear, plume, and thermal characterization. In Phase I, LIF measurements along the discharge chamber centerline found that upstream retraction of the thruster's peak magnetic field does result in an upstream shift of the acceleration zone but the magnitude of the shift does not correspond one-to-one to the shift in the location of the peak radial magnetic field magnitude. Phase II test segment results found at a normalized thruster magnetic field setting of 1, the thruster performance was very similar for all configurations. Discharge current waveforms ripple data indicated that configurations Bo, B1, and B2 has similar oscillatory profiles with the B2 configuration transitioning to higher oscillatory mode at 400V instead of the 450V observed for configurations B0 and B1. Configuration B4 ripple data indicates that the thruster was operating in a very oscillatory mode above 350V. Inner front pole cover erosion rates for configuration B1 were approximately 65% relative to B0 and the erosion rates for configuration B2 were 40% relative to B0. Analysis of the performance, stability, and wear test results of configurations B0, B1, B2, and B4 indicates that configuration B2 may present an alternative option for B0 because it mostly maintained the same performance and stability as configuration B0 but with 40% of the inner front pole cover erosion rate.

Kamhawi, Hani↗
DOE OSTI2023

Desolvation Processes in Channel Solvates of Niclosamide

The antiparasitic drug niclosamide (NCL) is notable for its ability to crystallize in multiple 1:1 channel solvate forms, none of which are isostructural. Here, using a combination of time-resolved synchrotron powder X-ray diffraction and thermogravimetry, the process-induced desolvation mechanisms of methanol and acetonitrile solvates are investigated. Structural changes in both solvates follow a complicated molecular-level trajectory characterized by a sudden shift in lattice parameters several degrees below the temperature where the desolvated phase first appears. Model fitting of kinetic data obtained under isothermal heating conditions suggests that the desolvation is rate-limited by the nucleation of the solvent-free product. The desolvation pathways identified in these systems stand in contrast to previous investigations of the NCL channel hydrate, where water loss by diffusion initially yields an anhydrous isomorph that converts to the thermodynamic polymorph at significantly higher temperatures. Taking the view that each solvate lattice is a unique “pre-organized” precursor, a comparison of the pathways from different starting topologies to the same final product provides the opportunity to reevaluate assumptions of how various factors (e.g., solvent binding strength, density) influence solid-state desolvation processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Performance analysis of distributed symmetric sparse matrix vector multiplication algorithm for multi-core architectures

In this article, sparse matrix vector multiply (SpMVM) is an important kernel that frequently arises in high performance computing applications. Due to its low arithmetic intensity, several approaches have been proposed in literature to improve its scalability and efficiency in large scale computations. In this paper, our target systems are high end multi-core architectures and we use messaging passing interface + open multiprocessing hybrid programming model for parallelism. We analyze the performance of recently proposed implementation of the distributed symmetric SpMVM, originally developed for large sparse symmetric matrices arising in ab initio nuclear structure calculations. We also study important features of this implementation and compare with previously reported implementations that do not exploit underlying symmetry. Our SpMVM implementations leverage the hybrid paradigm to efficiently overlap expensive communications with computations. Our main comparison criterion is the "CPU core hours" metric, which is the main measure of resource usage on supercomputers. We analyze the effects of topology-aware mapping heuristic using simplified network load model. Furthermore, we have tested the different SpMVM implementations on two large clusters with 3D Torus and Dragonfly topology. Our results show that the distributed SpMVM implementation that exploits matrix symmetry and hides communication yields the best value for the "CPU core hours" metric and significantly reduces data movement overheads.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Reverse-mode differentiation in arbitrary tensor network format: with application to supervised learning.

This paper describes an efficient reverse-mode differentiation algorithm for contraction operations of tensor networks that may have arbitrary and unconventional network topologies. The approach leverages the tensor contraction tree of Evenbly and Pfeifer (2014), which provides an instruction set for the contraction sequence of a network. We show that this tree can be efficiently leveraged for differentiation of a full tensor network contraction using a recursive scheme that exploits (1) the bilinear property of contraction and (2) the property that trees have single path from root to leaves. While differentiation of tensor-tensor contraction is already possible in most automatic differentiation packages, we show that exploiting these two additional properties in the specific context of contraction sequences can improve efficiency. Following a description of the algorithm and computational complexity analysis, we investigate its utility for gradient-based supervised learning for low-rank function recovery and for fitting real-world unstructured datasets. We demonstrate improved performance over alternating least-squares optimization approaches and the capability to handle heterogeneous and arbitrary tensor network formats. When compared to alternating minimization algorithms, we find that the gradient-based approach requires a smaller oversampling ratio (number of samples compared to number model parameters) for recovery. This increased efficiency extends to fitting unstructured data of varying dimensionality and when employing a variety of tensor network formats. Here, we show improved learning using the hierarchical Tucker method over the tensor-train in high-dimensional settings on a number of benchmark problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Network analysis of transcriptomics expands regulatory landscapes in Synechococcus sp. PCC 7002

Cyanobacterial regulation of gene expression must contend with a genome organization that lacks apparent functional context, as the majority of cellular processes and metabolic pathways are encoded by genes found at disparate locations across the genome. In addition, the fact that coordinated regulation of cyanobacterial cellular machinery takes place with significantly fewer transcription factors, compared to other Eubacteria, suggests the involvement of post-transcriptional mechanisms and regulatory adaptations which are not fully understood. Global transcript abundance from model cyanobacterium Synechococcus sp. PCC 7002 grown under 42 different conditions was analyzed using context-likelihood of relatedness. The resulting 903-gene network, which was organized into 11 modules, not only allowed classification of cyanobacterial responses to specific environmental variables but provided insight into the transcriptional network topology and led to the expansion of predicted regulons. When used in conjunction with genome sequence, the global transcript abundance allowed identification of putative post-transcriptional changes in expression as well as novel potential targets of both DNA binding proteins and asRNA regulators. The results offer a new perspective into the multi-level regulation that governs cellular adaptations of fast-growing physiologically robust cyanobacterium Synechococcus sp. PCC 7002 to changing environmental variables. It also extends a methodological knowledge-based framework for studying multi-scale regulatory mechanisms that operate in cyanobacteria. Finally, it provides valuable context for integrating systems-level data to enhance evidence-driven genomic annotation, especially in organisms where traditional context analyses cannot be implemented due to lack of operon-based functional organization.

Environmental Molecular Sciences Laboratory↗
NASA NTRS1995

Managing Complexity in Multidisciplinary Visualization

As high performance computing technology progresses, computational simulations are becoming more advanced in their capabilities. In the computational aerosciences domain, single discipline steady-state simulations computed on a single grid are far from the state-of-the-art. In their place are complex, time-dependent multidisciplinary simulations that attempt to model a given geometry more realistically. The product of these multidisciplinary simulations is a massive amount of data stored in different formats, grid topologies, units of measure, etc., as a result of the differences in the simulated physical domains. In addition to the challenges posed by setting up and performing the simulation, additional challenges exist in analyzing computational results. Visualization plays an important role in the advancement of multidisciplinary simulations. To date, visualization has been used to aid in the interpretation of large amounts of simulation data. Because the human visual system is effective in digesting a large amount of information presented graphically, visualization has helped simulation scientists to understand complex simulation results. As these simulations become even more complex, integrating several different physical domains, visualization will be critical to digest the massive amount of information. Another important role for visualization is to provide a common communication medium from which the domain scientists can use to develop, debug, and analyze their work. Multidisciplinary analyses are the next step in simulation technology, not only in computational aerosciences, but in many other areas such as global climate modeling. Visualization researchers must understand and work towards the challenges posed by multidisciplinary simulation scenarios. This paper addresses some of these challenges, describing technologies that must be investigated to create a useful visualization analysis tool for domain scientists.

Miceli, Kristina D.↗
DOE OSTI2013

In silico analysis of the three-dimensional structures of the homodimer of uridine phosphorylase from Yersinia Pseudotuberculosis in the ligand-free state and in a complex with 5-fluorouracil

Pseudotuberculosis is an acute infectious disease characterized by a lesion of the gastrointestinal tract. A positive therapeutic effect can be achieved by selectively suppressing the activity of uridine phosphorylase from the causative agent of the disease Yersinia pseudotuberculosis. The synergistic effect of a combination of the chemotherapeutic agent 5-fluorouracil and antimicrobial drugs, which block the synthesis of pyrimidine bases, on the cells of pathogenic protozoa and bacteria is described in the literature. The three-dimensional structures of uridine phosphorylase from Yersinia pseudotuberculosis (YptUPh) both in the ligand-free state and in complexes with pharmacological agents are unknown, which hinders the search for and design of selective inhibitors of YptUPh. The three-dimensional structure of the ligand-free homodimer of YptUPh was determined by homology-based molecular modeling. The three-dimensional structure of the subunit of the YptUPh molecule belongs to {alpha}/{beta} proteins, and its topology is a three-layer {alpha}/{beta}/{alpha} sandwich. The subunit monomer of the YptUPh molecule consists of 38% helices and 24% {beta} strands. A model of the homodimer structure of YptUPh in a complex with 5-FU was obtained by the molecular docking. The position of 5-FU in the active site of the molecule is very consistent with the known data on the X-ray diffraction structures of other bacterial uridine phosphorylases (the complex of uridine phosphorylase from Salmonella typhimurium (StUPh) with 5-FU, ID PDB: 4E1V and the complex of uridine phosphorylase from Escherichia coli (EcUPh) with 5-FU and ribose 1-phosphate, ID PDB: 1RXC).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Distribution of Bound and Free Water in Anatomical Fractions of Pine Residues and Corn Stover as a Function of Biological Degradation

Biomass quality is influenced by water’s abundance, distribution, and status in relation to other chemical species within the polymer matrix. Water interacts with polymers that make up the cell walls, and these interactions govern the physical and chemical changes that occur during the storage and preprocessing of biomass feedstocks. Time-domain nuclear magnetic resonance (TD-NMR) was employed to explore variations in the physical constraints of water within the lignocellulosic microstructure in distinct anatomical fractions of biomass and as a function of biological degradation. The Carr–Purcell–Meiboom–Gill sequence, when combined with knowledge of the chemical composition and physical structure of pine residues and corn stover anatomical fractions, gives an accurate measurement of the bound and free water. In this work, the impacts of storage and biological degradation were investigated to elucidate changes in the status and distribution of water within distinct plant tissues. We also investigate how degradation during storage affects water interactions in different pine residues (e.g., bark, branch, and needle) and corn stover (e.g., cob, leaf, and stalk) anatomical fractions using transverse relaxation times (T2). As demonstrated herein, TD-NMR provides quantitative data on lignocellulosic biomass–water interactions within anatomical fractions, which can further aid in the investigation of preprocessing effects on feedstock quality. Our findings suggest that biological heating enhances biomass–water interactions at the cellular and macromolecular scale. In addition, analysis of three-dimensional scanning electron microscopy reconstructions indicates that surface roughness wavelengths align with microscale roughness, suggesting that pine forestry residue and corn stover particles have primarily hydrophobic exterior surfaces. Furthermore, this study offers multiscale insights into understanding the microstructure, wettability, and chemical environment that dictate diffusion, enzyme access, and recalcitrance of lignocellulosic biomass.

09 BIOMASS FUELS↗
DOE OSTI2026

Characterizing skyrmion flow phases with principal component analysis

Principal component analysis (PCA) is a powerful method that can identify patterns in large, complex data sets by constructing low-dimensional order parameters from higher-dimensional feature vectors. There are increasing efforts to use space-and-time-dependent PCA to detect transitions in nonequilibrium systems that are difficult to characterize with equilibrium methods. Here, we demonstrate that feature vectors incorporating the position and velocity information of driven skyrmions moving through random disorder permit PCA to resolve different types of disordered skyrmion motion as a function of driving force and the ratio of the Magnus force to the dissipation. Since the Magnus force creates gyroscopic motion and a finite Hall angle, skyrmions can exhibit a greater range of flow phases than what is observed in overdamped driven systems with quenched disorder. We show that in addition to identifying previously known skyrmion flow phases, PCA detects several additional phases, including different types of channel flow, moving fluids, and partially ordered states. Guided by the PCA analysis, we further characterize the disordered flow phases to elucidate the different microscopic dynamics and show that the changes in the PCA-derived order parameters can be connected to features in bulk transport measures, including the transverse and longitudinal velocity-force curves, differential conductivity, topological defect density, and changes in the skyrmion Hall angle as a function of drive. We discuss how asymmetric feature vectors can be used to improve the resolution of the PCA analysis, and how this technique can be extended to find disordered phases in other nonequilibrium systems with time-dependent dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2012

MAGNETIC FIELDS OF AN ACTIVE REGION FILAMENT FROM FULL STOKES ANALYSIS OF Si I 1082.7 nm AND He I 1083.0 nm

Vector magnetic fields of an active region filament in the photosphere and upper chromosphere are obtained from spectro-polarimetric observations recorded with the Tenerife Infrared Polarimeter (TIP II) at the German Vacuum Tower Telescope. We apply Milne-Eddington inversions on full Stokes vectors of the photospheric Si I 1082.7 nm and the upper chromospheric He I triplet at 1083.0 nm to obtain the magnetic field vector and velocity maps in two atmosphere layers. We find that (1) a complete filament was already present in H{alpha} at the beginning of the TIP II data acquisition. Only a partially formed one, composed of multiple small threads, was present in He I. (2) The AR filament comprises two sections. One shows strong magnetic field intensities, about 600-800 G in the upper chromosphere and 800-1000 G in the photosphere. The other exhibits only comparatively weak magnetic field strengths in both layers. (3) The Stokes V signal is indicative of a dip in the magnetic field strength close to the chromospheric PIL. (4) In the chromosphere, consistent upflows are found along the PIL flanked by downflows. (5) The transversal magnetic field is nearly parallel to the PIL in the photosphere and inclined by 20 Degree-Sign -30 Degree-Sign in the chromosphere. (6) The chromospheric magnetic field around the filament is found to be in normal configuration, while the photospheric field presents a concave magnetic topology. The observations are consistent with the emergence of a flux rope with a subsequent formation of a filament.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

What does fault tolerant Deep Learning need from MPI?

Deep Learning (DL) algorithms have become the {\em de facto} Machine Learning (ML) algorithm for large scale data analysis. DL algorithms are computationally expensive -- even distributed DL implementations which use MPI require days of training (model learning) time on commonly studied datasets. Long running DL applications become susceptible to faults -- requiring development of a fault tolerant system infrastructure, in addition to fault tolerant DL algorithms. This raises an important question: {\em What is needed from MPI for designing fault tolerant DL implementations?} In this paper, we address this problem for permanent faults. We motivate the need for a fault tolerant MPI specification by an in-depth consideration of recent innovations in DL algorithms and their properties, which drive the need for specific fault tolerance features. We present an in-depth discussion on the suitability of different parallelism types (model, data and hybrid); a need (or lack thereof) for check-pointing of any critical data structures; and most importantly, consideration for several fault tolerance proposals (user-level fault mitigation (ULFM), Reinit) in MPI and their applicability to fault tolerant DL implementations. We leverage a distributed memory implementation of Caffe, currently available under the Machine Learning Toolkit for Extreme Scale (MaTEx). We implement our approaches by extending MaTEx-Caffe for using ULFM-based implementation. Our evaluation using the ImageNet dataset and AlexNet neural network topology demonstrates the effectiveness of the proposed fault tolerant DL implementation using OpenMPI based ULFM.

Amatya, Vinay C.↗
DOE OSTI2012

PHOTOSPHERIC FLUX CANCELLATION AND THE BUILD-UP OF SIGMOIDAL FLUX ROPES ON THE SUN

In this study we explore the scenario of photospheric flux cancellation being the primary formation mechanism of sigmoidal flux ropes in decaying active regions. We analyze magnetogram and X-ray observations together with data-driven non-linear force-free field (NLFFF) models of observed sigmoidal regions to test this idea. We measure the total and canceled fluxes in the regions from MDI magnetograms, as well as the axial and poloidal flux content of the modeled NLFFF flux ropes for three sigmoids-2007 February, 2007 December, and 2010 February. We infer that the sum of the poloidal and axial flux in the flux ropes for most models amounts to about 60%-70% of the canceled flux and 30%-50% of the total flux in the regions. The flux measurements and the analysis of the magnetic field structure show that the sigmoids first develop a strong axial field manifested as a sheared arcade and then, as flux cancellation proceeds, form long S-shaped field lines that contribute to the poloidal flux. In addition, the dips in the S-shaped field lines are located at the sites of flux cancellation that have been identified from the MDI magnetograms. We find that the line-of-sight-integrated free energy is also concentrated at these locations for all three regions, which can be liberated in the process of eruption. Flare-associated brightenings and flare loops coincide with the location of the X-line topology that develops at the site of most vigorous flux cancellation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

Pioneer 10 and 11 Data Analysis

This report finishes the work of NASA Grant NAS2-153, which supported data analysis for the UCSD instruments on Pioneers 10 and 11. The data analyzed under this grant span 22 years of interplanetary measurements in the inner and outer heliosphere. The UCSD instruments made their mark in cosmic ray research based upon their high energy thresholds, directional responses, and reliable data streams. one of their primary scientific objectives concerns the size, configuration, and time behavior of the heliosphere. The size scale is inferred from the radial intensity gradient, which is measured between the two spacecraft and extrapolated to interstellar intensity levels at the cosmic ray modulation boundary. This boundary still eludes us, and its position, motion, and the best method of extrapolation are ongoing problems. Current projections place the boundary beyond 100 AU, which may be beyond the termination shock, and raises the question of possible modulation in the heliosheath. Probably our only hope of seeing this region in the immediate future rides on the possibility that the boundary will move inward. Our instruments have recorded many Forbush, or transient, decreases in the outer heliosphere. These observations led us to a model that attributes many of the decreases to solar wind stream-stream interactions, and relates the cosmic ray variations to the locally observed magnetic field magnitude. As the cosmic ray variations in this model result only from topological changes in the modulation integral, the model is a tool for studying the possibility that the 11 year cosmic ray modulation cycle can be accounted for by a superposition of Forbush decreases. The cosmic ray angular distribution function is measurable, given a good telemetry rate, by the UCSD Cerenkov detector which counts particles of energy greater than 500 MeV/n. We obtained statistically significant samples from 1 to 9 AU, at 13 AU, and at 34 AU. The anisotropy tends to be a few tenths of a per cent at all radial distances. A quasiperiodic variation in the east-west anisotropy with period of about 50 days remains unexplained.

Fillius, Walker↗
DOE OSTI2016

Grizzly: A Multi-scale and Multi-Physics Tool to Model Aging of Nuclear Power Plant Components

Grizzly is a multi-scale multi-physics tool being developed at Idaho National Laboratory (INL) as part of the US Department of Energy's Light Water Reactor Sustainability program to provide improved safety assessments of systems, components, and structures in nuclear power plants subjected to age-related degradation. Its goal is to provide an improved scientific basis for decisions surrounding license renewal, which would permit operation of commercial nuclear power plants beyond 60 years. Grizzly is based on INL's MOOSE framework, which enables multi-physics simulations in a massively parallel computing environment. The reactor pressure vessel (RPV) has been chosen as the initial application for Grizzly. The primary concern for the RPV is that it becomes increasingly susceptible to fracture initiation at the sites of pre-existing flaws over time. This is because exposure to irradiation and high temperatures embrittles the steel. Assessing the safety of an RPV under transient conditions involves simulations at multiple scales. At the engineering scale, strongly coupled equations of heat conduction and solid mechanics are solved to simulate the global response of the RPV to accident conditions. Sub-models are used to represent regions with pre-existing flaws and the stress intensity factors are evaluated from the domain integrals. Crack propagation is then determined using a model that provides the temperature and neutron fluence dependent fracture toughness of the RPV steel. The ability of Grizzly to provide accurate engineering scale estimates of RPV safety beyond 60 years of operation strongly relies on the predictability of the embrittlement model. Multiple models for embrittlement have been proposed, and the physically based empirical model is implemented in Grizzly. However, the formation of late blooming phases greatly reduces the confidence on the extrapolated values obtained from embrittlement model, which is further exacerbated by the absence of real field data for validation. Experimentally validated lower length scale models can significantly improve the predictability of the embrittlement model beyond the calibrated range and the development of such models is another ongoing effort in Grizzly. The modeling activity in Grizzly spans across different length scales, and can be divided into: (i) Irradiation induced microstructure evolution at the atomistic and single crystal length scale using molecular dynamics, atomistic kinetic Monte-Carlo, phase-field and rate theory; (ii) Effect of the irradiation induced defects on the flow stress and fracture processes at the polycrystalline scale using crystal plasticity. Subsequently, the prediction from these models can be used to obtain irradiation induced ductile-to-brittle transition temperature (DBTT) shifts from specimen-scale simulations; (iii) Global thermo-mechanical response of the RPV under normal and accident scenarios, and probabilistic risk analyses at the engineering scale using the DBTT shifts obtained from the lower length scale models. The item (i) has been discussed in a previous report. In this extended abstract an overview of items (ii) and (iii) is provided. The extended abstract summarizes the current status of the modeling effort under Grizzly. At the engineering scale, reduced order model development is underway using the statistical tool RAVEN. This will allow probabilistic analysis considering a wide range of flaw topologies in a computationally efficient manner. Coupling the irradiation induced defect growth predictions to crystal plasticity is also in progress. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1985

Deformation of the heliospheric current sheet

A kinematic analysis of the previously neglected effect of velocity inhomogeneity on the topology of the heliospheric current sheet in a radially flowing solar wind shows how the originally smooth current sheet becomes 'ruffled'. In the highly idealized case of a totally uniform, radial, steady solar wind, the shape of the current sheet is independent of distance from the sun. However, the real solar wind is inhomogeneous; the velocity varies from point to point along the current sheet, causing a distortion in the current sheet of progressively greater amplitude with increasing distance from the sun. This is true even for purely radial flow. Significant and observable distortion is produced by relatively small gradients in velocity; thus to predict or understand the shape of the heliospheric current sheet it is essential to know the solar wind in which the current sheet is embedded. Examples are given of mild velocity gradients which demonstrate the principles, the magnitude, and the character of the effect; deformation of the actual heliospheric current sheet in the highly variable solar wind is expected to be of far greater amplitude and complexity than in the simplified, tutorial examples. A new expression for the inclination of the current sheet as a function of velocity inhomogeneity and distance from the sun that is easily applied to the interpretation of solar wind data is also derived.

Suess, S. T.↗