### Abstract

This paper deals with nonlinear diffusion problems including the Stefan problem, the porous medium equation and cross-diffusion systems. A linear discrete-time scheme was proposed by Berger, Brezis and Rogers [RAIRO Anal. Numér.13 (1979) 297–312] for degenerate parabolic equations and was extended to cross-diffusion systems by Murakawa [Math. Mod. Numer. Anal.45 (2011) 1141–1161]. There is a constant stability parameter μ in the linear scheme. In this paper, we propose a linear discrete-time scheme replacing the constant μ with given functions depending on time, space and species. After discretizing the scheme in space, we obtain an easy-to-implement numerical method for the nonlinear diffusion problems. Convergence rates of the proposed discrete-time scheme with respect to the time increment are analyzed theoretically. These rates are the same as in the case where μ is constant. However, actual errors in numerical computation become significantly smaller if varying μ is employed. Our scheme has many advantages even though it is very easy-to-implement, e.g., the ensuing linear algebraic systems are symmetric, it requires low computational cost, the accuracy is comparable to that of the well-studied nonlinear schemes, the computation is much faster than the nonlinear schemes to obtain the same level of accuracy.

Original language | English |
---|---|

Pages (from-to) | 71-101 |

Number of pages | 31 |

Journal | Japan Journal of Industrial and Applied Mathematics |

Volume | 35 |

Issue number | 1 |

DOIs | |

Publication status | Published - Mar 1 2018 |

### Fingerprint

### All Science Journal Classification (ASJC) codes

- Engineering(all)
- Applied Mathematics

### Cite this

**An efficient linear scheme to approximate nonlinear diffusion problems.** / Murakawa, Hideki.

Research output: Contribution to journal › Article

*Japan Journal of Industrial and Applied Mathematics*, vol. 35, no. 1, pp. 71-101. https://doi.org/10.1007/s13160-017-0279-3

}

TY - JOUR

T1 - An efficient linear scheme to approximate nonlinear diffusion problems

AU - Murakawa, Hideki

PY - 2018/3/1

Y1 - 2018/3/1

N2 - This paper deals with nonlinear diffusion problems including the Stefan problem, the porous medium equation and cross-diffusion systems. A linear discrete-time scheme was proposed by Berger, Brezis and Rogers [RAIRO Anal. Numér.13 (1979) 297–312] for degenerate parabolic equations and was extended to cross-diffusion systems by Murakawa [Math. Mod. Numer. Anal.45 (2011) 1141–1161]. There is a constant stability parameter μ in the linear scheme. In this paper, we propose a linear discrete-time scheme replacing the constant μ with given functions depending on time, space and species. After discretizing the scheme in space, we obtain an easy-to-implement numerical method for the nonlinear diffusion problems. Convergence rates of the proposed discrete-time scheme with respect to the time increment are analyzed theoretically. These rates are the same as in the case where μ is constant. However, actual errors in numerical computation become significantly smaller if varying μ is employed. Our scheme has many advantages even though it is very easy-to-implement, e.g., the ensuing linear algebraic systems are symmetric, it requires low computational cost, the accuracy is comparable to that of the well-studied nonlinear schemes, the computation is much faster than the nonlinear schemes to obtain the same level of accuracy.

AB - This paper deals with nonlinear diffusion problems including the Stefan problem, the porous medium equation and cross-diffusion systems. A linear discrete-time scheme was proposed by Berger, Brezis and Rogers [RAIRO Anal. Numér.13 (1979) 297–312] for degenerate parabolic equations and was extended to cross-diffusion systems by Murakawa [Math. Mod. Numer. Anal.45 (2011) 1141–1161]. There is a constant stability parameter μ in the linear scheme. In this paper, we propose a linear discrete-time scheme replacing the constant μ with given functions depending on time, space and species. After discretizing the scheme in space, we obtain an easy-to-implement numerical method for the nonlinear diffusion problems. Convergence rates of the proposed discrete-time scheme with respect to the time increment are analyzed theoretically. These rates are the same as in the case where μ is constant. However, actual errors in numerical computation become significantly smaller if varying μ is employed. Our scheme has many advantages even though it is very easy-to-implement, e.g., the ensuing linear algebraic systems are symmetric, it requires low computational cost, the accuracy is comparable to that of the well-studied nonlinear schemes, the computation is much faster than the nonlinear schemes to obtain the same level of accuracy.

UR - http://www.scopus.com/inward/record.url?scp=85033388081&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85033388081&partnerID=8YFLogxK

U2 - 10.1007/s13160-017-0279-3

DO - 10.1007/s13160-017-0279-3

M3 - Article

AN - SCOPUS:85033388081

VL - 35

SP - 71

EP - 101

JO - Japan Journal of Industrial and Applied Mathematics

JF - Japan Journal of Industrial and Applied Mathematics

SN - 0916-7005

IS - 1

ER -