
We present a computational study of the hydrodynamic coarsening in 3D of a
critical mixture using the CahnHilliard/NavierStokes model. The topology of
the resulting intricate bicontinuous microstructure is analyzed through the
principal curvatures to prove selfsimilar morphological evolution. We find
that the self similarity exists for both systems: isoviscous 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 (longtime) grain size distribution (LGSD) that
characterizes the late stage grain coarsening in twodimensional and
quasitwodimensional polycrystalline systems. While experiments and the
phasefield crystal (PFC) model (a simple dynamical density functional theory)
indicate a lognormal distribution, other works including theoretical studies
based on conventional phasefield simulations that rely on coarse grained
fields, like the multiphasefield (MPF) and orientation field (OF) models,
yield significantly different distributions. In a recent work, we have shown
that the coarse grained phasefield 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 advectiondiffusion equations for the environment. The model takes
into account two different bacterial phenotypes. In the motile state, bacteria
swim and perform a runandtumble motion that is biased toward regions of high
oxygen concentration (aerotaxis). In the matrixproducer 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
orientationfield 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
freeenergy functionals of existing models generally contain a term
proportional to the modulus of the orientation gradient, which complicates
their mathematical analysis and induces artificial longrange 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 twodimensional examples are also presented.

The general problem of twophase transport in phasefield 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
sharpinterface 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 phasefield model for the HeleShaw flow of nonNewtonian fluids is
developed. It extends a previous model for Newtonian fluids to a wide range of
sheardependent 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 shearthinning regime at high shear rate. The width selection
of a single steadystate finger is studied in detail for a 2plateaux
shearthinning law (Carreau law) in both its weakly and strongly shearthinning
limits, and the results are related to previous analyses. In the strongly
shearthinning regime a rescaling is found for powerlaw (OstwalddeWaehle)
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
sizedependence dynamics can lead to a novel arrhythmogenic mechanism.
Furthermore, we show that the essential phenomena can be captured in a much
simpler coupledmap 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 LuoRudy 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
earlyafterdepolarizations 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 earlyafterdepolarizations 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 modeI branching). In addition,
we test our model against recent experimental findings showing the presence of
oscillating cracks under biaxial 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 fluctuationinduced nucleation of waves in excitable
media. We utilize a stochastic FitzhughNagumo model for this study, a
spatiallyextended nonpotential 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 nonequilibrium patternforming processes.

Scroll waves are threedimensional analogs of spiral waves. The linear
stability spectrum of untwisted and twisted scroll waves is computed for a
twovariable reactiondiffusion 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 longwavelength
deformations of scroll waves.