A numerical approach to determine the sufficiency of given boundary data sets for uniquely estimating interior elastic properties

We consider an inverse elasticity problem in which forces and displacements are known on the boundary and the material property distribution inside the body is to be found. In other words, we need to estimate the distribution of constitutive properties using the finite boundary data sets. Uniqueness of the solution to this problem is proved in the literature only under certain assumptions for a given complete Dirichlet-to-Neumann map. Another complication in the numerical solution of this problem is that the number of boundary data sets needed to establish uniqueness is not known even under the restricted cases where uniqueness is proved theoretically. In this paper, we present a numerical technique that can assess the sufficiency of given boundary data sets by computing the rank of a sensitivity matrix that arises in the Gauss-Newton method used to solve the problem. Numerical experiments are presented to illustrate the method.

A QUADRATIC C degrees INTERIOR PENALTY METHOD FOR LINEAR FOURTH ORDER BOUNDARY VALUE PROBLEMS WITH BOUNDARY CONDITIONS OF THE CAHN-HILLIARD TYPE

We develop a quadratic C degrees interior penalty method for linear fourth order boundary value problems with essential and natural boundary conditions of the Cahn-Hilliard type. Both a priori and a posteriori error estimates are derived. The performance of the method is illustrated by numerical experiments.

Fidelity susceptibility of one-dimensional models with twisted boundary conditions

Recently it has been shown that the fidelity of the ground state of a quantum many-body system can be used todetect its quantum critical points (QCPs). If g denotes the parameter in the Hamiltonian with respect to which the fidelity is computed, we find that for one-dimensional models with large but finite size, the fidelity susceptibility chi(F) can detect a QCP provided that the correlation length exponent satisfies nu < 2. We then show that chi(F) can be used to locate a QCP even if nu >= 2 if we introduce boundary conditions labeled by a twist angle N theta, where N is the system size. If the QCP lies at g = 0, we find that if N is kept constant, chi(F) has a scaling form given by chi(F) similar to theta(-2/nu) f (g/theta(1/nu)) if theta << 2 pi/N. We illustrate this both in a tight-binding model of fermions with a spatially varying chemical potential with amplitude h and period 2q in which nu = q, and in a XY spin-1/2 chain in which nu = 2. Finally we show that when q is very large, the model has two additional QCPs at h = +/- 2 which cannot be detected by studying the energy spectrum but are clearly detected by chi(F). The peak value and width of chi(F) seem to scale as nontrivial powers of q at these QCPs. We argue that these QCPs mark a transition between extended and localized states at the Fermi energy. DOI: 10.1103/PhysRevB.86.245424

Temperature dependent structural and dielectric investigations of PbZr0.5Ti0.5O3 solid solution at the morphotropic phase boundary

PbZr1-xTixO3 ceramics synthesised by low temperature calcination followed by sintering at 1280 degrees C show a Morphotropic Phase Boundary (MPB) for compositions of x=0.44-0.51. The morphotropic phase boundary is wider for samples with smaller grain sizes due to the synthesis route. A Rietveld analysis is performed on a composition of x=0.5 composition to quantify the phase fractions of the tetragonal and monoclinic phases present in the PZT system. Temperature dependent X-ray diffraction and dielectric studies of PbZr0.5Ti0.5O3 composition demonstrated a phase transformation from monoclinic to tetragonal at 270 degrees C followed by a ferroelectric tetragonal to a paraelectric cubic transition at 370 degrees C. Thus, the poling of these ceramics should be performed below 270 degrees C to benefit from the presence of a monoclinic phase. (C) 2012 Elsevier Ltd and Techna Group S.r.l. All rights reserved.

Interdiffusion and the phase boundary compositions in the Co-To system

Diffusion parameters such as the interdiffusion coefficients and the ratio of the tracer diffusion coefficients are calculated in the Co2Ta Laves phase. The activation energy for the interdiffusion coefficients is calculated as 186 +/- 29 kJ/mol. The ratio of tracer diffusion coefficients indicates that Co has higher diffusion rate than that of Ta. This is explained with the help of possible point defects and the crystal structure of the phase: The phase boundary compositions measured in this study is different from the compositions published previously. (C) 2012 Elsevier Ltd. All rights reserved.

Optimal control problem for the time-dependent Kirchhoff-Love plate in a domain with rough boundary

In this paper, we study the asymptotic behavior of an optimal control problem for the time-dependent Kirchhoff-Love plate whose middle surface has a very rough boundary. We identify the limit problem which is an optimal control problem for the limit equation with a different cost functional.

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.

Analytical solutions for diffuse fluorescence spectroscopy/imaging in biological tissues. Part I: zero and extrapolated boundary conditions

The mathematical model for diffuse fluorescence spectroscopy/imaging is represented by coupled partial differential equations (PDEs), which describe the excitation and emission light propagation in soft biological tissues. The generic closed-form solutions for these coupled PDEs are derived in this work for the case of regular geometries using the Green's function approach using both zero and extrapolated boundary conditions. The specific solutions along with the typical data types, such as integrated intensity and the mean time of flight, for various regular geometries were also derived for both time-and frequency-domain cases. (C) 2013 Optical Society of America

Performance evaluation of chemistry transport models over India

Using continuous and near-real time measurements of the mass concentrations of black carbon (BC) aerosols near the surface, for a period of 1 year (from January to December 2006) from a network of eight observatories spread over different environments of India, a space-time synthesis is generated. The strong seasonal variations observed, with a winter high and summer low, are attributed to the combined effects of changes in synoptic air mass types, modulated strongly by the atmospheric boundary layer dynamics. Spatial distribution shows much higher BC concentration over the Indo-Gangetic Plain (IGP) than the peninsular Indian stations. These were examined against the simulations using two chemical transport models, GOCART (Goddard Global Ozone Chemistry Aerosol Radiation and Transport) and CHIMERE for the first time over Indian region. Both the model simulations significantly deviated from the measurements at all the stations; more so during the winter and pre-monsoon seasons and over mega cities. However, the CHIMERE model simulations show better agreement compared with the measurements. Notwithstanding this, both the models captured the temporal variations; at seasonal and subseasonal timescales and the natural variabilities (intra-seasonal oscillations) fairly well, especially at the off-equatorial stations. It is hypothesized that an improvement in the atmospheric boundary layer (ABL) parameterization scheme for tropical environment might lead to better results with GOCART.

Ages and relative sizes of pre-2004 tsunamis in the Bay of Bengal inferred from geologic evidence in the Andaman and Nicobar Islands

Geologic evidence along the northern part of the 2004 Aceh-Andaman rupture suggests that this region generated as many as five tsunamis in the prior 2000years. We identify this evidence by drawing analogy with geologic records of land-level change and the tsunami in 2004 from the Andaman and Nicobar Islands (A&N). These analogs include subsided mangrove swamps, uplifted coral terraces, liquefaction, and organic soils coated by sand and coral rubble. The pre-2004 evidence varies in potency, and materials dated provide limiting ages on inferred tsunamis. The earliest tsunamis occurred between the second and sixth centuries A.D., evidenced by coral debris of the southern Car Nicobar Island. A subsequent tsunami, probably in the range A.D. 770-1040, is inferred from deposits both in A&N and on the Indian subcontinent. It is the strongest candidate for a 2004-caliber earthquake in the past 2000years. A&N also contain tsunami deposits from A.D. 1250 to 1450 that probably match those previously reported from Sumatra and Thailand, and which likely date to the 1390s or 1450s if correlated with well-dated coral uplift offshore Sumatra. Thus, age data from A&N suggest that within the uncertainties in estimating relative sizes of paleo-earthquakes and tsunamis, the 1000year interval can be divided in half by the earthquake or earthquakes of A.D. 1250-1450 of magnitude >8.0 and consequent tsunamis. Unlike the transoceanic tsunamis generated by full or partial rupture of the subduction interface, the A&N geology further provides evidence for the smaller-sized historical tsunamis of 1762 and 1881, which may have been damaging locally.

Exact internal controllability for a hyperbolic problem in a domain with highly oscillating boundary

In this paper, by using the Hilbert Uniqueness Method (HUM), we study the exact controllability problem described by the wave equation in a three-dimensional horizontal domain bounded at the bottom by a smooth wall and at the top by a rough wall. The latter is assumed to consist in a plane wall covered with periodically distributed asperities whose size depends on a small parameter epsilon > 0, and with a fixed height. Our aim is to obtain the exact controllability for the homogenized equation. In the process, we study the asymptotic analysis of wave equation in two setups, namely solution by standard weak formulation and solution by transposition method.

Grain growth and grain boundary dynamics in colloidal polycrystals

Grain boundary dynamics and grain growth play a pivotal role in the fabrication of functional polycrystalline materials. However, not much is known about the delicate interplay between various microscopic processes that drive grain boundary motion which eventually culminates in the desired grain morphology. Colloidal systems are ideally suited to bridge the gap between the microscopic and macroscopic processes underlying grain growth, since their dynamics can be followed in real space and real time with single-particle resolution. The present review aims at highlighting contributions from colloid experiments that have led to a holistic understanding of grain growth in polycrystalline materials.

Rotating beams and non-rotating beams with shared eigenpair for pinned-free boundary condition

In this paper we look for a rotating beam, with pinned-free boundary conditions, whose eigenpair (frequency and mode-shape) is same as that of a uniform non-rotating beam for a particular mode. It is seen that for any given mode, there exists a flexural stiffness function (FSF) for which the ith mode eigenpair of a rotating beam with uniform mass distribution, is identical to that of a corresponding non-rotating beam with same length and mass distribution. Inserting these derived FSF's in a finite element code for a rotating pinned-free beam, the frequencies and mode shapes of a non-rotating pinned-free beam are obtained. For the first mode, a physically realistic equivalent rotating beam is possible, but for higher modes, the FSF has internal singularities. Strategies for addressing these singularities in the FSF for finite element analysis are provided. The proposed functions can be used as test functions for rotating beam codes and also for targeted destiffening of rotating beams.

Metastable morphotropic phase boundary state in the multiferroic BiFeO3-PbTiO3

Temperature-time study of the magnetoelectric multiferroic (1-x)BiFeO3-(x)PbTiO3 by x-ray and electron diffraction on the reported morphotropic phase boundary (MPB) compositions revealed that this MPB does not correspond to the equilibrium state. The MPB like state is rather of metastable nature and arise due to kinetic arrest of metastable rhombohedral (R3c) phase, along with the equilibrium tetragonal (P4mm) phase. The life time of the metastable R3c nuclei is very sensitive to composition and temperature, and nearly diverges at x -> 0.27. The MPB like state appears only if the system is cooled before the metastable R3c nuclei could vanish. These findings resolve the long standing controversy with regard to seemingly erratic phase formation behaviour reported by different groups and provides a rational basis for developing genuine equilibrium MPB compositions in this system for better piezoelectric properties. (C) 2013 AIP Publishing LLC.

A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs

In this work, first a Fortran code is developed for three dimensional linear elastostatics using constant boundary elements; the code is based on a MATLAB code developed by the author earlier. Next, the code is parallelized using BLACS, MPI, and ScaLAPACK. Later, the parallelized code is used to demonstrate the usefulness of the Boundary Element Method (BEM) as applied to the realtime computational simulation of biological organs, while focusing on the speed and accuracy offered by BEM. A computer cluster is used in this part of the work. The commercial software package ANSYS is used to obtain the `exact' solution against which the solution from BEM is compared; analytical solutions, wherever available, are also used to establish the accuracy of BEM. A pig liver is the biological organ considered. Next, instead of the computer cluster, a Graphics Processing Unit (GPU) is used as the parallel hardware. Results indicate that BEM is an interesting choice for the simulation of biological organs. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature. Also, a serial MATLAB code, and both serial and parallel versions of a Fortran code, which can solve three dimensional (3D) linear elastostatic problems using constant boundary elements, are provided as supplementary files that can be freely downloaded.