Search NASA⌕ Search

SEARCH · Search NASA

Results for “matrix diffusion”

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

Mesoscale shock structure in particulate composites

Multiscale experiments in heterogeneous materials and the knowledge of their physics under shock compression are limited. Here, this study examines the multiscale shock response of particulate composites comprised of soda-lime glass particles in a PMMA matrix using full-field high speed digital image correlation (DIC) for the first time. Normal plate impact experiments, and complementary numerical simulations, are conducted at stresses ranging from 1.1 - 3.1 GPa to elucidate the mesoscale mechanisms responsible for the distinct shock structure observed in particulate composites. The particle velocity from the macroscopic measurement at continuum scale shows a relatively smooth velocity profile, with shock thickness decreasing with an increase in shock stress, and the composite exhibits strain rate scaling as the second power of the shock stress. In contrast, the mesoscopic response was highly heterogeneous, which led to a rough shock front and the formation of a train of weak shocks traveling at different velocities. Additionally, the normal shock was seen to diffuse the momentum in the transverse direction, affecting the shock rise and the rounding-off observed at the continuum scale measurements. The numerical simulations indicate that the reflections at the interfaces, wave scattering, and interference of these reflected waves are the primary mechanisms for the observed rough shock fronts.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Microstructure evolution and mechanical performance of nickel based superalloy C1023 at elevated temperatures

Highlights: • Growth and strengthening mechanisms of carbide in superalloy C1023 were investigated. • Carbide dendrite primary arm grows selectively based on the vertex of blocky carbonitride. • Monolayer or multilayer carbide dendrite arm has certain crystal relationship with γ′ and γ. • Plastic deformation of γ hinders interior crack propagating from MC to matrix at 650 °C. • Oxidation of GBs significantly weakens the tensile strength of the superalloy at 850 °C. - Abstract: Microstructure evolution of polycrystalline nickel based superalloy C1023 casting was investigated by comparing as cast superalloy and heat treated superalloy in the paper. Influences of strengtheners and carbides on the tensile behaviors of heat treated superalloy at different temperatures (20 °C, 650 °C, and 850 °C) were studied. Resulting data from deep corrosion suggested that growth mode of scripted carbide included nucleation on the vertex of blocky carbonitride, expansion of dendrite arm in C-rich area, and enlargement of dendrite cap. Decomposition of MC carbide generated numerous octahedral M{sub 23}C{sub 6} carbide as well as nanoscale Cr-Mo-rich nuclei. Diffusion of Mo from γ to precipitates such as γ', MC, and M{sub 23}C{sub 6}, occurred during the heat treatment while irregular γ' (4 nm–500 nm) transformed into cuboidal particles. At 650 °C, no significant deterioration of tensile strength was found while plastic deformation ability of γ base was improved. Oxidation of grain boundaries (GBs) dominated the crack behavior of the superalloy at 850 °C.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

Rationalization of Reduced Penetration of Drugs through Ceramide Gel Phase Membrane

Since computing resources have advanced enough to allow routine molecular simulation studies of drug molecules interacting with biologically relevant membranes, a considerable amount of work has been carried out with fluid phospholipid systems. However, there is very little work in the literature on drug interactions with gel phase lipids. This poses a significant limitation for understanding permeation through the stratum corneum where the primary pathway is expected to be through a highly ordered lipid matrix. To address this point, we analyzed the interactions of p-aminobenzoic acid (PABA) and its ethyl (benzocaine) and butyl (butamben) esters with two membrane bilayers, which differ in their fluidity at ambient conditions. We considered a dioleoylphosphatidylcholine (DOPC) bilayer in a fluid state and a ceramide 2 (CER2, ceramide NS) bilayer in a gel phase. We carried out unbiased (100 ns long) and biased z-constraint molecular dynamics simulations and calculated the free energy profiles of all molecules along the bilayer normal. The free energy profiles converged significantly slower for the gel phase. While the compounds have comparable affinities for both membranes, they exhibit penetration barriers almost 3 times higher in the gel phase CER2 bilayer. This elevated barrier and slower diffusion in the CER2 bilayer, which are caused by the high ordering of CER2 lipid chains, explain the low permeability of the gel phase membranes. We also compared the free energy profiles from MD simulations with those obtained from COSMOmic. This method provided the same trends in behavior for the guest molecules in both bilayers; however, the penetration barriers calculated by COSMOmic did not differ between membranes. Finally, we show how membrane fluid properties affect the interaction of drug-like molecules with membranes.

chemistry↗
DOE OSTI2019

Random walks on jammed networks: Spectral properties

Using random walk analyses we explore diffusive transport on networks obtained from contacts between isotropically compressed, monodisperse, frictionless sphere packings generated over a range of pressures in the vicinity of the jamming transition p → 0 . For conductive particles in an insulating medium, conduction is determined by the particle contact network with nodes representing particle centers and edges contacts between particles. The transition rate is not homogeneous, but is distributed inhomogeneously due to the randomness of packing and concomitant disorder of the contact network, e.g., the distribution of the coordination number. A narrow escape time scale is used to write a Markov process for random walks on the particle contact network. This stochastic process is analyzed in terms of spectral density of the random, sparse, Euclidean and real, symmetric, positive, semidefinite transition rate matrix. Findings reveal network structures derived from jammed particles have properties similar to ordered, euclidean lattices but also some unique properties that distinguish them from other structures that are in some sense more homogeneous. Specifically, the distribution of eigenvalues of the transition rate matrix follow a power law with spectral dimension 3. Yet, quantitative details of the statistics of the eigenvectors show subtle differences with homogeneous lattices and allow us to distinguish between topological and geometric sources of disorder in the network.

42 ENGINEERING↗
DOE OSTI2017

SiC-Based Composite Materials Obtained by Siliconizing Carbon Matrices

We have developed a method for fabrication of parts of complicated configuration from composite materials based on SiC ceramics, which employs the interaction of silicon melt with the carbon matrix having a certain composition and porosity. For elevating the operating temperatures of ceramic components, we have developed a method for depositing protective silicon-carbide coatings that is based on the interaction of the silicon melt and vapor with carbon obtained during thermal splitting of hydrocarbon molecules. The new structural ceramics are characterized by higher operating temperatures; chemical stability; mechanical strength; thermal shock, wear and radiation resistance; and parameters stability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Initial precipitation and hardening mechanism during non-isothermal aging in an Al–Mg–Si–Cu 6005A alloy

The characterization of precipitation and hardening mechanism during non-isothermal aging had been investigated using high resolution transmission electron microscopy for an Al–Mg–Si–Cu 6005A alloy. It was proposed that the needle-shaped β″ precipitates with a three-dimension coherency strain-field and an increased number density in the Al matrix provided the maximum strengthening effect for the Al–Mg–Si–Cu 6005A alloy. Simultaneously, it was also found that the formation and evolution of clusters in the early precipitation were associated with the vacancy binding energy, during which Si atoms played an important role in controlling the numbers density of Mg/Si co-clusters, and the excess Si atoms provided the increased number of nucleation sites for the subsequent precipitates to strengthen and improve the precipitation rate. Finally, based on the experimental observation and theoretical analysis, the precipitation sequence during the early precipitation in the Al–Mg–Si–Cu 6005A alloy was proposed as: supersaturated solid solution → Si-vacancy pairs, Mg-vacancy pairs and Mg clusters → Si clusters, and dissolution of Mg clusters → Mg atoms diffusion into the existing Si clusters → Mg/Si co-clusters → GP zone. - Highlights: • β″ precipitates provide the maximum strengthening effect for the 6005A alloy. • Si atoms play an important role in controlling the numbers of Mg/Si co-clusters. • The early aging sequence is deduced based on the solute-vacancy binding energy.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Water confinement in nanoporous silica materials

The influence of the surface polarity of cylindrical silica nanopores and the presence of Na{sup +} ions as compensating charges on the structure and dynamics of confined water has been investigated by molecular dynamics simulations. A comparison between three different matrixes has been included: a protonated nanopore (PP, with SiOH groups), a deprotonated material (DP, with negatively charged surface groups), and a compensated-charge framework (CC, with sodium cations compensating the negative surface charge). The structure of water inside the different pores shows significant differences in terms of layer organization and hydrogen bonding network. Inside the CC pore the innermost layer is lost to be replaced by a quasi bulk phase. The electrostatic field generated by the DP pore is felt from the surface to the centre of pore leading to a strong orientation of water molecules even in the central part of the pore. Water dynamics inside both the PP and DP pores shows significant differences with respect to the CC pore in which the sub-diffusive regime of water is lost for a superdiffusive regime.

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

Holography, Gravity and Condensed Matter

Over the five years of funding from this grant, I produced 26 publications. These include a book-long monograph on "Holographic Quantum Matter" that is currently in press with MIT press. The remainder were mostly published in Physical Review Letters, the Journal of High Energy Physics, Nature Physics, Classical and Quantum Gravity and Physical Review B. Over this period, the field of holography applied to condensed matter physics developed from a promising theoretical approach to a mature conceptual and practical edifice, whose ideas were realized in experiments. My own work played a central role in this development. In particular, in the final year of this grant, I co-authored two experimental papers in which ideas that I had developed in earlier years were shown to usefully describe transport in strongly correlated materials — these papers were published in Science and in the Proceedings of the National Academy of Sciences (obviously my contribution to these papers was theoretical). My theoretical work in this period developed several new directions of research that have proven to be influential. These include (i) The construction of highly inhomogeneous black hole event horizons, realizing disordered fixed points and describing new regimes of classical gravity, (ii) The conjecture of a bound on diffusivities that could underpin transport in strongly interacting media — an idea which may be proven in the near future and has turned out to be intimately connected to studies of quantum chaos in black holes and strongly correlated media, (iii) The characterization of new forms of hydrodynamic transport, e.g. with phase-disordered order parameters. These studies pertain to key open questions in our understanding of how non-quasiparticle, intrinsically strongly interacting systems can behave. In addition to the interface between holography and strongly interacting condensed matter systems, I made several advances on understanding the role of entanglement in quantum gravity. These included the first computation of holographic entanglement beyond the bulk classical limit as well understanding short distance entanglement in the emergent spacetime of the c=1 matrix quantum mechanics. The objective here is ultimately to understanding how a priori non-local degrees of freedom can re-arrange themselves quantum mechanically to support emergent local dynamics. Much of work funded by this grant involved collaboration with postdocs and graduate students, several of which were directly funded by the grant. These students have now successfully graduated to postdoctoral positions and in one case to high tech industry. The ideas developed in this work have directly fed into my current research in which I am aiming to prove fundamental bounds on entropy production and transport from quantum mechanics and statistical physics. As often, as with much of my previous work, black hole physics can be an inspiration for extreme dynamics such as fundamental bounds, but ultimately one hopes to prove them using more general tools of quantum field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Characterisation of nano-assemblies inside mesopores using neutron scattering

Adsorption of molecular and nanoscale matter in mesoporous materials is important in filtration and chromatography as well as membrane processes. However, the morphology and distribution of self-assembled structures of adsorbate formed within nanoconfined environments is largely unknown. Small Angle Neutron Scattering (SANS) with porous matrix matching the neutron scattering length density of the solvent has the potential to provide detailed information on the self-assemblies formed in pore spaces. However, analysis of such SANS profiles remains a challenge. In this study, we extend the SANS analysis method previously developed by Findenegg and co-workers to include interparticle correlations, providing a qualitative characterisation of the adsorption state of surfactants and nanoparticles in the cylindrical pores of silica nanomaterials. We find that the pore filling fraction and the self-assembled state of the adsorbate governs the scattering profile. We apply the model to two materials namely, triethyleneglycol monohexyl ether (C 6 E 3 ) surfactant and gold nanoparticles adsorbed in SBA-15 mesoporous silica. In contrast to the C 6 E 3 system, Bragg scattering dominates the diffuse scattering in gold nanoparticle system, which we attribute to the immobile internal structure of the nanoparticles. In conclusion, this new SANS analysis method has potential to elucidate structures of soft matter nanoassemblies inside mesopores.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Evaluating quantitative 3-D image analysis as a design tool for low enriched uranium fuel compacts for the transient reactor test facility: A preliminary study

In this study, 3-D image analysis when combined with a non-destructive examination technique such as X-ray computed tomography (CT) provides a highly quantitative tool for the investigation of a material’s structure. In this investigation 3-D image analysis and X-ray CT were combined to analyze the microstructure of a preliminary subsized fuel compact for the Transient Reactor Test Facility’s low enriched uranium conversion program to assess the feasibility of the combined techniques for use in the optimization of the fuel compact fabrication process. The quantitative image analysis focused on determining the size and spatial distribution of the surrogate fuel particles and the size, shape, and orientation of voids within the compact. Additionally, the maximum effect of microstructural features on heat transfer through the carbonaceous matrix of the preliminary compact was estimated. The surrogate fuel particles occupied 0.8% of the compact by volume with a log-normal distribution of particle sizes with a mean diameter of 39 μm and a standard deviation of 16 μm. Roughly 39% of the particles had a diameter greater than the specified maximum particle size of 44 μm suggesting that the particles agglomerate during fabrication. The local volume fraction of particles also varies significantly within the compact although uniformities appear to be evenly dispersed throughout the analysed volume. The voids produced during fabrication were on average plate-like in nature with their major axis oriented perpendicular to the compaction direction of the compact. Finally, the microstructure, mainly the large preferentially oriented voids, may cause a small degree of anisotropy in the thermal diffusivity within the compact. α∥/α⊥, the ratio of thermal diffusivities parallel to and perpendicular to the compaction direction are expected to be no less than 0.95 with an upper bound of 1.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Post-irradiation examination of the Sirius-1 nuclear thermal propulsion fuel test

Nuclear Thermal Propulsion (NTP) systems hold promise in reducing transit times for exploration of the solar system by both crewed and uncrewed missions. NTP systems currently under investigation include a once-through high-temperature gas-cooled fission reactor to provide thermal energy to heat the coolant, which also serves as the propellant. The fuel system consists of angular UN fuel particles dispersed in a matrix of W/Re, creating a ceramic and metallic composite or cermet, which has been irradiated in Idaho National Laboratory's (INL) Transient Reactor Test Facility (TREAT). This enabled evaluation of these materials under representative nuclear heating rates (~95 K/s) and peak surface temperatures (~2527 K). Further, these surface conditions were indicative of achievement of the target peak internal temperature of approximately 2850 K. These tests named Sirius-1 (UN–W/Re), have been irradiated and this paper will present post-irradiation examination results. The Sirius-1 test produced minor, stable cracks in the fuel specimen and spalling of surface material. Volatilized uranium ‘soot’ was found deposited on the inner wall of the irradiation capsule indicating loss of some fissile material from the fuel specimen. Spalling from the surfaces was also noted upon visual inspection. Uranium diffusion from the fuel particles resulted in the formation of U/Re phases and the production of a laminar microstructure on the edge of the fuel system. Overall, the specimen was stable and remained within a coolable geometry upon conclusion of multiple thermal cycles to prototypical operating temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Holistic comparison of environmental barrier coating material candidates through design of a figure of merit

The first figure of merit for environmental barrier coating (EBC) materials was designed through a ranking system for material properties pertaining to established EBC failure modes in service. Seven past, present, and novel EBC candidate materials were used in the figure of merit design: SiO 2 , Ba 0.75 Sr 0.25 Al 2 Si 2 O 8 (BSAS), HfSiO 4 , Yb 2 Si 2 O 7 , Yb 2 SiO 5 , Yb 2 O 3 , and YbPO 4 . Utilizing compiled data from the literature, the presented figure of merit verified Yb 2 Si 2 O 7 as the state-of-the-art candidate with optimal EBC properties and inferred that YbPO 4 should be considered as a potentially viable EBC material candidate. Further, the figure of merit allows for a holistic comparison of EBC candidates and informs experimental and computational search efforts for next-generation complex EBCs. Clear knowledge gaps found through this work include CaO–MgO–Al 2 O 3 –SiO 2 (CMAS)-resistant coatings, EBC lifetimes before delamination, and oxidant diffusion rates in relevant EBC microstructures. It was shown that while some materials show promise for solving a single key failure mode for EBCs (i.e., CMAS reactivity), a community-wide goal should be placed on materials development to achieve acceptable resistance against all major failure modes, which are interconnected. Novel compositionally complex EBC materials, in addition to layered EBC architectures, show promise for the optimization of material properties for long-lifetime EBCs in combustion environments.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Reassessment of uncertainties for the PU-MET-INTER-002 benchmark - 14432

Uncertainty for the PU-MET-INTER-002 benchmark of ZPR-6/10 was unusually large for a ZPR assembly and was dominated by the uncertainties in matrix tube pitch and composition uncertainties for plutonium, graphite and stainless steel plates. The original uncertainty analysis for PU-MET-INTER-002 was limited by computer resources, software and the nuclear data available in 2002, all of which have been improved since. An as-built Monte Carlo model of ZPR-6/10 was used to quantify components of the transformation from the detailed model to a simplified RZ benchmark in 2002. Compositional uncertainties were estimated based upon TWODANT calculations of the RZ benchmark model (20 group, S{sub 8} - P{sub 1}). MC{sup 2}-2 was used with END/B-V.2 data to generate homogenized cross section data for TWODANT. The MC{sup 2}-2 calculations did not account for lattice cell spatial heterogeneity which was driven by the manual process required to generate ∼80 modifications to the deterministic model and by the software limitations of MC{sup 2}-2. Improvements to DIF3D-VARIANT, development of the MC{sup 2}-3 and BuildZPRModel codes, and better computer resources make accurate 3D transport calculations of reactivity changes practical today. BuildZPRModel imports the as-built geometry and creates a 'drawer' homogenized DIF3D model where each radial DIF3D mesh is the exact size of a single ZPR matrix tube and drawer. A MC{sup 2}-3 input is constructed for each drawer to produce homogenized cross section data where the spatial cell heterogeneity is captured by solving a 1D MOC transport calculation on the ultrafine energy group level. Collapsing 3D drawers into 1D MC{sup 2}-3 drawer models follows conventional procedures and is done by BuildZPRModel. MCNP6 was used to re-compute relevant geometry components of the uncertainty using ENDF/B-VII.0 data with the as-built model. Changes in the geometry components compared with earlier results are insignificant. DIF3D calculations of the compositional uncertainties were performed using 33, 70, 116, 230, 425, and 703 group energy structures and selectively combined with P{sub 3}-P{sub 7} transport approximations and P{sub 1} and P{sub 3} scattering kernels. 3D diffusion theory calculations were also performed for comparison. Many composition uncertainty components were essentially unchanged with the updated methodology. However, the dominant experimental uncertainty for PU-MET-INTER-002 is the composition of stainless steel. The new calculations reduce the stainless steel uncertainty from ∼115 pcm to ∼70 pcm which significantly reduces the total uncertainty for PU-MET-INTER-002. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

The Schwinger-Keldysh coset construction

The coset construction is a tool for systematically building low energy effective actions for Nambu-Goldstone modes. This technique is typically used to compute time-ordered correlators appropriate for S-matrix computations for systems in their ground state. In this paper, we extend this technique to the Schwinger-Keldysh formalism, which enables one to calculate a wider variety of correlators and applies also to systems in a mixed state. We focus our attention on internal symmetries and demonstrate that, after identifying the appropriate symmetry breaking pattern, Schwinger-Keldysh effective actions for Nambu-Goldstone modes can be constructed using the standard rules of the coset construction. Particular emphasis is placed on the thermal state and ensuring that correlators satisfy the KMS relation. We also discuss explicitly the power counting scheme underlying our effective actions. We comment on the similarities and differences between our approach and others that have previously appeared in the literature. In particular, our prescription does not require the introduction of additional “diffusive” symmetries and retains the full non-linear structure generated by the coset construction. We conclude with a series of explicit examples, including a computation of the finite-temperature two-point functions of conserved spin currents in non-relativistic paramagnets, antiferromagnets, and ferromagnets. Along the way, we also clarify the discrete symmetries that set antiferromagnets apart from ferromagnets, and point out that the dynamical KMS symmetry must be implemented in different ways in these two systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Microstructure and mechanical properties of NiCoCrAlYTa alloy processed by press and sintering route

Nickel-based superalloys such as NiCoCrAlY are widely used in high-temperature applications, such as gas turbine components in the energy and aerospace industries, due to their strength, high elastic modulus, and high-temperature oxidation resistance. However, the processing of these alloys is complex and costly, and the alloys are currently used as a bond coat in thermal barrier coatings. In this work, the effect of cold press and sintering processing parameters on the microstructure and mechanical properties of NiCoCrAlY alloy were studied using the powder metallurgy route as a new way to obtain NiCoCrAlYTa samples from a gas atomized prealloyed powder feedstock. High mechanical strength and adequate densification up to 98% were achieved. The most suitable compaction pressure and sintering temperature were determined for NiCoCrAlYTa alloy through microstructure characterization. Scanning electron microscopy (SEM), transmission electron microscopy (TEM), X-ray diffraction (XRD), and energy dispersive spectroscopy microanalysis (EDS) were performed to confirm the expected γ-Ni matrix and β-NiAl phase distribution. Additionally, the results demonstrated the unexpected presence of carbides and Ni–Y-rich zones in the microstructure due to the powder metallurgy processing parameters used. Thus, microhardness, nanoindentation and uniaxial compression tests were conducted to correlate the microstructure of the alloy samples with their mechanical properties under the different studied conditions. The results show that the compaction pressure did not significantly affect the mechanical properties of the alloy samples. In this work, the compaction pressures of 400, 700 and 1000 MPa were used. The sintering temperature of 1200 °C for NiCoCrAlYTa alloy was preferred; above this temperature, the improvement in mechanical properties is not significant due to grain coarsening, whereas a lower temperature produces a decrease in mechanical properties due to high porosity and poor solid-state diffusion. - Graphical abstract: Display Omitted - Highlights: • We made NiCoCrAlYTa alloy by a conventional powder metallurgy route. • High densification and adequate strength were observed. • The presence of unexpected carbides found along γ/γ and γ/β grain boundaries was detected. • The effect of cold press and sintering processing parameters on the microstructure and mechanical properties were studied.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evidence for the existence and ecological relevance of fast-cycling mineral-associated organic matter

Longstanding theories and models classify mineral-associated organic matter as the large ( ~ 60%) but slow-cycling and persistent portion of soil organic matter. Strong physico-chemical interactions and diffusion limitations restrict the turnover of mineral-associated organic matter, allowing carbon and nitrogen bound therein to persist in soil for as long as centuries to millennia. However, mineral-associated organic matter is a chemically and functionally diverse pool with a substantial portion cycling at relatively fast (i.e., minutes to years) timescales. Despite a growing body of evidence for the heterogenous and multi-pool nature of mineral-associated organic matter, we lack consensus on how to conceptualize and directly quantify fast-cycling mineral-associated organic matter and its ecological significance. We demonstrate that the dynamic qualities of fast-cycling mineral-associated organic matter vary based on 1) the chemistry of the mineral particles and organic matter, 2) the complex set of interactions between organic matter and the mineral matrix, and 3) the presence and strength of destabilizing forces that lead to decomposition or loss of mineral-associated organic matter (i.e., plant-microbe interactions, agricultural intensification, and climate change). Finally, we discuss potential implications and research opportunities for how we measure, manage, and model the dynamic subfraction of this otherwise persistent pool of soil organic matter.

Biological and medical sciences↗
DOE OSTI2019

A phase-field model for hydride formation in polycrystalline metals: Application to δ-hydride in zirconium alloys

We report a phase-field model for simulating metal hydride formation involving large volume expansion in single- and polycrystals. As an example, we consider δ-hydride formation in α-zirconium (Zr), which involves both displacive crystallographic structural change and hydrogen diffusion process. Thermodynamic Gibbs energy functions are extracted from the available thermodynamic database based on the sublattice model for the interstitial solid solutions. Solute-grain boundary interactions and inhomogeneous elasticity of polycrystals are taken into consideration within the context of diffuse-interface description. The stress-free transformation strains of multiple variants for hcp -Zr (α) to fcc -hydride (δ) transformation are derived based on the well-established orientation relationship between the α and δ phases as well as the corresponding temperature-dependent lattice parameters. In particular, to account for the large volume expansion, we introduced the mixed interfacial coherency concept between those phases—basal planes are coherent and prismatic planes are semi-coherent in computing the strain energy contribution to the thermodynamics. We analyzed the morphological characteristics of hydrides involving multiple structural variants and their interactions with grain boundaries. Moreover, our simulation study allows for the exploration of the possible hydride re-orientation mechanisms when precipitating under applied tensile load, taking into account the variation in the interfacial coherency between hydrides and matrix, their elastic interactions with the applied stress, as well as their morphology-dependent interactions with grain boundaries. Here, the phase-field model presented here is generally applicable to hydride formation in any binary metal–hydrogen systems.

36 MATERIALS SCIENCE↗
DOE OSTI2016

On the Degradation of the Effectiveness of Nonlinear Diffusion Acceleration with Parallel Block Jacobi Splitting

MOOSE's radiation transport module Rattlesnake comprises a first order S N solver that can be accelerated with a flexible Nonlinear Diffusion Acceleration. At each NDA iteration, the streaming and collision operator of the S{sub N} equations are (approximately) inverted using a transport sweep. When executed on multiple processors parallel block Jacobi (PBJ) splitting is utilized lagging the processor sub-domain's incoming angular fluxes making the implemented method asynchronous, i.e. the convergence properties of NDA become dependent on the number of concurrent tasks and details of the partitioning. Within this work, we study the NDA convergence properties for the two-dimensional C5G7 benchmark problem PBJ with p = 1, 2, 4, 8, 12, 16, 24, 32, 48, 60, 72, 96 and 144 concurrent tasks. Rosa studied the convergence properties of PBJ splitting methods with and without transport synthetic acceleration (TSA) in slab geometry and two-dimensional Cartesian geometry, respectively. Rosa's iteration scheme lags both the incoming angular fluxes and the scattering source making it very similar to the method described here. Rosa found by means of Fourier analysis that convergence of unaccelerated PBJ can degrade for optically thin subdomains, i.e. for a large number of processors given a fixed domain size. Depending on the choice of parameters application of TSA can lead to an effective algorithm. In contrast to Rosa's work, this work focuses on PBJ splitting for NDA accelerated solution of the multigroup eigenvalue problem. Later, Zerr and Anistratov used their integral transport matrix method in conjunction with PBJ. This iteration scheme is different as it does not lag the scattering source. Nevertheless, Fourier analysis with Diamond Difference spatial discretization showed that the method degrades with increasing the overall size of the spatial domain while keeping the sub-domain size constant. In the case of an infinite medium, the observed spectral radius becomes unity. These references establish the iterative performance of PBJ for one-group fixed source problems discretized on rectangular cells, but for the Rattlesnake first order NDA method, an emphasis lies on multigroup eigenvalue problems. In this setting, the additional complication arises that the low order equations are not exactly inverted, but only to some iterative criterion. Further, as the low order problem is solved using Jacobian-Free Newton Krylov methods, a suitable preconditioner for the GMRES iterations is required. As the effectiveness of the preconditioners may decrease when more concurrent processes are used, we expect the low order solve to impact the overall effectiveness of the NDA method in parallel execution as well. Therefore, this work focuses on the question of how detrimental the asynchronicity of the high and low order iterative solutions is for the effectiveness of NDA solutions of realistic, unstructured grid, multigroup eigenvalue problems. The iterative convergence properties of Rattlesnake's NDA method for multigroup eigenvalue problems are studied for varying number of concurrent processes. Rattlesnake's NDA method is asynchronous because the high order transport sweep uses PBJ splitting and because the low order solve's preconditioners, algebraic multigrid and additive Schwarz, are asynchronous. We found that up to 144 processes the number of Picard iterations increases by a factor of 2-3 affecting but not destroying NDA's parallel performance for the C5G7 benchmark problem, a representative LWR (light water reactor) configuration. The solution of the NDA low order eigenvalue problem posed greater problems than expected for the solvers currently available in Rattlesnake mainly because of problems to effectively precondition the modified low-order problem. As a matter of fact, the low order problem becomes an advection-diffusion problem. Non-robustness was observed for the AMG preconditioned coarse mesh acceleration for more than 60 concurrent processes, which was fixed by increasing the number of permissible V-cycles per AMG call. For fine mesh acceleration AMG proved to be ineffective and we resorted to ASM. For both AMG and ASM, increasing number of linear iterations per nonlinear are observed when increasing the number of processes, but for ASM this problem appears to be worse. The results from this work also indicate that initial development efforts should be diverted from implementation of a high-order parallel sweeping algorithm to the effective and robust iterative solution of the low order NDA multigroup eigenvalue problem. (authors)

97 MATHEMATICS AND COMPUTING↗