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.

Sliding-Mode Guidance and Control for All-Aspect Interceptors with Terminal Angle Constraints

In this paper, sliding-mode-control-based guidance laws to intercept stationary, constant-velocity, and maneuvering targets at a desired impact angle are proposed. The desired impact angle, which is defined in terms of a desired line-of-sight angle, is achieved in finite time by selecting the missile's lateral acceleration to enforce terminal sliding mode on a switching surface designed using nonlinear engagement dynamics. The conditions for capturability are also presented. In addition, by considering a three-degree-of-freedom linear-interceptor dynamic model and by following the procedure used to design a dynamic sliding-mode controller, the interceptor autopilot is designed as a simple static controller to track the lateral acceleration generated by the guidance law. Numerical simulation results are presented to validate the proposed guidance laws and the autopilot design for different initial engagement geometries and impact angles.

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.

Effect of Interfacial Microstructure of Adsorbed Poly(ethylene glycol) Monooleate on Steel Substrate on Sliding Friction in Oil in Water Emulsion

The concentration of a nonionic surfactant and water pH were varied in an oil-in-water emulsion to minimize the friction coefficient between a steel ball sliding on a steel flat. At a surfactant concentration near the CMC (critical micelle concentration) the oil droplet size was found to be minimum. In this paper we study the microstructure of the surfactant molecules self-assembled on the steel substrate in water to comment on the ability of the surfactant aggregate to attract and retain oil. We find that a near semicylindrical hemimiceller microstructure with hydrocarbon tails projecting into bulk water as obtained at CMC in near neutral water is best able to capture and retain oil in yielding a low coefficient of friction.

Microstructure and dry sliding wear resistance evaluation of plasma nitrided austenitic stainless steel type AISI 316LN against different sliders

In this work, Plasma Nitriding was carried out at a temperature of 570 degrees C on nuclear grade austenitic stainless steel type AISI 316 LN (316LN SS) in a gas mixture of 20% N-2-80% H-2 to improve the surface hardness and thereby sliding wear resistance. The Plasma Nitride (PN) treated surface has been characterized by Vickers microhardness measurements, Scanning Electron Microscopic (SEM) examination, X-ray Diffraction (XRD) and sliding wear assessment. The average thickness of the PN layer was found to be 70 mu m. Microhardness measurements showed a significant increase in the hardness from 210 HV25g (unnitrided sample) to 1040 HV25g (Plasma Nitrided sample). The XRD reveals that PN layer consists of CrN, Fe4N and Fe3N phases along with austenite phase. The tribological parameters such as the friction coefficient and wear mechanism have been evaluated at ambient conditions for PN treated ring (PN ring) vs. ASTM A453 grade 660 pin (ASTM pin), PN ring vs. Nickel based alloy hard faced pin (Colmonoy pin), PN ring vs. 316LN SS pin and 316LN SS ring vs. 316LN SS pin. The wear tracks have been analyzed by SEM, Energy Dispersive X-ray Analysis (EDX) and Optical Profilometry. The untreated 316LN SS ring vs. 316LN SS pin produced severe wear and was characterized by a combination of delamination and adhesion wear mechanism, whereas wear mechanism of the PN rings reveals mild abrasion and a transfer layer from pin materials. (C) 2012 Elsevier B.V. All rights reserved.

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.