VAM applied to dimensional reduction of non-linear hyperelastic plates

This work aims at dimensional reduction of non-linear isotropic hyperelastic plates in an asymptotically accurate manner. The problem is both geometrically and materially non-linear. The geometric non-linearity is handled by allowing for finite deformations and generalized warping while the material non-linearity is incorporated through hyperelastic material model. The development, based on the Variational Asymptotic Method (VAM) with moderate strains and very small thickness to shortest wavelength of the deformation along the plate reference surface as small parameters, begins with three-dimensional (3-D) non-linear elasticity and mathematically splits the analysis into a one-dimensional (1-D) through-the-thickness analysis and a two-dimensional (2-D) plate analysis. Major contributions of this paper are derivation of closed-form analytical expressions for warping functions and stiffness coefficients and a set of recovery relations to express approximately the 3-D displacement, strain and stress fields. Consistent with the 2-D non-linear constitutive laws, 2-D plate theory and corresponding finite element program have been developed. Validation of present theory is carried out with a standard test case and the results match well. Distributions of 3-D results are provided for another test case. (c) 2012 Elsevier Ltd. All rights reserved.

A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs

In this work, first a Fortran code is developed for three dimensional linear elastostatics using constant boundary elements; the code is based on a MATLAB code developed by the author earlier. Next, the code is parallelized using BLACS, MPI, and ScaLAPACK. Later, the parallelized code is used to demonstrate the usefulness of the Boundary Element Method (BEM) as applied to the realtime computational simulation of biological organs, while focusing on the speed and accuracy offered by BEM. A computer cluster is used in this part of the work. The commercial software package ANSYS is used to obtain the `exact' solution against which the solution from BEM is compared; analytical solutions, wherever available, are also used to establish the accuracy of BEM. A pig liver is the biological organ considered. Next, instead of the computer cluster, a Graphics Processing Unit (GPU) is used as the parallel hardware. Results indicate that BEM is an interesting choice for the simulation of biological organs. 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. Also, a serial MATLAB code, and both serial and parallel versions of a Fortran code, which can solve three dimensional (3D) linear elastostatic problems using constant boundary elements, are provided as supplementary files that can be freely downloaded.

A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs

In this work, possibility of simulating biological organs in realtime using the Boundary Element Method (BEM) is investigated, with specific reference to the speed and the accuracy offered by BEM. 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. A pig liver is the biological organ considered. Results indicate that BEM is an interesting choice for the simulation of biological organs. 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 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.

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.

Lifshitz transition and modulation of electronic and transport properties of bilayer graphene by sliding and applied normal compressive strain

Using density functional theory (DFT) we investigate the changes in electronic and transport properties of graphene bilayer caused by sliding one of the layers. Change in stacking pattern breaks the lattice symmetry, which results in Lifshitz transition together with the modulation of the electronic structure. Going from AA to AB stacking by sliding along armchair direction leads to a drastic transition in electronic structure from linear to parabolic dispersion. Our transport calculations show a significant change in the overall transmission value for large sliding distances along zigzag direction. The increase in interlayer coupling with normal compressive strain increases the overlapping of conduction and valence band, which leads to further shift in the Dirac points and an enhancement in the Lifshitz transition. The ability to tune the topology of band structure by sliding and/or applying normal compressive strain will open doors for controlled tuning of many physical phenomenon such as Landau levels and quantum Hall effect in graphene. (C) 2015 Elsevier Ltd. All rights reserved.

Effect of applied pressure on densification of monolithic ZrCx ceramic by reactive hot pressing

The effect of applied pressure on reactive hot pressing (RHP) of zirconium (Zr):graphite (C) in molar ratios of 1:0.5, 1:0.67, 1:0.8, and 1:1 was studied at 1200 degrees C for 60 min. The relative density achievable increased with increasing pressure and ranged from 99% at 4 MPa for ZrC0.5 to 93% for stoichiometric ZrC at 100 MPa. The diminishing influence of pressure on the final density with increasing stoichiometry is attributed to two causes: the decreasing initial volume fraction of the plastically deforming Zr metal which leads to the earlier formation of a contiguous, stress shielding carbide skeleton and the larger molar volume shrinkage during reaction which leads to pore formation in the final stages. A numerical model of the creep densification of a dynamically evolving microstructure predicts densities that are consistent with observations and confirm that the availability of a soft metal is primarily responsible for the achievement of such elevated densification during RHP. The ability to densify nonstoichiometric compositions like ZrC0.5 at pressures as low as 4 MPa offers an alternate route to fabricating dense nonstoichiometric carbides.