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.

On the tensile behaviors of a hard chromium coating plated on a steel substrate with periodic subsurface inclusions

The tensile behaviors of a hard chromium coating plated on a steel substrate with periodic laser pre-quenched regions have been investigated by experimental and theoretic analysis. In the experiment, three specimens are adopted to study the differences between homogeneous and periodic inhomogeneous substrates as well as between periodic inhomogeneous substrate of relatively softer and stiffer materials. The unique characteristics have been observed in the specimen of periodic inhomogeneous substrate under quasi-static tension loading. With the periodic laser pre-quenched regions being treated as periodic subsurface inclusions (PSI), the unique stress/strain pattern of the specimen is obtained by analytical modeling and FEM analysis, and the mechanisms accounting for the experimental results is preliminarily illustrated.

Effect of dentin tubules on the mechanical properties of dentin - part III: Numerical, analysis

Imitating a real tooth and the periodontal supporting tissues, we have established a 2D finite element model and carried out a numerical analysis based on the inhomogeneous and anisotropic (IA) stress-strain relation and strength model of dentin proposed in the preceding Parts I and II, and the conventional homogeneous and isotropic (III) model, respectively. Quite a few cases of loadings for a non-defected and a defected tooth are considered. The numerical results show that the stress level predicted by the IA model is remarkably higher than that by the III model, revealing that the effect of the dentin tubules should be taken into a serious consideration from the viewpoint of biomechanics.

Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated.

The induction plasma chemical reactor: Part I. Equilibrium model

A mathematical model is presented for the numerical simulation of the flow, temperature, and concentration fields in an rf plasma chemical reactor. The simulation is performed assuming chemical equilibrium. The extent of validity of this assumption is discussed. The system considered is the reaction of SiCl4 and NH3 for the production of Si3N4.

The induction plasma chemical reactor: Part II. Kinetic model

A kinetic model has been developed for the prediction of the concentration gelds in an rf plasma reactor. A sample calculation for a SiCl4/H2 system is then performed. The model considers the mixing processes along with the kinetics of seven reactions involving the decomposition of these reactants. The results obtained are compared to those assuming chemical equilibrium. The predictions indicate that an equilibrium assumption will result in lower predicted temperature fields in the reactor. Furthermore, for the chemical system considered here, while differences exist between the concentration fields obtained by the two models, the differences are not substantial.

Propagation of Explosion Wave in Hard-Soft-Hard Layered Media

The paper presents a reasonable analysis for dynamic response and failure process of a plane multi-layered media, which are subjected to a blast loading. This blast loading is induced by a cylindric explosive put on the center of top surface of the layered media. With the help of numerical simulation technique provided by LS-DYNA software, the whole process of explosion wave propagation and attenuation can be revealed. The feature of local failure around the blasting site is also discussed in some detail. Our focus will be on the explosion wave attenuation for the hard-soft-hard sandwich layers. As seen in the paper, the computational results are delivered in a feasible way by comparing with experimental data.

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.

Analysis of micro-indentation problem of a soft film on a hard substrate

Two stages have been observed in micro-indentation experiment of a soft film on a hard substrate. In the first stage, the hardness of the thin film decreases with increasing depth of indentation when indentation is shallow; and in the second stage, the hardness of the film increases with increasing depth of indentation when the indenter tip approaches the hard substrate. In this paper, the new strain gradient theory is used to analyze the micro-indentation behavior of a soft film on a hard substrate. Meanwhile, the classic plastic theory is also applied to investigating the problem. Comparing two theoretical results with the experiment data, one can find that the strain gradient theory can describe the experiment data at both the shallow and deep indentation depths quite well, while the classic theory can't explain the experiment results.

Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Energy spectrum of ground state and excitation spectrum of quasi-particle for hard-core boson in optical lattices

We investigate the energy spectrum of ground state and quasi-particle excitation spectrum of hard-core bosons, which behave very much like spinless noninteracting fermions, in optical lattices by means of the perturbation expansion and Bogoliubov approach. The results show that the energy spectrum has a single band structure, and the energy is lower near zero momentum; the excitation spectrum gives corresponding energy gap, and the system is in Mott-insulating state at Tonks limit. The analytic result of energy spectrum is in good agreement with that calculated in terms of Green's function at strong correlation limit.

Effects of the frequency chirp on the fields of a chirped Gaussian pulse passing through a hard-edged aperture

Based on the Huygens-Fresnel diffraction integral and Fourier transform, propagation expression of a chirped Gaussian pulse passing through a hard-edged aperture is derived. Intensity distributions of the pulse with different frequency chirp in the near-field and far-field are analyzed in detail by numerical calculations. In the near-field, amplitudes of the intensity peaks generated by the modulation of the hard-edged aperture decrease with increasing the frequency chirp, which results in the improving of the beam uniformity. A physical explanation for the smoothing effect brought by increasing the frequency chirp is given. The smoothing effect is achieved not only in the pulse with Gaussian transverse profile but also in the pulse with Hermite-Gaussian transverse profile when the frequency chirp increases. (C) 2005 Elsevier B.V. All rights reserved.

Transverse intensity distributions of a broadband laser modulated by a hard-edged aperture

By means of the Huygens-Fresnel diffraction integral, the field representation of a laser beam modulated by a hard-edged aperture is derived. The near-field and far-field transverse intensity distributions of the beams with different bandwidths are analyzed by using the representation. The numerical calculation results indicate that the amplitudes and numbers of the intensity spikes decrease with increasing bandwidth, and beam smoothing is achieved when the bandwidth takes a certain value in the near field. In the far field, the radius of the transverse intensity distribution decreases as the bandwidth increases, and the physical explanation of this fact is also given. (c) 2005 Optical Society of America.