Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “slip transfer”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Boosting Photoelectric Conductivity in Porphyrin-Based MOFs Incorporating C 60

Electronic structure calculations show that guest C60 in the porphyrin-containing metal organic frameworks Zn 2 (TCPB)(DA-ZnP) (DA-MOF; H 4 TCPB = 1,2,4,5-tetrakis(4-carboxyphenyl)benzene, DA-ZnP = [5,15-bis[(4-pyridyl)ethynyl]-10,20-diphenylporphinato]zinc(II)) and Zn 2 (TCPB)(F-ZnP) (F-MOF; F-ZnP = [5,15-di(4-pyridyl)-10,20-bis(pentafluorophenyl)porphinato]zinc(II)) engenders high photoelectrical conductivity due to efficient donor–acceptor charge-transfer (CT) interactions. Structural modifications at the meso position of the porphyrin influence the preferred positions of C 60 within the frameworks, giving rise to host–guest interactions with different anisotropic structural, electronic, and optoelectronic properties. A preferred slipped-parallel π-stacked interaction of C 60 that is predicted for NH 2 -substituted DA-MOF and F-MOF fosters strong CT transitions and lowers band gaps by ~1.0 eV compared to the pristine DA-MOF and F-MOF. Hopping rates computed using Marcus theory are found to be anisotropic and accelerated by multiple orders of magnitude across π-stacked interfaces created by C 60 incorporation, a consequence of strong electronic coupling between initial and final diabatic states. Calculations indicate that photoinduced electron transfer, as well as direct CT from porphyrin to C 60 upon irradiation, triggers a charge separation process that leads to the formation of what should be long-lived electron-trapped states at the heterojunctions. Lastly, design principles revealed here for the control of photophysical and electron-transfer processes will be useful for constructing new MOF-derived visible- and infrared-based optoelectronics.

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

The feasibility of electron beam slip stacking in EIC Rapid Cycling Synchrotron

Announced in 2020, the Electron Ion Collider will be constructed at Brookhaven National Lab (BNL) by modifying the existing Relativistic Heavy Ion Collider (RHIC) facility. Among the alterations include the installation of an electron storage ring (ESR) and a Rapid Cycling Synchrotron (RCS) in the existing RHIC tunnel. The RCS will accelerate polarized electrons from 400MeV to 5GeV,10GeV and 18 GeV and transfer it to the ESR. The LINAC source produces 7nC per shot meanwhile the RCS requires 28nC per bunch of charge. Currently, it is planned to inject 4 bunches into adjacent RF buckets and then merge them into a single RF bucket using fields from several RF cavities. However, this requires a kicker with an exceptionally fast rise and fall time to inject into adjacent RF bucket without disturbing the existing circulating bunch and this requires harmonic RF kickers. Our study focuses on determining the feasibility of an alternate particle injection method that will utilize slip stacking, a method that accommodates the circulation of multiple beams with different momentum in the same azimuthal space of the beampipe, thus permitting the accumulation of charge from several shots into a single RF bucket. We find the momentum aperture required for slip stacking using current RCS beam parameters exceeds the RCS’s momentum deviation limit,dp/p = 1.5%. We determine the appropriate dimensions the beam must have to enable slip stacking and conduct a literature review of slip stacking schemes in other facilities. The literature review tentatively indicates that slip stacking has not been attempted on electron beams at a major accelerator facility.. The slip stacking approach has the potential to save the cost of the harmonic RF kicker and bunch merging cavities.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Design Considerations and Operational Features of the Collimators for the Fermilab Main Injector and Recycler

The Fermilab Main Injector system delivers 700 kW of 120 GeV Proton beam for neutrino experiments. Since 2013 this has been achieved using slip stacking accumulation in the Recycler with up to 12 batches from the Fermilab Booster per Main Injector Ramp Cycle. To control activation from beam loss, collimation systems in the Booster to Recycler transfer line, in the Recycler and in the Main Injector are employed. Residual radiation measurements around the ring with detailed studies at the collimators are required to maintain adequate loss control. We will review design considerations, operational parameters and activation results for more than ten years of operation. Simulations with MARS15 are used to explore the activation rates and the isotopic composition of the resulting activation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Flexoelectricity and New Phenomena

It is easy to miss the scientific implications of our recent work on Triboelectricity. Everyone knows that rubbing and contact can produce static electricity; less appreciated is that the thermodynamic driver has been an open question since static electricity was first observed by Thales of Miletus around 585 BC. People think they understand it, for instance one common explanation that can still be found in the current literature is that differences in the work function drives charge transfer, often called the Volta-Helmholtz hypothesis. As summarized in 1967 by Harper, this fails to explain many experimental observations, for instance that charging can occur when two pieces of the same material are rubbed against each other. We are the first to place triboelectricity on a solid foundation rooted in quantum mechanics – the flexoelectric effect. We were able to explain a significant number of previously unexplained phenomena: 1) Bipolar tribocurrents associated with stick-slip, due to the change in sign of the strain gradients. 2) A one-third power scaling of tribocurrents with indentation force. 3) Tribocharging when two identical materials are used, these being due to local variations in the asperities so there are usually local potential differences. 4) Inhomogeneous charging of insulators, related to the statistical nature of asperities. 5) An experimentally observed reversal in the sign of charge transfer for negative and positive curvature, which is related to a change in the sign of the strain gradient. Exploiting our DOE prior funded work on flexoelectricity, we obtained semi-quantitative matching to existing experimental measurements of the surface charge in triboelectric experiments. The work has been well received in the literature. The work has been the focus of a number of popular science press articles, and also formed the basis for a Podcast for children 6-10 “The Rise and Fall of Static Man” posted in December 2019 by NPR as part of their “Wow in the World” series. I was also briefly interviewed by the Chicago PBS station in January 2020. This work are significant for a wide range of energy applications; to quote from an independent source: "Triboelectric power has plenty of potential, says Wenzhuo Wu, an assistant professor of engineering at Purdue. If the basics of static electricity are better understood, we could maximize the efficiency of wind or wave power generators, Wu says. The body's own movement could be used to power internal medical devices. Imagine being able to create a roof shaped to harness the power of a raindrop — the friction of the rain passing over the surface — to generate triboelectricity, powering the building below it." This is the start of new science, some of which we already partially understand such as the role of band bending in charge transfer. We need to understand charge transfer combining elasticity, quantum mechanics, band bending and defect states. These directly involve several of the DOE Grand Challenges: "How do we control material processes at the level of electrons? How do remarkable properties of matter emerge from complex correlations of the atomic or electronic constituents and how can we control these properties? How do we characterize and control matter away—especially very far away— from equilibrium?" I will argue that this work truly falls into the class of disruptive science; it is not just a simple extension, linear science. Not everyone will accept the approach. Since we explain far more about triboelectricity than anyone before, the preponderance of evidence supports the model. The feedback I have received is that many agree with the work, to quote: "The model makes sense, says Michael McAlpine, a professor of engineering at the University of Minnesota. "It's such a simple explanation, I was surprised I didn't put my finger on that," McAlpine says." The proposal received strong reviews. It was also publicized on the Department of Energy Web Site.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Strengthening mechanisms, hardening/softening behavior, and microstructure evolution in an LPSO magnesium alloy at elevated temperatures

Here, rapid degradation of the mechanical properties of conventional Mg alloys with temperature precludes their wide application in industry. The LPSO (long period stacking ordered) Mg alloys, on the other hand, show considerable potential to achieve excellent high-temperature mechanical performance. In this study, a Mg-5.1Y-2.4Zn (wt.%) extrusion alloy with the microstructure comprising of α-Mg matrix and LPSO phase was prepared. It is shown that the yield strength of the alloy at 200 °C retained ~83% of the yield strength of 198 MPa at room temperature (RT). Importantly, the ductility also improved, owing to the increased non-basal slip activities at elevated temperatures. In contrast to RT and 200 °C, the alloy exhibited softening tensile behavior with an apparent loss of yield strength at 300 °C. In-situ lattice strain analysis revealed that load transfer from α-Mg to LPSO remains effective at temperatures up to 200 °C, governing the high yield strength of the alloy. The post-mortem analysis by means of electron microscopy confirmed that the LPSO phase is structurally stable at 200 °C. In particular, the high strain hardenability is rationalized by the sandwiched LPSO structures together with stacking faults serving as effective obstacles to the motion of non-basal dislocations. The load transfer effect disappears at the temperature of 300 °C, resulting in the fast deterioration of strength of the alloy. The cracking of the LPSO phase along with the recrystallizing of α-Mg grains were found to be the main causes for the softening behavior of the alloy at 300 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗
DOE OSTI2020

Energy storage in magnetic textures driven by vorticity flow

Here, an experimentally feasible energy-storage concept is formulated based on vorticity (hydro)dynamics within an easy-plane insulating magnet. The free energy associated with the magnetic winding texture is built up in a circular easy-plane magnetic structure by injecting a vorticity flow in the radial direction. The latter is accomplished by electrically induced spin-transfer torque, which pumps energy into the magnetic system in proportion to the vortex flux. The resultant magnetic metastable state with a finite winding number can be maintained indefinitely because the process of its relaxation via phase slips is exponentially suppressed when the temperature is brought well below the Curie temperature. We characterize the vorticity-current interaction underlying the energy-loading mechanism through its contribution to the effective electric inductance in the rf response. Our proposal may open an avenue for naturally powering spintronic circuits and nontraditional magnet-based neuromorphic networks.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Data-Driven Modeling and Correction of Vehicle Dynamics

We develop a data-driven framework for learning and correcting nonautonomous vehicle dynamics. Physics-based vehicle models are often simplified for tractability and therefore exhibit inherent model-form uncertainty, motivating the need for data-driven correction. Moreover, nonautonomous dynamics are governed by time-dependent control inputs, which pose challenges in learning predictive models directly from temporal snapshot data. To address these, we reformulate the vehicle dynamics via a local parameterization of the time-dependent inputs, yielding a modified system composed ofa sequence of local parametric dynamical systems. Here, we approximate these parametric systems using two complementary approaches. First, we employ the dimension reduction and interpolation in parameter space (DRIPS) methodology to construct efficient linear surrogate models, equipped with lifted observable spaces and manifold-based operator interpolation. This enables data-efficient learning of vehicle models whose dynamics admit accurate linear representations in the lifted spaces. Second, for more strongly nonlinear systems, we employ flow map learning (FML), a deep neural network (DNN) approach that approximates the parametric evolution map without requiring special treatment of nonlinearities. We further extend FML with a transfer-learning-based model correction procedure, enabling the correction of misspecified prior models using only a sparse set of high-fidelity or experimental measurements, without assuming a prescribed form for the correction term. Through a suite of numerical experiments on unicycle, simplified bicycle, and slip-based bicycle models, we demonstrate that DRIPS offers robust and highly data-efficient learning of nonautonomous vehicle dynamics, while FML provides expressive nonlinear modeling and effective correction of model-form errors under severe data scarcity.

data-driven modeling↗
DOE OSTI2025

Solvent Polarity Independent Symmetry-Breaking Charge Separation in a Slip-Stacked Covalent Terrylene Monoimide Dimer

The design of organic materials capable of efficient photoinduced charge separation in low-polarity environments is a critical challenge for advancing organic photovoltaics. Symmetry-breaking charge separation (SB-CS) offers a promising route to efficient charge separation with minimal energy loss; however, in conventional organic materials this process typically relies on polar solvents to stabilize the charge-separated state. Here, we investigate a slip-stacked terrylene monoimide dimer (TMI 2 ) engineered with intrinsic electronic asymmetry to circumvent this limitation. We demonstrate through a comprehensive suite of ultrafast spectroscopic techniques that TMI 2 undergoes remarkably efficient SB-CS regardless of solvent polarity enabled by the permanent dipole moment of TMI and the accessible intramolecular charge-transfer (ICT) character of its monomer units. Furthermore, analysis of vibronic coherences reveals that the initial state mixing is actively driven by a 193 cm –1 intermolecular mode, while the final SB-CS state is marked by a distinct vibrational fingerprint of the radical ion pair product. Our findings reveal how the intrinsic electronic asymmetry of the TMI monomers that constitute TMI 2 creates an efficient charge separation pathway that precludes the need for external stabilization from a polar solvent. Here, this work establishes a powerful molecular design principle for leveraging intrinsic monomeric asymmetry to achieve efficient charge separation in low-dielectric environments, with significant implications for the development of next-generation organic optoelectronic materials.

Aromatic compounds↗
DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗
DOE OSTI2022

Late Paleozoic Evolution of the Anadarko Basin: Implications for Laurentian Tectonics and the Assembly of Pangea

Abstract Interplay between marginal and intraplate late Paleozoic tectonics in southeastern Laurentia is a subject of debate. In southern Laurentia, intraplate deformation has been attributed to combinations of compressive forces from orogenic belts in the west, southwest, and southeast. While left‐lateral transpression is generally an accepted mechanism for formation of uplifts in the Southern Oklahoma Aulacogen (SOA), to date no studies have explained the driving mechanisms of these kinematics. Here, we analyze the Southern Oklahoma Transpressional System (SOTS), a ∼500 km long ∼50 km wide fault zone that extends from the southeastern Laurentian margin into the plate interior along Cambrian lithospheric weaknesses inherited from the SOA. We use a compilation of published studies to present detailed, time‐integrated sedimentary thickness maps of SOTS‐related basins. We combine these maps with a new interpretation of three 2D seismic reflection lines to constrain timing and kinematics of all major late Paleozoic SOTS fault zones. We show that the SOTS kinematics evolved from contraction to left‐lateral strike‐slip during Mississippian‐Pennsylvanian time. Tectonic activity and fault kinematics in the SOTS are explained by a model in which a remnant ocean basin closed diachronously from northeast to southwest along the Southern Appalachian and Ouachita‐Marathon orogens. Left‐lateral transpression during diachronous closure of the Rheic oceanic basin along the Ouachita‐Marathon margin is related to rotation of the stress field due to changes in slab‐pull and eastward drift of Laurentia toward Gondwana. We speculate that the SOTS acted in part as a STEP (subduction‐transfer edge propagator) fault that followed the pre‐existing SOA.

Hobbs, Noah. F.↗
DOE OSTI2024

Resolving crystallographic geometrically necessary dislocations in three dimensions in a hexagonal close packed titanium alloy

Abstract Geometrically necessary dislocation (GND) content is measured from mm 3 -scaled Ti7Al three-dimensional (3D) microstructural data using a theory extended for hexagonal close packed crystals, which accounts for basal, prismatic and pyramidal ⟨ c + a ⟩ type dislocation content. The Ti7Al samples have been mechanically pre-strained to two different strain levels, and will then be strained along the same axis in uniaxial tension during simulation. Both inter- and intragranular GNDs across the microstructures have been characterized, with a large contribution of pyramidal ⟨ c + a ⟩ GNDs, consistent with the relative slip activity involved in pre-straining. The spatially resolved crystallographic GND distributions within the 3D microstructures are used to instantiate a microstructure model for forward modeling deformation simulations by a dislocation density hardening elasto-viscoplastic fast Fourier transform framework. Coarsening the voxel resolution during the initial microstructure construction procedure is shown to strongly impact both the magnitude and spatial distribution of the GNDs and in turn the forward deformation response of the pre-strained material. This study indicates that the voxel resolution desired when transferring from measured to model microstructures need not only be proportionally scaled with the microstructure but also sufficiently fine to capture the subgranular orientation gradients that may already be present in the material.

Materials Science↗
DOE OSTI2022

Implement and Test 3D Mortar Contact in BISON

We leverage the extension of the generation of mortar segment meshes to three dimensions in MOOSE’s framework to extend thermomechanical modeling capabilities to problems with three dimensions. A modular approach to gap heat transfer physics using the mortar finite element method was created and documented, mechanical contact was extended to three dimensions—including frictional behavior, performance and ease of use were improved, and steps towards scalability of solid mechanics problems involving contact were taken. Many of these new developments are demonstrated in the simulation of 3D light-water reactor (LWR) problems, where the thermomechanical interface problem is solved using the mortar finite element method. Usage of the mortar framework has improved convergence in 2D problems and has enabled employing friction in 3D problems, of which we show results of a short, local stack of 3D pellets. Consequently, the benefits of mortar in terms of solution convergence and quality are extended to three dimensions. Section 2 discusses fundamental developments that enabled the simulation of practical mortar problems in three dimensions and other general improvements, including the reduction of the derivative container size, the modification of dual basis computations when edge dropping (lack of secondary element projection) takes place, the improvement of conditioning when employing the VCP in-edge dropping conditions, and code usability and quality improvements. These latter code enhancements include the migration of tests using “old” mortar contact constraints to using dual mortar with a semi-smooth Newton solution strategy and the reuse of lower dimensional domains for straightforwardly setting up a mortar thermomechanical LWR problem, i.e. the MOOSE action is employed for mechanical contact and the thermal LWR action is employed to capture the gas conductance, contact, and radiation components of gap heat transfer physics. Independently of the mortar LWR thermal action, we developed a modular approach to gap heat transfer that resides in MOOSE and can be leveraged, e.g., in metallic fuel problems. This approach, whose code design based on MOOSE’s user objects to model specific physics was proposed by the maintenance activity, is detailed in Section 3. Based on the dual mortar finite element method, the frictional contact constraints were extended to three dimensions. A block sheared in two directions in and out of contact with a rigid plane is employed in Section 4 to show the way the approach handles changes in frictional states (e.g. stick to slip) within a competitive number of Newton iterations. Equations and numerical results on the use of the VCP with Cartesian Lagrange multipliers, whose combination enables their direct condensation, are described in Section 5.3. Two-dimensional and three-dimensional BISON LWR simulations are discussed in Section 6. Particularly, a stack of five eccentric pellets with a surface defect is simulated and the effect of pellet-cladding friction is assessed. Finally, conclusions are outlined in Section 7.

42 ENGINEERING↗
DOE OSTI2023

Geomechanical in situ testing of fault reactivation in argillite repositories

Abstract. Pressurization of natural faults as a result of repository-induced effects can lead to their reactivation and permeability generation in case such features are present near disposal tunnels. Potential driving forces for such pressurization are the temperature increase caused by heat-producing high-level radioactive waste and the generation of hydrogen and other gases due to corrosion of engineered materials. We are primarily concerned about pressurization and fault activation in host rocks and faults that have very low natural permeability for pore pressure increases to dissipate, such as the argillite rocks currently investigated in Switzerland, France and other countries. This presentation discusses a series of in situ experiments of fault activation by fluid injection conducted in the argillite rock (Opalinus clay) at the Mont Terri underground research laboratory in Switzerland. A multi-model monitoring test bed was installed at Mont Terri that includes distributed fiber optics for strain, temperature, and acoustics; local fault pore pressure and three-dimensional displacement sensors; active seismic imaging; and passive seismic monitoring. The fault experiments (and their subsequent analysis via hydromechanical modeling) provide a new fundamental understanding of the coupling between pore pressure, fault deformation and permeability generation as a function of time and allow exploring how seismic and aseismic events may impact the integrity of faulted argillite host rock. Our experimental observations demonstrate that significant flow (and transport) can occur along the initially impermeable argillite fault when rupture is activated. However, the fault permeability decreases to almost its pre-activation value when fluid injection ceases and fluid pressure drops. Dilatant slip on the fault plane alone does not explain the observed increase in fault permeability; pressure-induced fault opening also plays a role, favored by the softness of the shale along with the fact that the structure of the fault zone prevents fluids from diffusing into the adjacent damage zone. Rupture initiation and permeability generation is initially aseismic, associated only with an increase in noise level and emerging tremors. Micro-earthquakes are initiated later in the experiments and typically occur away from the fluid-pressurized area. Hydromechanical models show that stress transferred from the initial aseismic deformation can build up to stress criticality and later induce seismic rupture. After presenting the experimental results, we will close the presentation with an outlook to future experimental campaigns using the fault test bed at Mont Terri. We are currently planning a controlled thermal stimulation of the fault, where instead of fluid injection as a trigger mechanism we will heat up the nearby rock volume and measure potential effects on fault stability. Lessons learned from our past and future experiments will help inform the safety assessment of geologic disposal in argillite host rock.

Birkholzer, Jens T.↗
DOE OSTI2021

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The effect of microstructure on the relationship between grain boundary sliding and slip transmission in high purity aluminum

The relationship between grain boundary sliding (GBS) and slip transmission is investigated experimentally at grain boundaries in 99.99% aluminum with a through-thickness, coarse grained microstructure deformed in tension at 190 °C. Using scanning electron microscope enabled digital image correlation (SEM-DIC) and electron backscatter diffraction (EBSD), high resolution strain fields and microstructural information were measured to examine the influence of microstructural neighborhoods on interactions between GBS and slip transmission and strain localization. Several distinct cases are presented that highlight important microstructural factors that govern deformation near grain boundaries. The findings include (1) direct transmission and GBS were anti-compatible and facilitated by opposing boundary types (low misorientation and high energy grain boundaries respectively); (2) increased GBS activity was correlated with decreased indirect transmission behavior; (3) GBS accommodation at triple junctions was enabled by intragranular plasticity; and (4) the local intragranular plastic strain discontinuity between grains determined the magnitude of GBS gradients. Finally, this work provides insight into the nature of these mechanisms and can be used to identify strain transfer criteria that can lead to improved GBS-sensitive crystal plasticity models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Deformation Mechanisms at Grain Boundaries in Polycrystals (Final Report)

Despite its significance in polycrystalline materials, there have been few experimental investigations of the activity of grain boundary sliding (GBS) and the relationship between GBS and slip transmission at grain boundaries. This work addressed this knowledge gap by an experimental and analytical investigation into the interactions between grain boundary sliding (GBS) and slip transmission plastic deformation mechanisms, and how these are impacted by local microstructure. Polycrystalline columnar aluminum of 99.99% purity was examined as a highly modellable system for characterizing grain boundary-local deformation mechanism interaction. High-resolution, large field-of-view strain fields and grain orientation information were captured and statistically analyzed to identify deformation activity with respect to the local microstructure, using a combination of scanning electron microscopy and digital image correlation (SEM-DIC) with electron backscatter diffraction (EBSD). Alignment between the lower-resolution electron backscatter diffraction (EBSD) data and the higher-resolution strain fields was significantly improved by clustering of strain data within an EBSD-defined boundary mantle. Although this work specifically focuses on characterizing the relationship between grain boundary sliding (GBS) and slip transmission, the developed experimental approach, multimodal data alignment process, and automated displacement discontinuity measurement process can be used to study a variety of mechanism interactions across a broad range of materials. Regarding GBS and slip transmission interaction, significant findings include that (1) grain boundaries that slide easily do not facilitate slip transmission; (2) substantial gradients in GBS magnitude could be found along the length of individual grain boundaries, and these were caused by strain discontinuities across the boundary related to deformation anisotropy between adjacent grains; (3) direct transmission and GBS were anti-compatible and facilitated by opposing boundary types (low misorientation and high energy grain boundaries respectively); (4) increased GBS activity was correlated with decreased indirect transmission behavior; (5) GBS accommodation at triple junctions was enabled by intragranular plasticity; and (6) the local intragranular plastic strain discontinuity between grains determined the magnitude of GBS gradients. The experimental strain fields were compared to that of a tangent crystal plasticity finite element (CPFE) model developed by collaborators (Mr. Ajey Venkataraman and Professor Michael Sangid, Purdue University). Implication of this work for improving CPFE models and deformation prediction are discussed. This work provides insight into the nature of these mechanisms and can be used to identify strain transfer criteria that can lead to improved GBS-sensitive crystal plasticity models.

36 MATERIALS SCIENCE↗