American Institute of Mathematical Sciences

August  2021, 41(8): 3985-4012. doi: 10.3934/dcds.2021025

Equilibria of an anisotropic nonlocal interaction equation: Analysis and numerics

 1 Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, United Kingdom 2 Mathematics Institute, University of Warwick, Zeeman Building, Coventry CV4 7AL, United Kingdom 3 Department of Applied Mathematics and Theoretical Physics (DAMTP), University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom

Received  December 2019 Revised  December 2020 Published  January 2021

In this paper, we study the equilibria of an anisotropic, nonlocal aggregation equation with nonlinear diffusion which does not possess a gradient flow structure. Here, the anisotropy is induced by an underlying tensor field. Anisotropic forces cannot be associated with a potential in general and stationary solutions of anisotropic aggregation equations generally cannot be regarded as minimizers of an energy functional. We derive equilibrium conditions for stationary line patterns in the setting of spatially homogeneous tensor fields. The stationary solutions can be regarded as the minimizers of a regularised energy functional depending on a scalar potential. A dimension reduction from the two- to the one-dimensional setting allows us to study the associated one-dimensional problem instead of the two-dimensional setting. We establish $\Gamma$-convergence of the regularised energy functionals as the diffusion coefficient vanishes, and prove the convergence of minimisers of the regularised energy functional to minimisers of the non-regularised energy functional. Further, we investigate properties of stationary solutions on the torus, based on known results in one spatial dimension. Finally, we prove weak convergence of a numerical scheme for the numerical solution of the anisotropic, nonlocal aggregation equation with nonlinear diffusion and any underlying tensor field, and show numerical results.

Citation: José A. Carrillo, Bertram Düring, Lisa Maria Kreusser, Carola-Bibiane Schönlieb. Equilibria of an anisotropic nonlocal interaction equation: Analysis and numerics. Discrete & Continuous Dynamical Systems, 2021, 41 (8) : 3985-4012. doi: 10.3934/dcds.2021025
References:

show all references

References:
Stationary solution to the anisotropic interaction equation (10), obtained with the numerical scheme (33) on a grid of size 50 in each spatial direction and different diffusion coefficients $\delta$ for the spatially homogeneous tensor field with $s = (0,1)$ and $l = (1,0)$ and uniformly distributed initial data on a disc on the computational domain $[-0.5,0.5]^2$
Stationary solution to the anisotropic interaction equation (10), obtained with the numerical scheme (33) on grids of sizes 50,100 and 200 in each spatial direction for the diffusion coefficient $\delta = 10^{-10}$ for the spatially homogeneous tensor field with $s = (0,1)$ and $l = (1,0)$ and uniformly distributed initial data on a disc on the computational domain $[-0.5,0.5]^2$
Stationary solution to the anisotropic interaction equation (10), obtained with the numerical scheme (33) on grids of sizes 100 and 200 in each spatial direction for the diffusion coefficient $\delta = 10^{-10}$ for the spatially homogeneous tensor field with $s = (0,1)$ and $l = (1,0)$ and uniformly distributed initial data on the computational domain $[-0.5,0.5]^2$
Stationary solution to the anisotropic interaction equation (10), obtained with the numerical scheme (33) on a grid of size 200 in each spatial direction and different diffusion coefficients $\delta$ for the spatially homogeneous tensor field with $s = (0,1)$ and $l = (1,0)$ and uniformly distributed initial data on the computational domain $[-0.5,0.5]^2$
Cross-section of stationary solution to the anisotropic interaction equation (10), obtained with the numerical scheme (33) on a grid of size 200 in each spatial direction and diffusion coefficient $\delta = 10^{-9}$ for the spatially homogeneous tensor field with $s = (0,1)$ and $l = (1,0)$ and uniformly distributed initial data on the computational domain $[-0.5,0.5]^2$
Stationary solution to the anisotropic interaction equation (10), obtained with the numerical scheme (33) on a grid of size 50 in each spatial direction and diffusion coefficient $\delta = 10^{-10}$ for different spatially inhomogeneous tensor fields from real fingerprint images and uniformly distributed initial data on the computational domain $[-0.5,0.5]^2$
Numerical solution to the anisotropic interaction equation (10) after $n$ iterations for different $n$, obtained with the numerical scheme (33) on a grid of size 50 in each spatial direction with diffusion coefficient $\delta = 10^{-10}$ for the spatially inhomogeneous tensor field of part of a fingerprint and uniformly distributed initial data on the computational domain $[-0.5,0.5]^2$
Stationary solution to the anisotropic interaction equation (10), obtained with the numerical scheme (33) on a grid of size 50 in each spatial direction for different values of the diffusion coefficient $\delta$ for a given spatially inhomogeneous tensor field and uniformly distributed initial data on the computational domain $[-0.5,0.5]^2$
Stationary solution to the anisotropic interaction equation (10), obtained with the numerical scheme (33) on a grid of size 50 in each spatial direction, diffusion coefficient $\delta = 10^{-10}$ and different force rescalings $\eta$ for a given spatially inhomogeneous tensor field and uniformly distributed initial data on the computational domain $[-0.5,0.5]^2$
 [1] Dong Li, Xiaoyi Zhang. On a nonlocal aggregation model with nonlinear diffusion. Discrete & Continuous Dynamical Systems, 2010, 27 (1) : 301-323. doi: 10.3934/dcds.2010.27.301 [2] R.A. Satnoianu, Philip K. Maini, F.S. Garduno, J.P. Armitage. Travelling waves in a nonlinear degenerate diffusion model for bacterial pattern formation. Discrete & Continuous Dynamical Systems - B, 2001, 1 (3) : 339-362. doi: 10.3934/dcdsb.2001.1.339 [3] Hong Wang, Aijie Cheng, Kaixin Wang. Fast finite volume methods for space-fractional diffusion equations. Discrete & Continuous Dynamical Systems - B, 2015, 20 (5) : 1427-1441. doi: 10.3934/dcdsb.2015.20.1427 [4] Martin Burger, Marco Di Francesco. Large time behavior of nonlocal aggregation models with nonlinear diffusion. Networks & Heterogeneous Media, 2008, 3 (4) : 749-785. doi: 10.3934/nhm.2008.3.749 [5] Yuan Lou, Wei-Ming Ni, Shoji Yotsutani. Pattern formation in a cross-diffusion system. Discrete & Continuous Dynamical Systems, 2015, 35 (4) : 1589-1607. doi: 10.3934/dcds.2015.35.1589 [6] Maxime Breden, Christian Kuehn, Cinzia Soresina. On the influence of cross-diffusion in pattern formation. Journal of Computational Dynamics, 2021, 8 (2) : 213-240. doi: 10.3934/jcd.2021010 [7] Lili Ju, Wensong Wu, Weidong Zhao. Adaptive finite volume methods for steady convection-diffusion equations with mesh optimization. Discrete & Continuous Dynamical Systems - B, 2009, 11 (3) : 669-690. doi: 10.3934/dcdsb.2009.11.669 [8] Mostafa Bendahmane, Mauricio Sepúlveda. Convergence of a finite volume scheme for nonlocal reaction-diffusion systems modelling an epidemic disease. Discrete & Continuous Dynamical Systems - B, 2009, 11 (4) : 823-853. doi: 10.3934/dcdsb.2009.11.823 [9] Shanshan Chen, Junping Shi, Guohong Zhang. Spatial pattern formation in activator-inhibitor models with nonlocal dispersal. Discrete & Continuous Dynamical Systems - B, 2021, 26 (4) : 1843-1866. doi: 10.3934/dcdsb.2020042 [10] Christian Kuehn, Pasha Tkachov. Pattern formation in the doubly-nonlocal Fisher-KPP equation. Discrete & Continuous Dynamical Systems, 2019, 39 (4) : 2077-2100. doi: 10.3934/dcds.2019087 [11] Fengqi Yi, Eamonn A. Gaffney, Sungrim Seirin-Lee. The bifurcation analysis of turing pattern formation induced by delay and diffusion in the Schnakenberg system. Discrete & Continuous Dynamical Systems - B, 2017, 22 (2) : 647-668. doi: 10.3934/dcdsb.2017031 [12] Alexandra Köthe, Anna Marciniak-Czochra, Izumi Takagi. Hysteresis-driven pattern formation in reaction-diffusion-ODE systems. Discrete & Continuous Dynamical Systems, 2020, 40 (6) : 3595-3627. doi: 10.3934/dcds.2020170 [13] Qi Wang, Ling Jin, Zengyan Zhang. Global well-posedness, pattern formation and spiky stationary solutions in a Beddington–DeAngelis competition system. Discrete & Continuous Dynamical Systems, 2020, 40 (4) : 2105-2134. doi: 10.3934/dcds.2020108 [14] Jitraj Saha, Nilima Das, Jitendra Kumar, Andreas Bück. Numerical solutions for multidimensional fragmentation problems using finite volume methods. Kinetic & Related Models, 2019, 12 (1) : 79-103. doi: 10.3934/krm.2019004 [15] Zhangxin Chen. On the control volume finite element methods and their applications to multiphase flow. Networks & Heterogeneous Media, 2006, 1 (4) : 689-706. doi: 10.3934/nhm.2006.1.689 [16] Rajesh Kumar, Jitendra Kumar, Gerald Warnecke. Convergence analysis of a finite volume scheme for solving non-linear aggregation-breakage population balance equations. Kinetic & Related Models, 2014, 7 (4) : 713-737. doi: 10.3934/krm.2014.7.713 [17] Qiang Du, Zhan Huang, Richard B. Lehoucq. Nonlocal convection-diffusion volume-constrained problems and jump processes. Discrete & Continuous Dynamical Systems - B, 2014, 19 (2) : 373-389. doi: 10.3934/dcdsb.2014.19.373 [18] José A. Cañizo, Alexis Molino. Improved energy methods for nonlocal diffusion problems. Discrete & Continuous Dynamical Systems, 2018, 38 (3) : 1405-1425. doi: 10.3934/dcds.2018057 [19] Joseph G. Yan, Dong-Ming Hwang. Pattern formation in reaction-diffusion systems with $D_2$-symmetric kinetics. Discrete & Continuous Dynamical Systems, 1996, 2 (2) : 255-270. doi: 10.3934/dcds.1996.2.255 [20] Phillip Colella. High-order finite-volume methods on locally-structured grids. Discrete & Continuous Dynamical Systems, 2016, 36 (8) : 4247-4270. doi: 10.3934/dcds.2016.36.4247

2019 Impact Factor: 1.338