Un método de captura de choques basado en las funciones de forma para Galerkin discontinuo de alto orden

This article presents a high-order Discontinuous Galerkin method for compressible flow problems, in which is very frequent the formation of shocks. The stabilization is introduced by a new basis functions. This base has the flexibility to vary locally (within each element) between continuous polynomial functions space and a space of piecewise polynomial functions. Thus, the proposed method provides a bridge between the standard methods of high-order Discontinuous Galerkin and classical Finite Volume methods, maintaining the locality and compactness of the scheme. The variation of basis functions is automatically set according to the regularity of the solution and the stabilization is introduced by the jump operator, standard in Discontinuous Galerkin methods. Unlike the classical methods of slope limiting, the strategy here presented is very local, robust, and applies to any order of approximation. Moreover, the proposed method does not require adaptive mesh refinement techniques and it can be used with any temporal integration scheme. Several applications of the Euler equations are shown, demonstrating the validity and effectiveness of the method, especially for high orders of approximation.     

HIGH-ORDER NONLINEAR GALERKIN METHODS

The nonlinear Galerkin methods are numerical schemes well adapted to the long-term integration of nonlinear evolution partial differential equations. In this paper, a class of high-order nonlinear Galerkin methods are provided. Moreover, convergence results with high-order spectral accuracy are derived for the schemes introduced.     

Comparison of generalized differential quadrature and Galerkin methods for the analysis of micro-electro-mechanical coupled systems

This paper presents a comprehensive comparison study between the generalized differential quadrature (GDQ) and the well-known global Galerkin method for analysis of pull-in behavior of nonlinear micro-electro-mechanical coupled systems. The nonlinear governing integro-differential equation for double clamped MEMS devices which was derived using variational principle by the authors [Sadeghian H, Rezazadeh G, Osterberg PM. Application of the generalized differential quadrature method to the study of pull-in phenomena of MEMS switches. J Microelectromech Syst 2007;16(6):1334–40] is discretized by applying Galerkin and GDQ methods. The divergence instability or pull-in phenomenon is analyzed. Obtained results are compared with the results of the pervious works. The Galerkin method is implemented with effect of number of used shape functions. Different types of trail functions on calculated pull-in voltage are examined.Furthermore, compare to one term and two terms truncation Galerkin method, it is observed that the GDQ with small number of grid points (non-uniform) performs accurate results for nonlinear micro-electro-mechanical coupled behavior which requires a large number of grid points at high-order approximation.     

车路耦合非线性振动高阶Galerkin截断研究

将移动车辆模型化为运动的两自由度质量-弹簧-阻尼系统,道路模型化为立方非线性黏弹性地基上的弹性梁,并将路面不平度设定为简谐函数．通过受力分析,建立车路非线性耦合振动高阶偏微分方程．采用高阶Galerkin截断结合数值方法求解耦合系统的动态响应．首次研究不同截断阶数对车路耦合非线性振动动态响应的影响,确定Galerkin截断研究车路耦合振动的收敛性．研究结果表明,对于软土地基的沥青路面,耦合振动的动态响应,需要150阶以上的截断才能达到收敛效果．并通过高阶收敛的Galerkin截断研究了系统参数对车路耦合非线性振动动态响应的影响.     

Application of the Chebyshev pseudospectral method to van der Waals fluids

In this paper, we consider a class of van der Waals flows with non-convex flux functions. In these flows, nonclassical under-compressive shock waves can develop. Such waves, which are characterized by kinetic functions, violate classical entropy conditions. We propose to use a Chebyshev pseudospectral method for solving the governing equations. A comparison of the results of this method with very high-order (up to 10th-order accurate) finite difference schemes is presented, which shows that the proposed method leads to a lower level of numerical oscillations than other high-order finite difference schemes and also does not exhibit fast-traveling packages of short waves which are usually observed in high-order finite difference methods. The proposed method can thus successfully capture various complex regimes of waves and phase transitions in both elliptic and hyperbolic regimes.     

A Novel Staggered Semi-implicit Space-Time Discontinuous Galerkin Method for the Incompressible Navier-Stokes Equations

A new high-order accurate staggered semi-implicit space-time discontinuous Galerkin(DG)method is presented for the simulation of viscous incompressible flows on...     

A high-order non-conforming discontinuous Galerkin method for time-domain electromagnetics

In this paper, we discuss the formulation, stability and validation of a high-order non-dissipative discontinuous Galerkin (DG) method for solving Maxwell’s equations on non-conforming simplex meshes. The proposed method combines a centered approximation for the numerical fluxes at inter element boundaries, with either a second-order or a fourth-order leap-frog time integration scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary-level hanging nodes. The method is proved to be stable and conserves a discrete counterpart of the electromagnetic energy for metallic cavities. Numerical experiments with high-order elements show the potential of the method.     

Numerical study of one-dimensional Vlasov-Poisson equations for infinite homogeneous stellar systems

We consider the gravitational Vlasov-Poisson (VP), or the so-called collisionless Boltzmann-Poisson equations for the self-gravitating collisionless stellar systems. We compute the solutions using a high-order discontinuous Galerkin method for the Vlasov equation, and the classical representation by Green’s function for the Poisson equation in the one-dimensional setting. We study both the case of damping and Jeans instability depending on the wavenumbers, which are taken to be greater than or less than the Jeans wavenumber, respectively. The method is shown to be stable, accurate and conservative. We report for the first time the BGK modes for the gravitational VP system, and we study the behavior of solutions associated with these various wavenumbers.     

A TVD Type Wavelet-Galerkin Method for Hamilton-Jacobi Equations

In this paper,we use Daubechies scaling functions as test functions for the Galerkin method,and discuss Wavelet-Galerkin solutions for the Hamilton-Jacobi equations.It can be proved that the schemesare TVD schemes.Numerical tests indicate that the schemes are suitable for the Hamilton-Jacobi equations.Furthermore,they have high-order accuracy in smooth regions and good resolution of singularities.     

曲边区域上的多边形网格间断有限元离散及其多重网格算法

本文利用多边形网格上的间断有限元方法离散二阶椭圆方程,在曲边区域上,采用多条直短边逼近曲边的以直代曲的策略,实现了高阶元在能量范数下的最优收敛.本文还将这一方法用于带曲边界面问题的求解,同样得到高阶元的最优收敛.此外我们还设计并分析了这一方法的\linebreakW-cycle和Variable V-cycle多重网格预条件方法,证明当光滑次数足够多时,多重网格预条件算法一致收敛.最后给出了数值算例,证实该算法的可行性并验证了理论分析的结果.     

Newton-Multigrid for Biological Reaction-Diffusion Problems with Random Coefficients

An algebraic Newton-multigrid method is proposed in order to efficiently solve systems of nonlinear reaction-diffusion problems with stochastic coefficients. These problems model the conversion of starch into sugars in growing apples. The stochastic system is first converted into a large coupled system of deterministic equations by applying a stochastic Galerkin finite element discretization. This method leads to high-order accurate stochastic solutions. A stable and high-order time discretization is obtained by applying a fully implicit Runge-Kutta method. After Newton linearization, a point-based algebraic multigrid solution method is applied. In order to decrease the computational cost, alternative multigrid preconditioners are presented. Numerical results demonstrate the convergence properties, robustness and efficiency of the proposed multigrid methods.     

Compact high-order finite-element method for elliptic transport problems with variable coefficients

We present a new conforming bilinear Petrov-Galerkin finite-element scheme for elliptic transport problems with variable coefficients. This scheme combines a generalized test function and artificial diffusion to achieve O(h⁴) grid-point accuracy on uniform stencils of 3 × 3 in two dimensions without resorting to the extended stencils of high-order elements. The method is compared with upwind and high-order finite-difference schemes and the standard Galerkin finite-element method for representative test problems. © 1994 John Wiley & Sons, Inc.     

High-Order Leap-Frog Based Discontinuous Galerkin Method for the Time-Domain Maxwell Equations on Non-Conforming Simplicial Meshes

A high-order leap-flog based non-dissipative discontinuous Galerkin timedomain method for solving Maxwell's equations is introduced and analyzed. The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements, with a Nth-order leap-frog time scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with highorder elements show the potential of the method.     

Comparison of light transport models based on finite element and the discrete ordinates methods in view of optical tomography applications

Numerical tests are used to evaluate the accuracy of two finite element formulations associated with the discrete ordinates method for solving the radiative transfer equation: the Least Square and the Discontinuous Galerkin finite element formulations. The results show that the use of a penalization method to set the Dirichlet boundary conditions leads to a more accurate solution than the weakly type setting where the Least Square method is seen to be more sensitive. Convergence in mesh size shows that, while both methods give accurate results, the Discontinuous Galerkin formulation uses five times more degrees of freedom than the Least Square formulation, which may lead to large systems to handle when the number of mesh elements is large. The comparison of both methods using the S_n and the T_n angular quadratures has shown that the Discontinuous Galerkin gives more accurate solutions, as expected, for problems with strong discontinuities, but may exhibit some oscillations due to the Galerkin procedure. A last test featuring a collimated irradiation shows that both methods give the same accuracy due to the separation of the radiative intensity into transmitted and scattered components, which removes the discontinuities in the implementation of the boundary conditions.     

High-order accurate thin layer approximations for time-domain electromagnetics. Part I: General metal backed coatings

We present a family of high-order accurate thin layer approximations for time-domain electromagnetics. The thin layer approximations are valid for metal backed coatings of general isotropic materials and certain classes of anisotropic materials on smooth curvilinear backgrounds. Both dielectric and magnetic materials can be considered. These models are non-trivial and we discuss their properties and implementation in the context of discontinuous Galerkin methods which are particularly well suited for these models. The range of validity, accuracy, and stability of the resulting schemes is demonstrated through one- and two-dimensional examples.     

Excitation of electromagnetic oscillations in a domain with chiral filling

The problem of excitation of electromagnetic oscillations by a given distribution of charges and currents in a domain with inhomogeneous chiral filling is examined. The domain in which the problem is considered may either be finite with a perfectly conducting boundary surface or be the complement of a perfectly conducting bounded body. A special functional space is defined on which a generalized initial-boundary value problem is formulated. The Galerkin method is used to prove the existence and uniqueness of a weak solution of this problem.     

A parallel spectral Fourier–Nonlinear Galerkin algorithm for simulation of turbulence

We present a high-order parallel algorithm, which requires only the minimum interprocessor communication dictated by the physical nature of the problem at hand. The parallelization is achieved by domain decomposition. The discretization in space is performed using the Local Fourier Basis method. The continuity conditions on the interfaces are enforced by adding homogeneous solutions. Such solutions often have fast decay properties, which can be utilized to minimize interprocessor communication. In effect, the predominant part of the computation is performed independently in the subdomains (processors) or using only local communication. A novel element of the present parallel algorithm is the incorporation of a Nonlinear Galerkin strategy to accelerate the computation and stabilize the time integration process. The basic idea of this approach consists of decomposition of the variables into large scale and small scale components with different treatment of these large and small scales. The combination of the Multidomain Fourier techniques with the Nonlinear Galerkin (NLG) algorithm is applied here to solve incompressible Navier–Stokes equations. Results are presented on direct numerical simulation of two-dimensional homogeneous turbulence using the NLG method. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 699–715, 1997     

Mesh Curving and Refinement Based on Cubic Bézier Surface for High-Order Discontinuous Galerkin Methods

Computational Mathematics and Mathematical Physics - In this paper, three-dimensional mesh curving and refinement methods are examined for high-order flow simulations with discontinuous Galerkin...     

An optimal rational basis function in the finite element method for convection-diffusion problems in one space variables

A new fixed grid Galerkin finite element technique using a rational basis function of order (1, 1) to solve convection-dominated phenomenon is described. A special feature of the method is that it employs rational basis functions, which are biased in the upstream direction and thus replace the equivalent upwinding techniques. Optimality of the rational basis function results in an accurate and stable numerical scheme which is used to avoid numerical oscillations at high Peclet numbers. Numerical examples are discussed which illustrate the efficiency of the method.     

Discontinuous Galerkin finite element heterogeneous multiscale method for advection–diffusion problems with multiple scales

A discontinuous Galerkin finite element heterogeneous multiscale method is proposed for advection–diffusion problems with highly oscillatory coefficients. The method is based on a coupling of a discontinuous Galerkin discretization for an effective advection–diffusion problem on a macroscopic mesh, whose a priori unknown data are recovered from micro finite element calculations on sampling domains within each macro element. The computational work involved is independent of the high oscillations in the problem at the smallest scale. The stability of our method (depending on both macro and micro mesh sizes) is established for both diffusion dominated and advection dominated regimes without any assumptions about the type of heterogeneities in the data. Fully discrete a priori error bounds are derived for locally periodic data. Numerical experiments confirm the theoretical error estimates.