Abstract:This paper introduces a second-order method for solving general elliptic partial differential equations (PDEs) on irregular domains using GPU acceleration, based on Ying's kernel-free boundary integral (KFBI) method. The method addresses limitations imposed by CFL conditions in explicit schemes and accuracy issues in fully implicit schemes for the Laplacian operator. To overcome these challenges, the paper employs a series of second-order time discrete schemes and splits the Laplacian operator into explicit and implicit components. Specifically, the Crank-Nicolson method discretizes the heat equation in the temporal dimension, while the implicit scheme is used for the wave equation. The Schrodinger equation is treated using the Strang splitting method. By discretizing the temporal dimension implicitly, the heat, wave, and Schrodinger equations are transformed into a sequence of elliptic equations. The Laplacian operator on the right-hand side of the elliptic equation is obtained from the numerical scheme rather than being discretized and corrected by the five-point difference method. A Cartesian grid-based KFBI method is employed to solve the resulting elliptic equations. GPU acceleration, achieved through a parallel Cartesian grid solver, enhances the computational efficiency by exploiting high degrees of parallelism. Numerical results demonstrate that the proposed method achieves second-order accuracy for the heat, wave, and Schrodinger equations. Furthermore, the GPU-accelerated solvers for the three types of time-dependent equations exhibit a speedup of 30 times compared to CPU-based solvers.

A GPU Acceleration Method for Direct Solutions of Electric Field Integral Equations

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

An Improved Implementation of Multilevel Green's Function Interpolation Method for Electric Field Integral Equation

Discontinuous Galerkin Time Domain Electric Field Integral Equation Method For Solving Non-Conformal Structures

A New Coulomb Gauge Based Electric Field Integral Equation Method

The Implementation of the Three-Dimensional Unified Gas-Kinetic Wave-Particle Method on Multiple Graphics Processing Units

GPU-accelerated Vector-Form Particle-Element Method for 3D Elastoplastic Contact of Structures

OpenMP-CUDA accelerated moment method for homogeneous dielectric objects

FROM CPU TO GPU: GPU-BASED ELECTROMAGNETIC COMPUTING (GPUECO)

GPU Based Two-Level CMFD Accelerating Two-Dimensional MOC Neutron Transport Calculation

GPU_PBTE: an Efficient Solver for Three and Four Phonon Scattering Rates on Graphics Processing Units

Rapid simulation of elastic problems based on GPU

Enhanced Solution to the Surface - Volume - Surface EFIE for Arbitrary Metal - Dielectric Composite Objects

A GPU Poisson-Fermi Solver for Ion Channel Simulations

GPU acceleration of an iterative scheme for gas-kinetic model equations with memory reduction techniques

GPU-accelerated smoothed particle finite element method for large deformation analysis in geomechanics

An adaptive finite element multigrid solver using GPU acceleration

Solving Electromagnetic Scattering Problems with Tens of Billions of Unknowns Using GPU Accelerated Massively Parallel MLFMA

Optimization and acceleration of flow simulations for CFD on CPU/GPU architecture

Efficient GPU-Based Electromagnetic Transient Simulation for Power Systems with Thread-Oriented Transformation and Automatic Code Generation

A GPU-accelerated Cartesian grid method is proposed for solving the heat, wave, and Schrodinger equations on irregular domains