763 research outputs found
A Fast Parallel Poisson Solver on Irregular Domains Applied to Beam Dynamic Simulations
We discuss the scalable parallel solution of the Poisson equation within a
Particle-In-Cell (PIC) code for the simulation of electron beams in particle
accelerators of irregular shape. The problem is discretized by Finite
Differences. Depending on the treatment of the Dirichlet boundary the resulting
system of equations is symmetric or `mildly' nonsymmetric positive definite. In
all cases, the system is solved by the preconditioned conjugate gradient
algorithm with smoothed aggregation (SA) based algebraic multigrid (AMG)
preconditioning. We investigate variants of the implementation of SA-AMG that
lead to considerable improvements in the execution times. We demonstrate good
scalability of the solver on distributed memory parallel processor with up to
2048 processors. We also compare our SAAMG-PCG solver with an FFT-based solver
that is more commonly used for applications in beam dynamics
Solving elliptic problems with discontinuities on irregular domains – the Voronoi Interface Method.
We introduce a simple method, dubbed the Voronoi Interface Method, to solve Elliptic problems with discontinuities across the interface of irregular domains. This method produces a linear system that is symmetric positive definite with only its right-hand-side affected by the jump conditions. The solution and the solution's gradients are second-order accurate and first-order accurate, respectively, in the L∞L∞ norm, even in the case of large ratios in the diffusion coefficient. This approach is also applicable to arbitrary meshes. Additional degrees of freedom are placed close to the interface and a Voronoi partition centered at each of these points is used to discretize the equations in a finite volume approach. Both the locations of the additional degrees of freedom and their Voronoi discretizations are straightforward in two and three spatial dimensions
PDE-Based Multidimensional Extrapolation of Scalar Fields over Interfaces with Kinks and High Curvatures
We present a PDE-based approach for the multidimensional extrapolation of
smooth scalar quantities across interfaces with kinks and regions of high
curvature. Unlike the commonly used method of [2] in which normal derivatives
are extrapolated, the proposed approach is based on the extrapolation and
weighting of Cartesian derivatives. As a result, second- and third-order
accurate extensions in the norm are obtained with linear and
quadratic extrapolations, respectively, even in the presence of sharp geometric
features. The accuracy of the method is demonstrated on a number of examples in
two and three spatial dimensions and compared to the approach of [2]. The
importance of accurate extrapolation near sharp geometric features is
highlighted on an example of solving the diffusion equation on evolving
domains.Comment: 17 pages, 13 figures, submitted to SIAM Journal of Scientific
Computin
Immersed Boundary Smooth Extension: A high-order method for solving PDE on arbitrary smooth domains using Fourier spectral methods
The Immersed Boundary method is a simple, efficient, and robust numerical
scheme for solving PDE in general domains, yet it only achieves first-order
spatial accuracy near embedded boundaries. In this paper, we introduce a new
high-order numerical method which we call the Immersed Boundary Smooth
Extension (IBSE) method. The IBSE method achieves high-order accuracy by
smoothly extending the unknown solution of the PDE from a given smooth domain
to a larger computational domain, enabling the use of simple Cartesian-grid
discretizations (e.g. Fourier spectral methods). The method preserves much of
the flexibility and robustness of the original IB method. In particular, it
requires minimal geometric information to describe the boundary and relies only
on convolution with regularized delta-functions to communicate information
between the computational grid and the boundary. We present a fast algorithm
for solving elliptic equations, which forms the basis for simple, high-order
implicit-time methods for parabolic PDE and implicit-explicit methods for
related nonlinear PDE. We apply the IBSE method to solve the Poisson, heat,
Burgers', and Fitzhugh-Nagumo equations, and demonstrate fourth-order pointwise
convergence for Dirichlet problems and third-order pointwise convergence for
Neumann problems
A Computational Study of the Weak Galerkin Method for Second-Order Elliptic Equations
The weak Galerkin finite element method is a novel numerical method that was
first proposed and analyzed by Wang and Ye for general second order elliptic
problems on triangular meshes. The goal of this paper is to conduct a
computational investigation for the weak Galerkin method for various model
problems with more general finite element partitions. The numerical results
confirm the theory established by Wang and Ye. The results also indicate that
the weak Galerkin method is efficient, robust, and reliable in scientific
computing.Comment: 19 page
Finite element methods for surface PDEs
In this article we consider finite element methods for approximating the solution of partial differential equations on surfaces. We focus on surface finite elements on triangulated surfaces, implicit surface methods using level set descriptions of the surface, unfitted finite element methods and diffuse interface methods. In order to formulate the methods we present the necessary geometric analysis and, in the context of evolving surfaces, the necessary transport formulae. A wide variety of equations and applications are covered. Some ideas of the numerical analysis are presented along with illustrative numerical examples
An efficient neural-network and finite-difference hybrid method for elliptic interface problems with applications
A new and efficient neural-network and finite-difference hybrid method is
developed for solving Poisson equation in a regular domain with jump
discontinuities on embedded irregular interfaces. Since the solution has low
regularity across the interface, when applying finite difference discretization
to this problem, an additional treatment accounting for the jump
discontinuities must be employed. Here, we aim to elevate such an extra effort
to ease our implementation by machine learning methodology. The key idea is to
decompose the solution into singular and regular parts. The neural network
learning machinery incorporating the given jump conditions finds the singular
solution, while the standard finite difference method is used to obtain the
regular solution with associated boundary conditions. Regardless of the
interface geometry, these two tasks only require supervised learning for
function approximation and a fast direct solver for Poisson equation, making
the hybrid method easy to implement and efficient. The two- and
three-dimensional numerical results show that the present hybrid method
preserves second-order accuracy for the solution and its derivatives, and it is
comparable with the traditional immersed interface method in the literature. As
an application, we solve the Stokes equations with singular forces to
demonstrate the robustness of the present method
- …