skip to main content


Title: The leading edge of a free boundary interacting with a line of fast diffusion
The goal of this work is to explain an unexpected feature of the expanding level sets of the solutions of a system where a half-plane in which reaction-diffusion phenomena take place exchanges mass with a line having a large diffusion of its own. The system was proposed by H. Berestycki, L. Rossi and the second author as a model of enhancement of biological invasions by a line of fast diffusion. It was observed numerically by A.-C. Coulon that the leading edge of the front, rather than being located on the line, was in the lower half-plane. We explain this behavior for a closely related free boundary problem. We construct travelling waves for this problem, and the analysis of their free boundary near the line confirms the predictions of the numerical simulations.  more » « less
Award ID(s):
2000041
NSF-PAR ID:
10232046
Author(s) / Creator(s):
Date Published:
Journal Name:
Algebra i analiz
Volume:
32
Issue:
3
ISSN:
0234-0852
Page Range / eLocation ID:
149-179
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In this work, we investigate the two-component modified Korteweg-de Vries (mKdV) equation, which is a complete integrable system, and accepts a generalization of 4 × 4 matrix Ablowitz–Kaup–Newell-Segur (AKNS)-type Lax pair. By using of the unified transform approach, the initial-boundary value (IBV) problem of the two-component mKdV equation associated with a 4 × 4 matrix Lax pair on the half-line will be analyzed. Supposing that the solution {u1(x, t), u2(x, t)} of the two-component mKdV equation exists, we will show that it can be expressed in terms of the unique solution of a 4 × 4 matrix Riemann–Hilbert problem formulated in the complex λ-plane. Moreover, we will prove that some spectral functions s(λ) and S(λ) are not independent of each other but meet the global relationship. 
    more » « less
  2. We address the problem of the potential motion of an ideal incompressible fluid with a free surface and infinite depth in a two-dimensional geometry. We admit the presence of gravity forces and surface tension. A time-dependent conformal mapping $z(w,t)$ of the lower complex half-plane of the variable $w$ into the area filled with fluid is performed with the real line of $w$ mapped into the free fluid’s surface. We study the dynamics of singularities of both $z(w,t)$ and the complex fluid potential $\unicode[STIX]{x1D6F1}(w,t)$ in the upper complex half-plane of $w$ . We show the existence of solutions with an arbitrary finite number $N$ of complex poles in $z_{w}(w,t)$ and $\unicode[STIX]{x1D6F1}_{w}(w,t)$ which are the derivatives of $z(w,t)$ and $\unicode[STIX]{x1D6F1}(w,t)$ over $w$ . We stress that these solutions are not purely rational because they generally have branch points at other positions of the upper complex half-plane. The orders of poles can be arbitrary for zero surface tension while all orders are even for non-zero surface tension. We find that the residues of $z_{w}(w,t)$ at these $N$ points are new, previously unknown, constants of motion, see also Zakharov & Dyachenko (2012, authors’ unpublished observations, arXiv:1206.2046 ) for the preliminary results. All these constants of motion commute with each other in the sense of the underlying Hamiltonian dynamics. In the absence of both gravity and surface tension, the residues of $\unicode[STIX]{x1D6F1}_{w}(w,t)$ are also the constants of motion while non-zero gravity $g$ ensures a trivial linear dependence of these residues on time. A Laurent series expansion of both $z_{w}(w,t)$ and $\unicode[STIX]{x1D6F1}_{w}(w,t)$ at each poles position reveals the existence of additional integrals of motion for poles of the second order. If all poles are simple then the number of independent real integrals of motion is $4N$ for zero gravity and $4N-1$ for non-zero gravity. For the second-order poles we found $6N$ motion integrals for zero gravity and $6N-1$ for non-zero gravity. We suggest that the existence of these non-trivial constants of motion provides an argument in support of the conjecture of complete integrability of free surface hydrodynamics in deep water. Analytical results are solidly supported by high precision numerics. 
    more » « less
  3. Abstract

    We study viscosity solutions to the classical one-phase problem and its thin counterpart. In low dimensions, we show that when the free boundary is the graph of a continuous function, the solution is the half-plane solution. This answers, in the salient dimensions, a one-phase free boundary analogue of Bernstein’s problem for minimal surfaces.As an application, we also classify monotone solutions of semilinear equations with a bump-type nonlinearity.

     
    more » « less
  4. Ru(0001) and Co(0001) films with thickness d ranging from 5 to 300 nm are sputter deposited onto Al2O3(0001) substrates in order to quantify and compare the resistivity size effect. Both metals form epitaxial single crystal layers with their basal planes parallel to the substrate surface and exhibit a root-mean-square roughness < 0.4 nm for Ru and < 0.9 nm for Co. Transport measurements on these layers have negligible resistance contributions from roughness and grain boundary scattering which allows direct quantification of electron surface scattering. The measured resistivity ρ vs d is well described by the classical Fuchs-Sondheimer model, indicating a mean free path for transport within the basal plane of λ = 6.7 ± 0.3 nm for Ru and λ = 19.5 ± 1.0 nm for Co. Bulk Ru is 36% more resistive than Co; in contrast, Ru(0001) layers with d ≤ 25 nm are more conductive than Co(0001) layers, which is attributed to the shorter λ for Ru. The determined λ-values are utilized in combination with the Fuchs-Sondheimer and Mayadas-Shatzkes models to predict and compare the resistance of polycrystalline interconnect lines, assuming a grain boundary reflection coefficient R = 0.4 and accounting for the thinner barrier/adhesion layers available to Ru and Co metallizations. This results in predicted 10 nm half-pitch line resistances for Ru, Co, and Cu of 1.0, 2.2, and 2.1 kΩ/µm, respectively. 
    more » « less
  5. In the phase-field description of moving contact line problems, the two-phase system can be described by free energies, and the constitutive relations can be derived based on the assumption of energy dissipation. In this work we propose a novel boundary condition for contact angle hysteresis by exploring wall energy relaxation, which allows the system to be in non-equilibrium at the contact line. Our method captures pinning, advancing and receding automatically without the explicit knowledge of contact line velocity and contact angle. The microscopic dynamic contact angle is computed as part of the solution instead of being imposed. Furthermore, the formulation satisfies a dissipative energy law, where the dissipation terms all have their physical origin. Based on the energy law, we develop an implicit finite element method that is second order in time. The numerical scheme is proven to be unconditionally energy stable for matched density and zero contact angle hysteresis, and is numerically verified to be energy dissipative for a broader range of parameters. We benchmark our method by computing pinned drops and moving interfaces in the plane Poiseuille flow. When the contact line moves, its dynamics agrees with the Cox theory. In the test case of oscillating drops, the contact line transitions smoothly between pinning, advancing and receding. Our method can be directly applied to three-dimensional problems as demonstrated by the test case of sliding drops on an inclined wall. 
    more » « less