Analytical-Solution Based Corner Correction for Transient Thermal Measurement

2015 ◽  
Vol 137 (11) ◽  
Author(s):  
H. Jiang ◽  
W. Chen ◽  
Q. Zhang ◽  
L. He

The one-dimensional (1D) conduction analytical approaches for a semi-infinite domain, widely adopted in the data processing of transient thermal experiments, can lead to large errors, especially near a corner of solid domain. The problems could be addressed by adopting 2D/3D numerical solutions (finite element analysis (FEA) or computational fluid dynamics (CFD)) of the solid field. In addition to needing the access to a conduction solver and extra computing effort, the numerical field solution based processing methods often require extra experimental efforts to obtain full thermal boundary conditions around corners. On a more fundamental note, it would be highly preferable that the experimental data processing is completely free of any numerical solutions and associated discretization errors, not least because it is often the case that the main purposes of many experimental measurements are exactly to validate the numerical solution methods themselves. In the present work, an analytical-solution based method is developed to enable the correction of the 2D conduction errors in a corner region without using any conduction solvers. The new approach is based on the recognition that a temperature time trace in a 2D corner situation is the result of the accumulated heat conductions in both the normal and lateral directions. An equivalent semi-infinite 1D conduction temperature trace for a correct heat transfer coefficient (HTC) can then be generated by reconstructing and removing the lateral conduction component at each time step. It is demonstrated that this simple correction technique enables the use of the standard 1D conduction analysis to get the correct HTC completely analytically without any aid of CFD or FEA solutions. In addition to a transient infrared (IR) thermal measurement case, two numerical test cases of practical interest with turbine blade tip heat transfer and film cooling are used for validation and demonstration. It has been consistently shown that the errors of the conventional 1D conduction analysis in the near corner regions can be greatly reduced by the new corner correction method.

Author(s):  
W. Chen ◽  
H. Jiang ◽  
Q. Zhang ◽  
L. He

The 1D conduction analytical and semi-analytical approaches for a semi-infinite domain have been widely adopted in the data processing of transient thermal experiments. The convective heat transfer coefficient (HTC) calculated by the 1D approach contains large errors when lateral conduction effects are significant, especially near a corner of solid domain. The problems could be addressed by alternative full 3D numerical conduction analyses, which tend to be complex as well as requiring extra experimental efforts to obtain the full thermal boundary conditions around corners, in addition to an access to a 3D conduction solver (CFD or FEA). In the present work, a simple and effective method is developed to correct such errors in the near-corner region without using any commercial tool. The present approach is based on the recognition that a temperature time trace in a 2D situation is the result of the accumulated heat conductions from the normal and lateral directions respectively and summatively. An equivalent semi-infinite 1D conduction temperature trace for a correct HTC can be generated by reconstructing and removing the lateral conduction effect at each discrete time step. This simple new correction procedure enables the standard 1D conduction analysis to be properly used to get the correct HTC, completely analytically without needing any aid of CFD or FEA solutions. Two test cases of practical interest with turbine blade tip heat transfer and film cooling are used for validation and demonstration. It has been consistently shown that the errors of the conventional 1D conduction analysis in the near corner regions can be greatly reduced by the new corner correction method. The demonstrated validity, the simplicity and robustness of the present method makes it a good candidate for future applications in transient thermal experimental studies.


2018 ◽  
Vol 941 ◽  
pp. 2313-2318
Author(s):  
Jerry E. Gould

Most welding methods in use today involve heating and subsequent cooling of the substrates for joining. Not surprisingly, understanding of associated thermal cycles implicit with the various processes has been a key facet of welding research. While the tools are available for sophisticated numerical solutions, much insight can be gained from simplified analytical approaches. A wide range of joining technologies in use today can be addressed by nominal one-dimensional heat transfer analyses. These include, for example, resistance spot, flash-butt, and linear friction welding. In addressing heat transfer problems, the mathematical constructs for heat transfer are analogous to those for mass (diffusion) transfer. Not surprisingly, one dimensional heat transfer problems can be greatly simplified by adapting the Zener approximation from mass transfer. The work described here employs the Zener approximation to address the direct spot welding of aluminum to steel. The Zener approximation is used to understand heat flow progressively from the steel into the aluminum and finally the copper electrodes. The results are used to understand weld morphology and implicit cooling behavior


2021 ◽  
Vol 413 ◽  
pp. 29-46
Author(s):  
Axel von der Weth ◽  
Daniela Piccioni Koch ◽  
Frederik Arbeiter ◽  
Till Glage ◽  
Dmitry Klimenko ◽  
...  

In this work, the strategy for numerical solutions in transport processes is investigated. Permeation problems can be solved analytically or numerically by means of the Finite Difference Method (FDM), while choosing the Euler forward explicit or Euler backwards implicit formalism. The first method is the easiest and most commonly used, while the Euler backwards implicit is not yet well established and needs further development. Hereafter, a possible solution of the Crank-Nicolson algorithm is presented, which makes use of matrix multiplication and inversion, instead of the step-by-step FDM formalism. If one considers the one-dimensional diffusion case, the concentration of the elements can be expressed as a time dependent vector, which also contains the boundary conditions. The numerically stable matrix inversion is performed by the Branch and Bound (B&B) algorithm [2]. Furthermore, the paper will investigate, whether a larger time step can be used for speeding up the simulations. The stability range is investigated by eigenvalue estimation of the Euler forward and Euler backward. In addition, a third solver is considered, referred to as Combined Solver, that is made up of the last two ones. Finally, the Crank-Nicolson solver [9] is investigated. All these results are compared with the analytical solution. The solver stability is analyzed by means of the Steady State Eigenvector (SSEV), a mathematical entity which was developed ad hoc in the present work. In addition, the obtained results will be compared with the analytical solution by Daynes [6,7].


Author(s):  
H. Ma ◽  
Z. Wang ◽  
L. Wang ◽  
Q. Zhang ◽  
Z. Yang ◽  
...  

The uncertainty associated with the convective heat transfer coefficient (HTC) obtained in transient thermal measurement is often high, especially in high speed flow. The present study demonstrates that the experimental accuracy could be much improved by an actively controlled ramp heating instead of the conventional step heating approach. A general design guideline for the proposed ramp heating method is derived theoretically and further demonstrated by simulation cases. This paper also presents a detailed experimental study for transonic turbine blade tip heat transfer. Repeatable, high-resolution tip HTC contour was obtained through transient IR measurement with the proposed ramp heating method. Detailed uncertainty analysis shows that the resulting HTC uncertainty level is much lower than the experimental data currently available in the open literature. The ramp heating approach is specially recommended to the high-speed heat transfer experimental research community to improve the accuracy of the transient thermal measurement technique.


1999 ◽  
Author(s):  
M. K. Alam ◽  
H. Pasic ◽  
K. Anagurthi ◽  
R. Zhong

Abstract Quench probes have been used to collect temperature data in controlled quenching experiments; the data is then used to deduce the heat transfer coefficients in the quenching medium. The process of determination of the heat transfer coefficient at the surface is the inverse heat conduction problem, which is extremely sensitive to measurement errors. This paper reports on an experimental and theoretical study of quenching carried out to determine the surface heat flux history during a quenching process by an inverse algorithm based on an analytical solution. The algorithm is applied to experimental data from a quenching experiment. The surface heat flux is then calculated, and the theoretical curve obtained from the analytical solution is compared with experimental results. The inverse calculation appears to produce fast, stable, but approximate results. These results can be used as the initial guess to improve the efficiency of iterative numerical solutions which are sensitive to the initial guess.


Author(s):  
Hao Wang ◽  
Suresh V. Garimella ◽  
Jayathi Y. Murthy

When a liquid wets a solid wall, the extended meniscus may be divided into three regions: a non-evaporating region where liquid is adsorbed on the wall; a thin-film region where effects of long-range molecular forces (disjoining pressure) are felt; and an intrinsic meniscus region where capillary forces dominate. In this work, a simplified model based on the augmented Young-Laplace equation is developed and an analytical solution is obtained to more easily evaluate the total heat transfer in the thin-film region. The results are consistent with previously published numerical solutions. The present work is valid for a much wider range of fluid thermal conductivity than a previous analytical solution by Schonberg et al, which is only applicable for fluids with very low conductivity.


2020 ◽  
Vol 401 ◽  
pp. 1-13
Author(s):  
Luyanda Partner Ndlovu ◽  
Raseelo Joel Moitsheki

In this article, heat transfer through a moving fin with convective and radiative heat dissipation is studied. The analytical solutions are generated using the two-dimensional Differential Transform Method (2D DTM) which is an analytical solution technique that can be applied to various types of differential equations. The accuracy of the analytical solution is validated by benchmarking it against the numerical solution obtained by applying the inbuilt numerical solver in MATLAB ($pdepe$). A good agreement is observed between the analytical and numerical solutions. The effects of thermo-physical parameters, such as the Peclet number, surface emissivity coefficient, power index of heat transfer coefficient, convective-conductive parameter, radiative-conductive parameter and non-dimensional ambient temperature on non-dimensional temperature is studied and explained. Since numerous parameters are studied, the results could be useful in industrial and engineering applications.


Author(s):  
Robert L. McMasters ◽  
Filippo de Monte ◽  
Giampaolo D'Alessandro ◽  
James V. Beck

Abstract A two-dimensional transient thermal conduction problem is examined and numerical solutions to the problem generated by ANSYS and Matlab, employing the finite element method, are compared against an analytical solution. Various different grid densities and time-step combinations are used in the numerical solutions, including some as recommended by default in the ANSYS software, including coarse, medium and fine spatial grids. The transient temperature solutions from the analytical and numerical schemes are compared at four specific locations on the body and time-dependent error curves are generated for each point. Additionally, tabular values of each solution are presented for a more detailed comparison. The errors found in the numerical solutions by comparing them directly with the analytical solution vary depending primarily on the time step size used. The errors are much larger if calculated using the analytical solution at a given time as a basis of the comparison between the two solutions as opposed to using the steady-state temperature as a basis. The largest errors appear in the early time steps of the problem, which is typically the regime wherein the largest errors occur in mathematical solutions to transient conduction problems. Conversely, errors at larger values of dimensionless time are extremely small and the numerical solutions agree within one tenth of one percent of the analytical solutions at even the worst locations. In addition to difficulties during the early time values of the problem, temperatures calculated on convective boundaries or prescribed-heat-flux boundaries are locations generating larger-magnitude errors. Corners are particularly difficult locations and require finer gridding and finer time steps in order to generate a very precise solution from a numerical code. These regions are compared, using several grid densities, against the analytical solutions. The analytical solutions are, in turn, intrinsically verified to eight significant digits by comparing similar analytical solutions against one another at very small values of dimensionless time. The solution developed using the Matlab differential equation solver was found to have errors of a similar magnitude to those generated using ANSYS. Two different test cases are examined for the various numerical solutions using the selected grid densities. The first case involves steady heating on a portion of one surface for a long duration, up to a dimensionless time of 30. The second test case involves constant heating for a dimensionless time of one, immediately followed by an insulated condition on that same surface for another duration of one dimensionless time unit. Although the errors at large times were extremely small, the errors found within the short duration test were more significant.


Author(s):  
Robert McMasters ◽  
Filippo de Monte ◽  
Giampaolo D'Alessandro ◽  
James Beck

Abstract A two-dimensional transient thermal conduction problem is examined and numerical solutions to the problem generated by ANSYS and Matlab, employing the finite element (FE) method, are compared against an 'intrinsically' verified analytical solution. Various grid densities and time-step combinations are used in the numerical solutions, including some as recommended by default in the ANSYS software, including coarse, medium and fine spatial grids. The transient temperature solutions from the analytical and numerical schemes are compared at four specific locations on the body and time-dependent error curves are generated for each point. Additionally, tabular values of each solution are presented for a more detailed comparison. Two different test cases are examined for the various numerical solutions using selected grid densities. The first case involves uniform constant heating on a portion of one surface for a long duration, up to a dimensionless time of 30. The second test case still involves uniform constant heating but for a dimensionless time of one, immediately followed by an insulated condition on that same surface for another duration of one dimensionless time unit. Although the errors at large times for both ANSYS and Matlab are extremely small, the errors found within the short-duration test are more significant, in particular when the heating is suddenly set 'on'. Surprisingly, very small errors occur when the heating is suddenly set 'off'.


1986 ◽  
Vol 108 (2) ◽  
pp. 299-304 ◽  
Author(s):  
Win Aung ◽  
G. Worku

Numerical results are presented for the effects of buoyancy on the hydrodynamic and thermal parameters in the laminar, vertically upward flow of a gas in a parallel-plate channel. Solutions of the governing parabolic equations are obtained by the use of an implicit finite difference technique coupled with a marching procedure. It is found that buoyancy dramatically increases the hydrodynamic entry length but diminishes the thermal development distance. At a fixed wall temperature difference ratio, buoyancy enhances the heat transfer on the hot wall but has little impact on the cool wall heat transfer. Flow reversal is observed in some cases. Based on an analytical solution for fully developed flow, criteria for the occurrence of flow reversal are presented. The present numerical solutions yield results that asymptotically approach those from the analytical solution.


Sign in / Sign up

Export Citation Format

Share Document