Abstract:We consider numerical methods for linear parabolic equations in one spatial dimension having piecewise constant diffusion coefficients defined by a one parameter family of interface conditions at the discontinuity. We construct immersed interface finite element methods for an alternative formulation of the original deterministic diffusion problem in which the interface condition is recast as a natural condition on the interfacial flux for which the given operator is self adjoint. An Euler-Maruyama method is developed for the stochastic differential equation corresponding to the alternative divergence formulation of the equation having a discontinuous coefficient and a one-parameter family of interface conditions. We then prove convergence estimates for the Euler scheme. The main goal is to develop numerical schemes that can accommodate specification of any one of the possible interface conditions, and to illustrate the implementation for each of the deterministic and stochastic formulations, respectively. The issues pertaining to speed-ups of the numerical schemes are left to future work.

What problem does this paper attempt to address?

The core problem that this paper attempts to solve is how to numerically solve the linear diffusion equation efficiently and accurately in the presence of interfaces. Specifically, the focus of the research is on the numerical methods for linear parabolic equations with piecewise - constant diffusion coefficients in one - dimensional space. These equations define a family of interface conditions controlled by a single parameter on the interface. ### Specific description of the problem 1. **Treatment of interface conditions**: - The equation has a discontinuous diffusion coefficient \( D(x) \) on the interface and a family of interface conditions controlled by the parameter \( \lambda \) are imposed on the interface. These conditions affect the continuity of the flux and the solution at the interface. - The forms of the interface conditions are: \[ u(t,0^+) = u(t,0^-), \quad \forall t \in J, \] \[ \lambda \frac{\partial u}{\partial x}(t,0^+) = (1-\lambda) \frac{\partial u}{\partial x}(t,0^-), \quad \forall t \in J. \] 2. **Selection of numerical methods**: - The paper proposes the Immersed Interface Finite Element Method (IFEM) to deal with such interface problems. IFEM can construct a discrete grid without depending on the interface position and modify the basis functions to satisfy the interface conditions. - For time discretization, implicit finite - difference methods such as the Backward Euler method and the Crank - Nicolson method are used. 3. **Application of stochastic differential equations**: - The paper also explores the Stochastic Differential Equation (SDE) corresponding to the above - mentioned deterministic diffusion equation and develops the Euler - Maruyama method for numerical solution. - These SDEs include a local time term, which is used to reflect the influence of the interface conditions. ### Main objectives - **Develop numerical schemes**: Construct numerical methods that can adapt to any specified interface conditions and ensure that these methods have high precision and stability when dealing with interface discontinuities. - **Verify theoretical results**: Verify the convergence and error estimates of theoretical derivations through numerical simulations. - **Scope of application**: These methods can be applied in fields such as environmental science, ecology, and biology, especially when dealing with diffusion problems between media with different physical properties. ### Key technical points - **Immersed Interface Finite Element Method (IFEM)**: Used for spatial discretization and can handle interface conditions without depending on the interface position. - **Euler - Maruyama method**: Used for solving stochastic differential equations for time discretization, especially suitable for SDEs containing local time terms. - **Transformation to self - adjoint form**: Transform the original problem into a self - adjoint form, thereby simplifying the numerical solution process. Through these methods, the paper aims to provide a general and efficient numerical framework for dealing with diffusion problems with complex interface conditions.

Numerical Methods for Linear Diffusion Equations in the Presence of an Interface

Lattice Boltzmann Method for Conjugate Heat and Mass Transfer with Interfacial Jump Conditions

Analysis of a numerical technique for a singularly perturbed parabolic convection–diffusion interface problem with time delay

Evolving finite elements for advection diffusion with an evolving interface

A fitted space-time finite element method for an advection-diffusion problem with moving interfaces

Difference potentials method for the nonlinear convection-diffusion equation with interfaces

Mixed finite element methods for nonlinear reaction–diffusion equations with interfaces

Mixed finite element methods for nonlinear reaction-diffusion equations with interfaces

A high-order diffused-interface approach for two-phase compressible flow simulations using a discontinuous Galerkin framework

Advection-Dispersion Across Interfaces

Convergence of numerical schemes for the solution of partial integro-differential equations used in heat transfer

Decoupling Numerical Method Based on Deep Neural Network for Nonlinear Degenerate Interface Problems

Numerical studies of a class of reaction–diffusion equations with Stefan conditions

A Nonconforming Extended Virtual Element Method for Elliptic Interface Problems

A discrete de Rham discretization of interface diffusion problems with application to the Leaky Dielectric Model

Immersed boundary-conformal isogeometric LaTIn method for multiple non-linear interfaces

Numerical methods and regularity properties for viscosity solutions of nonlocal in space and time diffusion equations

High-Order Numerical Method for 1D Non-local Diffusive Equation

Efficient Numerical Algorithms for Three-Dimensional Fractional Partial Differential Equations

Numerical Analysis for a Non-isothermal Incompressible Navier–Stokes–Allen–Cahn System