Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex solids”

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 667 records · Page 37

Effects of Mineral Impurities and Heat on Uranium(VI) Sorption onto Bentonite (Final Report)

Most nuclear waste disposal options include compacted bentonite, consisting of montmorillonite clay, as part of a barrier system to minimize contaminant mobility. Uranium (U) is the primary element in spent nuclear fuel, and from an environmental perspective, a potential contaminant of water resources. Furthermore, uranium may control nuclear fuel degradation rates and the consequent release of other radioactive contaminants based on its diffusive transport away from waste canisters. Uranium sorption onto clay and other mineral surfaces is expected to limit U(VI) mobility in these systems. However, at this point a prediction of U(VI) sorption and transport behavior in performance assessment (PA) models is complicated by a series of factors, such as: (1) the presence of various U(VI) solution species with different charges and sorption characteristics; (2) the complex microstructure of montmorillonite clay resulting in two types of clay porosities and multiple sorption sites; (3) the largely unknown effects of bentonite mineral impurities on pore water chemistry, U(VI) sorption and diffusion behavior; and (4) the potential impacts of heat, generated by the decay of spent fuel, on mineralogical and microstructural transformations, and any subsequent effects on radionuclide sorption and mobility. In this project, we investigated the effects of calcite impurities on U(VI) sorption and diffusion onto montmorillonite before and after mineral exposure to heat. Overall, our main goals were to: 1. Gain a mechanistic understanding of the effects of calcite impurities on U(VI) sorption and diffusion in montmorillonite/bentonite at variable chemical conditions, and before and after mineral exposure to heat; and to 2. Work towards the development of new surface complexation and diffusion models to predict U(VI) sorption and mobility in montmorillonite/bentonite while accounting for the impacts of calcite impurities and heat. In order to achieve these goals, we pursued the following research objectives and tasks: 1. Characterize the effects of calcite impurities on U(VI) sorption onto Na-montmorillonite, and evaluate if and how these effects change due to solids’ exposure to heat (Task 1). 2. Characterize U(VI) surface speciation in solid samples from Task 1, with a particular focus on U(VI) precipitates and co-precipitates (Task 2). 3. Test the relevance of observed changes in U(VI) sorption characteristics for U(VI) diffusive transport behavior in bentonite barriers (Task 3). 4. Evaluate the exclusion of anionic U(VI) solution complexes from clay nanopores (Task 4). 5. Building on Tasks 1-4, start the development of new surface complexation and reactive transport models that can simulate and predict U(VI) sorption and diffusion in complex bentonite systems (Task 5). A detailed description of all tasks, and their specific results and conclusions is provided in this report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Towards the Stable Chelation of Radioantimony(V) for Targeted Auger Theranostics

Antimony-119 ( 119 Sb) is one of the most attractive Auger-electron emitters identified to date, but it remains practically unexplored for targeted radiotherapy because no chelators have been identified to stably bind this metalloid in vivo. In a departure from current studies focused on chelator development for Sb(III), we explore the chelation chemistry of Sb(V) using the tris-catecholate ligand TREN-CAM. Through a combination of radiolabeling, spectroscopic, solid-state, and computational studies, the radiochemistry and structural chemistry of TREN-CAM with 1XX/nat Sb(V) were established. The resulting [ 1XX Sb]Sb–TREN-CAM complex remained intact for several days in human serum, signifying high stability under biological conditions. Finally, the first in vivo single photon emission computed tomography and positron emission tomography imaging studies were carried out using 117 Sb, the diagnostic analogue of 119 Sb. These studies revealed marked differences in the uptake and distribution of activity in mice administered unchelated [ 117 Sb]Sb(OH) 6 – versus [ 117 Sb]Sb–TREN-CAM, suggesting that 117 Sb is largely retained by TREN-CAM over the time course of the study. Collectively, these findings demonstrate the most physiologically stable complex of no-carrier-added 1XX Sb yet reported, offering new promise for the clinical implementation of radioantimony in nuclear medicine. In conclusion, our results also establish the feasibility of 117 Sb as an elementally matched partner to 119 Sb for theranostic applications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Stress interference in multilayer additive friction stir deposition of AA6061 aluminum

Due to the multilayer deposition nature of metal additive manufacturing processes, each layer being printed experiences the state of thermokinetic and thermomechanical stress that in turn interfere with the state of thermokinetics and thermomechanical stress of subsequently deposited layers. Especially, this multilayer interference significantly affects the resultant properties of the component fabricated using solid-state additive friction stir deposition due to evolution of asymmetric state of planar stress. Due to the lack of comprehensive and suitable in situ diagnosis technique, the complex interference of inter- and multi-layer stresses during additive friction stir deposition was studied in an integrated approach of numerical simulation of fluidic state and experimental probing of stress influenced ultrasonic elastography. The uni-directional and bi-directional layer deposition configurations adopted during additive friction stir deposition result in the generation of constructive and destructive interference of the interlayer stress and hence, asymmetric and symmetric dynamic elasticity distribution respectively within the subsequent layers. With subsequent deposition of additional layers, the odd and even numbers of deposited layers generate asymmetric and nearly symmetric dynamic elasticity distributions.

Yang, Teng↗

Pyrrophens: Pyrrole-Based Hexadentate Ligands Tailor-Made for Uranyl (UO 2 2+ ) Coordination and Molecular Recognition

Derivatives of a novel pyrrole-containing Schiff base ligand system (called “pyrrophen”) are presented which feature substituted phenylene linkers (R 1 = R 2 = H (H 2 L 1 ); R 1 = R 2 = CH 3 (H 2 L 2 )) and a binding pocket modeled after macrocyclic species. These ligands bind neutral CH 3 OH in the solid state through pyrrolic hydrogen-bonding. Here, the interaction of the uranyl cation (UO 2 2+ ) and H 2 L 1–2 yields planar hexagonal bipyramdial uranyl complexes, while the Cu 2+ and Zn 2+ complexes were found to self-assemble as dinuclear helicate complexes (M 2 L 2 ) with H 2 L 1 under identical conditions. The favorable binding of UO 2 2+ over Zn 2+ provides insight into the molecular recognition of uranyl over other metal species. Structural features of these complexes are examined with special attention to features of the UO 2 2+ coordination environment which distinguish them from other related salophen and porphyrinoid complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Breaking the Neptunyl Barrier: Direct Access to Neptunium(IV) in Aqueous Solution via Polyoxometalate-Mediated Reduction and Stabilization

Neptunium exhibits truly unique chemistry as its speciation is dominated by the neptunyl(V) ion (NpO 2 + ). Here, in this study, we describe the spontaneous destabilization and reduction of neptunyl(V) via complexation to the Keggin-type polyoxometalate (POM) ligand PW 11 O 39 7– . The POM-mediated reduction of NpO 2 + does not require any reducing agent and occurs within minutes, at room temperature, and in aqueous solution. The resulting [Np(PW 11 O 39 ) 2 ] 10- complex (Np(PW 11 ) 2 ) remains soluble, water-stable, and air-stable for weeks and persists over an extended acidity range. Single-crystal structure, solid-state Raman and UV–visible absorbance characterization of Np(PW 11 ) 2 revealed a mixed alpha/beta isomery of the Keggin ion, forming an unprecedented 50:50 mixture of Np(α-PW 11 ) 2 and Np(α-PW 11 )(β-PW 11 ). Experiments with other tetravalent ions (i.e., Zr 4+ , Hf 4+ , Ce 4+ , and Th 4+ ) indicate that the occurrence of the beta isomer is specific to Np 4+ and independent of the cation’s size. Solution-state characterization of the Np-PW 11 system via UV-visible-NIR absorbance, 31 P NMR, VT NMR, and relaxometry further elucidated the speciation. Moreover, comparative experiments with uranium revealed that the two types of actinyl ions (UO 2 2+ vs NpO 2 + ) undergo drastically different reactions in the presence of PW 11 . UO 2 2+ is not reduced but instead uses PW 11 as a phosphate reservoir and precipitates as uranyl phosphate.

Hastings, Ashley M. [Lawrence Livermore National L↗

Computational alchemy clarifies origins of alloy strengthening

Solid solution strengthening (SSS) is widely used to enhance mechanical properties of metals. Originally developed for dilute alloys, classical SSS theories are presently challenged by the rise of complex concentrated alloys (CCA) with nearly equiatomic compositions. Here, we propose and develop a method of “computational alchemy” in which interatomic interactions are modified to systematically vary two key physical parameters defining SSS - atomic size misfit and elastic stiffness misfit - over a maximally wide range of two misfits. The resulting alchemical alloys are subjected to massive (~10 8 atoms) molecular dynamics (MD) simulations reproducing full complexity of plastic strength response. At variance with prevailing views, stiffness misfit is observed to contribute to SSS on par if not more than size misfit. Furthermore, depending on exactly how two misfits are combined, they result in synergistic (amplification) or antagonistic (compensation) effect on alloy strengthening. Unlike real CCAs in which each component element comes with its own specific size and stiffness, our alchemical model alloys span the space of two misfits continuously revealing trends in alloy strengthening unrecognized so far. Our study demonstrates unique value of intentionally unrealistic models for gaining deep physical insights into material behaviors that are difficult to reveal otherwise.

36 MATERIALS SCIENCE↗

Dynamic experiments to study the α - ε phase transition in cerium

The ability to understand and predict the response of matter at extreme conditions requires knowledge of a material’s equation-of-state including the location of phase boundaries, transition kinetics, and the evolution of material strength. Cerium is a material with a complex phase diagram that continues to attract significant scientific interest. Recent dynamic experiments have provided information on the low-pressure γ–α phase transition, sound speed, and Hugoniot data for the higher-pressure α phase, as well as the incipient shock melt transition. Despite these efforts, there are still regions of the phase diagram that are largely unexplored dynamically, including the high-pressure region below the melt boundary. Along a room temperature isotherm, diamond anvil cell data report a transition to the ϵ phase between 13 and 17 GPa. At higher temperatures, similar diamond anvil cell data show significant disagreement regarding the existence, location, and slope of the ϵ-phase boundary. In this work, double-shock loading was used to access the α–ϵ region of the phase diagram to obtain equation-of-state information and to determine the location of the ϵ-phase boundary for shock loading.

74 ATOMIC AND MOLECULAR PHYSICS↗

A hybrid machine-learning approach for analysis of methane hydrate formation dynamics in porous media with synchrotron CT imaging

Fast multi-phase processes in methane hydrate bearing samples pose a challenge for quantitative micro-computed tomography study and experiment steering due to complex tomographic data analysis involving time-consuming segmentation procedures. This is because of the sample's multi-scale structure, which changes over time, low contrast between solid and fluid materials, and the large amount of data acquired during dynamic processes. Here, a hybrid approach is proposed for the automatic segmentation of tomographic data from time-resolved imaging of methane gas-hydrate formation in sandy granular media, which includes a deep-learning 3D U-Net model. To prepare a training dataset for the 3D U-Net, a technique to automate data labeling based on sample-specific information about the mineral matrix immobility and occasional fluid movement in pores is proposed. Automatic segmentation allowed for studying properties of the hydrate growth in pores, as well as dynamic processes such as incremental flow and redistribution of pore brine. Results of the quantitative analysis showed that for typical gas-hydrate stability parameters (100 bar methane pressure, 7°C temperature) the rate of formation is slow (less than 1% per hour), after which the surface area of contact between brine and gas increases, resulting in faster formation (2.5% per hour). Hydrate growth reaches the saturation point after 11 h of the experiment. Finally, the efficacy of the proposed segmentation scheme in on-the-fly automatic data analysis and experiment steering with zooming to regions of interest is demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Geant4 RApid Pair Production Application (GRAPPA) v1.0.0

GRAPPA is a Geant4 application that simulates the interaction of an incoming particle beam (typically electron or photons) with an High-Z material. Particles are generated at a certain energy and interact with a solid target. The product of this interaction is detected and saved on file. The application implements a basic target that is just a simple box, and it accepts complex target descriptions via GDML files. The user has the freedom to set detectors to capture the desired particles. This application enables studying processes such as antiparticle production, muon generation, and atomic decay.

Terzani, Davide [Lawrence Berkeley National Labora↗

Development of Additively Manufactured Complex Tools for Autoclave Cure Composites

IACMI Project 4.9, Tooling for Composites with Washout Features Produced by Additive Manufacturing, assembled a team including the industry lead, Ability Composites, NREL and Colorado State University (CSU). Ability Composites had originally expressed interest in alternate methods of producing tooling for composite parts. In follow-up discussions, it became clear that one of the bigger tooling challenges revolved around small production volume composite parts that were tooled on washout material due to the complex geometry. To build an understanding of the potential, both from a technology and a cost perspective, for replacing conventional washout tooling with 3D printed thermoplastic tooling, a number of commercially available dissolvable FDM printing materials were evaluated, leading to tooling representative of commercial articles of interest to Ability Composites. Ultimately, Ability Composites was able to directly compare autoclave processed prepreg composite parts produced on conventional washout tooling to composite parts molded on 3D printed dissolvable tooling produced at CSU. Small, laboratory test specimens were developed to investigate the structural performance of the candidate materials under autoclave processing conditions, which were nominally 121 °C (250 °F) and 345 kPa (50 psi). In addition, several internal structural configurations (infills) were evaluated under autoclave conditions using model materials. The results of these tests indicated that two materials, Stratasys ST 130 and Infinite Materials Solutions Aquasys 180 (AQ 180), were the best candidates, given the specified autoclave processing conditions. ST-130 was slightly more robust than AQ-180; however, the AQ-180 was carried forward as it was dissolvable in water, not requiring the basic solution needed to dissolve ST-130. Based on the preliminary material and 3D printed structures evaluations, larger tools with a truncated square pyramid geometry were created to produce prepreg composite test articles for 3D printed dissolvable tool evaluation under standard autoclave fabrication conditions. Two tools were manufactured using ST-130 and one tool using traditional ceramic washout tooling media. The tools were evaluated for geometric fidelity and surface roughness changes before and after carbon fiber/epoxy prepreg composites were manufactured on the tooling. The autoclave processing did not impact the geometry significantly and was completed at 121 °C and 345 kPa, indicating satisfactory tool performance. The results from surface roughness testing of both the resulting composite and the associated tooling indicated that an adequate surface resulted without the need for a surface sealing step, as was required for the conventional washout tooling. Based on results of the truncated pyramid tests as a basis, ST-130, AQ-120 and AQ-180 materials were carried forward to the tool geometry of interest to Ability Composites. These hollow rectangular bent ducts, which were complex in nature and not extractable after cure, were used to understand the impacts of tool material and thickness. One ST-130 tool was produced as a partially solid part, with an enclosed 40% dense infill region to reduce weight and material use. This was the same approach evaluated in the truncated pyramid portion of the study. This tool was to be envelope vacuum bagged and directly compared to a monolithic tool of conventional washout material. The traditional monolithic ceramic tool was manufactured by Ability composites using CNC-based subtractive methods. An additional five dissolvable polymer tools, manufactured from ST-130, AQ-120, and AQ-180, using a hollow design were 3D printed and used to produce carbon fiber/epoxy composite evaluation articles. These hollow dissolvable tools were expected to be less influenced by the autoclave conditions as the wall was solid and vacuum bagging was inside and outside the tool. This alternative geometry was also evaluated as an option in techno-economic modeling. Print times were reduced from in excess of 3 days to under 30 hours, while surface quality and and tool integrity were substantially improved in the transition from the partially solid tool to the hollow tooling concept. Ability Composites produced autoclave-cured prepreg ducts on each of these tools. The autoclave conditions utilized were more severe than those of the initial trials, reaching temperatures of 160 °C and a pressure of 414 kPa. Under these conditions, the partially solid 3D printed tool with skin and 40% dense infill crushed significantly; however, the thicker ST-130 hollow tool showed good promise, deforming only slightly. The thinner hollow tool walls were unsuccessful as were the other materials. Overall, the hollow tool manufacturing process saved significant amounts of time and material in manufacturing as compared to the solid ducts and produced composite surface quality improvements compared to the traditional washout tooling. The TEM was developed to allow direct comparisons between conventional washout tool manufacturing processes and those developed at CSU. It also allowed for two separate 3D printed tool geometries to be analyzed and compared. In this case, the square bent duct tool geometry was determined to be representative of common washout tools. This geometry was compared with a scaled-up version of it to assess differences in the two manufacturing processes based on tool size. The model was developed to make use of user input in the form of geometry details, process steps, manufacturing parameters, bulk material costs, capital equipment costs, and general costs to calculate overall labor, material, capital equipment, and energy costs per manufactured tool for the conventional and additive manufacturing processes for the two representative geometries. It was also able to estimate step-by-step process times for the manufacturing process and geometries. Based on significant input from Ability Composites and CSU from their knowledge gained from hands-on manufacturing of the 3D printed bent duct tool geometry, costs and process times were calculated for the two manufacturing processes. Results showed that the additive manufacturing techniques developed at CSU can substantially reduce the costs of tool manufacturing by reducing labor times and material usage. This is because additive manufacturing is a relatively hands-off process and allows for the tool design to be optimized to reduce material usage. The disadvantage, however, is that process times for additive manufacturing are significantly longer. The three-dimensional (3D) printing process is slow if tight tolerances are required, but the analysis did show that print times could be reduced with the hollow tool geometry. Also, further advances in additive manufacturing could expedite the process. Costs and process times for the tool washout process were calculated separately. They showed that costs are relatively insignificant when compared to the overall tool manufacturing processes, but with increases in tool size, costs for the conventional manufacturing approach are larger than for additive manufacturing. Again, the washout process for conventional tools is very hands-on, whereas for additively manufactured tools the print medium is dissolved in an automated detergent bath at the sacrifice of process time. The analysis showed that optimizing the additively manufactured tools may also reduce washout times. Overall, Project 4.9 demonstrated that commercially available dissolvable 3D printing materials exist that can be used to produce dissolvable tooling capable of surviving prepreg composites fabrication under autoclave conditions of 121 °C (250 °F) and 345 kPa (50 psi). An alternative hollow dissolvable tool design was developed which was structurally superior to the initial concept and was cost and time effective versus conventional washout tooling. The 3D printed sacrificial tool required no added surface sealing steps prior to composite part layup and cure, offering a significant advantage over the porous conventional washout tooling.

36 MATERIALS SCIENCE↗

A New Complex Borohydride LiAl(BH4)2Cl2

A new mixed alkali metal–aluminum borohydride LiAl(BH4)2Cl2 has been prepared via mechanochemical synthesis from the 2LiBH4–AlCl3 mixture. Structural characterization, performed using a combination of X-ray powder diffraction and solid-state NMR methods, indicates that the LiAl(BH4)2Cl2 phase adopts a unique 3D framework and crystallizes in an orthorhombic structure with the space group C2221, a = 11.6709(6) Å, b = 8.4718(4) Å, c = 7.5114(3) Å. The material shows excellent dehydrogenation characteristics, where hydrogen evolution starts at Tons = 70 °C, releasing approximately 2 wt.% of nearly pure (99.8 vol.%) hydrogen and a very small amount (~0.2 vol.%) of diborane. When compared to halide-free mixed alkali metal–aluminum borohydrides, the presence of Al‒Cl bonding in the LiAl(BH4)2Cl2 structure likely prevents the formation of Al(BH4)3 upon decomposition, thus suppressing the formation of diborane.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Defect Thermodynamics and Transport Properties of Proton Conducting Oxide BaZr 1–x Y x O 3–δ (x ≤ 0.1) Guided by Density Functional Theory Modeling

Density functional theory-based thermodynamic modeling was performed to determine the effect of humidity and H 2 /O 2 gas pressure on the defect chemistry and transport properties of proton conducting oxide BaZr 1-x Y x O 3-δ (x ≤ 0.1) in the temperature range of 800–1200K, relevant for solid oxide fuel/electrolysis cell applications. The first-principles charge defect analysis was carried out to obtain the defect energetic ensembles as a function of Fermi level as well as the dependence of the hydration and oxidation reaction energies on complex defect-dopant configurations. It is shown that oxygen vacancies introduced to compensate for the extra Ba vacancies may cause a slight increase in the concentration of hydroxyl and proton species upon hydration of Ba deficient BaZr 1-x Y x O 3-δ , while the overall proton diffusivity is predicted to decrease with increased proton hopping barriers upon trapping near the Ba vacancies. The developed defect model is further demonstrated to be able to describe the bulk defect chemistry and transport properties of BaZr 0.9 Y 0.1 O 3-δ under the solid oxide cell operating conditions. Finally, the theoretical framework developed in this work further allows the inclusion of configurational, energetic and electronic characteristics of the defects that can be used as a complementary tool to experimental measurements for testing various mechanistic pathways or to provide essential mechanistic data.

36 MATERIALS SCIENCE↗

Point containment algorithms for constructive solid geometry with unbounded primitives

Here, we present several algorithms for evaluating point containment in constructive solid geometry (CSG) trees with unbounded primitives. Three algorithms are presented based on postfix, prefix, and infix notations of the CSG binary expression tree. We show that prefix and infix notations enable short-circuiting logic, which reduces the number of primitives that must be checked during point containment. To evaluate the performance of the algorithms, each algorithm was implemented in the OpenMC Monte Carlo particle transport code, which relies on CSG to represent solid bodies through which subatomic particles travel. Two sets of tests were carried out. First, the execution time to generate a rasterized image of a 2D slice of three CSG models of varying complexity was measured. Use of both prefix and infix notations offered significant speedup over the postfix notation that has traditionally been used in particle transport codes, with infix resulting in a 6 x reduction in execution time relative to postfix for a model of a tokamak fusion device. We then measured the execution time of neutron transport simulations of the same three models using each of the algorithms. The results and performance improvements reveal the same trends as for the rasterization test, with a 5.52 x overall speedup using the infix notation relative to the original postfix notation in OpenMC for the tokamak model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Developing Lanthanide-Nitrate Cluster Chemistry toward Rare Earth Separations

Nitrate-decorated hexamers with a [Ln 6 (μ 6 -O)(μ 3 –OH) 8 ] 8+ core have been reported for nearly every lanthanide ion and are used as precursors for the assembly of functional metal–organic frameworks. Yet, few studies have examined the correlation between the solution and solid-state species, and the formation of mixed-metal clusters. Toward this end, a series of homo- and heterometal lanthanide nitrate hexamers was prepared via pH adjustment of aqueous lanthanide nitrate solutions. Examination of the homometallic europium solutions using Small Angle X-ray Scattering and nESI-MS showed that lower order complexes dominate lanthanide speciation in nitrate media. Yet, powder X-ray diffraction data of the precipitated phase confirmed the formation of [Ln 6 (μ 6 -O)(μ 3 -OH) 8 (NO 3 ) 6 (H 2 O) 12 ]·2(NO 3 )·n(H 2 O), Ln 6 , for Ln = Eu and Tb. For heterometal systems, analysis of the solid-state product by ICP–MS showed the selective incorporation of the heavier rare earths into Ln 6 . Selectivity was quantified by calculating an average separation factor, which is defined as the ratio of recovery factors of both metals. Further examination of the luminescence behavior of mixed metal [Tb 6–x Eu x (μ 6 -O)(μ 3 -OH) 8 (NO 3 ) 6 (H 2 O) 12 ]·2(NO 3 )·n(H 2 O), with x = 1.1–3.6, showed that the relative intensities of the peaks at 489 nm (terbium, 5 D 4 → 7 F 6 ) and 690 nm (europium, 5 D 0 → 7 F 4 ) trend with the percent incorporation of europium and terbium into the cluster.

anions↗

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗

Potential-pH diagrams considering complex oxide solution phases for understanding aqueous corrosion of multi-principal element alloys

The potential-pH diagram, a graphical representation of the thermodynamically predominant reaction products in aqueous corrosion, is originally proposed for the corrosion of pure metals. The original approach only leads to stoichiometric oxides and hydroxides as the oxidation products. However, numerous experiments show that non-stoichiometric oxide scales are prevalent in the aqueous corrosion of alloys. In the present study, a room temperature potential-pH diagram considering oxide solid solutions, as a generalization of the traditional potential-pH diagram with stoichiometric oxides, is constructed for an FCC single-phase multi-principal element alloy (MPEA) based on the CALculation of PHAse Diagram method. The predominant reaction products, the ions in aqueous solution, and the cation distribution in oxides are predicted. The oxide solid solution is stabilized by the mixing free energy (or mixing entropy) and the stabilizing effect becomes more significant as the temperature increases. Consequently, solid solution oxides are stable in large regions of the potential-pH diagram and the mixing free energy mostly affects the equilibrium composition of the stable oxides, while the shape of stable regions for oxides is mostly determined by the structure of the stable oxides. Agreements are found for Ni 2+ , Fe 2+ , and Mn 2+ between the atomic emission spectroelectrochemistry measurements and thermodynamic calculations, while deviations exist for Cr 3+ and Co 2+ possibly due to surface complexation with species such as Cl - and the oxide dissolution. By incorporating the solution models of oxides, the current work presents a general and more accurate way to analyze the reaction products during aqueous corrosion of MPEAs.

Materials Science↗

Ion Transport in Batteries with Polymer Electrolytes

We discuss polymer electrolytes for use in rechargeable lithium batteries. Polymer electrolytes have the potential to enable batteries with lithium metal anodes. These batteries have significantly higher theoretical energy densities than current lithium-ion batteries. We consider binary mixtures of polymers and salts. We also cover more complex systems such as polymer electrolytes swollen with a solvent (gel polymer electrolytes) and microphase separated polymer electrolytes. By covalently attaching the anions to the chains in a polymer solid, one obtains a single-ion conductor. We mainly focus on experiments wherein the polymer electrolyte is placed between two lithium metal electrodes. These experiments enable the determination of three transport parameters, ionic conductivity, salt diffusion coefficient, and transference number, and the thermodynamic factor. The properties of dry polymer electrolytes are contrasted with those of gel polymer electrolytes. The gel systems exhibit higher conductivity while the dry systems exhibit superior mechanical properties. We discuss interfacial impedance when lithium metal is contacted with polymer electrolytes and the importance of coulombic efficiency.

Yu, Xiaopeng↗

A hybrid smoothed-particle hydrodynamics model of oxide skins on molten aluminum

A computational model of aluminum melting is proposed which captures both the thermal fluid-solid phase transition and the mechanical effects of oxidation. The model hybridizes ideas from smoothed particle hydrodynamics and bonded particle models to simulate both hydrodynamic flows and solid elasticity. Oxidation is represented by dynamically adding and deleting spring-like bonds between surface fluid particles to represent the formation and rupture of the oxide skin. Various complex systems are simulated to demonstrate the adaptability of the method and to illustrate the significant impact of skin properties on material flow. As a result, initial comparison to experiments of a melting aluminum cantilever highlights that the computational model can reproduce key qualitative features of aluminum relocation.

Clemmer, Joel Thomas↗