• We present a computational study of the hydrodynamic coarsening in 3D of a critical mixture using the Cahn-Hilliard/Navier-Stokes model. The topology of the resulting intricate bicontinuous microstructure is analyzed through the principal curvatures to prove self-similar morphological evolution. We find that the self similarity exists for both systems: iso-viscous and with variable viscosity. However the two system have distinct topological character. Our simulations confirm that the predicted viscous growth regime exists in both cases. Moreover the coarsening rate is inversely proportional to an \textit{effective viscosity} that is the geometrical average of the viscosities of the two phases.
  • In the literature, contradictory results have been published regarding the form of the limiting (long-time) grain size distribution (LGSD) that characterizes the late stage grain coarsening in two-dimensional and quasi-two-dimensional polycrystalline systems. While experiments and the phase-field crystal (PFC) model (a simple dynamical density functional theory) indicate a lognormal distribution, other works including theoretical studies based on conventional phase-field simulations that rely on coarse grained fields, like the multi-phase-field (MPF) and orientation field (OF) models, yield significantly different distributions. In a recent work, we have shown that the coarse grained phase-field models (whether MPF or OF) yield very similar limiting size distributions that seem to differ from the theoretical predictions. Herein, we revisit this problem, and demonstrate in the case of OF models [by R. Kobayashi et al., Physica D 140, 141 (2000) and H. Henry et al. Phys. Rev. B 86, 054117 (2012)] that an insufficient resolution of the small angle grain boundaries leads to a lognormal distribution close to those seen in the experiments and the molecular scale PFC simulations. Our work indicates, furthermore, that the LGSD is critically sensitive to the details of the evaluation process, and raises the possibility that the differences among the LGSD results from different sources may originate from differences in the detection of small angle grain boundaries.
  • The evolution of a crack front under mixed mode loading (I+III) is studied using a phase field model in 3 dimensions with no stress boundary conditions. As previously observed experimentally in gels, there is a relaxation toward a geometry where $K_{III}=0$ without any front fragmentation even for high values of the initial mode mixity $K_{III}/K_{I}$. The effects of the initial condition is studied and it is shown that irregularities in the initial slit can lead to front fragmentation for smaller values of the ratio $K_{III}/K_{I}$ as is observed in experiments.
  • The bacterium {\em Bacilus subtilis} frequently forms biofilms at the interface between the culture medium and the air. We develop a mathematical model that couples a description of bacteria as individual discrete objects to the standard advection-diffusion equations for the environment. The model takes into account two different bacterial phenotypes. In the motile state, bacteria swim and perform a run-and-tumble motion that is biased toward regions of high oxygen concentration (aerotaxis). In the matrix-producer state they excrete extracellular polymers, which allows them to connect to other bacteria and to form a biofilm. Bacteria are also advected by the fluid, and can trigger bioconvection. Numerical simulations of the model reproduce all the stages of biofilm formation observed in laboratory experiments. Finally, we study the influence of various model parameters on the dynamics and morphology of biofilms.
  • The solidification of polycrystalline materials can be modelled by orientation-field models, which are formulated in terms of two continuous fields: a phase field that describes the thermodynamic state and an orientation field that indicates the local direction of the crystallographic axes. The free-energy functionals of existing models generally contain a term proportional to the modulus of the orientation gradient, which complicates their mathematical analysis and induces artificial long-range interactions between grain boundaries. We present an alternative model, in which only the square of the orientation gradient appears, but in which the phase and orientation fields are coupled by a singular function that diverges in the solid phase. We show that this model exhibits stable grain boundaries whose interactions decay exponentially with their distance. Furthermore, we demonstrate that the anisotropy of the surface energy can be included while preserving the variational structure of the model. Illustrative numerical simulations of two-dimensional examples are also presented.
  • The general problem of two-phase transport in phase-field models is analyzed: the flux of a conserved quantity is driven by the gradient of a potential through a medium that consists of domains of two distinct phases which are separated by diffuse interfaces. It is shown that the finite thickness of the interfaces induces two effects that are not present in the analogous sharp-interface problem: a surface excess current and a potential jump at the interfaces. It is shown that both effects can be eliminated simultaneously only if the coefficient of proportionality between flux and potential gradient (mobility) is allowed to become a tensor in the interfaces. This opens the possibility for precise and efficient simulations of transport problems with finite interface thickness.
  • The shape of a crack front propagating through a thin sample is studied using a phase field model. The model is shown to have a well defined sharp interface limit. The crack front is found to be an ellipse with large axis the width of the sample and small axis a function of the Poisson ratio and the width of the sample. Numerical results also indicate that the front shape is independent of the crack speed and of the sample extension perpendicular to its width.
  • A phase-field model for the Hele-Shaw flow of non-Newtonian fluids is developed. It extends a previous model for Newtonian fluids to a wide range of shear-dependent fluids. The model is applied to perform simulations of viscous fingering in shear- thinning fluids, and it is found to be capable of describing the complete crossover from the Newtonian regime at low shear rate to the strongly shear-thinning regime at high shear rate. The width selection of a single steady-state finger is studied in detail for a 2-plateaux shear-thinning law (Carreau law) in both its weakly and strongly shear-thinning limits, and the results are related to previous analyses. In the strongly shear-thinning regime a rescaling is found for power-law (Ostwald-de-Waehle) fluids that allows for a direct comparison between simulations and experiments without any adjustable parameters, and good agreement is obtained.
  • We study numerically the dynamics of conduction blocks using a detailed electrophysiological model. We find that this dynamics depends critically on the size of the paced region. Small pacing regions lead to stationary conduction blocks while larger pacing regions can lead to conduction blocks that travel periodically towards the pacing region. We show that this size-dependence dynamics can lead to a novel arrhythmogenic mechanism. Furthermore, we show that the essential phenomena can be captured in a much simpler coupled-map model.
  • In this numerical study, we investigate the role of intrinsic heterogeneities of cardiac tissue due to M cells in the generation and maintenance of reentrant excitations using the detailed Luo-Rudy dynamic model. This model has been extended to include a description of the long QT 3 syndrome, and is studied in both one dimension, corresponding to a cable traversing the ventricular wall, and two dimensions, representing a transmural slice. We focus on two possible mechanisms for the generation of reentrant events. We first investigate if early-after-depolarizations occurring in M cells can initiate reentry. We find that, even for large values of the long QT strength, the electrotonic coupling between neighboring cells prevents early-after-depolarizations from creating a reentry. We then study whether M cell domains, with their slow repolarization, can function as wave blocks for premature stimuli. We find that the inclusion of an M cell domain can result in some cases in reentrant excitations and we determine the lifetime of the reentry as a function of the size and geometry of the domain and of the strength of the long QT syndrome.
  • We present a phase field model of the propagation of fracture under plane strain. This model, based on simple physical considerations, is able to accurately reproduce the different behavior of cracks (the principle of local symmetry, the Griffith and Irwin criteria, and mode-I branching). In addition, we test our model against recent experimental findings showing the presence of oscillating cracks under bi-axial load. Our model again reproduces well observed supercritical Hopf bifurcation, and is therefore the first simulation which does so.
  • I present the numerical computation of speed and direction of the drift of a spiral wave in an excitable medium in the presence of an electric field. In contrast to earlier results, the drift speed presents a strong variation close to the parameter value where the drift speed component along the field changes direction. Using a simple phenomenological model and results from a numerical linear stability analysis of scroll waves, I show this behavior can be attributed to a resonance of the meander modes with the translation modes of the spiral wave. Extending this phenomenological model to scroll waves also clarifies the link between the drift and long wavelength instabilities of scroll waves.
  • Motivated by recent experiments on intracellular calcium dynamics, we study the general issue of fluctuation-induced nucleation of waves in excitable media. We utilize a stochastic Fitzhugh-Nagumo model for this study, a spatially-extended non-potential pair of equations driven by thermal (i.e. white) noise. The nucleation rate is determined by finding the most probable escape path via minimization of an action related to the deviation of the fields from their deterministic trajectories. Our results pave the way both for studies of more realistic models of calcium dynamics as well as of nucleation phenomena in other non-equilibrium pattern-forming processes.
  • Scroll waves are three-dimensional analogs of spiral waves. The linear stability spectrum of untwisted and twisted scroll waves is computed for a two-variable reaction-diffusion model of an excitable medium. Different bands of modes are seen to be unstable in different regions of parameter space. The corresponding bifurcations and bifurcated states are characterized by performing direct numerical simulations. In addition, computations of the adjoint linear stability operator eigenmodes are also performed and serve to obtain a number of matrix elements characterizing the long-wavelength deformations of scroll waves.