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

Directional grain growth from anisotropic kinetic roughening of grain boundaries in sheared colloidal crystals

The fabrication of functional materials via grain growth engineering implicitly relies on altering the mobilities of grain boundaries (GBs) by applying external fields. Although computer simulations have alluded to kinetic roughening as a potential mechanism for modifying GB mobilities, its implications for grain growth have remained largely unexplored owing to difficulties in bridging the widely separated length and time scales. Here, by imaging GB particle dynamics as well as grain network evolution under shear, we present direct evidence for kinetic roughening of GBs and unravel its connection to grain growth in driven colloidal polycrystals. The capillary fluctuation method allows us to quantitatively extract shear-dependent effective mobilities. Remarkably, our experiments reveal that for sufficiently large strains, GBs with normals parallel to shear undergo preferential kinetic roughening, resulting in anisotropic enhancement of effective mobilities and hence directional grain growth. Single-particle level analysis shows that the mobility anisotropy emerges from strain-induced directional enhancement of activated particle hops normal to the GB plane. We expect our results to influence materials fabrication strategies for atomic and block copolymeric polycrystals as well.

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

Competing structural phase transition scenarios in the giant tetragonality ferroelectric BiFeO3-PbTiO3: Isostructural vs multiphase transition

A systematic x-ray and neutron powder diffraction study of the giant tetragonality multiferroic (1-x) BiFeO3-(x) PbTiO3 have revealed that the compositions close to the morphotropic phase boundary present two different structural phase transition scenarios on cooling from the cubic phase: (i) cubic -> tetragonal (T-2) + tetragonal (T-1) -> tetragonal (T-1) and (ii) cubic -> tetragonal (T-2) + tetragonal (T-1) + rhombohedral (R3c) -> tetragonal (T-1) + rhombohedral (R3c). The comparatively larger tetragonality (c/a - 1) of the T-1 phase as compared to the coexisting isostructural T-2 phase is shown to be a result of significantly greater degree of overlap of the Pb/Bi-6s and Ti/Fe-3d with the O-2p orbitals as compared to that in the T-2 phase. The formation/suppression of the minor metastable rhombohedral phase seems to be governed by subtle play of local kinetic factors. In the scenario when the minor rhombohedral (R) phase is formed along with the tetragonal phases it is able to accommodate the large transformation stress in the system due to formation of the tetragonal phases, and prevent the solid from disintegration into powder after sintering. When the metastable rhombohedral phase is not formed, the large transformation strain ruptures the grain boundaries leading to fragmentation of the dense solid to powder. (C) 2013 American Institute of Physics. http://dx.doi.org/10.1063/1.4792215]

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.

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.

Deformation-induced grain growth and twinning in nanocrystalline palladium thin films

The microstructure and mechanical properties of nanocrystalline Pd films prepared by magnetron sputtering have been investigated as a function of strain. The films were deposited onto polyimide substrates and tested in tensile mode. In order to follow the deformation processes in the material, several samples were strained to defined straining states, up to a maximum engineering strain of 10%, and prepared for post-mortem analysis. The nanocrystalline structure was investigated by quantitative automated crystal orientation mapping (ACOM) in a transmission electron microscope (TEM), identifying grain growth and twinning/detwinning resulting from dislocation activity as two of the mechanisms contributing to the macroscopic deformation. Depending on the initial twin density, the samples behaved differently. For low initial twin densities, an increasing twin density was found during straining. On the other hand, starting from a higher twin density, the twins were depleted with increasing strain. The findings from ACOM-TEM were confirmed by results from molecular dynamics (MD) simulations and from conventional and in-situ synchrotron X-ray diffraction (CXRD, SXRD) experiments.

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.

Extension of the universal erosive burning law to partly symmetric propellant grain geometries

This paper aims at extending the universal erosive burning law developed by two of the present authors from axi-symmetric internally burning grains to partly symmetric burning grains. This extension revolves around three dimensional flow calculations inside highly loaded grain geometry and benefiting from an observation that the flow gradients normal to the surface in such geometries have a smooth behavior along the perimeter of the grain. These are used to help identify the diameter that gives the same perimeter the characteristic dimension rather than a mean hydraulic diameter chosen earlier. The predictions of highly loaded grains from the newly chosen dimension in the erosive burning law show better comparison with measured pressure-time curves while those with mean hydraulic diameter definitely over-predict the pressures. (c) 2013 IAA. Published by Elsevier Ltd. All rights reserved.

Effect of Texture and Grain Size on Bio-Corrosion Response of Ultrafine-Grained Titanium

The bio-corrosion response of ultrafine-grained commercially pure titanium processed by different routes of equal-channel angular pressing has been studied in simulated body fluid. The results indicate that the samples processed through route B-c that involved rotation of the workpiece by 90 deg in the same sense between each pass exhibited higher corrosion resistance compared to the ones processed by other routes of equal-channel angular pressing, as well as the coarse-grained sample. For a similar grain size, the higher corrosion resistance of the samples exhibiting off-basal texture compared to shear texture indicates the major role of texture in corrosion behavior. It is postulated that an optimum combination of microstructure and crystallographic texture can lead to high strength and excellent corrosion resistance.