Spiral-Wave Dynamics in a Mathematical Model of Human Ventricular Tissue with Myocytes and Fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as E-f, the fibroblast resting-membrane potential, the fibroblast conductance G(f), and the MF gap-junctional coupling G(gap). Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as G(gap), G(f), and E-f, and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity CV decreases as a function of G(gap), for zero-sided and one-sided couplings; however, for two-sided coupling, CV decreases initially and then increases as a function of G(gap), and, eventually, we observe that conduction failure occurs for low values of G(gap). In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling G(gap) or E-f. Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

Experimental and numerical studies on friction welding of thixocast A356 aluminum alloy

This paper highlights the role of globular microstructure on the weldability of semi-solid processed aluminum alloys via high temperature flow behavior. The investigation was carried out on the joining of thixocast A356 aluminum alloy components by friction welding. A thermomechanical model was developed to predict the temperature and stress distributions, as well as to identify the suitable and safe range of parameters. Good comparisons between numerical and experimental results were observed. In addition, metallographic examinations and hardness and tensile tests of the welded samples were carried out. It was found that the tensile strength of the joint is higher than the tensile strength of the parent material for the optimum set of parameters. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Smoothed Functional Algorithms for Stochastic Optimization Using q-Gaussian Distributions

Smoothed functional (SF) schemes for gradient estimation are known to be efficient in stochastic optimization algorithms, especially when the objective is to improve the performance of a stochastic system However, the performance of these methods depends on several parameters, such as the choice of a suitable smoothing kernel. Different kernels have been studied in the literature, which include Gaussian, Cauchy, and uniform distributions, among others. This article studies a new class of kernels based on the q-Gaussian distribution, which has gained popularity in statistical physics over the last decade. Though the importance of this family of distributions is attributed to its ability to generalize the Gaussian distribution, we observe that this class encompasses almost all existing smoothing kernels. This motivates us to study SF schemes for gradient estimation using the q-Gaussian distribution. Using the derived gradient estimates, we propose two-timescale algorithms for optimization of a stochastic objective function in a constrained setting with a projected gradient search approach. We prove the convergence of our algorithms to the set of stationary points of an associated ODE. We also demonstrate their performance numerically through simulations on a queuing model.

Fully developed thermal transport in combined pressure and electroosmotically driven flow of nanofluid in a microchannel under the effect of a magnetic field

This paper critically analyzes, for the first time, the effect of nanofluid on thermally fully developed magnetohydrodynamic flows through microchannel, by considering combined effects of externally applied pressure gradient and electroosmosis. The classical boundary condition of uniform wall heat flux is considered, and the effects of viscous dissipation as well as Joule heating have been taken into account. Closed-form analytical expressions for the pertinent velocity and temperature distributions and the Nusselt number variations are obtained, in order to examine the role of nanofluids in influencing the fully developed thermal transport in electroosmotic microflows under the effect of magnetic field. Fundamental considerations are invoked to ascertain the consequences of particle agglomeration on the thermophysical properties of the nanofluid. The present theoretical formalism addresses the details of the interparticle interaction kinetics in tune with the pertinent variations in the effective particulate dimensions, volume fractions of the nanoparticles, as well as the aggregate structure of the particulate system. It is revealed that the inclusion of nanofluid changes the transport characteristics and system irreversibility to a considerable extent and can have significant consequences in the design of electroosmotically actuated microfluidic systems.

Looking for BSM physics using top-quark polarization and decay-lepton kinematic asymmetries

We explore beyond-standard-model (BSM) physics signatures in the l + jets channel of the t (t) over bar pair production process at the Tevatron and the LHC. We study the effects of BSM physics scenarios on the top-quark polarization and on the kinematics of the decay leptons. To this end, we construct asymmetries using the lepton energy and angular distributions. Further, we find their correlations with the top polarization, net charge asymmetry and top forward-backward asymmetry. We show that when used together, these observables can help discriminate effectively between SM and different BSM scenarios, which can lead to varying degrees of top polarization at the Tevatron as well as the LHC. We use two types of colored mediator models to demonstrate the effectiveness of proposed observables, an s-channel axigluon and a u-channel diquark.

AN INVESTIGATION OF CUTTING MECHANISMS AND STRAIN FIELDS DURING ORTHOGONAL CUTTING IN CFRPs

Fiber-reinforced plastics (FRPs) are typically difficult to machine due to their highly heterogeneous and anisotropic nature and the presence of two phases (fiber and matrix) with vastly different strengths and stiffnesses. Typical machining damage mechanisms in FRPs include series of brittle fractures (especially for thermosets) due to shearing and cracking of matrix material, fiber pull-outs, burring, fuzzing, fiber-matrix debonding, etc. With the aim of understanding the influence of the pronounced heterogeneity and anisotropy observed in FRPs, ``Idealized'' Carbon FRP (I-CFRP) plates were prepared using epoxy resin with embedded equispaced tows of carbon fibers. Orthogonal cutting of these I-CFRPs was carried out, and the chip formation characteristics, cutting force signals and strain distributions obtained during machining were analyzed using the Digital Image Correlation (DIC) technique. In addition, the same procedure was repeated on Uni-Directional CFRPs (UD-CFRPs). Chip formation mechanisms in FRPs were found to depend on the depth of cut and fiber orientation with pure epoxy showing a pronounced ``size effect.'' Experimental results indicate that in-situ full field strain measurements from DIC coupled with force measurements using dynamometry provide an adequate measure of anisotropy and heterogeneity during orthogonal cutting.

Precipitation In Small Systems--II. Mean Field Equations More Effective Than Population Balance

Part I (Manjunath et al., 1994, Chem. Engng Sci. 49, 1451-1463) of this paper showed that the random particle numbers and size distributions in precipitation processes in very small drops obtained by stochastic simulation techniques deviate substantially from the predictions of conventional population balance. The foregoing problem is considered in this paper in terms of a mean field approximation obtained by applying a first-order closure to an unclosed set of mean field equations presented in Part I. The mean field approximation consists of two mutually coupled partial differential equations featuring (i) the probability distribution for residual supersaturation and (ii) the mean number density of particles for each size and supersaturation from which all average properties and fluctuations can be calculated. The mean field equations have been solved by finite difference methods for (i) crystallization and (ii) precipitation of a metal hydroxide both occurring in a single drop of specified initial supersaturation. The results for the average number of particles, average residual supersaturation, the average size distribution, and fluctuations about the average values have been compared with those obtained by stochastic simulation techniques and by population balance. This comparison shows that the mean field predictions are substantially superior to those of population balance as judged by the close proximity of results from the former to those from stochastic simulations. The agreement is excellent for broad initial supersaturations at short times but deteriorates progressively at larger times. For steep initial supersaturation distributions, predictions of the mean field theory are not satisfactory thus calling for higher-order approximations. The merit of the mean field approximation over stochastic simulation lies in its potential to reduce expensive computation times involved in simulation. More effective computational techniques could not only enhance this advantage of the mean field approximation but also make it possible to use higher-order approximations eliminating the constraints under which the stochastic dynamics of the process can be predicted accurately.

Turbulent mixed convection in a shallow enclosure with a series of heat generating components

Turbulent mixed convection flow and heat transfer in a shallow enclosure with and without partitions and with a series of block-like heat generating components is studied numerically for a range of Reynolds and Grashof numbers with a time-dependent formulation. The flow and temperature distributions are taken to be two-dimensional. Regions with the same velocity and temperature distributions can be identified assuming repeated placement of the blocks and fluid entry and exit openings at regular distances, neglecting the end wall effects. One half of such module is chosen as the computational domain taking into account the symmetry about the vertical centreline. The mixed convection inlet velocity is treated as the sum of forced and natural convection components, with the individual components delineated based on pressure drop across the enclosure. The Reynolds number is based on forced convection velocity. Turbulence computations are performed using the standard k– model and the Launder–Sharma low-Reynolds number k– model. The results show that higher Reynolds numbers tend to create a recirculation region of increasing strength in the core region and that the effect of buoyancy becomes insignificant beyond a Reynolds number of typically 5×105. The Euler number in turbulent flows is higher by about 30 per cent than that in the laminar regime. The dimensionless inlet velocity in pure natural convection varies as Gr1/3. Results are also presented for a number of quantities of interest such as the flow and temperature distributions, Nusselt number, pressure drop and the maximum dimensionless temperature in the block, along with correlations.

Long Wavelength Peristaltic Transport of Non-Newton Fluids

Solutions are obtained for the stream function and the pressure field for the flow of non-Newtonian fluids in a tube by long peristaltic waves of arbitrary shape. The axial velocity profiles and stress distributions on the wall are discussed for particular waves of some practical interest. The effect of non- Newtonian character of the fluid is examined.

Two-dimensional solidification in a cylindrical mold with imperfect mold contact

A two-dimensional axisymmetric problem of solidification of a superheated liquid in a long cylindrical mold has been studied in this paper by employing a new embedding technique. The mold and the melt has an imperfect contact and the heat transfer coefficient has been taken as a function of space and time. Short-time exact analytical solutions for the moving boundary and temperature distributions in the liquid, solid and mold have been obtained. The numerical results indicate that with the present solution, for some parameter values, substantial solidified thickness can be obtained. The method of solution is simple and straightforward, and consists of assuming fictitious initial temperatures for some suitable fictitious extensions of the actual regions. Sufficient conditions for the commencement of the solidification have been discussed.

Threshold photoproduction of charged pions on14N

Pion photoproduction processes14Ngs(gamma, pgr +)14C and14Ngs(gamma, pgr –)14O have been studied in the threshold region. These processes provide an excellent tool to study the corrections to soft pion theorems and Kroll-Ruderman limit as applied to nuclear processes. The agreement with the available experimental data for these processes is better with the empirical wave functions while the shell-model wave functions predict a much higher value. Detailed experimental studies of these reactions at threshold, it is shown, are expected to lead to a better understanding of the shell-model inputs and radial distributions in the 1p state. We thank Dr. S.C.K. Nair for a helpful discussion during the initial stages of this work. One of us (MVN) thanks Dr. J.M. Laget for sending some unpublished data on pion photoproduction. He is also thankful to Dr. J. Pasupathy and Dr. R. Rajaraman for their interest and encouragement.

Stationary crack tip fields in elastic-plastic solids: an overview of recent numerical simulations

In this paper, an overview of some recent numerical simulations of stationary crack tip fields in elastic-plastic solids is presented. First, asymptotic analyses carried out within the framework of 2D plane strain or plane stress conditions in both pressure insensitive and pressure sensitive plastic solids are reviewed. This is followed by discussion of salient results obtained from recent computational studies. These pertain to 3D characteristics of elastic-plastic near-front fields under mixed mode loading, mechanics of fracture and simulation of near-tip shear banding process of amorphous alloys and influence of crack tip constraint on the structure of near-tip fields in ductile single crystals. These results serve to illustrate several important features associated with stress and strain distributions near the crack tip and provide the foundation for understanding the operative failure mechanisms. The paper concludes by highlighting some of the future prospects for this field of study.

Effect of surface radiation on conjugate natural convection in a horizontal annulus driven by inner heat generating solid cylinder

This paper reports a numerical study of the laminar conjugate natural convection heat transfer with and without the interaction of the surface radiation in a horizontal cylindrical annulus formed between an inner heat generating solid circular cylinder and an outer isothermal circular boundary. Numerical solutions are obtained by solving the governing equations with a pressure correction method on a collocated (non-staggered) mesh. Steady-state results are presented for the flow and temperature distributions and Nusselt numbers for the heat generation based Grashof number ranging from 10(7) to 10(10), solid-to-fluid thermal conductivity ratios of 1, 5, 10, 50 and 100, radius ratios of 0.226 and 0.452 and surface emissivities of 0-0.8 with air as the working medium. It is observed that surface radiation reduces the convective heat transfer in the annulus compared to the pure natural convection case and enhances the overall Nusselt number.

W physics at the ILC with polarized beams as a probe of the Littlest Higgs model

We study the possibility of using W pair production and leptonic decay of one of the W's at the ILC with polarized beams as a probe of the Littlest Higgs Model. We consider cross-sections, polarization fractions of the W's, leptonic decay energy and angular distributions, and left-right polarization asymmetry as probes of the model. With parameter values allowed by present experimental constraints detectable effects on these observables at typical ILC energies of 500 GeV and 800 GeV will be present. Beam polarization is further found to enhance the sensitivity.

Behaviour of Bi-adhesive Joints

The use of relatively low modulus adhesive at the ends of overlap in a bi-adhesive bondline of a bonded joint can reduce the stress concentration significantly and, therefore, potentially lead to higher strength of the joint. This study presents the two-dimensional and three-dimensional nonlinear (geometric and material) finite element analyses of adhesively bonded single lap joints having modulus-graded bondline under monotonic loading conditions. The adhesives were modelled as an elasto-plastic multi-linear material, while the substrates were regarded as both linear elastic and bi-linear elasto-plastic material. The computational simulations have been performed to investigate the bondline behaviour by studying the stress and strain distributions both at the mid-plane as well as at the interface of the bondline. It has been observed that the static strength is higher for joints with bi-adhesive bondlines compared to those with single adhesives in bondline. Higher joint strength has also been observed for optimum bi-adhesive bondline ratio through parametric studies. Effects of load level, and bondline thickness on stress distribution in the bi-adhesive bondline have also been studied. 3D analysis results reveal the existence of complex multi-axial stress/strain state at the ends of the overlap in the bondline which cannot be observed in 2D plane strain analysis. About 1/3rd of the width of the joint from the free edge in the width direction has 3D stress state, especially in the compliant adhesive of the bondline. Magnitudes of longitudinal and lateral stress/strain components are comparable to peel stress/strain components. It has also been analytically shown that the in-plane global stiffness of the joint remains unaffected by modulus gradation of the bondline adhesive. (C) Koninklijke Brill NV, Leiden, 2010.