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...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.展开更多
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...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.展开更多
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 ...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.展开更多
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 ca...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.展开更多
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 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.展开更多
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...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.展开更多
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,...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.展开更多
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 ...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.展开更多
It will be shown that finding solutions from the Poisson and Klein-Gordon equations under Neumann conditions are equivalent to solving an integral equation, which can be treated as a generalized two-dimensional moment...It will be shown that finding solutions from the Poisson and Klein-Gordon equations under Neumann conditions are equivalent to solving an integral equation, which can be treated as a generalized two-dimensional moment problem over a domain that is considered rectangular. The method consists to solve the integral equation numerically using the two-dimensional inverse moments problem techniques. We illustrate the different cases with examples.展开更多
In this paper, the global existence of weak s olutions to the initial boundary value problem for Boltzmann-Poisson system is proved. The proof is based on the regularization and the stability of the veloci ty averages...In this paper, the global existence of weak s olutions to the initial boundary value problem for Boltzmann-Poisson system is proved. The proof is based on the regularization and the stability of the veloci ty averages and the compactness results on L 1-theory.展开更多
We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditi...We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition (without the incompressibility condition), which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n(n ≥ 3) with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.展开更多
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. ...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.展开更多
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 funct...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.展开更多
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 c...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.展开更多
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 elem...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.展开更多
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 b...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.展开更多
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...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.展开更多
Poisson's equation is solved numerically by two direct methods, viz. Block Cyclic Reduction (BCR) method and Fourier Method. Qualitative and quantitative comparison between the numerical solutions obtained by two ...Poisson's equation is solved numerically by two direct methods, viz. Block Cyclic Reduction (BCR) method and Fourier Method. Qualitative and quantitative comparison between the numerical solutions obtained by two methods indicates that BCR method is superior to Fourier method in terms of speed and accuracy. Therefore. BCR method is applied to solve (?)2(?)= ζ and (?)2X= D from observed vorticity and divergent values. Thereafter the rotational and divergent components of the horizontal monsoon wind in the lower troposphere are reconstructed and are com pared with the results obtained by Successive Over-Relaxation (SOR) method as this indirect method is generally in more use for obtaining the streamfunction ((?)) and velocity potential (X) fields in NWP models. It is found that the results of BCR method are more reliable than SOR method.展开更多
Some new sixth-order compact finite difference schemes for Poisson/Helmholtz equations on rectangular domains in both two-and three-dimensions are developed and analyzed.Different from a few sixth-order compact finite...Some new sixth-order compact finite difference schemes for Poisson/Helmholtz equations on rectangular domains in both two-and three-dimensions are developed and analyzed.Different from a few sixth-order compact finite difference schemes in the literature,the finite difference and weight coefficients of the new methods have analytic simple expressions.One of the new ideas is to use a weighted combination of the source term at staggered grid points which is important for grid points near the boundary and avoids partial derivatives of the source term.Furthermore,the new compact schemes are exact for 2D and 3D Poisson equations if the solution is a polynomial less than or equal to 6.The coefficient matrices of the new schemes are M-matrices for Helmholtz equations with wave number K≤0,which guarantee the discrete maximum principle and lead to the convergence of the new sixth-order compact schemes.Numerical examples in both 2D and 3D are presented to verify the effectiveness of the proposed schemes.展开更多
A computational fluid dynamics(CFD)solver for a GPU/CPU heterogeneous architecture parallel computing platform is developed to simulate incompressible flows on billion-level grid points.To solve the Poisson equation,t...A computational fluid dynamics(CFD)solver for a GPU/CPU heterogeneous architecture parallel computing platform is developed to simulate incompressible flows on billion-level grid points.To solve the Poisson equation,the conjugate gradient method is used as a basic solver,and a Chebyshev method in combination with a Jacobi sub-preconditioner is used as a preconditioner.The developed CFD solver shows good performance on parallel efficiency,which exceeds 90%in the weak-scalability test when the number of grid points allocated to each GPU card is greater than 2083.In the acceleration test,it is found that running a simulation with 10403 grid points on 125 GPU cards accelerates by 203.6x over the same number of CPU cores.The developed solver is then tested in the context of a two-dimensional lid-driven cavity flow and three-dimensional Taylor-Green vortex flow.The results are consistent with previous results in the literature.展开更多
基金supported by the Shandong Post-Doctoral Innovation Fund(Grant No.201303064)the Qingdao Post-Doctoral Application Research Project+1 种基金the National Basic Research(973) Program of China(Grant No.2012CB417402 and 2010CB950402)the National Natural Science Foundation of China(Grant No.41176017)
文摘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.
基金This research was funded by the National Natural Science Foundation of China(Grant No.52108452)the Science Fund for Creative Research Groups of the National Natural Science Foundation of China(Grant No.51921006)the Guangdong Science and Technology Department(Grant No.2020B1212030001).
文摘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.
文摘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.
基金supported by the National Key R&D Program of China(No.2018YFA0702505)the project of CNOOC Limited(Grant No.CNOOC-KJ GJHXJSGG YF 2022-01)+1 种基金R&D Department of China National Petroleum Corporation(Investigations on fundamental experiments and advanced theoretical methods in geophysical prospecting application,2022DQ0604-02)NSFC(Grant Nos.U23B20159,41974142,42074129,12001311)。
文摘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.
文摘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.
文摘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.
文摘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.
基金Project supported by the National Key R&D Program of China(No.2022YFA1004504)the National Natural Science Foundation of China(Nos.12171404 and 12201229)the Fundamental Research Funds for Central Universities of China(No.20720210037)。
文摘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.
文摘It will be shown that finding solutions from the Poisson and Klein-Gordon equations under Neumann conditions are equivalent to solving an integral equation, which can be treated as a generalized two-dimensional moment problem over a domain that is considered rectangular. The method consists to solve the integral equation numerically using the two-dimensional inverse moments problem techniques. We illustrate the different cases with examples.
文摘In this paper, the global existence of weak s olutions to the initial boundary value problem for Boltzmann-Poisson system is proved. The proof is based on the regularization and the stability of the veloci ty averages and the compactness results on L 1-theory.
基金supported in part by the National Science Foundation under Grants DMS-0807551, DMS-0720925, and DMS-0505473the Natural Science Foundationof China (10728101)supported in part by EPSRC grant EP/F029578/1
文摘We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition (without the incompressibility condition), which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n(n ≥ 3) with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.
基金Project supported by the National Natural Science Foundation of China (No. 50876114)
文摘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.
文摘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.
文摘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.
基金supported by the National Natural Science Foundation of China (Nos. 50805028 and 50875195)the Open Foundation of the State Key Laboratory of Structural Analysis for In-dustrial Equipment (No. GZ0815)
文摘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.
基金Foundation item: Supported by the National Natural Science Foundation of China(50608036)
文摘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.
文摘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.
文摘Poisson's equation is solved numerically by two direct methods, viz. Block Cyclic Reduction (BCR) method and Fourier Method. Qualitative and quantitative comparison between the numerical solutions obtained by two methods indicates that BCR method is superior to Fourier method in terms of speed and accuracy. Therefore. BCR method is applied to solve (?)2(?)= ζ and (?)2X= D from observed vorticity and divergent values. Thereafter the rotational and divergent components of the horizontal monsoon wind in the lower troposphere are reconstructed and are com pared with the results obtained by Successive Over-Relaxation (SOR) method as this indirect method is generally in more use for obtaining the streamfunction ((?)) and velocity potential (X) fields in NWP models. It is found that the results of BCR method are more reliable than SOR method.
基金supported by the National Natural Science Foundation of China(Grant No.42274101)and the Excellent Youth Foundation of Hunan Province of China(Grant No.2018JJ1042)Hongling Hu was supported by the National Natural Science Foundation of China(Grant No.12071128)the Natural Science Foundation of Hunan Province(Grant No.2021JJ30434).Zhilin Li was partially supported by a Simons Grant No.633724.
文摘Some new sixth-order compact finite difference schemes for Poisson/Helmholtz equations on rectangular domains in both two-and three-dimensions are developed and analyzed.Different from a few sixth-order compact finite difference schemes in the literature,the finite difference and weight coefficients of the new methods have analytic simple expressions.One of the new ideas is to use a weighted combination of the source term at staggered grid points which is important for grid points near the boundary and avoids partial derivatives of the source term.Furthermore,the new compact schemes are exact for 2D and 3D Poisson equations if the solution is a polynomial less than or equal to 6.The coefficient matrices of the new schemes are M-matrices for Helmholtz equations with wave number K≤0,which guarantee the discrete maximum principle and lead to the convergence of the new sixth-order compact schemes.Numerical examples in both 2D and 3D are presented to verify the effectiveness of the proposed schemes.
基金supported by the National Natural Science Foundation of China (NSFC)Basic Science Center Program for Multiscale Problems in Nonlinear Mechanics’(Grant No. 11988102)NSFC project (Grant No. 11972038)
文摘A computational fluid dynamics(CFD)solver for a GPU/CPU heterogeneous architecture parallel computing platform is developed to simulate incompressible flows on billion-level grid points.To solve the Poisson equation,the conjugate gradient method is used as a basic solver,and a Chebyshev method in combination with a Jacobi sub-preconditioner is used as a preconditioner.The developed CFD solver shows good performance on parallel efficiency,which exceeds 90%in the weak-scalability test when the number of grid points allocated to each GPU card is greater than 2083.In the acceleration test,it is found that running a simulation with 10403 grid points on 125 GPU cards accelerates by 203.6x over the same number of CPU cores.The developed solver is then tested in the context of a two-dimensional lid-driven cavity flow and three-dimensional Taylor-Green vortex flow.The results are consistent with previous results in the literature.