A numerical method for solving elliptic equations on real closed algebraic curves and surfaces

Wenrui Hao,Jonathan D. Hauenstein,Margaret H. Regan,Tingting Tang

2023-11-16

Abstract:There are many numerical methods for solving partial different equations (PDEs) on manifolds such as classical implicit, finite difference, finite element, and isogeometric analysis methods which aim at improving the interoperability between finite element method and computer aided design (CAD) software. However, these approaches have difficulty when the domain has singularities since the solution at the singularity may be multivalued. This paper develops a novel numerical approach to solve elliptic PDEs on real, closed, connected, orientable, and almost smooth algebraic curves and surfaces. Our method integrates numerical algebraic geometry, differential geometry, and a finite difference scheme which is demonstrated on several examples.

Numerical Analysis

What problem does this paper attempt to address?

### Problems Addressed by the Paper

The paper aims to solve the problem of solving elliptic partial differential equations (PDEs) on real closed algebraic curves and surfaces. Specifically, the paper focuses on linear elliptic PDEs defined on closed, connected, orientable, and almost smooth algebraic curves and surfaces. These curves and surfaces can be implicitly represented as the solution sets of polynomial equations.

### Main Challenges

1. **Singularity Issue**: When the domain contains singularities, traditional methods such as finite difference methods and finite element methods are difficult to handle because the solution at the singularity may be multivalued.
2. **Parameterization Issue**: For curves and surfaces without global parameterization, traditional methods are hard to apply.
3. **Numerical Stability**: When dealing with embedded manifolds in high-dimensional spaces, traditional methods may be inefficient and numerically unstable.

### Solution

The paper proposes a new numerical method that combines numerical algebraic geometry, differential geometry, and finite difference methods to address the above problems. The specific steps are as follows:

1. **Local Tangential Parameterization**: Use local tangential parameterization at smooth points to simplify the computation of coefficients.
2. **Singularity Handling**: Use numerical algebraic geometry methods to compute the local irreducible branches at singularities and allow the solution at the singularity to take different values on each local irreducible branch.
3. **Numerical Discretization**: Discretize the PDE using finite difference methods to ensure the convergence and stability of the numerical solution.

### Key Contributions

1. **Theoretical Results**: Prove the existence and uniqueness of weak solutions of linear elliptic PDEs on closed parameterized C1 curves under certain conditions.
2. **Numerical Method**: Propose a numerical method based on local tangential parameterization, suitable for curves and surfaces without global parameterization.
3. **Singularity Handling**: Provide a method for handling singularities, ensuring the continuity and multivalued nature of the solution at singularities.

### Application Examples

The paper demonstrates the effectiveness of the method through several specific examples, including solving PDEs on the Lemniscate curve and solving on ellipses under different parameters.

### Conclusion

The method proposed in the paper not only effectively solves the problem of solving elliptic PDEs on closed algebraic curves and surfaces but also provides new ideas for handling singularities, with significant theoretical and practical application value.