This paper addresses the synchronization of follower agents’state vectors with that of a leader in high-order nonlinear multi-agent systems.The proposed low-complexity control scheme employs high-gain observers to es...This paper addresses the synchronization of follower agents’state vectors with that of a leader in high-order nonlinear multi-agent systems.The proposed low-complexity control scheme employs high-gain observers to estimate higher-order synchronization errors,enabling the controller to rely solely on relative output measurements.This approach significantly reduces the dependence on full-state information,which is often infeasible or costly in practical engineering applications.An output feedback control strategy is developed to overcome these limitations while ensuring robust and effective synchronization.Simulation results are provided to demonstrate the effectiveness of the proposed approach and validate the theoretical findings.展开更多
This paper is dedicated to solving the problem of adaptive fuzzy fault-tolerant tracking control for a class of time-varying high-order uncertain nonlinear systems.The motivation comes from how to construct a compact ...This paper is dedicated to solving the problem of adaptive fuzzy fault-tolerant tracking control for a class of time-varying high-order uncertain nonlinear systems.The motivation comes from how to construct a compact set large enough in which the approximation of any unknown continuous function by a fuzzy logic system(FLS)is effective while compensating sensor/actuator faults and external disturbances.The difficulty is to verify the boundedness of closed-loop signals on the constructed compact set and to reduce the number of the variables of the fuzzy membership functions as many as possible.By a new lemma,linear/nonlinear terms are introduced in adaptive laws to dominate unknown residual terms.With adding a power integrator method,a unified fault-tolerant controller is designed to drive the tracking error to converge to a small compact set of the origin within a fixed time,regardless of whether the system suffers from faults and disturbances.Superior to the existing results,in the presence of time-varying factors the scheme of this paper clarifies the logical relationship between the compactness of the approximation and the boundedness of the state variables.Finally,the application of control strategy is demonstrated by numerical/practical examples.展开更多
In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is de...In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is developed using the trigonometric scheme,which is based on zero,first,and second moments,and the direct discontinuous Galerkin(DDG)flux is used to discretize the diffusion term.Moreover,the DDG method directly applies the weak form of the parabolic equation to each computational cell,which can better capture the characteristics of the solution,especially the discontinuous solution.Meanwhile,the third-order TVD-Runge-Kutta method is applied for temporal discretization.Finally,the effectiveness and stability of the method constructed in this paper are evaluated through numerical tests.展开更多
An efficient and accurate scalar auxiliary variable(SAV)scheme for numerically solving nonlinear parabolic integro-differential equation(PIDE)is developed in this paper.The original equation is first transformed into ...An efficient and accurate scalar auxiliary variable(SAV)scheme for numerically solving nonlinear parabolic integro-differential equation(PIDE)is developed in this paper.The original equation is first transformed into an equivalent system,and the k-order backward differentiation formula(BDF k)and central difference formula are used to discretize the temporal and spatial derivatives,respectively.Different from the traditional discrete method that adopts full implicit or full explicit for the nonlinear integral terms,the proposed scheme is based on the SAV idea and can be treated semi-implicitly,taking into account both accuracy and effectiveness.Numerical results are presented to demonstrate the high-order convergence(up to fourth-order)of the developed schemes and it is computationally efficient in long-time computations.展开更多
In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws...In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.展开更多
To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS),optimized flux difference schemes are proposed.The disadvantages in previous optimization routines,i.e.,reducing formal...To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS),optimized flux difference schemes are proposed.The disadvantages in previous optimization routines,i.e.,reducing formal orders,or extending stencil widths,are avoided in the new optimized schemes by utilizing fluxes from both cell-edges and cell-nodes.Optimizations are implemented with Fourier analysis for linear schemes and the approximate dispersion relation(ADR)for nonlinear schemes.Classical difference schemes are restored near discontinuities to suppress numerical oscillations with use of a shock sensor based on smoothness indicators.The results of several benchmark numerical tests indicate that the new optimized difference schemes outperform the classical schemes,in terms of accuracy and resolution for smooth wave and vortex,especially for long-time simulations.Using optimized schemes increases the total CPU time by less than 4%.展开更多
In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.T...In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.The unconditional stability of the LDG scheme is proved,and an a priori error estimate with O(h^(k+1)+At^(2))is derived,where k≥0 denotes the index of the basis function.Extensive numerical results with Q^(k)(k=0,1,2,3)elements are provided to confirm our theoretical results,which also show that the second-order convergence rate in time is not impacted by the changed parameter θ.展开更多
The Spalart-Allmaras （S-A） turbulence model, the shear-stress transport （SST） turbulence model and their compressibility corrections are revaluated for hypersonic compression comer flows by using high-order differ...The Spalart-Allmaras （S-A） turbulence model, the shear-stress transport （SST） turbulence model and their compressibility corrections are revaluated for hypersonic compression comer flows by using high-order difference schemes. The compressibility effect of density gradient, pressure dilatation and turbulent Mach number is accounted. In order to reduce confusions between model uncertainties and discretization errors, the formally fifth-order explicit weighted compact nonlinear scheme （WCNS-E-5） is adopted for convection terms, and a fourth-order staggered central difference scheme is applied for viscous terms. The 15° and 34° compression comers at Mach number 9.22 are investigated. Numerical results show that the original SST model is superior to the original S-A model in the resolution of separated regions and predictions of wall pressures and wall heat-flux rates. The capability of the S-A model can be largely improved by blending Catris＇ and Shur＇s compressibility corrections. Among the three corrections of the SST model listed in the present paper, Catris＇ modification brings the best results. However, the dissipation and pressure dilatation corrections result in much larger separated regions than that of the experiment, and are much worse than the original SST model as well as the other two corrections. The correction of turbulent Mach number makes the separated region slightly smaller than that of the original SST model. Some results of low-order schemes are also presented. When compared to the results of the high-order schemes, the separated regions are smaller, and the peak wall pressures and peak heat-flux rates are lower in the region of the reattachment points.展开更多
Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor.A new high-order scheme(we call it WENO-K)proposed by our research group is employed to minimize numerical d...Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor.A new high-order scheme(we call it WENO-K)proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors.WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction,and its advantage is to reduce dissipation in smooth regions of flow while preserving high-resolution around discontinuities.Here WENO-K scheme is adopted to reconstruct left and right state values within the Roe Riemann solver updating the inviscid fluxes on a structured dynamic overset grid.To minimize the accuracy loss for high-order reconstruction on artificial boundaries of overset grid,a method of multilayer fringes is proposed to carry out interpolation between background grid and blade grid.Massively parallel computing considering automatic load balance on averagely partitioned overset grid is developed to reduce the wall-clock time of an unsteady simulation.Numerical results for Caradonna-Tung(C-T)rotor in hover at the conditions of subsonic and transonic tip Mach numbers show that the thrust coefficient error for the result of WENO-K scheme is no more than 3%.Compared with WENO-JS scheme,WENO-K scheme achieves about 40%improvement on accuracy of predicting rotor thrust with only 4.1%extra computational cost.More importantly,WENO-K scheme can capture more sophisticated unsteady flow structures and resolve tip vortices to a larger wake age with an increment of about 270°compared to WENO-JS scheme.展开更多
A high-order splitting scheme for the advection-diffusion equation of pollutants is proposed in this paper. The multidimensional advection-diffusion equation is splitted into several one-dimensional equations that are...A high-order splitting scheme for the advection-diffusion equation of pollutants is proposed in this paper. The multidimensional advection-diffusion equation is splitted into several one-dimensional equations that are solved by the scheme. Only three spatial grid points are needed in each direction and the scheme has fourth-order spatial accuracy. Several typically pure advection and advection-diffusion problems are simulated. Numerical results show that the accuracy of the scheme is much higher than that of the classical schemes and the scheme can he efficiently solved with little programming effort.展开更多
A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed. The stability condition is r = △t/△x^2 = △t/△y^2 = △t/△z^2 = △t/△w^2 〈 3/8, and the truncatio...A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed. The stability condition is r = △t/△x^2 = △t/△y^2 = △t/△z^2 = △t/△w^2 〈 3/8, and the truncation error is O（△t^2 ＋ △x^4）.展开更多
Based on the Taylor series method and Li’s spatial differential method, a high-order hybrid Taylor–Li scheme is proposed.The results of a linear advection equation indicate that, using the initial values of the squa...Based on the Taylor series method and Li’s spatial differential method, a high-order hybrid Taylor–Li scheme is proposed.The results of a linear advection equation indicate that, using the initial values of the square-wave type, a result with thirdorder accuracy occurs. However, using initial values associated with the Gaussian function type, a result with very high precision appears. The study demonstrates that, when the order of the time integral is more than three, the corresponding optimal spatial difference order could be higher than six. The results indicate that the reason for why there is no improvement related to an order of spatial difference above six is the use of a time integral scheme that is not high enough. The author also proposes a recursive differential method to improve the Taylor–Li scheme’s computation speed. A more rapid and highprecision program than direct computation of the high-order space differential item is employed, and the computation speed is dramatically boosted. Based on a multiple-precision library, the ultrahigh-order Taylor–Li scheme can be used to solve the advection equation and Burgers’ equation.展开更多
A sampling approximation for a function defined on a bounded interval is proposed by combining the Coiflet-type wavelet expansion and the boundary extension technique. Based on such a wavelet approximation scheme, a G...A sampling approximation for a function defined on a bounded interval is proposed by combining the Coiflet-type wavelet expansion and the boundary extension technique. Based on such a wavelet approximation scheme, a Galerkin procedure is developed for the spatial discretization of the generalized nonlinear Schr6dinger （NLS） equa- tions, and a system of ordinary differential equations for the time dependent unknowns is obtained. Then, the classical fourth-order explicit Runge-Kutta method is used to solve this semi-discretization system. To justify the present method, several widely considered problems are solved as the test examples, and the results demonstrate that the proposed wavelet algorithm has much better accuracy and a faster convergence rate in space than many existing numerical methods.展开更多
In this paper, we consider a class of high-order nonlinear systems with unmodelled dynamics from the viewpoint of maintaining the desired control performance （e,g., asymptotical stability） and reducing the control e...In this paper, we consider a class of high-order nonlinear systems with unmodelled dynamics from the viewpoint of maintaining the desired control performance （e,g., asymptotical stability） and reducing the control effort. By introducing a new reseating transformation, adopting an effective reduced-order observer, and choosing an ingenious Lyapunov function and appropriate design parameters, this paper designs all improved output-feedback controller. The output-feedback controller guarantees the globally asymptotieal stability of the closed-loop system. Subsequently, taking a concrete system for an example, the smaller critical values for gain parameter and resealing transformation parameter are obtained to effectively reduce the control effort.展开更多
High-order schemes based on block-structured adaptive mesh refinement method are prepared to solve computational aeroacoustic （CAA） problems with an aim at improving computational efficiency. A number of numerical i...High-order schemes based on block-structured adaptive mesh refinement method are prepared to solve computational aeroacoustic （CAA） problems with an aim at improving computational efficiency. A number of numerical issues associated with high-order schemes on an adaptively refined mesh, such as stability and accuracy are addressed. Several CAA benchmark problems are used to demonstrate the feasibility and efficiency of the approach.展开更多
The weakly ionized plasma flows in aerospace are commonly simulated by the single-fluid model,which cannot describe certain nonequilibrium phenomena by finite collisions of particles,decreasing the fidelity of the sol...The weakly ionized plasma flows in aerospace are commonly simulated by the single-fluid model,which cannot describe certain nonequilibrium phenomena by finite collisions of particles,decreasing the fidelity of the solution.Based on an alternative formulation of the targeted essentially non-oscillatory(TENO)scheme,a novel high-order numerical scheme is proposed to simulate the two-fluid plasmas problems.The numerical flux is constructed by the TENO interpolation of the solution and its derivatives,instead of being reconstructed from the physical flux.The present scheme is used to solve the two sets of Euler equations coupled with Maxwell's equations.The numerical methods are verified by several classical plasma problems.The results show that compared with the original TENO scheme,the present scheme can suppress the non-physical oscillations and reduce the numerical dissipation.展开更多
In this paper, we investigate the problem of global stabilization for a general class of high-order and non-smoothly stabilizable nonlinear systems with both lower-order and higher-order growth conditions. The designe...In this paper, we investigate the problem of global stabilization for a general class of high-order and non-smoothly stabilizable nonlinear systems with both lower-order and higher-order growth conditions. The designed continuous state feedback controller is recursively constructed to guarantee the global strong stabilization of the closed-loop system.展开更多
A new type of high-order multi-resolution weighted essentially non-oscillatory(WENO)schemes(Zhu and Shu in J Comput Phys,375:659-683,2018)is applied to solve for steady-state problems on structured meshes.Since the cl...A new type of high-order multi-resolution weighted essentially non-oscillatory(WENO)schemes(Zhu and Shu in J Comput Phys,375:659-683,2018)is applied to solve for steady-state problems on structured meshes.Since the classical WENO schemes(Jiang and Shu in J Comput Phys,126:202-228,1996)might suffer from slight post-shock oscillations(which are responsible for the residue to hang at a truncation error level),this new type of high-order finite-difference and finite-volume multi-resolution WENO schemes is applied to control the slight post-shock oscillations and push the residue to settle down to machine zero in steady-state simulations.This new type of multi-resolution WENO schemes uses the same large stencils as that of the same order classical WENO schemes,could obtain fifth-order,seventh-order,and ninth-order in smooth regions,and could gradually degrade to first-order so as to suppress spurious oscillations near strong discontinuities.The linear weights of such new multi-resolution WENO schemes can be any positive numbers on the condition that their sum is one.This is the first time that a series of unequal-sized hierarchical central spatial stencils are used in designing high-order finitedifference and finite-volume WENO schemes for solving steady-state problems.In comparison with the classical fifth-order finite-difference and finite-volume WENO schemes,the residue of these new high-order multi-resolution WENO schemes can converge to a tiny number close to machine zero for some benchmark steady-state problems.展开更多
In this paper,the bipartite consensus problem is studied for a class of uncertain high-order nonlinear multi-agent systems.A signed digraph is presented to describe the collaborative and competitive interactions among...In this paper,the bipartite consensus problem is studied for a class of uncertain high-order nonlinear multi-agent systems.A signed digraph is presented to describe the collaborative and competitive interactions among agents.For each agent with lower triangular structure,a time-varying gain compensator is first designed by relative output information of neighboring agents.Subsequently,a distributed controller with dynamic event-triggered mechanism is proposed to drive the bipartite consensus error to zero.It is worth noting that an internal dynamic variable is introduced in triggering function,which plays an essential role in excluding the Zeno behavior and reducing energy consumption.Furthermore,the dynamic event-triggered control protocol is developed for upper triangular multi-agent systems to realize the bipartite consensus without Zeno behavior.Finally,simulation examples are provided to illustrate the effectiveness of the presented results.展开更多
The problem of finite-time stabilization for uncertain nonlinear systems is investigated.It is proved that a class of high-order nonlinear systems in the lower-triangular form is globally stabilized via non-Lipschitz ...The problem of finite-time stabilization for uncertain nonlinear systems is investigated.It is proved that a class of high-order nonlinear systems in the lower-triangular form is globally stabilized via non-Lipschitz continuous state feedback.By using the finite-time Lyapunov stability theorem and the method of non-smooth feedback design,a recursive design procedure is provided,which guarantees the finite-time stability of the closed-loop system.The simulation results show the effectiveness of the theoretical results.展开更多
文摘This paper addresses the synchronization of follower agents’state vectors with that of a leader in high-order nonlinear multi-agent systems.The proposed low-complexity control scheme employs high-gain observers to estimate higher-order synchronization errors,enabling the controller to rely solely on relative output measurements.This approach significantly reduces the dependence on full-state information,which is often infeasible or costly in practical engineering applications.An output feedback control strategy is developed to overcome these limitations while ensuring robust and effective synchronization.Simulation results are provided to demonstrate the effectiveness of the proposed approach and validate the theoretical findings.
基金supported by National Natural Science Foundation of China[grant number 62173208]Taishan Scholar Project of Shandong Province of China[grant number tsqn202103061]。
文摘This paper is dedicated to solving the problem of adaptive fuzzy fault-tolerant tracking control for a class of time-varying high-order uncertain nonlinear systems.The motivation comes from how to construct a compact set large enough in which the approximation of any unknown continuous function by a fuzzy logic system(FLS)is effective while compensating sensor/actuator faults and external disturbances.The difficulty is to verify the boundedness of closed-loop signals on the constructed compact set and to reduce the number of the variables of the fuzzy membership functions as many as possible.By a new lemma,linear/nonlinear terms are introduced in adaptive laws to dominate unknown residual terms.With adding a power integrator method,a unified fault-tolerant controller is designed to drive the tracking error to converge to a small compact set of the origin within a fixed time,regardless of whether the system suffers from faults and disturbances.Superior to the existing results,in the presence of time-varying factors the scheme of this paper clarifies the logical relationship between the compactness of the approximation and the boundedness of the state variables.Finally,the application of control strategy is demonstrated by numerical/practical examples.
基金The Natural Science Foundation of Xinjiang Uygur Autonomous Region of China“RBF-Hermite difference scheme for the time-fractional kdv-Burgers equation”(2024D01C43)。
文摘In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is developed using the trigonometric scheme,which is based on zero,first,and second moments,and the direct discontinuous Galerkin(DDG)flux is used to discretize the diffusion term.Moreover,the DDG method directly applies the weak form of the parabolic equation to each computational cell,which can better capture the characteristics of the solution,especially the discontinuous solution.Meanwhile,the third-order TVD-Runge-Kutta method is applied for temporal discretization.Finally,the effectiveness and stability of the method constructed in this paper are evaluated through numerical tests.
基金Supported by the National Natural Science Foundation of China(Grant Nos.12001210 and 12261103)the Natural Science Foundation of Henan(Grant No.252300420308)the Yunnan Fundamental Research Projects(Grant No.202301AT070117).
文摘An efficient and accurate scalar auxiliary variable(SAV)scheme for numerically solving nonlinear parabolic integro-differential equation(PIDE)is developed in this paper.The original equation is first transformed into an equivalent system,and the k-order backward differentiation formula(BDF k)and central difference formula are used to discretize the temporal and spatial derivatives,respectively.Different from the traditional discrete method that adopts full implicit or full explicit for the nonlinear integral terms,the proposed scheme is based on the SAV idea and can be treated semi-implicitly,taking into account both accuracy and effectiveness.Numerical results are presented to demonstrate the high-order convergence(up to fourth-order)of the developed schemes and it is computationally efficient in long-time computations.
基金Project supported by the National Natural Science Foundation of China(No.11571366)the Basic Research Foundation of National Numerical Wind Tunnel Project(No.NNW2018-ZT4A08)
文摘In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.
基金Project supported by the National Key Project(No.GJXM92579)the Defense Industrial Technology Development Program(No.C1520110002)the State Administration of Science,Technology and Industry for National Defence,China。
文摘To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS),optimized flux difference schemes are proposed.The disadvantages in previous optimization routines,i.e.,reducing formal orders,or extending stencil widths,are avoided in the new optimized schemes by utilizing fluxes from both cell-edges and cell-nodes.Optimizations are implemented with Fourier analysis for linear schemes and the approximate dispersion relation(ADR)for nonlinear schemes.Classical difference schemes are restored near discontinuities to suppress numerical oscillations with use of a shock sensor based on smoothness indicators.The results of several benchmark numerical tests indicate that the new optimized difference schemes outperform the classical schemes,in terms of accuracy and resolution for smooth wave and vortex,especially for long-time simulations.Using optimized schemes increases the total CPU time by less than 4%.
基金This work is supported by the National Natural Science Foundation of China(11661058,11761053)the Natural Science Foundation of Inner Mongolia(2017MS0107)the Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region(NJYT-17-A07).
文摘In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.The unconditional stability of the LDG scheme is proved,and an a priori error estimate with O(h^(k+1)+At^(2))is derived,where k≥0 denotes the index of the basis function.Extensive numerical results with Q^(k)(k=0,1,2,3)elements are provided to confirm our theoretical results,which also show that the second-order convergence rate in time is not impacted by the changed parameter θ.
基金Foundation items： National Basic Research Program of China （2009CB723801） National Natural Science Foundation of China （11072259）
文摘The Spalart-Allmaras （S-A） turbulence model, the shear-stress transport （SST） turbulence model and their compressibility corrections are revaluated for hypersonic compression comer flows by using high-order difference schemes. The compressibility effect of density gradient, pressure dilatation and turbulent Mach number is accounted. In order to reduce confusions between model uncertainties and discretization errors, the formally fifth-order explicit weighted compact nonlinear scheme （WCNS-E-5） is adopted for convection terms, and a fourth-order staggered central difference scheme is applied for viscous terms. The 15° and 34° compression comers at Mach number 9.22 are investigated. Numerical results show that the original SST model is superior to the original S-A model in the resolution of separated regions and predictions of wall pressures and wall heat-flux rates. The capability of the S-A model can be largely improved by blending Catris＇ and Shur＇s compressibility corrections. Among the three corrections of the SST model listed in the present paper, Catris＇ modification brings the best results. However, the dissipation and pressure dilatation corrections result in much larger separated regions than that of the experiment, and are much worse than the original SST model as well as the other two corrections. The correction of turbulent Mach number makes the separated region slightly smaller than that of the original SST model. Some results of low-order schemes are also presented. When compared to the results of the high-order schemes, the separated regions are smaller, and the peak wall pressures and peak heat-flux rates are lower in the region of the reattachment points.
基金co-supported by the National Natural Science Foundation of China(No.12072285)Shaanxi Science foundation for Distinguished Young Scholars,China(No.2020JC-13)。
文摘Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor.A new high-order scheme(we call it WENO-K)proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors.WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction,and its advantage is to reduce dissipation in smooth regions of flow while preserving high-resolution around discontinuities.Here WENO-K scheme is adopted to reconstruct left and right state values within the Roe Riemann solver updating the inviscid fluxes on a structured dynamic overset grid.To minimize the accuracy loss for high-order reconstruction on artificial boundaries of overset grid,a method of multilayer fringes is proposed to carry out interpolation between background grid and blade grid.Massively parallel computing considering automatic load balance on averagely partitioned overset grid is developed to reduce the wall-clock time of an unsteady simulation.Numerical results for Caradonna-Tung(C-T)rotor in hover at the conditions of subsonic and transonic tip Mach numbers show that the thrust coefficient error for the result of WENO-K scheme is no more than 3%.Compared with WENO-JS scheme,WENO-K scheme achieves about 40%improvement on accuracy of predicting rotor thrust with only 4.1%extra computational cost.More importantly,WENO-K scheme can capture more sophisticated unsteady flow structures and resolve tip vortices to a larger wake age with an increment of about 270°compared to WENO-JS scheme.
文摘A high-order splitting scheme for the advection-diffusion equation of pollutants is proposed in this paper. The multidimensional advection-diffusion equation is splitted into several one-dimensional equations that are solved by the scheme. Only three spatial grid points are needed in each direction and the scheme has fourth-order spatial accuracy. Several typically pure advection and advection-diffusion problems are simulated. Numerical results show that the accuracy of the scheme is much higher than that of the classical schemes and the scheme can he efficiently solved with little programming effort.
基金NSF of the Education Department of Henan Province(20031100010)
文摘A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed. The stability condition is r = △t/△x^2 = △t/△y^2 = △t/△z^2 = △t/△w^2 〈 3/8, and the truncation error is O（△t^2 ＋ △x^4）.
基金supported by the National Natural Sciences Foundation of China(Grant Nos.41375112 and 41530426)the Chinese Academy of Sciences Key Technology Talent Program
文摘Based on the Taylor series method and Li’s spatial differential method, a high-order hybrid Taylor–Li scheme is proposed.The results of a linear advection equation indicate that, using the initial values of the square-wave type, a result with thirdorder accuracy occurs. However, using initial values associated with the Gaussian function type, a result with very high precision appears. The study demonstrates that, when the order of the time integral is more than three, the corresponding optimal spatial difference order could be higher than six. The results indicate that the reason for why there is no improvement related to an order of spatial difference above six is the use of a time integral scheme that is not high enough. The author also proposes a recursive differential method to improve the Taylor–Li scheme’s computation speed. A more rapid and highprecision program than direct computation of the high-order space differential item is employed, and the computation speed is dramatically boosted. Based on a multiple-precision library, the ultrahigh-order Taylor–Li scheme can be used to solve the advection equation and Burgers’ equation.
基金supported by the National Natural Science Foundation of China(Nos.11502103 and11421062)the Open Fund of State Key Laboratory of Structural Analysis for Industrial Equipment of China(No.GZ15115)
文摘A sampling approximation for a function defined on a bounded interval is proposed by combining the Coiflet-type wavelet expansion and the boundary extension technique. Based on such a wavelet approximation scheme, a Galerkin procedure is developed for the spatial discretization of the generalized nonlinear Schr6dinger （NLS） equa- tions, and a system of ordinary differential equations for the time dependent unknowns is obtained. Then, the classical fourth-order explicit Runge-Kutta method is used to solve this semi-discretization system. To justify the present method, several widely considered problems are solved as the test examples, and the results demonstrate that the proposed wavelet algorithm has much better accuracy and a faster convergence rate in space than many existing numerical methods.
基金supported by National Natural Science Founda-tion of China (No. 60774010)Natural Science Foundation of JiangsuProvince, Jiangsu "Six Top Talents" (No. 07-A-020)+1 种基金Program for Fundamental Research of Natural Sciences in Universities of JiangsuProvince (No. 07KJB510114)Natural Science Foundation ofXuzhou Normal University (No. 08XLB20)
文摘In this paper, we consider a class of high-order nonlinear systems with unmodelled dynamics from the viewpoint of maintaining the desired control performance （e,g., asymptotical stability） and reducing the control effort. By introducing a new reseating transformation, adopting an effective reduced-order observer, and choosing an ingenious Lyapunov function and appropriate design parameters, this paper designs all improved output-feedback controller. The output-feedback controller guarantees the globally asymptotieal stability of the closed-loop system. Subsequently, taking a concrete system for an example, the smaller critical values for gain parameter and resealing transformation parameter are obtained to effectively reduce the control effort.
基金supported by the National Natural Science Foundation of China (11150110134)the Science Foundation of Aeronautics of China (20101271004)
文摘High-order schemes based on block-structured adaptive mesh refinement method are prepared to solve computational aeroacoustic （CAA） problems with an aim at improving computational efficiency. A number of numerical issues associated with high-order schemes on an adaptively refined mesh, such as stability and accuracy are addressed. Several CAA benchmark problems are used to demonstrate the feasibility and efficiency of the approach.
基金Project supported by the National Natural Science Foundation of China(Nos.12072246,11972272,11872286)the National Numerical Wind Tunnel Project of China(No.NNW2020ZT3-A23)。
文摘The weakly ionized plasma flows in aerospace are commonly simulated by the single-fluid model,which cannot describe certain nonequilibrium phenomena by finite collisions of particles,decreasing the fidelity of the solution.Based on an alternative formulation of the targeted essentially non-oscillatory(TENO)scheme,a novel high-order numerical scheme is proposed to simulate the two-fluid plasmas problems.The numerical flux is constructed by the TENO interpolation of the solution and its derivatives,instead of being reconstructed from the physical flux.The present scheme is used to solve the two sets of Euler equations coupled with Maxwell's equations.The numerical methods are verified by several classical plasma problems.The results show that compared with the original TENO scheme,the present scheme can suppress the non-physical oscillations and reduce the numerical dissipation.
基金supported by National Natural Science Foundation of China (Nos. 61273125 and 61104222)Specialized Research Fund for the Doctoral Program of Higher Education (No. 20103705110002)+3 种基金Program for the Scientific Research Innovation Team in Colleges and Universities of Shandong ProvinceShandong Provincial Natural Science Foundation of China (No. ZR2012FM018)Natural Science Foundation of Jiangsu Province (No. BK2011205)Natural Science Foundation of Jiangsu Normal University(No. 11XLR08)
文摘In this paper, we investigate the problem of global stabilization for a general class of high-order and non-smoothly stabilizable nonlinear systems with both lower-order and higher-order growth conditions. The designed continuous state feedback controller is recursively constructed to guarantee the global strong stabilization of the closed-loop system.
基金supported by the National Natural Science Foundation of China(Grant No.11872210)supported by the National Science Foundation(Grant No.DMS-1719410)
文摘A new type of high-order multi-resolution weighted essentially non-oscillatory(WENO)schemes(Zhu and Shu in J Comput Phys,375:659-683,2018)is applied to solve for steady-state problems on structured meshes.Since the classical WENO schemes(Jiang and Shu in J Comput Phys,126:202-228,1996)might suffer from slight post-shock oscillations(which are responsible for the residue to hang at a truncation error level),this new type of high-order finite-difference and finite-volume multi-resolution WENO schemes is applied to control the slight post-shock oscillations and push the residue to settle down to machine zero in steady-state simulations.This new type of multi-resolution WENO schemes uses the same large stencils as that of the same order classical WENO schemes,could obtain fifth-order,seventh-order,and ninth-order in smooth regions,and could gradually degrade to first-order so as to suppress spurious oscillations near strong discontinuities.The linear weights of such new multi-resolution WENO schemes can be any positive numbers on the condition that their sum is one.This is the first time that a series of unequal-sized hierarchical central spatial stencils are used in designing high-order finitedifference and finite-volume WENO schemes for solving steady-state problems.In comparison with the classical fifth-order finite-difference and finite-volume WENO schemes,the residue of these new high-order multi-resolution WENO schemes can converge to a tiny number close to machine zero for some benchmark steady-state problems.
基金This work was supported by the National Natural Science Foundation of China(Nos.61973189,62073190)the Research Fund for the Taishan Scholar Project of Shandong Province of China(No.ts20190905)the Natural Science Foundation of Shandong Province of China(No.ZR2020ZD25).
文摘In this paper,the bipartite consensus problem is studied for a class of uncertain high-order nonlinear multi-agent systems.A signed digraph is presented to describe the collaborative and competitive interactions among agents.For each agent with lower triangular structure,a time-varying gain compensator is first designed by relative output information of neighboring agents.Subsequently,a distributed controller with dynamic event-triggered mechanism is proposed to drive the bipartite consensus error to zero.It is worth noting that an internal dynamic variable is introduced in triggering function,which plays an essential role in excluding the Zeno behavior and reducing energy consumption.Furthermore,the dynamic event-triggered control protocol is developed for upper triangular multi-agent systems to realize the bipartite consensus without Zeno behavior.Finally,simulation examples are provided to illustrate the effectiveness of the presented results.
基金Sponsored by the National Natural Science Foundation of China (Grant No. 61174001)
文摘The problem of finite-time stabilization for uncertain nonlinear systems is investigated.It is proved that a class of high-order nonlinear systems in the lower-triangular form is globally stabilized via non-Lipschitz continuous state feedback.By using the finite-time Lyapunov stability theorem and the method of non-smooth feedback design,a recursive design procedure is provided,which guarantees the finite-time stability of the closed-loop system.The simulation results show the effectiveness of the theoretical results.
