Beyond classical dynamic structural plasticity using mesh-free modelling techniques

The problem of modelling the transient response of an elastic-perfectly-plastic cantilever beam, carrying an impulsively loaded tip mass, is,often referred to as the Parkes cantilever problem 25]; The permanent deformation of a cantilever struck transversely at its tip, Proc. R. Soc. A., 288, pp. 462). This paradigm for classical modelling of projectile impact on structures is re-visited and updated using the mesh-free method, smoothed particle hydrodynamics (SPH). The purpose of this study is to investigate further the behaviour of cantilever beams subjected to projectile impact at its tip, by considering especially physically real effects such as plastic shearing close to the projectile, shear deformation, and the variation of the shear strain along the length and across the thickness of the beam. Finally, going beyond macroscopic structural plasticity, a strategy to incorporate physical discontinuity (due to crack formation) in SPH discretization is discussed and explored in the context of tip-severance of the cantilever beam. Consequently, the proposed scheme illustrates the potency for a more refined treatment of penetration mechanics, paramount in the exploration of structural response under ballistic loading. The objective is to contribute to formulating a computational modelling framework within which transient dynamic plasticity and even penetration/failure phenomena for a range of materials, structures and impact conditions can be explored ab initio, this being essential for arriving at suitable tools for the design of armour systems. (C) 2014 Elsevier Ltd. All rights reserved.

An open source massively parallel solver for Richards equation: Mechanistic modelling of water fluxes at the watershed scale

In this paper we present a massively parallel open source solver for Richards equation, named the RichardsFOAM solver. This solver has been developed in the framework of the open source generalist computational fluid dynamics tool box OpenFOAM (R) and is capable to deal with large scale problems in both space and time. The source code for RichardsFOAM may be downloaded from the CPC program library website. It exhibits good parallel performances (up to similar to 90% parallel efficiency with 1024 processors both in strong and weak scaling), and the conditions required for obtaining such performances are analysed and discussed. These performances enable the mechanistic modelling of water fluxes at the scale of experimental watersheds (up to few square kilometres of surface area), and on time scales of decades to a century. Such a solver can be useful in various applications, such as environmental engineering for long term transport of pollutants in soils, water engineering for assessing the impact of land settlement on water resources, or in the study of weathering processes on the watersheds. (C) 2014 Elsevier B.V. All rights reserved.

Modelling of transient heat conduction with diffuse interface methods

We present a survey on different numerical interpolation schemes used for two-phase transient heat conduction problems in the context of interface capturing phase-field methods. Examples are general transport problems in the context of diffuse interface methods with a non-equal heat conductivity in normal and tangential directions to the interface. We extend the tonsorial approach recently published by Nicoli M et al (2011 Phys. Rev. E 84 1-6) to the general three-dimensional (3D) transient evolution equations. Validations for one-dimensional, two-dimensional and 3D transient test cases are provided, and the results are in good agreement with analytical and numerical reference solutions.

Strategy for Enhancing the Dielectric Constant of Organic Semiconductors Without Sacrificing Charge Carrier Mobility and Solubility

Current organic semiconductors for organic photovoltaics (OPV) have relative dielectric constants (relative permittivities, epsilon(r)) in the range of 2-4. As a consequence, Coulombically bound electron-hole pairs (excitons) are produced upon absorption of light, giving rise to limited power conversion efficiencies. We introduce a strategy to enhance epsilon(r) of well-known donors and acceptors without breaking conjugation, degrading charge carrier mobility or altering the transport gap. The ability of ethylene glycol (EG) repeating units to rapidly reorient their dipoles with the charge redistributions in the environment was proven via density functional theory (DFT) calculations. Fullerene derivatives functionalized with triethylene glycol side chains were studied for the enhancement of epsilon(r) together with poly(p-phenylene vinylene) and diketo-pyrrolopyrrole based polymers functionalized with similar side chains. The polymers showed a doubling of epsilon(r) with respect to their reference polymers in identical backbone. Fullerene derivatives presented enhancements up to 6 compared with phenyl-C-61-butyric acid methyl ester (PCBM) as the reference. Importantly, the applied modifications did not affect the mobility of electrons and holes and provided excellent solubility in common organic solvents.

Identifying regions for conservation of sloth bears through occupancy modelling in north-eastern Karnataka, India

In the absence of information on species in decline with contracting ranges, management should emphasize remaining populations and protection of their habitats. Threatened by anthropogenic pressure including habitat degradation and loss, sloth bears (Melursus ursinus) in India have become limited in range, habitat, and population size. We identified ecological and anthropogenic determinants of occurrence within an occupancy framework to evaluate habitat suitability of non-protected regions (with sloth bears) in northeastern Karnataka, India. We employed a systematic sampling methodology to yield presence absence data to examine a priori hypotheses of determinants that affected occupancy. These covariates were broadly classified as habitat or anthropogenic factors. Mean number of termite mounds and trees positively influenced sloth bear occupancy, and grazing pressure expounded by mean number of livestock dung affected it negatively. Also, mean percentage of shrub coverage had no impact on bear inhabitance. The best fitting model further predicted habitats in Bukkasagara, Agoli, and Benakal reserved forests to have 38%, 75%, and 88%, respectively, of their sampled grid cells with high occupancies (>0.70) albeit little or no legal protection. We recommend a conservation strategy that includes protection of vegetation stand-structure, maintenance of soil moisture, and enrichment of habitat for the long-term welfare of this species.

High Proton Mobility, Solvent Induced Single Crystal to Single Crystal Structural Transformation, and Related Studies on a Family of Compounds Formed from Mn-3 Oxo-Clusters

The reaction between 4,4'-sulfonyldibenzoic acid (H(2)SDBA) and manganese under mild conditions resulted in the isolation of two new three-dimensional compounds, Mn-4(C14H8O6S)(4)(DMA)(2)]center dot 3DMA, I, and Mn-3(C14H8O6S)(3)(DMA)(2)(MeOH)]center dot DMA, IIa. Both structures have Mn-3 trimer oxo cluster units. While the Mn-3 oxoclusters are connected through octahedral manganese forming one-dimensional Mn-O-Mn chains in I, the Mn-3 units are isolated in IIa. The SDBA units connect the Mn-O-Mn chains and the Mn-3 clusters giving rise to the three-dimensional structure. Both compounds have coordinated and free solvent molecules. In IIa, two different solvent molecules are coordinated, of which one solvent can be reversibly exchanged by a variety of other similar solvents via a solvent-mediated single crystal to single crystal (SCSC) transformation. The free lattice DMA solvent molecules in I can be exchanged by water molecules resulting in hydrophilic channels. Proton conductivity studies on I reveals a high proton mobility with conductivity values of similar to 0.87 x 10(-3) Omega(-1) cm(-1) at 34 degrees C and 98% RH, which is comparable to some of the good proton conductivity values observed in inorganic coordination polymers. We have also shown structural transformation of I to IIa through a possible dissolution and recrystallization pathway. In addition, both I and IIa appear to transform to two other manganese compounds H3O]Mn-3(mu(3)-OH)(C14H8O6S)(3)(H2O)](DMF)(5) and H3O](2)Mn-7(mu 3-OH)(4)(C14H8O6S)(6)(H2O)(4)](H2O)(2)(DMF)(8) under suitable reaction conditions. We have partially substituted Co in place of Mn in the Mn-3 trimer clusters forming CoMn2(C14H8O6S)(3)(DMA)(2)(EtOH)]center dot DMA, III, a structure that is closely related to IIa. All the compounds reveal antiferromagnetic behavior. On heating, the cobalt substituted phase (compound III) forms a CoMn2O4 spinel phase with particle sizes in the nanometer range.

Modelling of soldier fly halteres for gyroscopic oscillations

Nature has evolved a beautiful design for small-scale vibratory rategyro in the form of dipteran halteres that detect body rotations via Coriolis acceleration. In most Diptera, including soldier fly, Hermetia illucens, halteres are a pair of special organs, located in the space between the thorax and the abdomen. The halteres along with their connecting joint with the fly's body constitute a mechanism that is used for muscle-actuated oscillations of the halteres along the actuation direction. These oscillations lead to bending vibrations in the sensing direction (out of the haltere's actuation plane) upon any impressed rotation due to the resulting Coriolis force. This induced vibration is sensed by the sensory organs at the base of the haltere in order to determine the rate of rotation. In this study, we evaluate the boundary conditions and the stiffness of the anesthetized halteres along the actuation and the sensing direction. We take several cross-sectional SEM (scanning electron microscope) images of the soldier fly haltere and construct its three dimensional model to get the mass properties. Based on these measurements, we estimate the natural frequency along both actuation and sensing directions, propose a finite element model of the haltere's joint mechanism, and discuss the significance of the haltere's asymmetric cross-section. The estimated natural frequency along the actuation direction is within the range of the haltere's flapping frequency. However, the natural frequency along the sensing direction is roughly double the haltere's flapping frequency that provides a large bandwidth for sensing the rate of rotation to the soldier flies.

REDUCED ORDER MODELLING OF COMBUSTION INSTABILITY IN A BACKWARD FACING STEP COMBUSTOR

This paper develops a fully coupled time domain Reduced Order Modelling (ROM) approach to model unsteady combustion dynamics in a backward facing step combustor The acoustic field equations are projected onto the canonical acoustic eigenmodes of the systems to obtain a coupled system of modal evolution equations. The heat release response of the flame is modelled using the G-equation approach. Vortical velocity fluctuations that arise due to shear layer rollup downstream of the step are modelled using a simplified 1D-advection equation whose phase speed is determined from a linear, local, temporal stability analysis of the shear layer just downstream of the step. The hydrodynamic stability analysis reveals a abrupt change in the value of disturbance phase speed from unity for Re < Re-crit to 0.5 for Re > Re-crit, where Remit for the present geometry was found to be approximate to 10425. The results for self-excited flame response show highly wrinkled flame shapes that are qualitatively similar to those seen in prior experiments of acoustically forced flames. The effect of constructive and destructive interference between the two contributions to flame surface wrinkling results in high amplitude wrinkles for the case when K-c -> 1.

Room-temperature bandlike transport and Hall effect in a high-mobility ambipolar polymer

The advent of a new class of high-mobility semiconducting polymers opens up a window to address fundamental issues in electrical transport mechanism such as transport between localized states versus extended state conduction. Here, we investigate the origin of the ultralow degree of disorder (E-a similar to 16 meV) and the ``bandlike'' negative temperature (T) coefficient of the field effect electron mobility: mu(e)(FET) (T) in a high performance (mu(e)(FET) > 2.5 cm(2) V-1 s(-1)) diketopyrrolopyrrole based semiconducting polymer. Models based on the framework of mobility edge with exponential density of states are invoked to explain the trends in transport. The temperature window over which the system demonstrates delocalized transport was tuned by a systematic introduction of disorder at the transport interface. Additionally, the Hall mobility (mu(e)(Hall)) extracted from Hall voltage measurements in these devices was found to be comparable to field effect mobility (mu(e)(FET)) in the high T bandlike regime. Comprehensive studies with different combinations of dielectrics and semiconductors demonstrate the effectiveness of rationale molecular design, which emphasizes uniform-energetic landscape and low reorganization energy.

Hybrid Elements for Modelling Squeeze Film Effects Coupled with Structural Interactions in Vibratory MEMS Devices

We present a hybrid finite element based methodology to solve the coupled fluid structure problem of squeeze film effects in vibratory MEMS devices, such as gyroscopes, RF switches, and 2D resonators. The aforementioned devices often have a thin plate like structure vibrating normally to a fixed substrate, and are generally not perfectly vacuum packed. This results in a thin air film being trapped between the vibrating plate and the fixed substrate which behaves like a squeeze film offering both stiffness and damping. For accurate modelling of such devices the squeeze film effects must be incorporated. Extensive literature is available on squeeze film modelling, however only a few studies address the coupled fluid elasticity problem. The majority of the studies that account for the plate elasticity coupled with the fluid equation, either use approximate mode shapes for the plate or use iterative solution strategies. In an earlier work we presented a single step coupled methodology using only one type of displacement based element to solve the coupled problem. The displacement based finite element models suffer from locking issues when it comes to modelling very thin structures with the lateral dimensions much larger than the plate thickness as is typical in MEMS devices with squeeze film effects. In this work we present another coupled formulation where we have used hybrid elements to model the structural domain. The numerical results show a huge improvement in convergence and accuracy with coarse hybrid mesh as compared to displacement based formulations. We further compare our numerical results with experimental data from literature and find them to be in good accordance.

Modelling ternary effects on antiphase boundary energy of Ni3Al

The shearing of ordered gamma' precipitates by matrix dislocations results in the formation of antiphase boundaries (APB) in Ni-base superalloys. The APB energy is an important source of order-strengthening in disk and blade alloys where Ti and Ta substitute for Al in gamma'. While the importance of APB energy is well-acknowledged, the effect of alloying on APB energy is not fully understood. In the present study, the effect of Ti and Ta additions on the {111} and {010} APB energies was probed via electronic structure calculations. Results suggest that at low levels of Ti/Ta, APB energies on either plane increases with alloying. However, at higher Ti/Ta levels, the APB energies decrease with alloying. These trends understood by accounting for nearest neighbour violations about the APB and additionally, invoking the effect of precipitate composition on the energy penalty of the violations. We propose an Environment Dependent Nearest Neighbour Bond (EDNNB) model that predicts APB energies that are in close agreement to calculated values.

Modelling the influence of land-use changes on biophysical and biochemical interactions at regional and global scales

Land-use changes since the start of the industrial era account for nearly one-third of the cumulative anthropogenic CO2 emissions. In addition to the greenhouse effect of CO2 emissions, changes in land use also affect climate via changes in surface physical properties such as albedo, evapotranspiration and roughness length. Recent modelling studies suggest that these biophysical components may be comparable with biochemical effects. In regard to climate change, the effects of these two distinct processes may counterbalance one another both regionally and, possibly, globally. In this article, through hypothetical large-scale deforestation simulations using a global climate model, we contrast the implications of afforestation on ameliorating or enhancing anthropogenic contributions from previously converted (agricultural) land surfaces. Based on our review of past studies on this subject, we conclude that the sum of both biophysical and biochemical effects should be assessed when large-scale afforestation is used for countering global warming, and the net effect on global mean temperature change depends on the location of deforestation/afforestation. Further, although biochemical effects trigger global climate change, biophysical effects often cause strong local and regional climate change. The implication of the biophysical effects for adaptation and mitigation of climate change in agriculture and agroforestry sectors is discussed. center dot Land-use changes affect global and regional climates through both biochemical and biophysical process. center dot Climate effect from biophysical process depends on the location of land-use change. center dot Climate mitigation strategies such as afforestation/reforestation should consider the net effect of biochemical and biophysical processes for effective mitigation. center dot Climate-smart agriculture could use bio-geoengineering techniques that consider plant biophysical characteristics such as reflectivity and water use efficiency.

Modelling Proximity Effects in Transverse Magnetic Mode for Lightning Studies

Lightning strike to instrumented and communication towers can be a source of electromagnetic disturbance to the system connected. Long cables running on these towers can get significant induction to their sheath/core, which would then couple to the connected equipments. For a quantitative analysis of the situation, suitable theoretical analysis is necessary. Due to the dominance of the transverse magnetic mode during the fast rising portion of the stroke current, which is the period of significant induction, a full wave solution based on Maxwell's equations is necessary. Owing to the large geometric aspect ratio of tower lattice elements and for feasibility of a numerical solution, the thin-wire formulation for the electric field integral equation is generally adopted. However, the classical thin-wire formulation is not set for handling non-cylindrical conductors like tower lattice elements and the proximity of other conductors. The present work investigates further into a recently proposed method for handling such a situation and optimizes the numerical solution approach.

Integrating AlGaN/GaN high electron mobility transistor with Si: A comparative study of integration schemes

AlGaN/GaN high electron mobility transistor stacks deposited on a single growth platform are used to compare the most common transition, AlN to GaN, schemes used for integrating GaN with Si. The efficiency of these transitions based on linearly graded, step graded, interlayer, and superlattice schemes on dislocation density reduction, stress management, surface roughness, and eventually mobility of the 2D-gas are evaluated. In a 500 nm GaN probe layer deposited, all of these transitions result in total transmission electron microscopy measured dislocations densities of 1 to 3 x 10(9)/cm(2) and <1 nm surface roughness. The 2-D electron gas channels formed at an AlGaN-1 nm AlN/GaN interface deposited on this GaN probe layer all have mobilities of 1600-1900 cm(2)/V s at a carrier concentration of 0.7-0.9 x 10(13)/cm(2). Compressive stress and changes in composition in GaN rich regions of the AlN-GaN transition are the most effective at reducing dislocation density. Amongst all the transitions studied the step graded transition is the one that helps to implement this feature of GaN integration in the simplest and most consistent manner. (C) 2015 AIP Publishing LLC.