Gradient discretization of Hybrid Dimensional Darcy Flows in Fractured Porous Media with discontinuous pressures at the matrix fracture interfaces

We investigate the discretization of Darcy flow through fractured porous media on general meshes. We consider a hybrid dimensional model, invoking a complex network of planar fractures. The model accounts for matrix-fracture interactions and fractures acting either as drains or as barriers, i.e. we have to deal with pressure discontinuities at matrix-fracture interfaces. The numerical analysis is performed in the general framework of gradient discretizations which is extended to the model under consideration. Two families of schemes namely the Vertex Approximate Gradient scheme (VAG) and the Hybrid Finite Volume scheme (HFV) are detailed and shown to satisfy the gradient scheme framework, which yields, in particular, convergence. Numerical tests confirm the theoretical results. Gradient Discretization; Darcy Flow, Discrete Fracture Networks, Finite Volum

On the stabilization of ion sputtered surfaces

The classical theory of ion beam sputtering predicts the instability of a flat surface to uniform ion irradiation at any incidence angle. We relax the assumption of the classical theory that the average surface erosion rate is determined by a Gaussian response function representing the effect of the collision cascade and consider the surface dynamics for other physically-motivated response functions. We show that although instability of flat surfaces at any beam angle results from all Gaussian and a wide class of non-Gaussian erosive response functions, there exist classes of modifications to the response that can have a dramatic effect. In contrast to the classical theory, these types of response render the flat surface linearly stable, while imperceptibly modifying the predicted sputter yield vs. incidence angle. We discuss the possibility that such corrections underlie recent reports of a ``window of stability'' of ion-bombarded surfaces at a range of beam angles for certain ion and surface types, and describe some characteristic aspects of pattern evolution near the transition from unstable to stable dynamics. We point out that careful analysis of the transition regime may provide valuable tests for the consistency of any theory of pattern formation on ion sputtered surfaces

The fractional Keller-Segel model

The Keller-Segel model is a system of partial differential equations modelling chemotactic aggregation in cellular systems. This model has blowing up solutions for large enough initial conditions in dimensions d >= 2, but all the solutions are regular in one dimension; a mathematical fact that crucially affects the patterns that can form in the biological system. One of the strongest assumptions of the Keller-Segel model is the diffusive character of the cellular motion, known to be false in many situations. We extend this model to such situations in which the cellular dispersal is better modelled by a fractional operator. We analyze this fractional Keller-Segel model and find that all solutions are again globally bounded in time in one dimension. This fact shows the robustness of the main biological conclusions obtained from the Keller-Segel model

Development of Knife-Edge Ridges on Ion-Bombarded Surfaces

We demonstrate in both laboratory and numerical experiments that ion bombardment of a modestly sloped surface can create knife-edge like ridges with extremely high slopes. Small pre-fabricated pits expand under ion bombardment, and the collision of two such pits creates knife-edge ridges. Both laboratory and numerical experiments show that the pit propagation speed and the precise shape of the knife edge ridges are universal, independent of initial conditions, as has been predicted theoretically. These observations suggest a novel method of fabrication in which a surface is pre-patterned so that it dynamically evolves to a desired target pattern made of knife-edge ridges.Comment: 5 pages, 4 figure

The Two Fluid Drop Snap-off Problem: Experiments and Theory

We address the dynamics of a drop with viscosity $\lambda \eta$ breaking up inside another fluid of viscosity $\eta$. For $\lambda=1$, a scaling theory predicts the time evolution of the drop shape near the point of snap-off which is in excellent agreement with experiment and previous simulations of Lister and Stone. We also investigate the $\lambda$ dependence of the shape and breaking rate.Comment: 4 pages, 3 figure