This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredepend...This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.展开更多
This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node...This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.展开更多
Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridyna...Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.展开更多
To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’...To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.展开更多
The numerical approach for finding the solution of fractional order systems of boundary value problems (BPVs) is derived in this paper. The implementation of the weighted residuals such as Galerkin, Least Square, and ...The numerical approach for finding the solution of fractional order systems of boundary value problems (BPVs) is derived in this paper. The implementation of the weighted residuals such as Galerkin, Least Square, and Collocation methods are included for solving fractional order differential equations, which is broadened to acquire the approximate solutions of fractional order systems with differentiable polynomials, namely Legendre polynomials, as basis functions. The algorithm of the residual formulations of matrix form can be coded efficiently. The interpretation of Caputo fractional derivatives is employed here. We have demonstrated these methods numerically through a few examples of linear and nonlinear BVPs. The results in absolute errors show that the present method efficiently finds the numerical solutions of fractional order systems of differential equations.展开更多
In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled b...In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.展开更多
The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect ...The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect is considered. Detailed formulations and procedures are given for the first time. Several examples are analyzed and accurate results are obtained by the DQEM. Based on the results in this paper, one may conclude that the DQEM is a useful tool for obtaining solutions of structural elements. It can be seen that the shell shape may be theore tically controlled and snap through may occur when the applied voltage reaches a critical value even without mechanical load for certain geometric configurations.展开更多
This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeabi...This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.展开更多
It is a new attempt to extend the differential quadrature method(DQM) to stability analysis of the straight and curved centerlinepipes conveying fluid. Emphasis is placed on the study of theinfluences of several param...It is a new attempt to extend the differential quadrature method(DQM) to stability analysis of the straight and curved centerlinepipes conveying fluid. Emphasis is placed on the study of theinfluences of several parameters on the critical flow velocity.Compared to other methods, this method can more easily deal with thepipe with spring support at its boundaries and asks for much lesscomputing effort while giving ac- ceptable precision in the numericalresults.展开更多
The Hamilton-Jacobi method for solving ordinary differential equations is presented in this paper. A system of ordinary differential equations of first order or second order can be expressed as a Hamilton system under...The Hamilton-Jacobi method for solving ordinary differential equations is presented in this paper. A system of ordinary differential equations of first order or second order can be expressed as a Hamilton system under certain conditions. Then the Hamilton-Jacobi method is used in the integration of the Hamilton system and the solution of the original ordinary differential equations can be found. Finally, an example is given to illustrate the application of the result.展开更多
A differential equation of first order can be expressed by the equation of motion of a mechanical system. In this paper, three methods of analytical mechanics, i.e. the Hamilton-Noether method, the Lagrange-Noether me...A differential equation of first order can be expressed by the equation of motion of a mechanical system. In this paper, three methods of analytical mechanics, i.e. the Hamilton-Noether method, the Lagrange-Noether method and the Poisson method, are given to solve a differential equation of first order, of which the way may be called the mechanical methodology in mathematics.展开更多
The quasi-static and dynamic responses of a thermoviscoelastic Timoshenko beam subject to thermal loads are analyzed. First, based on the small geometric deformation assumption and Boltzmann constitutive relation, the...The quasi-static and dynamic responses of a thermoviscoelastic Timoshenko beam subject to thermal loads are analyzed. First, based on the small geometric deformation assumption and Boltzmann constitutive relation, the governing equations for the beam are presented. Second, an extended differential quadrature method(DQM)in the spatial domain and a differential method in the temporal domain are combined to transform the integro-partial-differential governing equations into the ordinary differential equations. Third, the accuracy of the present discrete method is verified by elastic/viscoelastic examples, and the effects of thermal load parameters, material and geometrical parameters on the quasi-static and dynamic responses of the beam are discussed. Numerical results show that the thermal function parameter has a great effect on quasi-static and dynamic responses of the beam. Compared with the thermal relaxation time, the initial vibrational responses of the beam are more sensitive to the mechanical relaxation time of the thermoviscoelastic material.展开更多
This paper presents a combined application of the finite element method (FEM) and the differential quadrature method (DQM) to vibration and buckling problems of rectangular plates. The proposed scheme combines the...This paper presents a combined application of the finite element method (FEM) and the differential quadrature method (DQM) to vibration and buckling problems of rectangular plates. The proposed scheme combines the geometry flexibility of the FEM and the high accuracy and efficiency of the DQM. The accuracy of the present method is demonstrated by comparing the obtained results with those available in the literature. It is shown that highly accurate results can be obtained by using a small number of finite elements and DQM sample points. The proposed method is suitable for the problems considered due to its simplicity and potential for further development.展开更多
This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence ...This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing partial differential equations (PDEs) are transformed into a system of coupled non-linear ordinary differential equations (ODEs) with appropriate boundary conditions for various physical parameters. The remaining system of ODEs is solved numerically using a differential transformation method (DTM). The effects of the porous parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and the Nusselt numbers are presented. Comparison of the obtained numerical results is made with previously published results in some special cases, with good agreement. The results obtained in this paper confirm the idea that DTM is a powerful mathematical tool and can be applied to a large class of linear and non-linear problems in different fields of science and engineering.展开更多
In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference me...In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.展开更多
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.展开更多
This paper deals with the numerical solution of initial value problems for systems of differential equations with a delay argument. The numerical stability of a linear multistep method is investigated by analysing the...This paper deals with the numerical solution of initial value problems for systems of differential equations with a delay argument. The numerical stability of a linear multistep method is investigated by analysing the solution of the lest equation y’(t)=Ay(t) + By(1-t),where A,B denote constant complex N×N-matrices,and t】0.We investigate carefully the characterization of the stability region.展开更多
An improved precise integration method (IPIM) for solving the differential Riccati equation (DRE) is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise in...An improved precise integration method (IPIM) for solving the differential Riccati equation (DRE) is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method (PIM) for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the展开更多
We extend the differential quadrature element method (DQEM) to the buckling analysis of uniformly in-plane loaded functionally graded (FG) plates fully or partially resting on the Pasternak model of elastic support. M...We extend the differential quadrature element method (DQEM) to the buckling analysis of uniformly in-plane loaded functionally graded (FG) plates fully or partially resting on the Pasternak model of elastic support. Material properties of the FG plate are graded in the thickness direction and assumed to obey a power law distribution of the volume fraction of the constituents. To set up the global eigenvalue equation, the plate is divided into sub-domains or elements and the generalized differential quadrature procedure is applied to discretize the governing, boundary and compatibility equations. By assembling discrete equations at all nodal points, the weighting coefficient and force matrices are derived. To validate the accuracy of this method, the results are compared with those of the exact solution and the finite element method. At the end, the effects of different variables and local elastic support arrangements on the buckling load factor are investigated.展开更多
In this article, a new application to find the exact solutions of nonlinear partial time-space fractional differential Equation has been discussed. Firstly, the fractional complex transformation has been implemented t...In this article, a new application to find the exact solutions of nonlinear partial time-space fractional differential Equation has been discussed. Firstly, the fractional complex transformation has been implemented to convert nonlinear partial fractional differential Equations into nonlinear ordinary differential Equations. Afterwards, the (G'/G)-expansion method has been implemented, to celebrate the exact solutions of these Equations, in the sense of modified Riemann-Liouville derivative. As application, the exact solutions of time-space fractional Burgers’ Equation have been discussed.展开更多
基金supported by a grant from the National Science and Technology Council of the Republic of China(Grant Number:MOST 112-2221-E-006-048-MY2).
文摘This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.
基金Anhui Provincial Natural Science Foundation(2308085QD124)Anhui Province University Natural Science Research Project(GrantNo.2023AH050918)The University Outstanding Youth Talent Support Program of Anhui Province.
文摘This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.
文摘Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.
文摘To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.
文摘The numerical approach for finding the solution of fractional order systems of boundary value problems (BPVs) is derived in this paper. The implementation of the weighted residuals such as Galerkin, Least Square, and Collocation methods are included for solving fractional order differential equations, which is broadened to acquire the approximate solutions of fractional order systems with differentiable polynomials, namely Legendre polynomials, as basis functions. The algorithm of the residual formulations of matrix form can be coded efficiently. The interpretation of Caputo fractional derivatives is employed here. We have demonstrated these methods numerically through a few examples of linear and nonlinear BVPs. The results in absolute errors show that the present method efficiently finds the numerical solutions of fractional order systems of differential equations.
文摘In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.
文摘The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect is considered. Detailed formulations and procedures are given for the first time. Several examples are analyzed and accurate results are obtained by the DQEM. Based on the results in this paper, one may conclude that the DQEM is a useful tool for obtaining solutions of structural elements. It can be seen that the shell shape may be theore tically controlled and snap through may occur when the applied voltage reaches a critical value even without mechanical load for certain geometric configurations.
文摘This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.
基金National Key Project of China (No.PD9521907)the National Science Foundation of China (No.19872025).
文摘It is a new attempt to extend the differential quadrature method(DQM) to stability analysis of the straight and curved centerlinepipes conveying fluid. Emphasis is placed on the study of theinfluences of several parameters on the critical flow velocity.Compared to other methods, this method can more easily deal with thepipe with spring support at its boundaries and asks for much lesscomputing effort while giving ac- ceptable precision in the numericalresults.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10272021, 10572021) and the Doctoral Program Foundation of Institution of Higher Education of China (Grant No 20040007022).
文摘The Hamilton-Jacobi method for solving ordinary differential equations is presented in this paper. A system of ordinary differential equations of first order or second order can be expressed as a Hamilton system under certain conditions. Then the Hamilton-Jacobi method is used in the integration of the Hamilton system and the solution of the original ordinary differential equations can be found. Finally, an example is given to illustrate the application of the result.
基金Project supported by the National Natural Science Foundation of China (Grant No 10272021) and the Doctorate Foundation of the State Education Ministry of China (Grant No 20040007022).
文摘A differential equation of first order can be expressed by the equation of motion of a mechanical system. In this paper, three methods of analytical mechanics, i.e. the Hamilton-Noether method, the Lagrange-Noether method and the Poisson method, are given to solve a differential equation of first order, of which the way may be called the mechanical methodology in mathematics.
基金supported by the National Natural Science Foundation of China(Nos.11772182 and90816001)
文摘The quasi-static and dynamic responses of a thermoviscoelastic Timoshenko beam subject to thermal loads are analyzed. First, based on the small geometric deformation assumption and Boltzmann constitutive relation, the governing equations for the beam are presented. Second, an extended differential quadrature method(DQM)in the spatial domain and a differential method in the temporal domain are combined to transform the integro-partial-differential governing equations into the ordinary differential equations. Third, the accuracy of the present discrete method is verified by elastic/viscoelastic examples, and the effects of thermal load parameters, material and geometrical parameters on the quasi-static and dynamic responses of the beam are discussed. Numerical results show that the thermal function parameter has a great effect on quasi-static and dynamic responses of the beam. Compared with the thermal relaxation time, the initial vibrational responses of the beam are more sensitive to the mechanical relaxation time of the thermoviscoelastic material.
文摘This paper presents a combined application of the finite element method (FEM) and the differential quadrature method (DQM) to vibration and buckling problems of rectangular plates. The proposed scheme combines the geometry flexibility of the FEM and the high accuracy and efficiency of the DQM. The accuracy of the present method is demonstrated by comparing the obtained results with those available in the literature. It is shown that highly accurate results can be obtained by using a small number of finite elements and DQM sample points. The proposed method is suitable for the problems considered due to its simplicity and potential for further development.
文摘This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing partial differential equations (PDEs) are transformed into a system of coupled non-linear ordinary differential equations (ODEs) with appropriate boundary conditions for various physical parameters. The remaining system of ODEs is solved numerically using a differential transformation method (DTM). The effects of the porous parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and the Nusselt numbers are presented. Comparison of the obtained numerical results is made with previously published results in some special cases, with good agreement. The results obtained in this paper confirm the idea that DTM is a powerful mathematical tool and can be applied to a large class of linear and non-linear problems in different fields of science and engineering.
基金heprojectissupportedbyNNSFofChina (No .1 9972 0 39) .
文摘In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.
基金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.
文摘This paper deals with the numerical solution of initial value problems for systems of differential equations with a delay argument. The numerical stability of a linear multistep method is investigated by analysing the solution of the lest equation y’(t)=Ay(t) + By(1-t),where A,B denote constant complex N×N-matrices,and t】0.We investigate carefully the characterization of the stability region.
基金Project supported by the National Natural Science Foundation of China(Nos.10902020 and 10721062)
文摘An improved precise integration method (IPIM) for solving the differential Riccati equation (DRE) is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method (PIM) for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the
文摘We extend the differential quadrature element method (DQEM) to the buckling analysis of uniformly in-plane loaded functionally graded (FG) plates fully or partially resting on the Pasternak model of elastic support. Material properties of the FG plate are graded in the thickness direction and assumed to obey a power law distribution of the volume fraction of the constituents. To set up the global eigenvalue equation, the plate is divided into sub-domains or elements and the generalized differential quadrature procedure is applied to discretize the governing, boundary and compatibility equations. By assembling discrete equations at all nodal points, the weighting coefficient and force matrices are derived. To validate the accuracy of this method, the results are compared with those of the exact solution and the finite element method. At the end, the effects of different variables and local elastic support arrangements on the buckling load factor are investigated.
文摘In this article, a new application to find the exact solutions of nonlinear partial time-space fractional differential Equation has been discussed. Firstly, the fractional complex transformation has been implemented to convert nonlinear partial fractional differential Equations into nonlinear ordinary differential Equations. Afterwards, the (G'/G)-expansion method has been implemented, to celebrate the exact solutions of these Equations, in the sense of modified Riemann-Liouville derivative. As application, the exact solutions of time-space fractional Burgers’ Equation have been discussed.