Navier-Stokes equations are well known in modelling of an incompressible Newtonian fluid, such as air or water. This system of equations is very complex due to the non-linearity term that characterizes it. After the linearization and the discretization parts, we get a descriptor system of index-2 described by a set of differential algebraic equations (DAEs). The two main parts we develop through this paper are focused firstly on constructing an efficient algorithm based on a projection technique onto an extended block Krylov subspace, that appropriately allows us to construct a reduced system of the original DAE system. Secondly, we solve a Linear Quadratic Regulator (LQR) problem based on a Riccati feedback approach. This approach uses numerical solutions of large-scale algebraic Riccati equations. To this end, we use the extended Krylov subspace method that allows us to project the initial large matrix problem onto a low order one that is solved by some direct methods. These numerical solutions are used to obtain a feedback matrix that will be used to stabilize the original system. We conclude by providing some numerical results to confirm the performances of our proposed method compared to other known methods.
翻译:在模拟空气或水等不压缩牛顿液的模型中,纳维-斯托克斯方程式是众所周知的。这个方程式系统由于非线性术语的特点而非常复杂。在线性化和离散部分之后,我们得到一套不同的代数方程式(DAEs)所描述的指数-2的描述系统。我们通过本文件开发的两个主要部分首先侧重于在扩展的Krylov区块子空间的投影技术基础上构建一个高效的算法,这适当使我们能够构建一个原DAE系统的缩小系统。第二,我们根据Riccati反馈方法解决了线性二次调控管(LQR)问题。这个方法使用了大规模代数变相立方程的数值解决方案。我们为此使用扩展的Krylov子空间方法,将最初的大矩阵问题投射到一个以某种直接方法解决的低顺序上。这些数字解决方案用来获取一个反馈矩阵,它将用来稳定原始系统。我们通过一些数字结果来确定我们的拟议方法的性能。我们通过提供一些数字结果来确认其他方法。