Time and Frequency Domain Finite Element Models for Axial Wave Analysis in Hyperelastic Rods

The propagation of axial waves in hyperelastic rods is studied using both time and frequency domain finite element models. The nonlinearity is introduced using the Murnaghan strain energy function and the equations governing the dynamics of the rod are derived assuming linear kinematics. In the time domain, the standard Galerkin finite element method, spectral element method, and Taylor-Galerkin finite element method are considered. A frequency domain formulation based on the Fourier spectral method is also developed. It is found that the time domain spectral element method provides the most efficient numerical tool for the problem considered.

GINZBURG-LANDAU LIKE THEORY OF HIGH TEMPERATURE SUPERCONDUCTIVITY IN THE CUPRATES: EMERGENT d-WAVE ORDER

High temperature superconductivity in the cuprates remains one of the most widely investigated, constantly surprising and poorly understood phenomena in physics. Here, we describe briefly a new phenomenological theory inspired by the celebrated description of superconductivity due to Ginzburg and Landau and believed to describe its essence. This posits a free energy functional for the superconductor in terms of a complex order parameter characterizing it. We propose that there is, for superconducting cuprates, a similar functional of the complex, in plane, nearest neighbor spin singlet bond (or Cooper) pair amplitude psi(ij). Further, we suggest that a crucial part of it is a (short range) positive interaction between nearest neighbor bond pairs, of strength J'. Such an interaction leads to nonzero long wavelength phase stiffness or superconductive long range order, with the observed d-wave symmetry, below a temperature T-c similar to zJ' where z is the number of nearest neighbors; d-wave superconductivity is thus an emergent, collective consequence. Using the functional, we calculate a large range of properties, e. g., the pseudogap transition temperature T* as a function of hole doping x, the transition curve T-c(x), the superfluid stiffness rho(s)(x, T), the specific heat (without and with a magnetic field) due to the fluctuating pair degrees of freedom and the zero temperature vortex structure. We find remarkable agreement with experiment. We also calculate the self-energy of electrons hopping on the square cuprate lattice and coupled to electrons of nearly opposite momenta via inevitable long wavelength Cooper pair fluctuations formed of these electrons. The ensuing results for electron spectral density are successfully compared with recent experimental results for angle resolved photo emission spectroscopy (ARPES), and comprehensively explain strange features such as temperature dependent Fermi arcs above T-c and the ``bending'' of the superconducting gap below T-c.

Automatic energy-momentum conserving time integrators for hyperelastic waves

An energy-momentum conserving time integrator coupled with an automatic finite element algorithm is developed to study longitudinal wave propagation in hyperelastic layers. The Murnaghan strain energy function is used to model material nonlinearity and full geometric nonlinearity is considered. An automatic assembly algorithm using algorithmic differentiation is developed within a discrete Hamiltonian framework to directly formulate the finite element matrices without recourse to an explicit derivation of their algebraic form or the governing equations. The algorithm is illustrated with applications to longitudinal wave propagation in a thin hyperelastic layer modeled with a two-mode kinematic model. Solution obtained using a standard nonlinear finite element model with Newmark time stepping is provided for comparison. (C) 2012 Elsevier B.V. All rights reserved.

Development of micro-shock wave assisted dry particle and fluid jet delivery system

Small quantity of energetic material coated on the inner wall of a polymer tube is proposed as a new method to generate micro-shock waves in the laboratory. These micro-shock waves have been harnessed to develop a novel method of delivering dry particle and liquid jet into the target. We have generated micro-shock waves with the help of reactive explosive compound high melting explosive (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) and traces of aluminium] coated polymer tube, utilising 9 J of energy. The detonation process is initiated electrically from one end of the tube, while the micro-shock wave followed by the products of detonation escape from the open end of the polymer tube. The energy available at the open end of the polymer tube is used to accelerate tungsten micro-particles coated on the other side of the diaphragm or force a liquid jet out of a small cavity filled with the liquid. The micro-particles deposited on a thin metal diaphragm (typically 100-mu m thick) were accelerated to high velocity using micro-shock waves to penetrate the target. Tungsten particles of 0.7 mu m diameter have been successfully delivered into agarose gel targets of various strengths (0.6-1.0 %). The device has been tested by delivering micro-particles into potato tuber and Arachis hypogaea Linnaeus (ground nut) stem tissue. Along similar lines, liquid jets of diameter 200-250 mu m (methylene blue, water and oils) have been successfully delivered into agarose gel targets of various strengths. Successful vaccination against murine salmonellosis was demonstrated as a biological application of this device. The penetration depths achieved in the experimental targets are very encouraging to develop a future device for biological and biomedical applications.

Nonlocal continuum mechanics based ultrasonic flexural wave dispersion characteristics of a monolayer graphene embedded in polymer matrix

Wave propagation in graphene sheet embedded in elastic medium (polymer matrix) has been a topic of great interest in nanomechanics of graphene sheets, where the equivalent continuum models are widely used. In this manuscript, we examined this issue by incorporating the nonlocal theory into the classical plate model. The influence of the nonlocal scale effects has been investigated in detail. The results are qualitatively different from those obtained based on the local/classical plate theory and thus, are important for the development of monolayer graphene-based nanodevices. In the present work, the graphene sheet is modeled as an isotropic plate of one-atom thick. The chemical bonds are assumed to be formed between the graphene sheet and the elastic medium. The polymer matrix is described by a Pasternak foundation model, which accounts for both normal pressure and the transverse shear deformation of the surrounding elastic medium. When the shear effects are neglected, the model reduces to Winkler foundation model. The normal pressure or Winkler elastic foundation parameter is approximated as a series of closely spaced, mutually independent, vertical linear elastic springs where the foundation modulus is assumed equivalent to stiffness of the springs. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. An ultrasonic type of flexural wave propagation model is also derived and the results of the wave dispersion analysis are shown for both local and nonlocal elasticity calculations. From this analysis we show that the elastic matrix highly affects the flexural wave mode and it rapidly increases the frequency band gap of flexural mode. The flexural wavenumbers obtained from nonlocal elasticity calculations are higher than the local elasticity calculations. The corresponding wave group speeds are smaller in nonlocal calculation as compared to local elasticity calculation. The effect of y-directional wavenumber (eta(q)) on the spectrum and dispersion relations of the graphene embedded in polymer matrix is also observed. We also show that the cut-off frequencies of flexural wave mode depends not only on the y-direction wavenumber but also on nonlocal scaling parameter (e(0)a). The effect of eta(q) and e(0)a on the cut-off frequency variation is also captured for the cases of with and without elastic matrix effect. For a given nanostructure, nonlocal small scale coefficient can be obtained by matching the results from molecular dynamics (MD) simulations and the nonlocal elasticity calculations. At that value of the nonlocal scale coefficient, the waves will propagate in the nanostructure at that cut-off frequency. In the present paper, different values of e(0)a are used. One can get the exact e(0)a for a given graphene sheet by matching the MD simulation results of graphene with the results presented in this article. (c) 2012 Elsevier Ltd. All rights reserved.

Study of terahertz wave propagation properties in nanoplates with surface and small-scale effects

In macroscopic and even microscopic structural elements, surface effects can be neglected and classical theories are sufficient. As the structural size decreases towards the nanoscale regime, the surface-to-bulk energy ratio increases and surface effects must be taken into account. In the present work, the terahertz wave dispersion characteristics of a nanoplate are studied with consideration of the surface effects as well as the nonlocal small-scale effects. Nonlocal elasticity theory of plate is used to derive the general differential equation based on equilibrium approach to include those scale effects. Scale and surface property dependent wave characteristic equations are obtained via spectral analysis. For the present study the material properties of an anodic alumina with crystallographic of < 111 > direction are considered. The present analysis shows that the effect of surface properties on the flexural waves of nanoplates is more significant. It can be found that the flexural wavenumbers with surface effects are high as compared to that without surface effects. The scale effects show that the wavenumbers of the flexural wave become highly non-linear and tend to infinite at certain frequency. After that frequency the wave will not propagate and the corresponding wave velocities tend to zero at that frequency (escape frequency). The effects of surface stresses on the wave propagation properties of nanoplate are also captured in the present work. (C) 2012 Elsevier Ltd. All rights reserved.

SCALE EFFECTS ON ULTRASONIC WAVE DISPERSION CHARACTERISTICS OF MONOLAYER GRAPHENE EMBEDDED IN AN ELASTIC MEDIUM

Ultrasonic wave propagation in a graphene sheet, which is embedded in an elastic medium, is studied using nonlocal elasticity theory incorporating small-scale effects. The graphene sheet is modeled as an one-atom thick isotropic plate and the elastic medium/substrate is modeled as distributed springs. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. After that, an ultrasonic type of wave propagation model is also derived. The explicit expressions for the cut-off frequencies are also obtained as functions of the nonlocal scaling parameter and the y-directional wavenumber. Local elasticity shows that the wave will propagate even at higher frequencies. But nonlocal elasticity predicts that the waves can propagate only up to certain frequencies (called escape frequencies), after which the wave velocity becomes zero. The results also show that the escape frequencies are purely a function of the nonlocal scaling parameter. The effect of the elastic medium is captured in the wave dispersion analysis and this analysis is explained with respect to both local and nonlocal elasticity. The simulations show that the elastic medium affects only the flexural wave mode in the graphene sheet. The presence of the elastic matrix increases the band gap of the flexural mode. The present results can provide useful guidance for the design of next-generation nanodevices in which graphene-based composites act as a major element.

Investigations on micro-blast wave assisted metal foil forming for biomedical applications

The deformation dynamics of metal foils (<0.25 mm thick) subjected to micro-blast wave are presented in this paper. The energy of micro-blast wave emanating from the open end of a polymer tube is used to deliver micro-particles for bio-medical applications. In these experiments metal foils are used to transfer the energy of the micro-blast wave to the micro-particles. Using cubic root scaling law the over pressure of the blast wave at the open end of the polymer tube is estimated and using this peak plate over pressure is estimated. The finite element analysis is used to estimate the velocity profile of the deforming metal foils. The finite element analysis results are compared with experimental results for the maximum deformation and deformed shape. Based on the deformation velocity, metal foil to be used for experiments is selected. Among the materials investigated 0.1 mm thick brass foil has the maximum velocity of 205 m/s and is used in the experiments. It is found from finite element analysis that the particles deposited within a radius of 0.5 mm will leave the foil with nearly equal velocity (error < 5%). The spray cone angle which is the angle of deviation of the path of particles from the axis of the polymer tube is also estimated and found to be less than 7 degrees up to a radius of 0.75 mm. Illustrative experiments are carried out to deliver micro particles (0.7 mu m diameter tungsten) into plant tissues. Particle penetration depth up to 460 mu m was achieved in ground tissue of potato tuber. (C) 2012 Elsevier Ltd. All rights reserved.

Efficient energy transport in photosynthesis: roles of coherence and entanglement

Recently it has been discovered---contrary to expectations of physicists as well as biologists---that the energy transport during photosynthesis, from the chlorophyll pigment that captures the photon to the reaction centre where glucose is synthesised from carbon dioxide and water, is highly coherent even at ambient temperature and in the cellular environment. This process and the key molecular ingredients that it depends on are described. By looking at the process from the computer science view-point, we can study what has been optimised and how. A spatial search algorithmic model based on robust features of wave dynamics is presented.

Noniterative inversion strategy for photoacoustic tomography: recovery of absorbed energy map from boundary pressure measurements

The inverse problem in photoacoustic tomography (PAT) seeks to obtain the absorbed energy map from the boundary pressure measurements for which computationally intensive iterative algorithms exist. The computational challenge is heightened when the reconstruction is done using boundary data split into its frequency spectrum to improve source localization and conditioning of the inverse problem. The key idea of this work is to modify the update equation wherein the Jacobian and the perturbation in data are summed over all wave numbers, k, and inverted only once to recover the absorbed energy map. This leads to a considerable reduction in the overall computation time. The results obtained using simulated data, demonstrates the efficiency of the proposed scheme without compromising the accuracy of reconstruction.

Reconstructing the properties of dark energy using standard sirens

Future space-based gravity wave (GW) experiments such as the Big Bang Observatory (BBO), with their excellent projected, one sigma angular resolution, will measure the luminosity distance to a large number of GW sources to high precision, and the redshift of the single galaxies in the narrow solid angles towards the sources will provide the redshifts of the gravity wave sources. One sigma BBO beams contain the actual source in only 68% of the cases; the beams that do not contain the source may contain a spurious single galaxy, leading to misidentification. To increase the probability of the source falling within the beam, larger beams have to be considered, decreasing the chances of finding single galaxies in the beams. Saini et al. T.D. Saini, S.K. Sethi, and V. Sahni, Phys. Rev. D 81, 103009 (2010)] argued, largely analytically, that identifying even a small number of GW source galaxies furnishes a rough distance-redshift relation, which could be used to further resolve sources that have multiple objects in the angular beam. In this work we further develop this idea by introducing a self-calibrating iterative scheme which works in conjunction with Monte Carlo simulations to determine the luminosity distance to GW sources with progressively greater accuracy. This iterative scheme allows one to determine the equation of state of dark energy to within an accuracy of a few percent for a gravity wave experiment possessing a beam width an order of magnitude larger than BBO (and therefore having a far poorer angular resolution). This is achieved with no prior information about the nature of dark energy from other data sets such as type Ia supernovae, baryon acoustic oscillations, cosmic microwave background, etc. DOI:10.1103/PhysRevD.87.083001

Influence of crack parameters on wavelet energy correlated damage index

This paper illustrates a Wavelet Coefficient based approach using experiments to understand the sensitivity of ultrasonic signals due to parametric variation of a crack configuration in a metal plate. A PZT patch sensor/actuator system integrated to a metal plate with through-thickness crack is used. The proposed approach uses piezoelectric patches, which can be used to both actuate and sense the ultrasonic signals. While this approach leads to more flexibility and reduced cost for larger scalability of the sensor/actuator network, the complexity of the signals increases as compared to what is encountered in conventional ultrasonic NDE problems using selective wave modes. A Damage Index (DI) has been introduced, which is function of wavelet coefficient. Experiments have been carried out for various crack sizes, crack orientations and band-limited tone-burst signal through FIR filter. For a 1 cm long crack interrogated with 20 kHz tone-burst signal, the Damage Index (DI) for the horizontal crack orientation increases by about 70% with respect to that for 135 degrees oriented crack and it increases by about 33% with respect to the vertically oriented crack. The detailed results reported in this paper is a step forward to developing computational schemes for parametric identification of damage using sensor/actuator network and ultrasonic wave.

Light transmission through and its complete stoppage in an ultra slow wave optical medium

Light wave transmission - its compression, amplification, and the optical energy storage in an ultra slow wave medium (USWM) is studied analytically. Our phenomenological treatment is based entirely on the continuity equation for the optical energy flux, and the well-known distribution-product property of Dirac delta-function. The results so obtained provide a clear understanding of some recent experiments on light transmission and its complete stoppage in an USWM.

Energy loss due to adhesion in longitudinal impact of elastic cylinders

Adhesive interaction between impacting bodies can cause energy loss, even in an otherwise elastic impact. Adhesion force induces tensile stress in the bodies, which modifies the stress wave profile and influences the restitution behavior. We investigate this effect by developing a finite element framework, which incorporates a Lennard-Jones-type potential for modeling the adhesive interaction between volume elements. With this framework, the classical problems in contact mechanics can be revisited without the restrictive surface-force approximation. In this paper, we study the longitudinal impact of an elastic cylinder on a rigid half-space with adhesion. In the absence of adhesion, this problem reduces to the impact between two identical cylinders in which there is no energy loss. Adhesion causes a fraction of energy in the stress waves to remain in the cylinder as residual stress waves. This apparent loss in kinetic energy is shown to be a unique function of maximum tensile strain energy. We have developed a 1-D model in terms of interaction force parameters, velocity and material properties to estimate the tensile stain energy. We show that this model can be used to predict practically important phenomena like capture wherein the impacting bodies stick together. (C) 2013 Elsevier Masson SAS. All rights reserved.

Nonlinear gravity-wave interactions in stratified turbulence

To investigate the dynamics of gravity waves in stratified Boussinesq flows, a model is derived that consists of all three-gravity-wave-mode interactions (the GGG model), excluding interactions involving the vortical mode. The GGG model is a natural extension of weak turbulence theory that accounts for exact three-gravity-wave resonances. The model is examined numerically by means of random, large-scale, high-frequency forcing. An immediate observation is a robust growth of the so-called vertically sheared horizontal flow (VSHF). In addition, there is a forward transfer of energy and equilibration of the nonzero-frequency (sometimes called ``fast'') gravity-wave modes. These results show that gravity-wave-mode interactions by themselves are capable of systematic interscale energy transfer in a stratified fluid. Comparing numerical simulations of the GGG model and the full Boussinesq system, for the range of Froude numbers (Fr) considered (0.05 a parts per thousand currency sign Fr a parts per thousand currency sign 1), in both systems the VSHF is hardest to resolve. When adequately resolved, VSHF growth is more vigorous in the GGG model. Furthermore, a VSHF is observed to form in milder stratification scenarios in the GGG model than the full Boussinesq system. Finally, fully three-dimensional nonzero-frequency gravity-wave modes equilibrate in both systems and their scaling with vertical wavenumber follows similar power-laws. The slopes of the power-laws obtained depend on Fr and approach -2 (from above) at Fr = 0.05, which is the strongest stratification that can be properly resolved with our computational resources.