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.

Spiral-wave dynamics in ionically realistic mathematical models for human ventricular tissue: the effects of periodic deformation

We carry out an extensive numerical study of the dynamics of spiral waves of electrical activation, in the presence of periodic deformation (PD) in two-dimensional simulation domains, in the biophysically realistic mathematical models of human ventricular tissue due to (a) ten-Tusscher and Panfilov (the TP06 model) and (b) ten-Tusscher, Noble, Noble, and Panfilov (the TNNPO4 model). We first consider simulations in cable-type domains, in which we calculate the conduction velocity theta and the wavelength lambda of a plane wave; we show that PD leads to a periodic, spatial modulation of theta and a temporally periodic modulation of lambda; both these modulations depend on the amplitude and frequency of the PD. We then examine three types of initial conditions for both TP06 and TNNPO4 models and show that the imposition of PD leads to a rich variety of spatiotemporal patterns in the transmembrane potential including states with a single rotating spiral (RS) wave, a spiral-turbulence (ST) state with a single meandering spiral, an ST state with multiple broken spirals, and a state SA in which all spirals are absorbed at the boundaries of our simulation domain. We find, for both TP06 and TNNPO4 models, that spiral-wave dynamics depends sensitively on the amplitude and frequency of PD and the initial condition. We examine how these different types of spiral-wave states can be eliminated in the presence of PD by the application of low-amplitude pulses by square- and rectangular-mesh suppression techniques. We suggest specific experiments that can test the results of our simulations.

Role of gas in supporting grand spiral structure

The density wave theory for the grand-design two-armed spiral pattern in galaxies is successful in explaining several observed features. However, the long-term persistence of this spiral structure is a serious problem since the group transport would destroy it within about a billion years as shown in a classic paper by Toomre. In this paper, we include the low-velocity dispersion component, namely gas, on an equal footing with stars in the formulation of the density wave theory, and obtain the dispersion relation for this coupled system. We show that the inclusion of gas makes the group transport slower by a factor of few, thus allowing the pattern to persist longer - for several billion years. Though still less than the Hubble time, this helps in making the spiral structure more long-lived. Further we show that addition of gas is essential to get a stable wave for the observed pattern speed for the Galaxy, which otherwise is not possible for a one-component stellar disc.

Effect of dark matter halo on global spiral modes in galaxies

Low surface brightness (LSB) galaxies form a major class of galaxies, and are characterized by low disc surface density and low star formation rate. These are known to be dominated by dark matter halo from the innermost regions. Here, we study the role of the dark matter halo on the grand-design, m = 2, spiral modes in a galactic disc by carrying out a global mode analysis in the WKB approximation. The Bohr-Sommerfeld quantization rule is used to determine how many discrete global spiral modes are permitted. First, a typical superthin, LSB galaxy UGC 7321 is studied by taking only the galactic disc, modelled as a fluid; and then the disc embedded in a dark matter halo. We find that both cases permit the existence of global spiral modes. This is in contrast to earlier results where the inclusion of dark matter halo was shown to nearly fully suppress local, swing-amplified spiral features. Although technically global modes are permitted in the fluid model as shown here, we argue that due to lack of tidal interactions, these are not triggered in LSB galaxies. For comparison, we carried out a similar analysis for the Galaxy, for which the dark matter halo does not dominate in the inner regions. We show that here too the dark matter halo has little effect, hence the disc embedded in a halo is also able to support global modes. The derived pattern speed of the global mode agrees fairly well with the observed value for the Galaxy.

A theoretical approach to the passive control of spiral vortex breakdown

Previous numerical simulations have shown that vortex breakdown starts with the formation of a steady axisymmetric bubble and that an unsteady spiralling mode then develops on top of this.We study how this spiral mode of vortex breakdown might be suppressed or promoted. We use a Lagrangian approach to identify regions of the flow which are sensitive to small open-loop steady and unsteady (harmonic) forces. We find these regions to be upstream of the vortex breakdown bubble. We investigate passive control using a small axisymmetric control ring. In this case, the steady and unsteady control forces are caused by the drag force on the control ring. We find a narrow region upstream of the bubble where the control ring will stabilise the flow and we verify this using numerical simulations. © 2012 IEEE.

Structural sensitivity of spiral vortex breakdown

Previous numerical simulations have shown that vortex breakdown starts with the formation of a steady axisymmetric bubble and that an unsteady spiralling mode then develops on top of this. We investigate this spiral mode with a linear global stability analysis around the steady bubble and its wake. We obtain the linear direct and adjoint global modes of the linearized Navier-Stokes equations and overlap these to obtain the structural sensitivity of the spiral mode, which identifies the wavemaker region. We also identify regions of absolute instability with a local stability analysis. At moderate swirls, we find that the m=-1 azimuthal mode is the most unstable and that the wavemaker regions of the m=-1 mode lie around the bubble, which is absolutely unstable. The mode is most sensitive to feedback involving the radial and azimuthal components of momentum in the region just upstream of the bubble. To a lesser extent, the mode is also sensitive to feedback involving the axial component of momentum in regions of high shear around the bubble. At an intermediate swirl, in which the bubble and wake have similar absolute growth rates, other researchers have found that the wavemaker of the nonlinear global mode lies in the wake. We agree with their analysis but find that the regions around the bubble are more influential than the wake in determining the growth rate and frequency of the linear global mode. The results from this paper provide the first steps towards passive control strategies for spiral vortex breakdown. © 2013 Cambridge University Press.

Layout Design and Optimization of RF Spiral Inductors on Silicon Substrate

The effects of key geometrical parameters on the performance of integrated spiral inductors are investigated with the 3D electromagnetic simulator HFSS. While varying geometrical parameters such as the number of turns （N）,the width of the metal traces （W）,the spacing between the traces （S）,and the inner diameter （ID）, changes in the performance of the inductors are analyzed in detail. The reasons for these changes in performance are presented. Simulation results indicate that the performance of an integrated spiral inductor can be improved by optimizing its layout. Some design rules are summarized.

Development and transition of spiral wave in the coupled Hindmarsh-Rose neurons in two-dimensional space

The dynamics and the transition of spiral waves in the coupled Hindmarsh-Rose (H-R) neurons in two-dimensional space are investigated in the paper. It is found that the spiral wave can be induced and developed in the coupled HR neurons in two-dimensional space, with appropriate initial values and a parameter region given. However, the spiral wave could encounter instability when the intensity of the external current reaches a threshold value of 1.945. The transition of spiral wave is found to be affected by coupling intensity D and bifurcation parameter r. The spiral wave becomes sparse as the coupling intensity increases, while the spiral wave is eliminated and the whole neuronal system becomes homogeneous as the bifurcation parameter increases to a certain threshold value. Then the coupling action of the four sub-adjacent neurons, which is described by coupling coefficient D', is also considered, and it is found that the spiral wave begins to breakup due to the introduced coupling action from the sub-adjacent neurons (or sites) and together with the coupling action of the nearest-neighbour neurons, which is described by the coupling intensity D.

Evolution of spiral wave and pattern formation in a vortical polarized electric field

In this paper, the evolution of the pattern transition induced by the vortical electric field (VEF) is investigated. Firstly, a scheme is suggested to generate the VEF by changing the spatial magnetic field. Secondly, the VEF is imposed on the whole medium, and the evolutions of the spiral wave and the spatiotemporal chaos are investigated by using the numerical simulation. The result confirms that the drift and the breakup of the spiral wave and the new net-like pattern are observed when different polarized fields are imposed on the whole medium respectively. Finally, the pattern transition induced by the polarized field is discussed theoretically.

Shear-induced spiral-like morphology of a main-chain liquid crystalline poly(aryl ether ketone)

The shear-induced spiral-like morphology of a main-chain thermotropic liquid crystalline poly(aryl ether ketone) is observed and characterized by means of polarizing light microscopy, atomic force microscopy, transmission electron microscopy and electron diffraction techniques. The spiral-like texture is formed during shearing in the temperature range of liquid crystalline to isotropic transition (335-340 degreesC), and dispersed discontinuously in the mosaic matrix. Electron diffraction results indicate that the spiral exhibits orthorhombic lateral packing of the crystals and homeotropic alignment of the molecules. The spiral formation process and possible affecting factors are discussed.