Modeling of terahertz heating effects in realistic tissues

Theterahertz (THz) propagation in real tissues causes heating as with any other electromagnetic radiation propagation. A finite-element (FE) model that provides numerical solutions to the heat conduction equation coupled with realistic models of tissues is employed in this study to compute the temperature raise due to THz propagation. The results indicate that the temperature raise is dependent on the tissue type and is highly localized. The developed FE model was validated through obtaining solutions for the steady-state case and showing that they were in good agreement with the analytical solutions. These types of models can also enable computation of specific absorption rates, which are very critical in planning/setting up experiments involving biological tissues.

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.

Bearing capacity of piles in soft clay underlaid by cohesive frictional soil

By using the axisymmetric quasi-lower bound finite-element limit analysis, the bearing capacity factors N-c(p) and N-gamma q(p) have been computed for axially loaded piles, with the shaft embedded in a fully cohesive soil medium and the tip placed over cohesive frictional soil strata. The results were obtained for various combinations of L/D, phi(l), and c(l)/c(u); the subscripts l and u refer to lower and upper soil strata, respectively. The factors N-c(p) and N-gamma q(p) increase continuously with increases in L/D and phi(l); the rate of increase of N-c(p) and N-gamma q(p) with L/D, however, decreases with an increase in L/D. For c(l)/c(u) > 100, the factor N-c(p) hardly depends on L/D.

Exchange-spring mechanism of soft and hard ferrite nanocomposites

The paper reports exchange-spring soft and hard ferrite nanocomposites synthesized by chemical co-precipitation with or without the application of ultrasonic vibration. The composites contained BaFe12O19 as the hard phase and CoFe2O4/MgFe2O4 as the soft phase. X-ray diffraction patterns of the samples in the optimum calcined condition indicated the presence of soft ferrites as face-centred cubic (fcc) and hard ferrites as hexagonal close packed (hcp) structure respectively. Temperature dependence of magnetization in the range of 20-700 degrees C demonstrated distinct presence of soft and hard ferrites as magnetic phases which are characterized by wide difference in magnetic anisotropy and coercivity. Exchange-spring mechanism led these nanocomposite systems to exchange-coupled, which ultimately produced convex hysteresis loops characteristic of a single-phase permanent magnet. Fairly high value of coercivity and maximum energy product were observed for the samples in the optimum calcined conditions with a maximum applied field of 1600 kA/m (2 T).

A multifold reduction in the transition Reynolds number, and ultra-fast mixing, in a micro-channel due to a dynamical instability induced by a soft wall

A dynamical instability is observed in experimental studies on micro-channels of rectangular cross-section with smallest dimension 100 and 160 mu m in which one of the walls is made of soft gel. There is a spontaneous transition from an ordered, laminar flow to a chaotic and highly mixed flow state when the Reynolds number increases beyond a critical value. The critical Reynolds number, which decreases as the elasticity modulus of the soft wall is reduced, is as low as 200 for the softest wall used here (in contrast to 1200 for a rigid-walled channel) The instability onset is observed by the breakup of a dye-stream introduced in the centre of the micro-channel, as well as the onset of wall oscillations due to laser scattering from fluorescent beads embedded in the wall of the channel. The mixing time across a channel of width 1.5 mm, measured by dye-stream and outlet conductance experiments, is smaller by a factor of 10(5) than that for a laminar flow. The increased mixing rate comes at very little cost, because the pressure drop (energy requirement to drive the flow) increases continuously and modestly at transition. The deformed shape is reconstructed numerically, and computational fluid dynamics (CFD) simulations are carried out to obtain the pressure gradient and the velocity fields for different flow rates. The pressure difference across the channel predicted by simulations is in agreement with the experiments (within experimental errors) for flow rates where the dye stream is laminar, but the experimental pressure difference is higher than the simulation prediction after dye-stream breakup. A linear stability analysis is carried out using the parallel-flow approximation, in which the wall is modelled as a neo-Hookean elastic solid, and the simulation results for the mean velocity and pressure gradient from the CFD simulations are used as inputs. The stability analysis accurately predicts the Reynolds number (based on flow rate) at which an instability is observed in the dye stream, and it also predicts that the instability first takes place at the downstream converging section of the channel, and not at the upstream diverging section. The stability analysis also indicates that the destabilization is due to the modification of the flow and the local pressure gradient due to the wall deformation; if we assume a parabolic velocity profile with the pressure gradient given by the plane Poiseuille law, the flow is always found to be stable.

Hydrodynamics of soft active matter

This review summarizes theoretical progress in the field of active matter, placing it in the context of recent experiments. This approach offers a unified framework for the mechanical and statistical properties of living matter: biofilaments and molecular motors in vitro or in vivo, collections of motile microorganisms, animal flocks, and chemical or mechanical imitations. A major goal of this review is to integrate several approaches proposed in the literature, from semimicroscopic to phenomenological. In particular, first considered are ``dry'' systems, defined as those where momentum is not conserved due to friction with a substrate or an embedding porous medium. The differences and similarities between two types of orientationally ordered states, the nematic and the polar, are clarified. Next, the active hydrodynamics of suspensions or ``wet'' systems is discussed and the relation with and difference from the dry case, as well as various large-scale instabilities of these nonequilibrium states of matter, are highlighted. Further highlighted are various large-scale instabilities of these nonequilibrium states of matter. Various semimicroscopic derivations of the continuum theory are discussed and connected, highlighting the unifying and generic nature of the continuum model. Throughout the review, the experimental relevance of these theories for describing bacterial swarms and suspensions, the cytoskeleton of living cells, and vibrated granular material is discussed. Promising extensions toward greater realism in specific contexts from cell biology to animal behavior are suggested, and remarks are given on some exotic active-matter analogs. Last, the outlook for a quantitative understanding of active matter, through the interplay of detailed theory with controlled experiments on simplified systems, with living or artificial constituents, is summarized.

Numerical analysis of geocell reinforced ballast overlying soft clay subgrade

Geotextiles and geogrids have been in use for several decades in variety of geo-structure applications including foundation of embankments, retaining walls, pavements. Geocells is one such variant in geosynthetic reinforcement of recent years, which provides a three dimensional confinement to the infill material. Although extensive research has been carried on geocell reinforced sand, clay and layered soil subgrades, limited research has been reported on the aggregates/ballast reinforced with geocells. This paper presents the behavior of a railway sleeper subjected to monotonic loading on geocell reinforced aggregates, of size ranging from 20 to 75 mm, overlying soft clay subgrades. Series of tests were conducted in a steel test tank of dimensions 700 mm x 300 mm x 700 mm. In addition to the laboratory model tests, numerical simulations were performed using a finite difference code to predict the behavior of geocell reinforced ballast. The results from numerical simulations were compared with the experimental data. The numerical and experimental results manifested the importance that the geocell reinforcement has a significant effect on the ballast behaviour. The results depicted that the stiffness of underlying soft clay subgrade has a significant influence on the behavior of the geocell-aggregate composite material in redistributing the loading system.

Anomalous phase transitions in soft colloid-polymer binary mixtures

We have shown earlier [1] that these PGNPs resemble star polymers or spherical brushes in terms of their morphology in the melt. However, these particles show dynamics in melt which is quite different from other soft colloidal particles. Since most of the work on soft colloidal particles have been performed in solutions we have now explored the phase behavior of the PGNPs in good solvent using microscopic structural and dynamical measurements on binary mixtures of homopolymers and soft colloids consisting of polymer grafted nanoparticles. We observe anomalous structural and dynamical phase transitions of these binary mixtures, including appearance of spontaneous orientational alignment and logarithmic structural relaxations, as a function of added homopolymers of different molecular weights. Our experiments points to the possibility of exploiting the phase space in density and homopolymer size, of such hybrid systems, to create new materials with unique properties.

Monte Carlo simulation of light transport in turbid medium with embedded object-spherical, cylindrical, ellipsoidal, or cuboidal objects embedded within multilayered tissues

Monte Carlo modeling of light transport in multilayered tissue (MCML) is modified to incorporate objects of various shapes (sphere, ellipsoid, cylinder, or cuboid) with a refractive-index mismatched boundary. These geometries would be useful for modeling lymph nodes, tumors, blood vessels, capillaries, bones, the head, and other body parts. Mesh-based Monte Carlo (MMC) has also been used to compare the results from the MCML with embedded objects (MCML-EO). Our simulation assumes a realistic tissue model and can also handle the transmission/reflection at the object-tissue boundary due to the mismatch of the refractive index. Simulation of MCML-EO takes a few seconds, whereas MMC takes nearly an hour for the same geometry and optical properties. Contour plots of fluence distribution from MCML-EO and MMC correlate well. This study assists one to decide on the tool to use for modeling light propagation in biological tissue with objects of regular shapes embedded in it. For irregular inhomogeneity in the model (tissue), MMC has to be used. If the embedded objects (inhomogeneity) are of regular geometry (shapes), then MCML-EO is a better option, as simulations like Raman scattering, fluorescent imaging, and optical coherence tomography are currently possible only with MCML. (C) 2014 Society of Photo-Optical Instrumentation Engineers (SPIE)

Experimental signatures of a nonequilibrium phase transition governing the yielding of a soft glass

We present direct experimental signatures of a nonequilibrium phase transition associated with the yield point of a prototypical soft solid-a binary colloidal glass. By simultaneously quantifying single-particle dynamics and bulk mechanical response, we identified the threshold for the onset of irreversibility with the yield strain. We extracted the relaxation time from the transient behavior of the loss modulus and found that it diverges in the vicinity of the yield strain. This critical slowing down is accompanied by a growing correlation length associated with the size of regions of high Debye-Waller factor, which are precursors to yield events in glasses. Our results affirm that the paradigm of nonequilibrium critical phenomena is instrumental in achieving a holistic understanding of yielding in soft solids.

Phthalate mediated hydrogelation of a pyrene based system: a novel scaffold for shape-persistent, self-healing luminescent soft material

Two-component super-hydrogelation triggered by the acid-base interaction of a L-histidine appended pyrenyl derivative (PyHis) and phthalic acid (PA) was reported. The use of isomeric isophthalic or terephthalic acid or other comparable acids in place of PA does not lead to salt formation and therefore hydrogelation is not observed. Excimer formation of the pyrenyl unit has not been detected although the PyHis : PA = 1: 1 system undergoes extensive self-assembly in aqueous solution. The synergistic effect of intermolecular H-bonding forces, pi-pi stacking, electrostatic interactions, etc. is found to be responsible for robust hydrogel formation. Development of chiral supramotecular assemblies has been verified through circular dichroism spectroscopy. Morphological investigations involving the PyHis : PA = 1: 1 system show vesicular nano-structures with a definite bilayer width at relatively low concentrations. The latter fuses to construct coiled-coil left-handed helical fibers upon increase in the concentrations of the gelators. The intertwining of the resultant helical fibers eventually results in hydrogel formation. The probable bilayer packing in the self-assembled structures has been probed using X-ray diffraction (XRD) studies and lanthanide sensitization, which suggests that the polar imidazolium hydrogen phthalate unit of the gelator forms the head group and faces the hydrophilic water environment while the hydrophobic pyrenyl units sit inside the hydrophobic core of the bilayer. The hydrogel exhibits multi-stimuli responsiveness including thixotropic behavior. In addition, shape-persistent as well as rapid self-healing behaviour of the hydrogel was established. Furthermore load-bearing characteristics of the hydrogel have also been demonstrated.

OFDM-Aided Differential Space-Time Shift Keying Using Iterative Soft Multiple-Symbol Differential Sphere Decoding

Soft-decision multiple-symbol differential sphere decoding (MSDSD) is proposed for orthogonal frequency-division multiplexing (OFDM)-aided differential space-time shift keying (DSTSK)-aided transmission over frequency-selective channels. Specifically, the DSTSK signaling blocks are generated by the channel-encoded source information and the space-time (ST) blocks are appropriately mapped to a number of OFDM subcarriers. After OFDM demodulation, the DSTSK signal is noncoherently detected by our soft-decision MSDSD detector. A novel soft-decision MSDSD detector is designed, and the associated decision rule is derived for the DSTSK scheme. Our simulation results demonstrate that an SNR reduction of 2 dB is achieved by the proposed scheme using an MSDSD window size of N-w = 4 over the conventional soft-decision-aided differential detection benchmarker, while communicating over dispersive channels and dispensing with channel estimation (CE).

Unraveling the nature of electric field- and stress- induced structural transformations in soft PZT by a new powder poling technique

A `powder-poling' technique was developed to study electric field induced structural transformations in ferroelectrics exhibiting a morphotropic phase boundary (MPB). The technique was employed on soft PZT exhibiting a large longitudinal piezoelectric response (d(33) similar to 650 pCN(-1)). It was found that electric poling brings about a considerable degree of irreversible tetragonal to monoclinic transformation. The same transformation was achieved after subjecting the specimen to mechanical stress, which suggests an equivalence of stress and electric field with regard to the structural mechanism in MPB compositions. The electric field induced structural transformation was also found to be accompanied by a decrease in the spatial coherence of polarization.

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.

Soft-spring wall based non-periodic boundary conditions for non-equilibrium molecular dynamics of dense fluids

Non-equilibrium molecular dynamics (MD) simulations require imposition of non-periodic boundary conditions (NPBCs) that seamlessly account for the effect of the truncated bulk region on the simulated MD region. Standard implementation of specular boundary conditions in such simulations results in spurious density and force fluctuations near the domain boundary and is therefore inappropriate for coupled atomistic-continuum calculations. In this work, we present a novel NPBC model that relies on boundary atoms attached to a simple cubic lattice with soft springs to account for interactions from particles which would have been present in an untruncated full domain treatment. We show that the proposed model suppresses the unphysical fluctuations in the density to less than 1% of the mean while simultaneously eliminating spurious oscillations in both mean and boundary forces. The model allows for an effective coupling of atomistic and continuum solvers as demonstrated through multiscale simulation of boundary driven singular flow in a cavity. The geometric flexibility of the model enables straightforward extension to nonplanar complex domains without any adverse effects on dynamic properties such as the diffusion coefficient. (c) 2015 AIP Publishing LLC.