A study on the usefulness of Support Vector Machines for the realtime computational simulation of soft biological organs

Realistic and realtime computational simulation of soft biological organs (e.g., liver, kidney) is necessary when one tries to build a quality surgical simulator that can simulate surgical procedures involving these organs. Since the realistic simulation of these soft biological organs should account for both nonlinear material behavior and large deformation, achieving realistic simulations in realtime using continuum mechanics based numerical techniques necessitates the use of a supercomputer or a high end computer cluster which are costly. Hence there is a need to employ soft computing techniques like Support Vector Machines (SVMs) which can do function approximation, and hence could achieve physically realistic simulations in realtime by making use of just a desktop computer. Present work tries to simulate a pig liver in realtime. Liver is assumed to be homogeneous, isotropic, and hyperelastic. Hyperelastic material constants are taken from the literature. An SVM is employed to achieve realistic simulations in realtime, using just a desktop computer. The code for the SVM is obtained from [1]. The SVM is trained using the dataset generated by performing hyperelastic analyses on the liver geometry, using the commercial finite element software package ANSYS. The methodology followed in the present work closely follows the one followed in [2] except that [2] uses Artificial Neural Networks (ANNs) while the present work uses SVMs to achieve realistic simulations in realtime. Results indicate the speed and accuracy that is obtained by employing the SVM for the targeted realistic and realtime simulation of the liver.

Theoretical description of coherent doublon creation via lattice modulation spectroscopy

Using a recently developed strong-coupling method, we present a comprehensive theory for doublon production processes in modulation spectroscopy of a three-dimensional system of ultracold fermionic atoms in an optical lattice with a trap. The theoretical predictions compare well to the experimental time traces of doublon production. For experimentally feasible conditions, we provide a quantitative prediction for the presence of a nonlinear ``two-photon'' excitation at strong modulation amplitudes.

Evaluating FuncSION: A software for automated synthesis of design solutions for stimulating ideation during mechanical conceptual design

The goal of the work reported in this paper is to use automated, combinatorial synthesis to generate alternative solutions to be used as stimuli by designers for ideation. FuncSION, a computational synthesis tool that can automatically synthesize solution concepts for mechanical devices by combining building blocks from a library, is used for this purpose. The objectives of FuncSION are to help generate a variety of functional requirements for a given problem and a variety of concepts to fulfill these functions. A distinctive feature of FuncSION is its focus on automated generation of spatial configurations, an aspect rarely addressed by other computational synthesis programs. This paper provides an overview of FuncSION in terms of representation of design problems, representation of building blocks, and rules with which building blocks are combined to generate concepts at three levels of abstraction: topological, spatial, and physical. The paper then provides a detailed account of evaluating FuncSION for its effectiveness in providing stimuli for enhanced ideation.

Mining Unit Tests for Discovery and Migration of Math APIs

Today's programming languages are supported by powerful third-party APIs. For a given application domain, it is common to have many competing APIs that provide similar functionality. Programmer productivity therefore depends heavily on the programmer's ability to discover suitable APIs both during an initial coding phase, as well as during software maintenance. The aim of this work is to support the discovery and migration of math APIs. Math APIs are at the heart of many application domains ranging from machine learning to scientific computations. Our approach, called MATHFINDER, combines executable specifications of mathematical computations with unit tests (operational specifications) of API methods. Given a math expression, MATHFINDER synthesizes pseudo-code comprised of API methods to compute the expression by mining unit tests of the API methods. We present a sequential version of our unit test mining algorithm and also design a more scalable data-parallel version. We perform extensive evaluation of MATHFINDER (1) for API discovery, where math algorithms are to be implemented from scratch and (2) for API migration, where client programs utilizing a math API are to be migrated to another API. We evaluated the precision and recall of MATHFINDER on a diverse collection of math expressions, culled from algorithms used in a wide range of application areas such as control systems and structural dynamics. In a user study to evaluate the productivity gains obtained by using MATHFINDER for API discovery, the programmers who used MATHFINDER finished their programming tasks twice as fast as their counterparts who used the usual techniques like web and code search, IDE code completion, and manual inspection of library documentation. For the problem of API migration, as a case study, we used MATHFINDER to migrate Weka, a popular machine learning library. Overall, our evaluation shows that MATHFINDER is easy to use, provides highly precise results across several math APIs and application domains even with a small number of unit tests per method, and scales to large collections of unit tests.

A general methodology for calculating mixed mode stress intensity factors and fracture toughness of solder joints with interfacial cracks

Solder joints in electronic packages undergo thermo-mechanical cycling, resulting in nucleation of micro-cracks, especially at the solder/bond-pad interface, which may lead to fracture of the joints. The fracture toughness of a solder joint depends on material properties, process conditions and service history, as well as strain rate and mode-mixity. This paper reports on a methodology for determining the mixed-mode fracture toughness of solder joints with an interfacial starter-crack, using a modified compact mixed mode (CMM) specimen containing an adhesive joint. Expressions for stress intensity factor (K) and strain energy release rate (G) are developed, using a combination of experiments and finite element (FE) analysis. In this methodology, crack length dependent geometry factors to convert for the modified CMM sample are first obtained via the crack-tip opening displacement (CTOD)-based linear extrapolation method to calculate the under far-field mode I and II conditions (f(1a) and f(2a)), (ii) generation of a master-plot to determine a(c), and (iii) computation of K and G to analyze the fracture behavior of joints. The developed methodology was verified using J-integral calculations, and was also used to calculate experimental fracture toughness values of a few lead-free solder-Cu joints. (C) 2014 Elsevier Ltd. All rights reserved.

Using formal reasoning on a model of tasks for FreeRTOS

FreeRTOS is an open-source real-time microkernel that has a wide community of users. We present the formal specification of the behaviour of the task part of FreeRTOS that deals with the creation, management, and scheduling of tasks using priority-based preemption. Our model is written in the Z notation, and we verify its consistency using the Z/Eves theorem prover. This includes a precise statement of the preconditions for all API commands. This task model forms the basis for three dimensions of further work: (a) the modelling of the rest of the behaviour of queues, time, mutex, and interrupts in FreeRTOS; (b) refinement of the models to code to produce a verified implementation; and (c) extension of the behaviour of FreeRTOS to multi-core architectures. We propose all three dimensions as benchmark challenge problems for Hoare's Verified Software Initiative.

Improved design methodology for the development of electrically actuated MEMS structures

A comprehensive design flow is proposed for the design of Micro Electro Mechanical Systems that are fabricated using SOIMUMPs process. Many of the designers typically do not model the temperature dependency of electrical conductivity, thermal conductivity and convection coefficient, as it is very cumbersome to create/incorporate the same in the existing FEM simulators. Capturing these dependencies is very critical particularly for structures that are electrically actuated. Lookup tables that capture the temperature dependency of electrical conductivity, thermal conductivity and convection coefficient are created. These look up tables are taken as inputs for a commercially available FEM simulator to model the semiconductor behavior. It is demonstrated that when temperature dependency for all the above mentioned parameters is not captured, then the error in estimation of the maximum temperature (for a given structure) could be as high as 30%. Error in the estimated resistance value under the same conditions is as high as 40%. When temperature dependency of the above mentioned parameters is considered then error w.r.t the measured values is less than 5%. It is evident that error in temperature estimates leads to erroneous results from mechanical simulations as well.

Key challenges for the creation and maintenance of specialist protein resources

As the volume of data relating to proteins increases, researchers rely more and more on the analysis of published data, thus increasing the importance of good access to these data that vary from the supplemental material of individual articles, all the way to major reference databases with professional staff and long-term funding. Specialist protein resources fill an important middle ground, providing interactive web interfaces to their databases for a focused topic or family of proteins, using specialized approaches that are not feasible in the major reference databases. Many are labors of love, run by a single lab with little or no dedicated funding and there are many challenges to building and maintaining them. This perspective arose from a meeting of several specialist protein resources and major reference databases held at the Wellcome Trust Genome Campus (Cambridge, UK) on August 11 and 12, 2014. During this meeting some common key challenges involved in creating and maintaining such resources were discussed, along with various approaches to address them. In laying out these challenges, we aim to inform users about how these issues impact our resources and illustrate ways in which our working together could enhance their accuracy, currency, and overall value. Proteins 2015; 83:1005-1013. (c) 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Design of symmetric TIM barrel proteins from first principles

Background: Computational protein design is a rapidly maturing field within structural biology, with the goal of designing proteins with custom structures and functions. Such proteins could find widespread medical and industrial applications. Here, we have adapted algorithms from the Rosetta software suite to design much larger proteins, based on ideal geometric and topological criteria. Furthermore, we have developed techniques to incorporate symmetry into designed structures. For our first design attempt, we targeted the (alpha/beta)(8) TIM barrel scaffold. We gained novel insights into TIM barrel folding mechanisms from studying natural TIM barrel structures, and from analyzing previous TIM barrel design attempts. Methods: Computational protein design and analysis was performed using the Rosetta software suite and custom scripts. Genes encoding all designed proteins were synthesized and cloned on the pET20-b vector. Standard circular dichroism and gel chromatographic experiments were performed to determine protein biophysical characteristics. 1D NMR and 2D HSQC experiments were performed to determine protein structural characteristics. Results: Extensive protein design simulations coupled with ab initio modeling yielded several all-atom models of ideal, 4-fold symmetric TIM barrels. Four such models were experimentally characterized. The best designed structure (Symmetrin-1) contained a polar, histidine-rich pore, forming an extensive hydrogen bonding network. Symmetrin-1 was easily expressed and readily soluble. It showed circular dichroism spectra characteristic of well-folded alpha/beta proteins. Temperature melting experiments revealed cooperative and reversible unfolding, with a T-m of 44 degrees C and a Gibbs free energy of unfolding (Delta G degrees) of 8.0 kJ/mol. Urea denaturing experiments confirmed these observations, revealing a C-m of 1.6 M and a Delta G degrees of 8.3 kJ/mol. Symmetrin-1 adopted a monomeric conformation, with an apparent molecular weight of 32.12 kDa, and displayed well resolved 1D-NMR spectra. However, the HSQC spectrum revealed somewhat molten characteristics. Conclusions: Despite the detection of molten characteristics, the creation of a soluble, cooperatively folding protein represents an advancement over previous attempts at TIM barrel design. Strategies to further improve Symmetrin-1 are elaborated. Our techniques may be used to create other large, internally symmetric proteins.

Synthesis of Instruction Extensions on HyperCell, a Reconfigurable Datapath

In this paper we present HyperCell as a reconfigurable datapath for Instruction Extensions (IEs). HyperCell comprises an array of compute units laid over a switch network. We present an IE synthesis methodology that enables post-silicon realization of IE datapaths on HyperCell. The synthesis methodology optimally exploits hardware resources in HyperCell to enable software pipelined execution of IEs. Exploitation of temporal reuse of data in HyperCell results in significant reduction of input/output bandwidth requirements of HyperCell.

Experimental studies on evaporation of fuel droplets under forced convection using spray in crossflow methodology

Experimental data on evaporation of droplets of decane, Jet-A1, and Jet-A1 surrogate are generated using a spray in crossflow configuration. The advantage of a crossflow configuration is that it enables us to study droplet evaporation under forced convective conditions involving droplet diameters of size relevant in practical combustors. Specifically, spray from an airblast atomizer is injected into a preheated crossflow of air and the resulting spray is characterized in terms of spray structure along with droplet size and velocity. An existing correlation for the spray trajectory is modified to incorporate the effect of elevated temperature, and is found to be in good agreement with the experimental data. Droplet sizes and velocities are measured at different locations along the crossflow direction to assess droplet evaporation. Specifically, droplets having size less than 25-mu m are selected for further analysis since these droplets are observed to exhibit velocities which are aligned with the crossflow. By comparing the droplet diameter profiles at upstream and downstream locations, the evaporation constant k for the d(2)-law is obtained iteratively. To assess the efficacy of the values of k obtained, the calculated droplet size distribution using the proposed k values at the downstream location is compared with the measured droplet size distribution at that location. A reasonably good match is found for all the three liquids confirming the validity of the analysis. (C) 2015 Elsevier Ltd. All rights reserved.

Plasmonic Lipid Bilayer Membranes for Enhanced Detection Sensitivity of Biolabeling Fluorophores

Plasmonics based sensing, using the surface plasmon resonance of metal nanoparticles, has been effectively demonstrated in various applications. Extending this methodology to cell and artificial lipid bilayer membranes is extremely beneficial in enhancing the sensitivity of the detection of binding and cellular transport of molecules across such membranes. Here, the creation of an artificial plasmonic biomembrane template is demonstrated and used to show the enhanced detection sensitivity of certain widely used biomarker molecules. The efficacy of these templates is explained in terms of the ability of the hydrophobic polymer grafted gold nanoparticles used to organize, penetrate, and fluidize the membranes. The enhancement of photoluminescence of the dye molecules used occurs over a reasonably large spectral range as compared to the plasmon resonance of gold nanoparticles. The results could, possibly, be extended to cellular membranes with relevant modifications, as well as to the detection of any other biological molecule appropriately labeled with fluorescent dye molecules, and demonstrate the versatility of these plasmonic bioinspired platforms as potential biochemical sensors.

Three-dimensional warping in strip-based composite four-bar mechanisms

This work intends to demonstrate the effect of geometrically non-linear cross-sectional analysis of certain composite beam-based four-bar mechanisms in predicting the three-dimensional warping of the cross-section. The only restriction in the present analysis is that the strains within each elastic body remain small (i.e., this work does not deal with materials exhibiting non-linear constitutive laws at the 3-D level). Here, all component bars of the mechanism are made of fiber-reinforced laminates. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. Each component of the mechanism is modeled as a beam based on geometrically non-linear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and non-linear 1-D analyses along the three beam reference curves. The splitting of the three-dimensional beam problem into two- and one-dimensional parts, called dimensional reduction, results in a tremendous savings of computational effort relative to the cost of three-dimensional finite element analysis, the only alternative for realistic beams. The analysis of beam-like structures made of laminated composite materials requires a much more complicated methodology. Hence, the analysis procedure based on Variational Asymptotic Method (VAM), a tool to carry out the dimensional reduction, is used here. The representative cross-sections of all component bars are analyzed using two different approaches: (1) Numerical Model and (2) Analytical Model. Four-bar mechanisms are analyzed using the above two approaches for Omega = 20 rad/s and Omega = pi rad/s and observed the same behavior in both cases. The noticeable snap-shots of the deformation shapes of the mechanism about 1000 frames are also reported using commercial software (I-DEAS + NASTRAN + ADAMS). The maximum out-of-plane warping of the cross-section is observed at the mid-span of bar-1, bar-2 and bar-3 are 1.5 mm, 250 mm and 1.0 mm, respectively, for t = 0:5 s. (C) 2015 Elsevier Ltd. All rights reserved.

Numerical Investigations for Leakage and Windage Heating in Straight-Through Labyrinth Seals

The ability to quantify leakage flow and windage heating for labyrinth seals with honeycomb lands is critical in understanding gas turbine engine system performance and predicting its component life. Variety of labyrinth seal configurations (number of teeth, stepped or straight, honeycomb cell size) are in use in gas turbines, and for each configuration, there are many geometric factors that can impact a seal's leakage and windage characteristics. This paper describes the development of a numerical methodology aimed at studying the effect of honeycomb lands on leakage and windage heating. Specifically, a three-dimensional computational fluid dynamics (CFD) model is developed utilizing commercial finite volume-based software incorporating the renormalization group (RNG) k-epsilon turbulence model with modified Schmidt number. The modified turbulence model is benchmarked and fine-tuned based on several experiments. Using this model, a broad parametric study is conducted by varying honeycomb cell size, pressure ratio (PR), and radial clearance for a four-tooth straight-through labyrinth seal. The results show good agreement with available experimental data. They further indicate that larger honeycomb cells predict higher seal leakage and windage heating at tighter clearances compared to smaller honeycomb cells and smooth lands. However, at open seal clearances larger honeycomb cells have lower leakage compared to smaller honeycomb cells.