scispace - formally typeset
Search or ask a question

Showing papers by "Joaquim Peiró published in 2021"


Journal ArticleDOI
TL;DR: In this paper, an implicit discontinuous Galerkin compressible flow solver is proposed for the spectral/h p element open-source software framework Nektar++, where the integration in time is carried out by a singly diagonally implicit Runge-Kutta method and the nonlinear system arising from the implicit time integration is iteratively solved by the Jacobian-free Newton Krylov (JFNK) method.
Abstract: At high Reynolds numbers the use of explicit in time compressible flow simulations with spectral/ h p element discretization can become significantly limited by time step. To alleviate this limitation we extend the capability of the spectral/ h p element open-source software framework, Nektar++, to include an implicit discontinuous Galerkin compressible flow solver. The integration in time is carried out by a singly diagonally implicit Runge–Kutta method. The non-linear system arising from the implicit time integration is iteratively solved by the Jacobian-free Newton Krylov (JFNK) method. A favorable feature of the JFNK approach is its extensive use of the explicit operators available from the previous explicit in time implementation. The functionalities of different building blocks of the implicit solver are analyzed from the point of view of software design and placed in appropriate hierarchical levels in the C++ libraries. In the detailed implementation, the contributions of different parts of the solver to computational cost, memory consumption and programming complexity are also analyzed. A combination of analytical and numerical methods is adopted to simplify the programming complexity in forming the preconditioning matrix. The solver is verified and tested using cases such as manufactured compressible Poiseuille flow, Taylor–Green vortex, turbulent flow over a circular cylinder at Re = 3900 and shock wave boundary-layer interaction. The results show that the implicit solver can speed-up the simulations while maintaining good simulation accuracy.

15 citations


Journal ArticleDOI
TL;DR: In this article, a smooth particle hydrodynamics (SPH) formulation was proposed to capture the sloshing frequencies of tanks and accurately simulate long-duration violent slo-hing in partially filled tanks of arbitrary shape.

13 citations


Journal ArticleDOI
TL;DR: In this paper, the authors present a successful deployment of high-fidelity large-eddy simulation (LES) technologies based on spectral/$hp$ element methods to industrial flow problems, which are characterized by high Reyn...
Abstract: We present a successful deployment of high-fidelity large-eddy simulation (LES) technologies based on spectral/$hp$ element methods to industrial flow problems, which are characterized by high Reyn...

12 citations


Journal ArticleDOI
26 Aug 2021
TL;DR: The results indicate that the adaptive time-stepping strategy can maintain that the discretization error is dominated by the spatial error and relatively high efficiency is obtained for unsteady and steady, well-resolving and under-resolved simulations.
Abstract: A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations. A proper relation between the spatial, temporal and iterative errors generated within one time step is constructed. With an estimate of temporal and spatial error using an embedded Runge-Kutta scheme and a higher order spatial discretization, an adaptive time-stepping strategy is proposed based on the idea that the time step should be the maximum without obviously influencing the total error of the discretization. The designed adaptive time-stepping strategy is then tested in various types of problems including isentropic vortex convection, steady-state flow past a flat plate, Taylor-Green vortex and turbulent flow over a circular cylinder at $${Re}=3\,900$$ . The results indicate that the adaptive time-stepping strategy can maintain that the discretization error is dominated by the spatial error and relatively high efficiency is obtained for unsteady and steady, well-resolved and under-resolved simulations.

4 citations