Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regular hexagon”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Vertex-Aligned Model for Packing 4-Hexagonal Clusters in a Regular Hexagonal Container

This paper deals with a problem the packing polyhex clusters in a regular hexagonal container. It is a common problem in many applications with various cluster shapes used, but symmetric polyhex is the most useful in engineering due to its geometrical properties. Hence, we concentrate on mathematical modeling in such an application, where using the “bee” tetrahex is chosen for the new Compact Muon Solenoid (CMS) design upgrade, which is one of four detectors used in Large Hadron Collider (LHC) experiment at European Laboratory for Particle Physics (CERN). We start from the existing hexagonal containers with hexagonal cells packed inside, and uniform clustering applied. We compare the center-aligned (CA) and vertex-aligned (VA) models, analyzing cluster rotations providing the increased packing efficiency. We formally describe the geometrical properties of clustering approaches and show that cluster sharing is inevitable at the container border with uniform clustering. In addition, we propose a new vertex-aligned model decreasing the number of shared clusters in the uniform scenario, but with a smaller number of clusters contained inside the container. Also, we describe a non-uniform tetrahex cluster packing scheme in the proposed container model. With the proposed cluster packing solution, it is accomplished that all clusters are contained inside the container region. Since cluster-sharing is completely avoided at the container border, the maximal packing efficiency is obtained compared to the existing models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Optimizing Fast Charging and Wetting in Lithium-Ion Batteries with Optimal Microstructure Patterns Identified by Genetic Algorithm

To sustain the high-rate current required for fast charging electric vehicle batteries, electrodes must exhibit sufficiently high effective ionic diffusion. Additionally, to reduce battery manufacturing costs, wetting time must decrease. Both of these issues can be addressed by structuring the electrodes with mesoscale pore channels. However, their optimal spatial distribution, or patterns, is unknown. Herein, a genetic algorithm has been developed to identify these optimal patterns using a CPU-cheap proxy distance-based model to evaluate the impact of the added pore networks. Both coin-cell and pouch cell form factors have been considered for the wetting analysis, with their respective electrolyte infiltration mode. Regular hexagonal and mud-crack-like patterns, respectively, for fast charging and fast wetting were found to be optimal and have been compared with pre-determined, easier to manufacture, patterns. The model predicts that using cylindrical channels arranged in a regular hexagonal pattern is ∼6.25 times more efficient for fast charging as compared to grooved lines with both structuring strategies being restricted to a 5% electrode total volume loss. The model also shows that only a very limited electrode volume loss (1%–2%) is required to dramatically improve the wetting (5–20 times) compared to an unstructured electrode.

25 ENERGY STORAGE↗
DOE OSTI2024

Hexagonal Lattices of HIV Capsid Proteins Explored by Simulations Based on a Thermodynamically Consistent Model

HIV capsid proteins (CAs) may self-assemble into a variety of shapes under in vivo and in vitro conditions. Here, we employed simulations based on a residue-level coarse-grained (CG) model with full conformational flexibility to investigate hexagonal lattices, which are the underlying structural pattern for CA aggregations. Facilitated by enhanced sampling simulations to rigorously calculate CA dimerization and polymerization affinities, we calibrated our model to reproduce the experimentally measured affinities. Here using the calibrated model, we performed unbiased simulations on several large systems consisting of 1512 CA subunits, allowing reversible binding and unbinding of the CAs in a thermodynamically consistent manner. In one simulation, a preassembled hexagonal CA sheet developed spontaneous curvatures reminiscent of those observed in experiments, and the edges of the sheet exhibited local curvatures larger than those of the interior. In other simulations starting with randomly distributed CAs at different concentrations, existing CA assemblies grew by binding free capsomeres to the edges and by merging with other assemblies. At high CA concentrations, rapid establishment of predominant aggregates was followed by much slower adjustments toward more regular hexagonal lattices, with increasing numbers of intact CA hexamers and pentamers being formed. Our approach of adapting a general CG model to specific systems by using experimental binding data represents a practical and effective strategy for simulating and elucidating intricate protein aggregations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

MOOSE Framework Meshing Enhancements to Support Reactor Analysis

MOOSE-based physics codes require an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The use of difficult-touse, external licensed software is often required to generate high quality meshes for reactor geometries. High-fidelity geometry modeling also requires elaborate tracking of groups of elements for material property assignment and output reporting which can be considerably complex for the user to identify and maintain. Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, several meshingrelated enhancements have been developed for the MOOSE framework to address user challenges in creating finite element meshes for advanced reactor geometries. MOOSE mesh generators have been developed to mesh hexagonal geometries (pins, ducted assemblies, and cores) commonly found in liquid-metal cooled fast reactor concepts. The mesh generator used for hexagonal pin cells is generic for regular polygons and therefore may also be used for Cartesian pin cells. Hexagonal pin cells can be stitched into ducted assemblies, and assemblies can be stitched together into a core. The user may specify region ids, region names, and other preferences on the mesh. This control is useful for later material mapping in the MOOSE-based physics codes input. A capability was also developed for meshing rotating control drums including determination of material volume fractions in each mesh element as a function of time. Control drum meshes may be stitched to other hexagonal assemblies to create a core configuration. Additional mesh generators were developed that wrap around the hexagonal meshing capabilities and utilize “extra element integer” ID values on each element. In regular Cartesian or hexagonal assemblies or cores, the bookkeeping of element groups for both material assignment and output reporting can now be automated through assignment of pin, assembly, core, axial and depletion id values stored as extra element integers. The extra element tags on the mesh greatly speed the reactor analyst’s efforts to map materials to meshes, track depletion zones, and parse output such as axial pin power distributions. At the highest level, pin, assembly, and core mesh generators (with this reactor terminology) have also been developed to easily generate regular Cartesian and hexagonal cores, including axial extrusion. These reactor geometry builders call upon the previously mentioned capabilities to produce analysis-ready 3D meshes including material assignments. Open source mesh triangulation capabilities were also investigated for integration into the MOOSE framework to address the need for meshing the core periphery region which extends from the irregular outer assembly border to a cylindrical boundary. Options are limited due to licensing constraints, and the recommendation is pursue building a native MOOSE Delaunay triangulator routine with full functionality. Finally, a series of verification problems were performed with NEAMS physics tools. All developed capabilities will be available in the new open-source “Reactor” module of the MOOSE framework, which is accessible to any MOOSE-based NEAMS physics tool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Single crystal yield criterion for chromium based on atomistic studies of isolated 1/2[111] screw dislocations

Using the recently developed Bond Order Potential for Cr, we investigate the mobility of isolated 1/2[111] screw dislocations under stress at 0 K. These studies provide the dependence of the critical resolved shear stress (i.e., the Peierls stress) on the orientation of the maximum resolved shear stress plane and on the shear stress perpendicular to the slip direction. Similar studies are made for loading in tension and compression to show that the onset of yielding in Cr single crystals can be described using a combination of the shear stresses parallel and perpendicular to the slip direction. Based on these results, we develop a single crystal yield criterion for Cr that includes all non-Schmid contributions to the Peierls stress. This criterion is used to predict the shape of the yield surface and compare it with the regular Tresca's hexagon arising from the Schmid law. For all directions of uniaxial load, the criterion is used to quantify the character and magnitude of the tension-compression asymmetry.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Toward controlling wetting hysteresis with nanostructured surfaces derived from block copolymer self-assembly

The synthesis of nanostructured surfaces via block copolymer (BCP) self-assembly enables a precise control of the surface feature shape within a range of dimensions of the order of tens of nanometers. Herein, this work studies how to exploit this ability to control the wetting hysteresis and liquid adhesion forces as the substrate undergoes chemical aging and changes in its intrinsic wettability. Via BCP self-assembly we fabricate nanostructured surfaces on silicon substrates with a hexagonal array of regular conical pillars having a fixed period (52 nm) and two different heights (60 and 200 nm), which results in substantially different lateral and top surface areas of the nanostructure. The wetting hysteresis of the fabricated surfaces is characterized using force–displacement measurements under quasistaic conditions and over sufficiently long periods of time for which the substrate chemistry and surface energy, characterized by the Young contact angle, varies significantly. The experimental results and theoretical analysis indicate that controlling the lateral and top area of the nanostructure not only controls the degree of wetting hysteresis but can also make the advancing and receding contact angles less susceptible to chemical aging. These results can help rationalize the design of nanostructured surfaces for different applications such as self-cleaning, enhanced heat transfer, and drag reduction in micro/nanofluidic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Evaporative Assembly: Dual-Scale Nanostructures via Evaporative Assembly (Adv. Mater. Interfaces 7/2020)

In article number 1901954, Sunita Srivastava, Surita R. Bhatia, and co-workers demonstrate that evaporative assembly of functionalized nanoparticles can yield striking dual-scale hierarchical structures. Regular microscale stripes of nanoparticle monolayers with hexagonal nanoscale order are obtained on physically and chemically homogeneous substrates through evaporation of a suspension of DNAfunctionalized nanoparticles with a charged shell. The stripe width, spacing and nanoparticle ordering can be controlled by varying nanoparticle concentration and can be described by a simple analytical model. The results indicate that the interplay between “stick-slip” motion of the droplet contact line and Coulombic and steric nanoparticle interactions control the formation of the observed structures. Finally, this work demonstrates a simple cost-effective mask-free method for fabricating nanostructured 2D materials and metasurfaces for applications ranging from energy conversion/storage to optoelectronics and nanophotonics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Application of a long short-term memory for deconvoluting conductance contributions at charged ferroelectric domain walls

Ferroelectric domain walls are promising quasi-2D structures that can be leveraged for miniaturization of electronics components and new mechanisms to control electronic signals at the nanoscale. Despite the significant progress in experiment and theory, however, most investigations on ferroelectric domain walls are still on a fundamental level, and reliable characterization of emergent transport phenomena remains a challenging task. Here, we apply a neural-network-based approach to regularize local I ( V )-spectroscopy measurements and improve the information extraction, using data recorded at charged domain walls in hexagonal (Er 0.99 ,Zr 0.01 )MnO 3 as an instructive example. Using a sparse long short-term memory autoencoder, we disentangle competing conductivity signals both spatially and as a function of voltage, facilitating a less biased, unconstrained and more accurate analysis compared to a standard evaluation of conductance maps. The neural-network-based analysis allows us to isolate extrinsic signals that relate to the tip-sample contact and separating them from the intrinsic transport behavior associated with the ferroelectric domain walls in (Er 0.99 ,Zr 0.01 )MnO 3 . Our work expands machine-learning-assisted scanning probe microscopy studies into the realm of local conductance measurements, improving the extraction of physical conduction mechanisms and separation of interfering current signals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Magnetic excitations from the hexagonal spin clusters in the 𝑆 = $\frac{1}{2}$ distorted honeycomb lattice antiferromagnet Cu 2 ⁢(pymca)⁢ 3 ⁢(ClO 4 )

Cu 2 ⁢(pymca) ⁢3 (ClO 4 ) (pymca: pyrimidine-2-carboxylate) consists of a slightly distorted honeycomb lattice of Cu 2+ spins, which shows no long-range magnetic order down to 0.6 K. A magnetization study revealed 1/3 and 2/3 plateau phases [A. Okutani et al., J. Phys. Soc. Jpn. 88, 013703 (2019)], which is not expected for regular honeycomb antiferromagnets. Inelastic neutron scattering experiments were performed using a powder sample to investigate the exchange interactions of this material. The spin excitations from the singlet ground state to the first three triplet states, predicted from the antiferromagnetic hexagonal spin cluster interacting with 3.9 meV, were observed. Using the exact diagonalization methods, the intercluster coupling was estimated from the excitation peak width to be about 20% of the intracluster interaction, which is consistent with the previously reported value. Finally, our exchange path model explains the anisotropic exchange interactions in the distorted honeycomb plane.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Amino Acids Essential for the Assembly of Cellulose Synthase Complexes

Over 40 years ago, Andrew Staehelin's group showed that plant cellulose microfibrils are synthesized by protein complexes arranged in hexagonal arrays called 'particle rosettes'. Plant cellulose synthases (CESAs) differ from their bacterial ancestral forms by insertion of three unique sequences responsible for the assembly into compact multimeric units that aggregate further into rosettes. We used site-directed mutagenesis to replace amino acids within these sequences predicted to be essential for assembly and developed an in vivo method to determine the ability of mutated CesA1 transgenes to complement an Arabidopsis temperature-sensitive root-swelling1 (rsw1) mutant. Replacement of a Cys residue in the Class-Specific Region (CSR) or Pro417 and Arg453 of the Plant-Conserved Region (P-CR) rendered an AtCesA1 transgene unable to complement the rsw1 mutation. Despite an expected role for Arg457 in the trimerization of CESA proteins, AtCesA1 transgenes with Arg457Ala mutations were able to fully restore the wild-type phenotype in rsw1. Staehelin observed a third order of assembly in Micrasterias dendiculata, where hundreds of rosettes form large hexagonal arrays during the synthesis of secondary wall cellulose macrofibrils. We found by SEM of cytosolic faces of Micrasterias cell fragments that these arrays are associated with fields of regularly spaced slime secretion pore complexes.

arabidopsis↗
DOE OSTI2024

Anisotropic physics-regularized interpretable machine learning of microstructure evolution

Anisotropic Physics-Regularized Interpretable Machine Learning Microstructure Evolution (APRIMME) is a general-purpose machine learning solution for grain growth simulations. In prior work, PRIMME employed a deep neural network to predict site-specific migration as a function of its neighboring sites to model normal, isotropic, grain growth behavior. This work aims to extend this method by incorporating grain boundary misorientation-based grain growth behavior. APRIMME is trained on anisotropic simulations created using the Monte Carlo-Potts (MCP) model. Furthermore, the results of this work are compared statistically using grain radius, number of sides per grain, mean neighborhood misorientations, and the standard deviation of triple junction dihedral angles, and are found to match in most cases. The exceptions are small and seem to be related to two causes: (1) the deterministic model of APRIMME is learning from the stochastic simulations of MCP, which seems to accentuate triple junction behaviors; and, (2) a bias against very small grains is made evident in a quicker decrease in grains than expected at the beginning of an APRIMME simulation. APRIMME is also evaluated for its general ability to capture anisotropic grain growth behavior by first investigating different test case initial conditions, including a circle grain, three grain, and hexagonal grain microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dual‐Scale Nanostructures via Evaporative Assembly

Abstract Dual‐scale hierarchical structures with regular microscale patterns and varying degree of nanoscale crystalline order are synthesized on physically and chemically homogeneous substrates by evaporative self‐assembly with a suspension of DNA‐functionalized nanoparticles (NPs) with a charged core shell. For a certain NP concentration range, periodic concentric rings in a stripe‐like micropattern are produced over macroscale surface areas by an NP monolayer with hexagonal lattice structure at the nanoscale. The stripe width, spacing, and nanoparticle ordering can be controlled by varying the NP concentration. The results indicate that the interplay between “stick‐slip” motion of the droplet contact line and coulombic and steric NP interactions control the formation of the observed structures. A simple analytical model is proposed to account for the experimental observations and guide the future design of different nanostructure morphologies. This work demonstrates a simple cost‐effective mask‐free method for fabricating large‐area nanostructured 2D materials and metasurfaces for applications ranging from energy conversion/storage to optoelectronics and nanophotonics.

Srivastava, Sunita↗
DOE OSTI2020

Influence of crystal structure on constitutive anisotropy of silica sand at particle-scale

In spite of recent advanced geotechnical research regarding the influence of particle-scale morphology and fabric on the constitutive behavior of silica sands, there is still no published literature that investigated the influence of crystal structure of sand particles on their constitutive behavior. The main constituent of silica sand particles is quartz mineral, which has a trigonal crystal system and hexagonal lattice. This paper investigates the influence of quartz crystal structure on the constitutive response of silica sand particles using 3D X-ray diffraction (3DXRD), synchrotron micro-computed tomography (SMT), and 3D finite element (FE) analysis. Here, unconfined uniaxial compression experiments on single particles of synthetic silica cubes exposed the notion of constitutive anisotropy caused by the crystal structure of quartz. Since the silica cubes have simple regular geometry (1 mm) 3 , they were loaded in different directions with respect to their crystal local planes. The crystal structure was also identified for a natural silica sand specimen composed of 2705 particles and loaded in confined uniaxial compression. Physical limitations of the experimental measurements made a strong case for the use of 3D FE analysis to study how the change in crystal local orientation of individual sand particles can fundamentally produce anisotropy in the constitutive response of natural sand specimens. The 3D FE analysis was validated to accurately model the crystal-based constitutive anisotropy in silica particles using the results of the silica cube experiments. The analysis was then executed on 3D FE meshes that closely matched the physical shape and fabric of the 2705 particles in the natural sand specimen. Both the experiments and 3D FE analysis showed that the crystal structure of quartz essentially causes a directional anisotropy in the constitutive behavior of silica particles.

58 GEOSCIENCES↗
DOE OSTI2021

Bilayer Wigner crystals in a transition metal dichalcogenide heterostructure

One of the first theoretically predicted manifestations of strong interactions in many-electron systems was the Wigner crystal, in which electrons crystallize into a regular lattice. The crystal can melt via either thermal or quantum fluctuations. Quantum melting of the Wigner crystal is predicted to produce exotic intermediate phases and quantum magnetism because of the intricate interplay of Coulomb interactions and kinetic energy. However, studying two-dimensional Wigner crystals in the quantum regime has often required a strong magnetic field or a moiré superlattice potential, thus limiting access to the full phase diagram of the interacting electron liquid. Here we report the observation of bilayer Wigner crystals without magnetic fields or moiré potentials in an atomically thin transition metal dichalcogenide heterostructure, which consists of two MoSe 2 monolayers separated by hexagonal boron nitride. We observe optical signatures of robust correlated insulating states at symmetric (1:1) and asymmetric (3:1, 4:1 and 7:1) electron doping of the two MoSe 2 layers at cryogenic temperatures. Here, we attribute these features to bilayer Wigner crystals composed of two interlocked commensurate triangular electron lattices, stabilized by inter-layer interaction. The Wigner crystal phases are remarkably stable, and undergo quantum and thermal melting transitions at electron densities of up to 6 × 10 12 per square centimetre and at temperatures of up to about 40 kelvin. Our results demonstrate that an atomically thin heterostructure is a highly tunable platform for realizing many-body electronic states and probing their liquid–solid and magnetic quantum phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Controlled Self-Assembly of Gold Nanotetrahedra into Quasicrystals and Complex Periodic Supracrystals

The self-assembly of shape-anisotropic nanocrystals into large-scale structures is a versatile and scalable approach to creating multifunctional materials. The tetrahedral geometry is ubiquitous in natural and manmade materials, yet regular tetrahedra present a formidable challenge in understanding their self-assembly behavior as they do not tile space. Here, we report diverse supracrystals from gold nanotetrahedra including the quasicrystal (QC) and the dimer packing predicted more than a decade ago and hitherto unknown phases. We solve the complex three-dimensional (3D) structure of the QC by a combination of electron microscopy, tomography, and synchrotron X-ray scattering. Nanotetrahedron vertex sharpness, surface ligands, and assembly conditions work in concert to regulate supracrystal structure. We also discover that the surface curvature of supracrystals can induce structural changes of the QC tiling and eventually, for small supracrystals with high curvature, stabilize a hexagonal approximant. Finally, our findings bridge the gap between computational design and experimental realization of soft matter assemblies and demonstrate the importance of accurate control over nanocrystal attributes and the assembly conditions to realize increasingly complex nanopolyhedron supracrystals.

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

Energetics of hydroxylbastnäsite solid solutions, La 1-$\chi$ Nd $\chi$ CO 3 OH

Bastnäsites (LnCO 3 (F,OH)) are a group of common rare earth elements (REE)-bearing minerals and are one of the primary global sources of REE. Due to the chemical similarities among REE, bastnӓsites tend to occur as solid solutions instead of end members in REE containing ores. To better understand the processes and the mechanisms of formation of such deposits, it is essential to determine the thermodynamic properties of bastnӓsites, including hydroxylbastnӓsite (LnCO 3 OH) solid solutions. In this work, we performed detailed structural and calorimetric investigations on synthetic hexagonal La–Nd hydroxylbastnӓsite (La 1-$\chi$ NdCO 3 OH, $\chi$ = 0, 0.25, 0.5, 0.75, 1) solid solutions. X-ray diffraction confirms the crystal structure of the solid solution series in the $\ P$6 space group, and pair distribution function (PDF) analysis reveals local bonding environments characterized by three different types of 9-coordinated metal-oxygen polyhedra. Unit cell parameters of La1–xNdxCO3OH exhibit a nearly linear relation with the Nd content x, suggesting a random distribution of La and Nd in the structure. Their standard enthalpies of formation (Δ$\ H$° f ) were determined by high temperature oxide melt drop solution calorimetry, from which the enthalpies of mixing (Δ$\ H$ mix ) were derived. The Δ$\ H$ mix can be fitted by a regular solution model with an interaction parameter of 12.58 ± 0.16 kJ/mol, suggesting enthalpic metastability of La1–xNdxCO3OH relative to the two endmembers. Combining entropy and enthalpy, we further estimated the Gibbs free energies of mixing (Δ$\ G$ mix ) at relevant temperatures, revealing favorable temperatures under which the intermediate La 1-$\chi$ NdCO 3 OH phases can be stabilized. Such entropy-driven stabilization, as is consistent with our geochemical modeling results, may explain the enhancement of thermal stability of the solid solutions in nature. Additionally, the temperature range constrained from this study may be used to estimate the thermal history of REE bastnӓsite deposit.

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

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗