The deferred correction(DeC)is an iterative procedure,characterized by increasing the accuracy at each iteration,which can be used to design numerical methods for systems of ODEs.The main advantage of such framework i...The deferred correction(DeC)is an iterative procedure,characterized by increasing the accuracy at each iteration,which can be used to design numerical methods for systems of ODEs.The main advantage of such framework is the automatic way of getting arbitrarily high order methods,which can be put in the Runge-Kutta(RK)form.The drawback is the larger computational cost with respect to the most used RK methods.To reduce such cost,in an explicit setting,we propose an efcient modifcation:we introduce interpolation processes between the DeC iterations,decreasing the computational cost associated to the low order ones.We provide the Butcher tableaux of the new modifed methods and we study their stability,showing that in some cases the computational advantage does not afect the stability.The fexibility of the novel modifcation allows nontrivial applications to PDEs and construction of adaptive methods.The good performances of the introduced methods are broadly tested on several benchmarks both in ODE and PDE contexts.展开更多
In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial deriv...In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.展开更多
A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illus...A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.展开更多
Higher order finite difference weighted essentially non-oscillatory(WENO)schemes have been constructed for conservation laws.For multidimensional problems,they offer a high order accuracy at a fraction of the cost of ...Higher order finite difference weighted essentially non-oscillatory(WENO)schemes have been constructed for conservation laws.For multidimensional problems,they offer a high order accuracy at a fraction of the cost of a finite volume WENO or DG scheme of the comparable accuracy.This makes them quite attractive for several science and engineering applications.But,to the best of our knowledge,such schemes have not been extended to non-linear hyperbolic systems with non-conservative products.In this paper,we perform such an extension which improves the domain of the applicability of such schemes.The extension is carried out by writing the scheme in fluctuation form.We use the HLLI Riemann solver of Dumbser and Balsara(J.Comput.Phys.304:275-319,2016)as a building block for carrying out this extension.Because of the use of an HLL building block,the resulting scheme has a proper supersonic limit.The use of anti-diffusive fluxes ensures that stationary discontinuities can be preserved by the scheme,thus expanding its domain of the applicability.Our new finite difference WENO formulation uses the same WENO reconstruction that was used in classical versions,making it very easy for users to transition over to the present formulation.For conservation laws,the new finite difference WENO is shown to perform as well as the classical version of finite difference WENO,with two major advantages:(i)It can capture jumps in stationary linearly degenerate wave families exactly.(i)It only requires the reconstruction to be applied once.Several examples from hyperbolic PDE systems with non-conservative products are shown which indicate that the scheme works and achieves its design order of the accuracy for smooth multidimensional flows.Stringent Riemann problems and several novel multidimensional problems that are drawn from compressible Baer-Nunziato multiphase flow,multiphase debris flow and twolayer shallow water equations are also shown to document the robustness of the method.For some test problems that require well-balancing we have even been able to apply the scheme without any modification and obtain good results.Many useful PDEs may have stiff relaxation source terms for which the finite difference formulation of WENO is shown to provide some genuine advantages.展开更多
A new modification of the Method of Lines is proposed for the solution of first order partial differential equations. The accuracy of the method is shown with the matrix analysis. The method is applied to a number of ...A new modification of the Method of Lines is proposed for the solution of first order partial differential equations. The accuracy of the method is shown with the matrix analysis. The method is applied to a number of test problems, on uniform grids, to compare the accuracy and computational efficiency with the standard method.展开更多
The leaderless and leader-following finite-time consensus problems for multiagent systems(MASs)described by first-order linear hyperbolic partial differential equations(PDEs)are studied.The Lyapunov theorem and the un...The leaderless and leader-following finite-time consensus problems for multiagent systems(MASs)described by first-order linear hyperbolic partial differential equations(PDEs)are studied.The Lyapunov theorem and the unique solvability result for the first-order linear hyperbolic PDE are used to obtain some sufficient conditions for ensuring the finite-time consensus of the leaderless and leader-following MASs driven by first-order linear hyperbolic PDEs.Finally,two numerical examples are provided to verify the effectiveness of the proposed methods.展开更多
High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of th...High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.展开更多
This paper focuses on linear-quadratic(LQ)optimal control for a class of systems governed by first-order hyperbolic partial differential equations(PDEs).Different from most of the previous works,an approach of discret...This paper focuses on linear-quadratic(LQ)optimal control for a class of systems governed by first-order hyperbolic partial differential equations(PDEs).Different from most of the previous works,an approach of discretization-then-continuousization is proposed in this paper to cope with the infinite-dimensional nature of PDE systems.The contributions of this paper consist of the following aspects:(1)The differential Riccati equations and the solvability condition of the LQ optimal control problems are obtained via the discretization-then-continuousization method.(2)A numerical calculation way of the differential Riccati equations and a practical design way of the optimal controller are proposed.Meanwhile,the relationship between the optimal costate and the optimal state is established by solving a set of forward and backward partial difference equations(FBPDEs).(3)The correctness of the method used in this paper is verified by a complementary continuous method and the comparative analysis with the existing operator results is presented.It is shown that the proposed results not only contain the classic results of the standard LQ control problem of systems governed by ordinary differential equations as a special case,but also support the existing operator results and give a more convenient form of computation.展开更多
Dear Editor,This letter focuses on the distributed cooperative regulation problem for a class of networked re-entrant manufacturing systems(RMSs).The networked system is structured with a three-tier architecture:the p...Dear Editor,This letter focuses on the distributed cooperative regulation problem for a class of networked re-entrant manufacturing systems(RMSs).The networked system is structured with a three-tier architecture:the production line,the manufacturing layer and the workshop layer.The dynamics of re-entrant production lines are governed by hyperbolic partial differential equations(PDEs)based on the law of mass conservation.展开更多
For first-order quasilinear hyperbolic systems with zero eigenvaiues, the author establishes the local exact controllability in a shorter time-period by means of internal controls acting on suitable domains. In partic...For first-order quasilinear hyperbolic systems with zero eigenvaiues, the author establishes the local exact controllability in a shorter time-period by means of internal controls acting on suitable domains. In particular, under certain special but reasonable hypotheses, the local exact controllability can be realized only by internal controls, and the control time can be arbitrarily small.展开更多
In this paper,we are concerned with the numerical solutions for the parabolic and hyperbolic partial differential equations with nonlocal boundary conditions.Thus,we presented a new iterative algorithm based on the Re...In this paper,we are concerned with the numerical solutions for the parabolic and hyperbolic partial differential equations with nonlocal boundary conditions.Thus,we presented a new iterative algorithm based on the Restarted Adomian Decomposition Method(RADM)for solving the two equations of different types involving dissimilar boundary and nonlocal conditions.The algorithm presented transforms the given nonlocal initial boundary value problem to a local Dirichlet one and then employs the RADM for the numerical treatment.Numerical comparisons were made between our proposed method and the Adomian Decomposition Method(ADM)to demonstrate the efficiency and performance of the proposed method.展开更多
文摘The deferred correction(DeC)is an iterative procedure,characterized by increasing the accuracy at each iteration,which can be used to design numerical methods for systems of ODEs.The main advantage of such framework is the automatic way of getting arbitrarily high order methods,which can be put in the Runge-Kutta(RK)form.The drawback is the larger computational cost with respect to the most used RK methods.To reduce such cost,in an explicit setting,we propose an efcient modifcation:we introduce interpolation processes between the DeC iterations,decreasing the computational cost associated to the low order ones.We provide the Butcher tableaux of the new modifed methods and we study their stability,showing that in some cases the computational advantage does not afect the stability.The fexibility of the novel modifcation allows nontrivial applications to PDEs and construction of adaptive methods.The good performances of the introduced methods are broadly tested on several benchmarks both in ODE and PDE contexts.
文摘In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.
文摘A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.
基金support via NSF grants NSF-19-04774,NSF-AST-2009776,NASA-2020-1241NASA grant 80NSSC22K0628.DSB+3 种基金HK acknowledge support from a Vajra award,VJR/2018/00129a travel grant from Notre Dame Internationalsupport via AFOSR grant FA9550-20-1-0055NSF grant DMS-2010107.
文摘Higher order finite difference weighted essentially non-oscillatory(WENO)schemes have been constructed for conservation laws.For multidimensional problems,they offer a high order accuracy at a fraction of the cost of a finite volume WENO or DG scheme of the comparable accuracy.This makes them quite attractive for several science and engineering applications.But,to the best of our knowledge,such schemes have not been extended to non-linear hyperbolic systems with non-conservative products.In this paper,we perform such an extension which improves the domain of the applicability of such schemes.The extension is carried out by writing the scheme in fluctuation form.We use the HLLI Riemann solver of Dumbser and Balsara(J.Comput.Phys.304:275-319,2016)as a building block for carrying out this extension.Because of the use of an HLL building block,the resulting scheme has a proper supersonic limit.The use of anti-diffusive fluxes ensures that stationary discontinuities can be preserved by the scheme,thus expanding its domain of the applicability.Our new finite difference WENO formulation uses the same WENO reconstruction that was used in classical versions,making it very easy for users to transition over to the present formulation.For conservation laws,the new finite difference WENO is shown to perform as well as the classical version of finite difference WENO,with two major advantages:(i)It can capture jumps in stationary linearly degenerate wave families exactly.(i)It only requires the reconstruction to be applied once.Several examples from hyperbolic PDE systems with non-conservative products are shown which indicate that the scheme works and achieves its design order of the accuracy for smooth multidimensional flows.Stringent Riemann problems and several novel multidimensional problems that are drawn from compressible Baer-Nunziato multiphase flow,multiphase debris flow and twolayer shallow water equations are also shown to document the robustness of the method.For some test problems that require well-balancing we have even been able to apply the scheme without any modification and obtain good results.Many useful PDEs may have stiff relaxation source terms for which the finite difference formulation of WENO is shown to provide some genuine advantages.
文摘A new modification of the Method of Lines is proposed for the solution of first order partial differential equations. The accuracy of the method is shown with the matrix analysis. The method is applied to a number of test problems, on uniform grids, to compare the accuracy and computational efficiency with the standard method.
基金the National Natural Science Foundation of China(Nos.11671282 and 12171339)。
文摘The leaderless and leader-following finite-time consensus problems for multiagent systems(MASs)described by first-order linear hyperbolic partial differential equations(PDEs)are studied.The Lyapunov theorem and the unique solvability result for the first-order linear hyperbolic PDE are used to obtain some sufficient conditions for ensuring the finite-time consensus of the leaderless and leader-following MASs driven by first-order linear hyperbolic PDEs.Finally,two numerical examples are provided to verify the effectiveness of the proposed methods.
文摘High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.
基金supported by the National Natural Science Foundation of China under Grant Nos.61821004 and 62250056the Natural Science Foundation of Shandong Province under Grant Nos.ZR2021ZD14 and ZR2021JQ24+1 种基金Science and Technology Project of Qingdao West Coast New Area under Grant Nos.2019-32,2020-20,2020-1-4,High-level Talent Team Project of Qingdao West Coast New Area under Grant No.RCTDJC-2019-05Key Research and Development Program of Shandong Province under Grant No.2020CXGC01208.
文摘This paper focuses on linear-quadratic(LQ)optimal control for a class of systems governed by first-order hyperbolic partial differential equations(PDEs).Different from most of the previous works,an approach of discretization-then-continuousization is proposed in this paper to cope with the infinite-dimensional nature of PDE systems.The contributions of this paper consist of the following aspects:(1)The differential Riccati equations and the solvability condition of the LQ optimal control problems are obtained via the discretization-then-continuousization method.(2)A numerical calculation way of the differential Riccati equations and a practical design way of the optimal controller are proposed.Meanwhile,the relationship between the optimal costate and the optimal state is established by solving a set of forward and backward partial difference equations(FBPDEs).(3)The correctness of the method used in this paper is verified by a complementary continuous method and the comparative analysis with the existing operator results is presented.It is shown that the proposed results not only contain the classic results of the standard LQ control problem of systems governed by ordinary differential equations as a special case,but also support the existing operator results and give a more convenient form of computation.
文摘Dear Editor,This letter focuses on the distributed cooperative regulation problem for a class of networked re-entrant manufacturing systems(RMSs).The networked system is structured with a three-tier architecture:the production line,the manufacturing layer and the workshop layer.The dynamics of re-entrant production lines are governed by hyperbolic partial differential equations(PDEs)based on the law of mass conservation.
文摘For first-order quasilinear hyperbolic systems with zero eigenvaiues, the author establishes the local exact controllability in a shorter time-period by means of internal controls acting on suitable domains. In particular, under certain special but reasonable hypotheses, the local exact controllability can be realized only by internal controls, and the control time can be arbitrarily small.
文摘In this paper,we are concerned with the numerical solutions for the parabolic and hyperbolic partial differential equations with nonlocal boundary conditions.Thus,we presented a new iterative algorithm based on the Restarted Adomian Decomposition Method(RADM)for solving the two equations of different types involving dissimilar boundary and nonlocal conditions.The algorithm presented transforms the given nonlocal initial boundary value problem to a local Dirichlet one and then employs the RADM for the numerical treatment.Numerical comparisons were made between our proposed method and the Adomian Decomposition Method(ADM)to demonstrate the efficiency and performance of the proposed method.
