Role Of Vertical Component Of Surface Tension Of The Droplet On The Elastic Deformation Of Pdms Membrane

Poly(dimethylsiloxane) (PDMS) has been widely used in lab-on-a-chip and micro- total analysis systems (mu-TAS), thus wetting and electrowetting behaviors of PDMS are of great importance in these devices. PDMS is a kind of soft polymer material, so the elastic deformation of PDMS membrane by a droplet cannot be neglected due to the vertical component of the interfacial tension between the liquid and vapor, and this vertical component of liquid-vapor surface tension is also balanced by the stress distribution within the PDMS membrane. Such elastic deformation and stress distribution not only affect the exact measurement of contact angle, but also have influence on the micro-fluidic behavior of the devices. Using ANSYS code, we simulated numerically the elastic deformation and stress distribution of PDMS membrane on a rigid substrate due to the liquid-vapor surface tension. It is found that the vertical elastic deformation of the PDMS membrane is on the order of several tens of nanometers due to the application of a droplet with a diameter of 2.31 mm, which is no longer negligible for lab-on-a-chip and mu-TAS. The vertical elastic deformation increases with the thickness of the PDMS membrane, and there exists a saturated membrane thickness, regarded as a semi-infinite membrane thickness, and the vertical elastic deformation reaches a limiting value when the membrane thickness is equal to or thicker than such saturated thickness. (C) Koninklijke Brill NV, Leiden, 2008.

Impact of difference accuracy on computational properties of vertical grids for a nonhydrostatic model

Starting from nonhydrostatic Boussinesq approximation equations, a general method is introduced to deduce the dispersion relationships. A comparative investigation is performed on inertia-gravity wave with horizontal lengths of 100, 10 and 1 km. These are examined using the second-order central difference scheme and the fourth-order compact difference scheme on vertical grids that are currently available from the perspectives of frequency, horizontal and vertical component of group velocity. These findings are compared to analytical solutions. The obtained results suggest that whether for the second-order central difference scheme or for the fourth-order compact difference scheme, Charny-Phillips and Lorenz ( L) grids are suitable for studying waves at the above-mentioned horizontal scales; the Lorenz time-staggered and Charny-Phillips time staggered (CPTS) grids are applicable only to the horizontal scales of less than 10 km, and N grid ( unstaggered grid) is unsuitable for simulating waves at any horizontal scale. Furthermore, by using fourth-order compact difference scheme with higher difference precision, the errors of frequency and group velocity in horizontal and vertical directions produced on all vertical grids in describing the waves with horizontal lengths of 1, 10 and 100 km cannot inevitably be decreased. So in developing a numerical model, the higher-order finite difference scheme, like fourth-order compact difference scheme, should be avoided as much as possible, typically on L and CPTS grids, since it will not only take many efforts to design program but also make the calculated group velocity in horizontal and vertical directions even worse in accuracy.

Deformation Of Pdms Membrane And Microcantilever By A Water Droplet: Comparison Between Mooney-Rivlin And Linear Elastic Constitutive Models

In this paper, we studied the role of vertical component Of Surface tension of a water droplet on the deformation of membranes and microcantilevers (MCLs) widely used in lab-on-a-chip and micro-and nano-electromechanical system (MEMS/NEMS). Firstly, a membrane made of a rubber-like material, poly(dimethylsiloxane) (PDMS), was considered. The deformation was investigated using the Mooney-Rivlin (MR) model and the linear elastic constitutive relation, respectively. By comparison between the numerical solutions with two different models, we found that the simple linear elastic model is accurate enough to describe such kind of problem, which would be quite convenient for engineering applications. Furthermore, based on small-deflection beam theory, the effect of a liquid droplet on the deflection of a MCL was also studied. The free-end deflection of the MCL was investigated by considering different cases like a cylindrical droplet, a spherical droplet centered on the MCL and a spherical droplet arbitrarily positioned on the MCL. Numerical simulations demonstrated that the deflection might not be neglected, and showed good agreement with our theoretical analyses. (C) 2008 Elsevier Inc. All rights reserved.

Elastic deformation of soft membrane with finite thickness induced by a sessile liquid droplet

In this paper, the role of vertical component of Surface tension of a droplet on the elastic deformation of a finite-thickness flexible membrane was theoretically analyzed using Hankel transformation. The vertical displacement at the Surface was derived and can be reduced to Lester's or Rusanov's solutions when the thickness is infinite. Moreover, some Simulations of the effect of a liquid droplet on a membrane with a finite thickness were made. The numerical results showed that there exists a saturated membrane thickness of the order of millimeter, when the thickness of a membrane is larger than such a value, the membrane can be regarded as a half-infinite body. Further numerical calculations for soft membrane whose thickness is far below the saturated thickness were made. By comparison between the maximum vertical displacement of an ultrathin soft membrane and a half-infinite body, we found that Lester's or Rusanov's solutions for a half-infinite body cannot correctly describe Such cases. In other words, the thickness of a soft membrane has great effect on the surface deformation of the ultrathin membrane induced by a liquid droplet. (C) 2009 Elsevier Inc. All rights reserved.

Spatial structure of internal tidal waves generated on weak topographies

The generation of internal gravity waves by barotropic tidal flow passing over a two-dimensional topography is investigated. Rather than calculating the conversion of tidal energy, this study focuses on delineating the geometric characteristics of the spatial structure of the resulting internal wave fields (i.e., the configurations of the internal beams and their horizontal projections) which have usually been ignored. it is found that the various possible wave types can be demarcated by three characteristic frequencies: the tidal frequency, wo; the buoyancy frequency, N; and the vertical component of the Coriolis vector or earth's rotation.f. When different possibilities arising from the sequence of these frequencies are considered, there occur 12 kinds of wave structures in the full 3D space in contrast to the 5 kinds identified by the 2D theory. The constant wave phase lines may form as ellipses or hyperbolic lines on the horizontal plane, provided the buoyancy frequency is greater or less than the tidal frequency. The effect that stems from the consideration of the basic flow is also found, which not only serves as the reason for the occurrence of higtter harmonics but also increases the wave strength in the direction of basic flow. (C) 2009 Elsevier B.V. All rights reserved.

Millstone Hill电离层F层等效风场的气候学研究

Neutral winds and electric fields in the ionospheric F layer play important roles in the variations of the ionosphere, and also affect the thermospheric circulation via the close coupling between the ionosphere and the thermosphere. By now, the neutral winds and electric drifts are generally observed with ground-based Fabry-Perot interferometers (FPI) and incoherent scatter radars (ISR), rockets, and satellite-borne instrument. Based on the servo theory, the ionospheric equivalent winds, which include the information of both the neutral winds and electric fields, can be derived from these characteristic parameters observed by ionosondes. This indirect derivation has potential values in climatological researches and space weather forecast. With the data set of the incoherent scatter radar observations at Millstone Hill, USA, from 1976 to 2006, we statistically analyzed the climatological variations of the vertical component of the equivalent winds (VEWs) over Millstone Hill, which are derived from the ionospheric key parameters (the peak electron number density and peak height of the F2 layer, NmF2 and hmF2) on the basis of the servo theory, Liu's method, and measurements from the ion line-of-sight velocity as well. The main results of this analysis are summarized as follows: (1) The values of VEWs over Millstone Hill during nighttime are stronger than in the daytime, and the upward drift dominates most of the day. In 1993, Hagan found that the component of the neutral winds in the magnetic meridion in daytime is weaker than during nighttime under both solar maximum and minimum conditions; he also found that the equatorward winds dominate most of the day. Both results suggest that the thermosphere in Millstone Hill is modulated by the aurorally driven high-latitude circulation cell; that is, during geomagnetic quiet periods, the average auroral activity is strong enough to drive thermospheric circulation equatorward for most of the day at Millstone Hill. Moreover, since ion drag is the strongest during daytime when F region densities are enhanced by photoionization, the wind speeds are smaller during the daytime than in the nighttime. (2) There is equinoctial symmetry in VEWs at Millstone Hill. The amplitudes and phases of VEWs in spring are quite similar to those in autumn. In contrast, the nighttime upward drift in winter is weaker than in summer and the difference becomes more significant with increasing solar activity. This solstice asymmetry indicates that, the aurorally driven circulation in the northern hemisphere at Millstone Hill latitude is weaker in winter due to arctic darkness, because the subsolar point is in the southern hemisphere. (3) The comparison of the VEWs derived from three methods, i.e., the servo theory, Liu's method, and the ISR ion line-of-sight velocity measurements, indicates that the amplitudes and main phase tendencies of these VEWs accord well with each other during nighttime hours. However, the case in the daytime is relatively worse. This daytime discrepancy can be explained in terms of the effects of photochemical processes and the choices of the servo constants. A larger servo constant gives a stronger plasma drift in daytime. Therefore, this study tells how important to choose a suitable constant for deriving VEWs at Millstone Hill.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

On the mechanism of the formation of horizontal cracks in a vertical column of saturated sand

Extended horizontal cracks have! been observed experimentally in a vertical column of saturated sand when a flow of water is forced to percolate upward through it. This paper provides a theory for this phenomenon. It will be shown that the presence of inhomogeneity in permeability along the length of the column is essential for such cracks to develop. It will also be shown that small initial inhomogeneity may be magnified through the transport of the finer component of the sand by percolation. Under certain conditions liquefaction takes place at a section of the sand column causing a crack to initiate and grow there. This theory is found to be in good qualitative agreement with the experimental findings.

Visible vertical cavity surface emitting laser

We have designed and fabricated the visible vertical-cavity surface-emitting lasers (VCSEL's) by using metalorganic vapor phase epitaxy (MOVPE). We use the 8 lambda optical cavities with 3 quantum wells in AlGaInP/AlGaAs red VCSEL's to reduce the drift leakage current and enhance the model gain in AlGaInP active region. The structure has a p-type stack with 36 DBR pairs on the top and an n-type with 55-1/2 pairs on the bottom. Using micro-area reflectance spectrum, we try to get a better concordance between the center wavelength of DBR and the emitting wavelength of the active region. We used a component graded layer of 0.05 lambda thick (x = 0.5 similar to 0.9) at the p-type DBR AlGaAs/AlAs interface to reduce the resistance of p-type DBR. We use selective oxidation to define the current injection path. Because the oxidation rate of a thick layer is faster than a thinner one, we grown a thick AlAs layer close to the active region. In this way, we got a smaller active region for efficient confinement of injected carriers (the aperture area is 3 x 3 mu m) to reduce the threshold and, at the same time, a bigger conductive area in the DBR layers to reduce the resistance. We employ Zn doping on the p-side of the junction to improve hole injection and control the Zn dopant diffusion to get proper p-i-n junction. At room temperature, pulse operation of the laser has been achieved with the low threshold current of 0.8mA; the wavelength is about 670nm.

Biological application of multi-component nanowires in hybrid devices powered by F-1-ATPase motors

In this paper, construction of hybrid device by integrating nanowires with F1-ATPase motors is described. The nickel nanowires and multi-segment nanowires, including gold and nickel, were fabricated by electrochemical deposition in nanoporous templates. The nickel nanowires functionalized by biotinylated peptide can be assembled directly onto F1-ATPase motors to act as the propellers. If the multicomponent nanowires, including gold and nickel, were selectively functionalized by the thiol group modified ssDNA and the synthetic peptide, respectively, the biotinylated F1- ATPase motors can be attached to the biotinylated peptide on nickel segment of the nanowires. Then, the multi-component nanowires can also be used as the propellers, and one may observe the rotations of the multi-component nanowires driven by F1-ATPase motors. Therefore, introduction of multiple segments along the length of a nanowire can lead to a variety of multiple chemical functionalities, which can be selectively bound to cells and special biomolecules. This method provides an insight for the construction of other hybrid devices with its controlling arrangement of different biomolecule on designed nanometer scale structures.

Instability Analysis of Nonlinear Surface Waves in a Circular Cylindrical Container Subjected to a Vertical Excitation

Singular perturbation theory of two-time scale expansions was developed both in inviscid and weak viscous fluids to investigate the motion of single surface standing wave in a liquid-filled circular cylindrical vessel, which is subject to a vertical periodical oscillation. Firstly, it is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear evolution equation of slowly varying complex amplitude, which incorporates cubic nonlinear term, external excitation and the influence of surface tension, was derived from solvability condition of high-order approximation. It shows that when forced frequency is low, the effect of surface tension on mode selection of surface wave is not important. However, when forced frequency is high, the influence of surface tension is significant, and can not be neglected. This proved that the surface tension has the function, which causes free surface returning to equilibrium location. Theoretical results much close to experimental results when the surface tension is considered. In fact, the damping will appear in actual physical system due to dissipation of viscosity of fluid. Based upon weakly viscous fluids assumption, the fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates damping term and external excitation, was derived from linearized Navier-Stokes equation. The analytical expression of damping coefficient was determined and the relation between damping and other related parameters (such as viscosity, forced amplitude and depth of fluid) was presented. The nonlinear amplitude equation and a dispersion, which had been derived from the inviscid fluid approximation, were modified by adding linear damping. It was found that the modified results much reasonably close to experimental results. Moreover, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent. Finally, instability of the surface wave is analyzed and properties of the solutions of the modified amplitude equation are determined together with phase-plane trajectories. A necessary condition of forming stable surface wave is obtained and unstable regions are illustrated. (c) 2005 Elsevier SAS. All rights reserved.

Component Assembling Model and Its Application to Quasi-Brittle Damage

Potential energy can be approximated by ‘‘pair-functional’’ potentials which is composed of pair potentials and embedding energy. Pair potentials are grouped according to discrete directions of atomic bonds such that each group is represented by an orientational component. Meanwhile, another kind of component, the volumetric one is derived from embedding energy. Damage and fracture are the changing and breaking of atomic bonds at the most fundamental level and have been reflected by the changing of these components’ properties. Therefore, material is treated as a component assembly, and its constitutive equations are formed by means of assembling these two kinds of components’ response functions. This material model is referred to as the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of physical explicitness and intrinsic induced anisotropy, etc.

Frequency characteristics and dynamical behaviors of self-modulation in vertical-cavity surface-emitting lasers

The frequency characteristics of a VCSEL with a quarter-wave plate (QWP) and an external reflector are investigated with the translation matrix of the vectorial field. Two series of eigenmode with a shift of half the free spectrum range are linearly polarized, respectively, along the neutral axes of QWP. We also numerically explore the polarization self-modulation phenomenon by using a vectorial laser equation and considering the inhomogeneous broadening of the gain medium. If the external cavity is so short that the shift is bigger than the homogeneous broadening, two stable longitudinal modes oscillate, respectively, on the neutral axes of QWP because they consume different carriers. With a long external cavity, the competition of the modes for the common carriers causes the intensity fluctuation of the modes with a period of one round-trip time of the external cavity.

Characterization of Air-Water Two-Phase Vertical Flow by Using Electrical Resistance Imaging

The characterization of air-water two-phase vertical flow in a 12 m flow loop with 1.5 m of vertical section is studied by using electrical resistance tomography (ERT). By applying a fast data collection to a dual-plane ERT sensor and an iterative image reconstruction algorithm, relevant information is gathered for implementation of flow characteristics, particularly for flow regime recognition. A cross-correlation method is also used to interpret the velocity distribution of the gas phase on the cross section. The paper demonstrates that ERT can now be deployed routinely for velocity measurements and this capability will increase as faster measurement systems evolve.

Liquefaction and Displacement of Saturated Sand Under Vertical Vibration Loading

In order to investigate the influence of the vertical vibration loading on the liquefaction of saturated sand, one dimensional model for the saturated sand with a vertical vibration is presented based on the two phase continuous media theory. The development of the liquefaction and the liquefaction region are analyzed. It is shown that the vertical vibration loading could induce liquefaction. The rate of the liquefaction increases with the increase of the initial limit strain or initial porosity or amplitude and frequency of loading, and increases with the decrease of the permeability or initial modulus. It is shown also that there is a phase lag in the sand column. When the sand permeability distribution is non-uniform, the pore pressure and the strain will rise sharply where the permeability is the smallest, and fracture might be induced. With the development of liquefaction, the strength of the soil foundation becomes smaller and smaller. In the limiting case, landslides or debris flows could occur.