Abstract:We propose numerical schemes for the approximate solution of problems defined on the edges of a one-dimensional graph. In particular, we consider linear transport and a drift-diffusion equations, and discretize them by extending Finite Volume schemes with upwind flux to domains presenting bifurcation nodes with an arbitrary number of incoming and outgoing edges, and implicit time discretization. We show that the discrete problems admit positive unique solutions, and we test the methods on the intricate geometry of an electrical treeing.
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is to develop a numerical method for solving the linear drift - diffusion equations and pure drift equations defined on the edges of one - dimensional graphs. Specifically, the author focuses on how to effectively simulate the evolution process of charged - particle concentrations on complex one - dimensional graph structures with branching nodes. Such problems are particularly important in the simulation of insulation defects in power systems, especially the defect propagation phenomenon in electrical treeing.
### Problem Background
1. **Motivation**:
- The main motivation of this research is to simulate defect propagation in the insulation system of power networks. These defects exist in the form of tree - like structures (called "electrical trees") and continue to develop over time.
- The partial discharges inside these tree - like structures cause the movement of electrons and ions, which in turn triggers an avalanche effect and increases the number of charged particles.
- Currently, although there are some semi - empirical models, models based on first - principles are still lacking, especially when dealing with fully developed complex tree - like geometries.
2. **Objectives**:
- Develop a numerical solver suitable for one - dimensional branching structures to describe the evolution of charged - particle concentrations.
- Ensure that the numerical method can meet the non - negativity condition, which is crucial for subsequent coupling with chemical effects (such as partial discharges and electron avalanche effects).
### Mathematical Models
The paper considers two main types of mathematical models:
1. **Pure Transport Equation**:
\[
\frac{\partial u}{\partial t}+\frac{\partial}{\partial s}(cu) = f,\quad\text{on }\Lambda\times(0, T]
\]
where \(c\) is the transport speed and \(f\) is the source term.
2. **Drift - Diffusion Equation**:
\[
\frac{\partial u}{\partial t}-\frac{\partial}{\partial s}\left(\nu\frac{\partial u}{\partial s}\right)+\frac{\partial}{\partial s}(cu) = f,\quad\text{on }\Lambda\times(0, T]
\]
where \(\nu\) is the diffusion coefficient, \(c\) is the transport speed, and \(f\) is the source term.
### Numerical Method
To solve the above equations, the author proposes a numerical method that combines the finite volume method (FVM) and implicit Euler time discretization:
1. **Spatial Discretization**:
- Use the finite volume method to discretize the space. In particular, use the upwind scheme to handle the flux at the branching nodes.
2. **Time Discretization**:
- Use the implicit Euler method for time discretization to ensure unconditional stability.
3. **Flux Conservation**:
- At each internal node, ensure that the inflow and outflow fluxes are balanced, thereby avoiding mass accumulation.
### Results and Verification
The author verified the effectiveness of this method through a series of tests, including geometries from simple straight lines to complex tree - like structures. The experimental results show that this method can provide a positive and unique solution and exhibits good performance in complex electrical tree - like structures.
### Summary
The main contribution of this paper is to propose a new numerical method that can effectively solve the drift - diffusion equation and the pure transport equation on complex one - dimensional graph structures, and is particularly suitable for simulating the propagation of electrical tree - like defects in power systems. This method not only ensures the non - negativity and uniqueness of the numerical solution but also shows good applicability in complex geometric structures.