It has been evidenced that shallow gas hydrate resources are abundant in deep oceans worldwide.Their geological back-ground,occurrence,and other characteristics differ significantly from deep-seated hydrates.Because o...It has been evidenced that shallow gas hydrate resources are abundant in deep oceans worldwide.Their geological back-ground,occurrence,and other characteristics differ significantly from deep-seated hydrates.Because of the high risk of well construction and low production efficiency,they are difficult to be recovered by using conventional oil production methods.As a result,this paper proposes an alternative design based on a combination of radial drilling,heat injection,and backfilling methods.Multi-branch holes are used to penetrate shallow gas hydrate reservoirs to expand the depressurization area,and heat injection is utilized as a supplement to improve gas production.Geotechnical information collected from an investigation site close to the offshore production well in the South China Sea is used to assess the essential components of this plan,including well construction stability and gas production behavior.It demonstrates that the hydraulic fracturing of the 60mbsf overburden layer can be prevented by regulating the drilling fluid densities.However,the traditional well structure is unstable,and the suction anchor is advised for better mechanical performance.The gas produc-tion rate can be significantly increased by combining hot water injection and depressurization methods.Additionally,the suitable produc-tion equipment already in use is discussed.展开更多
BACKGROUND Congenital duodenal webs are rare in adults and can lead to various symptoms such as nausea,vomiting,and postprandial fullness.The treatment for this disease is mostly surgical.Endoscopic treatment techniqu...BACKGROUND Congenital duodenal webs are rare in adults and can lead to various symptoms such as nausea,vomiting,and postprandial fullness.The treatment for this disease is mostly surgical.Endoscopic treatment techniques have been developed and attempted for this disease.Endoscopic radial incision and cutting(RIC)techniques are reportedly very effective in benign anastomotic stricture.This case report highlights the effectiveness and safety of endoscopic RIC as a minimally invasive treatment for adult congenital duodenal webs.CASE SUMMARY A 23-year-old female patient with indigestion was referred to a tertiary hospital.The patient complained of postprandial fullness in the epigastric region.Previous physical examinations or blood tests indicated no abnormalities.Computed tomography revealed an eccentric broad-based delayed-enhancing mass-like lesion in the second portion of the duodenum.Endoscopy showed an enlarged gastric cavity and a significantly dilated duodenal bulb;a very small hole was observed in the distal part of the second portion,and scope passage was not possible.Gastrografin upper gastrointestinal series was performed,revealing an intraduodenal barium contrast-filled sac with a curvilinear narrow radiolucent rim,a typical"windsock"sign.Endoscopic RIC was performed on the duodenal web.The patient recovered uneventfully.Follow-up endoscopy showed a patent duodenal lumen without any residual stenosis.The patient reported complete resolution of symptoms at the 18-month follow-up.CONCLUSION Endoscopic RIC may be an effective treatment for congenital duodenal webs in adults.展开更多
Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with...Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.展开更多
This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structu...This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.展开更多
The paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral ...The paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral formulation of the governing equation.Due to the existence of thermal stress,the domain integral term appears in the boundary integral equation.We resolve this problem by incorporating radial integration method into IGABEM which converts the domain integral to the boundary integral.In this way,IGABEM can maintain its advantages in dimensionality reduction and more importantly,seamless integration of CAD and numerical analysis based on boundary representation.The algorithm is verified by numerical examples.展开更多
With the consideration of volume constraint of launch vehicle and trafficability of rover vehicle on lunar regolith terrain, a new design of radially deployable wheel is presented. For the purpose of achieving the mes...With the consideration of volume constraint of launch vehicle and trafficability of rover vehicle on lunar regolith terrain, a new design of radially deployable wheel is presented. For the purpose of achieving the meso-mechanics and dynamical behavior of lunar soil particles as well as macro-parameters of tractive performance for radially deployable wheel, the interaction between two types of wheel configurations and lunar soil particles is analyzed by means of discrete element method. The network of contact forces, the displacement vector chart, and the deformation of lunar soil beneath wheels are plotted. The equations of soil thrust, motion resistance, drawbar pull and driven torque are derived in granular scale based on the coordinates transformation and algebraic summation. The calculated results show that there is sufficient traction for both 6-split and 12-split radially deployable wheels with 304 mm outspread diameter to negotiate lunar regolith terrain specified here; the value of drawbar pull enhances with the increase of split number of radially deployable wheel, however, the required driven torque increases simultaneously, therefore, the tractive efficiency decreases.展开更多
The smoothing thin plate spline (STPS) interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the ...The smoothing thin plate spline (STPS) interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method (FEM) as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline (TPS) radial basis functions.展开更多
In this paper, Radial point collocation method (RPCM), a kind of meshfree method, is applied to solve convectiondiffusion problem. The main feature of this approach is to use the interpolation schemes in local suppo...In this paper, Radial point collocation method (RPCM), a kind of meshfree method, is applied to solve convectiondiffusion problem. The main feature of this approach is to use the interpolation schemes in local supported domains based on radial basis functions. As a result, this method is local and hence the system matrix is banded which is very attractive for practical engineering problems. In the numerical examination, RPCM is applied to solve non-linear convection-diffusion 2D Burgers equations. The results obtained by RPCM demonstrate the accuracy and efficiency of the proposed method for solving transient fluid dynamic problems. A fictitious point scheme is adopted to improve the solution accuracy while Neumann boundary conditions exist. The meshfree feature of the nresent method is verv attractive in solving comnutational fluid nroblems.展开更多
Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and the...Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.展开更多
Investigations into the magnetohydrodynamics of viscous fluids have become more important in recent years,owing to their practical significance and numerous applications in astro-physical and geo-physical phenomena.In...Investigations into the magnetohydrodynamics of viscous fluids have become more important in recent years,owing to their practical significance and numerous applications in astro-physical and geo-physical phenomena.In this paper,the radial base function was utilized to answer fractional equation associated with fluid flow passing through two parallel flat plates with a magnetic field.The magnetohydrodynamics coupled stress fluid flows between two parallel plates,with the bottom plate being stationary and the top plate moving at a persistent velocity.We compared the radial basis function approach to the numerical method(fourth-order Range-Kutta)in order to verify its validity.The findings demonstrated that the discrepancy between these two techniques is quite negligible,indicating that this method is very reliable.The impact of the magnetic field parameter and Reynolds number on the velocity distribution perpendicular to the fluid flow direction is illustrated.Eventually,the velocity parameter is compared for diverse conditionsα,Reynolds and position(y),the maximum of which occurs atα=0.4.Also,the maximum velocity values occur inα=0.4 and Re=1000 and the concavity of the graph is less forα=0.8.展开更多
On the basis of similar structure of solutions of ordinary differential equation (ODE) boundary value problem, the similar construction method was put forward by solving problems of fluid flow in porous media through ...On the basis of similar structure of solutions of ordinary differential equation (ODE) boundary value problem, the similar construction method was put forward by solving problems of fluid flow in porous media through the homogeneous reservoir. It is indicate that the pressure distribution of dimensionless reservoir and bottom hole in Laplace space, which take on the radial flow, also shows similar structure, and the internal relationship between the above solutions were illustrated in detail.展开更多
In order to calculate the pressure distribution of radial grooved thrust bearing, analytical and numerical methods were applied respectively. Grooved region and land region were linked by u- sing the mass conservation...In order to calculate the pressure distribution of radial grooved thrust bearing, analytical and numerical methods were applied respectively. Grooved region and land region were linked by u- sing the mass conservations principle at the groove/land boundary in each method. The block-weight approach was implemented to deal with the non-coincidence of mesh and radial groove pattern in nu- merical method. It was observed that the numerical solutions had higher precision as mesh number exceed 70 x 70, and the relaxation iteration of differential scheme presented the fastest convergence speed when relaxation factor was close to 1.94.展开更多
The present work describes the application of the method of fundamental solutions (MFS) along with the analog equation method (AEM) and radial basis function (RBF) approximation for solving the 2D isotropic and ...The present work describes the application of the method of fundamental solutions (MFS) along with the analog equation method (AEM) and radial basis function (RBF) approximation for solving the 2D isotropic and anisotropic Helmholtz problems with different wave numbers. The AEM is used to convert the original governing equation into the classical Poisson's equation, and the MFS and RBF approximations are used to derive the homogeneous and particular solutions, respectively. Finally, the satisfaction of the solution consisting of the homogeneous and particular parts to the related governing equation and boundary conditions can produce a system of linear equations, which can be solved with the singular value decomposition (SVD) technique. In the computation, such crucial factors related to the MFS-RBF as the location of the virtual boundary, the differential and integrating strategies, and the variation of shape parameters in multi-quadric (MQ) are fully analyzed to provide useful reference.展开更多
Based on the local discontinuous Galerkin methods for time-dependent convection-diffusion systems newly developed by Corkburn and Shu,according to the form of the generalized convection-diffusion equations which model...Based on the local discontinuous Galerkin methods for time-dependent convection-diffusion systems newly developed by Corkburn and Shu,according to the form of the generalized convection-diffusion equations which model the radial porous flow with dispersion and adsorption,a local discontinuous Galerkin method for radial porous flow with dispersion and adsorption was developed,a high order accurary new scheme for radial porous flow is obtained.The presented method was applied to the numerical tests of two cases of radial porous,i.e., the convection-dispersion flow and the convection-dispersion-adsorption flow,the corresponding parts of the numerical results are in good agreement with the published solutions,so the presented method is reliable.Reckoning of the computational cost also shows that the method is practicable.展开更多
基金financially supported by the Natural Science Foundation of Shandong Province(No.ZR202011030013)the National Natural Science Foundation of China(No.41976205)+1 种基金the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(No.2021QNLM020002)the China Geological Survey Program(No.DD20221704).
文摘It has been evidenced that shallow gas hydrate resources are abundant in deep oceans worldwide.Their geological back-ground,occurrence,and other characteristics differ significantly from deep-seated hydrates.Because of the high risk of well construction and low production efficiency,they are difficult to be recovered by using conventional oil production methods.As a result,this paper proposes an alternative design based on a combination of radial drilling,heat injection,and backfilling methods.Multi-branch holes are used to penetrate shallow gas hydrate reservoirs to expand the depressurization area,and heat injection is utilized as a supplement to improve gas production.Geotechnical information collected from an investigation site close to the offshore production well in the South China Sea is used to assess the essential components of this plan,including well construction stability and gas production behavior.It demonstrates that the hydraulic fracturing of the 60mbsf overburden layer can be prevented by regulating the drilling fluid densities.However,the traditional well structure is unstable,and the suction anchor is advised for better mechanical performance.The gas produc-tion rate can be significantly increased by combining hot water injection and depressurization methods.Additionally,the suitable produc-tion equipment already in use is discussed.
文摘BACKGROUND Congenital duodenal webs are rare in adults and can lead to various symptoms such as nausea,vomiting,and postprandial fullness.The treatment for this disease is mostly surgical.Endoscopic treatment techniques have been developed and attempted for this disease.Endoscopic radial incision and cutting(RIC)techniques are reportedly very effective in benign anastomotic stricture.This case report highlights the effectiveness and safety of endoscopic RIC as a minimally invasive treatment for adult congenital duodenal webs.CASE SUMMARY A 23-year-old female patient with indigestion was referred to a tertiary hospital.The patient complained of postprandial fullness in the epigastric region.Previous physical examinations or blood tests indicated no abnormalities.Computed tomography revealed an eccentric broad-based delayed-enhancing mass-like lesion in the second portion of the duodenum.Endoscopy showed an enlarged gastric cavity and a significantly dilated duodenal bulb;a very small hole was observed in the distal part of the second portion,and scope passage was not possible.Gastrografin upper gastrointestinal series was performed,revealing an intraduodenal barium contrast-filled sac with a curvilinear narrow radiolucent rim,a typical"windsock"sign.Endoscopic RIC was performed on the duodenal web.The patient recovered uneventfully.Follow-up endoscopy showed a patent duodenal lumen without any residual stenosis.The patient reported complete resolution of symptoms at the 18-month follow-up.CONCLUSION Endoscopic RIC may be an effective treatment for congenital duodenal webs in adults.
文摘Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.
基金funded by National Natural Science Foundation of China(NSFC)under Grant Nos.11702238,51904202,and 11902212Nanhu Scholars Program for Young Scholars of XYNU.
文摘This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.
基金This study was funded by the National Natural Science Foundation of China(NSFC)(Grant Nos.11702238,51904202 and 11902212)and Nanhu Scholars Program for Young Scholars of XYNU.
文摘The paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral formulation of the governing equation.Due to the existence of thermal stress,the domain integral term appears in the boundary integral equation.We resolve this problem by incorporating radial integration method into IGABEM which converts the domain integral to the boundary integral.In this way,IGABEM can maintain its advantages in dimensionality reduction and more importantly,seamless integration of CAD and numerical analysis based on boundary representation.The algorithm is verified by numerical examples.
基金the Doctoral Program of Higher Education of China(No.20070006012)and Pre-research Project of China Academy of Space Technology.
文摘With the consideration of volume constraint of launch vehicle and trafficability of rover vehicle on lunar regolith terrain, a new design of radially deployable wheel is presented. For the purpose of achieving the meso-mechanics and dynamical behavior of lunar soil particles as well as macro-parameters of tractive performance for radially deployable wheel, the interaction between two types of wheel configurations and lunar soil particles is analyzed by means of discrete element method. The network of contact forces, the displacement vector chart, and the deformation of lunar soil beneath wheels are plotted. The equations of soil thrust, motion resistance, drawbar pull and driven torque are derived in granular scale based on the coordinates transformation and algebraic summation. The calculated results show that there is sufficient traction for both 6-split and 12-split radially deployable wheels with 304 mm outspread diameter to negotiate lunar regolith terrain specified here; the value of drawbar pull enhances with the increase of split number of radially deployable wheel, however, the required driven torque increases simultaneously, therefore, the tractive efficiency decreases.
基金supported by the Key Program of the National Natural Science Foundation of China (Grand No. 51138001)the China-German Cooperation Project (Grand No. GZ566)+1 种基金the Innovative Research Groups Funded by the National Natural Science Foundation of China (Grand No. 51121005)the Special Funds for the Basic Scientific Research Expenses for the Central University (Grant No. DUT13LK16)
文摘The smoothing thin plate spline (STPS) interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method (FEM) as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline (TPS) radial basis functions.
基金Project (No. 10572128) supported by the National Natural ScienceFoundation of China
文摘In this paper, Radial point collocation method (RPCM), a kind of meshfree method, is applied to solve convectiondiffusion problem. The main feature of this approach is to use the interpolation schemes in local supported domains based on radial basis functions. As a result, this method is local and hence the system matrix is banded which is very attractive for practical engineering problems. In the numerical examination, RPCM is applied to solve non-linear convection-diffusion 2D Burgers equations. The results obtained by RPCM demonstrate the accuracy and efficiency of the proposed method for solving transient fluid dynamic problems. A fictitious point scheme is adopted to improve the solution accuracy while Neumann boundary conditions exist. The meshfree feature of the nresent method is verv attractive in solving comnutational fluid nroblems.
基金Project(2010CB732103)supported by the National Basic Research Program of ChinaProject(51179092)supported by the National Natural Science Foundation of ChinaProject(2012-KY-02)supported by the State Key Laboratory of Hydroscience and Engineering,China
文摘Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.
文摘Investigations into the magnetohydrodynamics of viscous fluids have become more important in recent years,owing to their practical significance and numerous applications in astro-physical and geo-physical phenomena.In this paper,the radial base function was utilized to answer fractional equation associated with fluid flow passing through two parallel flat plates with a magnetic field.The magnetohydrodynamics coupled stress fluid flows between two parallel plates,with the bottom plate being stationary and the top plate moving at a persistent velocity.We compared the radial basis function approach to the numerical method(fourth-order Range-Kutta)in order to verify its validity.The findings demonstrated that the discrepancy between these two techniques is quite negligible,indicating that this method is very reliable.The impact of the magnetic field parameter and Reynolds number on the velocity distribution perpendicular to the fluid flow direction is illustrated.Eventually,the velocity parameter is compared for diverse conditionsα,Reynolds and position(y),the maximum of which occurs atα=0.4.Also,the maximum velocity values occur inα=0.4 and Re=1000 and the concavity of the graph is less forα=0.8.
文摘On the basis of similar structure of solutions of ordinary differential equation (ODE) boundary value problem, the similar construction method was put forward by solving problems of fluid flow in porous media through the homogeneous reservoir. It is indicate that the pressure distribution of dimensionless reservoir and bottom hole in Laplace space, which take on the radial flow, also shows similar structure, and the internal relationship between the above solutions were illustrated in detail.
基金Supported by the Ministerial Level Foundation(2220060029)
文摘In order to calculate the pressure distribution of radial grooved thrust bearing, analytical and numerical methods were applied respectively. Grooved region and land region were linked by u- sing the mass conservations principle at the groove/land boundary in each method. The block-weight approach was implemented to deal with the non-coincidence of mesh and radial groove pattern in nu- merical method. It was observed that the numerical solutions had higher precision as mesh number exceed 70 x 70, and the relaxation iteration of differential scheme presented the fastest convergence speed when relaxation factor was close to 1.94.
文摘The present work describes the application of the method of fundamental solutions (MFS) along with the analog equation method (AEM) and radial basis function (RBF) approximation for solving the 2D isotropic and anisotropic Helmholtz problems with different wave numbers. The AEM is used to convert the original governing equation into the classical Poisson's equation, and the MFS and RBF approximations are used to derive the homogeneous and particular solutions, respectively. Finally, the satisfaction of the solution consisting of the homogeneous and particular parts to the related governing equation and boundary conditions can produce a system of linear equations, which can be solved with the singular value decomposition (SVD) technique. In the computation, such crucial factors related to the MFS-RBF as the location of the virtual boundary, the differential and integrating strategies, and the variation of shape parameters in multi-quadric (MQ) are fully analyzed to provide useful reference.
文摘Based on the local discontinuous Galerkin methods for time-dependent convection-diffusion systems newly developed by Corkburn and Shu,according to the form of the generalized convection-diffusion equations which model the radial porous flow with dispersion and adsorption,a local discontinuous Galerkin method for radial porous flow with dispersion and adsorption was developed,a high order accurary new scheme for radial porous flow is obtained.The presented method was applied to the numerical tests of two cases of radial porous,i.e., the convection-dispersion flow and the convection-dispersion-adsorption flow,the corresponding parts of the numerical results are in good agreement with the published solutions,so the presented method is reliable.Reckoning of the computational cost also shows that the method is practicable.