A convergent finite difference-quadrature scheme for the porous medium equation with nonlocal pressure
Fuente:
arXiv
Guardado en:
| Autores principales: | , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2023
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866912771442999296 |
|---|---|
| author | del Teso, Félix Jakobsen, Espen R. |
| author_facet | del Teso, Félix Jakobsen, Espen R. |
| contents | We introduce and analyze a numerical approximation of the porous medium equation with fractional potential pressure introduced by Caffarelli and Vázquez: \[ \partial_t u = \nabla \cdot (u^{m-1}\nabla (-Δ)^{-σ}u) \qquad \text{for} \qquad m\geq2 \quad \text{and} \quad σ\in(0,1). \] Our scheme is for one space dimension and positive solutions $u$. It consists of solving numerically the equation satisfied by $v(x,t)=\int_{-\infty}^xu(x,t)dx$, the quasilinear non-divergence form equation \[ \partial_t v= -|\partial_x v|^{m-1} (- Δ)^{s} v \qquad \text{where} \qquad s=1-σ, \] and then computing $u=v_x$ by numerical differentiation. Using upwinding ideas in a novel way, we construct a new and simple, monotone and $L^\infty$-stable, approximation for the $v$-equation, and show local uniform convergence to the unique discontinuous viscosity solution. Using ideas from probability theory, we then prove that the approximation of $u$ converges weakly-$*$, or more precisely, up to normalization, in $C(0,T; P(\mathbb{R}))$ where $P(\mathbb{R})$ is the space of probability measures under the Rubinstein-Kantorovich metric.The analysis include also fundamental solutions where the initial data for $u$ is a Dirac mass. Numerical tests are included to confirm the results. Our scheme seems to be the first numerical scheme for this type of problems. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2303_05168 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | A convergent finite difference-quadrature scheme for the porous medium equation with nonlocal pressure del Teso, Félix Jakobsen, Espen R. Numerical Analysis Analysis of PDEs We introduce and analyze a numerical approximation of the porous medium equation with fractional potential pressure introduced by Caffarelli and Vázquez: \[ \partial_t u = \nabla \cdot (u^{m-1}\nabla (-Δ)^{-σ}u) \qquad \text{for} \qquad m\geq2 \quad \text{and} \quad σ\in(0,1). \] Our scheme is for one space dimension and positive solutions $u$. It consists of solving numerically the equation satisfied by $v(x,t)=\int_{-\infty}^xu(x,t)dx$, the quasilinear non-divergence form equation \[ \partial_t v= -|\partial_x v|^{m-1} (- Δ)^{s} v \qquad \text{where} \qquad s=1-σ, \] and then computing $u=v_x$ by numerical differentiation. Using upwinding ideas in a novel way, we construct a new and simple, monotone and $L^\infty$-stable, approximation for the $v$-equation, and show local uniform convergence to the unique discontinuous viscosity solution. Using ideas from probability theory, we then prove that the approximation of $u$ converges weakly-$*$, or more precisely, up to normalization, in $C(0,T; P(\mathbb{R}))$ where $P(\mathbb{R})$ is the space of probability measures under the Rubinstein-Kantorovich metric.The analysis include also fundamental solutions where the initial data for $u$ is a Dirac mass. Numerical tests are included to confirm the results. Our scheme seems to be the first numerical scheme for this type of problems. |
| title | A convergent finite difference-quadrature scheme for the porous medium equation with nonlocal pressure |
| topic | Numerical Analysis Analysis of PDEs |
| url | https://arxiv.org/abs/2303.05168 |