Interactions of Penny-Shaped Cracks in Three-Dimensional Solids

The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is developed in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.

Theory Of Phase Transformation And Reorientation In Single Crystalline Shape Memory Alloys

A constitutive model, based on an (n + 1)-phase mixture of the Mori-Tanaka average theory, has been developed for stress-induced martensitic transformation and reorientation in single crystalline shape memory alloys. Volume fractions of different martensite lattice correspondence variants are chosen as internal variables to describe microstructural evolution. Macroscopic Gibbs free energy for the phase transformation is derived with thermodynamics principles and the ensemble average method of micro-mechanics. The critical condition and the evolution equation are proposed for both the phase transition and reorientation. This model can also simulate interior hysteresis loops during loading/unloading by switching the critical driving forces when an opposite transition takes place.

An automatic step adjustment method for average power analysis technique used in fiber amplifiers

An automatic step adjustment (ASA) method for average power analysis (APA) technique used in fiber amplifiers is proposed in this paper for the first time. In comparison with the traditional APA technique, the proposed method has suggested two unique merits such as a higher order accuracy and an ASA mechanism, so that it can significantly shorten the computing time and improve the solution accuracy. A test example demonstrates that, by comparing to the APA technique, the proposed method increases the computing speed by more than a hundredfold under the same errors. By computing the model equations of erbium-doped fiber amplifiers, the numerical results show that our method can improve the solution accuracy by over two orders of magnitude at the same amplifying section number. The proposed method has the capacity to rapidly and effectively compute the model equations of fiber Raman amplifiers and semiconductor lasers. (c) 2006 Optical Society of America

Burst event detection in wall turbulence by WVITA method

Wavelet Variable Interval Time Average (WVITA) is introduced as a method incorporating burst event detection in wall turbulence. Wavelet transform is performed to unfold the longitudinal fluctuating velocity time series measured in the near wall region of a turbulent boundary layer using hot-film anemometer. This unfolding is both in time and in space simultaneously. The splitted kinetic of the longitudinal fluctuating velocity time series among different scales is obtained by integrating the square of wavelet coefficient modulus over temporal space. The time scale that related to burst events in wall turbulence passing through the fixed probe is ascertained by maximum criterion of the kinetic energy evolution across scales. Wavelet transformed localized variance of the fluctuating velocity time series at the maximum kinetic scale is put forward instead of localized short time average variance in Variable Interval Time Average (VITA) scheme. The burst event detection result shows that WVITA scheme can avoid erroneous judgement and solve the grouping problem more effectively which is caused by VITA scheme itself and can not be avoided by adjusting the threshold level or changing the short time average interval.

Thermal mismatch stress in SiC whisker reinforced aluminium composite: new measurement method by X-ray diffraction

A new X-ray diffraction method for characterising thermal mismatch stress (TMS) in SiC_w–Al composite has been developed. The TMS and thermal mismatch strain (TMSN) in SiC whiskers are considered to be axis symmetrical, and can be calculated by measuring the lattice distortion of the whiskers. Not only the average TMS in whiskers and matrix can be obtained, but the TMS components along longitudinal and radial directions in the SiC whiskers can also be deduced. Experimental results indicate that the TMS in SiC whiskers is compressive, and tensile in the aluminium matrix. The TMS and TMSN components along the longitudinal direction in the SiC whiskers are greater than those along the radial direction for a SiC_w–Al composite quenched at 500°C.

Shallow Water Model On Cubed-Sphere By Multi-Moment Finite Volume Method

A global numerical model for shallow water flows on the cubed-sphere grid is proposed in this paper. The model is constructed by using the constrained interpolation profile/multi-moment finite volume method (CIP/MM FVM). Two kinds of moments, i.e. the point value (PV) and the volume-integrated average (VIA) are defined and independently updated in the present model by different numerical formulations. The Lax-Friedrichs upwind splitting is used to update the PV moment in terms of a derivative Riemann problem, and a finite volume formulation derived by integrating the governing equations over each mesh element is used to predict the VIA moment. The cubed-sphere grid is applied to get around the polar singularity and to obtain uniform grid spacing for a spherical geometry. Highly localized reconstruction in CIP/MM FVM is well suited for the cubed-sphere grid, especially in dealing with the discontinuity in the coordinates between different patches. The mass conservation is completely achieved over the whole globe. The numerical model has been verified by Williamson's standard test set for shallow water equation model on sphere. The results reveal that the present model is competitive to most existing ones. (C) 2008 Elsevier Inc. All rights reserved.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

Information Preservation (IP) Method in Simulation

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands m ? long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( 1000 L m ? = ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( 5 L m ? = ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Multiscale Method for Turbulent-Flow Computations

A main method of predicting turbulent flows is to solve LES equations, which was called traditional LES method. The traditional LES method solves the motions of large eddies of size larger than filtering scale An while modeling unresolved scales less than Delta_n. Hughes et al argued that many shortcomings of the traditional LES approaches were associated with their inabilities to successfully differentiate between large and small scales. One may guess that a priori scale-separation would be better, because it can predict scale-interaction well compared with posteriori scale-separation. To this end, a multi-scale method was suggested to perform scale-separation computation. The primary contents of the multiscale method are l) A space average is used to differentiate scale. 2) The basic equations include the large scale equations and fluctuation equations. 3) The large-scale equations and fluctuation equations are coupled through turbulent stress terms. We use the multiscale equations of n=2, i.e., the large and small scale (LSS) equations, to simulate 3-D evolutions of a channel flow and a planar mixing layer flow Some interesting results are given.

Investigation on average void fraction for air/non-Newtonian power-law fluids two-phase flow in downward inclined pipes

The present work has been carried out to investigate on the average void fraction of gas/non-Newtonian fluids flow in downward inclined pipes. The influences of pipe inclination angle on the average void fraction were studied experimentally. A simple correlation, which incorporated the method of Vlachos et al. for gas/Newtonain fluid horizontal flow, the correction factor of Farooqi and Richardson and the pipe inclination angle, was proposed to predict the average void fraction of gas/non-Newtonian power-law stratified flow in downward inclined pipes. The correlation was based on 470 data points covering a wide range of flow rates for different systems at diverse angles. A good agreement was obtained between theory and data and the fitting results could describe the majority of the experimental data within ±20%.

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.

Mode properties of off-axis one-sided hybrid resonator for novel high average power dual-slab laser

A novel dual-slab laser with off-axis one-sided hybrid resonator is presented. The mode properties of the hybrid resonator are calculated using a fast Fourier transform method (FFT). The influence of wavefront distoration on the output beam quality is considered. Results indicate that the novel dual-slab laser is better than the normal dual-slab laser with off-axis one-sided hybrid resonator.

A novel optical subdivision method for dual-frequency interferometer

The theory of optical subdivision techniques of dual-frequency laser interferometers is stated. And a novel optical subdivision technique is proposed originally to enhance resolution of a commercial interferometer by adding some corner-cubes. Then the performance of the interferometer is tested. The interferometer resolution of 1.24 nm and the average error of below 2 nm are achieved by using the technique. The most novel of the optical subdivision technique is without lambda/4 plates. It is less sensitive to environmental changes, it has prodigious potential to improve resolution farther and it can reduce polarization mixing error. (C) 2007 Elsevier GmbH. All rights reserved.

Structural and optical properties of YAG : Ce3+ phosphors by sol-gel combustion method

High-quality Ce3+-doped Y3Al5O12 (YAG:Ce3+) phosphors were synthesized by a facile sol-gel combustion method. In this sol-gel combustion process, citric acid acts as a fuel for combustion, traps the constituent cations and reduces the diffusion length of the precursors. The XRD and FT-IR results show that YAG phase can form through sintering at 900 degrees C for 2 h. This temperature is much lower than that required to synthesize YAG phase via the solid-state reaction method. There were no intermediate phases such as YAlO3 (YAP) and Y4Al2O9 (YAM) observed in the sintering process. The average grain size of the phosphors sintered at 900-1100 degrees C is about 40 nm. With the increasing of sintering temperature, the emission intensity increases due to the improved crystalline and homogeneous distribution of Ce3+ ions. A blue shift has been observed in the Ce3+ emission spectrum of YAG:Ce3+ phosphors with increasing sintering temperatures from 900 to 1200 degrees C. It can be explained that the decrease of lattice constant affects the crystal field around Ce3+ ions. The emission intensity of 0.06Ce-doped YAG phosphors is much higher than that of the 0.04Ce and 0.02Ce ones. The red-shift at higher Ce3+ concentrations may be Ce-Ce interactions or variations in the unit cell parameters between YAG:Ce3+ and YAG. It can be concluded that the sol-gel combustion synthesis method provides a good distribution of Ce3+ activators at the molecular level in YAG matrix. (c) 2005 Elsevier B.V. All rights reserved.