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.

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.

Effect of equal channel angular pressing on grain refinement and texture evolution in a biomedical alloy Ti-13Nb-13Zr

Evolution of texture and concomitant grain refinement during Equal Channel Angular Pressing (ECAP) of Ti - 13Nb - 13Zr alloy has been presented. Sub-micron sized equiaxed grains with narrow grain size distribution could be achieved after eight pass at 873 K. A characteristic ECAP texture evolved in alpha phase till four passes while the evolution of characteristic ECAP texture in the beta phase could be observed only beyond the fourth pass. On increasing the deformation up to eight passes, the texture in alpha phase weakens while the beta phase shows an ideal ECAP texture. A weaker texture, low dislocation density and high crystallite size values in alpha phase suggest the occurrence of dynamic recrystallization. The absence of texture evolution in beta phase till four passes can be attributed to local lattice rotations. The characteristic ECAP texture in the eight pass deformed sample is attributed to delayed dynamic recrystallization in the beta phase. (C) 2013 Elsevier Inc. All rights reserved.

Compact polarization dependent EBG surface with fractal boundary patches

Patches with variants of fractal Minkowski curves as boundaries are used here to design a polarization dependent electromagnetic bandgap surface. Reflection phases of the proposed structure depends upon the polarization state of the incident wave and frequency. The phase difference between the x-polarized and y-polarized components of the reflected wave can be as high as 200 degrees and this is achieved without excessive increase in unit cell dimensions and vias. The performance of the surface is analyzed numerically using CST microwave studio. The potential applications of the surface are in polarization conversion surfaces, polarimetric radar calibration, and RCS reduction.

Codes for solving three dimensional linear elastostatic problems using constant boundary elements while ignoring body forces

Three codes, that can solve three dimensional linear elastostatic problems using constant boundary elements while ignoring body forces, are provided here. The file 'bemconst.m' contains a MATLAB code for solving three dimensional linear elastostatic problems using constant boundary elements while ignoring body forces. The file 'bemconst.f90' is a Fortran translation of the MATLAB code contained in the file 'bemconst.m'. The file 'bemconstp.f90' is a parallelized version of the Fortran code contained in the file 'bemconst.f90'. The file 'inbem96.txt' is the input file for the Fortran codes contained in the files 'bemconst.f90' and 'bemconstp.f90'. Author hereby declares that the present codes are the original works of the author. Further, author hereby declares that any of the present codes, in full or in part, is not a translation or a copy of any of the existing codes written by someone else. Author's institution (Indian Institute of Science) has informed the author in writing that the institution is not interested in claiming any copyright on the present codes. Author is hereby distributing the present codes under the MIT License; full text of the license is included in each of the files that contain the codes.

On the k(1)(-1) scaling in sink-flow turbulent boundary layers

Scaling of the streamwise velocity spectrum phi(11)(k(1)) in the so-called sink-flow turbulent boundary layer is investigated in this work. The present experiments show strong evidence for the k(1)(-1) scaling i.e. phi(11)(k(1)) = Lambda(1)U(tau)(2)k(1)(-1), where k(1)(-1) is the streamwise wavenumber and U-tau is the friction velocity. Interestingly, this k(1)(-1) scaling is observed much farther from the wall and at much lower flow Reynolds number (both differing by almost an order of magnitude) than what the expectations from experiments on a zero-pressure-gradient turbulent boundary layer flow would suggest. Furthermore, the coefficient A(1) in the present sink-flow data is seen to be non-universal, i.e. A(1) varies with height from the wall; the scaling exponent -1 remains universal. Logarithmic variation of the so-called longitudinal structure function, which is the physical-space counterpart of spectral k(1)(-1) scaling, is also seen to be non-universal, consistent with the non-universality of A(1). These observations are to be contrasted with the universal value of A(1) (along with the universal scaling exponent of 1) reported in the literature on zero-pressure-gradient turbulent boundary layers. Theoretical arguments based on dimensional analysis indicate that the presence of a streamwise pressure gradient in sink-flow turbulent boundary layers makes the coefficient A(1) non-universal while leaving the scaling exponent -1 unaffected. This effect of the pressure gradient on the streamwise spectra, as discussed in the present study (experiments as well as theory), is consistent with other recent studies in the literature that are focused on the structural aspects of turbulent boundary layer flows in pressure gradients (Harun etal., J. Flui(d) Mech., vol. 715, 2013, pp. 477-498); the present paper establishes the link between these two. The variability of A(1) accommodated in the present framework serves to clarify the ideas of universality of the k(1)(-1) scaling.