A Study of Speed of the Boundary Element Method as applied to the Realtime Computational Simulation of Biological Organs

In this work, possibility of simulating biological organs in realtime using the Boundary Element Method (BEM) is investigated. Biological organs are assumed to follow linear elastostatic material behavior, and constant boundary element is the element type used. First, a Graphics Processing Unit (GPU) is used to speed up the BEM computations to achieve the realtime performance. Next, instead of the GPU, a computer cluster is used. Results indicate that BEM is fast enough to provide for realtime graphics if biological organs are assumed to follow linear elastostatic material behavior. Although the present work does not conduct any simulation using nonlinear material models, results from using the linear elastostatic material model imply that it would be difficult to obtain realtime performance if highly nonlinear material models that properly characterize biological organs are used. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature.

A Kushner-Stratonovich Monte Carlo filter applied to nonlinear dynamical system identification

A Monte Carlo filter, based on the idea of averaging over characteristics and fashioned after a particle-based time-discretized approximation to the Kushner-Stratonovich (KS) nonlinear filtering equation, is proposed. A key aspect of the new filter is the gain-like additive update, designed to approximate the innovation integral in the KS equation and implemented through an annealing-type iterative procedure, which is aimed at rendering the innovation (observation prediction mismatch) for a given time-step to a zero-mean Brownian increment corresponding to the measurement noise. This may be contrasted with the weight-based multiplicative updates in most particle filters that are known to precipitate the numerical problem of weight collapse within a finite-ensemble setting. A study to estimate the a-priori error bounds in the proposed scheme is undertaken. The numerical evidence, presently gathered from the assessed performance of the proposed and a few other competing filters on a class of nonlinear dynamic system identification and target tracking problems, is suggestive of the remarkably improved convergence and accuracy of the new filter. (C) 2013 Elsevier B.V. All rights reserved.

Elliptical tracers in two-dimensional, homogeneous, isotropic fluid turbulence: The statistics of alignment, rotation, and nematic order

We study the statistical properties of orientation and rotation dynamics of elliptical tracer particles in two-dimensional, homogeneous, and isotropic turbulence by direct numerical simulations. We consider both the cases in which the turbulent flow is generated by forcing at large and intermediate length scales. We show that the two cases are qualitatively different. For large-scale forcing, the spatial distribution of particle orientations forms large-scale structures, which are absent for intermediate-scale forcing. The alignment with the local directions of the flow is much weaker in the latter case than in the former. For intermediate-scale forcing, the statistics of rotation rates depends weakly on the Reynolds number and on the aspect ratio of particles. In contrast with what is observed in three-dimensional turbulence, in two dimensions the mean-square rotation rate increases as the aspect ratio increases.

Flame surface statistics of constant-pressure turbulent expanding premixed flames

In this paper we investigate the local flame surface statistics of constant-pressure turbulent expanding flames. First the statistics of local length ratio is experimentally determined from high-speed planar Mie scattering images of spherically expanding flames, with the length ratio on the measurement plane, at predefined equiangular sectors, defined as the ratio of the actual flame length to the length of a circular-arc of radius equal to the average radius of the flame. Assuming isotropic distribution of such flame segments we then convolute suitable forms of the length-ratio probability distribution functions (pdfs) to arrive at the corresponding area-ratio pdfs. It is found that both the length ratio and area ratio pdfs are near log-normally distributed and shows self-similar behavior with increasing radius. Near log-normality and rather intermittent behavior of the flame-length ratio suggests similarity with dissipation rate quantities which stimulates multifractal analysis. (C) 2014 AIP Publishing LLC.

Statistics of the inverse-cascade regime in two-dimensional magnetohydrodynamic turbulence

We present a detailed direct numerical simulation of statistically steady, homogeneous, isotropic, two-dimensional magnetohydrodynamic turbulence. Our study concentrates on the inverse cascade of the magnetic vector potential. We examine the dependence of the statistical properties of such turbulence on dissipation and friction coefficients. We extend earlier work significantly by calculating fluid and magnetic spectra, probability distribution functions (PDFs) of the velocity, magnetic, vorticity, current, stream-function, and magnetic-vector-potential fields, and their increments. We quantify the deviations of these PDFs from Gaussian ones by computing their flatnesses and hyperflatnesses. We also present PDFs of the Okubo-Weiss parameter, which distinguishes between vortical and extensional flow regions, and its magnetic analog. We show that the hyperflatnesses of PDFs of the increments of the stream function and the magnetic vector potential exhibit significant scale dependence and we examine the implication of this for the multiscaling of structure functions. We compare our results with those of earlier studies.

BELIEF PROPAGATION FOR OPTIMAL EDGE COVER IN THE RANDOM COMPLETE GRAPH

We apply the objective method of Aldous to the problem of finding the minimum-cost edge cover of the complete graph with random independent and identically distributed edge costs. The limit, as the number of vertices goes to infinity, of the expected minimum cost for this problem is known via a combinatorial approach of Hessler and Wastlund. We provide a proof of this result using the machinery of the objective method and local weak convergence, which was used to prove the (2) limit of the random assignment problem. A proof via the objective method is useful because it provides us with more information on the nature of the edge's incident on a typical root in the minimum-cost edge cover. We further show that a belief propagation algorithm converges asymptotically to the optimal solution. This can be applied in a computational linguistics problem of semantic projection. The belief propagation algorithm yields a near optimal solution with lesser complexity than the known best algorithms designed for optimality in worst-case settings.

Probing a spin-glass state in SrRuO3 thin films through higher-order statistics of resistance fluctuations

The complex perovskite oxide SrRuO3 shows intriguing transport properties at low temperatures due to the interplay of spin, charge, and orbital degrees of freedom. One of the open questions in this system is regarding the origin and nature of the low-temperature glassy state. In this paper we report on measurements of higher-order statistics of resistance fluctuations performed in epitaxial thin films of SrRuO3 to probe this issue. We observe large low-frequency non-Gaussian resistance fluctuations over a certain temperature range. Our observations are compatible with that of a spin-glass system with properties described by hierarchical dynamics rather than with that of a simple ferromagnet with a large coercivity.

Role of applied potentials on bioleaching of chalcopyrite concentrate and growth of Acidithiobacillus ferrooxidans

The effect of applied DC potentials on the bioleaching of a chalcopyrite concentrate in the presence of Acidithiobacillus ferrooxidans is discussed. Copper dissolution was the highest at an applied potential of +600mV (SCE), while all the dissolved copper got cathodically deposited at a negative potential of -600mV (SCE). Electrobioleaching at an applied potential of +600mV (SCE) was established at different pulp densities as a function of time. The effect of applied potentials and electrolytic currents on the activity and growth of bacterial cells was assessed Preadaptation of bacterial cells to the concentrate slurry and electrolytic growth conditions significantly enhanced copper dissolution. Electrochemical and biochemical mechanisms involved in electrobioleaching are illustrated with respect to oxidative dissolution and biocatalysis of anodic oxidation.

Colored polydimethylsiloxane micropillar arrays for high throughput measurements of forces applied by genetic model organisms

Measuring forces applied by multi-cellular organisms is valuable in investigating biomechanics of their locomotion. Several technologies have been developed to measure such forces, for example, strain gauges, micro-machined sensors, and calibrated cantilevers. We introduce an innovative combination of techniques as a high throughput screening tool to assess forces applied by multiple genetic model organisms. First, we fabricated colored Polydimethylsiloxane (PDMS) micropillars where the color enhances contrast making it easier to detect and track pillar displacement driven by the organism. Second, we developed a semiautomated graphical user interface to analyze the images for pillar displacement, thus reducing the analysis time for each animal to minutes. The addition of color reduced the Young's modulus of PDMS. Therefore, the dye-PDMS composite was characterized using Yeoh's hyperelastic model and the pillars were calibrated using a silicon based force sensor. We used our device to measure forces exerted by wild type and mutant Caenorhabditis elegans moving on an agarose surface. Wild type C. elegans exert an average force of similar to 1 mu N on an individual pillar and a total average force of similar to 7.68 mu N. We show that the middle of C. elegans exerts more force than its extremities. We find that C. elegans mutants with defective body wall muscles apply significantly lower force on individual pillars, while mutants defective in sensing externally applied mechanical forces still apply the same average force per pillar compared to wild type animals. Average forces applied per pillar are independent of the length, diameter, or cuticle stiffness of the animal. We also used the device to measure, for the first time, forces applied by Drosophila melanogaster larvae. Peristaltic waves occurred at 0.4Hz applying an average force of similar to 1.58 mu N on a single pillar. Our colored microfluidic device along with its displacement tracking software allows us to measure forces applied by multiple model organisms that crawl or slither to travel through their environment. (C) 2015 AIP Publishing LLC.

Stability of the flow in a soft tube deformed due to an applied pressure gradient

A linear stability analysis is carried out for the flow through a tube with a soft wall in order to resolve the discrepancy of a factor of 10 for the transition Reynolds number between theoretical predictions in a cylindrical tube and the experiments of Verma and Kumaran J. Fluid Mech. 705, 322 (2012)]. Here the effect of tube deformation (due to the applied pressure difference) on the mean velocity profile and pressure gradient is incorporated in the stability analysis. The tube geometry and dimensions are reconstructed from experimental images, where it is found that there is an expansion and then a contraction of the tube in the streamwise direction. The mean velocity profiles at different downstream locations and the pressure gradient, determined using computational fluid dynamics, are found to be substantially modified by the tube deformation. The velocity profiles are then used in a linear stability analysis, where the growth rates of perturbations are calculated for the flow through a tube with the wall modeled as a neo-Hookean elastic solid. The linear stability analysis is carried out for the mean velocity profiles at different downstream locations using the parallel flow approximation. The analysis indicates that the flow first becomes unstable in the downstream converging section of the tube where the flow profile is more pluglike when compared to the parabolic flow in a cylindrical tube. The flow is stable in the upstream diverging section where the deformation is maximum. The prediction for the transition Reynolds number is in good agreement with experiments, indicating that the downstream tube convergence and the consequent modification in the mean velocity profile and pressure gradient could reduce the transition Reynolds number by an order of magnitude.

VAM Applied to Dimensional Reduction of Nonlinear Multifunctional Film-fabric Laminates

This work aims at asymptotically accurate dimensional reduction of non-linear multi-functional film-fabric laminates having specific application in design of envelopes for High Altitude Airships (HAA). The film-fabric laminate for airship envelope consists of a woven fabric core coated with thin films on each face. These films provide UV protection and Helium leakage prevention, while the core provides required structural strength. This problem is both geometrically and materially non-linear. To incorporate the geometric non-linearity, generalized warping functions are used and finite deformations are allowed. The material non-linearity is handled by using hyper-elastic material models for each layer. The development begins with three-dimensional (3-D) nonlinear elasticity and mathematically splits the analysis into a one-dimensional through-the-thickness analysis and a two-dimensional (2-D) plate analysis. The through-the-thickness analysis provides the 2-D constitutive law which is then given as an input to the 2-D reference surface analysis. The dimensional reduction is carried out using Variational Asymptotic Method (VAM) for moderate strains and very small thickness-to-wavelength ratio. It features the identification and utilization of additional small parameters such as ratio of thicknesses and stiffness coefficients of core and films. Closed form analytical expressions for warping functions and 2-D constitutive law of the film-fabric laminate are obtained.

A thermodynamic correlation between damage and fracture as applied to concrete fatigue

Fatigue damage in concrete is characterized by the simultaneous presence of micro and macrocracks. The theory of fracture mechanics conveniently handles the propagation of macrocracks, whereas damage mechanics precisely describes the state of microcracking. This paper provides a platform to correlate fracture mechanics and damage mechanics theories through an energy equivalence within a thermodynamic framework by equating the energy dissipated according to each theory. Through this correlation, damage corresponding to a given crack length could be obtained, and alternatively a discrete crack could be transformed into an equivalent damage zone. The results are validated using available experimental data on concrete fatigue including stiffness degradation and acoustic emission. (C) 2015 Elsevier Ltd. All rights reserved.

A thermodynamic correlation between damage and fracture as applied to concrete fatigue

Fatigue damage in concrete is characterized by the simultaneous presence of micro and macrocracics. The theory of fracture mechanics conveniently handles the propagation of macrocracks, whereas damage mechanics precisely describes the state of microcracking. This paper provides a platform to correlate fracture mechanics and damage mechanics theories through an energy equivalence within a thermodynamic framework by equating the energy dissipated according to each theory. Through this correlation, damage corresponding to a given crack length could be obtained, and alternatively a discrete crack could be transformed into an equivalent damage zone. The results are validated using available experimental data on concrete fatigue including stiffness degradation and acoustic emission. (C) 2015 Elsevier Ltd. All rights reserved.