Energy conversion in shape memory alloy heat engine - Part II: Simulation

In this paper, a theoretical model proposed in Part I (Zhu et al., 2001a) is used to simulate the behavior of a twin crank NiTi SMA spring based heat engine, which has been experimentally studied by Iwanaga et al. (1988). The simulation results are compared favorably with the measurements. It is found that (1) output torque and heat efficiency decrease as rotation speed increase; (2) both output torque and output power increase with the increase of hot water temperature; (3) at high rotation speed, higher water temperature improves the heat efficiency. On the contrary, at low rotation speed, lower water temperature is more efficient; (4) the effects of initial spring length may not be monotonic as reported. According to the simulation, output torque, output power and heat efficiency increase with the decrease of spring length only in the low rotation speed case. At high rotation speed, the result might be on the contrary.

Electron Beam Surface Treatment. Part I: Surface Hardening of AISI D3 Tool Steel

The effects of electron beam surface hardening treatment on the microstructure and hardness of AISI D3 tool steel have been investigated in this paper. The results showed that the microstructure of the hardened layer consisted of martensite, a dispersion

An elastoplastic constitutive relation of whisker-reinforced composite for meso damage: Part I - formulation

An elastoplastic constitutive relation is developed for meso damage of whisker-reinforced composites. A model is constructed that includes orientation distribution of whiskers and slip systems as well as interface and crystal sliding. Evolution of damage will be addressed. Given in Part I is the formulation while examples will be illustrated in Part II.

Energy Conversion in Shape Memory Alloy Heat Engine Part I: Theory

Shape Memory Alloy (SMA) can be easily deformed to a new shape by applying a small external load at low temperature, and then recovers its original configuration upon heating. This unique shape memory phenomenon has inspired many novel designs. SMA based heat engine is one among them. SMA heat engine is an environment-friendly alternative to extract mechanical energy from low-grade energies, for instance, warm wastewater, geothermal energy, solar thermal energy, etc. The aim of this paper is to present an applicable theoretical model for simulation of SMA-based heat engines. First, a micro-mechanical constitutive model is derived for SMAs. The volume fractions of austenite and martensite variants are chosen as internal variables to describe the evolution of microstructure in SMA upon phase transition. Subsequently, the energy equation is derived based on the first thermodynamic law and the previous SMA model. From Fourier’s law of heat conduction and Newton’s law of cooling, both differential and integral forms of energy conversion equation are obtained.

Electron Beam Surface Treatment. Part II: Microstructure Evolution of Stainless Steel and Aluminum Alloy During Electron Beam Rapid Solidification

Surface rapid solidification microstructures of AISI 321 austenitic stainless steel and 2024 aluminum alloy have been investigated by electron beam remelting process and optical microscopy observation. It is indicated that the morphologies of the melted layer of both stainless steel and aluminum alloy change dramatically compared to the original materials. Also, the microstructures were greatly refined after the electron beam irradiation.

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.

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.

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.

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.