Isotropic finite volume discretization

Finite volume methods traditionally employ dimension by dimension extension of the one-dimensional reconstruction and averaging procedures to achieve spatial discretization of the governing partial differential equations on a structured Cartesian mesh in multiple dimensions. This simple approach based on tensor product stencils introduces an undesirable grid orientation dependence in the computed solution. The resulting anisotropic errors lead to a disparity in the calculations that is most prominent between directions parallel and diagonal to the grid lines. In this work we develop isotropic finite volume discretization schemes which minimize such grid orientation effects in multidimensional calculations by eliminating the directional bias in the lowest order term in the truncation error. Explicit isotropic expressions that relate the cell face averaged line and surface integrals of a function and its derivatives to the given cell area and volume averages are derived in two and three dimensions, respectively. It is found that a family of isotropic approximations with a free parameter can be derived by combining isotropic schemes based on next-nearest and next-next-nearest neighbors in three dimensions. Use of these isotropic expressions alone in a standard finite volume framework, however, is found to be insufficient in enforcing rotational invariance when the flux vector is nonlinear and/or spatially non-uniform. The rotationally invariant terms which lead to a loss of isotropy in such cases are explicitly identified and recast in a differential form. Various forms of flux correction terms which allow for a full recovery of rotational invariance in the lowest order truncation error terms, while preserving the formal order of accuracy and discrete conservation of the original finite volume method, are developed. Numerical tests in two and three dimensions attest the superior directional attributes of the proposed isotropic finite volume method. Prominent anisotropic errors, such as spurious asymmetric distortions on a circular reaction-diffusion wave that feature in the conventional finite volume implementation are effectively suppressed through isotropic finite volume discretization. Furthermore, for a given spatial resolution, a striking improvement in the prediction of kinetic energy decay rate corresponding to a general two-dimensional incompressible flow field is observed with the use of an isotropic finite volume method instead of the conventional discretization. (C) 2014 Elsevier Inc. All rights reserved.

A perturbed martingale approach to global optimization

A new global stochastic search, guided mainly through derivative-free directional information computable from the sample statistical moments of the design variables within a Monte Carlo setup, is proposed. The search is aided by imparting to the directional update term additional layers of random perturbations referred to as `coalescence' and `scrambling'. A selection step, constituting yet another avenue for random perturbation, completes the global search. The direction-driven nature of the search is manifest in the local extremization and coalescence components, which are posed as martingale problems that yield gain-like update terms upon discretization. As anticipated and numerically demonstrated, to a limited extent, against the problem of parameter recovery given the chaotic response histories of a couple of nonlinear oscillators, the proposed method appears to offer a more rational, more accurate and faster alternative to most available evolutionary schemes, prominently the particle swarm optimization. (C) 2014 Elsevier B.V. All rights reserved.

Structural and magnetic properties of ultra-small scale eutectic CoFeZr alloys

Aiming to develop high mechanical strength and toughness by tuning ultrafine lamellar spacing of magnetic eutectic alloys, we report the mechanical and magnetic properties of the binary eutectic alloys Co90.5Zr9.5 and Fe90.2Zr9.8, as well as the pseudo-binary eutectic alloys Co82.4Fe8Zr9.6, Co78Fe12.4Zr9.6 and Co49.2Fe49.2Zr9.6 developed by suction-casting. The lower lamellar spacing around 100 nm of the eutectics Co49.2Fe49.2Zr9.6 yields a high hardness of 713(+/- 20) VHN. Magnetic measurements reveal high magnetic moment of 1.92 mu B (at 5 K) and 1.82 mu B (at 300 K) per formula unit for this composition. The magnetization vs. applied field data at 5 K show a directional preference to some extent and therefore smaller non-collinear magnetization behavior compared to Co11Zr2 reported in the literature due to exchange frustration and transverse spin freezing owing to the presence of smaller Zr content. The decay of magnetization as a function of temperature along the easy axis of magnetization of all the eutectic compositions can be described fairly well by the spin wave excitation equation Delta M/M(0) = BT3/2 + CT5/2. (C) 2014 Elsevier B.V. All rights reserved.

Capturability Analysis of a 3-D Retro-PN Guidance Law for Higher Speed Nonmaneuvering Targets

This brief presents the capturability analysis of a 3-D Retro-proportional navigation (Retro-PN) guidance law, which uses a negative navigation constant (as against the usual positive one), for intercepting targets having higher speeds than interceptors. This modification makes it possible to achieve collision conditions that were inaccessible to the standard PN law. A modified polar coordinate system, that makes the model more compact, is used in this brief for capturability analysis. In addition to the ratio of the target to interceptor speeds, the directional cosines of the interceptor, and target velocity vectors play a crucial role in the capturability. The existence of nontrivial capture zone of the Retro-PN guidance law and necessary and sufficient conditions, for capturing the target in finite time, are presented. A sufficient condition on the navigation constant is derived to ensure finiteness of the line-of-sight turn rate. The results are more extensive than those available for 2-D engagements, which can be obtained as special cases of this brief. Simulation results are given to support the analytical results.

Investigations on pressure and thickness profiles in carbon fiber-reinforced polymers during vacuum assisted resin transfer molding

Systematic experiments have been carried out by monitoring the in-situ pressure and thickness profiles for three different configurations, viz., flat plate, flat plate with a central circular hole, and an L-section using vacuum assisted resin transfer molding (VARTM) process. The effect of anisotropy on resin flow has been quantified by considering uni-directional carbon fiber preforms with 0 degrees and 90 degrees orientation to the flow direction for each configuration. A quasi-isotropic 45 degrees/0 degrees/-45 degrees/90 degrees](S) layup has also been included for flat plate case. Additionally, the study has been extended to understand the effect of using high permeability medium for each configuration. Fluid pressure profiles and thickness variation profiles have been obtained using an array of pressure sensors and linear variable differential transformers for each configuration. Experimental data reveal that anisotropy (due to changing fiber orientations), configuration, and gravity significantly change fluid pressure and displacement fields obtained during VARTM.

Interphase anisotropy effects on lamellar eutectics: A numerical study

In directional solidification of binary eutectics, it is often observed that two-phase lamellar growth patterns grow tilted with respect to the direction z of the imposed temperature gradient. This crystallographic effect depends on the orientation of the two crystal phases alpha and beta with respect to z. Recently, an approximate theory was formulated that predicts the lamellar tilt angle as a function of the anisotropy of the free energy of the solid(alpha)-solid(beta) interphase boundary. We use two different numerical methods-phase field (PF) and dynamic boundary integral (BI)-to simulate the growth of steady periodic patterns in two dimensions as a function of the angle theta(R) between z and a reference crystallographic axis for a fixed relative orientation of alpha and beta crystals, that is, for a given anisotropy function (Wulff plot) of the interphase boundary. For Wulff plots without unstable interphase-boundary orientations, the two simulation methods are in excellent agreement with each other and confirm the general validity of the previously proposed theory. In addition, a crystallographic ``locking'' of the lamellae onto a facet plane is well reproduced in the simulations. When unstable orientations are present in the Wulff plot, it is expected that two distinct values of the tilt angle can appear for the same crystal orientation over a finite theta(R) range. This bistable behavior, which has been observed experimentally, is well reproduced by BI simulations but not by the PF model. Possible reasons for this discrepancy are discussed.

Giant actuation in bulk carbon nanotubes under coupled electric and magnetic fields

We report an enhanced actuation in bulk carbon nanotubes (CNTs) under coupled electric and magnetic fields, which is much higher than that evaluated in the presence of individual fields. Coupled electric and magnetic fields induce a directional actuation demonstrating a transformation from polarity independent to dependent actuation behavior of CNTs. Both qualitative and quantitative analyses are performed to understand this transformation in the bulk CNTs. Moreover, actuations along radial and axial directions of CNTs have also demonstrated a similar directional behavior.

A comparative study of failure features in aerospace grade unidirectional and bidirectional woven CFRP composite laminates under four-point bend fatigue loads

Damage mechanisms in unidirectional (UD) and bi-directional (BD) woven carbon fiber reinforced polymer (CFRP) laminates subjected to four point flexure, both in static and fatigue loadings, were studied. The damage progression in composites was monitored by observing the slopes of the load vs. deflection data that represent the stiffness of the given specimen geometry over a number of cycles. It was observed that the unidirectional composites exhibit gradual loss in stiffness whereas the bidirectional woven composites show a relatively quicker loss during stage II of fatigue damage progression. Both, the static and the fatigue failures in unidirectional carbon fiber reinforced polymer composites originates due to generation of cracks on compression face while in bidirectional woven composites the damage ensues from both the compression and the tensile faces. These observations are supported by a detailed fractographic analysis.

Algorithm for Autonomous Reorganization of Mobile Wireless Camera Sensor Networks to Improve Coverage

In this paper, sensing coverage by wireless camera-embedded sensor networks (WCSNs), a class of directional sensors is studied. The proposed work facilitates the autonomous tuning of orientation parameters and displacement of camera-sensor nodes in the bounded field of interest (FoI), where the network coverage in terms of every point in the FoI is important. The proposed work is first of its kind to study the problem of maximizing coverage of randomly deployed mobile WCSNs which exploits their mobility. We propose an algorithm uncovered region exploration algorithm (UREA-CS) that can be executed in centralized and distributed modes. Further, the work is extended for two special scenarios: 1) to suit autonomous combing operations after initial random WCSN deployments and 2) to improve the network coverage with occlusions in the FoI. The extensive simulation results show that the performance of UREA-CS is consistent, robust, and versatile to achieve maximum coverage, both in centralized and distributed modes. The centralized and distributed modes are further analyzed with respect to the computational and communicational overheads.

Anisotropy in the mechanical properties of organic crystals: temperature dependence

The nanoindentation technique has recently been utilized for quantitative evaluation of the mechanical properties of molecular materials successfully, including their temperature (T) dependence. In this paper, we examine how the mechanical anisotropy varies with T in saccharin and L-alanine single crystals. Our results show that elastic modulus (E) decreases linearly in all the cases examined, with the T-dependence of E being anisotropic. Correspondence between directional dependence of the slopes of the E vs. T plots and the linear thermal expansion coefficients was found. The T-dependence of hardness (H), on the other hand, was found to be nonlinear and significant when (100) of saccharin and (001) of L-alanine are indented. While the anisotropies in E and H of saccharin and E of L-alanine enhance with T, the anisotropy in H of L-alanine was found to reduce with T. Possible mechanistic origins of these variations are discussed.

Halogen Bonding and Chalcogen Bonding in 4,7-Dibromo-5, 6-dinitro-2,1,3-benzothiadiazole

An organic solid, 4,7-dibromo-5,6-dinitro-2,1,3-benzothiadiazole, has been designed to serve as an illustrative example to quantitatively evaluate the relative merits of halogen and chalcogen bonding in terms of charge density features. The compound displays two polymorphic modifications, one crystallizing in a non-centrosymmetric space group (Z' = 1) and the other in a centrosymmetric space group with two molecules in the asymmetric unit (Z' = 2). Topological analysis based on QTAIM clearly brings out the dominance of the chalcogen bond over the halogen bond along with an indication that halogen bonds are more directional compared to chalcogen bonds. The cohesive energies calculated with the absence of both strong and weak hydrogen bonds as well as stacking interaction are indicative of the stabilities associated with the polymorphic forms.

Effect of Surface Energy Anisotropy on the Stability of Growth Fronts in Multiphase Alloys

Eutectic growth offers a variety of examples for pattern formation which are interesting both for theoreticians as well as experimentalists. One such example of patterns is ternary eutectic colonies which arise as a result of instabilities during growth of two solid phases. Here, in addition to the two major components being exchanged between the solid phases during eutectic growth, there is an impurity component which is rejected by both solid phases. During progress of solidification, there develops a boundary layer of the third impurity component ahead of the solidification front of the two solid phases. Similar to Mullins-Sekerka type instabilities, such a boundary layer tends to make the global solidification envelope unstable to morphological perturbations giving rise to two-phase cells. This phenomenon has been studied numerically in two dimensions for the conditions of directional solidification, by Plapp and Karma (Phys Rev E 66:061608, 2002) using phase-field simulations. While, in the work by Plapp and Karma (Phys Rev E 66:061608, 2002) all interfaces are isotropic, in our presentation, we extend the phase-field model by considering interfacial anisotropy in the solid-solid and solid-liquid interfaces and characterize the role of interfacial anisotropy on the stability of the growth front through phase-field simulations in two dimensions.

Pattern-Formation During Self-organization in Three-Phase Eutectic Solidification

Eutectic growth is an interesting example for exploring the topic of pattern-formation in multi-phase systems, where the growth of the phases is coupled with the diffusive transport of one or more components in the melt. While in the case of binary alloys, the number of possibilities are limited (lamellae, rods, labyrinth etc.), their number rapidly increases with the number of components and phases. In this paper, we will investigate pattern formation during three-phase eutectic solidification using a state-of-the art phase-field method based on the grand-canonical density formulation. The major aim of the study is to highlight the role of two properties, which are the volume fraction of the solid phases and the solid-liquid interfacial energies, in the self-organization of the solid phases during directional growth. Thereafter, we will show representative phase-field simulations of a micro-structure in a real alloy (Ag-Al-Cu) using an asymmetric phase diagram as well as interfacial properties.

Investigation on the origin of exchange bias in epitaxial, oriented and polycrystalline Fe3O4 thin films

We observe exchange bias (EB) in a single magnetic film Fe3O4 at temperature T < 200 K. Irrespective of crystallographic orientations of grown Fe3O4; they exhibit similar nature of EB for (100) epitaxial, (111) oriented and polycrystalline Fe3O4 thin films. Growth induced defects such as anti-phase boundaries (APBs) in epitaxial Fe3O4 thin film is known to have an influence on the magnetic interaction. But, it is noticed that according to the common consensus of APBs alone cannot explain the origin of EB. If majority of APBs end up with mainly anti-ferromagnetic interactions across these boundaries together with the internal ordering modification in Fe3O4, then EB can emerge at low temperatures. Hence, we propose the idea of directional anti-ferromagnetic APB induced EB in Fe3O4 triggered by internal ordering for T <= 200 K. Similar arguments are extended to (111) oriented as well as polycrystalline Fe3O4 films where the grain boundaries can impart same consequence as that of APBs. (C) 2015 Author(s).

AN INVESTIGATION OF CUTTING MECHANISMS AND STRAIN FIELDS DURING ORTHOGONAL CUTTING IN CFRPs

Fiber-reinforced plastics (FRPs) are typically difficult to machine due to their highly heterogeneous and anisotropic nature and the presence of two phases (fiber and matrix) with vastly different strengths and stiffnesses. Typical machining damage mechanisms in FRPs include series of brittle fractures (especially for thermosets) due to shearing and cracking of matrix material, fiber pull-outs, burring, fuzzing, fiber-matrix debonding, etc. With the aim of understanding the influence of the pronounced heterogeneity and anisotropy observed in FRPs, ``Idealized'' Carbon FRP (I-CFRP) plates were prepared using epoxy resin with embedded equispaced tows of carbon fibers. Orthogonal cutting of these I-CFRPs was carried out, and the chip formation characteristics, cutting force signals and strain distributions obtained during machining were analyzed using the Digital Image Correlation (DIC) technique. In addition, the same procedure was repeated on Uni-Directional CFRPs (UD-CFRPs). Chip formation mechanisms in FRPs were found to depend on the depth of cut and fiber orientation with pure epoxy showing a pronounced ``size effect.'' Experimental results indicate that in-situ full field strain measurements from DIC coupled with force measurements using dynamometry provide an adequate measure of anisotropy and heterogeneity during orthogonal cutting.