High-order discontinuous Galerkin computation of axisymmetric transonic flows in safety relief valves

This paper presents a discontinuous Galerkin (DG) discretization of the compressible RANS and k- ω turbulence model equations for two-dimensional axisymmetric flows. The developed code has been applied to investigate the transonic flow in safety relief valves.This new DG implementation has evolved from the DG method presented in [1]. An "exact" Riemann solver is used to compute the interface numerical inviscid flux while the viscous flux discterization relies on the BRMPS scheme [2,3]. Control of oscillations of high-order solutions around shocks is obtained by means of a shock-capturing technique developed and assessed within the EU ADIGMA project [4].The code has been applied to compute the flow in a spring loaded safety valve at several back pressures and different disk lifts. The predicted device flow capacity and the pressure inside its bonnet have been checked against experimental data. The CFD simulations allow to clarify the complex flow patterns occurring and to explain the measured trends.

[1]  Chi-Wang Shu,et al.  High order residual distribution conservative finite difference WENO schemes for convection-diffusion steady state problems on non-smooth meshes , 2007, J. Comput. Phys..

[2]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[3]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[4]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[5]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[6]  S. Rebay,et al.  GMRES Discontinuous Galerkin Solution of the Compressible Navier-Stokes Equations , 2000 .

[7]  G. Bergeles,et al.  Notes on Numerical Fluid Mechanics and Multidisciplinary Design , 2012 .

[8]  Vincenzo Dossena,et al.  On the influence of back pressure and size on the performance of safety valves , 2002 .

[9]  Norbert Kroll,et al.  ADIGMA: A European Project on the Development of Adaptive Higher Order Variational Methods for Aerospace Applications , 2010 .

[10]  Marco Luciano Savini,et al.  Discontinuous Galerkin solution of the Reynolds-averaged Navier–Stokes and k–ω turbulence model equations , 2005 .

[11]  Jérôme Jaffré,et al.  CONVERGENCE OF THE DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR HYPERBOLIC CONSERVATION LAWS , 1995 .

[12]  Douglas N. Arnold,et al.  Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems , 2001, SIAM J. Numer. Anal..

[13]  F. Menter.  Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

[14]  J. Peraire,et al.  Sub-Cell Shock Capturing for Discontinuous Galerkin Methods , 2006 .

[15]  Alessandro Colombo,et al.  Very High-Order Accurate Discontinuous Galerkin Computation of Transonic Turbulent Flows on Aeronautical Configurations , 2010 .

[16]  Robert E. Harris,et al.  High-Order Adaptive Quadrature-Free Spectral Volume Method on Unstructured Grids , 2008 .

[17]  Rémi Abgrall,et al.  On uniformly high-order accurate residual distribution schemes for advection-diffusion , 2008 .

[18]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[19]  Ralf Hartmann,et al.  Adjoint-based error estimation and adaptive mesh refinement for the RANS and k-ω turbulence model equations , 2011, J. Comput. Phys..