Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Packaging Demonstration”

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 163 records · Page 9

The fastVFP code for solution of the Vlasov–Fokker–Planck equation

Abstract We describe the fast VFP code for solution of the Vlasov–Fokker–Planck equation for non-local electron transport and the generation of magnetic field, especially for application to laser-produced plasmas. We describe the essential features of the code that make it fast and robust and suitable for inclusion as a transport package in a fluid simulation. We present a few sample results that demonstrate the abilities of the code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lipid Droplets: Packing Hydrophobic Molecules Within the Aqueous Cytoplasm

Lipid droplets, also known as oil bodies or lipid bodies, are plant organelles that compartmentalize neutral lipids as a hydrophobic matrix covered by proteins embedded in a phospholipid monolayer. Some of these proteins have been known for decades, such as oleosins, caleosins, and steroleosins, whereas a host of others have been discovered more recently with various levels of abundance on lipid droplets, depending on the tissue and developmental stage. In addition to a growing inventory of lipid droplet proteins, the subcellular machinery that contributes to the biogenesis and degradation of lipid droplets is being identified and attention is turning to more mechanistic questions regarding lipid droplet dynamics. While lipid droplets are mostly regarded as storage deposits for carbon and energy in lipid-rich plant tissues such as seeds, these organelles are present in essentially all plant cells, where they display additional functions in signaling, membrane remodeling, and the compartmentalization of a variety of hydrophobic components. Remarkable metabolic engineering efforts have demonstrated the plasticity of vegetative tissues such as leaves to synthesize and package large amounts of storage lipids, which enable future applications in bioenergy and the engineering of high-value lipophilic compounds. Here, we review the growing body of knowledge about lipid droplets in plant cells, describe the evolutionary similarity and divergence in their associated subcellular machinery, and point to gaps that deserve future attention.

60 APPLIED LIFE SCIENCES↗

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Expanded Energy Credits in Energy Codes (Technical Brief)

This technical brief provides additional energy efficiency measures that go beyond the current prescriptive commercial energy codes. It demonstrates relative savings for multiple measures and shows a base cost effective package by building type and climate zone. An advanced savings package is also shown. Potential building savings for the base package for different building types across climate zones range from 3.9% to 22% of building energy cost, while the advanced package building savings range from 8.4% to 44%. In addition to measure descriptions, this technical brief includes code language that can be adopted by local jurisdictions to implement these measures for new buildings.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

IACMI Project 4.7: Pultruded Textile Carbon Fiber for Spar Caps (Final Report)

The primary objective of this project was to demonstrate the potential to significantly reduce the cost of wind turbine blades with carbon fiber reinforced polymer (CFRP) structure. Applicability of textile carbon fibers (TCF) were evaluated for use in pultruded spar cap (SC) elements as a path to cost reduction for utility scale wind turbine blades. In earlier work for the Department of Energy (DOE) Wind Energy Technologies Office (WETO), a collaboration of Sandia National Laboratory (SNL), Oak Ridge National Laboratory (ORNL), and Montana State University has demonstrated potential for pultruded TCF to compete with infused fiberglass and commercially available carbon fiber pultruded sections for spar cap construction. In the design cases evaluated, the TCF sections fared well when compared on cost per unit composite stiffness and cost per unit composite compressive strength for those designs [1]. Both stiffness and compressive strength tend to be key factors in the design of blade composite Spar Cap which carry the bulk of the blade structural loads in bending. Spar Cap design tends to distribute largely symmetric tensile and compressive stresses to opposite sides of the spar structure, but since carbon fiber composite compressive strength is typically 20-50% lower than tensile strength, the compressive loading reaches failure levels well before the tensile loading. Stiffness is critical in containing the large tip deflection in high wind loading situations. However, materials and process development were very limited in the earlier study and the work in this project was expanded to make the comparative information more representative of what will be required in order to make further inroads towards implementation. Similar to that study, this project team confirmed that the primary materials of interest for pultruded spar cap elements should be thermoset (TS) resins reinforced by carbon fibers, utilizing as high a percentage of TCF as practical to benchmark cost and performance against commercial carbon fibers. To make the closest comparison possible and eliminate specific test article size, resin selection, and equipment/operational nuance effects, the team planned to pultrude sections with 100% commercially available carbon fiber (Panex 35 carbon fiber from Zoltek) as well as samples utilizing high fractions of TCF. The resin system chosen was based on formulations recommended by large wind industry supplier Hexion and consisted of Hexion resin RSL-4597, curing agent CCA-138, and internal mold release additive 117, along with common kaolin filler ASP400P from BASF. As commonly deployed in spar cap configurations, the team had a mold built to pultrude a rectangular spar cap element of 100mm width and 3mm thickness. The extremely limited number of samples produced for the earlier study were produced with a “generic” epoxy utilized for a variety of applications by the pultruder contracted to produce test articles for demonstration purposes. More importantly, those samples were produced at a fiber fraction only slightly over 50%. Based on feedback from our industrial advisory team for that project and strongly recommended by this project team, the consensus is that it is highly desirable to obtain fiber fractions of 65-68% for significant penetration in wind blade spar cap. Although this requirement has yet to be exhaustively confirmed in readily available information, this was established as a project goal and informally decided we needed to exceed 60% fiber fraction to gain serious industry consideration. Previous TCF pultrusion trials have been challenged by the lack of robust TCF packages, resulting in non-uniform tension across and between tows, as well as excess labor and waste for removal of interleaved paper. The non-uniform tension and associated intermingling of tows in textile acrylic fiber tows and associated difficulties created from broken filaments in carbon fiber conversion inhibit the ordered packing necessary to enhance fiber fraction elevation. (These “cross-overs” are not considered undesirable for textile applications and there is some sense that they might be advantageous for those applications). In addition to work that is ongoing at the acrylic fiber manufacturers to improve their formats, The Institute of Advanced Composites Manufacturing Innovation (IACMI) Project 6.12 (report PA16-0349-6.12-01) [2] has developed and demonstrated a more robust packaging and creeling approach that at least partially addresses these issues, thus improving control of the TCF feed into the pultrusion unit. It was hoped that these and other improvements currently being implemented would allow us to achieve fiber fractions at least approaching these fiber fraction targets. During this project, sections utilizing 100% commercially available carbon fiber reinforcement were produced as a baseline, as well as sections reinforced with about 94% TCF and the balance being commercially available fiber for comparison. The most important finding was that similar to results reported in the earlier WETO-funded project and results from tests of TCF reported at IACMI meetings, this work demonstrated that sections pultruded with TCF in an epoxy resin frequently utilized in actual spar cap production had stiffness and compressive strengths largely comparable to similar sections pultruded with a commercially available carbon fiber also frequently utilized in the wind industry. Although the amount of that data is limited, some of the tensile strength results were actually closer than would have been expected based on fiber strength results provided by the TCF and commercial fiber producers. The actual test data are reported and discussed in detail in Section 5. The pultruded sections dominated by TCF reinforcement were approximately 8-10% lower in fiber fraction than for the sections produced using commercial fiber alone, making direct comparison difficult. The COVID-19 project has provided significant insight into the current state-of-the-art with various TCF product forms. The data obtained in this project will guide the planned improvements at the precursor level, especially in attaining uniform tensioning and payout to facilitate enhanced fiber fractions and overall processability of the TCF composites. The project team is providing guidance to stakeholders concerning the attributes, needs, and potential demand for TCF in wind blade spar caps. Results achieved in this project are consistent with findings in the related work cited [1] and support this guidance and the high potential for this product type. TCF precursor-producing partners continue to express interest in enhancing their product forms and the team looks forward to working with these improved materials as they become available.

42 ENGINEERING↗

RH-TRU Waste Streams: HFDA Crucibles Create New Complex Waste Stream

This project focuses on the removal of the Hot Fuel Dissolution Apparatus (HFDA) crucibles from the Hot Fuel Examination Facility’s (HFEF) hot cells to create space for upcoming projects. The project team has created a crucible and pyro-salt non-debris waste stream that can be disposed of at Waste Isolation Pilot Plant (WIPP). Chemical Compatibility Examinations (CCEs) must be performed on candidate waste by reviewing the types of chemicals used in the experiments, information found in INL Process Knowledge Summary Reports (PKSRs), and Environmental Protection Agency guides on chemical compatibility. The CCEs identify which chemicals can be safely combined to ensure adverse reactions do not occur in the waste packages. The CCEs are then used to update the PKSR, Basis of Knowledge, and Acceptable Knowledge documents and demonstrate compliance with the WIPP Waste Acceptance Criteria. Once documents are finalized, they are reviewed and approved by the WIPP Central Characterization Project (CCP) to ensure proper waste stream characterization. After the documents are approved, the crucibles can be packaged and shipped to the Idaho Nuclear Technical and Engineering Center for storage awaiting certification. CCP will certify the waste using the approved documents and a series of confirmatory testing before being shipped and disposed in WIPP’s deep geological repository. Currently, the CCE is in progress and the PKSRs are being updated.

Basis of Knowledge↗

Demonstrating solarpilot ’s Python Application Programmable Interface Through Heliostat Optimal Aimpoint Strategy Use Case

solarpilot is a software package that generates solar field layouts and characterizes the optical performance of concentrating solar power (CSP) tower systems. solarpilot was developed by the National Renewable Energy Laboratory (NREL) as a stand-alone desktop application but has also been incorporated into NREL’s System Advisor Model (SAM) in a simplified format. Prior means for user interaction with solarpilot have included the application’s graphical interface, the SAM routines with limited configurability, and through a built-in scripting language called “LK.” This article presents a new, full-featured, python-based application programmable interface (API) for solarpilot, which we hereafter refer to as CoPylot. CoPylot enables python users to perform detailed CSP tower analysis utilizing either the Hermite expansion technique (analytical) or the SolTrace ray-tracing engine. CoPylot’s enables CSP researchers to perform analysis that was previously not possible through solarpilot’s existing interfaces. This article discusses the capabilities of CoPylot and presents a use case wherein we populate a model that obtains optimal solar field aiming strategies.

14 SOLAR ENERGY↗

Hydrogenolysis of Poly(Ethylene–co–Vinyl Alcohol) and Related Polymer Blends over Ruthenium Heterogeneous Catalysts

The hydrogenolysis of polymers is emerging as a promising approach to deconstruct plastic waste into valuable chemicals. Yet, the complexity of plastic waste, including multilayer packaging, is a significant barrier to handling realistic waste streams. Herein, we reveal fundamental insights into a new chemical route for transforming a previously unaddressed fraction of plastic waste – poly(ethylene–co­­–vinyl alcohol) (EVOH) and related polymer blends – into alkane products. Additionally, we report that Ru/ZrO2 is active for the concurrent hydrogenolysis, hydrogenation, and hydrodeoxygenation of EVOH and its thermal degradation products into alkanes (C1–C35) and water. Detailed reaction data, product analysis, and catalyst characterization reveal that the in–situ thermal degradation of EVOH forms aromatic intermediates that are detrimental to catalytic activity. Increased hydrogen pressure promotes hydrogenation of these aromatics, preventing catalyst deactivation and improving alkane product yields. Calculated apparent rates of C–C scission reveal that the hydrogenolysis of EVOH is slower than low–density polyethylene. We apply these findings to achieve hydrogenolysis of EVOH/polyethylene blends and elucidate the sensitivity of hydrogenolysis catalysts to such blends. Overall, we demonstrate progress towards efficient catalytic processes for the hydroconversion of waste multilayer film plastic packaging into valuable products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

All-electron APW+ lo calculation of magnetic molecules with the SIRIUS domain-specific package

We report APW+lo (augmented plane wave plus local orbital) density functional theory (DFT) calculations of large molecular systems using the domain specific SIRIUS multi-functional DFT package. The APW and FLAPW (full potential linearized APW) task and data parallelism options and the advanced eigen-system solver provided by SIRIUS can be exploited for performance gains in ground state Kohn–Sham calculations on large systems. This approach is distinct from our prior use of SIRIUS as a library backend to another APW+lo or FLAPW code. We benchmark the code and demonstrate performance on several magnetic molecule and metal organic framework systems. We show that the SIRIUS package in itself is capable of handling systems as large as a several hundred atoms in the unit cell without having to make technical choices that result in the loss of accuracy with respect to that needed for the study of magnetic systems.

Chemistry↗

Measuring 3D Profilometry of SAVY-4000 Nuclear Material Storage Containers: Novacam TubeInspect Capabilities Report

The SAVY-4000 container series is a general-purpose interim storage container for nuclear materials, developed and maintained by Los Alamos National Laboratory (LANL). It is the first vented, general-use nuclear material container to be demonstrated as meeting the requirements outlined in DOE M 441.1-1, the Nuclear Material Packaging Manual. Due to the challenging radiation, thermal, and corrosive storage conditions that the SAVY containers must endure, continuous surveillance techniques are employed to ensure the containers meet all safety standards and specifications. These inspections are typically performed by human operators, who check for issues such as corrosion, O-ring deterioration, corrosion, filter integrity, and potential manufacturing defects. However, human inspections alone are not sufficient, and automated inspection technologies, such as the ATIS system, as well as other automated systems are also utilized. The MicroCam TubeInspect, developed by Novacam Technologies Inc., is designed to address the challenges of understanding how manufacturing variations in the SAVY-4000 container series may affect performance. It is a 3D profilometry measurement system that enables detailed analysis of surface features, including defects, surface roughness, and manufacturing variations. This advanced tool significantly enhances rapid surveillance techniques for both pristine and used containers. In this study, container properties such as surface roughness, thickness, and geometric attributes like circularity are measured for SAVY-4000 containers. Artificially corroded or dented containers are examined to demonstrate the MicroCam's ability to quantify defects. A sensitivity analysis is also conducted, comparing the MicroCam results to those obtained using more precise instruments such as confocal microscopy. This comparison aims to provide valuable insights into container quality, durability, and potential improvements in manufacturing processes.

42 ENGINEERING↗

Sparse Matrix-Based HPC Tomography

Tomographic imaging has benefited from advances in X-ray sources, detectors and optics to enable novel observations in science, engineering and medicine. These advances have come with a dramatic increase of input data in the form of faster frame rates, larger fields of view or higher resolution, so high performance solutions are currently widely used for analysis. Tomographic instruments can vary significantly from one to another, including the hardware employed for reconstruction: from single CPU workstations to large scale hybrid CPU/GPU supercomputers. Furthermore, flexibility on the software interfaces and reconstruction engines are also highly valued to allow for easy development and prototyping. This paper presents a novel software framework for tomographic analysis that tackles all aforementioned requirements. The proposed solution capitalizes on the increased performance of sparse matrix-vector multiplication and exploits multi-CPU and GPU reconstruction over MPI. Furthermore, the solution is implemented in Python and relies on CuPy for fast GPU operators and CUDA kernel integration, and on SciPy for CPU sparse matrix computation. As opposed to previous tomography solutions that are tailor-made for specific use cases or hardware, the proposed software is designed to provide flexible, portable and high-performance operators that can be used for continuous integration at different production environments, but also for prototyping new experimental settings or for algorithmic development. The experimental results demonstrate how our implementation can even outperform state-of-the-art software packages used at advanced X-ray sources worldwide.

97 MATHEMATICS AND COMPUTING↗

Low Noise Frequency-Domain Multiplexing of TES Bolometers Using SQUIDs at Sub-Kelvin Temperature

Digital Frequency-Domain Multiplexing (DfMux) is a technique that uses MHz superconducting resonators and Superconducting Quantum Interference Device (SQUID) arrays to read out sets of transition edge sensors. DfMux has been used by several Cosmic Microwave Background experiments, including most recently POLARBEAR-2 and SPT-3 G with multiplexing factors as high as 68, and is the baseline readout technology for the planned satellite mission LiteBIRD. Here, we present recent work focused on improving DfMux readout noise, reducing parasitic impedance, and improving sensor operation. We have achieved a substantial reduction in stray impedance by integrating the sensors, resonators, and SQUID array onto a single-carrier board operated at 250 mK. This also drastically simplifies the packaging of the cryogenic components and leads to better-controlled crosstalk. We demonstrate a low readout noise level of 8.6 pA/Hz1 1/2 , which was made possible by operating the SQUID array at a reduced temperature and with a low dynamic impedance. This is a factor of two improvement compared to the achieved readout noise level in currently operating Cosmic Microwave Background experiments using DfMux and represents a critical step toward maturation of the technology for the next generation of instruments.

47 OTHER INSTRUMENTATION↗

Massively scalable workflows for quantum chemistry: BigChem and ChemCloud

Electronic structure theory, i.e., quantum chemistry, is the fundamental building block for many problems in computational chemistry. Here we present a new distributed computing framework (BigChem), which allows for an efficient solution of many quantum chemistry problems in parallel. BigChem is designed to be easily composable and leverages industry-standard middleware (e.g., Celery, RabbitMQ, and Redis) for distributed approaches to large scale problems. BigChem can harness any collection of worker nodes, including ones on cloud providers (such as AWS or Azure), local clusters, or supercomputer centers (and any mixture of these). BigChem builds upon MolSSI packages, such as QCEngine to standardize the operation of numerous computational chemistry programs, demonstrated here with Psi4, xtb, geomeTRIC, and TeraChem. BigChem delivers full utilization of compute resources at scale, offers a programable canvas for designing sophisticated quantum chemistry workflows, and is fault tolerant to node failures and network disruptions. We demonstrate linear scalability of BigChem running computational chemistry workloads on up to 125 GPUs. Finally, we present ChemCloud, a web API to BigChem and successor to TeraChem Cloud. ChemCloud delivers scalable and secure access to BigChem over the Internet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗

Adversarial Sampling-Based Motion Planning

In this report there are many scenarios in which a mobile agent may not want its path to be predictable. Examples include preserving privacy or confusing an adversary. However, this desire for deception can conflict with the need for a low path cost. Optimal plans such as those produced by RRT* may have low path cost, but their optimality makes them predictable. Similarly, a deceptive path that features numerous zig-zags may take too long to reach the goal. We address this trade-off by drawing inspiration from adversarial machine learning. We propose a new planning algorithm, which we title Adversarial RRT*. Adversarial RRT* attempts to deceive machine learning classifiers by incorporating a predicted measure of deception into the planner cost function. Adversarial RRT* considers both path cost and a measure of predicted deceptiveness in order to produce a trajectory with low path cost that still has deceptive properties. We demonstrate the performance of Adversarial RRT*, with two measures of deception, using a simulated Dubins vehicle. We show how Adversarial RRT* can decrease cumulative RNN accuracy across paths to 10%, compared to 46% cumulative accuracy on near-optimal RRT* paths, while keeping path length within 16% of optimal. We also present an example demonstration where the Adversarial RRT* planner attempts to safely deliver a high value package while an adversary observes the path and tries to intercept the package.

42 ENGINEERING↗

Summary of the Brine Availability Test in Salt (BATS), Including Extended Plan for Experiments at the Waste Isolation Pilot Plant (WIPP)

This document lays out a set of near-future investigations in salt, the third phase of BATS (BATS 3). This phase is planned to answer the few remaining issues from the first two phases of BATS (BATS 1 and BATS 2), and to prepare for a subsequent large-scale demonstration phase. The BATS experiments are the first part of a larger plan to conduct field experiments to answer specific technical questions, improve the technical basis for disposal of heat-generating radioactive waste in salt (Stauffer et al., 2015; SNL et al., 2020), and demonstrate readiness for disposal of radioactive waste in salt, including large, hot waste packages.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Riemannian Optimization Applied to AC Optimal Power Flow: Preprint

The nonlinear, nonconvex AC optimal power flow problem is of growing importance as the nature of the power grid evolves. This problem can be difficult to solve for interior point methods. However, the advent of optimization algorithms over smooth Riemannian manifolds presents an alternative approach. The nonlinear, nonconvex constraints in the AC power flow problem form an embedded submanifold of Euclidean space. In this paper, the authors explore the performance of Riemannian optimization algorithms for the ACOPF problem where the optimization is performed directly on the AC power flow manifold. They demonstrate that these are viable computational alternatives to interior point methods. This is done by using Julia and the packages PowerModels.jl and Manopt.jl.

manifold optimization↗

GBOpt: Grain boundary structure optimization using Monte Carlo and evolutionary algorithms

Polycrystalline materials are made of many small crystals separated by grain boundaries (GBs), whose atomic structure strongly influences material properties. Because the structure of a GB determines its properties, the optimal structure must be known in order to determine those impacts. There are many ways of placing atoms in the GB region, but the optimal structure is defined as the one that gives the lowest value of a target property (typically energy). GB structure optimization has been successfully demonstrated using stochastic and evolutionary methods, but no reusable, community-maintained open-source workflow has been developed. GBOpt (Grain Boundary Optimization) is an open-source Python package that creates that workflow, where we have presently implemented two approaches: Markov Chain Monte Carlo, and genetic algorithm based on elite selection. We demonstrate this capability by successfully reproducing the known optimal structures of a specific GB in two materials, and point interested readers to the GitHub repository for additional examples, including optimization for different properties. Both of the implemented approaches recovered the known structures, with the genetic algorithm approach finding the optimal structure faster on average.

99 - GENERAL AND MISCELLANEOUS↗