Feedforward control of thermal history in laser powder bed fusion: Toward physics-based optimization of processing parameters
Not Available
SEARCH · Search NASA
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.
Not Available
Here, this study investigates the fabrication of oxide dispersion strengthened (ODS) Fe-Cr-Al alloys via laser powder bed fusion (LPBF) with strategic additions of 0.3 wt% Zr and 0.3 wt% Y 2 O 3 for enhanced mechanical performance in nuclear applications. Systematic processing parameter optimization yielded three distinct conditions: one low-density product with significant defects and two near-full-density materials with improved consolidation. Comprehensive characterization confirmed single-phase α-ferrite matrix formation with successful incorporation of Y-, Zr-, O-, and C-rich precipitates characteristic of ODS alloys. However, precipitate density remained low (∼10 7 cm −3 ), resulting in sink strength values substantially below optimal levels for radiation resistance. Microhardness values (mid-200s HV) correlated inversely with grain size following the Hall-Petch relationship, indicating grain boundary strengthening as the dominant mechanism rather than precipitation strengthening. The optimized processing conditions achieved excellent mechanical properties with room temperature yield strength of approximately 500 MPa and 30 % elongation, demonstrating superior strength-ductility synergy compared to other additively manufactured ODS materials and performance consistent with literature values for LPBF-processed ODS-FeCrAl alloys. This investigation reveals both the potential and limitations of LPBF processing for ODS Fe-Cr-Al alloys. While successful defect-free fabrication was achieved, results highlight the critical need for systematic optimization of processing parameters and post-processing heat treatments to enhance precipitate density for effective dispersion strengthening and radiation resistance while maintaining additive manufacturing advantages.
Static structural finite element models of an aluminum-framed crystalline silicon (c-Si) photovoltaic (PV) module and a glass-glass thin film PV module were constructed and validated against experimental measurements of deflection under uniform pressure loading. Parametric analyses using Latin Hypercube Sampling (LHS) were performed to propagate simulation input uncertainties related to module material properties, dimensions, and manufacturing tolerances into expected uncertainties in simulated deflection predictions. This exercise verifies the applicability and validity of finite element modeling for predicting mechanical behavior of solar modules across architectures and enables computational models to be used with greater confidence in assessment of module mechanical stressors and design for reliability. Sensitivity analyses were also performed on the uncertainty quantification data sets using linear correlation coefficients to elucidate the key parameters influencing module deformation. This information has implications on which materials or parameters may be optimized to best increase module stiffness and reliability, whether the key optimization parameters change with module architecture or loading magnitudes, and whether parameters such as frame design and racking must be replicated in reduced-scale reliability studies to adequately capture full module mechanical behavior.
While abstract models of quantum computation assume a closed system of two-level states, practical quantum devices inevitably couple to the environment in some way, creating sources of noise. Understanding the tolerance to noise of specific quantum algorithms run on specific devices is important for determining the feasibility of quantum computing in the current noisy intermediate-scale quantum era. Of particular interest is understanding the noise sensitivity of these devices as more qubits are added to the system. Classical simulations are a useful tool to understand the effects of this noise, but direct classical simulations of open quantum systems are burdened by an exponentially growing cost in the number of qubits and a large local Hilbert space dimension. For onedimensional, shallow circuits, using tensor networks can replace this exponential cost with a linear one and simulate far wider systems than what would normally be available. In this paper, we describe a tensor network simulation of a neutral atom quantum system under the presence of noise, while introducing a purity-preserving truncation technique that compromises between the simplicity of the matrix product state and the positivity of the matrix product density operator. We apply this simulation to a near-optimized iteration of the quantum approximate optimization algorithm on a transverse field Ising model in order to investigate the influence of large system sizes on the performance of the algorithm. We find that while circuits with a large number of qubits fail more often under noise that depletes the qubit population, their outputs on a successful measurement are just as robust under Rydberg atom dissipation or qubit dephasing as smaller systems. However, such circuits might not perform as well under coherent multiqubit errors such as Rydberg atom crosstalk. We also find that the optimized parameters are especially robust to noise, suggesting that a noisier quantum system can be used to find the optimal parameters before switching to a cleaner system for measurements of observables.
ABSTRACT Agrobacteria are a diverse, polyphyletic group of prokaryotes with multipartite genomes capable of transferring DNA into the genomes of host plants, making them an essential tool in plant biotechnology. Despite their utility in plant transformation, genome-wide transcriptional regulation is not well understood across the three main lineages of agrobacteria. Transcription start sites (TSSs) are a necessary component of gene expression and regulation. In this study, we used differential RNA-seq and a TSS identification algorithm optimized on manually annotated TSS, then validated with existing TSS to identify thousands of TSS with nucleotide resolution for representatives of each lineage. We extend upon the 356 TSSs previously reported in Agrobacterium fabrum C58 by identifying 1,916 TSSs. In addition, we completed genomes and phenotyping of Rhizobium rhizogenes C16/80 and Allorhizobium vitis T60/94, identifying 2,650 and 2,432 TSSs, respectively. Parameter optimization was crucial for an accurate, high-resolution view of genome and transcriptional dynamics, highlighting the importance of algorithm optimization in genome-wide TSS identification and genomics at large. The optimized algorithm reduced the number of TSSs identified internal and antisense to the coding sequence on average by 90.5% and 91.9%, respectively. Comparison of TSS conservation between orthologs of the three lineages revealed differences in cell cycle regulation of ctrA as well as divergence of transcriptional regulation of chemotaxis-related genes when grown in conditions that simulate the plant environment. These results provide a framework to elucidate the mechanistic basis and evolution of pathology across the three main lineages of agrobacteria. IMPORTANCE Transcription start sites (TSSs) are fundamental for understanding gene expression and regulation. Agrobacteria, a group of prokaryotes with the ability to transfer DNA into the genomes of host plants, are widely used in plant biotechnology. However, the genome-wide transcriptional regulation of agrobacteria is not well understood, especially in less-studied lineages. Differential RNA-seq and an optimized algorithm enabled identification of thousands of TSSs with nucleotide resolution for representatives of each lineage. The results of this study provide a framework for elucidating the mechanistic basis and evolution of pathology across the three main lineages of agrobacteria. The optimized algorithm also highlights the importance of parameter optimization in genome-wide TSS identification and genomics at large.
Laser Direct Energy Deposition (DED) is a metal additive manufacturing technique with the ability to fabricate large and complex parts through deposition of metal powders. However, achieving high-density parts and targeted build heights using DED can be challenging due to the large number of highly sensitive process variables. This work proposes a robust fabrication parameter optimization framework to generate process maps for primary parameters in DED, including laser power, scan speed, mass flow rate, hatch spacing, and layer height. Simple single-track experiments were utilized to map out the parameter space, and a combination of geometric criteria for hatch spacing and layer height were proposed to determine parameter sets that achieve both targeted build heights and mitigate porosity formation. Using this framework, specimens with >99 % density and consistent mechanical properties were successfully fabricated over a wide range of process parameters for an Fe-9wt.%Cr (Fe9Cr) alloy, a surrogate for radiation damage-resistant reduced activation ferritic/martensitic (RAFM) steels. Processing these materials using DED is of particular interest in the development of plasma facing components for nuclear fusion applications. The microstructure and mechanical properties of as-printed Fe9Cr were characterized using optical and electron microscopy, X-ray diffraction, and uniaxial tensile tests. As-printed Fe9Cr displayed ~25 % elongation and ultimate tensile strengths of up to 475 MPa which is comparable to similar wrought alloys. Finally, the proposed framework will allow for accelerated DED parameter optimization for novel alloy systems, as well as open the possibility for local microstructure control while simultaneously mitigating defect formation.
We develop a gradient based optimization approach for the equations of compressible, Lagrangian hydrodynamics and demonstrate how it can be employed to automatically uncover strategies to control hydrodynamic instabilities arising from shock acceleration of density interfaces. Strategies for controlling the Richtmyer-Meshkov instability (RMI) are of great benefit for inertial confinement fusion (ICF) where shock interactions with many small imperfections in the density interface lead to instabilities which rapidly grow over time. These instabilities lead to mixing which, in the case of laser driven ICF, quenches the runaway fusion process ruining the potential for positive energy return. Here, we demonstrate that control of these instabilities can be achieved by optimization of initial conditions with ( > 100) parameters. Optimizing over a large parameter space like this is not possible with gradient-free optimization strategies. This requires computation of the gradient of the outputs of a numerical solution to the equations of Lagrangian hydrodynamics with respect to the inputs. We show that the efficient computation of these gradients is made possible via a judicious application of (i) adjoint methods, the exact formal representation of sensitivities involving partial differential equations, and (ii) automatic differentiation (AD), the algorithmic calculation of derivatives of functions. Careful regularization of multiple operators including artificial viscosity and timestep control is required. We perform design optimization of > 100 parameter energy field driving the Richtmyer Meshkov instability showing significant suppression while simultaneously enhancing the acceleration of the interface relative to a nominal baseline case.
This repository contains both simulated and experimental structure factor data for the data challenge involving the inverse scattering problem. The simulated data were generated during a machine-learning-assisted optimization routine described in ref[1]. The experimental structure factor was measured on a rare-earth oxide, Dy2Ti2O7 using diffuse neutron scattering from time-of-flight techniques on the CORELLI instrument at the Spallation Neutron Source, Oak Ridge National Laboratory. A Metropolis Monte Carlo code implemented to run in a High-performance computing setting was used to calculated simulated structure factors for the spin-ice Hamiltonian at 680 mK, which is the same temperature as for the experimental data. The total size of all the files in this repository is 5.12 GB. A detailed description of the files is given below. ExperimentalData_630mK.dat – A linearized version of 3-dimensional experimental data of size 61×81×21. This data was processed to remove an estimation of non-magnetic background, including nuclear scattering signal and instrumentation background. Parameters.dat – 6700 samples were evaluated over the 4-dimensional parameter space (J_1, J_2, J_3 and J_(3^' )). There is an additional parameter, D in the spin Hamiltonian to mimic the dipolar interaction between magnetic ions. However, this parameter, D was fixed to a value determined by prior work. This file contains five columns for the parameters J_1, J_2, J_3, J_(3^' ) and D respectively. 3D_Simulation_Data.dat – The simulated structure factor, S(Q) data are included in this file. Each raw contains a linearized array of 3D volumes of S(Q) calculated for the parameter set given in the corresponding row of the file Parameters.dat. The size of the volume data was matched to the experimental data. Qx(h,-h,0).dat, Qy(k,k,-2k).dat, Qz(l,l,l).dat – These files contain the h, k, and l values along with the reciprocal vectors [h,-h,0], [k,k,-2k] and [l,l,l] respectively.
Chemical state analysis of a sample surface through fitting bell-shaped curves to x-ray photoelectron spectroscopic polymer data is reviewed using nylon to introduce and discuss aspects of data analysis. Different strategies for modeling chemistry in nylon spectra are presented and in so doing, a case is made to include in published science the design logic and implementation in terms of line shapes and optimization parameter constraints between components in a peak model. Imperfections in line shape relative to the true shape for photoemission lines, when compensated for using constraints to optimization parameters, are shown to provide chemical state information about a sample that justify, for peak models constructed with these limitations, metrics for goodness-of-fit different from those expected for pulse-counted data.
Thermomechanical analysis (TMA) experiments conducted on high-density polyethylene (HDPE) show both reversible and irreversible dimensional changes. To further explore these reversible and irreversible processes, modulated thermomechanical analysis (MTMA) was used. Before reliable data on compression-molded HDPE was collected, a parameter optimization was performed to obtain a suitable MTMA method. Once a suitable method was obtained, several MTMA experiments were conducted on compression-molded HDPE. This work highlights the steps taken during the MTMA parameter optimization and the results obtained from MTMA experiments conducted on pristine compression-molded HDPE samples.
CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.
Work proposed in this project focused on understanding the effect of martensitic transformation in the steel on the potential for cracking during seamless induction hardening (SIH) as a function of process conditions to allow the process to optimally scale up. Large-scale, three-dimensional phase-field simulations of martensitic transformation were performed using MEUMAPPS-SS (Microstructure Evolution Using Massively Parallel Phase-field Simulations – Solid State) code developed at Oak Ridge National Laboratory. The simulations were guided by location-specific thermal history generated by experimental measurements of time-temperature history generated at The Timken Company. The simulations were able to capture the morphological evolution of the martensite variants in an Fe-1.0C-1.5Cr steel based on the Nishiyama-Wasserman (NW) orientation relationship. The simulations were also able to quantify the stress-state at the interface between impinging martensite variants. The simulations indicated that the magnitude of the various stress and strain components were dependent on the sizes of the impinging plates with a reduction in these quantities with reduced plate size in agreement with experimental findings. The results obtained from the simulations will be used to guide the optimization of the alloy thermal conditions to eliminate quench cracking during SIH of bearing steels.
Work proposed in this project focused on understanding the effect of martensitic transformation in the steel on the potential for cracking during seamless induction hardening (SIH) as a function of process conditions to allow the process to optimally scale up. Large-scale, three-dimensional phase-field simulations of martensitic transformation were performed using MEUMAPPS-SS (Microstructure Evolution Using Massively Parallel Phase-field Simulations – Solid State) code developed at Oak Ridge National Laboratory. The simulations were guided by location-specific thermal history generated by experimental measurements of time-temperature history generated at The Timken Company. The simulations were able to capture the morphological evolution of the martensite variants in an Fe-1.0C-1.5Cr steel based on the Nishiyama-Wasserman (NW) orientation relationship. The simulations were also able to quantify the stress-state at the interface between impinging martensite variants. The simulations indicated that the magnitude of the various stress and strain components were dependent on the sizes of the impinging plates with a reduction in these quantities with reduced plate size in agreement with experimental findings. The results obtained from the simulations will be used to guide the optimization of the alloy thermal conditions to eliminate quench cracking during SIH of bearing steels.
For real-world installations of photovoltaic and other renewable energy resources, the critical design metric is the levelized cost of energy (LCOE); however, many power electronics design optimizations are performed with efficiency and power density as the primary design goals. Recent work has shown that a new LCOE-focused optimization approach can yield improved system designs balancing cost and energy generation. This paper expands the LCOE optimization approach by considering comprehensive optimization parameters, adding new modeling of inductor cost, extending the semiconductor model to include effects of losses on housing cost, and implementing a genetic algorithm to improve computation efficiency.
Commercial-scale generation of carbon-containing chemicals and fuels by means of electrochemical CO 2 reduction (CO 2 R) requires electrolyzers operating at high current densities and product selectivities. Membrane-electrode assemblies (MEAs) have been shown to be suitable for this purpose. In such devices, the cathode catalyst layer controls both the rate of CO 2 R and the distribution of products. In this study, we investigate how the ionomer-to-catalyst ratio (I:Cat), catalyst loading, and catalyst-layer thickness influence the performance of a cathode catalyst layer containing Ag nanoparticles supported on carbon. In this paper, we explore how these parameters affect the cell performance and establish the role of the exchange solution (water vs. CsHCO 3 ) behind the anode catalyst layer in cell performance. We show that a high total current density is best achieved using an I:Cat ratio of 3 at a Ag loading of 0.01–0.1 mg Ag /cm 2 and with a 1.0 M solution of CsHCO 3 circulated behind the anode catalyst layer. For these conditions, the optimal CO partial current density depends on the voltage applied to the MEA. The work also reveals that the performance of the cathode catalyst layer is limited by a combination of the electrochemically active surface area and the degree to which mass transfer of CO 2 to the surface of the Ag nanoparticles and the transport of OH – anions away from it limit the overall catalyst activity. Hydration of the ionomer in the cathode catalyst layer is found not to be an issue when using an exchange solution. The insights gained allowed for a Ag CO 2 R MEA that operates between 200 mA/cm 2 and 1 A/cm 2 with CO faradaic efficiencies of 78–91%, and the findings and understanding gained herein should be applicable to a broad range of CO 2 R MEA-based devices.
This study proposes a graded structure via additive manufacturing for divertor and first wall blanket applications in fusion reactors. Materials were selected based on thermodynamic calculations to operate from 1100 °C at the plasma-facing level to 550 °C at the structural steel level. Conventional joining methods often lead to failures due to discrete reaction layers with significant mechanical property differences. Using laser beam-directed energy deposition (LB-DED), this study demonstrates the fabrication of a VCrTi-Gr91 steel functionally graded component through a novel process parameter optimization framework. A systematic approach included powder characterization, single-track depositions, and construction of printability maps. Near full-density specimens of each interlayer were additively manufactured, and a transition from V-based alloys to reduced activation ferritic martensitic steels was achieved. Computational material selection of interlayer alloys and thermodynamic/diffusion kinetics simulations prevented most interface incompatibilities. A brittle intermetallic formed at one interface, causing cracking, which was not predicted by current thermodynamic models. Transition alloy design approach was updated with a more recent database and a mitigation strategy has been proposed to eliminate the formation of deleterious intermetallic phases. Ultimately, LB-DED has proven effective for producing multi-material graded systems for fusion applications, with the demonstrated process parameter optimization framework applicable to various materials.
Many stationary and heavy-duty on-road natural gas fueled engines today operate under stoichiometric conditions with a three-way catalyst. The disadvantage of stoichiometric natural gas engines compared to lean-burn natural gas and diesel engines is lower efficiency, resulting primarily from lower power density and compression ratio. Exhaust gas recirculation (EGR) coupled with advanced combustion controls can enable operation with higher compression ratio and power density, which yields higher efficiency. This also results in engine operation between the limits of knock and misfire. Operation between these limits has been named controlled end gas auto-ignition (C-EGAI) and can be used to improve the brake efficiency of the engine. Various methods of cylinder pressure-based knock quantification were explored to implement C-EGAI. The indicated quantification methods are used for the implementation of a control scheme for C-EGAI with a relation to the fractional heat release due to auto-ignition. A custom EGR system was built and the effect of EGR on the performance of a stochiometric, spark ignited, natural gas engine is evaluated. C-EGAI is implemented and the optimal parameters are determined for peak performance under EGR and C-EGAI conditions. In this study, knock detection is used for the recognition, magnitude, and location of the auto-ignition events. Cylinder pressure-based knock detection was the primary method for determining the occurrence and location of knock but was also used for implementing the ignition control scheme for controlled end gas auto-ignition. The combustion intensity metric (CIM) enabled parametric ignition timing control which allowed for the creation of a relationship between fractional heat release due to auto-ignition and CIM. Many stationary and heavy-duty on-road natural gas fueled engines today operate under stoichiometric conditions with a three-way catalyst. The disadvantage of stoichiometric natural gas engines compared to lean-burn natural gas and diesel engines is lower efficiency, resulting primarily from lower power density and compression ratio. Exhaust gas recirculation (EGR) coupled with advanced combustion controls can enable operation with higher compression ratio and power density, which yields higher efficiency. This also results in engine operation between the limits of knock and misfire. Operation between these limits has been named controlled end gas auto-ignition (C-EGAI) and can be used to improve the brake efficiency of the engine. Various methods of cylinder pressure-based knock quantification were explored to implement C-EGAI. The indicated quantification methods are used for the implementation of a control scheme for C-EGAI with a relation to the fractional heat release due to auto-ignition. A custom EGR system was built and the effect of EGR on the performance of a stochiometric, spark ignited, natural gas engine is evaluated. C-EGAI is implemented and the optimal parameters are determined for peak performance under EGR and C-EGAI conditions. In this study, knock detection is used for the recognition, magnitude, and location of the auto-ignition events. Cylinder pressure-based knock detection was the primary method for determining the occurrence and location of knock but was also used for implementing the ignition control scheme for controlled end gas auto-ignition. The combustion intensity metric (CIM) enabled parametric ignition timing control which allowed for the creation of a relationship between fractional heat release due to auto-ignition and CIM. Both exhaust gas recirculation and controlled end gas auto-ignition were analyzed with a cooperative fuel research (CFR) engine modified for boosted fuel/air intake. The data was interpreted to provide a proper evaluation of unique analytical methods to quantify the results of C_EGAI and characterize the live auto-ignition events. The control variables for this method of C-EGAI were optimized with EGR conditions to generate the point of peak performance on the CFR engine under stoichiometric, spark ignited, natural gas conditions.
Not Available