In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introduc...In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introducing only one additional evaluation of the function in each step.Furthermore,some new efficient methods with a higher-order of convergence are obtained by using only a single matrix inversion in each iteration.Analyses of convergence properties and computational efficiency of these new methods are made and testified by several numerical problems.By comparison,the new schemes are more efficient than the corresponding existing ones,particularly for large problem sizes.展开更多
In this paper,we construct a new sixth order iterative method for solving nonlinear equations.The local convergence and order of convergence of the new iterative method is demonstrated.In order to check the validity o...In this paper,we construct a new sixth order iterative method for solving nonlinear equations.The local convergence and order of convergence of the new iterative method is demonstrated.In order to check the validity of the new iterative method,we employ several chemical engineering applications and academic test problems.Numerical results show the good numerical performance of the new iterative method.Moreover,the dynamical study of the new method also supports the theoretical results.展开更多
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.展开更多
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.展开更多
Thymus quinquecostatus Celak.,a traditional aromatic edible plant from Lamiaceae,is widely used as food additive,food condiment,spice,and herbal teas.Polyphenol-rich fraction of T.quinquecostatus(PRF)has been proven t...Thymus quinquecostatus Celak.,a traditional aromatic edible plant from Lamiaceae,is widely used as food additive,food condiment,spice,and herbal teas.Polyphenol-rich fraction of T.quinquecostatus(PRF)has been proven to be effective protective effect for cerebral ischemia reperfusion injury(CIRI)in our previous study.In this study,we developed a novel“Gut flora-Compound-Target-Pathway”(GCTP)network based on network pharmacology coupled with gastrointestinal metabolism for screening bio-active components,key targets and gut floras through the classical technique for order preference by similarity to ideal solution(TOPSIS).This compensates for the lack of gut floras and gastrointestinal metabolism in network pharmacology.Firstly,four incubation models covering simulated gastric juice,simulated intestinal juice,gut floras of normal and transient middle cerebral artery occlusion(tMCAO)rat in vitro were applied to PRF.The 109 proto-components and 64 metabolites were elucidated by ultra-high performance liquid chromatography Q exactive orbitrap-mass spectrometry(UPLC-QE-Orbitrap-MS).Then,the key targets of matrix metalloproteinase 9(MMP9),prostaglandin-endoperoxide synthase 2(PTGS2),tyrosine-protein kinase fyn(FYN),estrogen receptor 1(ESR1),amyloid precursor protein(APP),and epidermal growth factor receptor(EGFR),and gut floras of Enterococcus avium LY1 were selected.Moreover,the selected key proto components were rosmarinic acid,daidzein,quercetin,luteolin,apigenin,methyl rosmarinate,kaempferol,luteoloside,and caffeic acid,and the key metabolites were isokaempferide,isorhamnetin,isoquercetin,and mangiferin.Binding of compounds to the key proteins was analyzed by molecular docking,and also verified though an 2,2'-azobis(2-amidinopropane)dihydrochloride(AAPH)induced oxidative stress zebrafish model and real-time quantitative polymerase chain reaction(RT-qPCR)assays.This study provides a new idea and a better understanding of PRF for its protective effects on CIRI and its underlying mechanisms.展开更多
In this paper, we study the solutions for variable-order time-fractional diffusion equations. A three-point combined compact difference (CCD) method is used to discretize the spatial variables to achieve sixth-order a...In this paper, we study the solutions for variable-order time-fractional diffusion equations. A three-point combined compact difference (CCD) method is used to discretize the spatial variables to achieve sixth-order accuracy, while the exponential-sum-approximation (ESA) is used to approximate the variable-order Caputo fractional derivative in the temporal direction, and a novel spatial sixth-order hybrid ESA-CCD method is implemented successfully. Finally, the accuracy of the proposed method is verified by numerical experiments.展开更多
In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident bou...In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident boundary in the HOS method. Based on the numerical model, the effects of wave parameters, such as the assumed focused amplitude, the central frequency, the frequency bandwidth, the wave amplitude distribution and the directional spreading on the surface elevation of the focused wave, the maximum generated wave crest, and the shifting of the focusing point, are numerically investigated. Especially, the effects of the wave directionality on the focused wave properties are emphasized. The numerical results show that the shifting of the focusing point and the maximum crest of the wave group are dependent on the amplitude of the focused wave, the central frequency, and the wave amplitude distribution type. The wave directionality has a definite effect on multidirectional focused waves. Generally, it can even out the difference between the simulated wave amplitude and the amplitude expected from theory and reduce the shifting of the focusing points, implying that the higher order interaction has an influence on wave focusing, especially for 2D wave. In 3D wave groups, a broader directional spreading weakens the higher nonlinear interactions.展开更多
The modified AOR method for solving linear complementarity problem(LCP(M,p))was proposed in literature,with some convergence results.In this paper,we considered the MAOR method for generalized-order linear complementa...The modified AOR method for solving linear complementarity problem(LCP(M,p))was proposed in literature,with some convergence results.In this paper,we considered the MAOR method for generalized-order linear complementarity problem(ELCP(M,N,p,q)),where M,N are nonsingular matrices of the following form:M=[D11H1K1D2],N=[D12H2K2D22],D11,D12,D21 and D22 are square nonsingular diagonal matrices.展开更多
In this paper, we present and analyze a family of fifth-order iterative methods free from second derivative for solving nonlinear equations. It is established that the family of iterative methods has convergence order...In this paper, we present and analyze a family of fifth-order iterative methods free from second derivative for solving nonlinear equations. It is established that the family of iterative methods has convergence order five. Numerical examples show that the new methods are comparable with the well known existing methods and give better results in many aspects.展开更多
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.展开更多
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.展开更多
A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method o...A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method on triangular grids. The developed HWENO methodology utilizes high-order derivative information to keep WENO re- construction stencils in the von Neumann neighborhood. A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils, making higher-order scheme stable and simplifying the reconstruction process at the same time. The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement. Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy, the designed HWENO limiters can simul- taneously obtain uniform high order accuracy and sharp, es- sentially non-oscillatory shock transition.展开更多
Owing to the Benjamin-Feir instability, the Stokes wave train experiences a modulation-demodulation process, and presents a recurrence characteristics. Stiassnie and Shemer researched the unstable evolution process an...Owing to the Benjamin-Feir instability, the Stokes wave train experiences a modulation-demodulation process, and presents a recurrence characteristics. Stiassnie and Shemer researched the unstable evolution process and provided a theoretical formulation for the recurrence period in 1985 on the basis of the nonlinear cubic Schrodinger equation （NLS）. However, NLS has limitations on the narrow band and the weak nonlinearity. The recurrence period is re-investigated in this paper by using a highly efficient High Order Spectral （HOS） method, which can be applied for the direct phase- resolved simulation of the nonlinear wave train evolution. It is found that the Stiassnie and Shemer＇s formula should be modified in the cases with most unstable initial conditions, which is important for such topics as the generation mechanisms of freak waves. A new recurrence period formula is presented and some new evolution characteristics of the Stokes wave train are also discussed in details.展开更多
Based on the newly-developed element energy projection （EEP） method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite ele...Based on the newly-developed element energy projection （EEP） method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method （FEM） is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.展开更多
In this study,the problem of time-optimal reconnaissance trajectory design for the aeroassisted vehicle is considered.Different from most works reported previously,we explore the feasibility of applying a high-order a...In this study,the problem of time-optimal reconnaissance trajectory design for the aeroassisted vehicle is considered.Different from most works reported previously,we explore the feasibility of applying a high-order aeroassisted vehicle dynamic model to plan the optimal flight trajectory such that the gap between the simulated model and the real system can be narrowed.A highly-constrained optimal control model containing six-degree-of-freedom vehicle dynamics is established.To solve the formulated high-order trajectory planning model,a pipelined optimization strategy is illustrated.This approach is based on the variable order Radau pseudospectral method,indicating that the mesh grid used for discretizing the continuous system experiences several adaption iterations.Utilization of such a strategy can potentially smooth the flight trajectory and improve the algorithm convergence ability.Numerical simulations are reported to demonstrate some key features of the optimized flight trajectory.A number of comparative studies are also provided to verify the effectiveness of the applied method as well as the high-order trajectory planning model.展开更多
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.展开更多
Routh order reduction method of the relativistic Birkhoffian equations is studied. For a relativistic Birkhoffian system, the cyclic integrals can be found by using the perfect differential method. Through these cycli...Routh order reduction method of the relativistic Birkhoffian equations is studied. For a relativistic Birkhoffian system, the cyclic integrals can be found by using the perfect differential method. Through these cyclic integrals, the order of the system can be reduced. If the relativistic Birkhoffian system has a cyclic integral, then the Birkhoffian equations can be reduced at least by two degrees and the Birkhoffian form can be kept. The relations among the relativistic Birkhoffian mechanics, the relativistic Hamiltonian mechanics, and the relativistic Lagrangian mechanics are discussed, and the Routh order reduction method of the relativistic Lagrangian system is obtained. And an example is given to illustrate the application of the result.展开更多
In this paper, we establish two new iterative methods of order four and five by using modified homotopy perturbation technique. We also present the convergence analysis of these iterative methods. To assess the validi...In this paper, we establish two new iterative methods of order four and five by using modified homotopy perturbation technique. We also present the convergence analysis of these iterative methods. To assess the validity and performance of these iterative methods, we have applied to solve some nonlinear problems.展开更多
A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be ...A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be applied to a variety of velocity potential problems. The present paper, however, focused on its application to the problem of water entry of a wedge with varying speed. The continuity of the velocity achieved herein is particularly important for this kind of nonlinear free surface flow problem, because when the time stepping method is used, the free surface is updated through the velocity obtained at each node and the accuracy of the velocity is therefore crucial. Calculation was made for a case when the distance S that the wedge has travelled and time t follow the relationship s=Dtα, where D and α are constants, which is found to lead to a self similar flow field when the effect due to gravity is ignored.展开更多
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.展开更多
基金Supported by the National Natural Science Foundation of China(12061048)NSF of Jiangxi Province(20232BAB201026,20232BAB201018)。
文摘In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introducing only one additional evaluation of the function in each step.Furthermore,some new efficient methods with a higher-order of convergence are obtained by using only a single matrix inversion in each iteration.Analyses of convergence properties and computational efficiency of these new methods are made and testified by several numerical problems.By comparison,the new schemes are more efficient than the corresponding existing ones,particularly for large problem sizes.
基金supported by the National Natural Science Foundation of China (No.12271518)the Key Program of the National Natural Science Foundation of China (No.62333016)。
文摘In this paper,we construct a new sixth order iterative method for solving nonlinear equations.The local convergence and order of convergence of the new iterative method is demonstrated.In order to check the validity of the new iterative method,we employ several chemical engineering applications and academic test problems.Numerical results show the good numerical performance of the new iterative method.Moreover,the dynamical study of the new method also supports the theoretical results.
文摘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.
文摘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.
基金supported by the Key Research and Development Project of Ningxia Hui Autonomous Region(2020BFG03007).
文摘Thymus quinquecostatus Celak.,a traditional aromatic edible plant from Lamiaceae,is widely used as food additive,food condiment,spice,and herbal teas.Polyphenol-rich fraction of T.quinquecostatus(PRF)has been proven to be effective protective effect for cerebral ischemia reperfusion injury(CIRI)in our previous study.In this study,we developed a novel“Gut flora-Compound-Target-Pathway”(GCTP)network based on network pharmacology coupled with gastrointestinal metabolism for screening bio-active components,key targets and gut floras through the classical technique for order preference by similarity to ideal solution(TOPSIS).This compensates for the lack of gut floras and gastrointestinal metabolism in network pharmacology.Firstly,four incubation models covering simulated gastric juice,simulated intestinal juice,gut floras of normal and transient middle cerebral artery occlusion(tMCAO)rat in vitro were applied to PRF.The 109 proto-components and 64 metabolites were elucidated by ultra-high performance liquid chromatography Q exactive orbitrap-mass spectrometry(UPLC-QE-Orbitrap-MS).Then,the key targets of matrix metalloproteinase 9(MMP9),prostaglandin-endoperoxide synthase 2(PTGS2),tyrosine-protein kinase fyn(FYN),estrogen receptor 1(ESR1),amyloid precursor protein(APP),and epidermal growth factor receptor(EGFR),and gut floras of Enterococcus avium LY1 were selected.Moreover,the selected key proto components were rosmarinic acid,daidzein,quercetin,luteolin,apigenin,methyl rosmarinate,kaempferol,luteoloside,and caffeic acid,and the key metabolites were isokaempferide,isorhamnetin,isoquercetin,and mangiferin.Binding of compounds to the key proteins was analyzed by molecular docking,and also verified though an 2,2'-azobis(2-amidinopropane)dihydrochloride(AAPH)induced oxidative stress zebrafish model and real-time quantitative polymerase chain reaction(RT-qPCR)assays.This study provides a new idea and a better understanding of PRF for its protective effects on CIRI and its underlying mechanisms.
文摘In this paper, we study the solutions for variable-order time-fractional diffusion equations. A three-point combined compact difference (CCD) method is used to discretize the spatial variables to achieve sixth-order accuracy, while the exponential-sum-approximation (ESA) is used to approximate the variable-order Caputo fractional derivative in the temporal direction, and a novel spatial sixth-order hybrid ESA-CCD method is implemented successfully. Finally, the accuracy of the proposed method is verified by numerical experiments.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51309050 and 51221961)the National Basic Research Program of China(973 Program,Grant Nos.2013CB036101 and 2011CB013703)
文摘In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident boundary in the HOS method. Based on the numerical model, the effects of wave parameters, such as the assumed focused amplitude, the central frequency, the frequency bandwidth, the wave amplitude distribution and the directional spreading on the surface elevation of the focused wave, the maximum generated wave crest, and the shifting of the focusing point, are numerically investigated. Especially, the effects of the wave directionality on the focused wave properties are emphasized. The numerical results show that the shifting of the focusing point and the maximum crest of the wave group are dependent on the amplitude of the focused wave, the central frequency, and the wave amplitude distribution type. The wave directionality has a definite effect on multidirectional focused waves. Generally, it can even out the difference between the simulated wave amplitude and the amplitude expected from theory and reduce the shifting of the focusing points, implying that the higher order interaction has an influence on wave focusing, especially for 2D wave. In 3D wave groups, a broader directional spreading weakens the higher nonlinear interactions.
文摘The modified AOR method for solving linear complementarity problem(LCP(M,p))was proposed in literature,with some convergence results.In this paper,we considered the MAOR method for generalized-order linear complementarity problem(ELCP(M,N,p,q)),where M,N are nonsingular matrices of the following form:M=[D11H1K1D2],N=[D12H2K2D22],D11,D12,D21 and D22 are square nonsingular diagonal matrices.
文摘In this paper, we present and analyze a family of fifth-order iterative methods free from second derivative for solving nonlinear equations. It is established that the family of iterative methods has convergence order five. Numerical examples show that the new methods are comparable with the well known existing methods and give better results in many aspects.
基金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.
文摘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.
基金supported by the National Basic Research Program of China (2009CB724104)the National Natural Science Foundation of China (90716010)
文摘A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method on triangular grids. The developed HWENO methodology utilizes high-order derivative information to keep WENO re- construction stencils in the von Neumann neighborhood. A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils, making higher-order scheme stable and simplifying the reconstruction process at the same time. The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement. Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy, the designed HWENO limiters can simul- taneously obtain uniform high order accuracy and sharp, es- sentially non-oscillatory shock transition.
基金supported by the National Natural Science Foundation of China (Grant No. 41106001)the Specialized Research Fund for the Doctoral Program of Higher Education (Grant No. 20100094110016)+1 种基金the Special Research Funding of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering (Grant No. 2009585812)the Priority Academic Program Development of Jiangsu Higher Education Institutions (Coastal Development and Conservancy)
文摘Owing to the Benjamin-Feir instability, the Stokes wave train experiences a modulation-demodulation process, and presents a recurrence characteristics. Stiassnie and Shemer researched the unstable evolution process and provided a theoretical formulation for the recurrence period in 1985 on the basis of the nonlinear cubic Schrodinger equation （NLS）. However, NLS has limitations on the narrow band and the weak nonlinearity. The recurrence period is re-investigated in this paper by using a highly efficient High Order Spectral （HOS） method, which can be applied for the direct phase- resolved simulation of the nonlinear wave train evolution. It is found that the Stiassnie and Shemer＇s formula should be modified in the cases with most unstable initial conditions, which is important for such topics as the generation mechanisms of freak waves. A new recurrence period formula is presented and some new evolution characteristics of the Stokes wave train are also discussed in details.
基金the National Natural Science Foundation of China(No.50678093)Program for Changjiang Scholars and Innovative Research Team in University(No.IRT00736)
文摘Based on the newly-developed element energy projection （EEP） method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method （FEM） is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.
文摘In this study,the problem of time-optimal reconnaissance trajectory design for the aeroassisted vehicle is considered.Different from most works reported previously,we explore the feasibility of applying a high-order aeroassisted vehicle dynamic model to plan the optimal flight trajectory such that the gap between the simulated model and the real system can be narrowed.A highly-constrained optimal control model containing six-degree-of-freedom vehicle dynamics is established.To solve the formulated high-order trajectory planning model,a pipelined optimization strategy is illustrated.This approach is based on the variable order Radau pseudospectral method,indicating that the mesh grid used for discretizing the continuous system experiences several adaption iterations.Utilization of such a strategy can potentially smooth the flight trajectory and improve the algorithm convergence ability.Numerical simulations are reported to demonstrate some key features of the optimized flight trajectory.A number of comparative studies are also provided to verify the effectiveness of the applied method as well as the high-order trajectory planning model.
文摘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.
基金The project supported by National Natural Science Foundation of China under Grant Nos, 10372053 and 10472040, the Natural Science Foundation of Hunan Province under Grant No. 03JJY3005, the Scientific Research Foundation of Eduction Burean of Hunan Province under Grant No. 02C033 and the 0utstanding Young Talents Training Fund of Liaoning Province under Grant No. 3040005
文摘Routh order reduction method of the relativistic Birkhoffian equations is studied. For a relativistic Birkhoffian system, the cyclic integrals can be found by using the perfect differential method. Through these cyclic integrals, the order of the system can be reduced. If the relativistic Birkhoffian system has a cyclic integral, then the Birkhoffian equations can be reduced at least by two degrees and the Birkhoffian form can be kept. The relations among the relativistic Birkhoffian mechanics, the relativistic Hamiltonian mechanics, and the relativistic Lagrangian mechanics are discussed, and the Routh order reduction method of the relativistic Lagrangian system is obtained. And an example is given to illustrate the application of the result.
文摘In this paper, we establish two new iterative methods of order four and five by using modified homotopy perturbation technique. We also present the convergence analysis of these iterative methods. To assess the validity and performance of these iterative methods, we have applied to solve some nonlinear problems.
文摘A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be applied to a variety of velocity potential problems. The present paper, however, focused on its application to the problem of water entry of a wedge with varying speed. The continuity of the velocity achieved herein is particularly important for this kind of nonlinear free surface flow problem, because when the time stepping method is used, the free surface is updated through the velocity obtained at each node and the accuracy of the velocity is therefore crucial. Calculation was made for a case when the distance S that the wedge has travelled and time t follow the relationship s=Dtα, where D and α are constants, which is found to lead to a self similar flow field when the effect due to gravity is ignored.
文摘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.
