Effects of superimposed pulse bias on TiN coating in cathodic arc deposition

Orthogonal designs are used to investigate the main factors when doing experiments in which pulse bias is superimposed on d.c. bias during cathodic are deposition of TiN. Pulse peak, duty cycle, frequency, direct voltage, are current and pressure all are investigated when coating TiN on HSS substrates. Roughness, surface micrograph, microhardness and thickness are tested. By analysis of variance, it is shown that pressure and frequency are the main factors. R-a and droplet density of the film with (d.c. + pulse) bias decrease. A simple explanation for the result is suggested.

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.

Practical matching of principal stress field geometries in composite components

A simple composite design methodology has been developed from the basic principles of composite component failure. This design approach applies the principles of stress field matching to develop suitable reinforcement patterns around three-dimensional details such as lugs in mechanical components. The resulting patterns are essentially curvilinear orthogonal meshes, adjusted to meet the restrictions imposed by geometric restraints and the intended manufacturing process. Whilst the principles behind the design methodology can be applied to components produced by differing manufacturing processes, the results found from looking at simple generic example problems suggest a realistic and practical generic manufacturing approach. The underlying principles of the design methodology are described and simple analyses are used to help illustrate both the methodology and how such components behave. These analyses suggest it is possible to replace high-strength steel lugs with composite components whose strength-to-weight ratio is some 4-5 times better. © 1998 Elsevier Science Ltd. All rights reserved.

The ensemble statistics of the energy of a random system subjected to harmonic excitation

This paper is concerned with the ensemble statistics of the response to harmonic excitation of a single dynamic system such as a plate or an acoustic volume. Random point process theory is employed, and various statistical assumptions regarding the system natural frequencies are compared, namely: (i) Poisson natural frequency spacings, (ii) statistically independent Rayleigh natural frequency spacings, and (iii) natural frequency spacings conforming to the Gaussian orthogonal ensemble (GOE). The GOE is found to be the most realistic assumption, and simple formulae are derived for the variance of the energy of the system under either point loading or rain-on-the-roof excitation. The theoretical results are compared favourably with numerical simulations and experimental data for the case of a mass loaded plate. © 2003 Elsevier Ltd. All rights reserved.

Static Study of Cantilever Beam stiction under Electrostatic Force Influence

The model and analysis of the cantilever beam adhesion problem under the action of electrostatic force are given. Owing to the nonlinearity of electrostatic force, the analytical solution for this kind of problem is not available. In this paper, a systematic method of generating polynomials which are the exact beamsolutions of the loads with different distributions is provided. The polynomials are used to approximate the beam displacement due to electrostatic force. The equilibrium equation offers an answer to how the beam deforms but no information about the unstuck length. The derivative of the functional with respect to the unstuck length offers such information. But to compute the functional it is necessary to know the beam deformation. So the problem is iteratively solved until the results are converged. Galerkin and Newton-Raphson methods are used to solve this nonlinear problem. The effects of dielectric layer thickness and electrostatic voltage on the cantilever beamstiction are studied.The method provided in this paper exhibits good convergence. For the adhesion problem of cantilever beam without electrostatic voltage, the analytical solution is available and is also exactly matched by the computational results given by the method presented in this paper.

Low-Dimensional Dynamical Systems for a Wake-Type Shear Flow with Vortex Dislocations

Low-dimensional systems are constructed to investigate dynamics of vortex dislocations in a wake-type shear flow. High-resolution direct numerical simulations are employed to obtain flow snapshots from which the most energetic modes are extracted using proper orthogonal decomposition (POD). The first 10 modes are classified into two groups. One represents the general characteristics of two-dimensional wake-type shear flow, and the other is related to the three-dimensional properties or non-uniform characteristics along the span. Vortex dislocations are generated by these two kinds of coherent structures. The results from the first 20 three-dimensional POD modes show that the low- dimensional systems have captured the basic properties of the wake-type shear flow with vortex dislocation, such as two incommensurable frequencies and their beat frequency.

Square wine glass mode resonator with quality factor of 4 million

We report on the experimental characterization of a single crystal silicon square-plate microresonator. The resonator is excited in the square wine glass (SWG) mode at a mechanical resonance frequency of 2.065 MHz. The resonator displays quality factor of 9660 in air and an ultra-high quality factor of Q = 4.05 × 106 in 12 mtorr vacuum. The SWG mode may be described as a square plate that contracts along one axis in the fabrication plane, while simultaneously extending along an orthogonal axis in the same plane. The resonant structure is addressed in a 2-terminal configuration by utilizing equal and opposite drive polarities on surrounding capacitor electrodes, thereby decreasing the motional resistance of the resonator. The resonant micromechanical device has been fabricated in a commercial silicon-on-insulator process through the MEMSCAP foundry utilising a minimum electrostatic gap of 2 μm. © 2008 IEEE.

Dynamic equations for curved submerged floating tunnel

In virtue of reference Cartesian coordinates, geometrical relations of spatial curved structure are presented in orthogonal curvilinear coordinates. Dynamic equations for helical girder are derived by Hamilton principle. These equations indicate that four generalized displacements are coupled with each other. When spatial structure degenerates into planar curvilinear structure, two generalized displacements in two perpendicular planes are coupled with each other. Dynamic equations for arbitrary curvilinear structure may be obtained by the method used in this paper.

Composite materials dynamic fracture studies by generalized Shmuely difference algorithm

The generalized Shmuely Difference Algorithm (GSDA) is presented here to analyze the dynamic fracture performance of orthogonal-anisotropic composite materials, such as glass fibre reinforced phenolplast. The difference recurrence Formulae and boundary condition difference extrapolation formulae are derived and programmed. The dynamic stress intensity factors (DSIF) of the isotropic and anisotropic centrally cracked plates are computed respectively using GSDA and compared with that published previously. GSDA is proved effective and reliable. Copyright (C) 1996 Elsevier Science Ltd.

An optimal theory for an expansion of flow quantities to capture the flow structures

An optimal theory on how database analysis to capture the flow structures has been developed in this paper, which include the POD method as its special case. By means of the remainder minimization method in the Sobolev space, for more general optimal conditions the new theory has the potential to overcome an inherent limitation of the POD method, i.e., it cannot be used to the situations in which the optimal condition is other than the inner product global one. As an example, using the new theory, the database of a two-dimensional flow over a backward-facing step is analyzed in detail, with velocity and vorticity bases.

Bending solution of high-order refined shear deformation-theory for rectangular composite plates

A new high-order refined shear deformation theory based on Reissner's mixed variational principle in conjunction with the state- space concept is used to determine the deflections and stresses for rectangular cross-ply composite plates. A zig-zag shaped function and Legendre polynomials are introduced to approximate the in-plane displacement distributions across the plate thickness. Numerical results are presented with different edge conditions, aspect ratios, lamination schemes and loadings. A comparison with the exact solutions obtained by Pagano and the results by Khdeir indicates that the present theory accurately estimates the in-plane responses.

3D Finite Volume Scheme for Czochralski Crystal Growth

A general three-dimensional model is developed for simulation of the growth process of silicon single crystals by Czochralski technique. The numerical scheme is based on the curvilinear non-orthogonal finite volume discretization. Numerical solutions show that the flow and temperature fields in the melt are asymmetric and unsteady for 8’’ silicon growth. The effects of rotation of crystal on the flow structure are studied. The rotation of crystal forms the Ekman layer in which the temperature gradient along solid/melt surface is small.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

Daily scale wintertime sea surface temperature and IPC-Navidad variability in the southern Bay of Biscay from 1981 to 2010

The combination of remotely sensed gappy Sea surface temperature (SST) images with the missing data filling DINEOF (data interpolating empirical orthogonal functions) technique, followed by a principal component analysis of the reconstructed data, has been used to identify the time evolution and the daily scale variability of the wintertime surface signal of the Iberian Poleward Current (IPC), or Navidad, during the 1981-2010 period. An exhaustive comparison with the existing bibliography, and the vertical temperature and salinity profiles related to its extremes over the Bay of Biscay area, show that the obtained time series accurately reflect the IPC-Navidad variability. Once a time series for the evolution of the SST signal of the current over the last decades is well established, this time series is used to propose a physical mechanism in relation to the variability of the IPC-Navidad, involving both atmospheric and oceanic variables. According to the proposed mechanism, an atmospheric circulation anomaly observed in both the 500 hPa and the surface levels generates atmospheric surface level pressure, wind-stress and heat-flux anomalies. In turn, those surface level atmospheric anomalies induce mutually coherent SST and sea level anomalies over the North Atlantic area, and locally, in the Bay of Biscay area. These anomalies, both locally over the Bay of Biscay area and over the North Atlantic, are in agreement with several mechanisms that have separately been related to the variability of the IPC-Navidad, i.e. the south-westerly winds, the joint effect of baroclinicity and relief (JEBAR) effect, the topographic beta effect and a weakened North Atlantic gyre.

Downscaling of surface moisture flux and precipitation in the Ebro Valley (Spain) using analogues and analogues followed by random forests and multiple linear regression

In this paper, reanalysis fields from the ECMWF have been statistically downscaled to predict from large-scale atmospheric fields, surface moisture flux and daily precipitation at two observatories (Zaragoza and Tortosa, Ebro Valley, Spain) during the 1961-2001 period. Three types of downscaling models have been built: (i) analogues, (ii) analogues followed by random forests and (iii) analogues followed by multiple linear regression. The inputs consist of data (predictor fields) taken from the ERA-40 reanalysis. The predicted fields are precipitation and surface moisture flux as measured at the two observatories. With the aim to reduce the dimensionality of the problem, the ERA-40 fields have been decomposed using empirical orthogonal functions. Available daily data has been divided into two parts: a training period used to find a group of about 300 analogues to build the downscaling model (1961-1996) and a test period (19972001), where models' performance has been assessed using independent data. In the case of surface moisture flux, the models based on analogues followed by random forests do not clearly outperform those built on analogues plus multiple linear regression, while simple averages calculated from the nearest analogues found in the training period, yielded only slightly worse results. In the case of precipitation, the three types of model performed equally. These results suggest that most of the models' downscaling capabilities can be attributed to the analogues-calculation stage.