Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

A dynamic coupling model for hybrid atomistic-continuum computations

A dynamic coupling model is developed for a hybrid atomistic-continuum computation in micro- and nano-fluidics. In the hybrid atomistic-continuum computation, a molecular dynamics (MD) simulation is utilized in one region where the continuum assumption breaks down and the Navier-Stokes (NS) equations are used in another region where the continuum assumption holds. In the overlapping part of these two regions, a constrained particle dynamics is needed to couple the MD simulation and the NS equations. The currently existing coupling models for the constrained particle dynamics have a coupling parameter, which has to be empirically determined. In the present work, a novel dynamic coupling model is introduced where the coupling parameter can be calculated as the computation progresses rather than inputing a priori. The dynamic coupling model is based on the momentum constraint and exhibits a correct relaxation rate. The results from the hybrid simulation on the Couette flow and the Stokes flow are in good agreement with the data from the full MD simulation and the solutions of the NS equations, respectively. (c) 2007 Elsevier Ltd. All rights reserved.

An improvement on the theoretical model for the large-scale coherent structures in the outer region of a turbulent boundary layer

According to the experimental results, there exist large-scale coherent structures in the outer region of a turbulent boundary layer, which have been studied by many authors.As experimental results, Antonia (1990) showed the phase- aver aged streamlines and isovorticity lines of the large-scale coherent structures in a turbulent boundary layer for different Reynolds numbers. Based on the hydrodynamic stability theory, the 2-D theoretical model for the large-scale structures was proposed by Luo and Zhou, in which the eddy viscosity was defined as a complex function of the position in the normal direction. The theoretical results showed in ref. were in agreement with those in ref. However, there were two problems in the results. One is that in the experimental results, there were divergent focuses between two saddle points in the streamlines, but in the theoretical results, there were centers. The other is that the stretched parts of the isovorticity lines appear at the location of centers in the theoretical results, while in the experimental results they located somewhere between the focuses and saddle points. The reason is that the computations were based on a 2-D model.

Triple-region structure for turbulent-flow in a square duct - a finite-element approach

The Reynolds-averaged Navier-Stokes equations for describing the turbulent flow in a straight square duct are formulated with two different turbulence models. The governing equations are then expanded as a multi-deck structure in a plane perpendicular to the streamwise direction, with each deck characterized by its dominant physical forces as commonly carried out in analytical work using triple-deck expansion. The resulting equations are numerically integrated using higher polynomial (H-P) finite element technique for each cross-sectional plane to be followed by finite difference representation in the streamwise direction until a fully developed state is reached. The computed results using the two different turbulence models show fair agreement with each other, and concur with the vast body of available experimental data. There is also general agreement between our results and the recent numerical works anisotropic k-epsilon turbulence model.

Theoretical and experimental researches of post-micro-buckling for fiber-reinforced composites

Post-microbuckling is a fundamental feature of compressive failure process for the unidirectional-fiber-reinforced composites and laminated composites. The post-microbuckling behavior of composites under compression in the light of the Kevlar49-reinforced 648/BF3.400 (brittle epoxy) and EP (flexible epoxy) is studied, theoretically and experimentally. Analytical results of compressive strength are in good agreement with experimental results, qualitatively and quantitatively. By the experimental research, the post-microbuckling feature of the advancing kink band model is clearly displayed.

Evolution of ideal micro-crack system

The ideal micro-cracks are treated with the number-density function; the characteristics of their evolution are investigated; a deterministic model is applied to the discussion of their extension. It is discovered that under certain conditions saturation may occur in the number-density. The main features of the statistical formulation are illustrated by several examples and compared with those observed in experiments.

Excitation mechanism of thermocapillary oscillatory convection

It is suggested that the oscillation of thermocapillary convection may be excited by the buoyancy instability. By means of numerical simulation of the finite-element method, the temperature distributions in the liquid bridge are qualitatively analyzed. The temperature gradient in a certain flow region of liquid bridge may turn to be parallel to the direction of gravity when the temperature difference △T between two boundary rods of liquid bridge is larger than the critical value. The buoyancy instability may be excited, and then the thermocapillary oscillatory convection appears, as the temperature difference increases further. The distribution of the critical Marangoni number in the micro-gravity environment is derived from the data on the ground experiments. The results show that the onset of thermocapillary oscillatory convection is delayed in the case of smaller typical scale of liquid bridge and lower gravity environment.

YAG激光毛化技术进展

二十世纪八十年代，当比利时冶金研究中心(CRM)开发出CO_2激光毛化冷轧辊技术后，尝试用YAG激光进行轧辊毛化一直吸引着众多的研究者，这是因为YAG(1.06μm)激光波长比CO_2(10.6μm)激光波长短一个量级，材料对YAG激光有更高的吸收率，并用YAG激光可以聚焦到更小的光斑尺寸，同时使用电信号驱动的声光开关技术便于对毛化分布进行可设定控制。但是用传统声光调制的YAG激光虽然可以碇以很高的脉冲频率(>30kHz)，但单脉冲有量仅为10mJ左右，难以达到辊面毛化粗糙度的要求，因此人们认为YAG激光用于毛化的主要困难是脉冲能量太小。

Formation and Structure of Composite Coating of HAD and Micro-Plasma Oxidation on A3 Steel

Composite coatings were obtained on A3 steel by hot dipping aluminum(HAD) at 720 degreesC for 6 min and micro-plasma oxidation (MPO) in alkali electrolyte. The surface morphology, element distribution and interface structure of composite coatings were studied by means of XRD, SEM and EDS. The results show that the composite coatings obtained through HAD/MPO on A3 steel consist of four layers. From the surface to the substrate, the layer is loose Al2O3 ceramic, compact Al2O3 ceramic, At and FeAl intermetallic compound layer in turn. The adhesions among all the layers are strengthened because the ceramic layer formed at the At surface originally, FeAl intermetallic compound layer and substrate are combined in metallurgical form through mutual diffusion during HAD process. Initial experiment results disclose that the anti-corrosion performance and wear resistance of composite coating are obviously improved through HAD/MPO treatment.

Edge and Size Effects in Micro Composite Structure

Because of the load transfer effect of interface layer, the stress distribution inside the composite structure of film/substrate can be very different from the Timoshenko's model. In this paper, we give the derivation and analysis of such load transfer effect of shear-lag (S-L) model. The micro-structure size (boundary conditions) effect together with interface load transfer effect becomes more and more important as the microstructure size including the three dimensions of thickness, width and length shrinks. The microstructure size is also responsible for the so-called edge-induced stress. The edge effect and difference of S-L model and Timoshenko model are also demonstrated.

The Droplet Group Micro-Explosions, Viscosity and Atomization Characteristics of W/O Emulsions

The spray of emulsified fuel, composed of diesel fuel, water and methanol can make micro-explosion under high temperature conditions, and the viscosity and the atomization characteristics of emulsion have significant effects on the micro- explosion of emulsions. To clarify the combustion mechanism of water-in-oil emulsion sprays, combustion bomb experiments were carried out, and the droplet group micro- explosions in W/O fuel emulsion sprays in a high-pressure, high-temperature bomb were observed clearly by a multi-pulsed, off-axis, image-plane ruby laser holocamera and continuously by a high-speed CCD camera.The viscosity and atomization characteristics of emulsions were also studied experimentally. The experimental results show that the higher concentration of the aqueous phase (water-methanol) (<50%) increases the viscosity of the emulsions, especially for higher agent concentration, and higher aqueous phase concentration and higher viscosity results in lager Sauter Mean Diameter (SMD). The experiment results also show that the different kinds of emulsifying agents, with different Hydrophile-Lipophile Balance (HLB) values, have significant influence on the viscosity of the emulsions.

Monte Carlo Modeling of Micro Scale Gas Flows

An information preservation (IP) method has been used to simulate many micro scale gas flows. It may efficiently reduce the statistical scatter inherent in conventional particle approaches such as the direct simulation Monte Carlo (DSMC) method. This paper reviews applications of IP to some benchmark problems. Comparison of the IP results with those given by experiment, DSMC, and the linearized Boltzmann equation, as well as the Navier-Stokes equations with a slip boundary condition, and the lattice Boltzmann equation, shows that the IP method is applicable to micro scale gas flows over the entire flow regime from continuum to free molecular.