Arbitrary High Order Of Accuracy Research Articles

A symplectic analytical singular element for V-notched analyses in elastic and viscoelastic plane problems

The linear elastic and viscoelastic problems with V-shaped notches under in-plane loading are investigated in this paper. We consider two different sets of boundary conditions on the radial edges: Free-Free case and Clamped-Clamped case. A novel singular finite element (SASE) with arbitrary high-order accuracy for the notched problems is proposed. Firstly, via Laplace transform, the original viscoelastic problem is transformed into a corresponding elastic one. Then we construct the SASE by using elastic symplectic eigen solutions with higher order expanding terms. The SASE can depict the characteristics of displacement fields and singular stress fields in the vicinity of the notch vertex having arbitrary opening angle. By taking advantage of the symplectic eigen solutions in the SASE, Mode I and/or Mode II notch stress intensity factors can be determined directly without any post-processing. Fine finite element meshes are not required. Numerical examples are provided to illustrate the validity and accuracy of the present method.

A new fifth-order finite difference well-balanced multi-resolution WENO scheme for solving shallow water equations

In this paper, a new fifth-order finite difference well-balanced multi-resolution weighted essentially non-oscillatory (WENO) scheme is designed to solve for one-dimensional and two-dimensional shallow water equations with or without source terms on structured meshes. We only use the information defined on a hierarchy of nested central spatial stencils without introducing any equivalent multi-resolution representation once again. For the shallow flow problems with smooth or discontinuous bed, we combine with the well-balanced procedure developed by Xing and Shu (2005) for balancing the flux gradients and the source terms, and then the new fifth-order well-balanced multi-resolution WENO scheme (Zhu and Shu, 2018) could satisfy the exact C-property for still stationary solutions, maintain the fifth-order accuracy in smooth regions, and keep essentially non-oscillatory property in non-smooth regions. Compared with the classical well-balanced WENO schemes (Xing and Shu, 2005), the new features of this finite difference well-balanced multi-resolution WENO scheme is its simplicity and hierarchical structure in obtaining higher-order accuracy and its linear weights could be arbitrarily chosen with one requirement that their summation is one. It is the first time that a series of unequal-sized hierarchical central spatial stencils are used in designing finite difference well-balanced WENO scheme for solving shallow water equations. This new well-balanced WENO scheme is simple to construct and can be easily implemented to arbitrary high-order accuracy and in higher dimensions. Some benchmark numerical examples are performed to illustrate the good performances of this new well-balanced WENO scheme which could get high-order accuracy, keep exact C-property, sustain good convergence property, and obtain sharp shock transitions in the whole computational field.

An efficient high-order gas-kinetic scheme (I): Euler equations

In this paper, an efficient high-order gas-kinetic scheme (EHGKS) is proposed to solve the Euler equations for compressible flows. We re-investigate the underlying mechanism of the high-order gas-kinetic scheme (HGKS) and find a new strategy to improve its efficiency. The main idea of the new scheme contains two parts. Firstly, inspired by the state-of-art simplifications on the third-order HGKS, we extend the HGKS to the case of arbitrary high-order accuracy and eliminate its unnecessary high-order dissipation terms. Secondly, instead of computing the derivatives of particle distribution function and their complex moments, we introduce a Lax-Wendroff procedure to compute the high-order derivatives of macroscopic quantities directly. The new scheme takes advantage of both HGKS and the Lax-Wendroff procedure, so that it can be easily extended to the case of arbitrary high-order accuracy with practical significance. Typical numerical tests are carried out by EHGKS, with the third, fifth and seventh order of accuracy. The presence of good resolution of the discontinuities and flow details, together with the optimal CFL numbers, validates the high accuracy and strong robustness of EHGKS. To compare the efficiency, we present the results computed by the EHGKS, the original HGKS and Runge-Kutta-WENO-GKS. This further demonstrates the advantages of EHGKS.

High-order Runge-Kutta discontinuous Galerkin methods with a new type of multi-resolution WENO limiters on triangular meshes

In this paper, high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (WENO) limiters are designed for solving hyperbolic conservation laws on triangular meshes. These multi-resolution WENO limiters are new extensions of the associated multi-resolution WENO finite volume schemes [49,50] which serve as limiters for RKDG methods from structured meshes [47] to triangular meshes. Such new WENO limiters use information of the DG solution essentially only within the troubled cell itself which is identified by a new modified version of the original KXRCF indicator [24], to build a sequence of hierarchical L2 projection polynomials from zeroth degree to the highest degree of the RKDG method. The second-order, third-order, and fourth-order RKDG methods with associated multi-resolution WENO limiters are developed as examples, which could maintain the original order of accuracy in smooth regions and keep essentially non-oscillatory property near strong shocks or contact discontinuities by gradually degrading from the highest order to the first order. The linear weights inside the procedure of the new multi-resolution WENO limiters can be any positive numbers on the condition that their sum equals one. This is the first time that a series of polynomials of different degrees within the troubled cell itself are applied in a WENO fashion to modify the DG solutions in the troubled cell on triangular meshes. These new WENO limiters are very simple to construct, and can be easily implemented to arbitrary high-order accuracy and in higher dimensions on unstructured meshes. Such spatial reconstruction methodology improves the robustness in the simulation on the same compact spatial stencil of the original DG methods on triangular meshes. Extensive one-dimensional (run as two-dimensional problems on triangular meshes) and two-dimensional tests are performed to demonstrate the effectiveness of these RKDG methods with the new multi-resolution WENO limiters.

H(div) finite elements based on nonaffine meshes for 3D mixed formulations of flow problems with arbitrary high order accuracy of the divergence of the flux

SummaryEffects of nonaffine elements on the accuracy of 3D H(div)‐conforming finite elements are studied. Instead of convergence order k+1 for the flux and the divergence of the flux obtained with Raviart‐Thomas or Nédélec spaces with normal traces of degree k, based on affine hexahedra or triangular prisms, reduced orders k for the flux and k−1 for the divergence of the flux may occur for distorted elements. To improve this scenario, a hierarchy of enriched flux approximations is considered, by adding internal shape functions up to a higher degree k+n, n>0, while keeping the original normal traces of degree k. The resulting enriched approximations, using multilinear transformations, keep the original flux accuracy (of order k+1 with affine elements or reduced order k otherwise), but enhanced divergence (of order k+n+1, in the affine case, or k+n−1 otherwise) can be reached. The reduced flux accuracy due to quadrilateral face distortions cannot be corrected by including higher order internal functions. The enriched spaces are applied to the mixed finite element formulation of Darcy's model. The computational cost of matrix assembly increases with n, but the condensed system to be solved has the same dimension and structure as the original scheme.

A Unified Framework for the Solution of Hyperbolic PDE Systems Using High Order Direct Arbitrary-Lagrangian–Eulerian Schemes on Moving Unstructured Meshes with Topology Change

In this work, we review the family of direct Arbitrary-Lagrangian–Eulerian (ALE) finite vlume (FV) and discontinuous Galerkin (DG) schemes on moving meshes that at each time step are rearranged by explicitly allowing topology changes, in order to guarantee a robust mesh evolution even for high shear flow and very long evolution times. Two different techniques are presented: a local nonconforming approach for dealing with sliding lines, and a global regeneration of Voronoi tessellations for treating general unpredicted movements. Corresponding elements at consecutive times are connected in space-time to construct closed space-time control volumes, whose bottom and top faces may be polygons with a different number of nodes, with different neighbors, and even degenerate space-time sliver elements. Our final ALE FV-DG scheme is obtained by integrating, over these arbitrary shaped space-time control volumes, the space-time conservation formulation of the governing hyperbolic PDE system: so, we directly evolve the solution in time avoiding any remapping stage, being conservative and satisfying the GCL by construction. Arbitrary high order of accuracy in space and time is achieved through a fully discrete one-step predictor–corrector ADER approach, also integrated with well balancing techniques to further improve the accuracy and to maintain exactly even at discrete level many physical invariants of the studied system. A large set of different numerical tests has been carried out in order to check the accuracy and the robustness of our methods for both smooth and discontinuous problems, in particular in the case of vortical flows.

A new type of increasingly high-order multi-resolution trigonometric WENO schemes for hyperbolic conservation laws and highly oscillatory problems

In this paper, we investigate designing a new type of high-order finite difference multi-resolution trigonometric weighted essentially non-oscillatory (TWENO) schemes for solving hyperbolic conservation laws and some benchmark highly oscillatory problems. We only use the information defined on a hierarchy of nested central spatial stencils in a trigonometric polynomial reconstruction framework without introducing any equivalent multi-resolution representations. These new finite difference trigonometric WENO schemes use the same large stencils as the classical WENO schemes (Jiang and Wu, 1996; Shu, 2009), could obtain the optimal order of accuracy in smooth regions, and simultaneously suppress spurious oscillations near strong discontinuities. The linear weights of such multi-resolution trigonometric WENO schemes can be any positive numbers on condition that their summation is one. This is the first time that a series of unequal-sized hierarchical central spatial stencils are used in designing high-order finite difference trigonometric WENO schemes. These new trigonometric WENO schemes are simple to construct and can be easily implemented to arbitrary high-order accuracy in multi-dimensions. Some benchmark examples including some highly oscillatory problems are given to demonstrate the robustness and good performance of these new trigonometric WENO schemes.

High-order Runge-Kutta discontinuous Galerkin methods with a new type of multi-resolution WENO limiters

In this paper, a new type of multi-resolution weighted essentially non-oscillatory (WENO) limiters for high-order Runge-Kutta discontinuous Galerkin (RKDG) methods is designed. This type of multi-resolution WENO limiters is an extension of the multi-resolution WENO finite volume and finite difference schemes developed in [43]. Such new limiters use information of the DG solution essentially only within the troubled cell itself, to build a sequence of hierarchical L2 projection polynomials from zeroth degree to the highest degree of the RKDG method. The second-order, third-order, fourth-order, and fifth-order RKDG methods with these multi-resolution WENO limiters have been developed as examples, which could maintain the original order of accuracy in smooth regions and could simultaneously suppress spurious oscillations near strong discontinuities. The linear weights of such new multi-resolution WENO limiters can be any positive numbers on the condition that their sum equals one. This is the first time that a series of polynomials of different degrees within the troubled cell itself are applied in a WENO fashion to modify the DG solutions in the troubled cell. These new WENO limiters are very simple to construct, and can be easily implemented to arbitrary high-order accuracy and in higher dimensions. Such spatial reconstruction methodology improves the robustness in the numerical simulation on the same compact spatial stencil of the original DG methods. Benchmark examples are given to demonstrate the good performance of these RKDG methods with the associated multi-resolution WENO limiters.

An arbitrary high order series-based time-marching method with discontinuous Galerkin method for 2D convection problem

This paper presents an arbitrary high-order series-based time-marching (SBTM) method to match with the high-order discontinuous Galerkin (DG) method in spatial discretization. And the detail of this hybrid method, SBTM-DG, is presented by solving a 2D convection problem where the exact solution of the semi-discrete system derived by DG method can be obtained by SBTM method. The SBTM-DG method is explicit single-step scheme of arbitrary high order of accuracy in time and with all advantages of DG method in space. The high order of accuracy in temporal direction is obtained by multiplications between small-scale matrixes and vectors in each cells with dimension K=(k+1)(k+2)/2, where k is the degree of basic function of DG method, which makes the SBTM-DG method memory-saving and efficient. With the help of projection, the error estimation between the exact solution and the full-discrete solution is O(h)k+1+O(Δt)S, where S is the degree of polynomial of time variable. By energy method we can prove the stability. By matrix method, we obtain the stability condition where time step Δt can be larger than mesh size h. Numerical examples confirm the analysis and show the good performance of SBTM-DG method. Furthermore, the convergence orders are around 20 about k, which makes SBTM-DG method significant.

Robust high‐order semi‐implicit backward difference formulas for Navier‐Stokes equations

SummaryTaylor‐Hood finite elements provide a robust numerical discretization of Navier‐Stokes equations (NSEs) with arbitrary high order of accuracy in space. To match the accuracy of the lowest degree Taylor‐Hood element, we propose a very efficient time‐stepping methods for unsteady flows, which are based on high‐order semi‐implicit backward difference formulas (SBDF) and the inclusion of grad‐div term in the NSE. To mitigate the impact on the numerical accuracy (in time) of the extrapolation of the nonlinear term in SBDF, several variants of nonlinear extrapolation formulas are investigated. The first approach is based on an extrapolation of the nonlinear advection term itself. The second formula uses the extrapolation of the velocity prior to the evaluation of the nonlinear advection term as a whole. The third variant is constructed such that it produces similar error on both velocity and pressure to that with fully implicit backward difference formulas methods at a given order of accuracy. This can be achieved by fixing one‐order higher than usually done in the extrapolation formula for the nonlinear advection term, while keeping the same extrapolation formula for the time derivative. The resulting truncation errors (in time) between these formulas are identified using Taylor expansions. These truncation error formulas are shown to properly represent the error seen in numerical tests using a 2D manufactured solution. Lastly, we show the robustness of the proposed semi‐implicit methods by solving test cases with high Reynolds numbers using one of the nonlinear extrapolation formulas, namely, the 2D flow past circular cylinder at Re=300 and Re = 1000 and the 2D lid‐driven cavity at Re = 50 000 and Re = 100 000. Our numerical solutions are found to be in a good agreement with those obtained in the literature, both qualitatively and quantitatively.

Default mode network is mechanism for semantic color-emotional interdomain integration

Progressing toward more accurate and more efficient numerical codes for the simulation of transport and turbulence in the edge plasma of tokamaks, we propose in this work a new hybrid discontinuous Galerkin solver. Based on 2D advection–diffusion conservation equations for the ion density and the particle flux in the direction parallel to the magnetic field, the code simulates plasma transport in the poloidal section of tokamaks, including the open field lines of the Scrape-off Layer (SOL) and the closed field lines of the core region. The spatial discretization is based on a high-order hybrid DG scheme on unstructured meshes, which provides an arbitrary high-order accuracy while reducing considerably the number of coupled degrees of freedom with a local condensation process. A discontinuity sensor is employed to identify critical elements and regularize the solution with the introduction of artificial diffusion. Based on a finite-element discretization, not constrained by a flux-aligned mesh, the code is able to describe plasma facing components of any complex shape using Bohm boundary conditions and to simulate the plasma in versatile magnetic equilibria, possibly extended up to the center. Numerical tests using a manufactured solution show appropriate convergence orders when varying independently the number of elements or the degree of interpolation. Validation is performed by benchmarking the code with the well-referenced edge transport code SOLEDGE2D (Bufferand et al., 2013, 2015 [1], [2]) in the WEST geometry. Final numerical experiments show the capacity of the code to deal with low-diffusion solutions.

A new type of multi-resolution WENO schemes with increasingly higher order of accuracy

In this paper, a new type of high-order finite difference and finite volume multi-resolution weighted essentially non-oscillatory (WENO) schemes is presented for solving hyperbolic conservation laws. We only use the information defined on a hierarchy of nested central spatial stencils and do not introduce any equivalent multi-resolution representation. These new WENO schemes use the same large stencils as the classical WENO schemes in [25,45], could obtain the optimal order of accuracy in smooth regions, and could simultaneously suppress spurious oscillations near discontinuities. The linear weights of such WENO schemes can be any positive numbers on the condition that their sum equals one. This is the first time that a series of unequal-sized hierarchical central spatial stencils are used in designing high-order finite difference and finite volume WENO schemes. These new WENO schemes are simple to construct and can be easily implemented to arbitrary high order of accuracy and in higher dimensions. Benchmark examples are given to demonstrate the robustness and good performance of these new WENO schemes.

A hybrid discontinuous Galerkin method for tokamak edge plasma simulations in global realistic geometry

Progressing toward more accurate and more efficient numerical codes for the simulation of transport and turbulence in the edge plasma of tokamaks, we propose in this work a new hybrid discontinuous Galerkin solver. Based on 2D advection–diffusion conservation equations for the ion density and the particle flux in the direction parallel to the magnetic field, the code simulates plasma transport in the poloidal section of tokamaks, including the open field lines of the Scrape-off Layer (SOL) and the closed field lines of the core region. The spatial discretization is based on a high-order hybrid DG scheme on unstructured meshes, which provides an arbitrary high-order accuracy while reducing considerably the number of coupled degrees of freedom with a local condensation process. A discontinuity sensor is employed to identify critical elements and regularize the solution with the introduction of artificial diffusion. Based on a finite-element discretization, not constrained by a flux-aligned mesh, the code is able to describe plasma facing components of any complex shape using Bohm boundary conditions and to simulate the plasma in versatile magnetic equilibria, possibly extended up to the center. Numerical tests using a manufactured solution show appropriate convergence orders when varying independently the number of elements or the degree of interpolation. Validation is performed by benchmarking the code with the well-referenced edge transport code SOLEDGE2D (Bufferand et al., 2013, 2015 [1,2]) in the WEST geometry. Final numerical experiments show the capacity of the code to deal with low-diffusion solutions.

Arbitrary high order accurate space–time discontinuous Galerkin finite element schemes on staggered unstructured meshes for linear elasticity

In this paper we propose a new high order accurate space–time discontinuous Galerkin (DG) finite element scheme for the solution of the linear elastic wave equations in first order velocity-stress formulation in two and three-space dimensions on staggered unstructured triangular and tetrahedral meshes. The method reaches arbitrary high order of accuracy in both space and time via the use of space–time basis and test functions. Within the staggered mesh formulation, we define the discrete velocity field in the control volumes of a primary mesh, while the discrete stress tensor is defined on a face-based staggered dual mesh. The space–time DG formulation leads to an implicit scheme that requires the solution of a linear system for the unknown degrees of freedom at the new time level. The number of unknowns is reduced at the aid of the Schur complement, so that in the end only a linear system for the degrees of freedom of the velocity field needs to be solved, rather than a system that involves both stress and velocity. Thanks to the use of a spatially staggered mesh, the stencil of the final velocity system involves only the element and its direct neighbors and the linear system can be efficiently solved via matrix-free iterative methods. Despite the necessity to solve a linear system, the numerical scheme is still computationally efficient. The chosen discretization and the linear nature of the governing PDE system lead to an unconditionally stable scheme, which allows large time steps even for low quality meshes that contain so-called sliver elements. The fully discrete staggered space–time DG method is proven to be energy stable for any order of accuracy, for any mesh and for any time step size. For the particular case of a simple Crank–Nicolson time discretization and homogeneous material, the final velocity system can be proven to be symmetric and positive definite and in this case the scheme is also exactly energy preserving. The new scheme is applied to several test problems in two and three space dimensions, providing also a comparison with high order explicit ADER-DG schemes.

High-order approximation of chromatographic models using a nodal discontinuous Galerkin approach

A nodal high-order discontinuous Galerkin finite element (DG-FE) method is presented to solve the equilibrium-dispersive model of chromatography with arbitrary high-order accuracy in space. The method can be considered a high-order extension to the total variation diminishing (TVD) framework used by Javeed et al. (2011a,b, 2013) with an efficient quadrature-free implementation. The framework is used to simulate linear and non-linear multicomponent chromatographic systems. The results confirm arbitrary high-order accuracy and demonstrate the potential for accuracy and speed-up gains obtainable by switching from low-order methods to high-order methods. The results reproduce an analytical solution and are in excellent agreement with numerical reference solutions already published in the literature.

HFVS: An arbitrary high order approach based on flux vector splitting

In this paper, a new scheme of arbitrary high order accuracy in both space and time is proposed to solve hyperbolic conservative laws. The basic idea in the construction is that, based on the idea of the flux vector splitting (FVS), we split all the spatial and time derivatives in the Taylor expansion of the numerical flux into two parts: one part with positive eigenvalues, another with negative eigenvalues. According to a Lax–Wendroff procedure, all the time derivatives are then replaced by spatial derivatives, which are evaluated by using WENO reconstruction polynomials. One of the most significant advantages of the current scheme is very easy to implement. In addition, it is found that the higher spatial and time derivatives produced in the construction of the numerical flux can be regarded as a building block, in the sense that they can be coupled with any extact/approximate Riemann solvers to extend a first-order scheme to very high order accuracy in both space and time. Numerous numerical tests for linear and nonlinear hyperbolic conservative laws are carried out, and the numerical results demonstrate that the proposed scheme is robust and can be of high order accuracy in both space and time.

Finite difference elastic wave modeling with an irregular free surface using ADER scheme

In numerical modeling of seismic wave propagation in the earth, we encounter two important issues: the free surface and the topography of the surface (i.e. irregularities). In this study, we develop a 2D finite difference solver for the elastic wave equation that combines a 4th- order ADER scheme (Arbitrary high-order accuracy using DERivatives), which is widely used in aeroacoustics, with the characteristic variable method at the free surface boundary. The idea is to treat the free surface boundary explicitly by using ghost values of the solution for points beyond the free surface to impose the physical boundary condition. The method is based on the velocity-stress formulation. The ultimate goal is to develop a numerical solver for the elastic wave equation that is stable, accurate and computationally efficient. The solver treats smooth arbitrary-shaped boundaries as simple plane boundaries. The computational cost added by treating the topography is negligible compared to flat free surface because only a small number of grid points near the boundary need to be computed. In the presence of topography, using 10 grid points per shortest shear-wavelength, the solver yields accurate results. Benchmark numerical tests using several complex models that are solved by our method and other independent accurate methods show an excellent agreement, confirming the validity of the method for modeling elastic waves with an irregular free surface.

High order operator splitting methods based on an integral deferred correction framework

Integral deferred correction (IDC) methods have been shown to be an efficient way to achieve arbitrary high order accuracy and possess good stability properties. In this paper, we construct high order operator splitting schemes using the IDC procedure to solve initial value problems (IVPs). We present analysis to show that the IDC methods can correct for both the splitting and numerical errors, lifting the order of accuracy by r with each correction, where r is the order of accuracy of the method used to solve the correction equation. We further apply this framework to solve partial differential equations (PDEs). Numerical examples in two dimensions of linear and nonlinear initial–boundary value problems are presented to demonstrate the performance of the proposed IDC approach.

A posteriori subcell limiting of the discontinuous Galerkin finite element method for hyperbolic conservation laws

The purpose of this work is to propose a novel a posteriori finite volume subcell limiter technique for the Discontinuous Galerkin finite element method for nonlinear systems of hyperbolic conservation laws in multiple space dimensions that works well for arbitrary high order of accuracy in space and time and that does not destroy the natural subcell resolution properties of the DG method. High order time discretization is achieved via a one-step ADER approach that uses a local space–time discontinuous Galerkin predictor method to evolve the data locally in time within each cell.Our new limiting strategy is based on the so-called MOOD paradigm, which a posteriori verifies the validity of a discrete candidate solution against physical and numerical detection criteria after each time step. Here, we employ a relaxed discrete maximum principle in the sense of piecewise polynomials and the positivity of the numerical solution as detection criteria. Within the DG scheme on the main grid, the discrete solution is represented by piecewise polynomials of degree N. For those troubled cells that need limiting, our new limiter approach recomputes the discrete solution by scattering the DG polynomials at the previous time step onto a set of Ns=2N+1 finite volume subcells per space dimension. A robust but accurate ADER-WENO finite volume scheme then updates the subcell averages of the conservative variables within the detected troubled cells. The recomputed subcell averages are subsequently gathered back into high order cell-centered DG polynomials on the main grid via a subgrid reconstruction operator. The choice of Ns=2N+1 subcells is optimal since it allows to match the maximum admissible time step of the finite volume scheme on the subgrid with the maximum admissible time step of the DG scheme on the main grid, minimizing at the same time also the local truncation error of the subcell finite volume scheme. It furthermore provides an excellent subcell resolution of discontinuities.Our new approach is therefore radically different from classical DG limiters, where the limiter is using TVB or (H)WENO reconstruction based on the discrete solution of the DG scheme on the main grid at the new time level. In our case, the discrete solution is recomputed within the troubled cells from the old time level using a different and more robust numerical scheme on a subgrid level.We illustrate the performance of the new a posteriori subcell ADER-WENO finite volume limiter approach for very high order DG methods via the simulation of numerous test cases run on Cartesian grids in two and three space dimensions, using DG schemes of up to tenth order of accuracy in space and time (N=9). The method is also able to run on massively parallel large scale supercomputing infrastructure, which is shown via one 3D test problem that uses 10 billion space–time degrees of freedom per time step.

Runge–Kutta Discontinuous Galerkin Method for Traffic Flow Model on Networks

We propose a bound-preserving Runge---Kutta (RK) discontinuous Galerkin (DG) method as an efficient, effective and compact numerical approach for numerical simulation of traffic flow problems on networks, with arbitrary high order accuracy. Road networks are modeled by graphs, composed of a finite number of roads that meet at junctions. On each road, a scalar conservation law describes the dynamics, while coupling conditions are specified at junctions to define flow separation or convergence at the points where roads meet. We incorporate such coupling conditions in the RK DG framework, and apply an arbitrary high order bound preserving limiter to the RK DG method to preserve the physical bounds on the network solutions (car density). We showcase the proposed algorithm on several benchmark test cases from the literature, as well as several new challenging examples with rich solution structures. Modeling and simulation of Cauchy problems for traffic flows on networks is notorious for lack of uniqueness or (Lipschitz) continuous dependence. The discontinuous Galerkin method proposed here deals elegantly with these problems, and is perhaps the only realistic and efficient high-order method for network problems.