Abstract:The high Weissenberg number problem (HWNP) has been regarded as a key challenge in computational rheology. Developing stable and flexible numerical methods for simulating highly elastic fluid flows has become the major research subject of scholars. In this work, we present a new stabilized high-order discontinuous Galerkin (DG) method to simulate viscoelastic flows at high Weissenberg (Wi) numbers. To decouple the velocity and pressure, and circumvent the limitation of the Ladyzhenskaya–Babuška–Brezzi condition, the dual splitting scheme based on the implicit backward differentiation formula and explicit extrapolation scheme is performed for the momentum and mass equations. In addition, integration by parts of the pressure gradient term in projection step and the velocity divergence term in pressure Poisson equation is required to enhance the numerical stability of the method at small time steps. To guarantee the positive definiteness of the conformation tensor and improve the computational stability at high Wi numbers, the square-root-conformation representation (SRCR) approach is applied to reconstruct the constitutive equation. Meanwhile, the second-order Runge–Kutta scheme is performed for temporal discretization and the local Lax–Friedrichs flux is employed for the nonlinear convective term to tackle the convection-dominated problems. The resulting numerical method is free from iterative errors, and does not require any additional stabilization term, which gives an efficient and easy-to-be-implemented solver for viscoelastic flows. The numerical study considers different viscoelastic flow problems, including the Poiseuille flow in a channel, the lid-driven cavity flow, the flow past a fixed cylinder and the flow around an oscillating cylinder. The numerical results reveal the high-order accuracy and robustness of our method at high Wi numbers. Compared with the scheme in the literature, the presented method is more flexible in complex irregular regions and much easier to construct the high-order format. This is the first attempt to apply SRCR scheme in framework of DG method, which provides a new efficient high-order method for simulating complex viscoelastic flows at high Wi numbers.

A Comparative Study of the Discontinuous Galerkin and Continuous SUPG Finite Element Methods for Computation of Viscoelastic Flows

Coupling of Finite Element Method and Discontinuous Galerkin Method to Simulate Viscoelastic Flows

A Coupled Continuous and Discontinuous Finite Element Method for the Incompressible Flows

Galerkin/least-square Finite-Element Methods for Steady Viscoelastic Flows

Investigation on Implicit Discontinuous Galerkin Method in Unsteady Flow

Some continuous and discontinuous Galerkin methods and structure preservation for incompressible flows

A Symmetric Direct Discontinuous Galerkin Method for the Compressible Navier-Stokes Equations

A Unified Study of Continuous and Discontinuous Galerkin Methods

Adaptive Discontinuous Galerkin Methods for Solving an Incompressible Stokes Flow Problem with Slip Boundary Condition of Frictional Type

A Numerical Comparison of the Lax–Wendroff Discontinuous Galerkin Method Based on Different Numerical Fluxes

Study on discontinuous Galerkin method for Navier-Stokes equations

Comparison of the high-order Runge-Kutta discontinuous Galerkin method and gas-kinetic scheme for inviscid compressible flow simulations

Discontinuous Galerkin Methods for a Stationary Navier–Stokes Problem with a Nonlinear Slip Boundary Condition of Friction Type

A Physical-Constraint-Preserving Discontinuous Galerkin Method for Weakly Compressible Two-Phase Flows.

A Reconstructed Discontinuous Galerkin Method for Incompressible Flows on Arbitrary Grids

A Multi-Dimensional High-Order Discontinuous Galerkin Method Based on Gas Kinetic Theory for Viscous Flow Computations

An ultraweak DPG method for viscoelastic fluids

On the dissipation of conforming and discontinuous Galerkin schemes for the incompressible Navier–Stokes equations

A Consistent Projection-Based Supg/Pspg Xfem for Incompressible Two-Phase Flows

Comparing Discontinuous Galerkin Shock-Capturing Techniques Applied to Inviscid Three-Dimensional Hypersonic Flows

A High-Order SRCR-DG Method for Simulating Viscoelastic Flows at High Weissenberg Numbers