Numerical simulation of aluminum holding furnace with fluid-solid coupled heat transfer

To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.

Calculation of Mass Concrete Temperature and Creep Stress under the Influence of Local Air Heat Transfer

Temperature-induced cracking during the construction of mass concrete is a significant concern.Numerical simulations of concrete temperature have primarily assumed that the concrete is placed in an open environment.The problem of heat transfer between the air and concrete has been simplified to the concrete’s heat dissipation boundary.However,in the case of tubular concrete structures,where air inlet and outlet are relatively limited,the internal air temperature does not dissipate promptly to the external environment as it rises.To accurately simulate the temperature and creep stress in tubular concrete structures with enclosed air spaces during construction,we establish an air–concrete coupled heat transfer model according to the principles of conjugate heat transfer,and the accuracy of the model is verified through experiments.Furthermore,we conduct a case study to analyze the impact of airflow within the ship lock corridor on concrete temperature and creep stress.The results demonstrate that enhancing airflow within the corridor can significantly reduce the maximum concrete temperature.Compared with cases in which airflow within the corridor is neglected,the maximum concrete temperature and maximum tensile stress can be reduced by 12.5℃ and 0.7 MPa,respectively,under a wind speed of 4 m/s.The results of the traditional calculation method are relatively close to those obtained at a wind speed of 1 m/s.However,the temperature reduction process in the traditional method is faster,and the method yields greater tensile stress values for the corridor location.

Conjugate Heat Transfer Analysis of an Ultrasonic Molten Metal Treatment System

In piezoceramic ultrasonic devices,the piezoceramic stacks may fail permanently or function improperly if their working temperatures overstep the Curie temperature of the piezoceramic material.While the end of the horn usually serves near the melting point of the molten metal and is enclosed in an airtight chamber,so that it is difficult to experimentally measure the temperature of the transducer and its variation with time,which bring heavy difficulty to the design of the ultrasonic molten metal treatment system.To find a way out,conjugate heat transfer analysis of an ultrasonic molten metal treatment system is performed with coupled fluid and heat transfer finite element method.In modeling of the system,the RNG model and the SIMPLE algorithm are adopted for turbulence and nonlinear coupling between the momentum equation and the energy equation.Forced air cooling as well as natural air cooling is analyzed to compare the difference of temperature evolution.Numerical results show that,after about 350 s of working time,temperatures in the surface of the ceramic stacks in forced air cooling drop about 7 K compared with that in natural cooling.At 240 s,The molten metal surface emits heat radiation with a maximum rate of about 19 036 W/m2,while the heat insulation disc absorbs heat radiation at a maximum rate of about 7922 W/m2,which indicates the effectiveness of heat insulation of the asbestos pad.Transient heat transfer film coefficient and its distribution,which are difficult to be measured experimentally are also obtained through numerical simulation.At 240 s,the heat transfer film coefficient in the surface of the transducer ranges from–17.86 to 20.17 W/（m2？K）.Compared with the trial and error method based on the test,the proposed research provides a more effective way in the design and analysis of the temperature control of the molten metal treatment system.

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.

Numerical Simulation and Control of Two-Phase Flow with Evaporation in a Vertical Tube Submitted to a Conjugate Heat Transfer

A better understanding of two-phase flows with evaporation allows leading to an optimal design of evaporators. For that purpose, numerical simulations are very useful. In this paper, a numerical study has been carried out in order to model and simulate the combination of a two-phase flow with evaporation in a vertical tube. The VOF （volume-of-fluid） multiphase flow method and a phase-change model for the mass transfer have been used. For an accurate modeling, the effect of axial conduction has been also taken into account using a conjugate heat transfer model. Since thermal oscillations are undesirable as they can lead to the failure of the tube, flow instabilities have also been analyzed, using FFT （fast Fourier transforms）, in order to comprehend their behavior and influence. A control study of the flow instabilities in the tube is also presented. For that purpose tube inlet temperature has been varied using a gain control parameter.

CONJUGATE MODEL FOR HEAT AND MASS TRANSFER OF POROUS WALL IN THE HIGH TEMPERATURE GAS FLOW

Heat and mass transfer of a porous permeable wall in a high temperature gas dynamical flow is considered. Numerical simulation is conducted on the ground of the conjugate mathematical model which includes filtration and heat transfer equations in a porous body and boundary layer equations on its surface. Such an approach enables one to take into account complex interaction between heat and mass transfer in the gasdynamical flow and in the structure subjected to this flow. The main attention is given to the impact of the intraporous heat transfer intensity on the transpiration cooling efficiency.

Investigation of convection cooling guide vane with conjugate heat transfer method

This paper studied a certain blade with ten radial cooling holes which employed conjugate heat transfer method. The cooling air entered the cooling channel from the bottom of the blade and went out from the top, it was not ejected into the main flow. This paper used different numerical conditions including different turbulence models,turbulence intensities,thermal conduction coefficients and the influence on fluid property via temperature variation. The temperature distribution and pressure distribution of the blade were compared with experimental data. The results show that the numerical results using different turbulence models are almost identical to experimental data even little deviation occurs at shock wave location. The trends of temperature distribution under different numerical conditions are coincident to experimental data,especially Reynolds stress turbulence model. It can be concluded that anisotropic turbulence models can simulate the transition from laminar to turbulence,and the influence of turbulence intensity on laminar region and transition region is more than that on developed turbulent region.

Numerical simulation on conjugate heat transfer of turbine cascade

Numerical simulation on conjugate heat transfer of an internal cooled turbine vane was carried out. Numerical techniques employed included the third-order accuracy TVD scheme, multi-block structured grids and the technique of arbitrary curved mesh. Comparison between results of commercial CFD codes with several turbulence models and those of this code shows that it is incorrect of commercial CFD codes to predict the thermal boundary layer with traditional turbulence models, and that turbulence models considering transition lead to more accurate heat transfer in thermal boundary layer with some reliability and deficiency yet. The results of this code are close to those of CFX with transition model.

A new algorithm of global tightly-coupled transient heat transfer based on quasi-steady flow to the conjugate heat transfer problem

Concerning the specific demand on solving the long-term conjugate heat transfer （CHT） problem, a new algorithm of the global tightly-coupled transient heat transfer based on the quasi-steady flow field is further put forward. Compared to the traditional loosely-coupled algorithm, the computational efficiency is further improved with the greatly reduced update frequency of the flow field, and moreover the update step of the flow field can be reasonably determined by using the engineering empirical formula of the Nusselt number based on the changes of the inlet and outlet boundary conditions. Taking a duct heated by inner forced air flow heating process as an example, the comparing results to the tightly-coupled transient calculation by Fluent software shows that the new algorithm can significantly improve the computational efficiency with a reasonable accuracy on the transient temperature distribution, such as the computing time is reduced to 22,8% and 40% while the duct wall temperature deviation are 7% and 5% respectively using two flow update time step of 100 s and 50 s on the variable inlet-flow rate conditions.

Numerical Study of Conjugate Heat Transfer for Cooling the Circuit Board

In this paper, a 3D model of a flat circuit board with a heat generating electronic chip mounted on it has been studied numerically. The conjugate heat transfer including the conduction in the chip and convection with the surrounding fluid has been investigated numerically. Computational fluid dynamics using the finite volume method has been used for modeling the conjugate heat transfer through the chip and the circuit board. Conjugate heat transfer has broad applications in engineering and industrial applications in design of cooling off electronic components. Effects of various inlet velocities have been studied on the heat transfer variation and temperature of the circuit board. Numerical results show that the temperature of the chip reduces as the velocity of the inlet fluid flow increases.

Lattice Boltzmann simulation of the effects of cavity structures and heater thermal conductivity on nucleate boiling heat transfer

The boiling heat transfer technology with cavity surfaces can provide higher heat flux under lower wall superheat,which is of great significance for the cooling of electronic chips and microelectromechanical devices.In this paper,the boiling characteristics of the cavity surfaces are investigated based on the lattice Boltzmann(LB)method,focusing on the effects of cavity shapes,sizes,and heater thermal conductivity on the heat transfer performance.The results show that the triangular cavity has the best boiling performance since it has less residual vapor and higher bubble departure frequency than those of the trapezoidal and rectangular cavities.As the cavity size increases,the enhancement of heat transfer by the cavity mouth is suppressed by the heat accumulation effect at the heater bottom.The liquid rewetting process during bubble departure is the reason for the fluctuation of the space-averaged heat flux,and the heater thermal conductivity determines the fluctuation amplitude.The evaporation of liquid in the cavity with high thermal conductivity walls is more intense,resulting in shorter waiting time and higher bubble departure frequency.

Second Law Analysis of Forced Convective Cooling in a Channel with a Heated Wall Mounted Obstacle

The present work details a numerical simulation of forced convective laminar flow in a channel with a heated obstacle attached to one wall. The second law analysis is employed to investigate the distribution of entropy generation in the flow domain to demonstrate the rate of irreversibilities in thermal system. The conjugate problem including the convection heat transfer in the fluid flow and conduction one inside the obstacle is solved numerically to obtain the velocity and temperature fields in both gas and solid phases. To reach this goal, the set of governing equations including momentum and energy equations for the gas phase and conduction equation for the obstacle are solved by CFD technique to determine the hydrodynamic and thermal behaviors of the fluid flow around the obstacle and the temperature distribution in the solid element. An attempt is made to detail the local Nusselt number distribution and mean Nusselt number and also the local entropy generation distribution for the individual exposed obstacle faces. A good consistency is found between the present numerical results with experiment.

延伸冲击-扇形气膜孔复合结构耦合换热实验与数值研究

针对冲击-气膜复合冷却结构,采用实验与数值方法,对比了基于圆形气膜孔的传统冲击、延伸冲击、小距离冲击、基于扇形气膜孔的传统冲击、延伸冲击5种复合冷却系统的耦合换热和流动特性。采用红外热成像技术,获得了5种复合冷却结构在吹风比分别为0.6、1.0、1.5时外壁面的综合冷却有效度,并通过数值计算进一步揭示了内部冷却的流动和换热特征。研究结果表明:在不产生额外气动损失的前提下,延伸冲击孔结构可提升内部冷却的换热系数,进而小幅度提升壁面的综合冷却有效度,幅度为1.2%~4.6%,但随着冷气量的增大综合冷却有效度提升幅度有所减小;减小冲击距离能够提升内部冲击换热效果,但不会对综合冷却有效度产生明显影响;采用扇形气膜孔可大幅度提升外部气膜冷却性能,且提升幅度大于采用延伸冲击内部改进结构的。相较于传统的圆形气膜孔-冲击复合冷却结构,在相同冷气量条件下,基于扇形孔的延伸冲击改进方案可将壁面的面平均综合冷却有效度提高7.6%~8.5%,并将系统的流量系数提升30%以上。

基于Workbench的高压圆盘气体轴承共轭传热研究

高压圆盘气体轴承流道间隙内高速气流的对流换热与轴承圆盘内部热传导紧密耦合在一起,是一个典型的共轭传热问题。基于ANSYS Workbench工作平台的Fluid Flow(Fluent)模块对高压圆盘气体轴承进行共轭传热数值模拟,获得轴承流道间隙内的速度和压力分布、流体域与固体域的温度分布以及共轭传热时流固耦合壁面的热流密度分布,并将其与非共轭传热恒温壁面条件下的计算结果进行对比,得到高压圆盘气体轴承共轭传热的一些基本特性。结果表明:2种情况下的计算结果存在较大差异,非共轭传热恒温壁面条件下,间隙内的气体只吸热,流体域耦合壁面上的热流密度均为正值;而共轭传热条件下流体域耦合壁面热流密度存在正负值,间隙内气体的吸热和放热同时存在,显示出轴承圆盘的热传导与间隙内气体的对流换热具有复杂的共轭作用机制;相比之下,采用共轭传热模型可以得到更为符合实际的结果。研究结果为该类轴承的设计和制造提供了有益的指导。

基于参数化建模的涡轮叶片气热耦合分析方法研究

气冷涡轮的出现使得航空发动机涡轮前温度可以进一步提高,促进了高性能航空发动机的飞速发展。然而随着复杂的内部冲击、扰流强化换热、对流气膜、致密气膜喷淋等冷却形式在涡轮部件上得到应用,冷气与燃气之间的相互作用也越发复杂。为对复杂冷却涡轮叶片的流动和传热耦合影响进行快速、精细化评估,本文建立了一种参数化建模方法,针对GE E3高压涡轮一级导向器进行高保真几何建模和气热耦合仿真。分析结果表明,该参数化建模方法可实现复杂气冷涡轮叶片的高保真几何建模,获得较高的气热耦合计算精度。同时通过详细分析获得了大量的涡轮叶片内部流动换热细节,可支撑高温涡轮叶片的精细化设计。

水下大功率LED灯散热数值模拟和试验研究

为进行水下大功率LED灯散热结构的多方案优化设计,基于CFD方法对某一型自主设计、加工和制造的LED灯进行自然对流工况下的流-固共轭传热数值模拟;为验证仿真结果准确性,进行了水箱试验验证,得到的热平衡状态下测点温度值与试验结果平均相对误差为2.9%;为预测实际水下环境中LED灯的内部工作温度,进行了水下大空间中LED灯的热平衡数值模拟,比较了不同海水流速下,横向和竖向翅片对应测点的温度变化趋势,为进一步优化特定作业环境下的LED灯散热结构提供参考。

Conjugate Heat Transfer Analysis of Film Cooling Flows

The objective of this study is to evaluate the potential of various grids to satisfactorily simulate the development of a cooling film, using a coupled computation that takes into account the full geometry. Detailed computations of a single row of 30 degrees round holes on a flat plate are presented for blowing ratios of 0.764, 1.01 and 1.54. The simulation results are compared well with experimental data. The two-layer model gave more accurate results but consumed much more computational time than the standard wall functions. The k-ε turbulence model with wall functions with appropriate values of y^＋ is suitable for practical use. The results show the importance of the conjugate calculation for accurately describing the influence of the heat transfer within the cooling film.

Retrofit design of composite cooling structure of a turbine blade by fluid networks and conjugate heat transfer methods

This article discusses the development of the numerical methods of gas flow coupled with heat transfer,and introduces the fluid net-works method for rapid prediction of the performance of the composite cooling structures in turbine blade.The reliability of these methods is verified by comparing experimental data.For a HPT rotor blade,a rapid prediction on the internal cooling structures is first made by using the fluid network analysis,then an assessment of aerodynamic and heat transfer characteristics is conducted.Based on the network analysis results,three ways to improve the design of the cooling structures are tested,i.e.,adjusting the cooling gas flow mass ratios for different inner cooling cavities,reducing the flow resistances of the channel turning structures,and improving the local internal cooling structure geometries with high temperature distribution.Through the verification of full three-dimensional gas/solid/coolant conjugate heat transfer calculation,we conclude that the modified design can make the overall temperature distribution more even by significantly reducing the highest temperature of the blade surface,and reasonably matching the parameters of different coolant inlets.The results show that the proposed calculation methods can remarkably reduce the design cycle of complex turbine blade cooling structure.

BEM/FDM Conjugate Heat Transfer Analysis of a Two-dimensional Air-cooled Turbine Blade Boundary Layer

A coupled boundary element method （BEM） and finite difference method （FDM） are applied to solve conjugate heat transfer problem of a two-dimensional air-cooled turbine blade boundary layer. A loosely coupled strategy is adopted, in which each set of field equations is solved to provide boundary conditions for the other. The Navier-Stokes equations are solved by HIT-NS code. In this code, the FDM is adopted and is used to resolve the convective heat transfer in the fluid region. The BEM code is used to resolve the conduction heat transfer in the solid region. An iterated convergence criterion is the continuity of temperature and heat flux at the fluid-solid interface. The numerical results from the BEM adopted in this paper are in good agreement with the results of analytical solution and the results of commercial code, such as Fluent 6.2. The BEM avoids the complicated mesh needed in other computation method and saves the computation time. The results prove that the BEM adopted in this paper can give the same precision in numerical results with less boundary points. Comparing the conjugate results with the numerical results of an adiabatic wall flow solution, it reveals a significant difference in the distribution of metal temperatures. The results from conjugate heat transfer analysis are more accurate and they are closer to realistic thermal environment of turbines.

涡轮叶片气热耦合壁温及压力分布计算与分析

针对气冷涡轮叶片建立了耦合分析数值模拟计算平台,采用气热耦合的方法,对径向气冷MARKⅡ型叶片进行三维气热耦合数值模拟。分析了该内冷涡轮叶片的多场耦合特性,并将计算结果与试验值进行对比。对比结果表明:湍流模型的选择影响叶片表面的温度分布,但对叶片表面压力分布影响较小;在选择的6种湍流模型中,SSTk-ω湍流模型对该流场状态模拟的效果较好,叶片表面温度和压力与试验结果最接近,满足工程上的计算要求。考虑涡轮进口总温径向不均匀时,会在叶片尾缘形成局部高温区,增加了叶身的温度梯度,可适当改进冷却方式,以提高叶片强度。