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 ...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.展开更多
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 ...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.展开更多
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...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.展开更多
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 metho...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.展开更多
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 crit...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.展开更多
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 solve...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.展开更多
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 hig...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.展开更多
Electromagnetic field distribution in the vertical metal organic chemical vapour deposition (MOCVD) reactor is simulated by using the finite element method (FEM). The effects of alternating current frequency, inte...Electromagnetic field distribution in the vertical metal organic chemical vapour deposition (MOCVD) reactor is simulated by using the finite element method (FEM). The effects of alternating current frequency, intensity, coil turn number and the distance between the coil turns on the distribution of the Joule heat are analysed separately, and their relations to the value of Joule heat are also investigated. The temperature distribution on the susceptor is also obtained. It is observed that the results of the simulation are in good agreement with previous measurements.展开更多
An adaptive heat source mode is proposed to account for the keyhole effect and the characteristics of volumetric distribution along the direction of the workpiece thickness. Finite element analysis of the temperature ...An adaptive heat source mode is proposed to account for the keyhole effect and the characteristics of volumetric distribution along the direction of the workpiece thickness. Finite element analysis of the temperature field in keyhole plasma arc welding is conducted and the weld geometry is obtained. The predicted results are in agreement with the measured ones.展开更多
This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structu...This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.展开更多
This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials (FGMs). First, a Laplace transform approach is used to handle the time variable. Then, a f...This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials (FGMs). First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest's algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.展开更多
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 interpolat...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.展开更多
The present work couples isogeometric analysis(IGA)and boundary element methods(BEM)for three dimensional steady heat conduction problems with variable coefficients.The Computer-Aided Design(CAD)geometries are built b...The present work couples isogeometric analysis(IGA)and boundary element methods(BEM)for three dimensional steady heat conduction problems with variable coefficients.The Computer-Aided Design(CAD)geometries are built by subdivision surfaces,and meantime the basis functions of subdivision surfaces are employed to discretize the boundary integral equations for heat conduction analysis.Moreover,the radial integration method is adopted to transform the additional domain integrals caused by variable coefficients to the boundary integrals.Several numerical examples are provided to demonstrate the correctness and advantages of the proposed algorithm in the integration of CAD and numerical analysis.展开更多
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 el...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.展开更多
Heat conduction through conventional and interlocking building bricks with cavities was studied in this work. Heat transfer analysis was carried out using MATLAB? partial differential equation toolbox. Regular and sta...Heat conduction through conventional and interlocking building bricks with cavities was studied in this work. Heat transfer analysis was carried out using MATLAB? partial differential equation toolbox. Regular and staggered hole arrangements were studied. Results showed that four staggered holed interlocking bricks were effective in thermal resistance into the bricks and increasing the holes beyond four did not give any thermal resistance advantage. For the conventional bricks staggered holes did not give any thermal resistance advantage but the four-holed bricks were also adjudged to be effective in thermal resistance into the brick surface. Increasing the number of holes beyond four in conventional bricks did give some thermal resistivity advantage but very minimal. Structural strengths of holed bricks were not considered in this study.展开更多
In this paper, a numerical technique is proposed to obtain the solution for transient heat conduction equation of Copper. The copper element is characterized by many characteristics;the most important of which is its ...In this paper, a numerical technique is proposed to obtain the solution for transient heat conduction equation of Copper. The copper element is characterized by many characteristics;the most important of which is its high ability to conduct heat and electrical conductivity, in addition to being a flexible and malleable metal that is easy to form without being broken, making it one of the basic minerals that humans have benefited from for thousands of years, it is one of the first minerals. That has been discovered and extracted, and still plays a major role in the development of societies. The obtained solutions are compared with the available exact solutions and the obtained solutions using the finite difference method. The results indicate that the finite difference method is a highly effective method for obtaining approximate solutions for the thermal conductivity equation for copper. It is also clear from the numerical results from copper in the high conductivity of heat and electricity.展开更多
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 ...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 (<em>Rot</em> = 0<span style="font-family:Verdana, Helvetica, Arial;white-space:normal;background-color:#FFFFFF;">°</span> - 90<span style="font-family:Verdana, Helvetica, Arial;white-space:normal;background-color:#FFFFFF;">°</span>) located somewhere inside the enclosure. The boundary top wall of the enclosure is moving through uniform speed <em>U</em><sub>0</sub>. 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 ≤ <em>Pr</em> ≤ 2), Reynolds number (60 ≤ <em>Re</em> ≤ 120), Rayleigh number (<em>Ra</em> = 10<sup>3</sup>) and Hartmann number (<em>Ha</em> = 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 <em>Nu</em><sub>av</sub> of hot bottom wall and <em>θ</em><sub>av</sub> in the enclosure is demonstrated here to show the characteristics of heat transfer in the enclosure.展开更多
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 develo...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.展开更多
A wide range of welding and surface treatment processes involve the use of a heat source which is moving at a constant speed over the component. The numerical simulation of such processes implies a transient analysis ...A wide range of welding and surface treatment processes involve the use of a heat source which is moving at a constant speed over the component. The numerical simulation of such processes implies a transient analysis using a very refined mesh in order to follow properly the path of the heat source. The 3D-mesh size can be very large if one consider the welds length or the heat-treated surface size in industrial components. To reduce the computational time to acceptable values, several techniques have been investigated. The first type is to use analytical methods such as Rosenthal equations. The second type of solutions consists in performing a transient analysis using adaptive meshing. But, for a large proportion of the involved processes, practical experience demonstrates the existence of quasi steady state conditions over the major part of the heat source path. Numerical algorithms have therefore been developed to directly compute the steady temperature, metallurgical phase proportion and stress distributions. This paper gives a general overview of the different numerical methods used to simulate welding and surface treatment processes with a special emphasis on the steady state calculation. The benefits and limitations of each of them are discussed and applications are presented.展开更多
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 met...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.展开更多
基金Project supported by the National Natural Science Foundation of China (Nos. 12102043, 12072375U2241240)the Natural Science Foundation of Hunan Province of China (Nos. 2023JJ40698 and 2021JJ40710)。
文摘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.
文摘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.
基金supported by the Innovation Training Project for Students in NUAA(No.2016C-X0010-129)the Key Laboratory of Aircraft Environment Control and Life Support(NUAA),Ministry of Industry and Information Technology
文摘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.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Special Fund for Basic Scientific Research of Central Colleges of Chang’an University, China (Grant No. CHD2011JC080)
文摘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.
文摘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.
文摘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.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10574073)
文摘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.
基金Project supported by the State Key Program of National Natural Science Foundation of China (Grant No 60736033)the National Natural Science Fund of China (Grant No 60676048)
文摘Electromagnetic field distribution in the vertical metal organic chemical vapour deposition (MOCVD) reactor is simulated by using the finite element method (FEM). The effects of alternating current frequency, intensity, coil turn number and the distance between the coil turns on the distribution of the Joule heat are analysed separately, and their relations to the value of Joule heat are also investigated. The temperature distribution on the susceptor is also obtained. It is observed that the results of the simulation are in good agreement with previous measurements.
文摘An adaptive heat source mode is proposed to account for the keyhole effect and the characteristics of volumetric distribution along the direction of the workpiece thickness. Finite element analysis of the temperature field in keyhole plasma arc welding is conducted and the weld geometry is obtained. The predicted results are in agreement with the measured ones.
基金funded by National Natural Science Foundation of China(NSFC)under Grant Nos.11702238,51904202,and 11902212Nanhu Scholars Program for Young Scholars of XYNU.
文摘This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.
文摘This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials (FGMs). First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest's algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.
文摘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.
文摘The present work couples isogeometric analysis(IGA)and boundary element methods(BEM)for three dimensional steady heat conduction problems with variable coefficients.The Computer-Aided Design(CAD)geometries are built by subdivision surfaces,and meantime the basis functions of subdivision surfaces are employed to discretize the boundary integral equations for heat conduction analysis.Moreover,the radial integration method is adopted to transform the additional domain integrals caused by variable coefficients to the boundary integrals.Several numerical examples are provided to demonstrate the correctness and advantages of the proposed algorithm in the integration of CAD and numerical analysis.
文摘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.
文摘Heat conduction through conventional and interlocking building bricks with cavities was studied in this work. Heat transfer analysis was carried out using MATLAB? partial differential equation toolbox. Regular and staggered hole arrangements were studied. Results showed that four staggered holed interlocking bricks were effective in thermal resistance into the bricks and increasing the holes beyond four did not give any thermal resistance advantage. For the conventional bricks staggered holes did not give any thermal resistance advantage but the four-holed bricks were also adjudged to be effective in thermal resistance into the brick surface. Increasing the number of holes beyond four in conventional bricks did give some thermal resistivity advantage but very minimal. Structural strengths of holed bricks were not considered in this study.
文摘In this paper, a numerical technique is proposed to obtain the solution for transient heat conduction equation of Copper. The copper element is characterized by many characteristics;the most important of which is its high ability to conduct heat and electrical conductivity, in addition to being a flexible and malleable metal that is easy to form without being broken, making it one of the basic minerals that humans have benefited from for thousands of years, it is one of the first minerals. That has been discovered and extracted, and still plays a major role in the development of societies. The obtained solutions are compared with the available exact solutions and the obtained solutions using the finite difference method. The results indicate that the finite difference method is a highly effective method for obtaining approximate solutions for the thermal conductivity equation for copper. It is also clear from the numerical results from copper in the high conductivity of heat and electricity.
文摘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 (<em>Rot</em> = 0<span style="font-family:Verdana, Helvetica, Arial;white-space:normal;background-color:#FFFFFF;">°</span> - 90<span style="font-family:Verdana, Helvetica, Arial;white-space:normal;background-color:#FFFFFF;">°</span>) located somewhere inside the enclosure. The boundary top wall of the enclosure is moving through uniform speed <em>U</em><sub>0</sub>. 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 ≤ <em>Pr</em> ≤ 2), Reynolds number (60 ≤ <em>Re</em> ≤ 120), Rayleigh number (<em>Ra</em> = 10<sup>3</sup>) and Hartmann number (<em>Ha</em> = 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 <em>Nu</em><sub>av</sub> of hot bottom wall and <em>θ</em><sub>av</sub> in the enclosure is demonstrated here to show the characteristics of heat transfer in the enclosure.
基金financially supported by the Program for New Century Excellent Talents in University(No.NCET-13-0229,NCET-09-0396)the National Science & Technology Key Projects of Numerical Control(No.2012ZX04010-031,2012ZX0412-011)the National High Technology Research and Development Program("863"Program)of China(No.2013031003)
文摘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.
文摘A wide range of welding and surface treatment processes involve the use of a heat source which is moving at a constant speed over the component. The numerical simulation of such processes implies a transient analysis using a very refined mesh in order to follow properly the path of the heat source. The 3D-mesh size can be very large if one consider the welds length or the heat-treated surface size in industrial components. To reduce the computational time to acceptable values, several techniques have been investigated. The first type is to use analytical methods such as Rosenthal equations. The second type of solutions consists in performing a transient analysis using adaptive meshing. But, for a large proportion of the involved processes, practical experience demonstrates the existence of quasi steady state conditions over the major part of the heat source path. Numerical algorithms have therefore been developed to directly compute the steady temperature, metallurgical phase proportion and stress distributions. This paper gives a general overview of the different numerical methods used to simulate welding and surface treatment processes with a special emphasis on the steady state calculation. The benefits and limitations of each of them are discussed and applications are presented.
文摘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.