Search

Information transmission delay can be widely observed in various systems. Here, we study the dynamics of scroll waves with time-delay propagation among slices in excitable media. Weak time delay induces scroll waves to meander. Through increasing the time delay, we find a series of dynamical transitions. Firstly, the straight filament of a scroll wave becomes twisted. Then, the scroll wave breaks and forms interesting patterns. With long time delay, loosed scroll waves are maintained while their period are greatly decreased. Also, cylinder waves appears. The influences of diffusively coupling strength on the time-delay-induced scroll waves are studied. It is found that the critical time delay characterizing those transitions decreases as the coupling strength is increased. A phase diagram in the diffusive coupling-time delay plane is presented.

In this paper, a novel lattice Boltzmann (LB) model based on the Allen-Cahn phase-field theory is proposed for simulating axisymmetric multiphase flows. The most striking feature of the model is that it enables to handle multiphase flows with large density ratio, which are unavailable in all previous axisymmetric LB models. The present model utilizes two LB evolution equations, one of which is used to solve fluid interface, and another is adopted to solve hydrodynamic properties. To simulate axisymmetric multiphase flows effectively, the appropriate source term and equilibrium distribution function are introduced into the LB equation for interface tracking, and simultaneously, a simple and efficient forcing distribution function is also delicately designed in the LB equation for hydrodynamic properties. Unlike many existing LB models, the source and forcing terms of the model arising from the axisymmetric effect include no additional gradients, and consequently, the present model contains only one non-local phase field variable, which in this regard is much simpler. We further conducted the Chapman-Enskog analysis to demonstrate the consistencies of our present MRT-LB model with the axisymmetric Allen-Cahn equation and hydrodynamic equations. A series of numerical examples, including static droplet, oscillation of a viscous droplet, breakup of a liquid thread, and bubble rising in a continuous phase, are used to test the performance of the proposed model. It is found that the present model can generate relatively small spurious velocities and can capture interfacial dynamics with higher accuracy than the previously improved axisymmetric LB model. Besides, it is also found that our present numerical results show excellent agreement with analytical solutions or available experimental data for a wide range of density ratios, which highlights the strengths of the proposed model., Comment: 50 pages, 7 figure

In this paper, a wetting boundary scheme used to describe the interactions among ternary fluids and solid is proposed in the framework of the lattice Boltzmann method. This scheme for three-phase fluids can preserve the reduction consistency property with the diphasic situation such that it could give physically relevant results. Combining this wetting boundary scheme and the lattice Boltzmann (LB) ternary fluid model based on the multicomponent phase-field theory, we simulated several ternary fluid flow problems involving solid substrate, including the spreading of binary drops on the substrate, the spreading of a compound drop on the substrate, and the shear of a compound liquid drop on the substrate. The numerical results are found to be good agreement with the analytical solutions or some available results. Finally, as an application, we use the LB model coupled with the present wetting boundary scheme to numerically investigate the impact of a compound drop on a solid circular cylinder. It is found that the dynamics of a compound drop can be remarkably influenced by the wettability of the solid surface and the dimensionless Weber number., Comment: 32 pages, 14 figures

In this paper, we present a simple and accurate lattice Boltzmann (LB) model for immiscible two-phase flows, which is able to deal with large density contrasts. This model utilizes two LB equations, one of which is used to solve the conservative Allen-Cahn equation, and the other is adopted to solve the incompressible Navier-Stokes equations. A novel forcing distribution function is elaborately designed in the LB equation for the Navier-Stokes equations, which make it much simpler than the existing LB models. In addition, the proposed model can achieve superior numerical accuracy compared with previous Allen-Cahn type of LB models. Several benchmark two-phase problems, including static droplet, layered Poiseuille flow, and Spinodal decomposition are simulated to validate the present LB model. It is found that the present model can achieve relatively small spurious velocities in the LB community, and the obtained numerical results also show good agreement with the analytical solutions or some available results. At last, we use the present model to investigate the droplet impact on a thin liquid film with a large density ratio of 1000 and the Reynolds number ranging from 20 to 500. The fascinating phenomenon of droplet splashing is successfully reproduced by the present model and the numerically predicted spreading radius exhibits to obey the power law reported in the literature., Comment: 31 pages, 8 figures

We study the behavior of spiral wave under the driving of a rotating electric field. The rotating electric field can drive a spiral wave to be synchronous, depending on four factors: its frequency and amplitude, chirality, and polarized modes. Rotation-synchronization characterized by the rotating direction is focused on. We discuss the behavior of synchronization, such as the dependence of angle-differences between the spiral tip and the electric field on ratio of frequency, the influences of different polarized modes of the electric field, the radius of synchronous spiral wave, and so on. A circularly polarized electric fields (CPEF) can suppress meandering spiral to rigid one and prevent breakup of spiral in medium with low excitability. The phase diagram describing the controllable region in excitability-period plane is presented. The influences of polarized modes of electric field on minimum excitability of medium are also studied.

The dynamics of a scroll wave in an excitable medium with gradient excitability is studied in detail. Three parameter regimes can be distinguished by the degree of gradient. For a small gradient, the system reaches a simple rotating synchronization. In this regime, the rigid rotating velocity of spiral waves is maximal in the layers with the highest filament twist. As the excitability gradient increases, the scroll wave evolutes into a meandering synchronous state. This transition is accompanied by a variation in twisting rate. Filament twisting may prevent the breakup of spiral waves in the bottom layers with a low excitability with which a spiral breaks in a 2D medium. When the gradient is large enough, the twisted filament breaks up, which results in a semi-turbulent state where the lower part is turbulent while the upper part contains a scroll wave with a low twisting filament.

We study pattern formation induced by a spiral wave developing from heterogeneities in an excitable medium. Turbulence can be suppressed by a spiral wave from the heterogeneity, forming multiple coexistent systems of regular geometrical patterns. We find that the types of these patterns depend critically on the degree of heterogeneity. The underlying mechanism is due to dispersion relation which is characterized by excitability.

Influences of periodic mechanical deformation (PMD) on spiral breakup that results from Doppler instability in excitable media are investigated. We present a new effect: a high degree of homogeneous PMD is favored to prevent the low-excitability-induced breakup of spiral waves. The frequency and amplitude of PMD are also significant for achieving this purpose. The underlying mechanism of successful control is also discussed, which is believed to be related to the increase of the minimum temporal period of the meandering spiral when the suitable PMD is applied.

Control of turbulence in two kinds of typical heterogeneous excitable media by applying a combined method is investigated. It is found that local-low-amplitude and high-frequency pacing (LHP) is effective to suppress turbulence if the deviation of the heterogeneity is minor. However, LHP is invalid when the deviation is large. Studies show that an additional radial electric field can greatly increase the efficiency of LHP. The underlying mechanisms of successful control in the two kinds of cases are different and are discussed separately. Since the developed strategy of combining LHP with a radial electric field can terminate turbulence in excitable media with a high degree of inhomogeneity, it has the potential contribution to promote the practical low-amplitude defibrillation approach.

The motion of spiral waves in excitable media driven by a weak pacing around the spiral tip is investigated numerically as well as theoretically. We presented a Bifurcations diagram containing four types of the spiral motion induced by different frequencies of pacing: rigidly rotating, inward-petal meandering, resonant drift, and outward-petal meandering spiral. Simulation shows that the spiral resonantly drifts when the frequency of pacing is close to that of the spiral rotation. We also find that the speed and direction of the drift can be efficiently controlled by means of the strength and phase of the local pacing, which is consistent with analytical results based on the framework of the weak deformation approximation.

The control of spiral breakup due to Doppler instability is investigated. It is found that applying an alternating advective field with suitable amplitude and period can prevent the breakup of spiral waves. Further numerical simulations show that the growing meandering behavior of a spiral tip caused by decreasing the excitability of the medium can be efficiently suppressed by the alternating advective field, which inhibits the breakup of spiral waves eventually.

In a generic model of excitable media, we study the behavior of spiral waves interacting with obstacles and their dynamics under the influences of simple periodic mechanical deformation (PMD). Depending on the characteristics of the obstacles, i.e., size and excitability, the rotation of a pinned spiral wave shows different scenarios, e.g., embedding into or anchoring on an obstacle. Three different drift phenomena induced by PMD are observed: scattering on small partial-excitable obstacles, meander-induced unpinning on big partial-excitable obstacles, and drifting around small unexcitable obstacles. Their underlying mechanisms are discussed. The dependence of the threshold amplitude of PMD on the characteristics of the obstacles to successfully remove pinned spiral waves on big partial-excitable obstacles is studied.

Spiral waves may be pinned to anatomical heterogeneities in the cardiac tissue, which leads to monomorphic ventricular tachycardia. Wave emission from heterogeneities (WEH) induced by electric pulses in one direction (EP) is a promising method for liberating such waves by using heterogeneities as internal virtual pacing sites. Here, based on the WEH effect, a new mechanism of liberation by means of a rotating electric pulse (REP) is proposed in a generic model of excitable media. Compared with the EP, the REP has the advantage of opening wider time window to liberate pinned spiral. The influences of rotating direction and frequency of the REP, and the radius of the obstacles on this new mechanism are studied. We believe this strategy may improve manipulations with pinned spiral waves in heart experiments.