Markov chain splitting methods in structural reliability integral estimation

Monte Carlo simulation methods involving splitting of Markov chains have been used in evaluation of multi-fold integrals in different application areas. We examine in this paper the performance of these methods in the context of evaluation of reliability integrals from the point of view of characterizing the sampling fluctuations. The methods discussed include the Au-Beck subset simulation, Holmes-Diaconis-Ross method, and generalized splitting algorithm. A few improvisations based on first order reliability method are suggested to select algorithmic parameters of the latter two methods. The bias and sampling variance of the alternative estimators are discussed. Also, an approximation to the sampling distribution of some of these estimators is obtained. Illustrative examples involving component and series system reliability analyses are presented with a view to bring out the relative merits of alternative methods. (C) 2015 Elsevier Ltd. All rights reserved.

Simultaneous perturbation methods for adaptive labor staffing in service systems

We consider the problem of optimizing the workforce of a service system. Adapting the staffing levels in such systems is non-trivial due to large variations in workload and the large number of system parameters do not allow for a brute force search. Further, because these parameters change on a weekly basis, the optimization should not take longer than a few hours. Our aim is to find the optimum staffing levels from a discrete high-dimensional parameter set, that minimizes the long run average of the single-stage cost function, while adhering to the constraints relating to queue stability and service-level agreement (SLA) compliance. The single-stage cost function balances the conflicting objectives of utilizing workers better and attaining the target SLAs. We formulate this problem as a constrained parameterized Markov cost process parameterized by the (discrete) staffing levels. We propose novel simultaneous perturbation stochastic approximation (SPSA)-based algorithms for solving the above problem. The algorithms include both first-order as well as second-order methods and incorporate SPSA-based gradient/Hessian estimates for primal descent, while performing dual ascent for the Lagrange multipliers. Both algorithms are online and update the staffing levels in an incremental fashion. Further, they involve a certain generalized smooth projection operator, which is essential to project the continuous-valued worker parameter tuned by our algorithms onto the discrete set. The smoothness is necessary to ensure that the underlying transition dynamics of the constrained Markov cost process is itself smooth (as a function of the continuous-valued parameter): a critical requirement to prove the convergence of both algorithms. We validate our algorithms via performance simulations based on data from five real-life service systems. For the sake of comparison, we also implement a scatter search based algorithm using state-of-the-art optimization tool-kit OptQuest. From the experiments, we observe that both our algorithms converge empirically and consistently outperform OptQuest in most of the settings considered. This finding coupled with the computational advantage of our algorithms make them amenable for adaptive labor staffing in real-life service systems.

Texture evolution during hot deformation of Moly-TZM

Moly-TZM was deformed at constant strain rate of 1.0 s(-1) to investigate the high strain rate deformation behaviour by microstructural and stress response change within a temperature range of 1400-1700 degrees C. To correlate the deformation behaviour with orientational change, recrystallization and recovery of the material, the microstructural investigation was undertaken using scanning electron microscopy (SEM) of electron back scattered diffraction (EBSD). Depending on the grain size and orientation spread recrystallized grains were identified and texture was calculated. Change in grain boundary characteristics with increasing temperature was determined by the misorientation angle distribution for the deformed and recrystallized grains. Subgrain coalescence and increase in subgrain size with increasing temperature was observed, indicating recrystallization not only occurred from the nucleation of the dislocation free grains in grain boundaries but also from the subgrain rotation and merging of the subgrains by annihilation of the low angle grain boundaries. Detailed studies on the evolution of texture of recrystallized grains showed continuous increase in <001> fiber texture in recrystallised grains, in contrast to a mixed fiber <001> +<111> for the deformed grains.

Study of Multilevel Sinusoidal PWM methods for Cascaded H-Bridge Multilevel Inverters

This paper is a study of Multilevel Sinusoidal Pulse Width Modulation (MSPWM) methods; Phase Disposition (PD), Alternate Phase Opposition Disposition (APOD), Phase Opposition Disposition (POD) on a single phase Cascaded H-Bridge Multilevel inverter. Various factors such as amplitude modulation index (Ma), frequency modulation index (M-f), phase angle between carrier and reference modulating wave (phi) have been considered for simulation. Variation in these factors and their effect on inverter performance is evaluated. Factors such as DC bus utilization, output r.m.s voltage, total harmonic distortion (%THD), dominant harmonic order, switching losses are evaluated based on simulation results.

Role of deformation temperature on the evolution and heterogeneity of texture during equal channel angular pressing of magnesium

Investigations on texture evolution and through-thickness texture heterogeneity during equal channel angular pressing (ECAP) of pure magnesium at 200 degrees C, 150 degrees C and room temperature (RT) was carried out by neutron, high energy synchrotron X-ray and electron back-scatter diffraction. Irrespective of the ECAP temperature, a distinctive basal (B) and pyramidal (C-2) II type of fibers forms. The texture differs in the bottom 1 mm portion, where the B-fiber is shifted similar to 55 degrees due to negative shear attributed to friction. (C) 2015 Elsevier Inc. All rights reserved.

Smaller is Plastic: Polymorphic Structures and Mechanism of Deformation in Nanoscale hcp Metals

Using first-principles calculations, we establish the existence of highly-stable polymorphs of hcp metals (Ti, Mg, Be, La and Y) with nanoscale structural periodicity. They arise from heterogeneous deformation of the hcp structure occurring in response to large shear stresses localized at the basal planes separated by a few nanometers. Through Landau theoretical analysis, we show that their stability derives from nonlinear coupling between strains at different length scales. Such multiscale hyperelasticity and long-period structures constitute a new mechanism of size-dependent plasticity and its enhancement in nanoscale hcp metals.

A Remark on the A Posteriori Error Analysis of Discontinuous Galerkin Methods for the Obstacle Problem

We revisit the a posteriori error analysis of discontinuous Galerkin methods for the obstacle problem derived in 25]. Under a mild assumption on the trace of obstacle, we derive a reliable a posteriori error estimator which does not involve min/max functions. A key in this approach is an auxiliary problem with discrete obstacle. Applications to various discontinuous Galerkin finite element methods are presented. Numerical experiments show that the new estimator obtained in this article performs better.

Volume conservation issues in incompressible smoothed particle hydrodynamics

A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates. (C) 2015 Elsevier Inc. All rights reserved.

Surface motions and intraplate continental deformation in Alaska driven by mantle flow

The degree to which the lithosphere and mantle are coupled and contribute to surface deformation beneath continental regions remains a fundamental question in the field of geodynamics. Here we use a new approach with a surface deformation field constrained by GPS, geologic, and seismicity data, together with a lithospheric geodynamic model, to solve for tractions inferred to be generated by mantle convection that (1) drive extension within interior Alaska generating southward directed surface motions toward the southern convergent plate boundary, (2) result in accommodation of the relative motions between the Pacific and North America in a comparatively small zone near the plate boundary, and (3) generate the observed convergence within the North American plate interior in the Mackenzie mountains in northwestern Canada. The evidence for deeper mantle influence on surface deformation beneath a continental region suggests that this mechanism may be an important contributing driver to continental plate assemblage and breakup.

Deformation field evolution in indentation of a porous brittle solid

An experimental study of plane strain wedge indentation of a model porous brittle solid has been made to understand the effect of indentation parameters on the evolution of the deformation field and the accompanying volume change. A series of high-speed, high-resolution images of the indentation region and simultaneous measurements of load response were captured through the progression of the indentation process. Particle image velocimetry analysis of the images facilitated in situ measurement of the evolution of the resulting plastic zone in terms of incremental material displacement (velocity), strain rate, strain and volume change (e.g., local pore collapse). These measurements revealed initiation and propagation of flow localizations and fractures, as well as enabled estimate of volume changes occurring in the deformation zone. The results were directly compared with theoretical estimates of indentation pressure and deformation zone geometry and were used to validate a modified cavity expansion solution that incorporates effects of volume changes in the plastic zone. (C) 2015 Elsevier Ltd. All rights reserved.

De-DUFing the DUFs: Deciphering distant evolutionary relationships of Domains of Unknown Function using sensitive homology detection methods

Background: In the post-genomic era where sequences are being determined at a rapid rate, we are highly reliant on computational methods for their tentative biochemical characterization. The Pfam database currently contains 3,786 families corresponding to ``Domains of Unknown Function'' (DUF) or ``Uncharacterized Protein Family'' (UPF), of which 3,087 families have no reported three-dimensional structure, constituting almost one-fourth of the known protein families in search for both structure and function. Results: We applied a `computational structural genomics' approach using five state-of-the-art remote similarity detection methods to detect the relationship between uncharacterized DUFs and domain families of known structures. The association with a structural domain family could serve as a start point in elucidating the function of a DUF. Amongst these five methods, searches in SCOP-NrichD database have been applied for the first time. Predictions were classified into high, medium and low-confidence based on the consensus of results from various approaches and also annotated with enzyme and Gene ontology terms. 614 uncharacterized DUFs could be associated with a known structural domain, of which high confidence predictions, involving at least four methods, were made for 54 families. These structure-function relationships for the 614 DUF families can be accessed on-line at http://proline.biochem.iisc.ernet.in/RHD_DUFS/. For potential enzymes in this set, we assessed their compatibility with the associated fold and performed detailed structural and functional annotation by examining alignments and extent of conservation of functional residues. Detailed discussion is provided for interesting assignments for DUF3050, DUF1636, DUF1572, DUF2092 and DUF659. Conclusions: This study provides insights into the structure and potential function for nearly 20 % of the DUFs. Use of different computational approaches enables us to reliably recognize distant relationships, especially when they converge to a common assignment because the methods are often complementary. We observe that while pointers to the structural domain can offer the right clues to the function of a protein, recognition of its precise functional role is still `non-trivial' with many DUF domains conserving only some of the critical residues. It is not clear whether these are functional vestiges or instances involving alternate substrates and interacting partners. Reviewers: This article was reviewed by Drs Eugene Koonin, Frank Eisenhaber and Srikrishna Subramanian.

Effect of Particle Concentration on Shape Deformation and Secondary Atomization Characteristics of a Burning Nanotitania Dispersion Droplet

Secondary atomization characteristics of burning bicomponent (ethanol-water) droplets containing titania nanoparticles (NPs) in dilute (0.5% and 1 wt.%) and dense concentrations (5% and 7.5 wt.%) are studied experimentally at atmospheric pressure under normal gravity. It is observed that both types of nanofuel droplets undergo distinct modes of secondary breakup, which are primarily responsible for transporting particles from the droplet domain to the flame zone. For dilute nanosuspensions, disruptive response is characterized by low intensity atomization modes that cause small-scale localized flame distortion. In contrast, the disruption behavior at dense concentrations is governed by high intensity bubble ejections, which result in severe disruption of the flame envelope.

Joint Strength and Wall Deformation Characteristics of a Single-Cell Geocell Subjected to Uniaxial Compression

Geocells are three-dimensional expandable panels with a wide range of applications in geotechnical engineering. A geocell is made up of many internally connected single cells. The current study discusses the joint strength and the wall deformation characteristics of a single cell when it is subjected to uniaxial compression. The study helps to understand the causes for the failure of the single cell in a cellular confinement system. Experimental studies were conducted on single cells with cell pockets filled up with three different infill materials, namely silty clay, sand, and the aggregates. The results of the experimental study revealed that the deformation of the geocell wall decreases with the increase in the friction angle of the infill material. Experimental results were also validated using numerical simulations carried out using Lagrangian analysis software. The experiment and the numerical results were found to be in good agreement with each other. A simple analytical model based on the theory of thin cylinders is also proposed to calculate the accumulated strain of the geocell wall. This model operates under a simple elastic solution framework. The proposed model slightly overestimates the strains as compared with experimental and numerical values. (C) 2014 American Society of Civil Engineers.

Kinematic flow patterns in slow deformation of a dense granular material

The kinematic flow pattern in slow deformation of a model dense granular medium is studied at high resolution using in situ imaging, coupled with particle tracking. The deformation configuration is indentation by a flat punch under macroscopic plane-strain conditions. Using a general analysis method, velocity gradients and deformation fields are obtained from the disordered grain arrangement, enabling flow characteristics to be quantified. The key observations are the formation of a stagnation zone, as in dilute granular flow past obstacles; occurrence of vortices in the flow immediately underneath the punch; and formation of distinct shear bands adjoining the stagnation zone. The transient and steady state stagnation zone geometry, as well as the strength of the vortices and strain rates in the shear bands, are obtained from the experimental data. All of these results are well-reproduced in exact-scale non-smooth contact dynamics simulations. Full 3D numerical particle positions from the simulations allow extraction of flow features that are extremely difficult to obtain from experiments. Three examples of these, namely material free surface evolution, deformation of a grain column below the punch and resolution of velocities inside the primary shear band, are highlighted. The variety of flow features observed in this model problem also illustrates the difficulty involved in formulating a complete micromechanical analytical description of the deformation.

Influence of mode of deformation on microstructural heterogeneities in Ni subjected to large strain deformation

In the present work, the effect of deformation mode (uniaxial compression, rolling and torsion) on the microstructural heterogeneities in a commercial purity Ni is reported. For a given equivalent von Mises strain, samples subjected to torsion have shown higher fraction of high-angle boundaries, kernel average misorientation and recrystallization nuclei when compared to uniaxially compressed and rolled samples. This is attributed to the differences in the slip system activity under different modes of deformation.