This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by...This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by a finite difference scheme,and thus the governing equation of transient heat transfer is transformed into a non-homogeneous modified Helmholtz equation.Secondly,the solution of the non-homogeneous modified Helmholtz equation is decomposed into a particular solution and a homogeneous solution.And then,the DRM and LKM are used to solve the particular solution of the non-homogeneous equation and the homogeneous solution of the modified Helmholtz equation,respectively.The LKM is a recently proposed local radial basis function collocationmethod with themerits of being simple,accurate,and free ofmesh and integration.Compared with the traditional domain-type and boundary-type schemes,the present coupling algorithm could be treated as a really good alternative for the analysis of transient heat conduction on high-dimensional and complicated domains.Numerical experiments,including two-and three-dimensional heat transfer models,demonstrated the effectiveness and accuracy of the new methodology.展开更多
In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is pres...In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.展开更多
Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the movi...Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.展开更多
Multiresolution analysis of wavelet theory can give an effective way to describe the information at various levels of approximations or different resolutions, based on spline wavelet analysis,so weight function is ort...Multiresolution analysis of wavelet theory can give an effective way to describe the information at various levels of approximations or different resolutions, based on spline wavelet analysis,so weight function is orthonormally projected onto a sequence of closed spline subspaces, and is viewed at various levels of approximations or different resolutions. Now, the useful new way to research weight function is found, and the numerical result is given.展开更多
A new direct method for solving unsymmetrical sparse linear systems(USLS) arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin (MLPG) method ...A new direct method for solving unsymmetrical sparse linear systems(USLS) arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin (MLPG) method need to solve large USLS. The proposed solution method for unsymmetrical case performs factorization processes symmetrically on the upper and lower triangular portion of matrix, which differs from previous work based on general unsymmetrical process, and attains higher performance. It is shown that the solution algorithm for USLS can be simply derived from the existing approaches for the symmetrical case. The new matrix factorization algorithm in our method can be implemented easily by modifying a standard JKI symmetrical matrix factorization code. Multi-blocked out-of-core strategies were also developed to expand the solution scale. The approach convincingly increases the speed of the solution process, which is demonstrated with the numerical tests.展开更多
Numerical quadrature is an important ingredient of Galerkin meshless methods. A new numerical quadrature technique, partition of unity quadrature (PUQ),for Galerkin meshless methods was presented. The technique is b...Numerical quadrature is an important ingredient of Galerkin meshless methods. A new numerical quadrature technique, partition of unity quadrature (PUQ),for Galerkin meshless methods was presented. The technique is based on finite covering and partition of unity. There is no need to decompose the physical domain into small cell. It possesses remarkable integration accuracy. Using Element-free Galerkin methods as example, Galerkin meshless methods based on PUQ were studied in detail. Meshing is always not required in the procedure of constitution of approximate function or numerical quadrature, so Galerkin meshless methods based on PUQ are “truly” meshless methods.展开更多
The element-free Galerkin (EFG) method for numerically solving the compound Korteweg-de Vries-Burgers (KdVB) equation is discussed in this paper. The Galerkin weak form is used to obtain the discrete equation and ...The element-free Galerkin (EFG) method for numerically solving the compound Korteweg-de Vries-Burgers (KdVB) equation is discussed in this paper. The Galerkin weak form is used to obtain the discrete equation and the essential boundary conditions are enforced by the penalty method. The effectiveness of the EFG method of solving the compound Korteweg-de Vries-Burgers (KdVB) equation is illustrated by three numerical examples.展开更多
An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with ...An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with smoothed ones, which are projected with Taylor series. Voronoi Cells are introduced not only for integration of potential energy but also for guidance of refinement. New nodes are placed within those cells with high estimated error. At the end of the paper, two numerical examples with severe stress gradient are analyzed. Through adaptive analysis accurate results are obtained at critical subdomains, which validates the efficiency of the method.展开更多
A meshless method integrated with linear elastic fracture mechanics(LEFM)is presented for 2D mixed-mode crack propagation analysis.The domain is divided automatically into sub-domains based on Voronoi cells,which are ...A meshless method integrated with linear elastic fracture mechanics(LEFM)is presented for 2D mixed-mode crack propagation analysis.The domain is divided automatically into sub-domains based on Voronoi cells,which are used for quadrature for the potential energy. The continuous crack propagation is simulated with an incremental crack-extension method which assumes a piecewise linear discretization of the unknown crack path.For each increment of the crack extension,the meshless method is applied to carry out a stress analysis of the cracked structure.The J-integral,which can be decomposed into mode Ⅰ and mode Ⅱ for mixed-mode crack,is used for the evaluation of the stress intensity factors(SIFs).The crack-propagation direction,predicted on an incremental basis, is computed by a criterion defined in terms of the SIFs. The flowchart of the proposed procedure is presented and two numerical problems are analyzed with this method.The meshless results agree well with the experimental ones,which validates the accuracy and efficiency of the method.展开更多
In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equi...In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equivalent functional of two-dimensional advection-diffusion problems is formed,the variation method is used to obtain the equation system,and the penalty method is employed to impose the essential boundary conditions.The difference method for twopoint boundary value problems is used to obtain the discrete equations.Then the corresponding formulas of the ICVMM for advection-diffusion problems are presented.Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper.It is shown that ICVMM is very effective for advection-diffusion problems,and has a good convergent character,accuracy,and computational efficiency.展开更多
Based on the improved version of the meshless singular boundary method(ISBM)in multi domain(MD),a numerical method is proposed in this paper to study the interaction of submerged permeable breakwaters and regular wave...Based on the improved version of the meshless singular boundary method(ISBM)in multi domain(MD),a numerical method is proposed in this paper to study the interaction of submerged permeable breakwaters and regular waves at normal incidence.To account for fluid flow inside the porous breakwaters,the conventional model of Sollitt and Cross for porous media is adopted.Both single and dual trapezoidal breakwaters are examined.The physical problem is formulated in the context of the linear potential wave theory.The domain decomposition method(DDM)is employed,in which the full computational domain is decomposed into separate domains,that is,the fluid domain and the domains of the breakwaters.Respectively,appropriate mixed type boundary and continuity conditions are applied for each subdomain and at the interfaces between domains.The solution is approximated in each subdomain by the ISBM.The discretized algebraic equations are combined,resulting in an overdetermined full system that is solved using a least-square solution procedure.The numerical results are presented in terms of the hydrodynamic quantities of reflection,transmission,and wave-energy dissipation.The relevance of the results of the present numerical procedure is first validated against data of previous studies,and then selected computations are discussed for various structural conditions.The proposed method is demonstrated to be highly accurate and computationally efficient.展开更多
This paper presents a meshless method for the nonlinear generalized regularized long wave (GRLW) equation based on the moving least-squares approximation. The nonlinear discrete scheme of the GRLW equation is obtain...This paper presents a meshless method for the nonlinear generalized regularized long wave (GRLW) equation based on the moving least-squares approximation. The nonlinear discrete scheme of the GRLW equation is obtained and is solved using the iteration method. A theorem on the convergence of the iterative process is presented and proved using theorems of the infinity norm. Compared with numerical methods based on mesh, the meshless method for the GRLW equation only requires the.scattered nodes instead of meshing the domain of the problem. Some examples, such as the propagation of single soliton and the interaction of two solitary waves, are given to show the effectiveness of the meshless method.展开更多
It is weN-known that the standard Galerkin is not ideally suited to deal with the spatial discretization of convection-dominated problems. In this paper, several techniques are proposed to overcome the instabilitY iss...It is weN-known that the standard Galerkin is not ideally suited to deal with the spatial discretization of convection-dominated problems. In this paper, several techniques are proposed to overcome the instabilitY issues in convection-dominated problems in the simulation with a meshless method. These stable techniques included nodal refinement, enlargement of the nodal influence domain, full upwind meshless technique and adaptive upwind meshless technique. Numerical results for sample problems show that these techniques are effective in solving convection-dominated problems, and the adaptive upwind meshless technique is the most effective method of all.展开更多
Numerical solutions of the second-order one-dimensional hyperbolic telegraph equations are presented using the radial basis functions.The purpose of this paper is to propose a simple novel direct meshless scheme for s...Numerical solutions of the second-order one-dimensional hyperbolic telegraph equations are presented using the radial basis functions.The purpose of this paper is to propose a simple novel direct meshless scheme for solving hyperbolic telegraph equations.This is fulfilled by considering time variable as normal space variable.Under this scheme,there is no need to remove time-dependent variable during the whole solution process.Since the numerical solution accuracy depends on the condition of coefficient matrix derived from the radial basis function method.We propose a simple shifted domain method,which can avoid the full-coefficient interpolation matrix easily.Numerical experiments performed with the proposed numerical scheme for several second-order hyperbolic telegraph equations are presented with some discussions.展开更多
The reservoir is the networked rock skeleton of an oil and gas trap,as well as the generic term for the fluid contained within pore fractures and karst caves.Heterogeneity and a complex internal pore structure charact...The reservoir is the networked rock skeleton of an oil and gas trap,as well as the generic term for the fluid contained within pore fractures and karst caves.Heterogeneity and a complex internal pore structure characterize the reservoir rock.By introducing the fractal permeability formula,this paper establishes a fractal mathematical model of oil-water two-phase flow in an oil reservoir with heterogeneity characteristics and numerically solves the mathematical model using the weighted least squares meshless method.Additionally,the method’s correctness is verified by comparison to the exact solution.The numerical results demonstrate that the fractal oil-water two-phase flow mathematical model developed using the meshless method is capable of more accurately and efficiently describing the flow characteristics of the oil-water two-phase migration process.In comparison to the conventional numerical model,this method achieves a greater degree of convergence and stability.This paper examines the effect of varying the initial viscosity of the oil,the initial formation pressure,and the production and injection ratios on daily oil production per well,water cut in the block,and accumulated oil in the block.For 10 and 60 cp initial crude oil viscosities,the water cut can be 0.62 and 0.80,with 3100 and 1900 m^(3)cumulative oil production.Initial pressures have little effect on production.In this case,the daily oil production of well PRO1 is 1.7 m^(3)at 7 and 10 MPa initial pressure.Block cumulative oil production is 3465.4 and 2149.9m^(3)when the production injection ratio is 1.4 and 0.8.The two-phase meshless method described in this paper is essential for a rational and effective study of production dynamics patterns in complex reservoirs and the development of reservoir simulations of oil-water flow in heterogeneous reservoirs.展开更多
In this paper a meshless method of lines is proposed for the numerical solution of time-dependent nonlinear coupled partial differential equations. Contrary to mesh oriented methods of lines using the finite-differenc...In this paper a meshless method of lines is proposed for the numerical solution of time-dependent nonlinear coupled partial differential equations. Contrary to mesh oriented methods of lines using the finite-difference and finite element methods to approximate spatial derivatives, this new technique does not require a mesh in the problem domain, and a set of scattered nodes provided by initial data is required for the solution of the problem using some radial basis functions. Accuracy of the method is assessed in terms of the error norms L2, L∞ and the three invariants C1, C2, C3. Numerical experiments are performed to demonstrate the accuracy and easy implementation of this method for the three classes of time-dependent nonlinear coupled partial differential equations.展开更多
In the present work, the improved version of the meshless singular boundary method(ISBM) is developed for analyzing the performance of bottom standing submerged permeable breakwaters in regular normally incident waves...In the present work, the improved version of the meshless singular boundary method(ISBM) is developed for analyzing the performance of bottom standing submerged permeable breakwaters in regular normally incident waves and in the proximity of a vertical wall. Both single and dual prismatic breakwaters of rectangular and trapezoidal forms are examined. The physical problem is cast in terms of the Laplace equation governing an irrotational flow and incompressible fluid motion with appropriate mixed type boundary conditions, and solved numerically using the ISBM. To model the permeability of the breakwaters fully absorbing boundary conditions are assumed. Numerical results are presented in terms of hydrodynamic quantities of the reflection coefficients. These are firstly validated against the results of a multi-domain boundary element method(BEM) developed independently for a previous study. The agreement between the results of the two methods is excellent. The coefficients of reflection are then computed and discussed for a variety of structural conditions including the breakwaters height, width, spacing, and absorbing permeability. Effects of the proximity of the vertical plane wall are also investigated. The breakwater's width is found to have only marginal effects compared with its height. Permeability tends to decrease the minimum reflections. These coefficients show periodic variations with the spacing relative to the wavelength. Trapezoidal breakwaters are found to be more cost-effective than the rectangular breakwaters. Dual breakwater systems are confirmed to perform much better than single structures.展开更多
To avoid mesh distortion and iterative remeshing in mesh-based numerical analysis,a meshless approach based on element free Galerkin (EFG) method is applied to the metal forming analysis of ring compression. Discrete ...To avoid mesh distortion and iterative remeshing in mesh-based numerical analysis,a meshless approach based on element free Galerkin (EFG) method is applied to the metal forming analysis of ring compression. Discrete equations are formulated upon the moving least-squares (MLS) approximation and modified Markov variational principles for rigid-plastic/ rigid-viscoplastic (RP/RVP) material models. The penalty function is used for the incompressible condition without volumetric locking. Based on the axisymmetric mechanical model,ring tests with different friction coefficients are studied. The deformed nodal configurations and shaded contours of equivalent strains are shown by developed meshless post processor. The comparison of meshless and finite element (FE) results validates the feasibility and accuracy for meshless method to simulate metal forming process.展开更多
In this work, an algorithm based on method of lines coupled with radial basis functions namely meshless method of lines (MMOL) is presented for the numerical solution of Kawahara, modified Kawahara and KdV Kawahara eq...In this work, an algorithm based on method of lines coupled with radial basis functions namely meshless method of lines (MMOL) is presented for the numerical solution of Kawahara, modified Kawahara and KdV Kawahara equations. The motion of a single solitary wave, interaction of two and three solitons and the phenomena of wave generation is discussed. The results are compared with the exact solution and with the results in the relevant literature to show the efficiency of the method.展开更多
基金supported by the NationalNatural Science Foundation of China (No.11802151)the Natural Science Foundation of Shandong Province of China (No.ZR2019BA008)the China Postdoctoral Science Foundation (No.2019M652315).
文摘This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by a finite difference scheme,and thus the governing equation of transient heat transfer is transformed into a non-homogeneous modified Helmholtz equation.Secondly,the solution of the non-homogeneous modified Helmholtz equation is decomposed into a particular solution and a homogeneous solution.And then,the DRM and LKM are used to solve the particular solution of the non-homogeneous equation and the homogeneous solution of the modified Helmholtz equation,respectively.The LKM is a recently proposed local radial basis function collocationmethod with themerits of being simple,accurate,and free ofmesh and integration.Compared with the traditional domain-type and boundary-type schemes,the present coupling algorithm could be treated as a really good alternative for the analysis of transient heat conduction on high-dimensional and complicated domains.Numerical experiments,including two-and three-dimensional heat transfer models,demonstrated the effectiveness and accuracy of the new methodology.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Shanghai Leading Academic Discipline Project,China(Grant No.S30106)the Innovation Fund for Graduate Student of Shanghai University of China (Grant No.SHUCX120125)
文摘In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072117)the Natural Science Foundation of Ningbo City,China(GrantNo.2013A610103)+2 种基金the Natural Science Foundation of Zhejiang Province,China(Grant No.Y6090131)the Disciplinary Project of Ningbo City,China(GrantNo.SZXL1067)the K.C.Wong Magna Fund in Ningbo University,China
文摘Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.
基金theNationalNaturalScienceFoundationofChina (No .50 40 90 0 8)
文摘Multiresolution analysis of wavelet theory can give an effective way to describe the information at various levels of approximations or different resolutions, based on spline wavelet analysis,so weight function is orthonormally projected onto a sequence of closed spline subspaces, and is viewed at various levels of approximations or different resolutions. Now, the useful new way to research weight function is found, and the numerical result is given.
基金Project supported by the National Natural Science Foundation of China (Nos. 10232040, 10572002 and 10572003)
文摘A new direct method for solving unsymmetrical sparse linear systems(USLS) arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin (MLPG) method need to solve large USLS. The proposed solution method for unsymmetrical case performs factorization processes symmetrically on the upper and lower triangular portion of matrix, which differs from previous work based on general unsymmetrical process, and attains higher performance. It is shown that the solution algorithm for USLS can be simply derived from the existing approaches for the symmetrical case. The new matrix factorization algorithm in our method can be implemented easily by modifying a standard JKI symmetrical matrix factorization code. Multi-blocked out-of-core strategies were also developed to expand the solution scale. The approach convincingly increases the speed of the solution process, which is demonstrated with the numerical tests.
文摘Numerical quadrature is an important ingredient of Galerkin meshless methods. A new numerical quadrature technique, partition of unity quadrature (PUQ),for Galerkin meshless methods was presented. The technique is based on finite covering and partition of unity. There is no need to decompose the physical domain into small cell. It possesses remarkable integration accuracy. Using Element-free Galerkin methods as example, Galerkin meshless methods based on PUQ were studied in detail. Meshing is always not required in the procedure of constitution of approximate function or numerical quadrature, so Galerkin meshless methods based on PUQ are “truly” meshless methods.
基金Project supported by the National Natural Science Foundation of China (Grant No.10871124)the Natural Science Foundation of Zhejiang Province of China (Grant No.Y6110007)
文摘The element-free Galerkin (EFG) method for numerically solving the compound Korteweg-de Vries-Burgers (KdVB) equation is discussed in this paper. The Galerkin weak form is used to obtain the discrete equation and the essential boundary conditions are enforced by the penalty method. The effectiveness of the EFG method of solving the compound Korteweg-de Vries-Burgers (KdVB) equation is illustrated by three numerical examples.
基金Project supported by the National Natural Science Foundation of China (No. 50175060).
文摘An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with smoothed ones, which are projected with Taylor series. Voronoi Cells are introduced not only for integration of potential energy but also for guidance of refinement. New nodes are placed within those cells with high estimated error. At the end of the paper, two numerical examples with severe stress gradient are analyzed. Through adaptive analysis accurate results are obtained at critical subdomains, which validates the efficiency of the method.
基金Project supported by the National Natural Science Foundation of China(Nos.59825117 and 50175060).
文摘A meshless method integrated with linear elastic fracture mechanics(LEFM)is presented for 2D mixed-mode crack propagation analysis.The domain is divided automatically into sub-domains based on Voronoi cells,which are used for quadrature for the potential energy. The continuous crack propagation is simulated with an incremental crack-extension method which assumes a piecewise linear discretization of the unknown crack path.For each increment of the crack extension,the meshless method is applied to carry out a stress analysis of the cracked structure.The J-integral,which can be decomposed into mode Ⅰ and mode Ⅱ for mixed-mode crack,is used for the evaluation of the stress intensity factors(SIFs).The crack-propagation direction,predicted on an incremental basis, is computed by a criterion defined in terms of the SIFs. The flowchart of the proposed procedure is presented and two numerical problems are analyzed with this method.The meshless results agree well with the experimental ones,which validates the accuracy and efficiency of the method.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project,China(Grant No. S30106)the Innovation Fund for Graduate Student of Shanghai University,China (Grant No. SHUCX120125)
文摘In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equivalent functional of two-dimensional advection-diffusion problems is formed,the variation method is used to obtain the equation system,and the penalty method is employed to impose the essential boundary conditions.The difference method for twopoint boundary value problems is used to obtain the discrete equations.Then the corresponding formulas of the ICVMM for advection-diffusion problems are presented.Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper.It is shown that ICVMM is very effective for advection-diffusion problems,and has a good convergent character,accuracy,and computational efficiency.
基金the Ministry of Higher Edu-cation and Scientific Research of Algeria(grant PRFU number A01L06UN310220200002).
文摘Based on the improved version of the meshless singular boundary method(ISBM)in multi domain(MD),a numerical method is proposed in this paper to study the interaction of submerged permeable breakwaters and regular waves at normal incidence.To account for fluid flow inside the porous breakwaters,the conventional model of Sollitt and Cross for porous media is adopted.Both single and dual trapezoidal breakwaters are examined.The physical problem is formulated in the context of the linear potential wave theory.The domain decomposition method(DDM)is employed,in which the full computational domain is decomposed into separate domains,that is,the fluid domain and the domains of the breakwaters.Respectively,appropriate mixed type boundary and continuity conditions are applied for each subdomain and at the interfaces between domains.The solution is approximated in each subdomain by the ISBM.The discretized algebraic equations are combined,resulting in an overdetermined full system that is solved using a least-square solution procedure.The numerical results are presented in terms of the hydrodynamic quantities of reflection,transmission,and wave-energy dissipation.The relevance of the results of the present numerical procedure is first validated against data of previous studies,and then selected computations are discussed for various structural conditions.The proposed method is demonstrated to be highly accurate and computationally efficient.
基金supported by the National Natural Science Foundation of China (Grant No. 10871124)the Innovation Program of the Shanghai Municipal Education Commission,China (Grant No. 09ZZ99)
文摘This paper presents a meshless method for the nonlinear generalized regularized long wave (GRLW) equation based on the moving least-squares approximation. The nonlinear discrete scheme of the GRLW equation is obtained and is solved using the iteration method. A theorem on the convergence of the iterative process is presented and proved using theorems of the infinity norm. Compared with numerical methods based on mesh, the meshless method for the GRLW equation only requires the.scattered nodes instead of meshing the domain of the problem. Some examples, such as the propagation of single soliton and the interaction of two solitary waves, are given to show the effectiveness of the meshless method.
基金the National Natural Science Foundation of China(No.10590353)theNatural Science Foundation of Shaanxi Province of China(No.2005A16)
文摘It is weN-known that the standard Galerkin is not ideally suited to deal with the spatial discretization of convection-dominated problems. In this paper, several techniques are proposed to overcome the instabilitY issues in convection-dominated problems in the simulation with a meshless method. These stable techniques included nodal refinement, enlargement of the nodal influence domain, full upwind meshless technique and adaptive upwind meshless technique. Numerical results for sample problems show that these techniques are effective in solving convection-dominated problems, and the adaptive upwind meshless technique is the most effective method of all.
基金The first author is supported by the Natural Science Foundation of Anhui Province(Project No.1908085QA09)the University Natural Science Research Project of Anhui Province(Project Nos.KJ2019A0591&KJ2020B06)。
文摘Numerical solutions of the second-order one-dimensional hyperbolic telegraph equations are presented using the radial basis functions.The purpose of this paper is to propose a simple novel direct meshless scheme for solving hyperbolic telegraph equations.This is fulfilled by considering time variable as normal space variable.Under this scheme,there is no need to remove time-dependent variable during the whole solution process.Since the numerical solution accuracy depends on the condition of coefficient matrix derived from the radial basis function method.We propose a simple shifted domain method,which can avoid the full-coefficient interpolation matrix easily.Numerical experiments performed with the proposed numerical scheme for several second-order hyperbolic telegraph equations are presented with some discussions.
基金The National Natural Science Foundation of China(Nos.51874044,51922007).
文摘The reservoir is the networked rock skeleton of an oil and gas trap,as well as the generic term for the fluid contained within pore fractures and karst caves.Heterogeneity and a complex internal pore structure characterize the reservoir rock.By introducing the fractal permeability formula,this paper establishes a fractal mathematical model of oil-water two-phase flow in an oil reservoir with heterogeneity characteristics and numerically solves the mathematical model using the weighted least squares meshless method.Additionally,the method’s correctness is verified by comparison to the exact solution.The numerical results demonstrate that the fractal oil-water two-phase flow mathematical model developed using the meshless method is capable of more accurately and efficiently describing the flow characteristics of the oil-water two-phase migration process.In comparison to the conventional numerical model,this method achieves a greater degree of convergence and stability.This paper examines the effect of varying the initial viscosity of the oil,the initial formation pressure,and the production and injection ratios on daily oil production per well,water cut in the block,and accumulated oil in the block.For 10 and 60 cp initial crude oil viscosities,the water cut can be 0.62 and 0.80,with 3100 and 1900 m^(3)cumulative oil production.Initial pressures have little effect on production.In this case,the daily oil production of well PRO1 is 1.7 m^(3)at 7 and 10 MPa initial pressure.Block cumulative oil production is 3465.4 and 2149.9m^(3)when the production injection ratio is 1.4 and 0.8.The two-phase meshless method described in this paper is essential for a rational and effective study of production dynamics patterns in complex reservoirs and the development of reservoir simulations of oil-water flow in heterogeneous reservoirs.
文摘In this paper a meshless method of lines is proposed for the numerical solution of time-dependent nonlinear coupled partial differential equations. Contrary to mesh oriented methods of lines using the finite-difference and finite element methods to approximate spatial derivatives, this new technique does not require a mesh in the problem domain, and a set of scattered nodes provided by initial data is required for the solution of the problem using some radial basis functions. Accuracy of the method is assessed in terms of the error norms L2, L∞ and the three invariants C1, C2, C3. Numerical experiments are performed to demonstrate the accuracy and easy implementation of this method for the three classes of time-dependent nonlinear coupled partial differential equations.
基金financially supported by the Direction Général des Enseignements et de la Formation Supérieure of Algeria(Grant CNEPRU No.G0301920140029)
文摘In the present work, the improved version of the meshless singular boundary method(ISBM) is developed for analyzing the performance of bottom standing submerged permeable breakwaters in regular normally incident waves and in the proximity of a vertical wall. Both single and dual prismatic breakwaters of rectangular and trapezoidal forms are examined. The physical problem is cast in terms of the Laplace equation governing an irrotational flow and incompressible fluid motion with appropriate mixed type boundary conditions, and solved numerically using the ISBM. To model the permeability of the breakwaters fully absorbing boundary conditions are assumed. Numerical results are presented in terms of hydrodynamic quantities of the reflection coefficients. These are firstly validated against the results of a multi-domain boundary element method(BEM) developed independently for a previous study. The agreement between the results of the two methods is excellent. The coefficients of reflection are then computed and discussed for a variety of structural conditions including the breakwaters height, width, spacing, and absorbing permeability. Effects of the proximity of the vertical plane wall are also investigated. The breakwater's width is found to have only marginal effects compared with its height. Permeability tends to decrease the minimum reflections. These coefficients show periodic variations with the spacing relative to the wavelength. Trapezoidal breakwaters are found to be more cost-effective than the rectangular breakwaters. Dual breakwater systems are confirmed to perform much better than single structures.
文摘To avoid mesh distortion and iterative remeshing in mesh-based numerical analysis,a meshless approach based on element free Galerkin (EFG) method is applied to the metal forming analysis of ring compression. Discrete equations are formulated upon the moving least-squares (MLS) approximation and modified Markov variational principles for rigid-plastic/ rigid-viscoplastic (RP/RVP) material models. The penalty function is used for the incompressible condition without volumetric locking. Based on the axisymmetric mechanical model,ring tests with different friction coefficients are studied. The deformed nodal configurations and shaded contours of equivalent strains are shown by developed meshless post processor. The comparison of meshless and finite element (FE) results validates the feasibility and accuracy for meshless method to simulate metal forming process.
文摘In this work, an algorithm based on method of lines coupled with radial basis functions namely meshless method of lines (MMOL) is presented for the numerical solution of Kawahara, modified Kawahara and KdV Kawahara equations. The motion of a single solitary wave, interaction of two and three solitons and the phenomena of wave generation is discussed. The results are compared with the exact solution and with the results in the relevant literature to show the efficiency of the method.