A Full Predictor-Corrector Finite Element Method for the One-Dimensional Heat Equation with Time-Dependent Singularities

The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.

Galerkin-based quasi-smooth manifold element(QSME)method for anisotropic heat conduction problems in composites with complex geometry

The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.

Heat Distribution in Rectangular Fins Using Efficient Finite Element and Differential Quadrature Methods

Finite element method (FEM) and differential quadrature method (DQM) are among important numerical techniques used in engineering analyses. Usually elements are sub-divided uniformly in FEM (conventional FEM, CFEM) to obtain temperature distribution behavior in a fin or plate. Hence, extra computational complexity is needed to obtain a fair solution with required accuracy. In this paper, non-uniform sub-elements are considered for FEM (efficient FEM, EFEM) solution to reduce the computational complex-ity. Then this EFEM is applied for the solution of one-dimensional heat transfer problem in a rectangular thin fin. The obtained results are compared with CFEM and efficient DQM (EDQM), with non-uniform mesh generation). It is found that the EFEM exhibit more accurate results than CFEM and EDQM showing its potentiality.

Finite element modeling of heating phenomena of cracks excited by high-intensity ultrasonic pulses

A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by high-intensity ultrasonic pulses. In the finite element model, the high-power ultrasonic transducer is modeled by using a piezoelectric thermal-analogy method, and the dynamical interaction between both crack faces is modeled using a contact-impact theory. In the simulations, the frictional heating taking place at the crack faces is quantitatively calculated by using finite element thermal-structural coupling analysis, especially, the influences of acoustic chaos to plate vibration and crack heating are calculated and analysed in detail. Meanwhile, the related ultrasonic infrared images are also obtained experimentally, and the theoretical simulation results are in agreement with that of the experiments. The results show that, by using the theoretical method, a good simulation of dynamic interaction and friction heating process of the crack faces under non-chaotic or chaotic sound excitation can be obtained.

Increment-Dimensional Scaled Boundary Finite Element Method for Solving Transient Heat Conduction Problem

An increment-dimensional scaled boundary finite element method (ID-SBFEM) is developed to solve the transient temperature field.To improve the accuracy of SBFEM,the effect of high frequency factor on dynamic stiffness is considered,and the first-order continued fraction technique is used.After the derivation,the SBFE equations are obtained,and the dimensions of thermal conduction,the thermal capacity matrix and the vector of the right side term in the equations are doubled.An example is presented to illustrate the feasibility and good accuracy of the proposed method.

Combining the complex variable reproducing kernel particle method and the finite element method for solving transient heat conduction problems

In this paper, the complex variable reproducing kernel particle （CVRKP） method and the finite element （FE） method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.

Nodeless variable finite element method for heat transfer analysis by means of flux-based formulation and mesh adaptation

Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method, The solution accuracy is further improved by implementing an adaptive meshing technique to generaie finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems tire： （a） a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and （b） a transient heat conduction analysis in a long plate subjected to moving heat source.

Modal and Thermal Analysis of a Modified Connecting Rod of an Internal Combustion Engine Using Finite Element Method

The connecting rod is one of the most important moving components in an internal combustion engine. The present work determined the possibility of using aluminium alloy 7075 material to design and manufacture a connecting rod for weight optimisation without losing the strength of the connecting rod. It considered modal and thermal analyses to investigate the suitability of the material for connecting rod design. The parameters that were considered under the modal analysis were: total deformation, and natural frequency, while the thermal analysis looked at the temperature distribution, total heat flux and directional heat flux of the four connecting rods made with titanium alloy, grey cast iron, structural steel and aluminium 7075 alloy respectively. The connecting rod was modelled using Autodesk inventor2017 software using the calculated parameters. The steady-state thermal analysis was used to determine the induced heat flux and directional heat flux. The study found that Aluminium 7075 alloy deformed more than the remaining three other materials but has superior qualities in terms of vibrational natural frequency, total heat flux and lightweight compared to structural steel, grey cast iron and titanium alloy.

Streamline upwind finite element method for conjugate heat transfer problems

This paper presents a combined finite element method for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow. The streamline upwind finite element method is used for the analysis of thermal viscous flow in the fluid region, whereas the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the proposed method is to consistently couple heat transfer along the fluid-solid interface. Three test cases, i.e. conjugate Couette flow problem in parallel plate channel, counter-flow in heat exchanger, and conjugate natural convection in a square cavity with a conducting wall, are selected to evaluate the efficiency of the present method.

Radiation heat transfer model for complex superalloy turbine blade in directional solidification process based on finite element method

For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is developed in this study. Key technologies, such as distinguishing boundaries automatically, local matrix and lumped heat capacity matrix, are also stated. In order to analyze the effect of withdrawing rate on DS process,the solidification processes of a complex superalloy turbine blade in the High Rate Solidification(HRS) process with different withdrawing rates are simulated; and by comparing the simulation results, it is found that the most suitable withdrawing rate is determined to be 5.0 mm·min^(-1). Finally, the accuracy and reliability of the radiation heat transfer model are verified, because of the accordance of simulation results with practical process.

Finite Element Analysis of Magnetohydrodynamic Mixed Convection in a Lid-Driven Trapezoidal Enclosure Having Heated Triangular Block

A numerical research on magnetohydrodynamic mixed convection flow in a lid-driven trapezoidal enclosure at non-uniform heating of bottom wall has been studied numerically. The enclosure consists of insulated top wall and cold side walls, too. It also contains a heated triangular block (Rot = 0° - 90°) located somewhere inside the enclosure. The boundary top wall of the enclosure is moving through uniform speed U0. The geometry of the model has been represented mathematically by coupled governing equations in accordance with proper boundary conditions and then a two-dimensional Galerkin finite element based numerical approach has been adopted to solve this paper. The numerical computations have been carried out for the wide range of parameters Prandtl number (0.5 ≤ Pr ≤ 2), Reynolds number (60 ≤ Re ≤ 120), Rayleigh number (Ra = 103) and Hartmann number (Ha = 20) taking with different rotations of heated triangular block. The results have been shown in the form of streamlines, temperature patterns or isotherms, average Nusselt number and average bulk temperature of the fluid in the enclosure at non-uniform heating of bottom wall. It is also indicated that both the streamlines, isotherm patterns strongly depend on the aforesaid governing parameters and location of the triangular block but the thermal conductivity of the triangular block has a noteworthy role on the isotherm pattern lines. Moreover, the variation of Nuav of hot bottom wall and θav in the enclosure is demonstrated here to show the characteristics of heat transfer in the enclosure.

Impact of a Magnetic Dipole on Heat Transfer in Non-Conducting Magnetic Fluid Flow over a Stretching Cylinder

The thermal behavior of an electrically non-conducting magnetic liquid flowing over a stretching cylinder under the influence of a magnetic dipole is considered.The governing nonlinear differential equations are solved numerically using a finite element approach,which is properly validated through comparison with earlier results available in the literature.The results for the velocity and temperature fields are provided for different values of the Reynolds number,ferromagnetic response number,Prandtl number,and viscous dissipation parameter.The influence of some physical parameters on skin friction and heat transfer on the walls of the cylinder is also investigated.The applicability of this research to heat control in electronic devices is discussed to a certain extent.

CRITERIA FOR FINITE ELEMENT ALGORITHM OF GENERALIZED HEAT CONDUCTION EQUATION

To eliminate oscillation and overbounding of finite element solutions of classical heat conduction equation, the author and Xiao have put forward two new concepts of monotonies and have derived and proved several criteria. This idea is borrowed here to deal with generalized conduction equation and finite element criteria for eliminating oscillation and overbounding are also presented. Some new and useful conclusions are drawn.

ANALYSIS OF THE THERMOPHYSICAL PARAMETERS OF MOIST WOOD PARTICLE MATERIAL IN A COUPLED HEAT AND MASS TRANSFER PROCESS OF FREEZING BY USING FINITE ELEMENT METHOD

The coupled heat and moisture transfer in a freezing process of wood particle material was mathematically modeled in the paper. The models were interactively solved by using the numerical method(the finite element method and the finite difference method). By matching the theoretical calculation to an experiment, the nonlinear problem was analyzed and the variable thermophysical parameters concerned was evaluated. The analysis procedure and the evaluation of the parameters were presented in detail. The result of the study showed that by using the method as described in the paper, it was possible to determine the variable (with respect to temperature, moisture content and freezing state) thermophysical parameters which were unknown or difficult to measure as long as the governing equations for a considered process were available. The method can significantly reduces the experiment efforts for determining thermophysical parameters which arc very complicated to measure. The determined variable of the effective heat conductivity of wood particle material was given in the paper. The error of the numerical calculation was also estimated by the comparison with a matched experiment.

Finite element analysis of welding residual stress of aero engine blisk by controlling heat input

In order to improve aero engine performance, it is necessary to reduce the welding residual stress of aero engine blisk. In this paper, finite element method was employed to simulate electron beam welding process of blisk, in accordance with the deducing formula（ p = kh ） , the heat input is changed with the weld depth to control welding residual stress of blisk. The calculation results show that welding residual stress of blisk can be controlled effectively by reducing the heat input on the conditions of meeting the demand of weld penetration and guaranteeing the welding quality, a new theoretical method and some numerical data are provided for controlling welding residual stress of blisk.

Analysis Process of Finite Element Method in Heat Transfer through Fabrics

According to heat transfer principle and the process of solving engineering problems by finite element method, examples were given to demonstrate how finite element analysis can be used to describe transient heat transfer through fabrics. Details were given to describe how conduction and convection affect temperature distribution and heat loss during heat transfer processes by taking advantage of the quick calculation of FEA software MSC.Marc. Experimental results show good agreement with the theoretical results.

A Finite Element Analysis on MHD Free Convection Flow in Open Square Cavity Containing Heated Circular Cylinder

The problem of Magnetohydrodynamic (MHD) free convection heat transfer in a square open cavity containing a heated circular cylinder at the centre has been investigated in this work. As boundary conditions of the cavity, the left vertical wall is kept at a constant heat flux, bottom and top walls are kept at different high and low temperature respectively. The remaining side wall is open. Finite element analysis based on Galerkin weighted Residual approach is used to visualize the temperature distribution and fluid flow solving two-dimensional governing mass, momentum and energy equations for steady state, natural convection flow in presence of magnetic field in side an open square cavity. A uniformly heated circular cylinder is located at the centre of the cavity. The object of this study is to describe the effects of MHD on the thermal fields and flow in presence of such heated circular cylinder by visualization of graph. The investigations are conducted for different values of Rayleigh number (Ra) and Hartmann number (Ha). The results show that the temperature field and flow pattern are significantly dependent on the above mentioned parameters.

A New Formulation of the Scaled Boundary Finite Element Method for Heterogeneous Media:Application to Heat Transfer Problems

The solution to heat transfer problems in two-dimensional heterogeneous media is attended based on the scaled boundary finite element method(SBFEM)coupled with equilibrated basis functions(EqBFs).The SBFEM reduces the model order by scaling the boundary solution onto the inner element.To this end,tri-lateral elements are emanated from a scaling center,followed by the development of a semi-analytical solution along the radial direction and a finite element solution along the circumferential/boundary direction.The discretization is thus limited to the boundaries of the model,and the semi-analytical radial solution is found through the solution of an eigenvalue problem,which restricts the methods’applicability to heterogeneous media.In this research,we first extracted the SBFEM formulation considering the heterogeneity of the media.Then,we replaced the semi-analytical radial solution with the EqBFs and removed the eigenvalue solution step from the SBFEM.The varying coefficients of the partial differential equation(PDE)resulting from the heterogeneity of the media are replaced by a finite series in the radial and circumferential directions of the element.A weighted residual approach is applied to the radial equation.The equilibrated radial solution series is used in the new formulation of the SBFEM.

FREEZE PROFILE AND HEAT BALANCE CALCULATION OF THE 160kA DRAINED CELL

A 2D full cell thermo-electric model of 160kA drained cell was set up using finiteelement code to calculate its freeze profile, then the drained cell model was modifiedaccording to the freeze profile computed and its heat balance was calculated. Comparedwith that of a 160kA conventional Hall-Heroult cell (H-H cell), though the melts vol-ume of the drained cell reduced greatly, the whole heat loss from it didn't drop downapparently, and an analysis was presented in the paper. On the other hand, the anode-cathode distance (ACD) of a drained cell was much less than that of a H-H cell, sothe voltage drop on it and heat produced decreased too, steps should be taken to keepa workable heat balance on a drained cell.

3D Finite Element Numerical Simulation of Residual Stresses on Electron Beam Welded BT20 Plates

A three-dimensional finite-element model (FEM) used for calculating electron beam (EB) welding temperature and stresses fields of thin plates of BT20 titanium has been developed in which the nonlinear thermophysical and thermo-mechanical properties of the material has been considered. The welding temperature field, the distributions of residual stresses in as-welded (AW) and electron beam local post-weld heat treatment (EBLPWHT) conditions have been successfully simulated. The results show that: (1) In the weld center, the maximum magnitude of residual tensile stresses of BT20 thin plates of Ti alloy is equal to 60%- 70% of its yield strength σs. (2) The residual tensile stresses in weld center can be even decreased after EBLPWHT and the longitudinal tensile stresses are decreased about 50% compared to joints in AW conditions. (3) The numerical calculating results of residual stresses by using FEM are basically in agreement with the experimental results. Combined with numerical calculating results, the effects of electron beam welding and EBLPWHT on the distribution of welding residual stresses in thin plates of BT20 have been analyzed in detail.