Constitutive modeling of shape memory alloy wire with non-local rate kinetics

A constitutive modeling approach for shape memory alloy (SMA) wire by taking into account the microstructural phase inhomogeneity and the associated solid-solid phase transformation kinetics is reported in this paper. The approach is applicable to general thermomechanical loading. Characterization of various scales in the non-local rate sensitive kinetics is the main focus of this paper. Design of SMA materials and actuators not only involve an optimal exploitation of the hysteresis loops during loading-unloading, but also accounts for fatigue and training cycle identifications. For a successful design of SMA integrated actuator systems, it is essential to include the microstructural inhomogeneity effects and the loading rate dependence of the martensitic evolution, since these factors play predominant role in fatigue. In the proposed formulation, the evolution of new phase is assumed according to Weibull distribution. Fourier transformation and finite difference methods are applied to arrive at the analytical form of two important scaling parameters. The ratio of these scaling parameters is of the order of 10(6) for stress-free temperature-induced transformation and 10(4) for stress-induced transformation. These scaling parameters are used in order to study the effect of microstructural variation on the thermo-mechanical force and interface driving force. It is observed that the interface driving force is significant during the evolution. Increase in the slopes of the transformation start and end regions in the stress-strain hysteresis loop is observed for mechanical loading with higher rates.

Multi-port network modeling for fractal shape microstrip antenna

In this paper we demonstrate the use of multi-port network modeling to analyze one such antenna with fractal shaped parts. Based on simulation and experimental studies, it has been demonstrated that model can accurately predict the input characteristics of antennas with Minkowski geometry replacing a side micro strip square ring.

3D shape extraction of human face in presence of facial hair: A profilometric approach

The increasing use of 3D modeling of Human Face in Face Recognition systems, User Interfaces, Graphics, Gaming and the like has made it an area of active study. Majority of the 3D sensors rely on color coded light projection for 3D estimation. Such systems fail to generate any response in regions covered by Facial Hair (like beard, mustache), and hence generate holes in the model which have to be filled manually later on. We propose the use of wavelet transform based analysis to extract the 3D model of Human Faces from a sinusoidal white light fringe projected image. Our method requires only a single image as input. The method is robust to texture variations on the face due to space-frequency localization property of the wavelet transform. It can generate models to pixel level refinement as the phase is estimated for each pixel by a continuous wavelet transform. In cases of sparse Facial Hair, the shape distortions due to hairs can be filtered out, yielding an estimate for the underlying face. We use a low-pass filtering approach to estimate the face texture from the same image. We demonstrate the method on several Human Faces both with and without Facial Hairs. Unseen views of the face are generated by texture mapping on different rotations of the obtained 3D structure. To the best of our knowledge, this is the first attempt to estimate 3D for Human Faces in presence of Facial hair structures like beard and mustache without generating holes in those areas.

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Symmetrical Bisbenzimidazoles with Benzenediyl Spacer: The Role of the Shape of the Ligand on the Stabilization and Structural Alterations in Telomeric G-Quadruplex DNA and Telomerase Inhibition

The extremities of chromosomes end in a G-rich single-stranded overhang that has been implicated in the onset of the replicate senescence. The repeated sequence forming a G-overhang is able to adopt a four-stranded DNA structure called G-quadruplex, which is a poor substrate for the enzyme telomerase. Small molecule based ligands that selectively stabilize the telomeric G-quadruplex DNA, induce telomere shortening eventually leading to cell death. Herein, we have investigated the G-quadruplex DNA interaction with two isomeric bisbenzimidazole-based compounds that differ in terms of shape (V-shaped angular vs linear).While the linear isomer induced some stabilization of the intramolecular G-quadruplex structure generated in the presence of Na+ the other, having V-shaped central planar core, caused a dramatic structural alteration of the latter, above a threshold concentration. This transition was evident from the pronounced changes observed in the circular dichroism spectra and from the get mobility shift assa involving the G-quadruples DNA. Notably, this angular isomer could also induce the G-quadruplex formation in the absence of any added cation. The ligand-quadruples complexes were investigated by computational molecular modeling, providing further information on structure-activity relationships. Finally, TRAP (telomerase repeat amplification protocol) experiments demonstrated that the angular isomer is selective toward the inhibition of telomerase activity.

The Reproducing Kernel DMS-FEM: 3D Shape Functions and Applications to Linear Solid Mechanics

We propose a family of 3D versions of a smooth finite element method (Sunilkumar and Roy 2010), wherein the globally smooth shape functions are derivable through the condition of polynomial reproduction with the tetrahedral B-splines (DMS-splines) or tensor-product forms of triangular B-splines and ID NURBS bases acting as the kernel functions. While the domain decomposition is accomplished through tetrahedral or triangular prism elements, an additional requirement here is an appropriate generation of knotclouds around the element vertices or corners. The possibility of sensitive dependence of numerical solutions to the placements of knotclouds is largely arrested by enforcing the condition of polynomial reproduction whilst deriving the shape functions. Nevertheless, given the higher complexity in forming the knotclouds for tetrahedral elements especially when higher demand is placed on the order of continuity of the shape functions across inter-element boundaries, we presently emphasize an exploration of the triangular prism based formulation in the context of several benchmark problems of interest in linear solid mechanics. In the absence of a more rigorous study on the convergence analyses, the numerical exercise, reported herein, helps establish the method as one of remarkable accuracy and robust performance against numerical ill-conditioning (such as locking of different kinds) vis-a-vis the conventional FEM.

Shape optimization of RC flexural members

A natural velocity field method for shape optimization of reinforced concrete (RC) flexural members has been demonstrated. The possibility of shape optimization by modifying the shape of an initially rectangular section, in addition to variation of breadth and depth along the length, has been explored. Necessary shape changes have been computed using the sequential quadratic programming (SQP) technique. Genetic algorithm (Goldberg and Samtani 1986) has been used to optimize the diameter and number of main reinforcement bars. A limit-state design approach has been adopted for the nonprismatic RC sections. Such relevant issues as formulation of optimization problem, finite-element modeling, and solution procedure have been described. Three design examples-a simply supported beam, a cantilever beam, and a two-span continuous beam, all under uniformly distributed loads-have been optimized. The results show a significant savings (40-56%) in material and cost and also result in aesthetically pleasing structures. This procedure will lead to considerable cost saving, particularly in cases of mass-produced precast members and a heavy cast-in-place member such as a bridge girder.

Modeling Indian Ocean circulation: Bay of Bengal fresh water plume and the Arabian Sea Mini Warm Pool

The Indian subcontinent divides the north Indian Ocean into two tropical basins, namely the Arabian Sea and the Bay of Bengal. The Arabian Sea has high salinity whereas the salinity of the Bay of Bengal is much lower due to the contrast in freshwater forcing of the two basins. The freshwater received by the Bay in large amounts during the summer monsoon through river discharge is flushed out annually by ocean circulation. After the withdrawal of the summer monsoon, the Ganga – Brahmaputra river plume flows first along the Indian coast and then around Sri Lanka into the Arabian Sea creating a low salinity pool in the southeastern Arabian Sea (SEAS). In the same region, during the pre-monsoon months of February – April, a warm pool, known as the Arabian Sea Mini Warm Pool (ASMWP), which is distinctly warmer than the rest of the Indian Ocean, takes shape. In fact, this is the warmest region in the world oceans during this period. Simulation of the river plume and its movement as well as its implications to thermodynamics has been a challenging problem for models of Indian Ocean. Here we address these issues using an ocean general circulation model – first we show that the model is capable of reproducing fresh plumes in the Bay of Bengal as well as its movement and then we use the model to determine the processes that lead to formation of the ASMWP. Hydrographic observations from the western Bay of Bengal have shown the presence of a fresh plume along the northern part of the Indian coast during summer monsoon. The Indian Ocean model when forced by realistic winds and climatological river discharge reproduces the fresh plume with reasonable accuracy. The fresh plume does not advect along the Indian coast until the end of summer monsoon. The North Bay Monsoon Current, which flows eastward in the northern Bay, separates the low salinity water from the more saline southern parts of the bay and thus plays an important role in the fresh water budget of the Bay of Bengal. The model also reproduces the surge of the fresh-plume along the Indian coast, into the Arabian Sea during northeast monsoon. Mechanisms that lead to the formation of the Arabian Sea Mini Warm Pool are investigated using several numerical experiments. Contrary to the existing theories, we find that salinity effects are not necessary for the formation of the ASMWP. The orographic effects of the Sahyadris (Western Ghats) and resulting reduction in wind speed leads to the formation of the ASMWP. During November – April, the SEAS behave as a low-wind heatdominated regime where the evolution of sea surface temperature is solely determined by atmospheric forcing. In such regions the evolution of surface layer temperature is not dependent on the characteristics of the subsurface ocean such as the barrier layer and temperature inversion.

Modeling of growth and motion of equiaxed dendrites in a convecting melt

A numerical micro-scale model is developed to study the behavior of dendrite growth in presence of melt convection. In this method, an explicit, coupled enthalpy model is used to simulate the growth of an equiaxed dendrite, while a Volume of Fluid (VOF) method is used to track the movement of the dendrite in the convecting melt in a two-dimensional Eulerian framework. Numerical results demonstrate the effectiveness of the enthalpy model in simulating the dendritic growth involving complex shape, and the accuracy of VOF method in conserving mass and preserving the complex dendritic shape during motion. Simulations are performed in presence of uniform melt flow for both fixed and moving dendrites, and the difference in dendrite morphology is shown.

Premixed flame response to equivalence ratio fluctuations: Comparison between reduced order modeling and detailed computations

Combustion instability events in lean premixed combustion systems can cause spatio-temporal variations in unburnt mixture fuel/air ratio. This provides a driving mechanism for heat-release oscillations when they interact with the flame. Several Reduced Order Modelling (ROM) approaches to predict the characteristics of these oscillations have been developed in the past. The present paper compares results for flame describing function characteristics determined from a ROM approach based on the level-set method, with corresponding results from detailed, fully compressible reacting flow computations for the same two dimensional slot flame configuration. The comparison between these results is seen to be sensitive to small geometric differences in the shape of the nominally steady flame used in the two computations. When the results are corrected to account for these differences, describing function magnitudes are well predicted for frequencies lesser than and greater than a lower and upper cutoff respectively due to amplification of flame surface wrinkling by the convective Darrieus-Landau (DL) instability. However, good agreement in describing function phase predictions is seen as the ROM captures the transit time of wrinkles through the flame correctly. Also, good agreement is seen for both magnitude and phase of the flame response, for large forcing amplitudes, at frequencies where the DL instability has a minimal influence. Thus, the present ROM can predict flame response as long as the DL instability, caused by gas expansion at the flame front, does not significantly alter flame front perturbation amplitudes as they traverse the flame. (C) 2012 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Experiments and zero D modeling studies using specific Wiebe coefficients for producer gas as fuel in spark-ignited engines

The paper addresses experiments and modeling studies on the use of producer gas, a bio-derived low energy content fuel in a spark-ignited engine. Producer gas, generated in situ, has thermo-physical properties different from those of fossil fuel(s). Experiments on naturally aspirated and turbo-charged engine operation and subsequent analysis of the cylinder pressure traces reveal significant differences in the heat release pattern within the cylinder compared with a typical fossil fuel. The heat release patterns for gasoline and producer gas compare well in the initial 50% but beyond this, producer gas combustion tends to be sluggish leading to an overall increase in the combustion duration. This is rather unexpected considering that producer gas with nearly 20% hydrogen has higher flame speeds than gasoline. The influence of hydrogen on the initial flame kernel development period and the combustion duration and hence on the overall heat release pattern is addressed. The significant deviations in the heat release profiles between conventional fuels and producer gas necessitates the estimation of producer gas-specific Wiebe coefficients. The experimental heat release profiles are used for estimating the Wiebe coefficients. Experimental evidence of lower fuel conversion efficiency based on the chemical and thermal analysis of the engine exhaust gas is used to arrive at the Wiebe coefficients. The efficiency factor a is found to be 2.4 while the shape factor m is estimated at 0.7 for 2% to 90% burn duration. The standard Wiebe coefficients for conventional fuels and fuel-specific coefficients for producer gas are used in a zero D model to predict the performance of a 6-cylinder gas engine under naturally aspirated and turbo-charged conditions. While simulation results with standard Wiebe coefficients result in excessive deviations from the experimental results, excellent match is observed when producer gas-specific coefficients are used. Predictions using the same coefficients on a 3-cylinder gas engine having different geometry and compression ratio(s) indicate close match with the experimental traces highlighting the versatility of the coefficients.

Morphing fixed wing MAV modeling using VAM

The design and implementation of a morphing Micro Air Vehicle (MAV) wing using a smart composite is attempted in this research work. Control surfaces actuated by traditional servos are difficult to instrument and fabricate on thin composite-wings of MAVs. Piezoelectric Fiber Reinforced Composites (PFRCs) are the chosen smart structural materials in the current work for incorporation onto fixed-wing MAVs to simultaneously perform the dual functions of structural load-bearing and actuation of flexure, torsion and/or extension for morphing. Further, PFRC use can be extended towards shape control of a “fixed” wing MAV to meet changing performance requirements. Wings that can warp into desired shapes and/or have variable camber are well-known to exhibit improved efficiency in aerodynamic control. During an entire flight cycle, there are multiple optimal configurations, each of which suits a particular phase of the flight regime. Widely proposed methods of wing morphing include changes in camber, twist, sweep and span. However, camber change during flight is already established, in terms of its potential, as a major factor in improving the aerofoil efficiency and flow separation behavior. Hence, for this work, morphing with camber change is adopted with the goal to better tailor aerodynamic properties.

Hot deformation behavior of Ni-Fe-Ga-based ferromagnetic shape memory alloy - A study using processing map

Ni-Fe-Ga-based alloys form a new class of ferromagnetic shape memory alloys (FSMAs) that show considerable formability because of the presence of a disordered fcc gamma-phase. The current study explores the deformation processing of this alloy using an off-stoichiometric Ni55Fe59Ga26 alloy that contains the ductile gamma-phase. The hot deformation behavior of this alloy has been characterized on the basis of its flow stress variation obtained by isothermal constant true strain rate compression tests in the 1123-1323 K temperature range and strain rate range of 10(-3)-10 s(-1) and using a combination of constitutive modeling and processing map. The dynamic recrystallization (DRX) regime for thermomechanical processing has been identified for this Heusler alloy on the basis of the processing maps and the deformed microstructures. This alloy also shows evidence of dynamic strain-aging (DSA) effect which has not been reported so far for any Heusler FSMAs. Similar effect is also noticed in a Ni-Mn-Ga-based Heusler alloy which is devoid of any gamma-phase. (C) 2014 Elsevier Ltd. All rights reserved.

Effects of Interface Treatment on the Fatigue Behavior of Shape Memory Alloy Reinforced Polymer Composites

Interfacial properties of Shape Memory Alloy (SMA) reinforced polymer matrix composites can be enhanced by improving the interfacial bonding. This paper focuses on studying the interfacial stresses developed in the SMA-epoxy interface due to various laser shot penning conditions. Fiber-pull test-setup is designed to understand the role of mechanical bias stress cycling and thermal actuation cycling. Phase transformation is tracked over mechanical and thermal fatigue cycles. A micromechanics based model developed earlier based on shear lag in SMA and energy based consistent homogenization is extended here to incorporate the stress-temperature phase diagram parameters for modeling fatigue.