In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Mill...In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.展开更多
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.展开更多
By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is propose...By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.展开更多
The EI Nino-southern oscillation (ENSO) is an interannual phenomenon involved in tropical Pacific ocean- atmosphere interactions. In this paper, we develop an asymptotic method of solving the nonlinear equation usin...The EI Nino-southern oscillation (ENSO) is an interannual phenomenon involved in tropical Pacific ocean- atmosphere interactions. In this paper, we develop an asymptotic method of solving the nonlinear equation using the ENSO model. Based on a class of the oscillator of the ENSO model, a approximate solution of the corresponding problem is studied employing the perturbation method.展开更多
In this article, we deal with a class of semilinear elliptic equations which are perturbations of the problems with the critical Hardy-Sobolev exponent. Some existence results are given via an abstract perturbation me...In this article, we deal with a class of semilinear elliptic equations which are perturbations of the problems with the critical Hardy-Sobolev exponent. Some existence results are given via an abstract perturbation method in critical point theory.展开更多
A third-order correction was recently suggested to improve the accuracy of the half-power bandwidth method in estimating the damping of single DOF systems.This paper analyzes the accuracy of the half-power bandwidth m...A third-order correction was recently suggested to improve the accuracy of the half-power bandwidth method in estimating the damping of single DOF systems.This paper analyzes the accuracy of the half-power bandwidth method with the third-order correction in damping estimation for multi-DOF linear systems.Damping ratios in a two-DOF linear system are estimated using its displacement and acceleration frequency response curves,respectively.A wide range of important parameters that characterize the shape of these response curves are taken into account.Results show that the third-order correction may greatly improve the accuracy of the half-power bandwidth method in estimating damping in a two-DOF system.In spite of this,the half-power bandwidth method may significantly overestimate the damping ratios of two-DOF systems in some cases.展开更多
The purpose of this work is to implement a discontinuous Galerkin(DG)method with a one-sided flux for a singularly perturbed Volterra integro-differential equation(VIDE)with a smooth kernel.First,the regularity proper...The purpose of this work is to implement a discontinuous Galerkin(DG)method with a one-sided flux for a singularly perturbed Volterra integro-differential equation(VIDE)with a smooth kernel.First,the regularity property and a decomposition of the exact solution of the singularly perturbed VIDE with the initial condition are provided.Then the existence and uniqueness of the DG solution are proven.Then some appropriate projection-type interpolation operators and their corresponding approximation properties are established.Based on the decomposition of the exact solution and the approximation properties of the projection type interpolants,the DG method achieves the uniform convergence in the L2 norm with respect to the singular perturbation parameter e when the space of polynomials with degree p is used.A numerical experiment validates the theoretical results.Furthermore,an ultra-convergence order 2p+1 at the nodes for the one-sided flux,uniform with respect to the singular perturbation parameter e,is observed numerically.展开更多
In this paper, using the interpolation perturbation method. the author seeks tosolve several nonlinear problems. Numerical examples show that the method Df thispaper has good accuracy.
We present a perturbation study of the ground-state energy of the beryllium atom by incorporating double parameters in the atom's Hamiltonian. The eigenvalue of the Hamiltonian is then solved with a double-fold pertu...We present a perturbation study of the ground-state energy of the beryllium atom by incorporating double parameters in the atom's Hamiltonian. The eigenvalue of the Hamiltonian is then solved with a double-fold perturbation scheme,where the spin-spin interaction of electrons from different shells of the atom is also considered. Calculations show that the obtained ground-state energy is in satisfactory agreement with experiment. It is found that the Coulomb repulsion of the inner-shell electrons enhances the effective nuclear charge seen by the outer-shell electrons, and the shielding effect of the outer-shell electrons to the nucleus is also notable compared with that of the inner-shell electrons.展开更多
The Homotopy Perturbation Method (HPM) is used to solve the Burgers-Huxley non-linear differential equations. Three case study problems of Burgers-Huxley are solved using the HPM and the exact solutions are obtained. ...The Homotopy Perturbation Method (HPM) is used to solve the Burgers-Huxley non-linear differential equations. Three case study problems of Burgers-Huxley are solved using the HPM and the exact solutions are obtained. The rapid convergence towards the exact solutions of HPM is numerically shown. Results show that the HPM is efficient method with acceptable accuracy to solve the Burgers-Huxley equation. Also, the results prove that the method is an efficient and powerful algorithm to construct the exact solution of non-linear differential equations.展开更多
This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matr...This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.展开更多
The finite_element_displacement_perturbation method (FEDPM)for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes (Ⅰ) was e...The finite_element_displacement_perturbation method (FEDPM)for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes (Ⅰ) was employed to calculate the stress distributions and the stiffness of the bellows. Firstly, by applying the first_order perturbation solution (the linear solution)of the FEDPM to the bellows, the obtained results were compared with those of the general solution and the initial parameter integration solution proposed by the present authors earlier, as well as of the experiments and the FEA by others.It is shown that the FEDPM is with good precision and reliability, and as it was pointed out in (Ⅰ) the abrupt changes of the meridian curvature of bellows would not affect the use of the usual straight element. Then the nonlinear behaviors of the bellows were discussed. As expected, the nonlinear effects mainly come from the bellows ring plate,and the wider the ring plate is, the stronger the nonlinear effects are. Contrarily, the vanishing of the ring plate, like the C_shaped bellows, the nonlinear effects almost vanish. In addition, when the pure bending moments act on the bellows, each convolution has the same stress distributions calculated by the linear solution and other linear theories, but by the present nonlinear solution they vary with respect to the convolutions of the bellows. Yet for most bellows, the linear solutions are valid in practice.展开更多
In this article, the application of variational homotopy perturbation method is applied to solve Benjamin-Bona-Mahony equation. Then, we obtain the numerical solution of BBM equation using the initial condition. Compa...In this article, the application of variational homotopy perturbation method is applied to solve Benjamin-Bona-Mahony equation. Then, we obtain the numerical solution of BBM equation using the initial condition. Comparison with Adomian’s decomposition method, homotopy perturbation method, and with the exact solution shows that VHPM is more effective and accurate than ADM and HPM, and is reliable and manageable for this type of equation.展开更多
Based on the first-order upwind and second-order central type of finite volume (UFV and CFV) scheme, upwind and central type of perturbation finite volume (UPFV and CPFV) schemes of the Navier-Stokes equations were de...Based on the first-order upwind and second-order central type of finite volume (UFV and CFV) scheme, upwind and central type of perturbation finite volume (UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from (0.3) to (0.8) and convergence perform excellent with Reynolds number variation from 10~2 to 10~4.展开更多
In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomi...In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomials [1]. The analytical results of examples are calculated in terms of convergent series with easily computed components [2].展开更多
The main aim of the paper is to examine the concentration of the longitudinal dispersion phenomenon arising in fluid flow through porous media. These phenomenon yields a partial differential equation namely Burger’s ...The main aim of the paper is to examine the concentration of the longitudinal dispersion phenomenon arising in fluid flow through porous media. These phenomenon yields a partial differential equation namely Burger’s equation, which is solved by mixture of the new integral transform and the homotopy perturbation method under suitable conditions and the standard assumption. This method provides an analytical approximation in a rapidly convergent sequence with in exclusive manner computed terms. Its rapid convergence shows that the method is trustworthy and introduces a significant improvement in solving nonlinear partial differential equations over existing methods. It is concluded that the behaviour of concentration in longitudinal dispersion phenomenon is decreases as distance x is increasing with fixed time t > 0 and slightly increases with time t.展开更多
This research paper deals with the boundary and initial value problems for the Bratu-type model by using the New Improved Variational Homotopy Perturbation Method. The New Method does not require discritization, linea...This research paper deals with the boundary and initial value problems for the Bratu-type model by using the New Improved Variational Homotopy Perturbation Method. The New Method does not require discritization, linearization or any restrictive assumption of any form in providing analytical or approximate solutions to linear and nonlinear equation without the integral related with nonlinear term. Theses virtues make it to be reliable and its efficiency is demonstrated with numerical examples.展开更多
In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We the...In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.展开更多
Recently application of functionally graded materials(FGMs) have attracted a great deal of interest. These materials are composed of various materials with different micro-structures which can vary spatially in FGMs...Recently application of functionally graded materials(FGMs) have attracted a great deal of interest. These materials are composed of various materials with different micro-structures which can vary spatially in FGMs. Such composites with varying thickness and non-uniform pressure can be used in the aerospace engineering. Therefore, analysis of such composite is of high importance in engineering problems. Thermoelastic analysis of functionally graded cylinder with variable thickness under non-uniform pressure is considered. First order shear deformation theory and total potential energy approach is applied to obtain the governing equations of non-homogeneous cylinder. Considering the inner and outer solutions, perturbation series are applied to solve the governing equations. Outer solution for out of boundaries and more sensitive variable in inner solution at the boundaries are considered. Combining of inner and outer solution for near and far points from boundaries leads to high accurate displacement field distribution. The main aim of this paper is to show the capability of matched asymptotic solution for different non-homogeneous cylinders with different shapes and different non-uniform pressures. The results can be used to design the optimum thickness of the cylinder and also some properties such as high temperature residence by applying non-homogeneous material.展开更多
This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order...This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order pertur- bation equation, which is solved approximately by resolv- ing the Hamiltonian coefficient matrix into a "major compo- nent" and a "high order small quantity" and using perturba- tion transformation technique, then the solution to the orig- inal equation of Hamiltonian system is determined through a series of inverse transform. Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes, the transfer matrix is a sym- plectic matrix; furthermore, the exponential matrices can be calculated accurately by the precise time integration method, so the method presented in this paper has fine accuracy, ef- ficiency and stability. The examples show that the proposed method can also give good results even though a large time step is selected, and with the increase of the perturbation or- der, the perturbation solutions tend to exact solutions rapidly.展开更多
基金sponsored by the Graduate Student Research and Innovation Fund of Xinyang Normal University under No.2024KYJJ012.
文摘In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.
基金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.
基金supported by the Natural Science Foundation of Zhejiang Province,China(Grant Nos.LY20A010021,LY19A010002,LY20G030025)the Natural Science Founda-tion of Ningbo City,China(Grant Nos.2021J147,2021J235).
文摘By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11071205 and 11101349), the “Strate- gic Priority Research Program-Climate Change: Carbon Budget and Relevant Issues” of the Chinese Academy of Sciences, China (Grant No. XDA01020304), the Natural Science Foundation from the Education Bureau of Anhui Province, China (Grant No. KJ2011A135), and the Natural Science Foundation of Jiangsu Province, China (Grant No. BK2011042).
文摘The EI Nino-southern oscillation (ENSO) is an interannual phenomenon involved in tropical Pacific ocean- atmosphere interactions. In this paper, we develop an asymptotic method of solving the nonlinear equation using the ENSO model. Based on a class of the oscillator of the ENSO model, a approximate solution of the corresponding problem is studied employing the perturbation method.
基金Supported by National Natural Science Foundation of China(11071198)
文摘In this article, we deal with a class of semilinear elliptic equations which are perturbations of the problems with the critical Hardy-Sobolev exponent. Some existence results are given via an abstract perturbation method in critical point theory.
基金National Natural Science Foundation under Grant No. 51179093National Basic Research Program of China under Grant No. 2011CB013602Program for New Century Excellent Talents in University under Grant No.NCET-10-0531
文摘A third-order correction was recently suggested to improve the accuracy of the half-power bandwidth method in estimating the damping of single DOF systems.This paper analyzes the accuracy of the half-power bandwidth method with the third-order correction in damping estimation for multi-DOF linear systems.Damping ratios in a two-DOF linear system are estimated using its displacement and acceleration frequency response curves,respectively.A wide range of important parameters that characterize the shape of these response curves are taken into account.Results show that the third-order correction may greatly improve the accuracy of the half-power bandwidth method in estimating damping in a two-DOF system.In spite of this,the half-power bandwidth method may significantly overestimate the damping ratios of two-DOF systems in some cases.
基金supported by the National Natural Science Foundation of China(12001189)supported by the National Natural Science Foundation of China(11171104,12171148)。
文摘The purpose of this work is to implement a discontinuous Galerkin(DG)method with a one-sided flux for a singularly perturbed Volterra integro-differential equation(VIDE)with a smooth kernel.First,the regularity property and a decomposition of the exact solution of the singularly perturbed VIDE with the initial condition are provided.Then the existence and uniqueness of the DG solution are proven.Then some appropriate projection-type interpolation operators and their corresponding approximation properties are established.Based on the decomposition of the exact solution and the approximation properties of the projection type interpolants,the DG method achieves the uniform convergence in the L2 norm with respect to the singular perturbation parameter e when the space of polynomials with degree p is used.A numerical experiment validates the theoretical results.Furthermore,an ultra-convergence order 2p+1 at the nodes for the one-sided flux,uniform with respect to the singular perturbation parameter e,is observed numerically.
文摘In this paper, using the interpolation perturbation method. the author seeks tosolve several nonlinear problems. Numerical examples show that the method Df thispaper has good accuracy.
基金Project supported by the National Natural Science Foundation of China(Grant No.11647071)the Natural Science Foundation of Jiangsu Province,China(Grant No.BK20160435)
文摘We present a perturbation study of the ground-state energy of the beryllium atom by incorporating double parameters in the atom's Hamiltonian. The eigenvalue of the Hamiltonian is then solved with a double-fold perturbation scheme,where the spin-spin interaction of electrons from different shells of the atom is also considered. Calculations show that the obtained ground-state energy is in satisfactory agreement with experiment. It is found that the Coulomb repulsion of the inner-shell electrons enhances the effective nuclear charge seen by the outer-shell electrons, and the shielding effect of the outer-shell electrons to the nucleus is also notable compared with that of the inner-shell electrons.
文摘The Homotopy Perturbation Method (HPM) is used to solve the Burgers-Huxley non-linear differential equations. Three case study problems of Burgers-Huxley are solved using the HPM and the exact solutions are obtained. The rapid convergence towards the exact solutions of HPM is numerically shown. Results show that the HPM is efficient method with acceptable accuracy to solve the Burgers-Huxley equation. Also, the results prove that the method is an efficient and powerful algorithm to construct the exact solution of non-linear differential equations.
基金Project supported by the National Natural Science Foundation of China(No.10672194)the China-Russia Cooperative Project(the National Natural Science Foundation of China and the Russian Foundation for Basic Research)(No.10811120012)
文摘This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.
文摘The finite_element_displacement_perturbation method (FEDPM)for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes (Ⅰ) was employed to calculate the stress distributions and the stiffness of the bellows. Firstly, by applying the first_order perturbation solution (the linear solution)of the FEDPM to the bellows, the obtained results were compared with those of the general solution and the initial parameter integration solution proposed by the present authors earlier, as well as of the experiments and the FEA by others.It is shown that the FEDPM is with good precision and reliability, and as it was pointed out in (Ⅰ) the abrupt changes of the meridian curvature of bellows would not affect the use of the usual straight element. Then the nonlinear behaviors of the bellows were discussed. As expected, the nonlinear effects mainly come from the bellows ring plate,and the wider the ring plate is, the stronger the nonlinear effects are. Contrarily, the vanishing of the ring plate, like the C_shaped bellows, the nonlinear effects almost vanish. In addition, when the pure bending moments act on the bellows, each convolution has the same stress distributions calculated by the linear solution and other linear theories, but by the present nonlinear solution they vary with respect to the convolutions of the bellows. Yet for most bellows, the linear solutions are valid in practice.
文摘In this article, the application of variational homotopy perturbation method is applied to solve Benjamin-Bona-Mahony equation. Then, we obtain the numerical solution of BBM equation using the initial condition. Comparison with Adomian’s decomposition method, homotopy perturbation method, and with the exact solution shows that VHPM is more effective and accurate than ADM and HPM, and is reliable and manageable for this type of equation.
文摘Based on the first-order upwind and second-order central type of finite volume (UFV and CFV) scheme, upwind and central type of perturbation finite volume (UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from (0.3) to (0.8) and convergence perform excellent with Reynolds number variation from 10~2 to 10~4.
文摘In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomials [1]. The analytical results of examples are calculated in terms of convergent series with easily computed components [2].
文摘The main aim of the paper is to examine the concentration of the longitudinal dispersion phenomenon arising in fluid flow through porous media. These phenomenon yields a partial differential equation namely Burger’s equation, which is solved by mixture of the new integral transform and the homotopy perturbation method under suitable conditions and the standard assumption. This method provides an analytical approximation in a rapidly convergent sequence with in exclusive manner computed terms. Its rapid convergence shows that the method is trustworthy and introduces a significant improvement in solving nonlinear partial differential equations over existing methods. It is concluded that the behaviour of concentration in longitudinal dispersion phenomenon is decreases as distance x is increasing with fixed time t > 0 and slightly increases with time t.
文摘This research paper deals with the boundary and initial value problems for the Bratu-type model by using the New Improved Variational Homotopy Perturbation Method. The New Method does not require discritization, linearization or any restrictive assumption of any form in providing analytical or approximate solutions to linear and nonlinear equation without the integral related with nonlinear term. Theses virtues make it to be reliable and its efficiency is demonstrated with numerical examples.
基金The project supported by the China NKBRSF(2001CB409604)
文摘In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.
文摘Recently application of functionally graded materials(FGMs) have attracted a great deal of interest. These materials are composed of various materials with different micro-structures which can vary spatially in FGMs. Such composites with varying thickness and non-uniform pressure can be used in the aerospace engineering. Therefore, analysis of such composite is of high importance in engineering problems. Thermoelastic analysis of functionally graded cylinder with variable thickness under non-uniform pressure is considered. First order shear deformation theory and total potential energy approach is applied to obtain the governing equations of non-homogeneous cylinder. Considering the inner and outer solutions, perturbation series are applied to solve the governing equations. Outer solution for out of boundaries and more sensitive variable in inner solution at the boundaries are considered. Combining of inner and outer solution for near and far points from boundaries leads to high accurate displacement field distribution. The main aim of this paper is to show the capability of matched asymptotic solution for different non-homogeneous cylinders with different shapes and different non-uniform pressures. The results can be used to design the optimum thickness of the cylinder and also some properties such as high temperature residence by applying non-homogeneous material.
基金supported by the National Natural Science Foun-dation of China (11172334)
文摘This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order pertur- bation equation, which is solved approximately by resolv- ing the Hamiltonian coefficient matrix into a "major compo- nent" and a "high order small quantity" and using perturba- tion transformation technique, then the solution to the orig- inal equation of Hamiltonian system is determined through a series of inverse transform. Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes, the transfer matrix is a sym- plectic matrix; furthermore, the exponential matrices can be calculated accurately by the precise time integration method, so the method presented in this paper has fine accuracy, ef- ficiency and stability. The examples show that the proposed method can also give good results even though a large time step is selected, and with the increase of the perturbation or- der, the perturbation solutions tend to exact solutions rapidly.