In this paper,we present the semi-implicit Euler(SIE)numerical solution for stochastic pantograph equations with jumps and prove that the SIE approximation solution converges to the exact solution in the mean-square...In this paper,we present the semi-implicit Euler(SIE)numerical solution for stochastic pantograph equations with jumps and prove that the SIE approximation solution converges to the exact solution in the mean-square sense under the Local Lipschitz condition.展开更多
Mineral exploration is done by different methods. Geophysical and geochemical studies are two powerful tools in this field. In integrated studies, the results of each study are used to determine the location of the dr...Mineral exploration is done by different methods. Geophysical and geochemical studies are two powerful tools in this field. In integrated studies, the results of each study are used to determine the location of the drilling boreholes. The purpose of this study is to use field geophysics to calculate the depth of mineral reserve. The study area is located 38 km from Zarand city called Jalalabad iron mine. In this study, gravimetric data were measured and mineral depth was calculated using the Euler method. 1314 readings have been performed in this area. The rocks of the region include volcanic and sedimentary. The source of the mineralization in the area is hydrothermal processes. After gravity measuring in the region, the data were corrected, then various methods such as anomalous map remaining in levels one and two, upward expansion, first and second-degree vertical derivatives, analytical method, and analytical signal were drawn, and finally, the depth of the deposit was estimated by Euler method. As a result, the depth of the mineral deposit was calculated to be between 20 and 30 meters on average.展开更多
To capitalize on the synergies between the Econometrics course and the Environmental Economics major,this paper aims to enhance students’ability to conduct empirical analysis and practical application using econometr...To capitalize on the synergies between the Econometrics course and the Environmental Economics major,this paper aims to enhance students’ability to conduct empirical analysis and practical application using econometric models.It also seeks to promote collaborative teaching through case studies and model research.The primary focus is on the hot research issues within the field of environmental economics,utilizing the econometric model as a vehicle for instruction.To achieve this,the paper proposes the development of a comprehensive case library specific to environmental economics.This resource will serve to optimize the case teaching approach,incorporating the use of econometric software,and fostering interactive teaching models between educators and students.By implementing these strategies,the paper outlines a path and mode for collaborative teaching that effectively bridges the gap between econometrics and environmental economics.展开更多
Based on the analysis methods of non-parametric Malmquist index and spatial econometrics as well as the provincial panel data in 2007-2010, this paper estimates the efficiency of fiscal expenditure from local governme...Based on the analysis methods of non-parametric Malmquist index and spatial econometrics as well as the provincial panel data in 2007-2010, this paper estimates the efficiency of fiscal expenditure from local governments in china in terms of reducing the income gap between urban and rural residents for the first time and evaluates the spatial correlation and heterogeneity of this efficiency. The results have shown that the fiscal expenditure of most provinces is of low efficiency in reducing the income gap between urban and rural residents, and the expenditure efficiency of local governments is not relevant to their levels of economic development. Besides, the efficiency on reducing the urban-rural income gap between different regions of China has a tendency of convergence. But this is mainly reflected inside the regional economic belt. There is significant difference between the efficiency of each economic belt. The central region has the highest efficiency in a rising trend, the western region has the lowest efficiency in a downward trend, while the eastern region is relatively stable.展开更多
A class of stochastic Besov spaces BpL^(2)(Ω;˙H^(α)(O)),1≤p≤∞andα∈[−2,2],is introduced to characterize the regularity of the noise in the semilinear stochastic heat equation du−Δudt=f(u)dt+dW(t),under the fol...A class of stochastic Besov spaces BpL^(2)(Ω;˙H^(α)(O)),1≤p≤∞andα∈[−2,2],is introduced to characterize the regularity of the noise in the semilinear stochastic heat equation du−Δudt=f(u)dt+dW(t),under the following conditions for someα∈(0,1]:||∫_(0)^(t)e−(t−s)^(A)dW(s)||L^(2)(Ω;L^(2)(O))≤C^(t^(α/2))and||∫_(0)^(t)e−(t−s)^(A)dW(s)||_B^(∞)L^(2)(Ω:H^(α)(O))≤C..The conditions above are shown to be satisfied by both trace-class noises(withα=1)and one-dimensional space-time white noises(withα=1/2).The latter would fail to satisfy the conditions withα=1/2 if the stochastic Besov norm||·||B∞L^(2)(Ω;˙H^(α)(O))is replaced by the classical Sobolev norm||·||L^(2)(Ω;˙Hα(O)),and this often causes reduction of the convergence order in the numerical analysis of the semilinear stochastic heat equation.In this paper,the convergence of a modified exponential Euler method,with a spectral method for spatial discretization,is proved to have orderαin both the time and space for possibly nonsmooth initial data in L^(4)(Ω;˙H^(β)(O))withβ>−1,by utilizing the real interpolation properties of the stochastic Besov spaces and a class of locally refined stepsizes to resolve the singularity of the solution at t=0.展开更多
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.展开更多
In the literature (Tan and Wang, 2010), Tan and Wang investigated the convergence of the split-step backward Euler (SSBE) method for linear stochastic delay integro-differential equations (SDIDEs) and proved the...In the literature (Tan and Wang, 2010), Tan and Wang investigated the convergence of the split-step backward Euler (SSBE) method for linear stochastic delay integro-differential equations (SDIDEs) and proved the mean-square stability of SSBE method under some condition. Unfortu- nately, the main result of stability derived by the condition is somewhat restrictive to be applied for practical application. This paper improves the corresponding results. The authors not only prove the mean-square stability of the numerical method but also prove the general mean-square stability of the numerical method. Furthermore, an example is given to illustrate the theory.展开更多
In this paper,the numerical methods for semi-linear stochastic delay integro-difFerential equations are studied.The uniqueness,existence and stability of analytic solutions of semi-linear stochastic delay integro-diff...In this paper,the numerical methods for semi-linear stochastic delay integro-difFerential equations are studied.The uniqueness,existence and stability of analytic solutions of semi-linear stochastic delay integro-differential equations are studied and some suitable conditions for the mean-square stability of the analytic solutions are also obtained.Then the numerical approximation of exponential Euler method for semi-linear stochastic delay integro-differential equations is constructed and the convergence and the stability of the numerical method are studied.It is proved that the exponential Euler method is convergent with strong order 1/2 and can keep the mean-square exponential stability of the analytical solutions under some restrictions on the step size.In addition,numerical experiments are presented to confirm the theoretical results.展开更多
The paper presents the comparative study on numerical methods of Euler method,Improved Euler method and fourth-order Runge-Kutta method for solving the engineering problems and applications.The three proposed methods ...The paper presents the comparative study on numerical methods of Euler method,Improved Euler method and fourth-order Runge-Kutta method for solving the engineering problems and applications.The three proposed methods are quite efficient and practically well suited for solving the unknown engineering problems.This paper aims to enhance the teaching and learning quality of teachers and students for various levels.At each point of the interval,the value of y is calculated and compared with its exact value at that point.The next interesting point is the observation of error from those methods.Error in the value of y is the difference between calculated and exact value.A mathematical equation which relates various functions with its derivatives is known as a differential equation.It is a popular field of mathematics because of its application to real-world problems.To calculate the exact values,the approximate values and the errors,the numerical tool such as MATLAB is appropriate for observing the results.This paper mainly concentrates on identifying the method which provides more accurate results.Then the analytical results and calculates their corresponding error were compared in details.The minimum error directly reflected to realize the best method from different numerical methods.According to the analyses from those three approaches,we observed that only the error is nominal for the fourth-order Runge-Kutta method.展开更多
Level Set interface treatment method is introduced into Euler method,which is employed for interface treatment method for multi-materials. Combined with the ghost fluid method,the moving interface is tracked. Fifth-or...Level Set interface treatment method is introduced into Euler method,which is employed for interface treatment method for multi-materials. Combined with the ghost fluid method,the moving interface is tracked. Fifth-order WENO spatial discretization and third-order TVD Runge-Kutta time discretization methods are used. Shock-wave action on bubble,implosion and velocity field Shock effect bubbles; implosion and velocity field are simulated by means of LS-MMIC3D programmed by C++. Nu-merical results show that the Level Set interface treatment method is effective and feasible for multi-material interface treatment in comparison with the WENO method.展开更多
In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume (RKCV) discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian co...In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume (RKCV) discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin (RKDG) method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.展开更多
Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions(symmetric and non-symmetric)are studied in this work.Reliability of the obtained result...Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions(symmetric and non-symmetric)are studied in this work.Reliability of the obtained results is verified by the finite difference method(FDM)and the finite element method(FEM)with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes(regular and non-regular).The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.展开更多
A preconditioned gridless method is developed for solving the Euler equations at low Mach numbers.The preconditioned system in a conservation form is obtained by multiplying apreconditioning matrix of the type of Weis...A preconditioned gridless method is developed for solving the Euler equations at low Mach numbers.The preconditioned system in a conservation form is obtained by multiplying apreconditioning matrix of the type of Weiss and Smith to the time derivative of the Euler equations,which are discretized using agridless technique wherein the physical domain is distributed by clouds of points.The implementation of the preconditioned gridless method is mainly based on the frame of the traditional gridless method without preconditioning,which may fail to converge for low Mach number simulations.Therefore,the modifications corresponding to the affected terms of preconditioning are mainly addressed.The numerical results show that the preconditioned gridless method still functions for compressible transonic flow simulations and additionally,for nearly incompressible flow simulations at low Mach numbers as well.The paper ends with the nearly incompressible flow over a multi-element airfoil,which demonstrates the ability of the method presented for treating flows over complicated geometries.展开更多
The numerical solutions of standing waves for Euler equations with the nonlinear free surface boundary condition in a two-dimensional (2D) tank are studied. The irregular tank is mapped onto a fixed square domain th...The numerical solutions of standing waves for Euler equations with the nonlinear free surface boundary condition in a two-dimensional (2D) tank are studied. The irregular tank is mapped onto a fixed square domain through proper mapping functions. A staggered mesh system is employed in a 2D tank to calculate the elevation of the transient fluid. A time-independent finite difference method, which is developed by Bang- fuh Chen, is used to solve the Euler equations for incompressible and inviscid fluids. The numerical results agree well with the analytic solutions and previously published results. The sloshing profiles of surge and heave motion with initial standing waves are presented. The results show very clear nonlinear and beating phenomena.展开更多
This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification o...This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional. Using this method, a rapid convergent sequence is produced which converges to the exact solutions of the problem. Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.展开更多
This paper develops the mean-square exponential input-to-state stability(exp-ISS) of the Euler-Maruyama(EM) method for stochastic delay control systems(SDCSs).The definition of mean-square exp-ISS of numerical m...This paper develops the mean-square exponential input-to-state stability(exp-ISS) of the Euler-Maruyama(EM) method for stochastic delay control systems(SDCSs).The definition of mean-square exp-ISS of numerical methods is established.The conditions of the exact and EM method for an SDCS with the property of mean-square exp-ISS are obtained without involving control Lyapunov functions or functional.Under the global Lipschitz coefficients and mean-square continuous measurable inputs,it is proved that the mean-square exp-ISS of an SDCS holds if and only if that of the EM method is preserved for a sufficiently small step size.The proposed results are evaluated by using numerical experiments to show their effectiveness.展开更多
In this work, we apply the Zhou’s method [1] or differential transformation method (DTM) for solving the Euler equidimensional equation. The Zhou’s method may be considered as alternative and efficient for finding t...In this work, we apply the Zhou’s method [1] or differential transformation method (DTM) for solving the Euler equidimensional equation. The Zhou’s method may be considered as alternative and efficient for finding the approximate solutions of initial values problems. We prove superiority of this method by applying them on the some Euler type equation, in this case of order 2 and 3 [2]. The power series solution of the reduced equation transforms into an approximate implicit solution of the original equations. The results agreed with the exact solution obtained via transformation to a constant coefficient equation.展开更多
基金Supported by the NSF of the Higher Education Institutions of Jiangsu Province(10KJD110006)Supported by the grant of Jiangsu Institute of Education(Jsjy2009zd03)Supported by the Qing Lan Project of Jiangsu Province(2010)
文摘In this paper,we present the semi-implicit Euler(SIE)numerical solution for stochastic pantograph equations with jumps and prove that the SIE approximation solution converges to the exact solution in the mean-square sense under the Local Lipschitz condition.
文摘Mineral exploration is done by different methods. Geophysical and geochemical studies are two powerful tools in this field. In integrated studies, the results of each study are used to determine the location of the drilling boreholes. The purpose of this study is to use field geophysics to calculate the depth of mineral reserve. The study area is located 38 km from Zarand city called Jalalabad iron mine. In this study, gravimetric data were measured and mineral depth was calculated using the Euler method. 1314 readings have been performed in this area. The rocks of the region include volcanic and sedimentary. The source of the mineralization in the area is hydrothermal processes. After gravity measuring in the region, the data were corrected, then various methods such as anomalous map remaining in levels one and two, upward expansion, first and second-degree vertical derivatives, analytical method, and analytical signal were drawn, and finally, the depth of the deposit was estimated by Euler method. As a result, the depth of the mineral deposit was calculated to be between 20 and 30 meters on average.
基金supported by the Ministry of Education of Humanities and Social Science Project(21YJC630009)the National Natural Science Foundation of China(No.72104116).
文摘To capitalize on the synergies between the Econometrics course and the Environmental Economics major,this paper aims to enhance students’ability to conduct empirical analysis and practical application using econometric models.It also seeks to promote collaborative teaching through case studies and model research.The primary focus is on the hot research issues within the field of environmental economics,utilizing the econometric model as a vehicle for instruction.To achieve this,the paper proposes the development of a comprehensive case library specific to environmental economics.This resource will serve to optimize the case teaching approach,incorporating the use of econometric software,and fostering interactive teaching models between educators and students.By implementing these strategies,the paper outlines a path and mode for collaborative teaching that effectively bridges the gap between econometrics and environmental economics.
基金Supported by National Science Fund for Distinguished Young Scholars(GrantNo.:70825003)Key Project of National Social Science Foundation(GrantNo.:07AJL002,12AGL008 and 12ASH004)+3 种基金Young Scholar Project of National Social Science Foundation(Grant No.:12CGL063 and 12CJY062)Key Project of Ministry of Education(Grant No.:DFA100209)Social Science Planning Fund of Ministry of Education (Grant No.:07JA790104)Foundation Project for Central Universities-Xiamen University(Grant No. :2009ZK1007)
文摘Based on the analysis methods of non-parametric Malmquist index and spatial econometrics as well as the provincial panel data in 2007-2010, this paper estimates the efficiency of fiscal expenditure from local governments in china in terms of reducing the income gap between urban and rural residents for the first time and evaluates the spatial correlation and heterogeneity of this efficiency. The results have shown that the fiscal expenditure of most provinces is of low efficiency in reducing the income gap between urban and rural residents, and the expenditure efficiency of local governments is not relevant to their levels of economic development. Besides, the efficiency on reducing the urban-rural income gap between different regions of China has a tendency of convergence. But this is mainly reflected inside the regional economic belt. There is significant difference between the efficiency of each economic belt. The central region has the highest efficiency in a rising trend, the western region has the lowest efficiency in a downward trend, while the eastern region is relatively stable.
基金supported by National Natural Science Foundation of China(Grant Nos.12071020,12131005 and U2230402)the Research Grants Council of Hong Kong(Grant No.Poly U15300519)an Internal Grant of The Hong Kong Polytechnic University(Grant No.P0038843,Work Programme:ZVX7)。
文摘A class of stochastic Besov spaces BpL^(2)(Ω;˙H^(α)(O)),1≤p≤∞andα∈[−2,2],is introduced to characterize the regularity of the noise in the semilinear stochastic heat equation du−Δudt=f(u)dt+dW(t),under the following conditions for someα∈(0,1]:||∫_(0)^(t)e−(t−s)^(A)dW(s)||L^(2)(Ω;L^(2)(O))≤C^(t^(α/2))and||∫_(0)^(t)e−(t−s)^(A)dW(s)||_B^(∞)L^(2)(Ω:H^(α)(O))≤C..The conditions above are shown to be satisfied by both trace-class noises(withα=1)and one-dimensional space-time white noises(withα=1/2).The latter would fail to satisfy the conditions withα=1/2 if the stochastic Besov norm||·||B∞L^(2)(Ω;˙H^(α)(O))is replaced by the classical Sobolev norm||·||L^(2)(Ω;˙Hα(O)),and this often causes reduction of the convergence order in the numerical analysis of the semilinear stochastic heat equation.In this paper,the convergence of a modified exponential Euler method,with a spectral method for spatial discretization,is proved to have orderαin both the time and space for possibly nonsmooth initial data in L^(4)(Ω;˙H^(β)(O))withβ>−1,by utilizing the real interpolation properties of the stochastic Besov spaces and a class of locally refined stepsizes to resolve the singularity of the solution at t=0.
文摘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.
基金supported by the Fundamental Research Funds for the Central Universities under Grant No. 2012089:31541111213China Postdoctoral Science Foundation Funded Project under Grant No.2012M511615the State Key Program of National Natural Science of China under Grant No.61134012
文摘In the literature (Tan and Wang, 2010), Tan and Wang investigated the convergence of the split-step backward Euler (SSBE) method for linear stochastic delay integro-differential equations (SDIDEs) and proved the mean-square stability of SSBE method under some condition. Unfortu- nately, the main result of stability derived by the condition is somewhat restrictive to be applied for practical application. This paper improves the corresponding results. The authors not only prove the mean-square stability of the numerical method but also prove the general mean-square stability of the numerical method. Furthermore, an example is given to illustrate the theory.
基金This research is supported by National Natural Science Foundation of China(Project No.11901173)by the Heilongjiang province Natural Science Foundation(LH2019A030)by the Heilongjiang province Innovation Talent Foundation(2018CX17).
文摘In this paper,the numerical methods for semi-linear stochastic delay integro-difFerential equations are studied.The uniqueness,existence and stability of analytic solutions of semi-linear stochastic delay integro-differential equations are studied and some suitable conditions for the mean-square stability of the analytic solutions are also obtained.Then the numerical approximation of exponential Euler method for semi-linear stochastic delay integro-differential equations is constructed and the convergence and the stability of the numerical method are studied.It is proved that the exponential Euler method is convergent with strong order 1/2 and can keep the mean-square exponential stability of the analytical solutions under some restrictions on the step size.In addition,numerical experiments are presented to confirm the theoretical results.
文摘The paper presents the comparative study on numerical methods of Euler method,Improved Euler method and fourth-order Runge-Kutta method for solving the engineering problems and applications.The three proposed methods are quite efficient and practically well suited for solving the unknown engineering problems.This paper aims to enhance the teaching and learning quality of teachers and students for various levels.At each point of the interval,the value of y is calculated and compared with its exact value at that point.The next interesting point is the observation of error from those methods.Error in the value of y is the difference between calculated and exact value.A mathematical equation which relates various functions with its derivatives is known as a differential equation.It is a popular field of mathematics because of its application to real-world problems.To calculate the exact values,the approximate values and the errors,the numerical tool such as MATLAB is appropriate for observing the results.This paper mainly concentrates on identifying the method which provides more accurate results.Then the analytical results and calculates their corresponding error were compared in details.The minimum error directly reflected to realize the best method from different numerical methods.According to the analyses from those three approaches,we observed that only the error is nominal for the fourth-order Runge-Kutta method.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10872085 and 10472042)Program for New Century Excellent Talents of Ministry of Education (Grant No. NCET-08-0043)+1 种基金National Basic Research Program of China (Grant No. 2010CB832700 6)Key Program of Numerical Simulation Key Laboratory of Sichuan Higher School (Grant No. 07NJZZ001)
文摘Level Set interface treatment method is introduced into Euler method,which is employed for interface treatment method for multi-materials. Combined with the ghost fluid method,the moving interface is tracked. Fifth-order WENO spatial discretization and third-order TVD Runge-Kutta time discretization methods are used. Shock-wave action on bubble,implosion and velocity field Shock effect bubbles; implosion and velocity field are simulated by means of LS-MMIC3D programmed by C++. Nu-merical results show that the Level Set interface treatment method is effective and feasible for multi-material interface treatment in comparison with the WENO method.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11261035 and 11171038)the Science Research Foundation of the Institute of Higher Education of Inner Mongolia Autonomous Region, China (Grant No. NJZZ12198)the Natural Science Foundation of Inner Mongolia Autonomous Region, China (Grant No. 2012MS0102)
文摘In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume (RKCV) discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin (RKDG) method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.
文摘Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions(symmetric and non-symmetric)are studied in this work.Reliability of the obtained results is verified by the finite difference method(FDM)and the finite element method(FEM)with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes(regular and non-regular).The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.
基金supported by the National Natural Science Foundation of China(No.11172134)
文摘A preconditioned gridless method is developed for solving the Euler equations at low Mach numbers.The preconditioned system in a conservation form is obtained by multiplying apreconditioning matrix of the type of Weiss and Smith to the time derivative of the Euler equations,which are discretized using agridless technique wherein the physical domain is distributed by clouds of points.The implementation of the preconditioned gridless method is mainly based on the frame of the traditional gridless method without preconditioning,which may fail to converge for low Mach number simulations.Therefore,the modifications corresponding to the affected terms of preconditioning are mainly addressed.The numerical results show that the preconditioned gridless method still functions for compressible transonic flow simulations and additionally,for nearly incompressible flow simulations at low Mach numbers as well.The paper ends with the nearly incompressible flow over a multi-element airfoil,which demonstrates the ability of the method presented for treating flows over complicated geometries.
文摘The numerical solutions of standing waves for Euler equations with the nonlinear free surface boundary condition in a two-dimensional (2D) tank are studied. The irregular tank is mapped onto a fixed square domain through proper mapping functions. A staggered mesh system is employed in a 2D tank to calculate the elevation of the transient fluid. A time-independent finite difference method, which is developed by Bang- fuh Chen, is used to solve the Euler equations for incompressible and inviscid fluids. The numerical results agree well with the analytic solutions and previously published results. The sloshing profiles of surge and heave motion with initial standing waves are presented. The results show very clear nonlinear and beating phenomena.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10771019 and 10826107)
文摘This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional. Using this method, a rapid convergent sequence is produced which converges to the exact solutions of the problem. Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.
基金Supported by National Natural Science Foundation of China(10571036)the Key Discipline Development Program of Beijing Municipal Commission (XK100080537)
基金supported by the National Natural Science Foundation of China(6127312660904032)the Natural Science Foundation of Guangdong Province(10251064101000008)
文摘This paper develops the mean-square exponential input-to-state stability(exp-ISS) of the Euler-Maruyama(EM) method for stochastic delay control systems(SDCSs).The definition of mean-square exp-ISS of numerical methods is established.The conditions of the exact and EM method for an SDCS with the property of mean-square exp-ISS are obtained without involving control Lyapunov functions or functional.Under the global Lipschitz coefficients and mean-square continuous measurable inputs,it is proved that the mean-square exp-ISS of an SDCS holds if and only if that of the EM method is preserved for a sufficiently small step size.The proposed results are evaluated by using numerical experiments to show their effectiveness.
文摘In this work, we apply the Zhou’s method [1] or differential transformation method (DTM) for solving the Euler equidimensional equation. The Zhou’s method may be considered as alternative and efficient for finding the approximate solutions of initial values problems. We prove superiority of this method by applying them on the some Euler type equation, in this case of order 2 and 3 [2]. The power series solution of the reduced equation transforms into an approximate implicit solution of the original equations. The results agreed with the exact solution obtained via transformation to a constant coefficient equation.