Search NASA⌕ Search

SEARCH · Search NASA

Results for “part”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems

In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Investigation of Mechanical Properties of Parts Fabricated with Gas- and Water-Atomized 304L Stainless Steel Powder in the Laser Powder Bed Fusion Process

In this report the use of gas-atomized powder as the feedstock material for the Laser Powder Bed Fusion (L-PBF) process is common in the Additive Manufacturing (AM) community. Although gas-atomization produces powder with high sphericity, its relatively expensive production cost is a downside for application in AM processes. Water atomization of powder may overcome this limitation due to its low cost relative to the gas-atomization process. In this work, gas- and water-atomized 304L stainless steel powders were morphologically characterized through Scanning Electron Microscopy (SEM). The water-atomized powder had a wider particle size distribution and exhibited less sphericity. Measuring powder flowability using the Revolution Powder Analyzer (RPA) indicated that the water-atomized powder had less flowablility than the gas-atomized powder. Through examining the mechanical properties of L-PBF fabricated parts using tensile tests, the gas-atomized powder had significantly higher yield tensile strength and elongation than the water-atomized powder, however, their ultimate tensile strengths were not significantly different.

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

Manufacturing Processes for Permanent Magnets: Part I—Sintering and Casting

Abstract Permanent magnets (PMs) produce magnetic fields and maintain the field even in the presence of an opposing magnetic field. Electrical machines using permanent magnets are more efficient than those without. Currently, all known strong magnets contain rare earth (RE) elements, and they are core components of a wide range of applications including electric vehicles and wind turbines. RE elements such as Nd and Dy have become critical materials due to the growing demand and constrained supply. Improving the manufacturing process is effective in mitigating the RE criticality issue by reducing waste and improving parts consistency. In this article, the state of the industry for PM is reviewed in detail considering both the technical and economic drivers. The importance of RE elements is discussed along with their economic importance to green energy. The conventional sintering and casting manufacturing processes for commercial magnets, including Nd-Fe-B, Sm-Co, Alnico, and ferrite, are described in detail.

36 MATERIALS SCIENCE↗
DOE OSTI2021

AFRL Additive Manufacturing Modeling Series: Challenge 1, Characterization of Residual Strain Distribution in Additively-Manufactured Metal Parts Using Energy-Dispersive Diffraction

A machine component fabricated by additive manufacturing (AM) of metal powders often possesses steep residual stress gradients with significant magnitudes due to large temperature gradients that transpire within a localized area during fabrication. These processing-induced residual stresses can cause distortion of the part, and if they are of significant magnitude, they could induce cracking of the component, degrade printability, and/or diminish subsequent mechanical performance. The ability to predict these residual stresses imparted by AM is an important step in permeating AM technology for advanced manufacturing. Calibration and validation of AM process models used for prediction are, therefore, a critical step in understanding the origin and mitigating the challenges associated with residual stresses inherent to the AM process. In the present work, the residual strain distributions in components with simple geometries fabricated by a laser powder bed fusion (LBPF) process were characterized non-destructively using energy-dispersive X-ray diffraction in support of the US Air Force Research Laboratory Additive Manufacturing Modeling Challenge Series Groeber et al. (JOM 70:441-448, 2018). The measurement setup and approach are described in detail so that the data can be used as a benchmark to calibrate and validate models for the prediction of macro-scale residual stresses due to the LPBF process.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predicting Sintering Window of Binder Jet Additively Manufactured Parts Using a Coupled Data Analytics and CALPHAD Approach

Batch-to-batch variation in powder compositions for binder jet additive manufacturing (BJAM) can significantly deter defining an “ideal” sintering window for a given alloy. One way to overcome the problem is by running sintering experiments at various temperatures for each batch of the powder. However, such an approach increases the time required to achieve large-scale production of parts. The predictive capabilities of computational thermodynamic tools like CALPHAD can be leveraged to overcome the challenge, especially for binder jet additive manufacturing, since the process occurs under near-equilibrium conditions. However, calculating the sintering window using CALPHAD can be computationally expensive, considering many possible feedstock compositions within “specification”. Here, we generate high throughput CALPHAD data for nickel-based superalloys to develop machine learning models to predict the sintering window rapidly. The predictive capability of the models has been validated using published results on BJAM of Inconel 718 and 625. Further, validated models are lightweight and can be deployed in an industrial setting to get sintering window in an accelerated manner.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part II: Entropy Stability

In the hyperbolic research community, there exists the strong belief that a continuous Galerkin scheme is notoriously unstable and additional stabilization terms have to be added to guarantee stability. In the first part of the series, the application of simultaneous approximation terms for linear problems is investigated where the boundary conditions are imposed weakly. By applying this technique, the authors demonstrate that a pure continuous Galerkin scheme is indeed linearly stable if the boundary conditions are imposed in the correct way. In this work, we extend this investigation to the nonlinear case and focus on entropy conservation. Here, by switching to entropy variables, we provide an estimation of the boundary operators also for nonlinear problems, that guarantee conservation. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

MURR LEU structural and thermal hydraulics analyses: Part II – Impacts of irradiation thermo-mechanical behavior on thermal hydraulics safety analyses

A series of structural analyses have been performed to support the conversion of the University of Missouri Research Reactor (MURR) from the use of highly enriched uranium (HEU; ≥20 wt% U-235) to low-enriched uranium (LEU; <20 wt% U-235) fuel. The irradiation thermo-mechanical analysis evaluated the effects of fuel swelling, irradiation creep, thermal expansion, as well as thermal resistance from the oxide layer growth for the MURR LEU element in prototypic thermal and irradiation conditions as presented in Part I of this article. Overall, this irradiation thermo-mechanical analysis predicts smaller gap thickness reductions in previously limiting regions, and larger reductions in the middle of the outermost end channels where power density is not typically a maximum. Due to substantial differences between the channel gap reductions assumed for the previous safety analyses and those predicted by the irradiation thermo-mechanical analysis, a need to evaluate their impact on the thermal hydraulics safety analyses arose. This article presents the results from the steady-state safety analysis for normal operation as well as the two most limiting accident scenarios. The calculation models were revised in order to account for the spatial and temporal variation of the channel gap thicknesses. The results show that sufficient safety margins are still maintained for normal operation as well as during the postulated accident transients. This work provides a methodology of incorporating the irradiation thermo-mechanical behavior of plate-type fuel into the thermal hydraulics safety analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

MURR LEU structural and thermal hydraulics analyses: Part I – Preliminary irradiation thermo-mechanical behavior

The University of Missouri Research Reactor (MURR) is expected to be converted from highly enriched uranium (HEU, ≥ 20 wt% U-235) U-Al x dispersion fuel to low-enriched uranium (LEU, < 20 wt% U-235) with U-10Mo monolithic fuel. This work introduces high-fidelity irradiation thermo–mechanical (T-M) analysis of the MURR LEU focusing on changes in coolant channel gap thickness. Three-dimensional (3D) finite element (FE) models were developed to simulate the irradiation T-M behavior of the MURR LEU element with all 23 curved fuel plates, the two side plates, and the combs. It was shown that channel gap thickness changes were influenced not only by plate thickness variations due to fuel swelling and creep but also by the radial displacement of consecutive MURR LEU plates. Modeling the fuel element assembly captured side plate displacements, which were shown to reduce radial fuel plate displacements towards the convex side. The maximum local radial displacement in the element was predicted at the end of life (EOL) as 23.7 mil (602.0 µm) on the lateral centerline of plate 23 towards the convex side. The maximum stripe-averaged reduction in channel gap thickness, particularly relevant for thermal hydraulics (TH) safety analysis, was calculated as 15.9 mil (403.9 µm) in single-side heated channel 24 (the outermost channel). These results account for the thermal resistance from the oxide build-up on cladding surfaces which was shown to be up to 0.82 mil (20.8 µm) thick. It was demonstrated that accounting for oxide layer thermal resistance led to a 10 °C higher peak fuel temperature and a 4.4 mil (111.8 µm) greater maximum local radial displacement. The impact of the calculated channel gap thickness changes on the MURR LEU TH safety analysis is evaluated in Part II.

Irradiation thermo-mechanical analysis↗
DOE OSTI2025

Analysis of contrasting aerosol indirect effects in liquid water clouds over the northern part of Arabian Sea

The extensive daily statistics of aerosol properties, cloud properties, and their mutual correlations provide crucial information for better assessing future climate change. Here, in this paper, 14 years (2010–2023) of data from the Moderate Resolution Imaging Spectroradiometer (MODIS) are analyzed over the northern part of Arabian Sea (Latitude: 21°–25° N and Longitude: 62°–68° E) to assess the characteristics of aerosols and clouds and their relationships under different meteorological conditions. When aerosol optical depth (AOD) is less than ~0.7, the observations exhibit a positive correlation between AOD and cloud droplet effective radius (CDR) but negative correlations between AOD and cloud droplet number concentration (CDNC), between AOD and cloud optical depth (COD), between AOD and cloud liquid water path (CLWP), and between AOD and cloud geometrical thickness (H). The corresponding aerosol-cloud correlations change signs when the AOD values are larger than 0.7. However, the single folded positive AOD-cloud fraction (CF) relationship is observed in both AOD regimes. Similar correlations are also observed between precipitable water vapor (PWV) and CDR, CDNC, COD, H, CF and CLWP, together with a positive correlation between PWV and AOD. Further isolation of the environmental effects from aerosol effects by stratifying AOD and cloud data into different LTS and PWV bins shows that the signature of the well-known Twomey effect is observed under high LTS-high PWV conditions, while an opposite effect (anti-Twomey) is observed under low PWV conditions, regardless of LTS values. Additionally, negative correlations between AOD and COD, AOD and CLWP, and AOD and H are observed under low LTS, regardless of PWV conditions, with a slight positive correlation when AOD >0.4 under high LTS and PWV conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

High-dimensional isotomics, part 2: Observations of over 100 constraints on methionine's isotome

The abundances of different isotopic forms of a compound, or isotopologues, will vary based on its physical and chemical history. The number of isotopologues increases combinatorically with the size of a molecule, and even small molecules such as amino acids have thousands of potentially observable isotopic variants. However, due to the analytical challenges of separating and observing isotopologues, only a few dimensions of isotopic diversity are routinely measured. Overcoming these challenges requires both an experimental method to observe many isotopic properties and a theoretical framework for interpreting these experiments. In Part 1, we presented such a theoretical framework; here, we demonstrate an experimental method, which we apply to methionine. Our approach uses a Q Exactive HF Orbitrap to perform several “M + N experiments”, where a sample is ionized, a subset of its isotopologues with cardinal mass N daltons greater than the unsubstituted isotopologue is selected and fragmented, and the proportions of all detectable isotopic forms of those fragment ions are quantified. We perform M + 1, M + 2, M + 3, and M + 4 experiments of a methionine sample and standard where the sample has a 100 ‰ enrichment of 13C at the methyl carbon relative to the natural 13C abundance at that position in the standard, and is otherwise identical to the standard. We observe isotopic forms of 8 fragment ion species for each version of the M + N experiment. With the assistance of a forward model of expected mass spectra, we identify isotopic peaks for each fragment ion based on observed mass and abundance, screen these for data quality, and quantify abundances for 146 unique isotopic peaks at precisions of ≈ 0.3–3 ‰. We present our direct observations and use them to reconstruct the concentrations of 19 individual singly, doubly, and triply-substituted isotopologues; doing so gives fewer constraints and broader error bars than working with the direct observations, but may be more interpretable for some applications. We also examine possibilities for measuring additional peaks, which are primarily limited by the detection limit of the Orbitrap-IRMS method. We then suggest some possible uses of our direct measurements for chemical forensics and hypothesis testing. Furthermore, our results demonstrate the diversity of isotopic constraints currently observable and interpretable for organic molecules.

58 GEOSCIENCES↗
DOE OSTI2023

Correcting distortions of thin-walled machined parts by machine hammer peening

Thin-walled aerostructural components frequently get distorted after the machining process. Reworking to correct distortions or eventually rejecting parts significantly increases the cost. This paper proposes a new approach to correct distortions in thin-walled components by strategically applying hammer peening on target surfaces of a machined component. Aluminium alloy 7475-T7351 was chosen for this research. The study was divided in two stages. First, the residual stresses (RS) induced by four different pneumatic hammer peening conditions (modifying the stepover distance and initial offset) were characterised in a test coupon, and one of the conditions was selected for the next stage. In the second stage, a FEM model was used to predict distortions caused by machining in a representative workpiece. Then, the RS induced by hammer peening were included in an FEM model to define two hammer peening strategies (varying the coverage area) to analyse the capability to reduce distortions. Two workpieces were machined and then treated with the simulated hammer peening strategies for experimental validation. Results in the test coupon showed that pneumatic hammer peening can generate high compressive RS (-50 to –350 MPa) up to 800 μm depth, with their magnitude increasing with a reduced stepover distance. Application of hammer peening over 4% of the surface of the representative workpiece reduced the machining-induced distortions by 37%, and a coverage area of 100% led to and overcorrection by a factor of five. This confirms that hammer peening can be strategically applied (in target areas and changing the percentage of coverage) to correct low or severe distortions.

42 ENGINEERING↗
DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 2

A wide variety of high temperature experimental data obtained in this study complement the data on the oxidation of the two di-isobutylene isomers presented in Part I and offers a basis for an extensive validation of the kinetic model developed in this study. Due to the increasing importance of unimolecular decomposition reactions in high-temperature combustion, we have investigated the di-isobutylene isomers in high dilution utilizing a pyrolysis microflow reactor and detected radical intermediates and stable products using vacuum ultraviolet (VUV) synchrotron radiation and photoelectron photoion coincidence (PEPICO) spectroscopy. Additional speciation data at oxidative conditions were also recorded utilizing a plug flow reactor at atmospheric pressure in the temperature range 725-1150 K at equivalence ratios of 1.0 and 3.0 and at residence times of 0.35 s and 0.22 s, respectively. Combustion products were analyzed using gas chromatography (GC) and mass spectrometry (MS). Ignition delay time measurements for di-isobutylene were performed at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in 'air' in the temperature range 900-1400 K using a high-pressure shock-tube facility. New measurements of the laminar burning velocities of di-isobutylene/air flames are also presented. The experiments were performed using the heat flux method at atmospheric pressure and initial temperatures of 298-358 K. Moreover, data consistency was assessed with the help of analysis of the temperature and pressure dependencies of laminar burning velocity measurements, which was interpreted using an empirical power-law expression. Electronic structure calculations were performed to compute the energy barriers to the formation of many of the product species formed. The predictions of the present mechanism were found to be in adequate agreement with the wide variety of experimental measurements performed.

09 BIOMASS FUELS↗
DOE OSTI2021

The sum is greater than the parts: exploiting microbial communities to achieve complex functions

We report multi-species microbial communities are ubiquitous in nature. The widespread prevalence of these communities is due to highly elaborated interactions among their members thereby accomplishing metabolic functions that are unattainable by individual members alone. Harnessing these communal capabilities is an emerging field in biotechnology. The rational intervention of microbial communities for the purpose of improved function has been facilitated in part by developments in multi-omics approaches, synthetic biology, and computational methods. Recent studies have demonstrated the benefits of rational interventions to human and animal health as well as agricultural productivity. Emergent technologies, such as in situ modification of complex microbial community and community metabolic modeling, represent an avenue to engineer sustainable microbial communities. In this opinion, we review relevant computational and experimental approaches to study and engineer microbial communities and discuss their potential for biotechnological applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Operation of helium sub-atmospheric multistage cryogenic centrifugal compressor trains: Part 2 – Transient modeling and pump-down path selection

Low-pressure conditions required for operation of helium cryogenic systems below the normal boiling point (i.e. 4.2 K) are established through a transient process, commonly referred to as ‘pump-down’. This process is defined as the transition from pressures above atmospheric conditions to the saturation pressure which corresponds to a specified operational temperature. The FRIB 2 K system consists of five cryogenic centrifugal compressors which are operated in series. Historically, the pump-down process path has been established through empirical methods and system operator experience. Investigation into the pump-down process at FRIB aimed to develop a pump-down methodology which relies on theoretical model predictions rather than empirically developed process paths. Ensuring stable operation during the pump-down process involved application of a centrifugal compressor performance prediction model, which is described in Part 1 of this paper. Compressor performance maps can be directly used to evaluate the stability of a selected pump-down path and anticipate the overall reliability of the selected path. In conjunction with the compressor performance maps, a system pressure model was developed to estimate the transient pressure response during the pump-down process. Lastly, an explicit equation was developed to establish a mass flow rate profile for the pump-down process. Implementation of the presented methodology (including the developed models) allows for the system operator to determine a continuous pump-down path which maintains compressor stability while conforming to overall system capabilities. Altogether, the methodology presented has resulted in simplification of transient pump-down operations and increased the reliability, stability and efficiency of the pump-down process.

Compressor train control↗
DOE OSTI2024

PINN surrogate of Li-ion battery models for parameter inference, Part II: Regularization and application of the pseudo-2D model

Bayesian parameter inference is useful to improve Li-ion battery diagnostics and can help formulate battery aging models. However, it is computationally intensive and cannot be easily repeated for multiple cycles, multiple operating conditions, or multiple replicate cells. To reduce the computational cost of Bayesian calibration, numerical solvers for physics-based models can be replaced with faster surrogates. A physics-informed neural network (PINN) is developed as a surrogate for the pseudo-2D (P2D) battery model calibration. For the P2D surrogate, additional training regularization was needed as compared to the PINN single-particle model (SPM) developed in Part I. Both the PINN SPM and P2D surrogate models are exercised for parameter inference and compared to data obtained from a direct numerical solution of the governing equations. A parameter inference study highlights the ability to use these PINNs to calibrate scaling parameters for the cathode Li diffusion and the anode exchange current density. By realizing computational speed-ups of ~2250x for the P2D model, as compared to using standard integrating methods, the PINN surrogates enable rapid state-of-health diagnostics. Finally, in the low-data availability scenario, the testing error was estimated to ~2 mV for the SPM surrogate and ~10 mV for the P2D surrogate which could be mitigated with additional data.

25 ENERGY STORAGE↗
DOE OSTI2023

Effects of knock intensity measurement technique and fuel chemical composition on the research octane number (RON) of FACE gasolines: Part 2 – Effects of spark timing

The Research and Motor Octane Number (RON, MON) characterize a fuel’s knock resistance by rating the knock intensity of a sample fuel relative to that of Primary Reference Fuels (PRF) in a Cooperative Fuel Research (CFR) Engine. A fuel’s octane number is regulated to prevent damage from autoignition leading to knocking combustion in spark-ignition engines. The operational differences between the standard RON rating and modern engine operation are explored in a three-part publication series. The previous study focused on the effects of lambda and knock characterization. This second study primarily focuses on the effects of spark timing on RON determination. Following the findings from the first publication, the knock intensity was captured by the knockmeter and by the maximum amplitude of pressure oscillations (MAPO) at the lambda of peak knock intensity and stoichiometry. Knock-limited spark advance tests were conducted for a set of seven Fuels for Advanced Combustion Engines (FACE) from the Coordinating Research Council (CRC) with varying chemical composition, PRFs, and Toluene Standardization Fuels (TSFs). For retarded spark timings, pre-spark low-temperature heat release was found for low RON PRFs. Low RON PRFs also showed knocking characteristics before reaching the center of combustion suggesting that the use of knock-limited spark advance (KLSA) was preferred over the knock-limited combustion phasing. Primarily paraffinic fuels tended towards increased pressure oscillations while dominantly aromatic fuels experienced higher pressure rise rates. A MAPO-based KLSA correlated best to Octane Index at a negative K-factor suggesting beyond RON operation despite being at otherwise RON conditions. At stoichiometry, the MAPO-based KLSA did neither correlate to RON nor Octane Index. Finally, good agreement was found between KLSA-based effective RON from this study to the MAPO-based effective RON from the first study.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Investigating relationship between surface topography and emissivity of metallic additively manufactured parts

Due to the direct relationship between thermal history and mechanical behavior, in situ thermal monitoring is key in gauging quality of parts produced with additive manufacturing (AM). Accurate monitoring of temperatures in an AM process requires knowledge of environment and object parameters including object emissivity. The emissivity is dependent on several variables, including: wavelength, material composition, temperature, and surface topography. Researchers have been concerned with the thermal emissivity dependence on temperature since large ranges are seen in metal powder bed processes, but there is also an extensive range of surfaces produced by AM. This work focused on discovering what roughness characteristics control thermal emissivity through investigation of prototypic 316 stainless steel AM samples produced with a range of build conditions on a laser powder bed fusion machine. Through experimental measurements of emissivity using hemispherical directional reflectance (HDR), guided by simulations using a finite-difference time-domain (FDTD) Maxwell solver, it was found that combinations of existing roughness parameters describing both height and slope of the surface correlate well with emissivity changes. Furthermore, these parameters work well due to their apt description of surface features encouraging internal reflection, which is the phenomenon that increases emissivity when a surface falls under the geometric optical region conditions.

36 MATERIALS SCIENCE↗