A novel numerical method is explored and named as mesh-free poly-cell Galerkin method. An improved moving least-square (MLS) scheme is presented, which can avoid the matrix inversion in standard MLS and can be used ...A novel numerical method is explored and named as mesh-free poly-cell Galerkin method. An improved moving least-square (MLS) scheme is presented, which can avoid the matrix inversion in standard MLS and can be used to construct shape functions possessing delta Kronecher property. A new type of local support is introduced to ensure the alignment of integral domains with the cells of the back-ground mesh, which will reduce the difficult in integration. An intensive numerical study is conducted to test the accuracy of the present method. It is observed that solutions with good accuracy can be obtained with the present method.展开更多
In this paper,we analyse the equal width(EW) wave equation by using the mesh-free reproducing kernel particle Ritz(kp-Ritz) method.The mesh-free kernel particle estimate is employed to approximate the displacement...In this paper,we analyse the equal width(EW) wave equation by using the mesh-free reproducing kernel particle Ritz(kp-Ritz) method.The mesh-free kernel particle estimate is employed to approximate the displacement field.A system of discrete equations is obtained through the application of the Ritz minimization procedure to the energy expressions.The effectiveness of the kp-Ritz method for the EW wave equation is investigated by numerical examples in this paper.展开更多
In this study,to simulate open channel flows,an explicit incompressible mesh-free method is employed in which the pressure field is obtained by explicitly solving the pressure Poisson equation.To capture the velocity ...In this study,to simulate open channel flows,an explicit incompressible mesh-free method is employed in which the pressure field is obtained by explicitly solving the pressure Poisson equation.To capture the velocity information in open channel flows,the source term in the pressure Poisson equation is modified while the spatial discretization of gradient and Laplacian models is based on the moving particle semi-implicit(MPS)method.The inflow boundary condition is treated by injecting fluid particles into the domain according to the inlet discharge,and the outflow condition is handled by prescribing the pressure distribution and removing the fluid particles beyond the domain.The explicit incompressible mesh-free method is then used to simulate open channel flows,including weir flow,hydraulic jump,and flow over an obstacle.In the simulations,velocity distribution and flow pattern are examined.The simulated results are compared to available experimental measurements and other numerical results.There is a good agreement between the simulated results and the experimental measurements.It shows that the explicit incompressible mesh-free method can reproduce the flow characteristics in the open channel flows.展开更多
A mesh-free method based on local Petrov-Galerkin formulation is presented to solve dynamic impact problems of hyperelastic material.In the present method,a simple Heaviside test function is chosen for simplifying dom...A mesh-free method based on local Petrov-Galerkin formulation is presented to solve dynamic impact problems of hyperelastic material.In the present method,a simple Heaviside test function is chosen for simplifying domain integrals.Trial function is constructed by using a radial basis function (RBF) coupled with a polynomial basis function,in which the shape function possesses the kronecker delta function property.So,additional treatment is not required for imposing essential boundary conditions.Governing equations of impact problems are established and solved node by node by using an explicit time integration algorithm in a local domain,which is very similar to that of the collocation method except that numerical integration can be implemented over local domain in the present method.Numerical results for several examples show that the present method performs well in dealing with the dynamic impact problem of hyperelastic material.展开更多
This paper is dedicated to the expansion of the framework of general interpolant observables introduced by Azouani,Olson,and Titi for continuous data assimilation of nonlinear partial differential equations.The main f...This paper is dedicated to the expansion of the framework of general interpolant observables introduced by Azouani,Olson,and Titi for continuous data assimilation of nonlinear partial differential equations.The main feature of this expanded framework is its mesh-free aspect,which allows the observational data itself to dictate the subdivision of the domain via partition of unity in the spirit of the so-called Partition of Unity Method by Babuska and Melenk.As an application of this framework,we consider a nudging-based scheme for data assimilation applied to the context of the two-dimensional Navier-Stokes equations as a paradigmatic example and establish convergence to the reference solution in all higher-order Sobolev topologies in a periodic,mean-free setting.The convergence analysis also makes use of absorbing ball bounds in higherorder Sobolev norms,for which explicit bounds appear to be available in the literature only up to H^(2);such bounds are additionally proved for all integer levels of Sobolev regularity above H^(2).展开更多
In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be r...In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.展开更多
文摘A novel numerical method is explored and named as mesh-free poly-cell Galerkin method. An improved moving least-square (MLS) scheme is presented, which can avoid the matrix inversion in standard MLS and can be used to construct shape functions possessing delta Kronecher property. A new type of local support is introduced to ensure the alignment of integral domains with the cells of the back-ground mesh, which will reduce the difficult in integration. An intensive numerical study is conducted to test the accuracy of the present method. It is observed that solutions with good accuracy can be obtained with the present method.
基金Project supported by the Natural Science Foundation of Zhejiang Province,China (Grant No. Y6110007)
文摘In this paper,we analyse the equal width(EW) wave equation by using the mesh-free reproducing kernel particle Ritz(kp-Ritz) method.The mesh-free kernel particle estimate is employed to approximate the displacement field.A system of discrete equations is obtained through the application of the Ritz minimization procedure to the energy expressions.The effectiveness of the kp-Ritz method for the EW wave equation is investigated by numerical examples in this paper.
基金This work was supported by the Key Research and Development Program of Zhejiang Province(Grant No.2020C03082)the Visiting Researcher Fund Program of State Key Laboratory of Water Resources and Hydropower Engineering Science,Wuhan University(Grant No.2021HLG01).
文摘In this study,to simulate open channel flows,an explicit incompressible mesh-free method is employed in which the pressure field is obtained by explicitly solving the pressure Poisson equation.To capture the velocity information in open channel flows,the source term in the pressure Poisson equation is modified while the spatial discretization of gradient and Laplacian models is based on the moving particle semi-implicit(MPS)method.The inflow boundary condition is treated by injecting fluid particles into the domain according to the inlet discharge,and the outflow condition is handled by prescribing the pressure distribution and removing the fluid particles beyond the domain.The explicit incompressible mesh-free method is then used to simulate open channel flows,including weir flow,hydraulic jump,and flow over an obstacle.In the simulations,velocity distribution and flow pattern are examined.The simulated results are compared to available experimental measurements and other numerical results.There is a good agreement between the simulated results and the experimental measurements.It shows that the explicit incompressible mesh-free method can reproduce the flow characteristics in the open channel flows.
基金supported by the National Natural Science Foundation of China(No.10902038)
文摘A mesh-free method based on local Petrov-Galerkin formulation is presented to solve dynamic impact problems of hyperelastic material.In the present method,a simple Heaviside test function is chosen for simplifying domain integrals.Trial function is constructed by using a radial basis function (RBF) coupled with a polynomial basis function,in which the shape function possesses the kronecker delta function property.So,additional treatment is not required for imposing essential boundary conditions.Governing equations of impact problems are established and solved node by node by using an explicit time integration algorithm in a local domain,which is very similar to that of the collocation method except that numerical integration can be implemented over local domain in the present method.Numerical results for several examples show that the present method performs well in dealing with the dynamic impact problem of hyperelastic material.
基金partially supported by the award PSC-CUNY64335-0052,jointly funded by The Professional Staff Congress and The City University of New York。
文摘This paper is dedicated to the expansion of the framework of general interpolant observables introduced by Azouani,Olson,and Titi for continuous data assimilation of nonlinear partial differential equations.The main feature of this expanded framework is its mesh-free aspect,which allows the observational data itself to dictate the subdivision of the domain via partition of unity in the spirit of the so-called Partition of Unity Method by Babuska and Melenk.As an application of this framework,we consider a nudging-based scheme for data assimilation applied to the context of the two-dimensional Navier-Stokes equations as a paradigmatic example and establish convergence to the reference solution in all higher-order Sobolev topologies in a periodic,mean-free setting.The convergence analysis also makes use of absorbing ball bounds in higherorder Sobolev norms,for which explicit bounds appear to be available in the literature only up to H^(2);such bounds are additionally proved for all integer levels of Sobolev regularity above H^(2).
基金supported by the National Natural Science Foundation of China (No.12172154)the 111 Project (No.B14044)+1 种基金the Natural Science Foundation of Gansu Province (No.23JRRA1035)the Natural Science Foundation of Anhui University of Finance and Economics (No.ACKYC20043).
文摘In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.