Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
131 tokens/sec
GPT-4o
10 tokens/sec
Gemini 2.5 Pro Pro
47 tokens/sec
o3 Pro
4 tokens/sec
GPT-4.1 Pro
38 tokens/sec
DeepSeek R1 via Azure Pro
28 tokens/sec
2000 character limit reached

Finite volume schemes for multilayer diffusion (1711.10052v2)

Published 28 Nov 2017 in math.NA

Abstract: This paper focusses on finite volume schemes for solving multilayer diffusion problems. We develop a finite volume method that addresses a deficiency of recently proposed finite volume/difference methods, which consider only a limited number of interface conditions and do not carry out stability or convergence analysis. Our method also retains second order accuracy in space while preserving the tridiagonal matrix structure of the classical single-layer discretisation. Stability and convergence analysis of the new finite volume method is presented for each of the three classical time discretisation methods: forward Euler, backward Euler and Crank-Nicolson. We prove that both the backward Euler and Crank-Nicolson schemes are always unconditionally stable. The key contribution of the work is the presentation of a set of sufficient stability conditions for the forward Euler scheme. Here, we find that to ensure stability of the forward Euler scheme it is not sufficient that the time step $\tau$ satisfies the classical constraint of $\tau\leq h_{i}2/(2D_{i})$ in each layer (where $D_{i}$ is the diffusivity and $h_{i}$ is the grid spacing in the $i$th layer) as more restrictive conditions can arise due to the interface conditions. The paper concludes with some numerical examples that demonstrate application of the new finite volume method, with the results presented in excellent agreement with the theoretical analysis.

Summary

We haven't generated a summary for this paper yet.