An official website of the United States government
We consider a model of congestion dynamics with chemotaxis, where the density of cells follows the chemical signal it generates, while observing an incompressibility constraint (incompressible parabolic-elliptic Patlak-Keller-Segel model). We show that when the chemical diffuses slowly and attracts the cells strongly, then the dynamics of the congested cells is well approximated by a surface-tension driven free boundary problem. More precisely, we rigorously establish the convergence of the solution to the characteristic function of a set whose evolution is determined by the classical Hele-Shaw free boundary problem with surface tension. The problem is set in a bounded domain, which leads to an interesting analysis on the limiting boundary conditions. Namely, we prove that the assumption of Robin boundary conditions for the chemical potential leads to a contact angle condition for the free interface (in particular Neumann boundary conditions lead to an orthogonal contact angle condition, while Dirichlet boundary conditions lead to a tangential contact angle condition).
more »
« less
A density-constrained model for chemotaxis
Abstract We consider a model of congestion dynamics with chemotaxis: the density of cells follows a chemical signal it generates, while subject to an incompressibility constraint. The incompressibility constraint results in the formation of patches , describing regions where the maximal density has been reached. The dynamics of these patches can be described by either Hele-Shaw or Richards equation type flow (depending on whether we consider the model with diffusion or the model with pure advection). Our focus in this paper is on the construction of weak solutions for this problem via a variational discrete time scheme of JKO type. We also establish the uniqueness of these solutions. In addition, we make more rigorous the connection between this incompressible chemotaxis model and the free boundary problems describing the motion of the patches in terms of the density and associated pressure variable. In particular, we obtain new results characterising the pressure variable as the solution of an obstacle problem and prove that in the pure advection case the dynamic preserves patches.
more »
« less
- PAR ID:
- 10425073
- Date Published:
- Journal Name:
- Nonlinearity
- Volume:
- 36
- Issue:
- 2
- ISSN:
- 0951-7715
- Page Range / eLocation ID:
- 1082 to 1119
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract This work presents a stabilized formulation for phase‐field fracture of hyperelastic materials near the limit of incompressibility. At this limit, traditional mixed displacement and pressure formulations must satisfy the inf‐sup condition for solution stability. The mixed formulation coupled with the damage field can lead to an inhibition of crack opening as volumetric changes are severely penalized effectively creating a pressure‐bubble. To overcome this bottleneck, we utilize a mixed formulation with a perturbed Lagrangian formulation which enforces the incompressibility constraint in the undamaged material and reduces the pressure effect in the damaged material. A mesh‐dependent stabilization technique based on the residuals of the Euler–Lagrange equations multiplied with a differential operator acting on the weight space is used, allowing for linear interpolation of all field variables of the elastic subproblem. This formulation was validated with three examples at finite deformations: a plane‐stress pure‐shear test, a two‐dimensional geometry in plane‐stress, and a three‐dimensional notched sample. In the last example, we incorporate a hybrid formulation with an additive strain energy decomposition to account for different behaviors in tension and compression. The results show close agreement with analytical solutions for crack tip opening displacements and performs well at the limit of incompressibility.more » « less
-
D. M. Kilgour et al. (Ed.)We consider a Keller-Segel type chemotaxis model with logarithmic sensitivity and density-dependent production/consumption rate. It is a 2 × 2 reaction-diffusion system describing the interaction of cells and a chemical signal. We study Cauchy problem for the original system and its transformed system, which is one of hyperbolic-parabolic conservation laws. In both cases of diffusive and non-diffusive chemical,we obtain optimal L^2 time decay rates for the solution. Our results improve those in Li et al. (Nonlinearity 28:2181-2210, 2015 [5]), Martinez et al. (Indiana Univ Math J 67:1383-1424, 2018 [7]).more » « less
-
Rhenium diboride (ReB2) exhibits high differential strain due to its puckered boron sheets that impede shear deformation. Here, we demonstrate the use of solid solution formation to enhance the Vickers hardness and differential strain of ReB2. ReB2-structured solid solutions (Re0.98Os0.02B2 and Re0.98Ru0.02B2, noted as “ReOsB2” and “ReRuB2”) were synthesized via arc-melting from the pure elements. In-situ high-pressure radial x-ray diffraction was performed in the diamond anvil cell to study the incompressibility and lattice strain of ReOsB2 and ReRuB2 up to ∼56 GPa. Both solid solutions exhibit higher incompressibility and differential strain than pure ReB2. However, while all lattice planes are strengthened by doping osmium (Os) into the ReB2 structure, only the weakest ReB2 lattice plane is enhanced with ruthenium (Ru). These results are in agreement with the Vickers hardness measurements of the two systems, where higher hardness was observed in ReOsB2. The combination of high-pressure studies with experimentally observed hardness data provides lattice specific information about the strengthening mechanisms behind the intrinsic hardness enhancement of the ReB2 system.more » « less
-
We consider a 2D free boundary model of cell motility, inspired by the 1D contraction-driven cell motility model due to P. Recho, T. Putelat, and L. Truskinovsky [Phys. Rev. Lett. 111 (2013), p. 108102]. The key ingredients of the model are the Darcy law for overdamped motion of the acto-myosin network, coupled with the advection-diffusion equation for myosin density. These equations are supplemented with the Young-Laplace equation for the pressure and no-flux condition for the myosin density on the boundary, while evolution of the boundary is subject to the acto-myosin flow at the edge. The focus of the work is on stability analysis of stationary solutions and translationally moving traveling wave solutions. We study stability of radially symmetric stationary solutions and show that at some critical radius a pitchfork bifurcation occurs, resulting in emergence of a family of traveling wave solutions. We perform linear stability analysis of these latter solutions with small velocities and reveal the type of bifurcation (sub- or supercritical). The main result of this work is an explicit asymptotic formula for the stability determining eigenvalue in the limit of small traveling wave velocities.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1361-6544/acad5f
