A higher order boundary element method(HOBEM)is presented for inviscid flow passing cylinders in bounded or unbounded domain.The traditional boundary integral equation is established with respect to the velocity poten...A higher order boundary element method(HOBEM)is presented for inviscid flow passing cylinders in bounded or unbounded domain.The traditional boundary integral equation is established with respect to the velocity potential and its normal derivative.In present work,a new integral equation is derived for the tangential velocity.The boundary is discretized into higher order elements to ensure the continuity of slope at the element nodes.The velocity potential is also expanded with higher order shape functions,in which the unknown coefficients involve the tangential velocity.The expansion then ensures the continuities of the velocity and the slope of the boundary at element nodes.Through extensive comparison of the results for the analytical solution of cylinders,it is shown that the present HOBEM is much more accurate than the conventional BEM.展开更多
The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quad...The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quadratic Newton-Raphson method.These six methods do not require higher order derivatives to achieve a higher convergence rate.Six algorithms are developed to use the higher order methods in place of the Newton-Raphson method to solve the nonlinear equilibrium equations in geometrically nonlinear analysis of structures.The higher order methods are applied to both continuum and discrete problems(spherical shell and dome truss).The computational cost and the sensitivity of the higher order solution methods and the Newton-Raphson method with respect to the load increment size are comparatively investigated.The numerical results reveal that the higher order methods require a lower number of iterations that the Newton-Raphson method to converge.It is also shown that these methods are less sensitive to the variation of the load increment size.As it is indicated in numerical results,the average residual reduces in a lower number of iterations by the application of the higher order methods in the nonlinear analysis of structures.展开更多
The higher order displacement discontinuity method(HODDM) utilizing special crack tip elements has been used in the solution of linear elastic fracture mechanics(LEFM) problems. The paper has selected several example ...The higher order displacement discontinuity method(HODDM) utilizing special crack tip elements has been used in the solution of linear elastic fracture mechanics(LEFM) problems. The paper has selected several example problems from the fracture mechanics literature(with available analytical solutions) including center slant crack in an infinite and finite body, single and double edge cracks, cracks emanating from a circular hole. The numerical values of Mode Ⅰ and Mode Ⅱ SIFs for these problems using HODDM are in excellent agreement with analytical results(reaching up to 0.001% deviation from their analytical results). The HODDM is also compared with the XFEM and a modified XFEM results. The results show that the HODDM needs a considerably lower computational effort(with less than 400 nodes) than the XFEM and the modified XFEM(which needs more than 10000 nodes) to reach a much higher accuracy. The proposed HODDM offers higher accuracy and lower computation effort for a wide range of problems in LEFM.展开更多
In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs...In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.展开更多
One-dimensional heat equation was solved for different higher-order finite difference schemes, namely, forward time and fourth-order centered space explicit method, backward time and fourth-order centered space implic...One-dimensional heat equation was solved for different higher-order finite difference schemes, namely, forward time and fourth-order centered space explicit method, backward time and fourth-order centered space implicit method, and fourth-order implicit Crank-Nicolson finite difference method. Higher-order schemes have complexity in computing values at the neighboring points to the boundaries. It is required there a specification of the values of field variables at some points exterior to the domain. The complexity was incorporated using Hicks approximation. The convergence and stability analysis was also computed for those higher-order finite difference explicit and implicit methods in case of solving a one dimensional heat equation. The obtained numerical results were compared with exact solutions. It is found that backward time and fourth-order centered space implicit scheme along with Hicks approximation performed well over the other mentioned higher-order approaches.展开更多
This paper is devoted to find the numerical solutions of one dimensional general nonlinear system of third-order boundary value problems (BVPs) for the pair of functions using Galerkin weighted residual method. We der...This paper is devoted to find the numerical solutions of one dimensional general nonlinear system of third-order boundary value problems (BVPs) for the pair of functions using Galerkin weighted residual method. We derive mathematical formulations in matrix form, in detail, by exploiting Bernstein polynomials as basis functions. A reasonable accuracy is found when the proposed method is used on few examples. At the end of the study, a comparison is made between the approximate and exact solutions, and also with the solutions of the existing methods. Our results converge monotonically to the exact solutions. In addition, we show that the derived formulations may be applicable by reducing higher order complicated BVP into a lower order system of BVPs, and the performance of the numerical solutions is satisfactory. .展开更多
Considered under their standard form, the fifth-order KdV equations are a sort of reading table on which new prototypes of higher order solitary waves residing there, have been uncovered and revealed to broad daylight...Considered under their standard form, the fifth-order KdV equations are a sort of reading table on which new prototypes of higher order solitary waves residing there, have been uncovered and revealed to broad daylight. The mathematical tool that made it possible to explore and analyze this equation is the Bogning-Djeumen Tchaho-Kofané method extended to the new implicit Bogning' functions. The analytical form of the solutions chosen in this manuscript is particular in the sense that it contains within its bosom, a package of solitary waves made up of three solitons, especially, the bright type soliton, the hybrid soliton and the dark type soliton which we estimate capable in their interactions of generating new hybrid or multi-form solitons. Existence conditions of the obtained solitons have been determined. It emerges that, these existence conditions of the chosen ansatz could open the way to other new varieties of fifth-order KdV equations including to which it will be one of the solutions. Some of the obtained solitons are exact solutions. Intense numerical simulations highlighted numerical stability and confirmed the hybrid character of the obtained solutions. These results will help to model new nonlinear wave phenomena, in plasma media and in fluid dynamics, especially, on the shallow water surface.展开更多
An approximate analysis for free vibration of a laminated curved panel(shell)with four edges simply supported(SS2),is presented in this paper.The transverse shear deformation is considered by using a higher-order shea...An approximate analysis for free vibration of a laminated curved panel(shell)with four edges simply supported(SS2),is presented in this paper.The transverse shear deformation is considered by using a higher-order shear deformation theory.For solving the highly coupled partial differential governing equations and associated boundary conditions,a set of solution functions in the form of double trigonometric Fourier series,which are required to satisfy the geometry part of the considered boundary conditions,is assumed in advance.By applying the Galerkin procedure both to the governing equations and to the natural boundary conditions not satisfied by the assumed solution functions,an approximate solution,capable of providing a reliable prediction for the global response of the panel,is obtained.Numerical results of antisymmetric angle-ply as well as symmetric cross-ply and angle-ply laminated curved panels are presented and discussed.展开更多
The boundary-layer method is used to study a wide moving jam to a class of higher-order viscous models. The equations for characteristic parameters are derived to determine the asymptotic solution. The sufficient and ...The boundary-layer method is used to study a wide moving jam to a class of higher-order viscous models. The equations for characteristic parameters are derived to determine the asymptotic solution. The sufficient and essential conditions for the wide moving jam formation are discussed in detail, respectively, and then used to prove or disprove the existence of the wide moving jam solutions to many well-known higher-order models. It is shown that the numerical results agree with the analytical results.展开更多
By using the extended homogeneous balance method, the localized coherent structures are studied. A nonlinear transformation was first established, and then the linearization form was obtained based on the extended hom...By using the extended homogeneous balance method, the localized coherent structures are studied. A nonlinear transformation was first established, and then the linearization form was obtained based on the extended homogeneous balance method for the higher order (2 + 1)-dimensional Broer-Kaup equations. Starting from this linearization form equation, a variable separation solution with the entrance of some arbitrary functions and some arbitrary parameters was constructed. The quite rich localized coherent structures were revealed. This method, which can be generalized to other (2 + I) -dimensional nonlinear evolution equation, is simple and powerful.展开更多
The initial boundary value problems for a class of high order Kirchhoff type equations with nonlinear strongly damped terms are considered. We establish the existence and uniqueness of the global solution of the probl...The initial boundary value problems for a class of high order Kirchhoff type equations with nonlinear strongly damped terms are considered. We establish the existence and uniqueness of the global solution of the problem by using prior estimates and Galerkin’s method under proper assumptions for the rigid term. Then the compact method is used to prove the existence of a compact family of global attractors in the solution semigroup generated by the problem. Finally, the Frechet differentiability of the operator semigroup and the decay of the volume element of linearization problem are proved, and the Hausdorff dimension and Fractal dimension of the family of global attractors are obtained.展开更多
A differential quadrature (DQ) method for orthotropic plates was proposed based on Reddy' s theory of plates with the effects of the higher-order transverse shear deformations. Wang-Bert's DQ approach was also...A differential quadrature (DQ) method for orthotropic plates was proposed based on Reddy' s theory of plates with the effects of the higher-order transverse shear deformations. Wang-Bert's DQ approach was also further extended to handle the boundary conditions of plates. The computational convergence was studied, and the numerical results were obtained for different grid spacings and compared with the existing results. The results show that the DQ method is fairly reliable and effective.展开更多
Through a higher-order boundary element method based on NURBS （Non-uniform Rational B-splines）, the calculation of second-order low-frequency forces and slow drift motions is conducted for floating bodies. In the fl...Through a higher-order boundary element method based on NURBS （Non-uniform Rational B-splines）, the calculation of second-order low-frequency forces and slow drift motions is conducted for floating bodies. In the floating body＇s inner domain, an auxiliary equation is obtained by applying a Green function which satisfies the solid surface condition. Then, the auxiliary equation and the velocity potential equation are combined in the fluid domain to remove the solid angle coefficient and the singularity of the double layer potentials in the integral equation. Thus, a new velocity potential integral equation is obtained. The new equation is extended to the inner domain to reheve the irregular frequency effects; on the basis of the order analysis, the comparison is made about the contribution of all integral terms with the result in the second-order tow-frequency problem; the higher-order boundary element method based on NURBS is apphed to calculate the geometric position and velocity potentials; the slow drift motions are calculated by the spectrum analysis method. Removing the solid angle coefficient can apply NURBS technology to the hydrodynamic calculation of floating bodies with complex surfaces, and the extended boundary integral method can reduce the irregular frequency effects. Order analysis shows that free surface integral can be neglected, and the numerical results can also prove the correctness of order analysis. The results of second-order low-frequency forces and slow drift motions and the comparison with the results from references show that the application of the NURBS technology to the second-order low-frequency problem is of high efficiency and credible results.展开更多
This paper studies integration of a higher-order differential equation which can be reduced to a second-order ordinary differential equation. The solution of the second-order equation can be obtained by the Noether me...This paper studies integration of a higher-order differential equation which can be reduced to a second-order ordinary differential equation. The solution of the second-order equation can be obtained by the Noether method and the Poisson method. Then the solution of the higher-order equation can be obtained by integrating the solution of the second-order equation.展开更多
基金financially supported by the National Natural Science Foundation of China (Grant Nos.52271276,52271319,and 52201364)the Natural Science Foundation of Jiangsu Province (Grant No.BK20201006)。
文摘A higher order boundary element method(HOBEM)is presented for inviscid flow passing cylinders in bounded or unbounded domain.The traditional boundary integral equation is established with respect to the velocity potential and its normal derivative.In present work,a new integral equation is derived for the tangential velocity.The boundary is discretized into higher order elements to ensure the continuity of slope at the element nodes.The velocity potential is also expanded with higher order shape functions,in which the unknown coefficients involve the tangential velocity.The expansion then ensures the continuities of the velocity and the slope of the boundary at element nodes.Through extensive comparison of the results for the analytical solution of cylinders,it is shown that the present HOBEM is much more accurate than the conventional BEM.
文摘The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quadratic Newton-Raphson method.These six methods do not require higher order derivatives to achieve a higher convergence rate.Six algorithms are developed to use the higher order methods in place of the Newton-Raphson method to solve the nonlinear equilibrium equations in geometrically nonlinear analysis of structures.The higher order methods are applied to both continuum and discrete problems(spherical shell and dome truss).The computational cost and the sensitivity of the higher order solution methods and the Newton-Raphson method with respect to the load increment size are comparatively investigated.The numerical results reveal that the higher order methods require a lower number of iterations that the Newton-Raphson method to converge.It is also shown that these methods are less sensitive to the variation of the load increment size.As it is indicated in numerical results,the average residual reduces in a lower number of iterations by the application of the higher order methods in the nonlinear analysis of structures.
文摘The higher order displacement discontinuity method(HODDM) utilizing special crack tip elements has been used in the solution of linear elastic fracture mechanics(LEFM) problems. The paper has selected several example problems from the fracture mechanics literature(with available analytical solutions) including center slant crack in an infinite and finite body, single and double edge cracks, cracks emanating from a circular hole. The numerical values of Mode Ⅰ and Mode Ⅱ SIFs for these problems using HODDM are in excellent agreement with analytical results(reaching up to 0.001% deviation from their analytical results). The HODDM is also compared with the XFEM and a modified XFEM results. The results show that the HODDM needs a considerably lower computational effort(with less than 400 nodes) than the XFEM and the modified XFEM(which needs more than 10000 nodes) to reach a much higher accuracy. The proposed HODDM offers higher accuracy and lower computation effort for a wide range of problems in LEFM.
文摘In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.
文摘One-dimensional heat equation was solved for different higher-order finite difference schemes, namely, forward time and fourth-order centered space explicit method, backward time and fourth-order centered space implicit method, and fourth-order implicit Crank-Nicolson finite difference method. Higher-order schemes have complexity in computing values at the neighboring points to the boundaries. It is required there a specification of the values of field variables at some points exterior to the domain. The complexity was incorporated using Hicks approximation. The convergence and stability analysis was also computed for those higher-order finite difference explicit and implicit methods in case of solving a one dimensional heat equation. The obtained numerical results were compared with exact solutions. It is found that backward time and fourth-order centered space implicit scheme along with Hicks approximation performed well over the other mentioned higher-order approaches.
文摘This paper is devoted to find the numerical solutions of one dimensional general nonlinear system of third-order boundary value problems (BVPs) for the pair of functions using Galerkin weighted residual method. We derive mathematical formulations in matrix form, in detail, by exploiting Bernstein polynomials as basis functions. A reasonable accuracy is found when the proposed method is used on few examples. At the end of the study, a comparison is made between the approximate and exact solutions, and also with the solutions of the existing methods. Our results converge monotonically to the exact solutions. In addition, we show that the derived formulations may be applicable by reducing higher order complicated BVP into a lower order system of BVPs, and the performance of the numerical solutions is satisfactory. .
文摘Considered under their standard form, the fifth-order KdV equations are a sort of reading table on which new prototypes of higher order solitary waves residing there, have been uncovered and revealed to broad daylight. The mathematical tool that made it possible to explore and analyze this equation is the Bogning-Djeumen Tchaho-Kofané method extended to the new implicit Bogning' functions. The analytical form of the solutions chosen in this manuscript is particular in the sense that it contains within its bosom, a package of solitary waves made up of three solitons, especially, the bright type soliton, the hybrid soliton and the dark type soliton which we estimate capable in their interactions of generating new hybrid or multi-form solitons. Existence conditions of the obtained solitons have been determined. It emerges that, these existence conditions of the chosen ansatz could open the way to other new varieties of fifth-order KdV equations including to which it will be one of the solutions. Some of the obtained solitons are exact solutions. Intense numerical simulations highlighted numerical stability and confirmed the hybrid character of the obtained solutions. These results will help to model new nonlinear wave phenomena, in plasma media and in fluid dynamics, especially, on the shallow water surface.
文摘An approximate analysis for free vibration of a laminated curved panel(shell)with four edges simply supported(SS2),is presented in this paper.The transverse shear deformation is considered by using a higher-order shear deformation theory.For solving the highly coupled partial differential governing equations and associated boundary conditions,a set of solution functions in the form of double trigonometric Fourier series,which are required to satisfy the geometry part of the considered boundary conditions,is assumed in advance.By applying the Galerkin procedure both to the governing equations and to the natural boundary conditions not satisfied by the assumed solution functions,an approximate solution,capable of providing a reliable prediction for the global response of the panel,is obtained.Numerical results of antisymmetric angle-ply as well as symmetric cross-ply and angle-ply laminated curved panels are presented and discussed.
基金Project supported by the National Natural Science Foundation of China(No.11602128)the Natural Science Foundation of Fujian Province of China(No.2016J01679)
文摘The boundary-layer method is used to study a wide moving jam to a class of higher-order viscous models. The equations for characteristic parameters are derived to determine the asymptotic solution. The sufficient and essential conditions for the wide moving jam formation are discussed in detail, respectively, and then used to prove or disprove the existence of the wide moving jam solutions to many well-known higher-order models. It is shown that the numerical results agree with the analytical results.
文摘By using the extended homogeneous balance method, the localized coherent structures are studied. A nonlinear transformation was first established, and then the linearization form was obtained based on the extended homogeneous balance method for the higher order (2 + 1)-dimensional Broer-Kaup equations. Starting from this linearization form equation, a variable separation solution with the entrance of some arbitrary functions and some arbitrary parameters was constructed. The quite rich localized coherent structures were revealed. This method, which can be generalized to other (2 + I) -dimensional nonlinear evolution equation, is simple and powerful.
文摘The initial boundary value problems for a class of high order Kirchhoff type equations with nonlinear strongly damped terms are considered. We establish the existence and uniqueness of the global solution of the problem by using prior estimates and Galerkin’s method under proper assumptions for the rigid term. Then the compact method is used to prove the existence of a compact family of global attractors in the solution semigroup generated by the problem. Finally, the Frechet differentiability of the operator semigroup and the decay of the volume element of linearization problem are proved, and the Hausdorff dimension and Fractal dimension of the family of global attractors are obtained.
基金key Project of the Municipal Commission of Science and Technology of Shanghai
文摘A differential quadrature (DQ) method for orthotropic plates was proposed based on Reddy' s theory of plates with the effects of the higher-order transverse shear deformations. Wang-Bert's DQ approach was also further extended to handle the boundary conditions of plates. The computational convergence was studied, and the numerical results were obtained for different grid spacings and compared with the existing results. The results show that the DQ method is fairly reliable and effective.
文摘Through a higher-order boundary element method based on NURBS （Non-uniform Rational B-splines）, the calculation of second-order low-frequency forces and slow drift motions is conducted for floating bodies. In the floating body＇s inner domain, an auxiliary equation is obtained by applying a Green function which satisfies the solid surface condition. Then, the auxiliary equation and the velocity potential equation are combined in the fluid domain to remove the solid angle coefficient and the singularity of the double layer potentials in the integral equation. Thus, a new velocity potential integral equation is obtained. The new equation is extended to the inner domain to reheve the irregular frequency effects; on the basis of the order analysis, the comparison is made about the contribution of all integral terms with the result in the second-order tow-frequency problem; the higher-order boundary element method based on NURBS is apphed to calculate the geometric position and velocity potentials; the slow drift motions are calculated by the spectrum analysis method. Removing the solid angle coefficient can apply NURBS technology to the hydrodynamic calculation of floating bodies with complex surfaces, and the extended boundary integral method can reduce the irregular frequency effects. Order analysis shows that free surface integral can be neglected, and the numerical results can also prove the correctness of order analysis. The results of second-order low-frequency forces and slow drift motions and the comparison with the results from references show that the application of the NURBS technology to the second-order low-frequency problem is of high efficiency and credible results.
基金Project supported by the National Natural Science Foundation of China(Grant No10572021)Doctoral Programme Foundation of Institution of Higher Education of China(Grant No20040007022)
文摘This paper studies integration of a higher-order differential equation which can be reduced to a second-order ordinary differential equation. The solution of the second-order equation can be obtained by the Noether method and the Poisson method. Then the solution of the higher-order equation can be obtained by integrating the solution of the second-order equation.
