Numerical Estimation of Nonlinear Thermal Conductivity of SAE 1020 Steel

Abstract

Thermally characterizing high-thermal conductivity materials is challenging, especially considering high temperatures. However, the modeling of heat transfer processes requires specific material information. The present study addresses an inverse approach to estimate the thermal conductivity of SAE 1020 relative to temperature during an autogenous LASER Beam Welding (LBW) experiment. The temperature profile during LBW is computed with the aid of an in-house CUDA-C algorithm. Here, the governing three-dimensional heat diffusion equation is discretized through the Finite Volume Method (FVM) and solved using the Successive Over-Relaxation (SOR) parallelized iterative solver. With temperature information, one may employ a minimization procedure to assess thermal properties or process parameters. In this work, the Quadrilateral Optimization Method (QOM) is applied to perform estimations because it allows for the simultaneous optimization of variables with no quantity restriction and renders the assessment of parameters in unsteady states valid, thereby preventing the requirement for steady-state experiments. We extended QOM’s prior applicability to account for more parameters concurrently. In Case I, the optimization of the three parameters that compose the second-degree polynomial function model of thermal conductivity is performed. In Case II, the heat distribution model’s gross heat rate (Ω) is also estimated in addition to the previous parameters. Ω [W] quantifies the power the sample receives and is related to the process’s efficiency. The method’s suitability for estimating the parameters was confirmed by investigating the reduced sensitivity coefficients, while the method’s stability was corroborated by performing the estimates with noisy data. There is a good agreement between the reference and estimated values. Hence, this study introduces a proper methodology for estimating a temperature-dependent thermal property and an LBW parameter. As the performance of the present algorithm is increased using parallel computation, a pondered solution between estimation reliability and computational cost is achieved.