Interface delamination of the thermal barrier coating subjected to local heating

To investigate the possible failure modes of the thermal barrier coating (TBC) used to protect the scramjet combustion chamber, the local heating via laser beam irradiation was utilized to simulate the service condition of high thermal flux and high temperature gradient. Firstly, the experimental method and process were described and the typical fracture morphology of the TBC under test were provided. Then, the theoretical and finite element modeling were carried out to study the temperature, deformation and stresses of the specimen when the top ceramic coat was subjected to local heating, and to demonstrate the mechanism on the failure of the TBC. It is revealed that the interface delamination shall appear and ultimately lead to the failure of the TBC under such thermal loading of local quick heating. According to the outcome of this study, the driving force of the interface delamination is influenced greatly by the key structural parameters and performance matching. Moreover, by utilizing the rules of the effects of these parameters on the fracture driving force, there is some possibility for the designer to optimize the performances of the TBC.

Modeling Nonlinear Peeling of Ductile Thin Films-Critical Assessment of Analytical Bending Models Using FE Simulations

Three analytical double-parameter criteria based on a bending model and a two-dimensional finite element analysis model are presented for the modeling of ductile thin film undergoing a nonlinear peeling process. The bending model is based on different governing parameters: (1) the interfacial fracture toughness and the separation strength, (2) the interfacial fracture toughness and the crack tip slope angle, and (3) the interfacial fracture toughness and the critical Mises effective strain of the delaminated thin film at the crack tip. Thin film nonlinear peeling under steady-state condition is solved with the different governing parameters. In addition, the peeling test problem is simulated by using the elastic-plastic finite element analysis model. A critical assessment of the three analytical bending models is made by comparison of the bending model solutions with the finite element analysis model solutions. Furthermore, through analyses and comparisons for solutions based on both the bending model and the finite element analysis model, some connections between the bending model and the finite element analysis model are developed. Moreover, in the present research, the effect of different selections for cohesive zone shape on the ductile film peeling solutions is discussed.

Anomalous Thermal Stability of Nd-Fe-Co-Al Bulk Metallic Glass

The thermal stability of Nd60Fe20Co10Al10 bulk metallic glass (BMG) has been studied by differential scanning calorimetry (DSC), dynamic mechanical thermal analysis (DMTA), isochronal dilatation and compression tests. The results show that the glass transition of the BMG takes place quite gradually between about 460 and 650 K at a heating rate of 0.17 K/s. Several transformation processes are observed during continuous heating with the first crystallization process beginning at about 460 K, while massive crystallization takes place near the solidus temperature of the alloy. The positive heat of mixing between the two major constituents, Nd and Fe, and, consequently, a highly inhomogeneous composition of the attained amorphous phase are responsible for the anomalous thermal stability in this system. (C) 2002 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved.

Modeling and mesoscopic damage constitutive relation of brittle short-fiber-reinforced composites

Aimed at brittle composites reinforced by randomly distributed short-fibers with a relatively large aspect ratio, stiffness modulus and strength, a mesoscopic material model was proposed. Based on the statistical description, damage mechanisms, damage-induced anisotropy, damage rate effect and stress redistribution, the constitutive relation were derived. By taking glass fiber reinforced polypropylene polymers as an example, the effect of initial orientation distribution of fibers, damage-induced anisotropy, and damage-rate effect on macro-behaviors of composites were quantitatively analyzed. The theoretical predictions compared favorably with the experimental results.

Modeling Study on the Characteristics of Laminar and Turbulent Argon Plasma Jets Impinging Normally upon a Flat Plate in Ambient Air

Modeling study is performed to compare the flow and heat transfer characteristics of laminar and turbulent argon thermal-plasma jets impinging normally upon a flat plate in ambient air. The combined-diffusion-coefficient method and the turbulence-enhanced combined-diffusion-coefficient method are employed to treat the diffusion of argon in the argon-air mixture for the laminar and the turbulent cases, respectively. Modeling results presented include the flow, temperature and argon concentration fields, the air mass flow-rates entrained into the impinging plasma jets, and the distributions of the heat flux density on the plate surface. It is found that the formation of a radial wall jet on the plate surface appreciably enhances the mass flow rate of the ambient air entrained into the laminar or turbulent plasma impinging-jet. When the plate standoff distance is comparatively small, there exists a significant difference between the laminar and turbulent plasma impinging-jets in their flow fields due to the occurrence of a large closed recirculation vortex in the turbulent plasma impinging-jet, and no appreciable difference is found between the two types of jets in their maximum values and distributions of the heat flux density at the plate surface. At larger plate standoff distances, the effect of the plate on the jet flow fields only appears in the region near the plate, and the axial decaying-rates of the plasma temperature, axial velocity and argon mass fraction along the axis of the laminar plasma impinging-jet become appreciably less than their turbulent counterparts.

Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

Nonequilibrium microstructures and their evolution in a Fe-Cr-W-Ni-C laser clad coating

The laser-solidified microstructural and compositional characterization and phase evolution during tempering at 963 K were investigated using an analytical transmission electron microscope with energy dispersive X-ray analysis. The cladded alloy, a powder mixture of Fe, Cr, W, Ni, and C with a weight ratio of 10:5:1:1:1, was processed with a 3 kW continuous wave CO2 laser. The processing parameters were 16 mm/s beam scanning speed, 3 mm beam diameter. 2 kW laser power, and 0.3 g/s feed rate. The coating was metallurgically bonded to the substrate, with a maximum thickness of 730 mu m, a microhardness of about 860 Hv and a volumetric dilution ratio of about 6%. Microanalyses revealed that the cladded coating possessed the hypoeutectic microstructure comprising the primary dendritic gamma-austenite and interdendritic eutectic consisted of gamma-austenite and M7C3 carbide. The gamma-austenite was a non-equilibrium phase with extended solid solution of alloying elements and a great deal of defect structures, i.e. a high density of dislocations, twins, and stacking faults existed in gamma phase. During high temperature aging, in situ carbide transformation occurred of M7C3 to M23C6 and M6C. The precipitation of M23C6, MC and M2C carbides from austenite was also observed.

Optical Diagnostic and Modeling Solution Growth Process of Sodium Chlorate Crystals

Both a real time optical interferometric experiment and a numerical simulation of two-dimension non-steady state model were employed to study the growth process of aqueous sodium chlorate crystals. The parameters such as solution concentration distribution, crystal dimensions, growth rate and velocity field were obtained by both experiment and numerical simulation. The influence of earth gravity during crystal growth process was analyzed. A reasonable theory model corresponding to the present experiment is advanced. The thickness of concentration boundary layer was investigated especially. The results from the experiment and numerical simulation match well.

A Dynamic Loading Device for Suction Foundations in Centrifuge Modeling

Suction bucket foundations are widely used in the offshore platform for the exploitation of the offshore petroleum and natural gas resources. During winter seasons, ice sheets formed in Bohai Bay will impose strong impact and result in strong vibration on the platform. This paper describes a dynamic loading device developed on the geotechnical centrifuge and its application in modeling suction bucket foundation under the equivalent ice-induced vibration loadings. Some experimental results are presented. It is shown that when the loading amplitude is over a critical value, the sand at the upper part around the bucket softens or even liquefies. The excess pore pressure decreases from the upper part to the lower part of the sand foundation in vertical direction while decreases from near to far away from the bucket's side wall in the horizontal direction. Large settlements of the bucket and the sand around the bucket occur under the horizontal dynamic loading. The dynamic responses of the bucket with smaller size are heavier.

Growth of Silicon Carbide Bulk Crystals by Physical Vapor Transport Method and Modeling Efforts in the Process Optimization

Silicon carbide bulk crystals were grown in an induction-heating furnace using the physical vapor transport method. Crystal growth modeling was performed to obtain the required inert gas pressure and temperatures for sufficiently large growth rates. The SiC crystals were expanded by designing a growth chamber having a positive temperature gradient along the growth interface. The obtained 6H-SiC crystals were cut into wafers and characterized by Raman scattering spectroscopy and X-ray diffraction, and the results showed that most parts of the crystals had good crystallographic structures.

A Discussion on Modeling Shape Memory Alloy Embedded in A Composite Laminate as Axial Force and Elastic Foundation

In this paper, the possible error sources of the composite natural frequencies due to modeling the shape memory alloy (SMA) wire as an axial force or an elastic foundation and anisotropy are discussed. The great benefit of modeling the SMA wire as an axial force and an elastic foundation is that the complex constitutive relation of SMA can be avoided. But as the SMA wire and graphite-epoxy are rigidly bonded together, such constraint causes the re-distribution of the stress in the composite. This, together with anisotropy, which also reduces the structural stiffness can cause the relatively large error between the experimental data and theoretical results.

Effects of Subgrid-Scale Modeling on Time Correlations in Large Eddy Simulation

The effects of the unresolved subgrid-scale (SGS) motions on the energy balance of the resolved scales in large eddy simulation (LES) have been investigated actively because modeling the energy transfer between the resolved and unresolved scales is crucial to constructing accurate SGS models. But the subgrid scales not only modify the energy balance, they also contribute to temporal decorrelation of the resolved scales. The importance of this effect in applications including the predictability problem and the evaluation of sound radiation by turbulent flows motivates the present study of the effect of SGS modeling on turbulent time correlations. This paper compares the two-point, two-time Eulerian velocity correlation in isotropic homogeneous turbulence evaluated by direct numerical simulation (DNS) with the correlations evaluated by LES using a standard spectral eddy viscosity. It proves convenient to express the two-point correlations in terms of spatial Fourier decomposition of the velocity field. The LES fields are more coherent than the DNS fields: their time correlations decay more slowly at all resolved scales of motion and both their integral scales and microscales are larger than those of the DNS field. Filtering alone is not responsible for this effect: in the Fourier representation, the time correlations of the filtered DNS field are identical to those of the DNS field itself. The possibility of modeling the decorrelating effects of the unresolved scales of motion by including a random force in the model is briefly discussed. The results could have applications to the problem of computing sound sources in isotropic homogeneous turbulence by LES