Elasto-Electrostatic Analysis of Circular Microplates Used in Capacitive Micromachined Ultrasonic Transducers

The active structural component of a capacitive micromachined ultrasonic transducer (CMUT) is the top plate which vibrates under the influence of a time-varying electrostatic force thereby producing ultrasound waves of the desired frequency in the surrounding medium. Analysis of MEMS devices which rely on electrostatic actuation is complicated due to the fact that the structural deformations alter the electrostatic forces, which redistribute and modify the applied loads. Hence, it becomes imperative to consider the electrostatics-structure coupling aspect in the design of these devices. This paper presents an approximate analytical solution for the static deflection of a thin, clamped circular plate caused by electrostatic forces which are inherently nonlinear. Traditionally, finite element simulations using some commercial software such as ANSYS are employed to determine the structural deflections caused by electrostatic forces. Since the structural deformation alters the electrostatic field, a coupled-field simulation is required wherein the electrostatic mesh is continuously updated to coincide with the deflection of the structure. Such simulations are extremely time consuming, in addition to being nontransparent and somewhat hard to implement. We employ the classical thin-plate theory which is adequate when the ratio of the diameter to thickness of the plate is very large, a situation commonly prevalent in many MEMS devices, especially the CMUTs. We solve the thin-plate electrostatic-elastic equation using the Galerkin-weighted residual technique, under the assumption that the deflections are small in comparison to the thickness of the plate. The evaluation of the electrostatic force between the two plates is simplified due to the fact that the electrostatic gap is much smaller than the lateral dimensions of the device. The results obtained are compared to those found from ANSYS simulations and an excellent agreement is observed between the two. The pull-in voltage predicted by our model is close to the value predicted by ANSYS simulations.

Convergence feedback and unstable low frequency oscillations in a simple coupled ocean‐atmosphere model

The role of convergence feedback on the stability of a coupled ocean‐atmosphere system is studied using model III of Hirst (1986). It is shown that the unstable coupled mode found by Hirst is greatly modified by the convergence feedback. If the convergence feedback strength exceeds a critical value, several new unstable intraseasonal modes are also introduced. These modes have very weak dependence on the wave number. These results may explain the behaviour of some coupled models and to some extent provide a mechanism for the observed aperiodicity of the El‐Nino and Southern Oscillation (ENSO) events.

Structure and representation of correlation functions and the density matrix for a statistical wave field in optics

A systematic structure analysis of the correlation functions of statistical quantum optics is carried out. From a suitably defined auxiliary two‐point function we are able to identify the excited modes in the wave field. The relative simplicity of the higher order correlation functions emerge as a byproduct and the conditions under which these are made pure are derived. These results depend in a crucial manner on the notion of coherence indices and of unimodular coherence indices. A new class of approximate expressions for the density operator of a statistical wave field is worked out based on discrete characteristic sets. These are even more economical than the diagonal coherent state representations. An appreciation of the subtleties of quantum theory obtains. Certain implications for the physics of light beams are cited.

A real transformation for modal analysis of flexible damped spacecraft

Modal approach is widely used for the analysis of dynamics of flexible structures. However, space analysts yet lack an intimate modal analysis of current spacecraft which are rich with flexibility and possess both structural and discrete damping. Mathematical modeling of such spacecraft incapacitates the existing real transformation procedure, for it cannot include discrete damping, demands uncomputable inversion of a modal matrix inaccessible due to its overwhelming size and does not permit truncation. On the other hand, complex transformation techniques entail more computational time and cannot handle structural damping. This paper presents a real transformation strategy which averts inversion of the associated real transformation matrix, allows truncation and accommodates both forms of damping simultaneously. This is accomplished by establishing a key relation between the real transformation matrix and its adjoint. The relation permits truncation of the matrices and leads to uncoupled pairs of coupled first order equations which contain a number of adjoint eigenvectors. Finally these pairs are solved to obtain a literal modal response of forced gyroscopic damped flexibile systems at arbitrary initial conditions.

A mathematical analysis of nonlinear waves in a fluid filled visco-elastic tube

Our investigations in this paper are centred around the mathematical analysis of a ldquomodal waverdquo problem. We have considered the axisymmetric flow of an inviscid liquid in a thinwalled viscoelastic tube under certain simplifying assumptions. We have first derived the propagation space equations in the long wave limit and also given a general procedure to derive these equations for arbitrary wave length, when the flow is irrotational. We have used the method of operators of multiple scales to derive the nonlinear Schrödinger equation governing the modulation of periodic waves and we have elaborated on the ldquolong modulated wavesrdquo and the ldquomodulated long wavesrdquo. We have also examined the existence and stability of Stokes waves in this system. This is followed by a discussion of the progressive wave solutions of the long wave equations. One of the most important results of our paper is that the propagation space equations are no longer partial differential equations but they are in terms of pseudo-differential operators.Die vorliegenden Untersuchungen beziehen sich auf die mathematische Behandlung des ldquorModalwellenrdquo-Problems. Die achsensymmetrische Strömung einer nichtviskosen Flüssigkeit in einem dünnwandigen viskoelastischen Rohr, unter bestimmten vereinfachenden Annahmen, wird betrachtet. Zuerst werden die Gleichungen des Ausbreitungsraumes im Langwellenbereich abgeleitet und eine allgemeine Methode zur Herleitung dieser Gleichungen für beliebige Wellenlängen bei nichtrotierender Strömung angegeben. Eine Operatorenmethode mit multiplem Maßstab wird verwendet zur Herleitung der nichtlinearen Schrödinger-Gleichung für die Modulation der periodischen Wellen, und die ldquorlangmodulierten Wellenrdquo sowie die ldquormodulierten Langwellenrdquo werden aufgezeigt. Weiters wird die Existenz und die Stabilität der Stokes-Wellen im System untersucht. Anschließend werden die progressiven Wellenlösungen der Langwellengleichungen diskutiert. Eines der wichtigsten Ergebnisse dieser Arbeit ist, daß die Gleichungen des Ausbreitungsraumes keine partiellen Differentialgleichungen mehr sind, sondern Ausdrücke von Pseudo-Differentialoperatoren.

Spectral element approach to wave propagation in uncertain beam structures

This paper presents a study on the uncertainty in material parameters of wave propagation responses in metallic beam structures. Special effort is made to quantify the effect of uncertainty in the wave propagation responses at high frequencies. Both the modulus of elasticity and the density are considered uncertain. The analysis is performed using a Monte Carlo simulation (MCS) under the spectral finite element method (SEM). The randomness in the material properties is characterized by three different distributions, the normal, Weibull and extreme value distributions. Their effect on wave propagation in beams is investigated. The numerical study shows that the CPU time taken for MCS under SEM is about 48 times less than for MCS under a conventional one-dimensional finite element environment for 50 kHz loading. The numerical results presented investigate effects of material uncertainties on high frequency modes. A study is performed on the usage of different beam theories and their uncertain responses due to dynamic impulse load. These studies show that even for a small coefficient of variation, significant changes in the above parameters are noticed. A number of interesting results are presented, showing the true effects of uncertainty response due to dynamic impulse load.

Development of a New Finite Element for the Analysis of Sandwich Beams with Soft Core

A new super convergent sandwich beam finite element formulation is presented in this article. This element is a two-nodded, six degrees of freedom (dof) per node (3 dof u(0), w, phi for top and bottom face sheets each), which assumes that all the axial and flexural loads are taken by face sheets, while the core takes only the shear loads. The beam element is formulated based on first-order shear deformation theory for the face sheets and the core displacements are assumed to vary linearly across the thickness. A number of numerical experiments involving static, free vibration, and wave propagation analysis examples are solved with an aim to show the super convergent property of the formulated element. The examples presented in this article consider both metallic and composite face sheets. The formulated element is verified in most cases with the results available in the published literature.

Analysis of an Aquifer-water table aquitard system

The problem of pumping an aquifer in an aquifer-water table aquitard system is considered, accounting for the elastic properties of both the aquifer and the aquitard, the gravity drainage in the aquitard and treating the water table as an unknown boundary. The coupled partial differential equations are nondimensionalised, yielding three principal parameters governing the problem. The numerical solution of these equations is obtained for a wide range of parameter values. Type curves are generated and their use is illustrated through a field application.

Analysis of Nonlinear Optical Properties of Photonic Crystal Beam-splitters

The paper analyses electromagnetic wave propagation through nonlinear photonic crystal beam-splitters. Different lattice configurations of Y-junction beam-splitters are simulated and propagation properties are investigated with introducing nonlinearity with varying the rod size in crystal lattice. It is seen that nonlinear photonic crystal shows a considerable band-gap even at low refractive contrast. The division of power in both arms of beam-splitters can be controlled by varying the nonlinearity.

Electrostatic Field Analysis of Electrooptic Devices

We describe a Finite Difference Method for the determination of the electrostatic field in a multilayered electrooptic device. The Laplace equation is solved, assuming a suitable closed area, by taking into account the different permittivities of the various layers. The effect of a higher permittivity in the guiding layer has been explicitly considered. As a practical example, we calculate the phase shift of a guided optical wave within an electrooptic modulator. A review of the various methods in use for the field analysis is given. Some criteria for the selection of the appropriate method are also mentioned.

Mechanisms of variability and predictability of the tropical coupled ocean-atmosphere system

A conceptual model is proposed to explain the observed aperiodicity in the short term climate fluctuations of the tropical coupled ocean-atmosphere system. This is based on the evidence presented here that the tropical coupled ocean-atmosphere system sustains a low frequency inter-annual mode and a host of higher frequency intra-seasonal unstable modes. At long wavelengths, the low frequency mode is dominant while at short wavelengths, the high frequency modes are dominant resulting in the co-existence of a long wave low frequency mode with some short wave intra-seasonal modes in the tropical coupled system. It is argued that due to its long wavelength, the low frequency mode would behave like a linear oscillator while the higher frequency short wave modes would be nonlinear. The conceptual model envisages that an interaction between the low frequency linear oscillator and the high frequency nonlinear oscillations results in the observed aperiodicity of the tropical coupled system. This is illustrated by representing the higher frequency intra-seasonal oscillations by a nonlinear low order model which is then coupled to a linear oscillator with a periodicity of four years. The physical mechanism resulting in the aperiodicity in the low frequency oscillations and implications of these results on the predictability of the coupled system are discussed.

Surface acoustic wave distribution and acousto-optic interaction in proton exchanged LiNbO3 waveguides

The efficiency of acoustooptic (AO) interaction in YZ-cut proton exchanged (PE) LiNbO3 waveguides is theoretically analysed by determining the overlap between the optical and acoustic field distributions. The present analysis takes into account the perturbed SAW field distribution due to the presence of the PE layer on the LiNbO3 substrate determined by the rigorous layered medium approach. The overlap is found to be significant upto very high acoustic frequencies of the order of 5 GHz, whereas in the earlier analysis by vonHelmolt and Schaffer [6] for diffused waveguides, it was shown that the overlap integral rolls down to nearly zero at this high frequency range.

Strong nonlocalization induced by small scale parameter on terahertz flexural wave dispersion characteristics of a monolayer graphene

This paper presents the strong nonlocal scale effect on the flexural wave propagation in a monolayer graphene sheet. The graphene is modeled as an isotropic plate of one atom thick. Nonlocal governing equation of motion is derived and wave propagation analysis is performed using spectral analysis. The present analysis shows that the flexural wave dispersion in graphene obtained by local and nonlocal elasticity theories is quite different. The nonlocal elasticity calculation shows that the wavenumber escapes to infinite at certain frequency and the corresponding wave velocity tends to zero at that frequency indicating localization and stationary behavior. This behavior is captured in the spectrum and dispersion curves. The cut-off frequency of flexural wave not only depend on the axial wavenumber but also on the nonlocal scaling parameter. The effect of axial wavenumber on the wave behavior in graphene is also discussed in the present manuscript. (C) 2010 Elsevier B.V. All rights reserved.

Effect of open joint on stress wave propagation

Stress wave characteristics are drastically altered by joints and other inhomogenities. This paper addresses the effect of an open joint on stress wave transmission. An elastodynamic analysis is developed to supplement and explain some recent observations by Fourney and Dick(1995) on open as well as filled joints. The analytical model developed here assuming spherical symmetry can be extended to filled joints between dissimilar media, but results are presented only for open joints separating identical materials. As a special case, stress wave transmission across a joint with no gap is also addressed.

Novel geodesic splitting on paraboloid of revolution with applications to ray-theoretic analysis

Computer simulations have shown a novel geodesic splitting on the paraboloid of revolution leading to a multiplicity of surface ray paths. Such a phenomenon would have wide ramifications for wave propagation problems in general, besides applications in target-detection problems and the computational requirements of ray-theoretic formulations such as the UTD, in computing the antenna characteristics in the high-frequency domain.