The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectra...The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.展开更多
The spectral/hp element method combines the geometric flexibility of the classical h-type finite element technique with the desirable numerical properties of spectral methods, employing high-degree piecewise polynomia...The spectral/hp element method combines the geometric flexibility of the classical h-type finite element technique with the desirable numerical properties of spectral methods, employing high-degree piecewise polynomial basis functions on coarse finite element-type meshes. The spatial approximation is based upon orthogonal polynomials, such as Legendre or Chebychev polynomials,modified to accommodate a C~0-continuous expansion. Computationally and theoretically, by increasing the polynomial order p,high-precision solutions and fast convergence can be obtained and, in particular, under certain regularity assumptions an exponential reduction in approximation error between numerical and exact solutions can be achieved. This method has now been applied in many simulation studies of both fundamental and practical engineering flows. This paper briefly describes the formulation of the spectral/hp element method and provides an overview of its application to computational fluid dynamics. In particular, it focuses on the use of the spectral/hp element method in transitional flows and ocean engineering. Finally, some of the major challenges to be overcome in order to use the spectral/hp element method in more complex science and engineering applications are discussed.展开更多
This paper investigates an optimal control problem governed by an elliptic equation with integral control and state constraints.The control problem is approxi-mated by the hp spectral element method with high accuracy...This paper investigates an optimal control problem governed by an elliptic equation with integral control and state constraints.The control problem is approxi-mated by the hp spectral element method with high accuracy and geometricflexibil-ity.Optimality conditions of the continuous and discrete optimal control problems are presented,respectively.The a posteriori error estimates both for the control and state variables are established in detail.In addition,illustrative numerical examples are carried out to demonstrate the accuracy of theoretical results and the validity of the proposed method.展开更多
基金Research was supported in part by NSF grant DMS-0800612Research was supported by Applied Mathematics program of the US DOE Office of Advanced Scientific Computing ResearchThe Pacific Northwest National Laboratory is operated by Battelle for the U.S. Department of Energy under Contract DE-AC05-76RL01830
文摘The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.
文摘The spectral/hp element method combines the geometric flexibility of the classical h-type finite element technique with the desirable numerical properties of spectral methods, employing high-degree piecewise polynomial basis functions on coarse finite element-type meshes. The spatial approximation is based upon orthogonal polynomials, such as Legendre or Chebychev polynomials,modified to accommodate a C~0-continuous expansion. Computationally and theoretically, by increasing the polynomial order p,high-precision solutions and fast convergence can be obtained and, in particular, under certain regularity assumptions an exponential reduction in approximation error between numerical and exact solutions can be achieved. This method has now been applied in many simulation studies of both fundamental and practical engineering flows. This paper briefly describes the formulation of the spectral/hp element method and provides an overview of its application to computational fluid dynamics. In particular, it focuses on the use of the spectral/hp element method in transitional flows and ocean engineering. Finally, some of the major challenges to be overcome in order to use the spectral/hp element method in more complex science and engineering applications are discussed.
基金supported by the State Key Program of National Natural Science Foundation of China(No.11931003)National Natural Science Foundation of China(Nos.41974133 and 11971410)+1 种基金Hunan Provincial Innovation Foundation for Postgraduate(Grant No.CX20190463)Research Innovation Foundation for Postgraduate of Xiangtan University(Grant No.XDCX2019B055).
文摘This paper investigates an optimal control problem governed by an elliptic equation with integral control and state constraints.The control problem is approxi-mated by the hp spectral element method with high accuracy and geometricflexibil-ity.Optimality conditions of the continuous and discrete optimal control problems are presented,respectively.The a posteriori error estimates both for the control and state variables are established in detail.In addition,illustrative numerical examples are carried out to demonstrate the accuracy of theoretical results and the validity of the proposed method.
