Inverse problem for the wave equation : partial data and novel boundary sources

An inverse problem for the wave equation is a mathematical formulation of the problem to convert measurements of sound waves to information about the wave speed governing the propagation of the waves. This doctoral thesis extends the theory on the inverse problems for the wave equation in cases with partial measurement data and also considers detection of discontinuous interfaces in the wave speed. A possible application of the theory is obstetric sonography in which ultrasound measurements are transformed into an image of the fetus in its mother's uterus. The wave speed inside the body can not be directly observed but sound waves can be produced outside the body and their echoes from the body can be recorded. The present work contains five research articles. In the first and the fifth articles we show that it is possible to determine the wave speed uniquely by using far apart sound sources and receivers. This extends a previously known result which requires the sound waves to be produced and recorded in the same place. Our result is motivated by a possible application to reflection seismology which seeks to create an image of the Earth s crust from recording of echoes stimulated for example by explosions. For this purpose, the receivers can not typically lie near the powerful sound sources. In the second article we present a sound source that allows us to recover many essential features of the wave speed from the echo produced by the source. Moreover, these features are known to determine the wave speed under certain geometric assumptions. Previously known results permitted the same features to be recovered only by sequential measurement of echoes produced by multiple different sources. The reduced number of measurements could increase the number possible applications of acoustic probing. In the third and fourth articles we develop an acoustic probing method to locate discontinuous interfaces in the wave speed. These interfaces typically correspond to interfaces between different materials and their locations are of interest in many applications. There are many previous approaches to this problem but none of them exploits sound sources varying freely in time. Our use of more variable sources could allow more robust implementation of the probing.

Homonuclear NOE in spatially periodic spin systems

The nuclear Overhauser effect equations are solved analytically for a homonuclear group of spins whose sites are periodically arranged, including the special cases where the spins lie at the vertices of a regular polygon and on a one-dimensional lattice. t is shown that, for long correlation times, the equations governing magnetization transfer resemble a diffusion equation. Furthermore the deviation from exact diffusion is quantitatively related to the molecular tumbling correlation time. Equations are derived for the range of magnetization travel subsequent to the perturbation of a single spin in a lattice for both the case of strictly dipolar relaxation and the more general situation where additional T1 mechanisms may be active. The theory given places no restrictions on the delay (or mixing) times, and it includes all the spins in the system. Simulations are presented to confirm the theory.

Magnetoresistance of FexNi80-xCr20 (50

The authors have measured longitudinal and transverse magnetoresistance (MR) of crystalline pseudo-binary alloys FexNi80-xCr20 (50

Meridional propagation of large-scale monsoon convective zones

Observational studies indicate that the convective activity of the monsoon systems undergo intraseasonal variations with multi-week time scales. The zone of maximum monsoon convection exhibits substantial transient behavior with successive propagating from the North Indian Ocean to the heated continent. Over South Asia the zone achieves its maximum intensity. These propagations may extend over 3000 km in latitude and perhaps twice the distance in longitude and remain as coherent entities for periods greater than 2-3 weeks. Attempts to explain this phenomena using simple ocean-atmosphere models of the monsoon system had concluded that the interactive ground hydrology so modifies the total heating of the atmosphere that a steady state solution is not possible, thus promoting lateral propagation. That is, the ground hydrology forces the total heating of the atmosphere and the vertical velocity to be slightly out of phase, causing a migration of the convection towards the region of maximum heating. Whereas the lateral scale of the variations produced by the Webster (1983) model were essentially correct, they occurred at twice the frequency of the observed events and were formed near the coastal margin, rather than over the ocean. Webster's (1983) model used to pose the theories was deficient in a number of aspects. Particularly, both the ground moisture content and the thermal inertia of the model were severely underestimated. At the same time, the sea surface temperatures produced by the model between the equator and the model's land-sea boundary were far too cool. Both the atmosphere and the ocean model were modified to include a better hydrological cycle and ocean structure. The convective events produced by the modified model possessed the observed frequency and were generated well south of the coastline. The improved simulation of monsoon variability allowed the hydrological cycle feedback to be generalized. It was found that monsoon variability was constrained to lie within the bounds of a positive gradient of a convective intensity potential (I). The function depends primarily on the surface temperature, the availability of moisture and the stability of the lower atmosphere which varies very slowly on the time scale of months. The oscillations of the monsoon perturb the mean convective intensity potential causing local enhancements of the gradient. These perturbations are caused by the hydrological feedbacks, discussed above, or by the modification of the air-sea fluxes caused by variations of the low level wind during convective events. The final result is the slow northward propagation of convection within an even slower convective regime. The ECMWF analyses show very similar behavior of the convective intensity potential. Although it is considered premature to use the model to conduct simulations of the African monsoon system, the ECMWF analysis indicates similar behavior in the convective intensity potential suggesting, at least, that the same processes control the low frequency structure of the African monsoon. The implications of the hypotheses on numerical weather prediction of monsoon phenomenon are discussed.

Analogues of the Wiener-Tauberian and Schwartz theorems for radial functions on symmetric spaces

We prove a Wiener Tauberian theorem for the L-1 spherical functions on a semisimple Lie group of arbitrary real rank. We also establish a Schwartz-type theorem for complex groups. As a corollary we obtain a Wiener Tauberian type result for compactly supported distributions.

E. C. G. Sudarshan and Symmetry in Classical Dynamics, Optics and Quantum Mechanics

We review work initiated and inspired by Sudarshan in relativistic dynamics, beam optics, partial coherence theory, Wigner distribution methods, multimode quantum optical squeezing, and geometric phases. The 1963 No Interaction Theorem using Dirac's instant form and particle World Line Conditions is recalled. Later attempts to overcome this result exploiting constrained Hamiltonian theory, reformulation of the World Line Conditions and extending Dirac's formalism, are reviewed. Dirac's front form leads to a formulation of Fourier Optics for the Maxwell field, determining the actions of First Order Systems (corresponding to matrices of Sp(2,R) and Sp(4,R)) on polarization in a consistent manner. These groups also help characterize properties and propagation of partially coherent Gaussian Schell Model beams, leading to invariant quality parameters and the new Twist phase. The higher dimensional groups Sp(2n,R) appear in the theory of Wigner distributions and in quantum optics. Elegant criteria for a Gaussian phase space function to be a Wigner distribution, expressions for multimode uncertainty principles and squeezing are described. In geometric phase theory we highlight the use of invariance properties that lead to a kinematical formulation and the important role of Bargmann invariants. Special features of these phases arising from unitary Lie group representations, and a new formulation based on the idea of Null Phase Curves, are presented.

Robust Semi-Blind Estimation for Beamforming Based MIMO Wireless Communication

In this paper, we present robust semi-blind (SB) algorithms for the estimation of beamforming vectors for multiple-input multiple-output wireless communication. The transmitted symbol block is assumed to comprise of a known sequence of training (pilot) symbols followed by information bearing blind (unknown) data symbols. Analytical expressions are derived for the robust SB estimators of the MIMO receive and transmit beamforming vectors. These robust SB estimators employ a preliminary estimate obtained from the pilot symbol sequence and leverage the second-order statistical information from the blind data symbols. We employ the theory of Lagrangian duality to derive the robust estimate of the receive beamforming vector by maximizing an inner product, while constraining the channel estimate to lie in a confidence sphere centered at the initial pilot estimate. Two different schemes are then proposed for computing the robust estimate of the MIMO transmit beamforming vector. Simulation results presented in the end illustrate the superior performance of the robust SB estimators.

High-Rate Space-Time Coded Large MIMO Systems: Low-Complexity Detection and Performance

Large MIMO systems with tens of antennas in each communication terminal using full-rate non-orthogonal space-time block codes (STBC) from Cyclic Division Algebras (CDA) can achieve the benefits of both transmit diversity as well as high spectral efficiencies. Maximum-likelihood (ML) or near-ML decoding of these large-sized STBCs at low complexities, however, has been a challenge. In this paper, we establish that near-ML decoding of these large STBCs is possible at practically affordable low complexities. We show that the likelihood ascent search (LAS) detector, reported earlier by us for V-BLAST, is able to achieve near-ML uncoded BER performance in decoding a 32x32 STBC from CDA, which employs 32 transmit antennas and sends 32(2) = 1024 complex data symbols in 32 time slots in one STBC matrix (i.e., 32 data symbols sent per channel use). In terms of coded BER, with a 16x16 STBC, rate-3/4 turbo code and 4-QAM (i.e., 24 bps/Hz), the LAS detector performs close to within just about 4 dB from the theoretical MIMO capacity. Our results further show that, with LAS detection, information lossless (ILL) STBCs perform almost as good as full-diversity ILL (FD-ILL) STBCs. Such low-complexity detectors can potentially enable implementation of high spectral efficiency large MIMO systems that could be considered in wireless standards.

Low-Complexity Detection and Performance in Multi-Gigabit High Spectral Efficiency Wireless Systems

Recently, we reported a low-complexity likelihood ascent search (LAS) detection algorithm for large MIMO systems with several tens of antennas that can achieve high spectral efficiencies of the order of tens to hundreds of bps/Hz. Through simulations, we showed that this algorithm achieves increasingly near SISO AWGN performance for increasing number of antennas in Lid. Rayleigh fading. However, no bit error performance analysis of the algorithm was reported. In this paper, we extend our work on this low-complexity large MIMO detector in two directions: i) We report an asymptotic bit error probability analysis of the LAS algorithm in the large system limit, where N-t, N-r -> infinity keeping N-t = N-r, where N-t and N-r are the number of transmit and receive antennas, respectively. Specifically, we prove that the error performance of the LAS detector for V-BLAST with 4-QAM in i.i.d. Rayleigh fading converges to that of the maximum-likelihood (ML) detector as N-t, N-r -> infinity keeping N-t = N-r ii) We present simulated BER and nearness to capacity results for V-BLAST as well as high-rate non-orthogonal STBC from Division Algebras (DA), in a more realistic spatially correlated MIMO channel model. Our simulation results show that a) at an uncoded BER of 10(-3), the performance of the LAS detector in decoding 16 x 16 STBC from DA with N-t = = 16 and 16-QAM degrades in spatially correlated fading by about 7 dB compared to that in i.i.d. fading, and 19) with a rate-3/4 outer turbo code and 48 bps/Hz spectral efficiency, the performance degrades by about 6 dB at a coded BER of 10(-4). Our results further show that providing asymmetry in number of antennas such that N-r > N-t keeping the total receiver array length same as that for N-r = N-t, the detector is able to pick up the extra receive diversity thereby significantly improving the BER performance.

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.

SiO Maser Emission and the Intrinsic Properties of Mira Variables

This paper presents observations of SiO maser emission from 161 Mira variables distributed over a wide range of intrinsic parameters like spectral type, bolometric magnitude and amplitude of pulsation. The observations were made at 86.243 GHz, using the 10.4 m millimeter-wave telescope of the Raman Research Institute at Bangalore, India. These are the first observations made using this telescope. From these observations, we have established that the maser emission is restricted to Miras having mean spectral types between M6 and M10. The infrared period-luminosity relation for Mira variables is used to calculate their distances and hence estimate their maser luminosities from the observed fluxes. The maser luminosity is found to be correlated with the bolometric magnitude of the Mira variable. On an H-R diagram, the masing Mira variables are shown to lie in a region distinct from that for the non-masing ones.

Variational principles in evolution

For a one-locus selection model, Svirezhev introduced an integral variational principle by defining a Lagrangian which remained stationary on the trajectory followed by the population undergoing selection. It is shown here (i) that this principle can be extended to multiple loci in some simple cases and (ii) that the Lagrangian is defined by a straightforward generalization of the one-locus case, but (iii) that in two-locus or more general models there is no straightforward extension of this principle if linkage and epistasis are present. The population trajectories can be constructed as trajectories of steepest ascent in a Riemannian metric space. A general method is formulated to find the metric tensor and the surface-in the metric space on which the trajectories, which characterize the variations in the gene structure of the population, lie. The local optimality principle holds good in such a space. In the special case when all possible linkage disequilibria are zero, the phase point of the n-locus genetic system moves on the surface of the product space of n higher dimensional unit spheres in a certain Riemannian metric space of gene frequencies so that the rate of change of mean fitness is maximum along the trajectory. In the two-locus case the corresponding surface is a hyper-torus.

Structure, Spectra and Redox Behaviour of Copper(11) Complexes of Bis( benzimidazoly1)diamine Ligands

The linear quadridentate ligand N,N'-bis(benzimidazoI-2-ylethyl)ethane-l,2-diamine (L') and its 1 - methylbenzimidazole analogue (L2) and homologues form 1 : 1 complexes with Cu(CIO,),; L' also forms complexes of the types CuL'X, where X = NO,, PF,, Br or CI and CuL'(X)Y where X = CI or Br and Y = CIO, or Br. Deep blue CuL1Br,*2H20 crystallizes in the monoclinic space group C2/c with Z = 4, a = 9.91 9(2), b = 16.626(3), c = 14.1 02(3) le\ and p = 94.39(2)". The structure was solved by Patterson and Fourier difference methods and refined by the least-squares technique to R = 0.064 for 2195 independent reflections with / > 1.50(/). The molecule lies on a two-fold axis symmetrically around Cu". The co-ordination around Cu" is found to be square planar with two amino nitrogens and two benzimidazole nitrogens forming the equatorial plane [CU-N 1.983(3) and 2.037(4) A]. The bromides are at longer distances [3.349(1) A] in axial sites. Ligand field and EPR spectra indicate that one bromide or chloride ion is axially co-ordinated to Cu" in [CuL1l2+. This ion exhibits quasi-reversible redox behaviour. Electrochemical studies of the dihalides in methanol have established the presence of [CuL'X,], [CuL'(X)]+ and [CuL'I2+ in equilibrium. In complexes with 565 [CuL4I2+ [L4 = N,Nbis( benzimidazol-2-ylmethyl)ethane-l,2-diamine] and 555 [CuL3] [L3 = N,N'-bis(1 -methylbenzimidazol- 2-ylmethyl)propane-l,3-diamine] chelate rings, Cull does not seem to lie in the N, square plane, as revealed by their low A values and irreversible electrochemical behaviour. The Cu"-Cu' redox potentials in methanol are in the order [CuL1I2+ < [CuL3I2+ < [CuL4I2+; this illustrates that sixmembered chelate rings are suitable to stabilize Cu", when CU-N 0 interactions are favourable.

Vibration and buckling of hybrid laminated curved panels

Vibration and buckling of curved plates, made of hybrid laminated composite materials, are studied using first-order shear deformation theory and Reissner's shallow shell theory. For an initial study, only simply-supported boundary conditions are considered. The natural frequencies and critical buckling loads are calculated using the energy method (Lagrangian approach) by assuming a combination of sine and cosine functions in the form of double Fourier series. The effects of curvature, aspect ratio, stacking sequence and ply-orientation are studied. The non-dimensional frequencies and critical buckling load of a hybrid laminate lie in between the values for laminates made of all plies of higher strength and lower strength fibres. Curvature enhances natural frequencies and it is more predominant for a thin panel than a thick one.

Plate-shaped transformation products in zirconium-base alloys

Plate-shaped products resulting from martensitic, diffusional, and mixed mode transformations in zirconium-base alloys are compared. in the present study. These alloys are particularly suitable for the comparison in view of the fact that the lattice correspondence between the parent beta (bcc) and the product alpha (hcp) or gamma-hydride (fct) phases are remarkably similar for different types of transformations. Crystallographic features such as orientation relations, habit planes, and interface structures associated with these transformations have been compared:, with a view toward examining whether the transformation mechanisms have characteristic imprints on these experimental observables. Martensites exhibiting dislocated lath, internally twinned plate, and self-accommodating three-plate cluster morphologies have been encountered in Zr-2.5Nb alloy. Habit planes corresponding to all these morphologies have been found to be consistent with the predictions based on the invariant plane strain (IFS) criterion. Different morphologies have been found to reflect the manner in which the neighboring martensite variants are assembled. Lattice-invariant shears (LISs) for all these cases have been identified to be either {10 (1) over bar 1}(alpha) ((1) over bar 123)(alpha) slip or twinning on (10 (1) over bar 1)(alpha) planes. Widmanstatten alpha precipitates, forming in a step-quenching treatment, have been shown to have a lath morphology, the alpha/beta interface being decorated with a periodic array of (c + a) dislocations at a spacing of 8 to 10 nm. The line vectors of these dislocations are nearly parallel to the invariant lines. The alpha precipitates, forming in the retained beta phase on aging, exhibit an internally twinned structure with a zigzag habit plane. Average habit planes for the morphologies have been found to lie near the {103}(beta) - {113}(beta) poles, which are close to the specific variant of the {112}(beta) plane, which transforms into a prismatic plane of the type {1 (1) over bar 00}(alpha). The crystallography of the formation of the gamma-hydride phase (fct) from both the alpha and beta phases is seen to match the IFS predictions. While the beta-gamma transformation can be treated approximately as a simple shear on the basal plane involving a change in the stacking sequence, the alpha-gamma transformation call be conceptually broken into a alpha --> beta transformation following the Burgers correspondence and the simple beta-gamma shear process. The active eutectoid decomposition in the Zr-Cu system, beta --> alpha + beta', has been described in terms of cooperative growth of the alpha phase from the beta phase through the Burgers correspondence and of the partially ordered beta' (structurally similar to the equilibrium Zr2Cu phase) through an ordering process. Similarities and differences in crystallographic features of these transformations have been discussed. and the importance of the invariant line vector in deciding the geometry of the corresponding habit planes has been pointed out.