Solvation dynamics in a Brownian dipolar lattice. Comparison between computer simulation and various molecular theories of solvation dynamics

Several recent theoretical and computer simulation studies have considered solvation dynamics in a Brownian dipolar lattice which provides a simple model solvent for which detailed calculations can be carried out. In this article a fully microscopic calculation of the solvation dynamics of an ion in a Brownian dipolar lattice is presented. The calculation is based on the non‐Markovian molecular hydrodynamic theory developed recently. The main assumption of the present calculation is that the two‐particle orientational correlation functions of the solid can be replaced by those of the liquid state. It is shown that such a calculation provides an excellent agreement with the computer simulation results. More importantly, the present calculations clearly demonstrate that the frequency‐dependent dielectric friction plays an important role in the long time decay of the solvation time correlation function. We also find that the present calculation provides somewhat better agreement than either the dynamic mean spherical approximation (DMSA) or the Fried–Mukamel theory which use the simulated frequency‐dependent dielectric function. It is found that the dissipative kernels used in the molecular hydrodynamic approach and in the Fried–Mukamel theory are vastly different, especially at short times. However, in spite of this disagreement, the two theories still lead to comparable results in good agreement with computer simulation, which suggests that even a semiquantitatively accurate dissipative kernel may be sufficient to obtain a reliable solvation time correlation function. A new wave vector and frequency‐dependent dissipative kernel (or memory function) is proposed which correctly goes over to the appropriate expressions in both the single particle and the collective limits. This form is expected to lead to better results than all the existing descriptions.

A fuzzy logic based dynamic wave model inversion algorithm for canal regulation

A fuzzy logic based centralized control algorithm for irrigation canals is presented. Purpose of the algorithm is to control downstream discharge and water level of pools in the canal, by adjusting discharge release from the upstream end and gates settings. The algorithm is based on the dynamic wave model (Saint-Venant equations) inversion in space, wherein the momentum equation is replaced by a fuzzy rule based model, while retaining the continuity equation in its complete form. The fuzzy rule based model is developed on fuzzification of a new mathematical model for wave velocity, the derivational details of which are given. The advantages of the fuzzy control algorithm, over other conventional control algorithms, are described. It is transparent and intuitive, and no linearizations of the governing equations are involved. Timing of the algorithm and method of computation are explained. It is shown that the tuning is easy and the computations are straightforward. The algorithm provides stable, realistic and robust outputs. The disadvantage of the algorithm is reduced precision in its outputs due to the approximation inherent in the fuzzy logic. Feed back control logic is adopted to eliminate error caused by the system disturbances as well as error caused by the reduced precision in the outputs. The algorithm is tested by applying it to water level control problem in a fictitious canal with a single pool and also in a real canal with a series of pools. It is found that results obtained from the algorithm are comparable to those obtained from conventional control algorithms.

New light on the optical equivalence theorem and a new type of discrete diagonal coherent state representation

In many instances we find it advantageous to display a quantum optical density matrix as a generalized statistical ensemble of coherent wave fields. The weight functions involved in these constructions turn out to belong to a family of distributions, not always smooth functions. In this paper we investigate this question anew and show how it is related to the problem of expanding an arbitrary state in terms of an overcomplete subfamily of the overcomplete set of coherent states. This provides a relatively transparent derivation of the optical equivalence theorem. An interesting by-product is the discovery of a new class of discrete diagonal representations.

Modeling and analysis of a tunable piezoelectric structure for transverse shear wave generation

An analytical investigation of the transverse shear wave mode tuning with a resonator mass (packing mass) on a Lead Zirconium Titanate (PZT) crystal bonded together with a host plate and its equivalent electric circuit parameters are presented. The energy transfer into the structure for this type of wave modes are much higher in this new design. The novelty of the approach here is the tuning of a single wave mode in the thickness direction using a resonator mass. First, a one-dimensional constitutive model assuming the strain induced only in the thickness direction is considered. As the input voltage is applied to the PZT crystal in the thickness direction, the transverse normal stress distribution induced into the plate is assumed to have parabolic distribution, which is presumed as a function of the geometries of the PZT crystal, packing mass, substrate and the wave penetration depth of the generated wave. For the PZT crystal, the harmonic wave guide solution is assumed for the mechanical displacement and electric fields, while for the packing mass, the former is solved using the boundary conditions. The electromechanical characteristics in terms of the stress transfer, mechanical impedance, electrical displacement, velocity and electric field are analyzed. The analytical solutions for the aforementioned entities are presented on the basis of varying the thickness of the PZT crystal and the packing mass. The results show that for a 25% increase in the thickness of the PZT crystal, there is ~38% decrease in the first resonant frequency, while for the same change in the thickness of the packing mass, the decrease in the resonant frequency is observed as ~35%. Most importantly the tuning of the generated wave can be accomplished with the packing mass at lower frequencies easily. To the end, an equivalent electric circuit, for tuning the transverse shear wave mode is analyzed.

A new impulse peak voltmeter

The paper describes a simple instrument for the measurement of the peak amplitudes of impulse voltages up to 250 v with an accuracy to ±3%. The response of the instrument is fast enough to read the peak amplitude of a 0.5/10 μsec impulse wave and its response remains the same for impulses of longer duration. Its favourable response has been obtained by charging a capacitor through a thyratron and measuring the voltage across it by an inverted triode voltmeter. The discharge time constant of the instrument is 5000 sec so that the reading can be taken at leisure. It can be used for the measurement of peak amplitudes of repetitive impulse and power frequency voltages also

New type of high field electrooptic response in nematics

We have made concurrent measurements of ionic current and optical transmission between crossed polarisers on several nematics with positive dielectric anisotropy under the action of applied low frequency (< 1KHz) square wave voltages. When the field E is low, the measured current is linear in E and there is no electrooptic response. Beyond some value of the field (E(0)similar to 100 esu), the current becomes independent of the field (phenomenon of limiting current). Further an electrooptic signal is measured at twice the frequency of the applied voltage, which exhibits a peak as a function of the field. The width of the peak is 3 to 4 times the value of E-0, and the signal level at the peak decreases as the frequency is increased. These measurements have been made on three highly polar compounds with cyano end groups. Careful observations do not show any evidence of electrohydrodynamic instabilities in the sample. It is argued that the observations can be understood if at the onset of the phenomenon of the limiting current, a strong electric field gradient is established near one of the electrodes due to the sweeping of an ionic species with high mobility. The field gradient produces a flexoelectric deformation of the director field, which in turn gives rise to the electrooptic effect. At higher fields, the stabilising dielectric torque takes over to suppress this instability.

Rapid localization of damage using a circular sensor array and Lamb wave based triangulation

A damage detection and imaging methodology based on symmetry of neighborhood sensor path and similarity of signal patterns with respect to radial paths in a circular array of sensors has been developed It uses information regarding Limb wave propagation along with a triangulation scheme to rapidly locate and quantify the severity of damage without using all of the sensor data. In a plate like structure, such a scheme can be effectively employed besides full field imaging of wave scattering pattern from the damage, if present in the plate. This new scheme is validated experimentally. Hole and corrosion type damages have been detected and quantified using the proposed scheme successfully. A wavelet based cumulative damage index has been studied which shows monotonic sensitivity against the severity of the damage. which is most desired in a Structural Health Monitoring system. (C) 2010 Elsevier Ltd. All rights reserved.

Topological defects in crystals: A density-wave theory

A new approach for describing dislocations and other topological defects in crystals, based on the density wave theory of Ramakrishnan and Yussouff is presented. Quantitative calculations are discussed in brief for the order parameter profiles, the atomic configuration and the free energy of a screw dislocation with Burgers vector b = (a/2, a/2,a/2 ) in a bcc solid. Our results for the free energy of the dislocation in a crystal of sizeR, when expressed as (λb 2/4π) ln (αR/|b|) whereλ is the shear elastic constant, yield, for example, the valueα ⋍ 1·85 for sodium at its freezing temperature (371°K). The density distribution in the presence of the dislocation shows that the dislocation core has a columnar character. To our knowledge, this study represents the first calculation of dislocation structure, including the core, within the framework of an order parameter theory incorporating thermal effects.

Microscopic expression for frequency and wave vector dependent dielectric constant of a dipolar liquid

A microscopic expression for the frequency and wave vector dependent dielectric constant of a dense dipolar liquid is derived starting from the linear response theory. The new expression properly takes into account the effects of the translational modes in the polarization relaxation. The longitudinal and the transverse components of the dielectric constant show vastly different behavior at the intermediate values of the wave vector k. We find that the microscopic structure of the dense liquid plays an important role at intermediate wave vectors. The continuum model description of the dielectric constant, although appropriate at very small values of wave vector, breaks down completely at the intermediate values of k. Numerical results for the longitudinal and the transverse dielectric constants are obtained by using the direct correlation function from the mean‐spherical approximation for dipolar hard spheres. We show that our results are consistent with all the limiting expressions known for the dielectric function of matter.

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.

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.

PVU and wave-particle splitting schemes for Euler equations of gas dynamics

A new way of flux-splitting, termed as the wave-particle splitting is presented for developing upwind methods for solving Euler equations of gas dynamics. Based on this splitting, two new upwind methods termed as Acoustic Flux Vector Splitting (AFVS) and Acoustic Flux Difference Splitting (AFDS) methods are developed. A new Boltzmann scheme, which closely resembles the wave-particle splitting, is developed using the kinetic theory of gases. This method, termed as Peculiar Velocity based Upwind (PVU) method, uses the concept of peculiar velocity for upwinding. A special feature of all these methods that the unidirectional and multidirectional parts of the flux vector are treated separately. Extensive computations done using these schemes demonstrate the soundness of the ideas.

A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

Diaphragmless shock wave generators for industrial applications of shock waves

The prime focus of this study is to design a 50 mm internal diameter diaphragmless shock tube that can be used in an industrial facility for repeated loading of shock waves. The instantaneous rise in pressure and temperature of a medium can be used in a variety of industrial applications. We designed, fabricated and tested three different shock wave generators of which one system employs a highly elastic rubber membrane and the other systems use a fast acting pneumatic valve instead of conventional metal diaphragms. The valve opening speed is obtained with the help of a high speed camera. For shock generation systems with a pneumatic cylinder, it ranges from 0.325 to 1.15 m/s while it is around 8.3 m/s for the rubber membrane. Experiments are conducted using the three diaphragmless systems and the results obtained are analyzed carefully to obtain a relation between the opening speed of the valve and the amount of gas that is actually utilized in the generation of the shock wave for each system. The rubber membrane is not suitable for industrial applications because it needs to be replaced regularly and cannot withstand high driver pressures. The maximum shock Mach number obtained using the new diaphragmless system that uses the pneumatic valve is 2.125 +/- 0.2%. This system shows much promise for automation in an industrial environment.

A new beam finite element for the analysis of functionally graded materials

A new beam element is developed to study the thermoelastic behavior of functionally graded beam structures. The element is based on the first-order shear deformation theory and it accounts for varying elastic and thermal properties along its thickness. The exact solution of static part of the governing differential equations is used to construct interpolating polynomials for the element formulation. Consequently, the stiffness matrix has super-convergent property and the element is free of shear locking. Both exponential and power-law variations of material property distribution are used to examine different stress variations. Static, free vibration and wave propagation problems are considered to highlight the behavioral difference of functionally graded material beam with pure metal or pure ceramic beams. (C) 2003 Elsevier Science Ltd. All rights reserved.