Solar-powered aircraft have attracted great attention owing to their potential for longendurance flight and wide application prospects.Due to the particularity of energy system,flight strategy optimization is a signif...Solar-powered aircraft have attracted great attention owing to their potential for longendurance flight and wide application prospects.Due to the particularity of energy system,flight strategy optimization is a significant way to enhance the flight performance for solar-powered aircraft.In this study,a flight strategy optimization model for high-altitude long-endurance solar-powered aircraft was proposed.This model consists of three-dimensional kinematic model,aerodynamic model,energy collection model,energy store model and energy loss model.To solve the nonlinear optimal control problem with process constraints and terminal constraints,Gauss pseudo-spectral method was employed to discretize the state equations and constraint equations.Then a typical mission flying from given initial point to given final point within a time interval was considered.Results indicate that proper changes of the attitude angle contribute to increasing the energy gained by photovoltaic cells.Utilization of gravitational potential energy can partly take the role of battery pack.Integrating these two measures,the optimized flight strategy can improve the final state of charge compared with current constant-altitude constant-velocity strategy.The optimized strategy brings more profits on condition of lower sunlight intensity and shorter daytime.展开更多
Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. ...Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. In this paper, an efficient numerical method for solving the fractional Advection-dispersion equation (ADE) is considered. The fractional derivative is described in the Caputo sense. The method is based on Chebyshev approximations. The properties of Chebyshev polynomials are used to reduce ADE to a system of ordinary differential equations, which are solved using the finite difference method (FDM). Moreover, the convergence analysis and an upper bound of the error for the derived formula are given. Numerical solutions of ADE are presented and the results are compared with the exact solution.展开更多
With more applications of seismic exploration in metal ore exploration,forward modelling of seismic wave has become more important in metal ore. Finite difference method and pseudo-spectral method are two important me...With more applications of seismic exploration in metal ore exploration,forward modelling of seismic wave has become more important in metal ore. Finite difference method and pseudo-spectral method are two important methods of wave-field simulation. Results of previous studies show that both methods have distinct advantages and disadvantages: Finite difference method has high precision but its dispersion is serious; pseudospectral method considers both computational efficiency and precision but has less precision than finite-difference. The authors consider the complex structural characteristics of the metal ore,furthermore add random media in order to simulate the complex effects produced by metal ore for wave field. First,the study introduced the theories of random media and two forward modelling methods. Second,it compared the simulation results of two methods on fault model. Then the authors established a complex metal ore model,added random media and compared computational efficiency and precision. As a result,it is found that finite difference method is better than pseudo-spectral method in precision and boundary treatment,but the computational efficiency of pseudospectral method is slightly higher than the finite difference method.展开更多
In this study, the numerical solution for the Modified Equal Width Wave (MEW) equation is presented using Fourier spectral method that use to discretize the space variable and Leap-frog method scheme for time dependen...In this study, the numerical solution for the Modified Equal Width Wave (MEW) equation is presented using Fourier spectral method that use to discretize the space variable and Leap-frog method scheme for time dependence. Test problems including the single soliton wave motion, interaction of two solitary waves and interaction of three solitary waves will use to validate the proposed method. The three invariants of the motion are evaluated to determine the conservation properties of the generated scheme. Finally, a Maxwellian initial condition pulse is then studied. The L<sub>2</sub> and L<sub>∞</sub> error norms are computed to study the accuracy and the simplicity of the presented method.展开更多
This paper presents a numerical scheme for space fractional diffusion equations (SFDEs) based on pseudo-spectral method. In this approach, using the Guass-Lobatto nodes, the unknown function is approximated by orthogo...This paper presents a numerical scheme for space fractional diffusion equations (SFDEs) based on pseudo-spectral method. In this approach, using the Guass-Lobatto nodes, the unknown function is approximated by orthogonal polynomials or interpolation polynomials. Then, by using pseudo-spectral method, the SFDE is reduced to a system of ordinary differential equations for time variable t. The high order Runge-Kutta scheme can be used to solve the system. So, a high order numerical scheme is derived. Numerical examples illustrate that the results obtained by this method agree well with the analytical solutions.展开更多
This paper deals with the numerical simulation of incompressible turbulent boundary flow of a flat plate with the pseudo-spectral matrix method. In order to appear more than 10 nodes in the turbulent base-stratum and ...This paper deals with the numerical simulation of incompressible turbulent boundary flow of a flat plate with the pseudo-spectral matrix method. In order to appear more than 10 nodes in the turbulent base-stratum and transition of 43×43 computational grids,a coordinate transformation is put up from physical panel to computational panel. Several zero turbulent models are computed comparatively. The results are credible when comparing with the previous methods.展开更多
In this paper, we apply the Legendre spectral-collocation method to obtain approximate solutions of nonlinear multi-order fractional differential equations (M-FDEs). The fractional derivative is described in the Caput...In this paper, we apply the Legendre spectral-collocation method to obtain approximate solutions of nonlinear multi-order fractional differential equations (M-FDEs). The fractional derivative is described in the Caputo sense. The study is conducted through illustrative example to demonstrate the validity and applicability of the presented method. The results reveal that the proposed method is very effective and simple. Moreover, only a small number of shifted Legendre polynomials are needed to obtain a satisfactory result.展开更多
This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructe...This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructed to find the minimum number of Gaussian quadrature points satisfying the given accuracy.The constructed model is trained by using a large amount of data calculated in the traditional boundary element method and the optimal network architecture is selected.The two-dimensional potential problem of a circular structure is tested and analyzed based on the determined model,and the accuracy of the model is about 90%.Finally,by incorporating the predicted Gaussian quadrature points into the boundary element analysis,we find that the numerical solution and the analytical solution are in good agreement,which verifies the robustness of the proposed method.展开更多
To make the quantitative results of nuclear magnetic resonance(NMR) transverse relaxation(T;) spectrums reflect the type and pore structure of reservoir more directly, an unsupervised clustering method was developed t...To make the quantitative results of nuclear magnetic resonance(NMR) transverse relaxation(T;) spectrums reflect the type and pore structure of reservoir more directly, an unsupervised clustering method was developed to obtain the quantitative pore structure information from the NMR T;spectrums based on the Gaussian mixture model(GMM). Firstly, We conducted the principal component analysis on T;spectrums in order to reduce the dimension data and the dependence of the original variables. Secondly, the dimension-reduced data was fitted using the GMM probability density function, and the model parameters and optimal clustering numbers were obtained according to the expectation-maximization algorithm and the change of the Akaike information criterion. Finally, the T;spectrum features and pore structure types of different clustering groups were analyzed and compared with T;geometric mean and T;arithmetic mean. The effectiveness of the algorithm has been verified by numerical simulation and field NMR logging data. The research shows that the clustering results based on GMM method have good correlations with the shape and distribution of the T;spectrum, pore structure, and petroleum productivity, providing a new means for quantitative identification of pore structure, reservoir grading, and oil and gas productivity evaluation.展开更多
We present a method for derivation of the density matrix of an arbitrary multi-mode continuous variable Gaussian entangled state from its phase space representation.An explicit computer algorithm is given to reconstru...We present a method for derivation of the density matrix of an arbitrary multi-mode continuous variable Gaussian entangled state from its phase space representation.An explicit computer algorithm is given to reconstruct the density matrix from Gaussian covariance matrix and quadrature average values.As an example,we apply our method to the derivation of three-mode symmetric continuous variable entangled state.Our method can be used to analyze the entanglement and correlation in continuous variable quantum network with multi-mode quantum entanglement states.展开更多
The Gaussian weighted trajectory method （GWTM） is a practical implementation of classical S matrix theory （CSMT） in the random phase approximation, CSMT being the first and simplest semi-classical approach of mole...The Gaussian weighted trajectory method （GWTM） is a practical implementation of classical S matrix theory （CSMT） in the random phase approximation, CSMT being the first and simplest semi-classical approach of molecular collisions, developped in the early seventies. Though very close in spirit to the purely classical description, GWTM accounts to some extent for the quantization of the different degrees-of-freedom involved in the processes. While CSMT may give diverging final state distributions, in relation to the rainbow effect of elastic scattering theory, GWTM has never led to such a mathematical catastrophe. The goal of the present note is to explain this finding.展开更多
Aiming at the missile avoidance problem of the unmanned aerial vehicle(UAV)in complex obstacle environments,this work proposes a collision-avoidance method based on receding horizon optimization.The proposed method ge...Aiming at the missile avoidance problem of the unmanned aerial vehicle(UAV)in complex obstacle environments,this work proposes a collision-avoidance method based on receding horizon optimization.The proposed method generated a specific trajectory for the UAV to effectively induce the proportional navigation missile to successfully intercept the obstacle,thereby accomplishing the evasive maneuver.The evasive maneuver was divided into two distinct stages,namely the collision-inducing phase and the fast departure phase.The obstacle potential field-based target selection algorithm was employed to identify the most appropriate target obstacle,while the induced trajectory was determined through a combination of receding horizon optimization and the hp-adaptive pseudo-spectral method.Simulation experiments were carried out under three different types of obstacle environments and one multiobstacle environment,and the simulation results show that the method proposed in this paper greatly improves the success rate of UAV evasive maneuvers,proving the effectiveness of this method.展开更多
Fluid-conveying pipes generally face combined excitations caused by periodic loads and random noises.Gaussian white noise is a common random noise excitation.This study investigates the random vibration response of a ...Fluid-conveying pipes generally face combined excitations caused by periodic loads and random noises.Gaussian white noise is a common random noise excitation.This study investigates the random vibration response of a simply-supported pipe conveying fluid under combined harmonic and Gaussian white noise excitations.According to the generalized Hamilton’s principle,the dynamic model of the pipe conveying fluid under combined harmonic and Gaussian white noise excitations is established.Subsequently,the averaged stochastic differential equations and Fokker–Planck–Kolmogorov(FPK)equations of the pipe conveying fluid subjected to combined excitations are acquired by the modified stochastic averaging method.The effectiveness of the analysis results is verified through the Monte Carlo method.The effects of fluid speed,noise intensity,amplitude of harmonic excitation,and damping factor on the probability density functions of amplitude,displacement,as well as velocity are discussed in detail.The results show that with an increase in fluid speed or noise intensity,the possible greatest amplitude for the fluid-conveying pipe increases,and the possible greatest displacement and velocity also increase.With an increase in the amplitude of harmonic excitation or damping factor,the possible greatest amplitude for the pipe decreases,and the possible greatest displacement and velocity also decrease.展开更多
The numerical simulation of the fluid flow and the flexible rod(s)interaction is more complicated and has lower efficiency due to the high computational cost.In this paper,a semi-resolved model coupling the computatio...The numerical simulation of the fluid flow and the flexible rod(s)interaction is more complicated and has lower efficiency due to the high computational cost.In this paper,a semi-resolved model coupling the computational fluid dynamics and the flexible rod dynamics is proposed using a two-way domain expansion method.The gov-erning equations of the flexible rod dynamics are discretized and solved by the finite element method,and the fluid flow is simulated by the finite volume method.The interaction between fluids and solid rods is modeled by introducing body force terms into the momentum equations.Referred to the traditional semi-resolved numerical model,an anisotropic Gaussian kernel function method is proposed to specify the interactive forces between flu-ids and solid bodies for non-circle rod cross-sections.A benchmark of the flow passing around a single flexible plate with a rectangular cross-section is used to validate the algorithm.Focused on the engineering applications,a test case of a finite patch of cylinders is implemented to validate the accuracy and efficiency of the coupled model.展开更多
In order to enhance the reliability of the moving target detection, an adaptive moving target detection algorithm based on the Gaussian mixture model is proposed. This algorithm employs Gaussian mixture distributions ...In order to enhance the reliability of the moving target detection, an adaptive moving target detection algorithm based on the Gaussian mixture model is proposed. This algorithm employs Gaussian mixture distributions in modeling the background of each pixel. As a result, the number of Gaussian distributions is not fixed but adaptively changes with the change of the pixel value frequency. The pixels of the difference image are divided into two parts according to their values. Then the two parts are separately segmented by the adaptive threshold, and finally the foreground image is obtained. The shadow elimination method based on morphological reconstruction is introduced to improve the performance of foreground image＇s segmentation. Experimental results show that the proposed algorithm can quickly and accurately build the background model and it is more robust in different real scenes.展开更多
文摘Solar-powered aircraft have attracted great attention owing to their potential for longendurance flight and wide application prospects.Due to the particularity of energy system,flight strategy optimization is a significant way to enhance the flight performance for solar-powered aircraft.In this study,a flight strategy optimization model for high-altitude long-endurance solar-powered aircraft was proposed.This model consists of three-dimensional kinematic model,aerodynamic model,energy collection model,energy store model and energy loss model.To solve the nonlinear optimal control problem with process constraints and terminal constraints,Gauss pseudo-spectral method was employed to discretize the state equations and constraint equations.Then a typical mission flying from given initial point to given final point within a time interval was considered.Results indicate that proper changes of the attitude angle contribute to increasing the energy gained by photovoltaic cells.Utilization of gravitational potential energy can partly take the role of battery pack.Integrating these two measures,the optimized flight strategy can improve the final state of charge compared with current constant-altitude constant-velocity strategy.The optimized strategy brings more profits on condition of lower sunlight intensity and shorter daytime.
文摘Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. In this paper, an efficient numerical method for solving the fractional Advection-dispersion equation (ADE) is considered. The fractional derivative is described in the Caputo sense. The method is based on Chebyshev approximations. The properties of Chebyshev polynomials are used to reduce ADE to a system of ordinary differential equations, which are solved using the finite difference method (FDM). Moreover, the convergence analysis and an upper bound of the error for the derived formula are given. Numerical solutions of ADE are presented and the results are compared with the exact solution.
基金Supported by the National"863"Project(No.2014AA06A605)
文摘With more applications of seismic exploration in metal ore exploration,forward modelling of seismic wave has become more important in metal ore. Finite difference method and pseudo-spectral method are two important methods of wave-field simulation. Results of previous studies show that both methods have distinct advantages and disadvantages: Finite difference method has high precision but its dispersion is serious; pseudospectral method considers both computational efficiency and precision but has less precision than finite-difference. The authors consider the complex structural characteristics of the metal ore,furthermore add random media in order to simulate the complex effects produced by metal ore for wave field. First,the study introduced the theories of random media and two forward modelling methods. Second,it compared the simulation results of two methods on fault model. Then the authors established a complex metal ore model,added random media and compared computational efficiency and precision. As a result,it is found that finite difference method is better than pseudo-spectral method in precision and boundary treatment,but the computational efficiency of pseudospectral method is slightly higher than the finite difference method.
文摘In this study, the numerical solution for the Modified Equal Width Wave (MEW) equation is presented using Fourier spectral method that use to discretize the space variable and Leap-frog method scheme for time dependence. Test problems including the single soliton wave motion, interaction of two solitary waves and interaction of three solitary waves will use to validate the proposed method. The three invariants of the motion are evaluated to determine the conservation properties of the generated scheme. Finally, a Maxwellian initial condition pulse is then studied. The L<sub>2</sub> and L<sub>∞</sub> error norms are computed to study the accuracy and the simplicity of the presented method.
文摘This paper presents a numerical scheme for space fractional diffusion equations (SFDEs) based on pseudo-spectral method. In this approach, using the Guass-Lobatto nodes, the unknown function is approximated by orthogonal polynomials or interpolation polynomials. Then, by using pseudo-spectral method, the SFDE is reduced to a system of ordinary differential equations for time variable t. The high order Runge-Kutta scheme can be used to solve the system. So, a high order numerical scheme is derived. Numerical examples illustrate that the results obtained by this method agree well with the analytical solutions.
文摘This paper deals with the numerical simulation of incompressible turbulent boundary flow of a flat plate with the pseudo-spectral matrix method. In order to appear more than 10 nodes in the turbulent base-stratum and transition of 43×43 computational grids,a coordinate transformation is put up from physical panel to computational panel. Several zero turbulent models are computed comparatively. The results are credible when comparing with the previous methods.
文摘In this paper, we apply the Legendre spectral-collocation method to obtain approximate solutions of nonlinear multi-order fractional differential equations (M-FDEs). The fractional derivative is described in the Caputo sense. The study is conducted through illustrative example to demonstrate the validity and applicability of the presented method. The results reveal that the proposed method is very effective and simple. Moreover, only a small number of shifted Legendre polynomials are needed to obtain a satisfactory result.
基金The authors thank the financial support of National Natural Science Foundation of China(NSFC)under Grant(No.11702238).
文摘This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructed to find the minimum number of Gaussian quadrature points satisfying the given accuracy.The constructed model is trained by using a large amount of data calculated in the traditional boundary element method and the optimal network architecture is selected.The two-dimensional potential problem of a circular structure is tested and analyzed based on the determined model,and the accuracy of the model is about 90%.Finally,by incorporating the predicted Gaussian quadrature points into the boundary element analysis,we find that the numerical solution and the analytical solution are in good agreement,which verifies the robustness of the proposed method.
基金Supported by the National Natural Science Foundation of China (42174142)National Science and Technology Major Project (2017ZX05039-002)+2 种基金Operation Fund of China National Petroleum Corporation Logging Key Laboratory (2021DQ20210107-11)Fundamental Research Funds for Central Universities (19CX02006A)Major Science and Technology Project of China National Petroleum Corporation (ZD2019-183-006)。
文摘To make the quantitative results of nuclear magnetic resonance(NMR) transverse relaxation(T;) spectrums reflect the type and pore structure of reservoir more directly, an unsupervised clustering method was developed to obtain the quantitative pore structure information from the NMR T;spectrums based on the Gaussian mixture model(GMM). Firstly, We conducted the principal component analysis on T;spectrums in order to reduce the dimension data and the dependence of the original variables. Secondly, the dimension-reduced data was fitted using the GMM probability density function, and the model parameters and optimal clustering numbers were obtained according to the expectation-maximization algorithm and the change of the Akaike information criterion. Finally, the T;spectrum features and pore structure types of different clustering groups were analyzed and compared with T;geometric mean and T;arithmetic mean. The effectiveness of the algorithm has been verified by numerical simulation and field NMR logging data. The research shows that the clustering results based on GMM method have good correlations with the shape and distribution of the T;spectrum, pore structure, and petroleum productivity, providing a new means for quantitative identification of pore structure, reservoir grading, and oil and gas productivity evaluation.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11574400 and 11204379the Beijing Institute of Technology Research Fund Program for Young Scholarsthe NSFC-ICTP Proposal under Grant No 11981240356
文摘We present a method for derivation of the density matrix of an arbitrary multi-mode continuous variable Gaussian entangled state from its phase space representation.An explicit computer algorithm is given to reconstruct the density matrix from Gaussian covariance matrix and quadrature average values.As an example,we apply our method to the derivation of three-mode symmetric continuous variable entangled state.Our method can be used to analyze the entanglement and correlation in continuous variable quantum network with multi-mode quantum entanglement states.
文摘The Gaussian weighted trajectory method （GWTM） is a practical implementation of classical S matrix theory （CSMT） in the random phase approximation, CSMT being the first and simplest semi-classical approach of molecular collisions, developped in the early seventies. Though very close in spirit to the purely classical description, GWTM accounts to some extent for the quantization of the different degrees-of-freedom involved in the processes. While CSMT may give diverging final state distributions, in relation to the rainbow effect of elastic scattering theory, GWTM has never led to such a mathematical catastrophe. The goal of the present note is to explain this finding.
基金Natural Science Foundation of Heilongjiang Province of China(Grant No.YQ2022F012)the Fundamental Research Funds for the Central Universities(Grant No.HIT.OCEF.2023010)to provide fund for conducting experiments.
文摘Aiming at the missile avoidance problem of the unmanned aerial vehicle(UAV)in complex obstacle environments,this work proposes a collision-avoidance method based on receding horizon optimization.The proposed method generated a specific trajectory for the UAV to effectively induce the proportional navigation missile to successfully intercept the obstacle,thereby accomplishing the evasive maneuver.The evasive maneuver was divided into two distinct stages,namely the collision-inducing phase and the fast departure phase.The obstacle potential field-based target selection algorithm was employed to identify the most appropriate target obstacle,while the induced trajectory was determined through a combination of receding horizon optimization and the hp-adaptive pseudo-spectral method.Simulation experiments were carried out under three different types of obstacle environments and one multiobstacle environment,and the simulation results show that the method proposed in this paper greatly improves the success rate of UAV evasive maneuvers,proving the effectiveness of this method.
基金supported by the National Natural Science Foundation of China(Nos.12272211 and 12072181).
文摘Fluid-conveying pipes generally face combined excitations caused by periodic loads and random noises.Gaussian white noise is a common random noise excitation.This study investigates the random vibration response of a simply-supported pipe conveying fluid under combined harmonic and Gaussian white noise excitations.According to the generalized Hamilton’s principle,the dynamic model of the pipe conveying fluid under combined harmonic and Gaussian white noise excitations is established.Subsequently,the averaged stochastic differential equations and Fokker–Planck–Kolmogorov(FPK)equations of the pipe conveying fluid subjected to combined excitations are acquired by the modified stochastic averaging method.The effectiveness of the analysis results is verified through the Monte Carlo method.The effects of fluid speed,noise intensity,amplitude of harmonic excitation,and damping factor on the probability density functions of amplitude,displacement,as well as velocity are discussed in detail.The results show that with an increase in fluid speed or noise intensity,the possible greatest amplitude for the fluid-conveying pipe increases,and the possible greatest displacement and velocity also increase.With an increase in the amplitude of harmonic excitation or damping factor,the possible greatest amplitude for the pipe decreases,and the possible greatest displacement and velocity also decrease.
基金supported by Shanghai 2021“Science and Technology Innovation Action Plan”:Social Development Science and Technology Research Project(Grant No.21DZ1202703).
文摘The numerical simulation of the fluid flow and the flexible rod(s)interaction is more complicated and has lower efficiency due to the high computational cost.In this paper,a semi-resolved model coupling the computational fluid dynamics and the flexible rod dynamics is proposed using a two-way domain expansion method.The gov-erning equations of the flexible rod dynamics are discretized and solved by the finite element method,and the fluid flow is simulated by the finite volume method.The interaction between fluids and solid rods is modeled by introducing body force terms into the momentum equations.Referred to the traditional semi-resolved numerical model,an anisotropic Gaussian kernel function method is proposed to specify the interactive forces between flu-ids and solid bodies for non-circle rod cross-sections.A benchmark of the flow passing around a single flexible plate with a rectangular cross-section is used to validate the algorithm.Focused on the engineering applications,a test case of a finite patch of cylinders is implemented to validate the accuracy and efficiency of the coupled model.
基金The National Natural Science Foundation of China (No.61172135,61101198)the Aeronautical Foundation of China (No.20115152026)
文摘In order to enhance the reliability of the moving target detection, an adaptive moving target detection algorithm based on the Gaussian mixture model is proposed. This algorithm employs Gaussian mixture distributions in modeling the background of each pixel. As a result, the number of Gaussian distributions is not fixed but adaptively changes with the change of the pixel value frequency. The pixels of the difference image are divided into two parts according to their values. Then the two parts are separately segmented by the adaptive threshold, and finally the foreground image is obtained. The shadow elimination method based on morphological reconstruction is introduced to improve the performance of foreground image＇s segmentation. Experimental results show that the proposed algorithm can quickly and accurately build the background model and it is more robust in different real scenes.
