Helmholtz decomposition with a scalar Poisson equation in elastic anisotropic media

P-and S-wave separation plays an important role in elastic reverse-time migration.It can reduce the artifacts caused by crosstalk between different modes and improve image quality.In addition,P-and Swave separation can also be used to better understand and distinguish wave types in complex media.At present,the methods for separating wave modes in anisotropic media mainly include spatial nonstationary filtering,low-rank approximation,and vector Poisson equation.Most of these methods require multiple Fourier transforms or the calculation of large matrices,which require high computational costs for problems with large scale.In this paper,an efficient method is proposed to separate the wave mode for anisotropic media by using a scalar anisotropic Poisson operator in the spatial domain.For 2D problems,the computational complexity required by this method is 1/2 of the methods based on solving a vector Poisson equation.Therefore,compared with existing methods based on pseudoHelmholtz decomposition operators,this method can significantly reduce the computational cost.Numerical examples also show that the P and S waves decomposed by this method not only have the correct amplitude and phase relative to the input wavefield but also can reduce the computational complexity significantly.

A Global Spectral Element Model for Poisson Equations and Advective Flow over a Sphere

A global spherical Fourier-Legendre spectral element method is proposed to solve Poisson equations and advective flow over a sphere. In the meridional direction, Legendre polynomials are used and the region is divided into several elements. In order to avoid coordinate singularities at the north and south poles in the meridional direction, Legendre-Gauss-Radau points are chosen at the elements involving the two poles. Fourier polynomials are applied in the zonal direction for its periodicity, with only one element. Then, the partial differential equations are solved on the longitude-latitude meshes without coordinate transformation between spherical and Cartesian coordinates. For verification of the proposed method, a few Poisson equations and advective flows are tested. Firstly, the method is found to be valid for test cases with smooth solution. The results of the Poisson equations demonstrate that the present method exhibits high accuracy and exponential convergence. High- precision solutions are also obtained with near negligible numerical diffusion during the time evolution for advective flow with smooth shape. Secondly, the results of advective flow with non-smooth shape and deformational flow are also shown to be reasonable and effective. As a result, the present method is proved to be capable of solving flow through different types of elements, and thereby a desirable method with reliability and high accuracy for solving partial differential equations over a sphere.

Essential consistency of pressure Poisson equation method and projection method on staggered grids

A new pressure Poisson equation method with viscous terms is established on staggered grids. The derivations show that the newly established pressure equation has the identical equation form in the projection method. The results show that the two methods have the same velocity and pressure values except slight differences in the CPU time.

ON SOLVABILITY OF A BOUNDARY VALUE PROBLEM FOR THE POISSON EQUATION WITH A NONLOCAL BOUNDARY OPERATOR

In this work, we investigate the solvability of the boundary value problem for the Poisson equation, involving a generalized Riemann-Liouville and the Caputo derivative of fractional order in the class of smooth functions. The considered problems are generalization of the known Dirichlet and Neumann oroblems with operators of a fractional order.

AN ACCURATE SOLUTION OF THE POISSON EQUATION BY THE FINITE DIFFERENCE-CHEBYSHEV-TAU METHOD

A new finite difference-Chebyshev-Tau method for the solution of the two-dimensional Poisson equation is presented. Some of the numerical results are also presented which indicate that the method is satisfactory and compatible to other methods.

Numerical solution of Poisson equation with wavelet bases of Hermite cubic splines on the interval

A new wavelet-based finite element method is proposed for solving the Poisson equation. The wavelet bases of Hermite cubic splines on the interval are employed as the multi-scale interpolation basis in the finite element analysis. The lifting scheme of the wavelet-based finite element method is discussed in detail. For the orthogonal characteristics of the wavelet bases with respect to the given inner product, the corresponding multi-scale finite element equation can be decoupled across scales, totally or partially, and suited for nesting approximation. Numerical examples indicate that the proposed method has the higher efficiency and precision in solving the Poisson equation.

A machine learning based solver for pressure Poisson equations

When using the projection method(or fractional step method)to solve the incompressible Navier-Stokes equations,the projection step involves solving a large-scale pressure Poisson equation(PPE),which is computationally expensive and time-consuming.In this study,a machine learning based method is proposed to solve the large-scale PPE.An machine learning(ML)-block is used to completely or partially(if not sufficiently accurate)replace the traditional PPE iterative solver thus accelerating the solution of the incompressible Navier-Stokes equations.The ML-block is designed as a multi-scale graph neural network(GNN)framework,in which the original high-resolution graph corresponds to the discrete grids of the solution domain,graphs with the same resolution are connected by graph convolution operation,and graphs with different resolutions are connected by down/up prolongation operation.The well trained MLblock will act as a general-purpose PPE solver for a certain kind of flow problems.The proposed method is verified via solving two-dimensional Kolmogorov flows(Re=1000 and Re=5000)with different source terms.On the premise of achieving a specified high precision(ML-block partially replaces the traditional iterative solver),the ML-block provides a better initial iteration value for the traditional iterative solver,which greatly reduces the number of iterations of the traditional iterative solver and speeds up the solution of the PPE.Numerical experiments show that the ML-block has great advantages in accelerating the solving of the Navier-Stokes equations while ensuring high accuracy.

Singular Hybrid Boundary Node Method for Solving Poisson Equation

As a boundary-type meshless method, the singular hybrid boundary node method（SHBNM） is based on the modified variational principle and the moving least square（MLS） approximation, so it has the advantages of both boundary element method（BEM） and meshless method. In this paper, the dual reciprocity method（DRM） is combined with SHBNM to solve Poisson equation in which the solution is divided into particular solution and general solution. The general solution is achieved by means of SHBNM, and the particular solution is approximated by using the radial basis function（RBF）. Only randomly distributed nodes on the bounding surface of the domain are required and it doesn＇t need extra equations to compute internal parameters in the domain. The postprocess is very simple. Numerical examples for the solution of Poisson equation show that high convergence rates and high accuracy with a small node number are achievable.

The Numerical Solution of Poisson Equation with Dirichlet Boundary Conditions

This work mainly focuses on the numerical solution of the Poisson equation with the Dirichlet boundary conditions. Compared to the traditional 5-point finite difference method, the Chebyshev spectral method is applied. The numerical results show the Chebyshev spectral method has high accuracy and fast convergence;the more Chebyshev points are selected, the better the accuracy is. Finally, the error of two numerical results also verifies that the algorithm has high precision.

Implementation of the Integrated Green’s Function Method for 3D Poisson’s Equation in a Large Aspect Ratio Computational Domain

The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.

Fractal dynamics and computational analysis of local fractional Poisson equations arising in electrostatics

In this paper,the local fractional natural decomposition method(LFNDM)is used for solving a local fractional Poisson equation.The local fractional Poisson equation plays a significant role in the study of a potential field due to a fixed electric charge or mass density distribution.Numerical examples with computer simulations are presented in this paper.The obtained results show that LFNDM is effective and convenient for application.

A NEW METHOD OF DETERMING SOURCE TERMS OF POISSON EQUATIONS FOR GRID GENERATION

This paper is concerned with a new method to determine the source terms P and Q of Poisson equations for grid generation, with which a satisfactory grid can be obtained. The interior grid distribution is controlled by the prior selection of grid point distribution along the boundary of the region, and the orthogonality condition in the neighborhood of the boundary is satisfied. This grid generation technique can be widely used in the numerical solution of 2 D flow in rivers, lakes, and shallow water regions.

A Lagrange-multiplier-based XFEM to solve pressure Poisson equations in problems with quasi-static interfaces

The XFEM(extended finite element method) has a lot of advantages over other numerical methods to resolve discontinuities across quasi-static interfaces due to the jump in fluidic parameters or surface tension.However,singularities corresponding to enriched degrees of freedom(DOFs) embedded in XFEM arise in the discrete pressure Poisson equations.In this paper,constraints on these DOFs are derived from the interfacial equilibrium condition and introduced in terms of stabilized Lagrange multipliers designed for non-boundary-fitted meshes to address this issue.Numerical results show that the weak and strong discontinuities in pressure with straight and circular interfaces are accurately reproduced by the constraints.Comparisons with the SUPG/PSPG(streamline upwind/pressure stabilizing Petrov-Galerkin) method without Lagrange multipliers validate the applicability and flexibility of the proposed constrained algorithm to model problems with quasi-static interfaces.

Physics-informed neural networks with residual/gradient-based adaptive sampling methods for solving partial differential equations with sharp solutions

We consider solving the forward and inverse partial differential equations(PDEs)which have sharp solutions with physics-informed neural networks(PINNs)in this work.In particular,to better capture the sharpness of the solution,we propose the adaptive sampling methods(ASMs)based on the residual and the gradient of the solution.We first present a residual only-based ASM denoted by ASMⅠ.In this approach,we first train the neural network using a small number of residual points and divide the computational domain into a certain number of sub-domains,then we add new residual points in the sub-domain which has the largest mean absolute value of the residual,and those points which have the largest absolute values of the residual in this sub-domain as new residual points.We further develop a second type of ASM(denoted by ASMⅡ)based on both the residual and the gradient of the solution due to the fact that only the residual may not be able to efficiently capture the sharpness of the solution.The procedure of ASMⅡis almost the same as that of ASMⅠ,and we add new residual points which have not only large residuals but also large gradients.To demonstrate the effectiveness of the present methods,we use both ASMⅠand ASMⅡto solve a number of PDEs,including the Burger equation,the compressible Euler equation,the Poisson equation over an Lshape domain as well as the high-dimensional Poisson equation.It has been shown from the numerical results that the sharp solutions can be well approximated by using either ASMⅠor ASMⅡ,and both methods deliver much more accurate solutions than the original PINNs with the same number of residual points.Moreover,the ASMⅡalgorithm has better performance in terms of accuracy,efficiency,and stability compared with the ASMⅠalgorithm.This means that the gradient of the solution improves the stability and efficiency of the adaptive sampling procedure as well as the accuracy of the solution.Furthermore,we also employ the similar adaptive sampling technique for the data points of boundary conditions(BCs)if the sharpness of the solution is near the boundary.The result of the L-shape Poisson problem indicates that the present method can significantly improve the efficiency,stability,and accuracy.

An Effective Indirect Trefftz Method for Solving Poisson Equation in 2D

In the solution domain, the inhomogeneous part of Poisson equation is approximated with the 5-order polynomial using Galerkin method, and the particular solution of the polynomial can be determined easily. Then, the solution of the Poisson equation is approximated by superposition of the particular solution and the Tcomplete functions related to the Laplace equation. Unknown parameters are determined by Galerkin method, so that the approximate solution is to satisfy the boundary conditions. Comparison with analogous results of others numerical method, the two calculating examples of the paper indicate that the accuracy of the method is very high, which also has a very fast convergence rate.

A New H^2 Regularity Condition of the Solution to Dirichlet Problem of the Poisson Equation and Its Applications

We study the regularity of the solution of Dirichlet problem of Poisson equations over a bounded domain.A new sufficient condition,uniformly positive reach is introduced.Under the assumption that the closure of the underlying domain of interest has a uniformly positive reach,the H^2 regularity of the solution of the Poisson equation is established.In particular,this includes all star-shaped domains whose closures are of positive reach,regardless if they are Lipschitz domains or non-Lipschitz domains.Application to the strong solution to the second order elliptic PDE in non-divergence form and the regularity of Helmholtz equations will be presented to demonstrate the usefulness of the new regularity condition.

Poisson Equation on Some Complete Noncompact Manifolds

We study the Poisson equation on some complete noncompact manifolds with asymptotically nonnegative curvature. We will also study the limiting behavior of the nonhomogeneous heat equation on some complete noncompact manifolds with nonnegative curvature.

Discrete Maximum Principle for Poisson Equation with Mixed Boundary Conditions Solved by hp-FEM

We present a proof of the discrete maximum principle(DMP)for the 1D Poisson equation−u"=f equipped with mixed Dirichlet-Neumann boundary conditions.The problem is discretized using finite elements of arbitrary lengths and polynomial degrees(hp-FEM).We show that the DMP holds on all meshes with no limitations to the sizes and polynomial degrees of the elements.

A Differential Algebraic Method for the Solution of the Poisson Equation for Charged Particle Beams

The design optimization and analysis of charged particle beam systems employing intense beams requires a robust and accurate Poisson solver.This paper presents a new type of Poisson solver which allows the effects of space charge to be elegantly included into the system dynamics.This is done by casting the charge distribution function into a series of basis functions,which are then integrated with an appropriate Green’s function to find a Taylor series of the potential at a given point within the desired distribution region.In order to avoid singularities,a Duffy transformation is applied,which allows singularity-free integration and maximized convergence region when performed with the help of Differential Algebraic methods.The method is shown to perform well on the examples studied.Practical implementation choices and some of their limitations are also explored.

ON THE EXISTENCE OF ALMOST GLOBAL WEAK SOLUTION TO MULTIDIMENSIONAL VLASOV-POISSON EQUATION

The authors prove the existence of almost global weak solution to multidimensional Vlasov Poisson equation with a class of Randon measure as initial data.