Wave propagation analysis in laminated composite plates with transverse cracks using the wavelet spectral finite element method

This paper presents a newly developed wavelet spectral finite element (WFSE) model to analyze wave propagation in anisotropic composite laminate with a transverse surface crack penetrating part-through the thickness. The WSFE formulation of the composite laminate, which is based on the first-order shear deformation theory, produces accurate and computationally efficient results for high frequency wave motion. Transverse crack is modeled in wavenumber-frequency domain by introducing bending flexibility of the plate along crack edge. Results for tone burst and impulse excitations show excellent agreement with conventional finite element analysis in Abaqus (R). Problems with multiple cracks are modeled by assembling a number of spectral elements with cracks in frequency-wavenumber domain. Results show partial reflection of the excited wave due to crack at time instances consistent with crack locations. (C) 2014 Elsevier B.V. All rights reserved.

Atom Lithography

We review the two kinds of forces that near-resonant light exerts on atoms the spontaneous force that is used for laser cooling, and the stimulated force that is used for coherent manipulation of atoms. We will discuss an experiment where laser cooling is used to collimate an atomic beam of sodium atoms, and the stimulated force within one period of a one-dimensional standing wave is used as a lens to focus the atoms to a narrow line about 20 nm wide. This kind of atom lithography is an example of the general field of atom optics in which light is used to manipulate atoms.

A spectral multiscale method for wave propagation analysis: Atomistic-continuum coupled simulation

In this paper, we present a new multiscale method which is capable of coupling atomistic and continuum domains for high frequency wave propagation analysis. The problem of non-physical wave reflection, which occurs due to the change in system description across the interface between two scales, can be satisfactorily overcome by the proposed method. We propose an efficient spectral domain decomposition of the total fine scale displacement along with a potent macroscale equation in the Laplace domain to eliminate the spurious interfacial reflection. We use Laplace transform based spectral finite element method to model the macroscale, which provides the optimum approximations for required dynamic responses of the outer atoms of the simulated microscale region very accurately. This new method shows excellent agreement between the proposed multiscale model and the full molecular dynamics (MD) results. Numerical experiments of wave propagation in a 1D harmonic lattice, a 1D lattice with Lennard-Jones potential, a 2D square Bravais lattice, and a 2D triangular lattice with microcrack demonstrate the accuracy and the robustness of the method. In addition, under certain conditions, this method can simulate complex dynamics of crystalline solids involving different spatial and/or temporal scales with sufficient accuracy and efficiency. (C) 2014 Elsevier B.V. All rights reserved.

Nonclassical Paths in Quantum Interference Experiments

In a double slit interference experiment, the wave function at the screen with both slits open is not exactly equal to the sum of the wave functions with the slits individually open one at a time. The three scenarios represent three different boundary conditions and as such, the superposition principle should not be applicable. However, most well-known text books in quantum mechanics implicitly and/or explicitly use this assumption that is only approximately true. In our present study, we have used the Feynman path integral formalism to quantify contributions from nonclassical paths in quantum interference experiments that provide a measurable deviation from a naive application of the superposition principle. A direct experimental demonstration for the existence of these nonclassical paths is difficult to present. We find that contributions from such paths can be significant and we propose simple three-slit interference experiments to directly confirm their existence.

Dynamic Site Characterization and Correlation of Shear Wave Velocity with Standard Penetration Test `N' Values for the City of Agartala, Tripura State, India

Seismic site characterization is the basic requirement for seismic microzonation and site response studies of an area. Site characterization helps to gauge the average dynamic properties of soil deposits and thus helps to evaluate the surface level response. This paper presents a seismic site characterization of Agartala city, the capital of Tripura state, in the northeast of India. Seismically, Agartala city is situated in the Bengal Basin zone which is classified as a highly active seismic zone, assigned by Indian seismic code BIS-1893, Indian Standard Criteria for Earthquake Resistant Design of Structures, Part-1 General Provisions and Buildings. According to the Bureau of Indian Standards, New Delhi (2002), it is the highest seismic level (zone-V) in the country. The city is very close to the Sylhet fault (Bangladesh) where two major earthquakes (M (w) > 7) have occurred in the past and affected severely this city and the whole of northeast India. In order to perform site response evaluation, a series of geophysical tests at 27 locations were conducted using the multichannel analysis of surface waves (MASW) technique, which is an advanced method for obtaining shear wave velocity (V (s)) profiles from in situ measurements. Similarly, standard penetration test (SPT-N) bore log data sets have been obtained from the Urban Development Department, Govt. of Tripura. In the collected data sets, out of 50 bore logs, 27 were selected which are close to the MASW test locations and used for further study. Both the data sets (V (s) profiles with depth and SPT-N bore log profiles) have been used to calculate the average shear wave velocity (V (s)30) and average SPT-N values for the upper 30 m depth of the subsurface soil profiles. These were used for site classification of the study area recommended by the National Earthquake Hazard Reduction Program (NEHRP) manual. The average V (s)30 and SPT-N classified the study area as seismic site class D and E categories, indicating that the city is susceptible to site effects and liquefaction. Further, the different data set combinations between V (s) and SPT-N (corrected and uncorrected) values have been used to develop site-specific correlation equations by statistical regression, as `V (s)' is a function of SPT-N value (corrected and uncorrected), considered with or without depth. However, after considering the data set pairs, a probabilistic approach has also been presented to develop a correlation using a quantile-quantile (Q-Q) plot. A comparison has also been made with the well known published correlations (for all soils) available in the literature. The present correlations closely agree with the other equations, but, comparatively, the correlation of shear wave velocity with the variation of depth and uncorrected SPT-N values provides a more suitable predicting model. Also the Q-Q plot agrees with all the other equations. In the absence of in situ measurements, the present correlations could be used to measure V (s) profiles of the study area for site response studies.

Comparative study of different nonconserving time integrators for wave propagation in hyperelastic waveguides

This paper addresses the formulation and numerical efficiency of various numerical models of different nonconserving time integrators for studying wave propagation in nonlinear hyperelastic waveguides. The study includes different nonlinear finite element formulations based on standard Galerkin finite element model, time domain spectral finite element model, Taylor-Galerkin finite element model, generalized Galerkin finite element model and frequency domain spectral finite element model. A comparative study on the computational efficiency of these different models is made using a hyperelastic rod model, and the optimal computational scheme is identified. The identified scheme is then used to study the propagation of transverse and longitudinal waves in a Timoshenko beam with Murnaghan material nonlinearity.

Efficient Adaptive Mesh Refinement for MoM-based Package-Board 3D Full-wave Extraction

The accurate solution of 3D full-wave Method of Moments (MoM) on an arbitrary mesh of a package-board structure does not guarantee accuracy, since the discretizations may not be fine enough to capture rapid spatial changes in the solution variable. At the same time, uniform over-meshing on the entire structure generates large number of solution variables and therefore requires an unnecessarily large matrix solution. In this work, a suitable refinement criterion for MoM based electromagnetic package-board extraction is proposed and the advantages of the adaptive strategy are demonstrated from both accuracy and speed perspectives.

Spiral-wave dynamics in ionically realistic mathematical models for human ventricular tissue: the effects of periodic deformation

We carry out an extensive numerical study of the dynamics of spiral waves of electrical activation, in the presence of periodic deformation (PD) in two-dimensional simulation domains, in the biophysically realistic mathematical models of human ventricular tissue due to (a) ten-Tusscher and Panfilov (the TP06 model) and (b) ten-Tusscher, Noble, Noble, and Panfilov (the TNNPO4 model). We first consider simulations in cable-type domains, in which we calculate the conduction velocity theta and the wavelength lambda of a plane wave; we show that PD leads to a periodic, spatial modulation of theta and a temporally periodic modulation of lambda; both these modulations depend on the amplitude and frequency of the PD. We then examine three types of initial conditions for both TP06 and TNNPO4 models and show that the imposition of PD leads to a rich variety of spatiotemporal patterns in the transmembrane potential including states with a single rotating spiral (RS) wave, a spiral-turbulence (ST) state with a single meandering spiral, an ST state with multiple broken spirals, and a state SA in which all spirals are absorbed at the boundaries of our simulation domain. We find, for both TP06 and TNNPO4 models, that spiral-wave dynamics depends sensitively on the amplitude and frequency of PD and the initial condition. We examine how these different types of spiral-wave states can be eliminated in the presence of PD by the application of low-amplitude pulses by square- and rectangular-mesh suppression techniques. We suggest specific experiments that can test the results of our simulations.

Response of shock wave deformation in AA5086 aluminum alloy

The AA5086 aluminum alloy sheets with different starting textures were subjected to shock wave deformation with an input impulse of similar to 0.2 Ns. Microstructural examination indicate no significant change in grain size; however, the evolution of substructure manifesting intra-granular misorientation was evident. The improvement in hardness indicates the absence of recovery and strain hardening during shock deformation. Shock deformed samples show characteristic texture evolution with high Brass {110}< 112 > component. The study demonstrates the viability of high velocity forming of AA5086 aluminum alloy sheet using shock wave. (C) 2014 Elsevier B.V. All rights reserved.

Wave propagation analysis in adhesively bonded composite joints using the wavelet spectral finite element method

A wavelet spectral finite element (WSFE) model is developed for studying transient dynamics and wave propagation in adhesively bonded composite joints. The adherands are formulated as shear deformable beams using the first order shear deformation theory (FSDT) to obtain accurate results for high frequency wave propagation. Equations of motion governing wave motion in the bonded beams are derived using Hamilton's principle. The adhesive layer is modeled as a line of continuously distributed tension/compression and shear springs. Daubechies compactly supported wavelet scaling functions are used to transform the governing partial differential equations from time domain to frequency domain. The dynamic stiffness matrix is derived under the spectral finite element framework relating the nodal forces and displacements in the transformed frequency domain. Time domain results for wave propagation in a lap joint are validated with conventional finite element simulations using Abaqus. Frequency domain spectrum and dispersion relation results are presented and discussed. The developed WSFE model yields efficient and accurate analysis of wave propagation in adhesively-bonded composite joints. (C) 2014 Elsevier Ltd. All rights reserved.

Remotely triggered micro-shock wave responsive drug delivery system for resolving diabetic wound infection and controlling blood sugar levels

A novel, micro-shock wave responsive spermidine and dextran sulfate microparticle was developed. Almost 90% of the drug release was observed when the particles were exposed to micro-shock waves 5 times. Micro-shock waves served two purposes; of releasing the antibiotic from the system and perhaps disrupting the S. aureus biofilm in the skin infection model. A combination of shock waves with ciprofloxacin loaded microparticles could completely cure the S. aureus infection lesion in a diabetic mouse model. As a proof of concept insulin release was triggered using micro-shock waves in diabetic mice to reduce the blood glucose level. Insulin release could be triggered for at least 3 days by exposing subcutaneously injected insulin loaded particles.

Wave propagation of phonon and phason displacement modes in quasicrystals: Determination of wave parameters

The paper presents the study of wave propagation in quasicrystals. Our interest is in the computation of the wavenumber (k(n)) and group speed (c(g)) of the phonon and phason displacement modes of one, two, and three dimensional quasicrystals. These wave parameter expressions are derived and computed using the elasto-hydrodynamic equations for quasicrystals. For the computation of the wavenumber and group speeds, we use Fourier transform approximation of the phonon and the phason displacement modes. The characteristic equations obtained are a polynomial equation of the wavenumber (k(n)), with frequency as a parameter. The corresponding group speeds (c(g)) for different frequencies are then computed from the wavenumber k(n). The variation of wavenumber and group speeds with frequency is plotted for the 1-D quasicrystal, 2-D decagonal Al-Ni-Co quasicrystals, and 3-D icosahedral Al-Pd-Mn and Zn-Mg-Sc quasicrystals. From the wavenumber and group speeds plots, we obtain the cut-off frequencies for different spatial wavenumber eta(m). The results show that for 1-D, 2-D, and 3-D quasicrystals, the phonon displacement modes are non-dispersive for low values of eta(m) and becomes dispersive for increasing values of eta(m). The cut-off frequencies are not observed for very low values of eta(m), whereas the cut-off frequency starts to appear with increasing eta(m). The group speeds of the phason displacement modes are orders of magnitude lower than that of the phonon displacement modes, showing that the phason modes do not propagate, and they are essentially the diffusive modes. The group speeds of the phason modes are also not influenced by eta(m). The group speeds for the 2-D quasicrystal at 35 kHz is also simulated numerically using Galerkin spectral finite element methods in frequency domain and is compared with the results obtained using wave propagation analysis. The effect of the phonon and phason elastic constants on the group speeds is studied using 3-D icosahedral Al-Pd-Mn and Zn-Mg-Sc quasicrystals. It is also shown that the phason elastic constants and the coupling coefficient do not affect the group speeds of the phonon displacement modes. (C) 2015 AIP Publishing LLC.

Optimal Pulse Width Modulation of Voltage-Source Inverter fed Motor Drives with Relaxation of Quarter Wave Symmetry Condition

Optimal switching angles for minimization of total harmonic distortion of line current (I-THD) in a voltage source inverter are determined traditionally by imposing half-wave symmetry (HWS) and quarter-wave symmetry (QWS) conditions on the pulse width modulated waveform. This paper investigates optimal switching angles with QWS relaxed. Relaxing QWS expands the solution space and presents the possibility of improved solutions. The optimal solutions without QWS are shown here to outperform the optimal solutions with QWS over a range of modulation index (M) between 0.82 and 0.94 for a switching frequency to fundamental frequency ratio of 5. Theoretical and experimental results are presented on a 2.3kW induction motor drive.

L-p estimates for the wave equation associated to the Grushin operator

We prove that the solution of the wave equation associated to the Grushin operator G = -Delta -vertical bar x vertical bar(2)partial derivative(2)(t) is bounded on L-P (Rn+1), with 1 < p < infinity, when vertical bar 1/p - 1/2 vertical bar < 1/n+2.

Anisotropic isometric fluctuation relations in experiment and theory on a self-propelled rod

The isometric fluctuation relation (IFR) P. I. Hurtado et al., Proc. Natl. Acad. Sci. USA 108, 7704 (2011)] relates the relative probability of current fluctuations of fixed magnitude in different spatial directions. We test its validity in an experiment on a tapered rod, rendered motile by vertical vibration and immersed in a sea of spherical beads. We analyze the statistics of the velocity vector of the rod and show that they depart significantly from the IFR of Hurtado et al. Aided by a Langevin-equation model we show that our measurements are largely described by an anisotropic generalization of the IFR R. Villavicencio et al., Europhys. Lett. 105, 30009 (2014)], with no fitting parameters, but with a discrepancy in the prefactor whose origin may lie in the detailed statistics of the microscopic noise. The experimentally determined large-deviation function of the velocity vector has a kink on a curve in the plane.