Embedded discontinuous Galerkin transport schemes with localised limiters

Motivated by finite element spaces used for representation of temperature in the compatible finite element approach for numerical weather prediction, we introduce locally bounded transport schemes for (partially-)continuous finite element spaces. The underlying high-order transport scheme is constructed by injecting the partially-continuous field into an embedding discontinuous finite element space, applying a stable upwind discontinuous Galerkin (DG) scheme, and projecting back into the partially-continuous space; we call this an embedded DG scheme. We prove that this scheme is stable in L2 provided that the underlying upwind DG scheme is. We then provide a framework for applying limiters for embedded DG transport schemes. Standard DG limiters are applied during the underlying DG scheme. We introduce a new localised form of element-based flux-correction which we apply to limiting the projection back into the partially-continuous space, so that the whole transport scheme is bounded. We provide details in the specific case of tensor-product finite element spaces on wedge elements that are discontinuous P1/Q1 in the horizontal and continuous P2 in the vertical. The framework is illustrated with numerical tests

On two-dimensional finite amplitude electro-convection in a dielectric liquid induced by a strong unipolar injection

The hydrodynamic stability of a dielectric liquid subjected to strong unipolar injection is numerically investigated. We determined the linear criterion Tc (T being the electric Rayleigh number) and finite amplitude one Tf over a wide range of the mobility parameter M. A noticeable discrepancy is shown for Tf between our numerical prediction and the value predicted by stability analysis, which is due to the velocity field used in stability analysis. Recent studies revealed a transition of the flow structure from one cell to two with an increase in T. We demonstrate that this transition results in a new subcritical bifurcationMinisterio de Ciencia y Tecnología FIS2011-25161Junta de Andalucía P10-FQM-5735Junta de Andalucía P09-FQM-458

A unified analysis of algebraic flux correction schemes for convection–diffusion equations

Recent results on the numerical analysis of algebraic flux correction (AFC) finite element schemes for scalar convection–diffusion equations are reviewed and presented in a unified way. A general form of the method is presented using a link between AFC schemes and nonlinear edge-based diffusion schemes. Then, specific versions of the method, that is, different definitions for the flux limiters, are reviewed and their main results stated. Numerical studies compare the different versions of the scheme

A unified analysis of Algebraic Flux Correction schemes for convection-diffusion equations

Recent results on the numerical analysis of Algebraic Flux Correction (AFC) finite element schemes for scalar convection-diffusion equations are reviewed and presented in a unified way. A general form of the method is presented using a link between AFC schemes and nonlinear edge-based diffusion scheme. Then, specific versions of the method, this is, different definitions for the flux limiters, are reviewed and their main results stated. Numerical studies compare the different versions of the scheme

Institute for Computational Mechanics in Propulsion (ICOMP)

The Institute for Computational Mechanics in Propulsion (ICOMP) is a combined activity of Case Western Reserve University, Ohio Aerospace Institute (OAI) and NASA Lewis. The purpose of ICOMP is to develop techniques to improve problem solving capabilities in all aspects of computational mechanics related to propulsion. The activities at ICOMP during 1991 are described

A refined r-factor algorithm for TVD schemes on arbitrary unstructured meshes

A refined r-factor algorithm for implementing TVD schemes on arbitrary unstructured meshes, referred to henceforth as FFISAM (a Face-perpendicular Far-upwind Interpolation Scheme for Arbitrary Meshes), is proposed based on an extensive review of the existing r-factor algorithms available in the literature. The design principles, as well as the respective advantages and disadvantages, of the existing algorithms are first systematically analyzed before presenting the FFISAM. The FFISAM is designed to combine the merits of various existing r-factor algorithms. The performance of the FFISAM, implemented in ten classical TVD schemes, is evaluated against four two-dimensional pure-advection benchmark test cases where analytical solutions are available. The numerical results clearly show that the FFISAM leads to a better overall performance than the existing algorithms in terms of accuracy and convergence on arbitrary unstructured meshes for the ten classical TVD schemes.The authors gratefully acknowledge the financial support provided by the National Natural Science Foundation of China (Grant No. 51279082) and the support from Australian Research Council through a Discovery Grant (Project ID: DP110105171).This is the author accepted manuscript. The final version is available from Wiley at http://onlinelibrary.wiley.com/doi/10.1002/fld.4073/abstract