Estimating the Hausdorff-Besicovitch dimension of boundary of basin of attraction in helicopter trim

Helicopter trim involves solution of nonlinear force equilibrium equations. As in many nonlinear dynamic systems, helicopter trim problem can show chaotic behavior. This chaotic behavior is found in the basin of attraction of the nonlinear trim equations which have to be solved to determine the main rotor control inputs given by the pilot. This study focuses on the boundary of the basin of attraction obtained for a set of control inputs. We analyze the boundary by considering it at different magnification levels. The magnified views reveal intricate geometries. It is also found that the basin boundary exhibits the characteristic of statistical self-similarity, which is an essential property of fractal geometries. These results led the authors to investigate the fractal dimension of the basin boundary. It is found that this dimension is indeed greater than the topological dimension. From all the observations, it is evident that the boundary of the basin of attraction for helicopter trim problem is fractal in nature. (C) 2012 Elsevier Inc. All rights reserved.

A MATLAB Code for Three Dimensional Linear Elastostatics using Constant Boundary Elements

Present work presents a code written in the very simple programming language MATLAB, for three dimensional linear elastostatics, using constant boundary elements. The code, in full or in part, is not a translation or a copy of any of the existing codes. Present paper explains how the code is written, and lists all the formulae used. Code is verified by using the code to solve a simple problem which has the well known approximate analytical solution. Of course, present work does not make any contribution to research on boundary elements, in terms of theory. But the work is justified by the fact that, to the best of author’s knowledge, as of now, one cannot find an open access MATLAB code for three dimensional linear elastostatics using constant boundary elements. Author hopes this paper to be of help to beginners who wish to understand how a simple but complete boundary element code works, so that they can build upon and modify the present open access code to solve complex engineering problems quickly and easily. The code is available online for open access (as supplementary file for the present paper), and may be downloaded from the website for the present journal.

Deformation-induced thermally activated grain growth in nanocrystalline nickel

Grain growth during indentation at low temperatures has been taken to imply that grain growth is largely stress induced and athermal in nanometals. Indentation experiments on electrodeposited nano-Ni indicate clearly that the load required for grain growth decreases with an increase in temperature, suggesting strongly that concurrent grain growth is thermally activated. (C) 2012 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Continuum theory of edge states of topological insulators: variational principle and boundary conditions

We develop a continuum theory to model low energy excitations of a generic four-band time reversal invariant electronic system with boundaries. We propose a variational energy functional for the wavefunctions which allows us to derive natural boundary conditions valid for such systems. Our formulation is particularly suited for developing a continuum theory of the protected edge/surface excitations of topological insulators both in two and three dimensions. By a detailed comparison of our analytical formulation with tight binding calculations of ribbons of topological insulators modelled by the Bernevig-Hughes-Zhang (BHZ) Hamiltonian, we show that the continuum theory with a natural boundary condition provides an appropriate description of the low energy physics.

Cumulus clouds and convective boundary layers: a tropical perspective on two turbulent shear flows

This paper is a review prepared for the second Marseille Colloquium on the mechanics of turbulence, held in 2011, 50 years after the first. The review covers recent developments in our understanding of the large-scale dynamics of cumulus cloud flows and of the atmospheric boundary layer in the low-wind convective regime that is often encountered in the tropics. It has recently been shown that a variety of cumulus cloud forms and life cycles can be experimentally realized in the laboratory, with the transient diabatic plume taken as the flow model for a cumulus cloud. The plume is subjected to diabatic heating scaled to be dynamically similar to heat release from phase changes in clouds. The experiments are complemented by exact numerical solutions of the Navier-Stokes-Boussinesq equations for plumes with scaled off-source heating. The results show that the Taylor entrainment coefficient first increases with heating, reaches a positive maximum and then drops rapidly to zero or even negative values. This reduction in entrainment is a consequence of structural changes in the flow, smoothing out the convoluted boundaries in the non-diabatic plume, including the tongues engulfing the ambient flow. This is accompanied by a greater degree of mixedness in the core flow because of lower dilution by the ambient fluid. The cloud forms generated depend strongly on the history of the diabatic heating profile in the vertical direction. The striking effects of heating on the flow are attributable to the operation of the baroclinic torque due to the temperature field. The mean baroclinic torque is shown to peak around a quasi-cylindrical sheet situated midway between the axis of the flow and the edges. This torque is shear-enhancing and folds down the engulfment tongues. The increase in mixedness can be traced to an explosive growth in the enstrophy, triggered by a strong fluctuating baroclinic torque that acts as a source, especially at the higher wave numbers, thus enhancing the mixedness. In convective boundary layers field measurements show that, under conditions prevailing in the tropics, the eddy fluxes of momentum and energy do not follow the Monin-Obukhov similarity. Instead, the eddy momentum flux is found to be linear in the wind speed at low winds; and the eddy heat flux is, to a first approximation, governed by free convection laws, with wind acting as a small perturbation on a regime of free convection. A new boundary layer code, based on heat flux scaling rather than wall-stress scaling, shows promising improvements in predictive skills of a general circulation model.

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]

Effect of resonator length and working fluid on the performance of twin thermoacoustic heat engine - experimental and simulation studies

Success in the advancement of thermoacoustic field led the researchers to develop the thermoacoustic engines which found its applications in various fields such as refrigeration, gas mixture separation, natural gas liquefaction, and cryogenics. The objective of this study is to design and fabricate the twin thermoacoustic heat engine (TAHE) producing the acoustic waves with high resonance frequencies which is used to drive a thermoacoustic refrigerator efficiently by the influence of geometrical parameters and working fluids. Twin TAHE has gained significant attention due to the production of high intensity acoustic waves than single TAHE. In order to drive an efficient thermoacoustic refrigerator, a twin thermoacoustic heat engine is built up and its performance are analysed by varying the resonator length and working fluid. The performance is measured in terms of onset temperature difference, resonance frequency and pressure amplitude of the oscillations generated from twin TAHE. The simulation is performed using free software DeltaEC, from LANL, USA. The simulated DeltaEC results are compared with experimental results and the deviations are found within +10%.