A Molecular Dynamics Simulation: the Effect of Finite Size on the Thermal Conductivity in a Single Crystal Silicon

Non-equilibrium molecular dynamics (NEMD) simulations are performed to calculate thermal conductivity. The environment-dependent interatomic potential (EDIP) potential on crystal silicon is adopted as a model system. The issues are related to nonlinear response, local thermal equilibrium and statistical averaging. The simulation results by non-equilibrium molecular dynamics show that the calculated thermal conductivity decreases almost linearly as the film thickness reduced at the nanometre scale. The effect of size on the thermal conductivity is also obtained by a theoretic analysis of the kinetic theory and formulas of the heat capacity. The analysis reveals that the contributions of phonon mean free path (MFP) and phonon number in a finite cell to thermal conductivity are very important.

Thermal Conductivity Of Composites With Nanoscale Inclusions And Size-Dependent Percolation

An analytical model for thermal conductivity of composites with nanoparticles in a matrix is developed based on the effective medium theory by introducing the intrinsic size effect of thermal conductivity of nanoparticles and the interface thermal resistance effect between two phases. The model predicts the percolation of thermal conductivity with the volume fraction change of the second phase, and the percolation threshold depends on the size and the shape of the nanoparticles. The theoretical predictions are in agreement with the experimental results.

Size-Dependent Interface Phonon Transmission And Thermal Conductivity Of Nanolaminates

An analytical model for size-dependent interface phonon transmission and thermal conductivity of nanolaminates is derived based on the improved acoustic mismatch theory and the Lindemann melting theory by considering the size effect of phonon velocity and the interface lattice mismatch effect. The model suggests that the interface phonon transmission is dominant for the cross-plane thermal conductivity of nanolaminates and superlattices, and the intrinsic variety of size effect of thermal conductivity for different systems is proposed based on the competition mechanism of size effect of phonon transport between two materials constituting the interfaces. The model's prediction for thermal conductivity of nanolaminates agrees with the experimental results. (C) 2008 American Institute of Physics.

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.

Semi-gas kinetics model for performance modeling of flowing chemical oxygen-iodine lasers (COIL)

A semi-gas kinetics (SGK) model for performance analyses of flowing chemical oxygen-iodine laser (COIL) is presented. In this model, the oxygen-iodine reaction gas flow is treated as a continuous medium, and the effect of thermal motions of particles of different laser energy levels on the performances of the COIL is included and the velocity distribution function equations are solved by using the double-parameter perturbational method. For a premixed flow, effects of different chemical reaction systems, different gain saturation models and temperature, pressure, yield of excited oxygen, iodine concentration and frequency-shift on the performances of the COIL are computed, and the calculated output power agrees well with the experimental data. The results indicate that the power extraction of the SGK model considering 21 reactions is close to those when only the reversible pumping reaction is considered, while different gain saturation models and adjustable parameters greatly affect the output power, the optimal threshold gain range, and the length of power extraction.

Thermal stress wave and spallation induced by an electron beam

Thermal stress wave and spallation in aluminium alloy exposed to a high fluency and low energy electron beams are studied theoretically. A simple model for the study of energy deposition of electrons in materials is presented on the basis of some empirical formulae. Under the stress wave induced by energy deposition, microcracks and/or microvoids may appear in target materials, and in this case, the inelastic volume deformation should not vanish. The viscoplastic model proposed by Bodner and Partom with corresponding Gurson's yield function requires modification for this situation. The new constitutive model contains a scalar field variable description of the material damage which is taken as the void volume fraction of the polycrystalline material. Incorporation of the damage parameter permits description of rate-dependent, compressible, inelastic deformation and ductile fracture. The melting phenomenon has been observed in the experiment, therefore one needs to take into account the melting process in the intermediate energy deposition range. A three-phase equation of state used in the paper provides a more detailed and thermodynamical description of metals, particularly, in the melting region. The computational results based on the suggested model are compared with the experimental test for aluminium alloy, which is subjected to a pulsed electron beam with high fluency and low energy. (C) 1997 Elsevier Science Ltd.

The effect of thermal-activation on dislocation processes at an atomistic crack-tip

The effects of thermal activation on the dislocation emission from an atomistic crack tip are discussed, Molecular dynamics simulations at different constant temperatures are carried out to investigate the thermal effects. The simulated results show that the processes of the partial dislocation generation and emission are temperature dependent. As the temperature increases, the incipient duration of the partial dislocation nucleation becomes longer, the critical stress intensity factor for partial dislocation emission is reduced and, at the same loading level, more dislocations are emitted. The dislocation velocity moving away from the crack tip and the separations of partial dislocations are apparently not temperature dependent. The simulated results also show that, as the temperature increases, the stress distribution along the crack increases slightly. Therefore stress softening at the crack tip induced by thermal activation does not exist in the present simulation. A simple model is proposed to evaluate the relation of the critical stress intensity factor versus temperature. The obtained relation is in good agreement with our molecular dynamics results.

Influence of inertial and thermal effects on the dynamic growth of voids in porous ductile materials

The influence of inertial, thermal and rate - sensitive effects on the void growth at high strain rate in a thermal - viscoplastic solid is investigated by means of a theoretical model presented in the present paper. Numerical analysis of the model suggests that inertial, thermal and rate - sensitive effects are three major factors which greatly influence the behavior of void growth in the high strain rate case. Comparison of the mathematical model proposed in the present work and Johnson's model shows that if the temperature - dependence is considered, material viscosity eta can take the experimentally measured values.

A new model for the floating potential of fine particles in plasma

In the plasma processing of ultrafine particles of material, the heat transfer and force are considerably affected by particle charging. In this communication a new model, including thermal electron emission and incorporating the effect of electric field near the particle surface, is developed for metallic spherical particles under the condition of a thin plasma sheath. Based on this model, the particle floating potential, and thus the heat transfer and force, can be detemined more accurately and more realistically than previously.

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-conductivity of periodic composite media with spherical inclusions

The thermal conductivity of periodic composite media with spherical inclusions embedded in a homogeneous matrix is discussed. Using Green's function, we show that the Rayleigh identity can be generalized to deal with the thermal properties of these systems. A technique for calculating effective thermal conductivities is proposed. Systems with cubic symmetries (including simple cubic, body centered cubic and face centered cubic symmetry) are investigated in detail, and useful formulae for evaluating effective thermal conductivities are derived.

The thermal and mechanical structure of a two-dimensional plume in the earths mantle

On the condition that the distribution of velocity and temperature at the mid-plane of a mantle plume has been obtained (pages 213–218, this issue), the problem of determining the lateral structure of the plume at a given depth is reduced to solving an eigenvalue problem of a set of ordinary differential equations with five unknown functions, with an eigenvalue being related to the thermal thickness of the plume at this depth. The lateral profiles of upward velocity, temperature and viscosity in the plume and the thickness of the plume at various depths are calculated for two sets of Newtonian rheological parameters. The calculations show that the precondition for the existence of the plume, δT/L 1 (L = the height of the plume, δT = lateral distance from the mid-plane), can be satisfied, except for the starting region of the plume or near the base of the lithosphere. At the lateral distance, δT, the upward velocity decreases to 0.1 – 50% of its maximum value at different depths. It is believed that this model may provide an approach for a quantitative description of the detailed structure of a mantle plume.

Numeric Calculation for Transient Thermal Stress Field in Wall-Shape Metal Part During Laser Direct Forming

In order to investigate the transient thermal stress field in wall-shape metal part during laser direct forming, a FEM model basing on ANSYS is established, and its algorithm is also dealt with. Calculation results show that while the wall-shape metal part is being deposited, in X direction, the thermal stress in the top layer of the wall-shape metal part is tensile stress and in the inner of the wall-shape metal part is compressive stress. The reason causing above-mentioned thermal stress status in the wall-shape metal part is illustrated, and the influence of the time and the processing parameters on the thermal stress field in wall-shape metal part is also studied. The calculation results are consistent with experimental results in tendency.

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.

Mechanisms of thermal shock failure for ultra-high temperature ceramic

The materials considered in our analysis were ZrB2 ceramic matrix composites. Effect of two different additives (graphite and AlN) on thermal shock stability for the materials was measured by water quench test. It showed that it may provide more stable thermal shock properties with additives of graphite. It was explained by different thermal properties and crack resistance of the two materials in detail. Surface oxidation was one of main reasons for strength degradation of ceramic with additives of graphite after quenched in water, and surface crack was one of main reasons for strength degradation of ceramic with additives of AlN after quenched in water. It was presented that it was a potential method for improving thermal shock stability of ZrB2 ceramic matrix composites by introducing proper quantities of graphite.