Thermoelastic Stresses in SiC Single Crystals Grown by the Physical Vapor Transport Method

A finite element-based thermoelastic anisotropic stress model for hexagonal silicon carbide polytype is developed for the calculation of thermal stresses in SiC crystals grown by the physical vapor transport method. The composite structure of the growing SiC crystal and graphite lid is considered in the model. The thermal expansion match between the crucible lid and SiC crystal is studied for the first time. The influence of thermal stress on the dislocation density and crystal quality is discussed.

Characters of Surface Deformation and Surface Wave in Thermal Capillary Convection

In the field of fluid mechanics, free surface phenomena is one of the most important physical processes. In the present research work, the surface deformation and surface wave caused by temperature difference of sidewalls in a rectangular cavity have been investigated. The horizontal cross-section of the container is 52 mmx42 mm, and there is a silicon oil layer of height 3.5 mm in the experimental cavity. Temperature difference between the two side walls of the cavity is increased gradually, and the flow on the liquid layer will develop from stable convection to un-stable convection. An optical diagnostic system consisting of a modified Michelson interferometer and image processor has been developed for study of the surface deformation and surface wave of thermal capillary convection. The Fourier transformation method is used to interferometer fringe analysis. The quantitative results of surface deformation and surface wave have been calculated from a serial of the interference fringe patterns.The characters of surface deformation and surface wave have been obtained. They are related with temperature gradient and surface tension. Surface deformation is fluctuant with time, which shows the character of surface wave. The cycle period of the wave is 4.8 s, and the amplitudes are from 0 to 0.55 mu m. The phase of the wave near the cool side of the cavity is opposite and correlative to that near the hot side. The present experiment proves that the surface wave of thermal capillary convection exists on liquid free surface, and it is wrapped in surface deformation.

Thermocapillary Effect in a Thermal Rheological Jet

The process of die swell in polymer jets is an important feature within polymer processing and can be explained through a study of its rheological effects. The existence of a thermocapillary effect, driven by the gradient of its surface tension, should be considered when examining a thermal jet that has a non-uniform temperature distribution on its free surface, as in various polymer processings. Both the rheological effect and thermocapillary effect on die swell can be studied numerically through a finite element method as used on a two-dimensional and unsteady model, in which a Coleman-Noll second-order fluid model is employed. The results show that the expanding angle depends on both the rheological property of the fluid and the pressure at the vessel exit. Although both the thermocapillary and the rheological effects contribute to the cross-section expansion of the fluid jet, the latter is more important in determining the expansion.

The Optical Methods on Measuring Surface Deformation and Surface Wave in The Thermal Capillary Convection

An optical diagnostic system consisting of the Michelson interferometer with the image processor has been developed for the study of the kinetics of the thermal capillary convection. The capillary convection, surface deformation, surface wave and the velocity field in a rectangular cavity with different temperature's sidewalls have been investigated by optical interference method and PIV technique. In order to calculate the surface deformation from the interference fringe, Fourier transformation is used to grating analysis. The quantitative results of the surface deformation and surface wave have been calculated from the interference fringe pattern.

A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.

Pressure and thermal effects on elastic properties of single crystal alpha-Al2O3 from Monte-Carlo simulation

A rectangular structural unit cell of a-Al2O3 is generated from its hexagonal one. For the rectangular structural crystal with a simple interatomic potential [Matsui, Mineral Mag. 58A, 571 (1994)], the relations of lattice constants to homogeneous pressure and temperature are calculated by using Monte-Carlo method at temperature 298K and 0 GPa, respectively. Both numerical results agree with experimental ones fairly well. By comparing pair distribution function, the crystal structure of a-Al2O3 has no phase transition in the range of systematic parameters. Based on the potential model, pressure dependence of isothermal bulk moduli is predicted. Under variation of general strains, which include of external and internal strains, elastic constants of a-Al2O3 in the different homogeneous load are determined. Along with increase of pressure, axial elastic constants increase appreciably, but nonaxial elastic constants are slowly changed.

Elastodynamic erosion of thermal barrier coatings

The finite element method is used to analyze the elastodynamic response of a columnar thermal barrier coating due to normal impact and oblique impact by an erosive particle. An assessment is made of the erosion by crack growth from preexisting flaws at the edge of each column: it is demonstrated that particle impacts can be sufficiently severe to give rise to columnar cracking. First, the transient stress state induced by the normal impact of a circular cylinder or a sphere is calculated in order to assess whether a 2D calculation adequately captures the more realistic 3D behavior. It is found that the transient stress states for the plane strain and axisymmetric models are similar. The sensitivity of response to particle diameter and to impact velocity is determined for both the cylinder and the sphere. Second, the transient stress state is explored for 2D oblique impact by a circular cylindrical particle and by an angular cylindrical particle. The sensitivity of transient tensile stress within the columns to particle shape (circular and angular), impact angle, impact location, orientation of the angular particle, and to the level of friction is explored in turn. The paper concludes with an evaluation of the effect of inclining the thermal barrier coating columns upon their erosion resistance. © 2011 The American Ceramic Society.

Thermal fatigue on pistons induced by shaped high power laser. Part I: Experimental study of transient temperature field and temperature oscillation

Thermal fatigue behavior is one of the foremost considerations in the design and operation of diesel engines. It is found that thermal fatigue is closely related to the temperature field and temperature fluctuation in the structure. In this paper, spatially shaped high power laser was introduced to simulate thermal loadings on the piston. The incident Gaussian beam was transformed into concentric multi-circular beam of specific intensity distribution with the help of diffractive optical element (DOE), and the transient temperature fields in the piston similar to those under working conditions could be achieved by setting up appropriate loading cycles. Simulation tests for typical thermal loading conditions, i.e., thermal high cycle fatigue (HCF) and thermal shock (or thermal low cycle fatigue, LCF) were carried out. Several important parameters that affect the transient temperature fields and/or temperature oscillations, including controlling mode, intensity distribution of shaped laser, laser power, temporal profile of laser pulse, heating time and cooling time in one thermal cycle, etc., were investigated and discussed. The results show that as a novel method, the shaped high power laser can simulate thermal loadings on pistons efficiently, and it is helpful in the study of thermal fatigue behavior in pistons. (C) 2007 Elsevier Ltd. All rights reserved.

Thermal fatigue on pistons induced by shaped high power laser. Part II: Design of spatial intensity distribution via numerical simulation

In the laser induced thermal fatigue simulation test on pistons, the high power laser was transformed from the incident Gaussian beam into a concentric multi-circular pattern with specific intensity ratio. The spatial intensity distribution of the shaped beam, which determines the temperature field in the piston, must be designed before a diffractive optical element (DOE) can be manufactured. In this paper, a reverse method based on finite element model (FEM) was proposed to design the intensity distribution in order to simulate the thermal loadings on pistons. Temperature fields were obtained by solving a transient three-dimensional heat conduction equation with convective boundary conditions at the surfaces of the piston workpiece. The numerical model then was validated by approaching the computational results to the experimental data. During the process, some important parameters including laser absorptivity, convective heat transfer coefficient, thermal conductivity and Biot number were also validated. Then, optimization procedure was processed to find favorable spatial intensity distribution for the shaped beam, with the aid of the validated FEM. The analysis shows that the reverse method incorporated with numerical simulation can reduce design cycle and design expense efficiently. This method can serve as a kind of virtual experimental vehicle as well, which makes the thermal fatigue simulation test more controllable and predictable. (C) 2007 Elsevier Ltd. All rights reserved.

Influence of thermal residual stresses on the composite macroscopic behavior

The influence of the thermal residual stress on the deformation behavior of a composite has been analyzed with a new micromechanical method. The method is based on secant moduli approximation and a new homogenized effective stress to characterize the plastic state of the matrix. It is found that the generated thermal residual stresses after cooling and their influence on the subsequent deformation behavior depends significantly on the aspect ratio of the inclusions. With prolate inclusions, the presence of thermal residual stresses generate a higher compressive hardening curves of the composite, but it is reversed with oblate inclusions. For particle reinforced composite, thermal residual stresses induce a tensile hardening curve higher than the compressive one and this is in agreement with experimental observations. (C) 1998 Elsevier Science Ltd.

Effects of contact resistance on thermal-conductivity of composite media with a periodic structure

The thermal conductivity of periodic composite media with spherical or cylindrical inclusions embedded in a homogeneous matrix is discussed. Using Green functions, we show that the Rayleigh identity can be generalized to deal with thermal properties ot these systems. A new calculating method for effective conductivity of composite media is proposed. Useful formulae for effective thermal conductivity are derived, and meanings of contact resistance in engineering problems are explained.

Thermal Residual Stresses in SiCw/Al Composite

<正> The SiCw/6061Al composites were fabricated by squeeze casting method. Varia-tions of thermal residual stresses with quenching temperature, cooling manner, aging time and thethermal-cold cycle process in thin specimens,and the distributions of thermal residual stresses alongthe distances from the surface and changes with heating temperatnres in thick specimens were stud-ied by means of X-ray diffraction (XRD). The effects of residual stresses on the mierostructure, di-mensional stability and age-hardening behavior were studied by SEM, TEM observations, and tensiletest. The results showed that there existed macrostress, microstress and thermal mismatch stress inSiCw/Al compo-site,and the presence of microstress and thermal mismatch stress has no influenceon the measurement of macrostress, but the macrostress can affect the measured value of thermalmismatch stress. Thermal res dual stress induced during the composite fabrication process, will be further in-creased when the composite were subjected to the e

Experimental Observation on Surface Wave Phenomena in Thermal Capillary Convection

An optical diagnostic system consisting of the Michelson interferometer with the image processor has been developed for studying of the surface wave in the thermal capillary convection in a rectangular cavity. In this paper, the capillary convection, surface deformation and surface wave due to the different temperature between the two sidewalls have been investigated. The cavity is 52mm?42mm in horizontal cross section and 4mm in height. The temperature difference is increased gradually and flow in liquid layer will change from steady convection to unstable convection. The optical interference method measures the surface deformation and the surface wave of the convection. The deformation of the interference fringes, which produced by the meeting of the reflected light from the liquid surface and the reference light has been captured, and the surface deformation appears when the steady convection is developed. The surface deformation is enhanced with the increasing of the temperature difference, and then several knaggy peeks in the interference fringes appear and move from the heated side to the cooled side, it demonstrates that the surface wave is existed. The surface deformation, the wavelength, the frequency, and the wave amplitude of the surface wave have been calculated.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Size- And Temperature-Dependent Thermal Expansion Coefficient Of A Nanofilm

The thermal expansion coefficient (TEC) of an ideal crystal is derived by using a method of Boltzmann statistics. The Morse potential energy function is adopted to show the dependence of the TEC on the temperature. By taking the effects of the surface relaxation and the surface energy into consideration, the dimensionless TEC of a nanofilm is derived. It is shown that with decreasing thickness, the TEC can increase or decrease, depending on the surface relaxation of the nanofilm.