Waves at a fluid-solid interface: Explicit versus implicit formulation of boundary conditions using a discontinuous Galerkin method.

The analysis of the pulsed laser generation of interfacial waves propagating along a plane solid-fluid interface and their interaction with a defect (notch) on the interface are presented. The solid-fluid interface motion due to laser-generated interfacial waves including leaky Rayleigh and Scholte wave is simulated by finite element method for two configurations, one is “hard solid-fluid” interface, like steel-water, the other one is “soft solid-fluid” interface, like Plexiglas-water. And then their interaction with a notch-like defect in the interface has also been studied. Two series defects are introduced into the interface. One is with same width but different depth, and the other one with the same depth but different width. For hard and soft solid-fluid configurations, leaky Rayleigh and Scholte wave is severely affected by the notch, respectively. In contrast, Scholte and leaky Rayleigh wave experiences almost no influence for hard and soft solid-fluid interface. The results may lead to an application for defect identification and nondestructive materials characterization.

A new discontinuous Galerkin finite element method for solving diffusion problems is introduced. Unlike the traditional local discontinuous Galerkin method, the scheme called the direct discontinuous Galerkin (DDG) method is based on the direct weak formulation for solutions of parabolic equations in each computational cell and lets cells communicate via the numerical flux $\widehat{u_x}$ only. We propose a general numerical flux formula for the solution derivative, which is consistent and conservative; and we then introduce a concept of admissibility to identify a class of numerical fluxes so that the nonlinear stability for both one-dimensional (1D) and multidimensional problems are ensured. Furthermore, when applying the DDG scheme with admissible numerical flux to the 1D linear case, kth order accuracy in an energy norm is proven when using kth degree polynomials. The DDG method has the advantage of easier formulation and implementation and efficient computation of the solution. A series of numerical examples are presented to demonstrate the high order accuracy of the method. In particular, we study the numerical performance of the scheme with different admissible numerical fluxes.

Split form ALE discontinuous Galerkin methods with applications to under-resolved turbulent low-Mach number flows

A discontinuous Galerkin finite element method for solving the Troesch’s problem

A dissipation-rate reserving DG method for wave catching-up phenomena in a nonlinearly elastic composite bar

Devising discontinuous Galerkin methods for non-linear hyperbolic conservation laws

We further investigate the high order positivity-preserving discontinuous Galerkin (DG) methods for linear hyperbolic and radiative transfer equations developed in Yuan et al. (SIAM J Sci Comput 38:A2987–A3019, 2016). The DG methods in Yuan et al. (2016) can maintain positivity and high order accuracy, but they rely both on the scaling limiter in Zhang and Shu (J Comput Phys 229:8918–8934, 2010) and a rotational limiter, the latter may alter cell averages of the unmodulated DG scheme, thereby affecting conservation. Even though a Lax–Wendroff type theorem is proved in Yuan et al. (2016), guaranteeing convergence to weak solutions with correct shock speed when such rotational limiter is applied, it would still be desirable if a conservative DG method without changing the cell averages can be obtained which has both high order accuracy and positivity-preserving capability. In this paper, we develop and analyze such a DG method for both linear hyperbolic equations and radiative transfer equations. In the one-dimensional case, the method uses traditional DG space $$P^k$$ of piecewise polynomials of degree at most k. A key result is proved that the unmodulated DG solver in this case can maintain positivity of the cell average if the inflow boundary value and the source term are both positive, therefore the positivity-preserving framework in Zhang and Shu (2010) can be used to obtain a high order conservative positivity-preserving DG scheme. Unfortunately, in two-dimensions this is no longer the case. We show that the unmodulated DG solver based either on $$P^k$$ or $$Q^k$$ spaces (piecewise kth degree polynomials or piecewise tensor-product kth degree polynomials) could generate negative cell averages. We augment the DG space with additional functions so that the positivity of cell averages from the unmodulated DG solver can be restored, thereby leading to high order conservative positivity-preserving DG scheme based on these augmented DG spaces following the framework in Zhang and Shu (2010). Computational results are provided to demonstrate the good performance of our DG schemes.

In this paper, a finite element formulation is defined in the framework of the discontinuous Galerkin method. Discontinuous Galerkin (dG) methods are classically used in fluid mechanics, however recently their application in solid mechanics has become more vivid among scientists. Of special interest is their application in elliptic problems with constraints such as incompressibility which leads to volumetric locking phenomenon and also in some structural models of shells, plates and beams with compatibility constraints, which brings about shear locking [1]. While classical standard Galerkin methods must be continuous, dG methods can be applied for discontinuities across element boundaries, where a jump of a value (displacement) can be observed. In the present work, a dG method is applied to a linear elastic bar, where a weak discontinuity is allowed in the bar. (© 2015 Wiley‐VCH Verlag GmbH & Co. KGaA, Weinheim)

The aim of this work is to provide a solver for viscoelastic multi-phase flows within the Bounded Support Spectral Solver (BoSSS) currently under development at the Chair of Fluid Dynamics at the Technical University of Darmstadt. The discretisation in BoSSS consists of a high-order Discontinuous Galerkin (DG) method for single-phase flow and a high-order eXtended Discontinuous Galerkin (XDG) method for the multi-phase purpose. The solver shall be used to investigate numerically the behaviour of viscoelastic droplets. The macroscopic Oldroyd B model which is used in a wide range of applications is chosen as the constitutive model. A detailed derivation of the system of equations including the modeling principles for the Oldroyd B model is presented. A DG discretisation of the system of equations including the Local Discontinuous Galerkin (LDG) method is presented after introducing the field of the DG method. The derivation of appropriate flux functions for the constitutive equations and the extra stress tensor are one of the key derivations of this scientific work. Difficulties arising in the numerical solution of viscoelastic flow problems for higher Weissenberg numbers for different discretisation methods are due to the convection dominated, mixed hyperbolic-elliptic-parabolic nature of the system of equations. Several strategies are presented which overcome these problems and are known from the literature. A key achievement of this scientific work is the application of the LDG method, originally developed for a hyperbolic system of equations for a Newtonian fluid, on the viscoelastic system of equations which renders methods for preserving ellipticity unnecessary. Furthermore, various strategiesn to enhance and to support convergence of the solution of the DG discretised system are presented. These are the Newton method with different approaches determining the Jacobian of the system, a homotopy continuation method based on the Weissenberg number for a better initial guess for the Newton method, and a troubled cell indicator combined with an artificial diffusion approach or an adaptive mesh refinement strategy, respectively. For the completeness of this work the XDG method is presented using a sharp interface approach with a signed distance level-set function as it is implemented in BoSSS. The single-phase solver is combined with these methods and appropriate flux functions for the interface are implemented to enable multi-phase applications for viscoelastic fluid. Several numerical experiments are conducted to verify and to validate the viscoelastic singlephase solver and to show the capability of the viscoelastic multi-phase solver to simulate viscoelastic droplets. Advantages and disadvantages of the implementation and an outlook for future research can be found in the conclusion.

Coupling local discontinuous and continuous galerkin methods for flow problems

Cavities are strong scattering parts for an aircraft. A discontinuous Galerkin (DG) method is presented to efficiently optimize cavity scattering. The DG method is shown to have ability of only remeshing the optimized parts and not remeshing the whole cavity in the computing process of scattering from the whole cavity. An effective and simple preconditioner, called distance sparse preconditioner (DSP), is applied to the DG method, different from the conventional block-diagonal preconditioner (BDP) and sparse approximate inverse preconditioner (SAI). The former one is usually recommended in the DG methods and the latter one is generally used in method of moment (MoM). Numerical results show that the DSP is more efficient than the BDP and SAI in the DG method for the cavity analysis. The effects of cavity opening shape, area, and transition shape are studied for scattering from cavities. Some interesting and useful numerical results are presented for cavity design.

Error analysis of a novel discontinuous Galerkin method for the two-dimensional Poisson's equation

Convective cooling channels are applied in a two-dimensional compressor vane to use the intercooling method to improve the efficiency of Brayton cycle and reduce the temperature of the vane. In this paper, we analyze the effect of coolant to the aerodynamic performance and heat transfer performance of the main stream and the vane. For the case of a two-dimensional compressor vane NACA65-(12A2I8b)10, the vane which has five convective cooling channels has been numerically simulated in different test conditions by discontinuous Galerkin (DG) method. The coolant is supercritical carbon dioxide whose pressure is 10MPa. Conjugate heat transfer method has been used in this paper. The numerical simulation result is similar to the experiment data and has been compared with the result of the vane without cooling channels to prove the effect of cooling channels. Cooling channels have large effect on the distribution of temperature and heat transfer coefficient. In addition, the relationship between Nu and Re on the fluid-solid interface has been analyzed and a suitable empirical equation has been obtained. This work analyzes the effect of intercooling system in the compressor and give several advice on future engineering applications in aero engines and gas turbines.

The discontinuous Galerkin (dG) method was introduced by Reed and Hill [73] in 1973 for steady-state neutron transport as an hyperbolic problem. This was followed by other studies; by Bassi and Rebay [12] for the compressible Navier-Stokes equations, Cockburn and Shu [31] developed the local discontinuous Galerkin (ldG) method for advection-diffusion equations, and Peraire and Persson [69] introduced the compact discontinuous Galerkin (cdG) method. Independent of the dG methods, interior penalty (IP) methods have been developed for elliptic and parabolic problems by Douglas and Dupont [40] and Wheeler [95]. Then, in the 1980’s, Arnold et al. [6] proposed a unified classification and analysis of various kinds of dG methods. Later on, the dG methods were developed for elliptic problems [8, 24, 77] and for problems with advection [7, 14, 51, 56].

Discontinuous Galerkin (DG) methods support elements of various types, nonmatching grid and varying polynomial order in each element. In DG methods continuity at element interfaces is weakly enforced with the addition of proper penalty terms on the variational formulation commonly referred to as numerical fluxes. An interior penalty approach to derive a DG method for solving the two first order Maxwell's equations in the time domain is presented. The proposed method is explicit and conditionally stable. In addition, a local time-stepping strategy is applied to increase efficiency and reduce the computational time.