Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced basis method”

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 181 records · Page 10

Iterative X-ray spectroscopic ptychography

Spectroscopic ptychography is a powerful technique to determine the chemical composition of a sample with high spatial resolution. In spectro-ptychography, a sample is rastered through a focused X-ray beam with varying photon energy so that a series of phaseless diffraction data are recorded. Each chemical component in the material under investigation has a characteristic absorption and phase contrast as a function of photon energy. Using a dictionary formed by the set of contrast functions of each energy for each chemical component, it is possible to obtain the chemical composition of the material from high-resolution multi-spectral images. This paper presents SPA (spectroscopic ptychography with alternating direction method of multipliers), a novel algorithm to iteratively solve the spectroscopic blind ptychography problem. First, a nonlinear spectro-ptychography model based on Poisson maximum likelihood is designed, and then the proposed method is constructed on the basis of fast iterative splitting operators. SPA can be used to retrieve spectral contrast when considering either a known or an incomplete (partially known) dictionary of reference spectra. By coupling the redundancy across different spectral measurements, the proposed algorithm can achieve higher reconstruction quality when compared with standard state-of-the-art two-step methods. It is demonstrated how SPA can recover accurate chemical maps from Poisson-noised measurements, and its enhanced robustness when reconstructing reduced-redundancy ptychography data using large scanning step sizes is shown.

47 OTHER INSTRUMENTATION↗

GIS-Based Modeling of Contaminated Soil Volumes at Multiple Sites in the Formerly Utilized Sites Remedial Action Program - 20149

The remediation of hazardous, toxic, and radioactive waste (HTRW) sites produces cost-related risks associated with the estimation of contaminated soil or debris volumes. Historical risk-management techniques include cost contingencies to cover volume uncertainties that affect project budgeting and decision-making. The Buffalo District teamed with project partners to lessen volume uncertainty and reduce project risks at multiple HTRW sites managed under the Formerly Utilized Sites Remedial Action Program (FUSRAP). Historical remedial investigations under FUSRAP commonly identified the presence of radiological material in site media, the associated human health risk, and then areas of remediation. To manage remedial execution and reduce risk, pre-design or remediation-phase sampling essentially 'chased' contamination, which was not conducive to efficient predictive budgeting derived from Feasibility Study (FS) cost analyses. The Buffalo District first optimized their approach to better understand volume uncertainty by utilizing the Argonne National Laboratory's Bayesian Approaches for Adaptive Spatial Sampling (BAASS) software [1]. BAASS processed soft data (e.g., gamma walk-over data) and spatial sampling data to estimate the lateral extent of contaminated soil irrespective of depth (i.e., gross contamination extent) and define areas of contaminant uncertainty. The software performed a binary transformation of contaminant concentrations at all sampling points based upon remedial action goals or a sum of ratios approach (i.e., clean, impacted, or range of impacts in soil). The model produced two-dimensional (horizontal) contaminant probability contours and statistical uncertainty in the sampling coverage and resulting contaminant extents. This method was translated vertically by partitioning the sampling data into depth brackets that produced a stacked representation of contaminant extents and uncertainty in the subsurface (i.e., similar to construction lifts). The results commonly led to a better understanding of project uncertainty and the need for sampling strategies that produce high-confidence soil volumes, which control costs. The BAASS-based delineations were eventually replaced by Empirical Bayesian Kriging (EBK) methods available in ArcGIS Spatial or 3D Analysts [2]. The EBK method calculates contaminant probability zones derived from user-controlled semivariograms of the spatial datasets. The resulting probability zones (e.g., 50% or 80% of contaminant probability) represent the two-dimensional surface delineation of the overall horizontal remedial area, similarly to BAASS. However, unlike BAASS, the vertical sampling data within these probability zones became vertical control points to contour a subterranean surface that connects subsurface points to the land-surface delineations of contamination. The resulting representation of horizontal and vertical impacts within an enclosed envelop (volume) of soil included uncertainty distributions that are used to plan uncertainty-reduction sampling. These data-driven and math-based models of three-dimensional sampling results produced well-bounded remedial volumes for project planning and better uncertainty predictions during project budgeting. The EBK method was applied to several FUSRAP sites managed by the Buffalo District and compared to less rigorously modeled sites previously remediated by the District. The comparison of modeled to actual remediated volumes provide a basis for validating the volume-estimation method. This comparison is important to ensure modeled volumes match physical boundaries of site remediation. FUSRAP sites with denser investigative sampling and lesser volume uncertainty proved useful in remedial planning and contracting. The Buffalo District noted that historical sites with sparser sampling arrays had greater disparity between estimated volumes and final remedial volumes. The benefit achieved over the cost of detailed soil sampling appears positive for FUSRAP projects, especially where impacts vary widely and appear unbounded by investigation-phase sampling. The subsequent Empirical Bayesian Kriging of contamination coupled with vertical contouring for soil estimations reduces uncertainty in soil volumes or indicates where sampling is required to reduce uncertainty, which together optimize remedial planning and budgeting. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗

Islet: interpolation semi-Lagrangian element-based transport

Abstract. Advection of trace species, or tracers, also called tracer transport, in models of the atmosphere and other physical domains is an important and potentially computationally expensive part of a model's dynamical core. Semi-Lagrangian (SL) advection methods are efficient because they permit a time step much larger than the advective stability limit for explicit Eulerian methods without requiring the solution of a globally coupled system of equations as implicit Eulerian methods do. Thus, to reduce the computational expense of tracer transport, dynamical cores often use SL methods to advect tracers. The class of interpolation semi-Lagrangian (ISL) methods contains potentially extremely efficient SL methods. We describe a finite-element ISL transport method that we call the interpolation semi-Lagrangian element-based transport (Islet) method, such as for use with atmosphere models discretized using the spectral element method. The Islet method uses three grids that share an element grid: a dynamics grid supporting, for example, the Gauss–Legendre–Lobatto basis of degree three; a physics parameterizations grid with a configurable number of finite-volume subcells per element; and a tracer grid supporting use of Islet bases with particular basis again configurable. This method provides extremely accurate tracer transport and excellent diagnostic values in a number of verification problems.

58 GEOSCIENCES↗

Canister Helium Leak Detection Based on Temperature Information Using RAMM (Final CRADA Report)

The performance of RAMM-TM in detecting gas leakage from a canister with a small, simulated chloride-induced stress corrosion crack has been demonstrated in a series of experiments conducted at CRIEPI by using a 1/4.5- scale model cask. Leakage of both helium and air from a canister was detected within hours after the start of the leakage. The change in surface temperatures at the top and bottom of the canister (ΔTBT) during gas leakage (depressurization) triggered automatic alarms, providing a sound basis for early detection of gas leakage from the canister. This methodology would enable consequence management through the implementation of mitigatory actions to continue effective aging management, to reduce risks to public safety and health, and to protect the environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Correlation Function Approach for Diffusion in Confined Geometries

This paper describes a formalism for extracting spatially varying transport coefficients from simulations of a molecular fluid in a nano channel. This approach is applied to self-diffusion of a Lennard-Jones fluid confined between two parallel surfaces. A numerical grid is laid over the domain confining the fluid, and fluid properties are projected onto the grid cells. The time correlation functions between properties in different grid cells are calculated and can be used as the basis for a fitting procedure for extracting spatially varying diffusion coefficients from the simulation. Results for the Lennard-Jones system show that transport behavior varies sharply near the liquid-solid boundary and that the changes depend on the details of the liquid-solid interaction. A quantitative difference between the reduced and detailed models is discussed. It is found that the difference could be associated with assumptions about the form of the transport equations at molecular scales in lieu of problems with the method itself. The study suggests that this approach to fitting molecular simulations to continuum equations may guide the development of appropriate coarse-grained equations to model transport phenomena at nanometer scales.

Nanoscale flow, Transport, molecular simulation↗

Benchmark study of a new simplified DFN model for shearing of intersecting fractures and faults

It is challenging to quantitatively predict shearing of intersecting fractures/faults because of dynamic frictional contacts accompanied by possible nonlinear rock deformation. To address such challenges, a new conceptual model—the simplified DFN model—was proposed and validated by Hu et al. 46 to use major paths (MPs) to represent complicated DFNs for calculation of shearing. In this work, we conducted a benchmark study for three examples that involve different levels of complexity of intersecting fractures, and correspondingly different numbers of MPs. The codes and software that were used in the benchmark cover a range of continuum, discontinuum and hybrid numerical methods: NMM (LBNL), FLAC3D (LBNL), GBDEM (KIGAM), FRACOD (DynaFrax), and CASRock (CAS). The general consistency between DFN and MP cases as predicted by all the codes/software demonstrates that major paths can be used to simplify the geometry of DFNs in a wide range of software. Disagreement in results made by some software and potential future improvements are discussed. We show that (1) shearing of one or multiple major fractures can be reduced if there are multiple smaller intersecting fractures in that area, which is a useful basis for understanding and controlling induced seismicity and merits further analysis, and (2) the agreement achieved in the benchmark examples provide confidence that the simplified DFN model is a promising conceptual model that can be used for different types of numerical approaches and software for simplifying the analysis of the shearing of intersecting fractures and faults.

58 GEOSCIENCES↗

A Two-Step Control Approach for Torque Ripple and Vibration Reduction in Switched Reluctance Motor Drives

The objective of this work is to reduce the torque ripple and vibration in Switched Reluctance Motor (SRM) drives. To achieve such a goal, a two-step approach is proposed. In the first step, the reference current is profiled using an optimization method, aimed to reduce the torque pulsation. In the second, the optimum current profile is applied to the SRM through an adaptive hysteresis band controller to actively cancel the radial vibration. The proposed technique uses an optimization procedure based on the Field Reconstruction Method (FRM) and Non-Derivative Optimization Method to find the optimal current profile that mitigates the torque pulsation, however, the investigation has shown that the proposed excitation leads to a higher radial vibration, which has been mitigated by using an adaptive hysteresis band controller. The theoretical basis for the proposed approach is initially presented, followed by system modeling and simulation. In addition, experimental verification has been conducted using an 8/6 SRM to show the feasibility of the proposed approach. The proposed method alleviates the acoustic noise and torque pulsation in SRM and offers an effective step towards the prevalent use of SRM in a wide range of applications.

42 ENGINEERING↗

Butterfly, Vinylidene-Like, Monobridged and Trans Structures of Si 2 H 2 + : Comparison to the Well-Characterized Neutral Si 2 H 2

This paper investigates four of the constitutional isomers of Si 2 H 2 + , namely the butterfly, vinylidene-like, monobridged, and trans structures. These isomer geometries were all studied using the CCSD(T) method with basis sets as large as cc-pV5Z. Higher level methods CCSDT and CCSDT(Q) were used for final energetics. It is found that the butterfly isomer has the lowest energy, followed by vinylidene-like and monobridged structures; the trans structure has the highest energy of the four. All structures were compared with their neutral counterparts. Partial charges are computed and it is found that all isomers share the positive charge between the silicon atoms. NBO analysis shows that the cation Si-Si bond order is reduced by 0.5 relative to the neutral versions of the same isomer in the case of the butterfly, vinylidene-like, and monobridged; and reduced by 0.6 for the trans isomer. Vibrational frequencies, infrared intensities, and dipole moments are predicted to encourage the spectroscopic identification of Si 2 H 2 + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reactor and fuel cycle performance of light water reactor fuel with 235 U enrichments above 5%

Recent advances in nuclear fuel materials research, particularly on the topic of accident-tolerant fuels, have brought up potential opportunities for expanding the operating envelope of existing light water reactors. As many of the performance improvements offered by these technologies may be most fully realized by increasing fuel enrichment beyond the standard 5% limit, this paper examines the potential reactor performance and fuel cycle performance of low-enriched uranium oxide fueled light water reactors by generically considering pressurized water reactors with 235 U enrichment from 5 to 7%. Advanced cladding, including accident-tolerant cladding, has the potential to increase fuel burnup limits related to hydrogen in the cladding that coincide with those limits associated with end-of-life reactivity. Therefore, higher enrichment will be necessary in order to realize the higher fuel burnups. This work includes evaluation of the fuel cycle length, discharge burnup, reactivity coefficients, and fuel cycle performance, including radioactive waste and environmental impact metrics per unit energy generated. The analysis was performed using the evaluation metrics from the US Department of Energy Office of Nuclear Energy Fuel Cycle Evaluation and Screening Study. The reactor performance and safety analysis show that enrichments between 5 and 7% would have similar fuel temperature and moderator temperature coefficients. However, the soluble boron coefficient would decrease in magnitude, requiring more corrosive boric acid in the coolant or other methods of reactivity control during the fuel cycle. At these higher enrichments the maximum burnup at the rim of the fuel pellet would increase by almost a factor of two, which is expected to impact the formation of high-burnup structure in the fuel and the corresponding thermo-mechanical fuel properties. The fuel cycle performance assessment shows that increasing enrichment reduces the quantity of high-level waste disposed per unit energy generated, but it increases the natural resource requirements normalized to a gigawatt-electricity-per-year basis. Another impact is the slightly higher discharge burnup, resulting in somewhat different activity levels of the spent nuclear fuel and high-level waste radioactivity at 100 and 100,000 years after fuel discharge. The environmental impacts—including land use, water use, carbon emission, and radiological exposure—are of the same magnitude per unit energy generated. Finally, however, the impacts are distributed differently. Less than 5% enrichment has marginally more impact on the back-end of the fuel cycle, and greater than 5% enrichment has marginally more impact on the front-end of the fuel cycle. Ultimately, no neutronic or reactor safety hindrances to employing light water reactor fuel with enrichments greater than 5% are identified; given the achievable reactor performance benefits with advanced fuels, further practical exploration of increased enrichment fuel is recommended.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Generalized representative structures for atomistic systems

A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.

atomic cluster expansion↗

Electronic structure theory with molecular point group symmetries on quantum annealers

Quantum computation has the potential to revolutionize quantum chemistry through major speedups in computation times and an exponential reduction in computational resources. Here, we combine the symmetry-adapted Jordan–Wigner encoding based on the full Boolean symmetry group $\mathbb{Z}$$^{k}_{2}$ with our new implementation of the Xia–Bian–Kais (XBK) method for improving the efficiency of electronic structure theory calculations on quantum annealers, particularly by reducing the number of qubits needed to achieve the same accuracy. By providing a more extensive symmetry-adapted encoding (SAE) than previous work, we are able to simulate molecules larger than those previously reported that have been studied using methods developed for quantum annealers and without using an active space. We calculated the potential energy surfaces of H 2 , LiH, He 2 , H 2 O, O 2 , N 2 , Li 2 , F 2 , CO, BH 3 , NH 3 , and CH 4 , with the largest molecule in the STO-6G basis set requiring 16 qubits with our SAE, and compared them with full configuration interaction results. The application of SAE to the XBK method provides an exponential reduction in the size of the Hilbert space and scales well with the size of the problem. It does not introduce significant additional errors for even or large values of a key variational parameter that determines the number of ancilla qubits used in the XBK method’s Hamiltonian embedding, or for certain molecules such as He 2 and H 2 O. Here, we provide an explanation for this behavior and a recommendation on the usage of our method. In addition, we briefly discuss the potential of extracting electronic excited states from our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Coupled loads analysis of a novel shared-mooring floating wind farm

Shared mooring lines are a method to reduce the stationkeeping cost of floating offshore wind farms; however, their impacts on floating array dynamics are not yet well understood. Here this article presents advancements to floating wind farm simulation capabilities that allow the dynamic modeling of shared-mooring floating wind arrays - specifically, platform coupling through shared lines and proper phasing of wave loads across the array. These capabilities are demonstrated on a 10-turbine shared-mooring array featuring two staggered rows of wind turbines. A baseline array design with individually moored turbines is used as a basis for comparison. Coupled analysis using FAST.Farm shows that the shared-mooring design has smaller fluctuations in surge when compared to the baseline design. The mean and maximum platform offsets are highest in normal operating conditions, but the maximum mooring line tensions are highest in 50-year storm conditions. Power spectral density analysis of platform surge and sway motions shows that the shared mooring configuration does not introduce any significant inter-platform resonances. The tower-base bending moment statistics and the line tension damage equivalent loads indicate no consequential differences in turbine loading or mooring fatigue life. Anchor load analysis shows that the total required anchor capacity for the shared-mooring array is decreased by approximately 25% when shared anchors are used. Dynamic simulations of the shared mooring system under line failures show that remaining line tensions stay within bounds and that offsets are significantly smaller when compared to the baseline design. Overall, the results show that shared moorings do not introduce any dynamic response concerns in this floating wind array design.

17 WIND ENERGY↗

Active Learning Meets Foundation Models: Fast Remote Sensing Data Annotation for Object Detection

Object detection in remote sensing demands extensive, high-quality annotations—a process that is both labor-intensive and time-consuming. In this work, we introduce a real-time active learning and semi-automated labeling framework that leverages foundation models to streamline dataset annotation for object detection in remote sensing imagery. For example, by integrating a Segment Anything Model (SAM), our approach generates mask-based bounding boxes that serve as the basis for dual sampling: (a) uncertainty estimation to pinpoint challenging samples, and (b) diversity assessment to ensure broad data coverage. Furthermore, our Dynamic Box Switching Module (DBS) addresses the well-known cold start problem for object detection models by replacing its suboptimal initial predictions with SAM-derived masks, thereby enhancing early-stage localization accuracy. Extensive evaluations on multiple remote sensing datasets plus a real-world user study, demonstrate that our framework not only reduces annotation effort, but also significantly boosts detection performance compared to traditional active learning sampling methods. The code for training and the user interface will be made available.

Burges, Marvin [ORNL] (ORCID:0000000312690769)↗

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↗

Optimized Auxiliary Functions for Robust Mitigation of Finite-Size Errors in Periodic Hybrid Density Functional Theory

When calculating properties of periodic systems at the thermodynamic limit (TDL), the dominant source of finite size error (FSE) arises from the long-range Coulomb interaction, and can manifest as a slowly converging quadrature error when approximating an integral in the reciprocal space by a finite sum. The singularity subtraction (SS) method offers a systematic approach for reducing this quadrature error and thus the FSE. Here, in this work, we first investigate the performance of the SS method in the simplest setting, aiming at reducing the FSE in exact exchange calculations by subtracting the Coulomb contribution with a single, adjustable Gaussian auxiliary function. We demonstrate that a simple fitting method can robustly estimate the optimal Gaussian width and leads to rapid convergence toward the TDL. Furthermore, we suggest new forms of the auxiliary function, whose optimal parameters could also be determined through least-squares fitting. For a range of semiconductors and insulators, the proposed auxiliary functions achieve robust, millihartree-level accuracy in hybrid density functional theory calculations, including cases with sparse k-meshes and large basis sets.

Quiton, Stephen Jon [University of California, Ber↗