Multidimensional method-of-lines transport for atmospheric flows over steep terrain using arbitrary meshes

Including terrain in atmospheric models gives rise to mesh distortions near the lower boundary that can degrade accuracy and challenge the stability of transport schemes. Multidimensional transport schemes avoid splitting errors on distorted, arbitrary meshes, and method-of-lines schemes have low computational cost because they perform reconstructions at fixed points. This paper presents a multidimensional method-of-lines finite volume transport scheme, “cubicFit”, which is designed to be numerically stable on arbitrary meshes. Stability conditions derived from a von Neumann analysis are imposed during model initialisation to obtain stability and improve accuracy in distorted regions of the mesh, and near steeply-sloping lower boundaries. Reconstruction calculations depend upon the mesh only, needing just one vector multiply per face per time-stage irrespective of the velocity field. The cubicFit scheme is evaluated using three, idealised numerical tests. The first is a variant of a standard horizontal transport test on severely distorted terrain-following meshes. The second is a new test case that assesses accuracy near the ground by transporting a tracer at the lower boundary over steep terrain on terrain-following meshes, cut-cell meshes, and new, slanted-cell meshes that do not suffer from severe time-step constraints associated with cut cells. The third, standard test deforms a tracer in a vortical flow on hexagonal-icosahedral meshes and cubed-sphere meshes. In all tests, cubicFit is stable and largely insensitive to mesh distortions, and cubicFit results are more accurate than those obtained using a multidimensional linear upwind transport scheme. The cubicFit scheme is second-order convergent regardless of mesh distortions

A mimetic, semi-implicit, forward-in-time, finite volume shallow water model: comparison of hexagonal–icosahedral and cubed-sphere grids

A new algorithm is presented for the solution of the shallow water equations on quasi-uniform spherical grids. It combines a mimetic finite volume spatial discretization with a Crank–Nicolson time discretization of fast waves and an accurate and conservative forward-in-time advection scheme for mass and potential vorticity (PV). The algorithm is implemented and tested on two families of grids: hexagonal–icosahedral Voronoi grids, and modified equiangular cubed-sphere grids. <br><br> Results of a variety of tests are presented, including convergence of the discrete scalar Laplacian and Coriolis operators, advection, solid body rotation, flow over an isolated mountain, and a barotropically unstable jet. The results confirm a number of desirable properties for which the scheme was designed: exact mass conservation, very good available energy and potential enstrophy conservation, consistent mass, PV and tracer transport, and good preservation of balance including vanishing &nabla; &times; &nabla;, steady geostrophic modes, and accurate PV advection. The scheme is stable for large wave Courant numbers and advective Courant numbers up to about 1. <br><br> In the most idealized tests the overall accuracy of the scheme appears to be limited by the accuracy of the Coriolis and other mimetic spatial operators, particularly on the cubed-sphere grid. On the hexagonal grid there is no evidence for damaging effects of computational Rossby modes, despite attempts to force them explicitly

Three-phase AC-AC hexagonal chopper system with heterodyne modulation for power flow control enhancement

This paper proposes a three-phase AC chopper system for the interconnection of various distributed generation (DG) farms or main utilities to enhance the active and reactive power flow control. The absence of large energy storage component in direct AC-AC converter makes the system footprint small and reliable. As the interface for different AC sources, the presented converter can be configured as star or delta. However, delta connection is preferred as it can trap the potential zero-sequence current and reduce the current rating of the switching devices. In this way, the proposed converter resembles the hexagonal chopper, and it offers an inherent degree of freedom for output voltage phase-shifting. Considering the scalability in high voltage applications, a new version of the hexagonal chopper with half-bridge cell modular multilevel structure is developed. The modular multilevel AC hexagonal chopper (M2AHC) is operated in quasi-2-level mode to suppress the electro-magnetic interference (EMI) caused by high voltage switching. Quasi-2-level operation divides the voltage level transition into multi-steps, diminishing the voltage rising and falling rates (dv/dt) in high voltage condition. Then, heterodyne modulation is adopted for the presented chopper system, supplying a new degree of freedom to decouple the phase and amplitude regulation. Based on this idea, system control strategy is developed in synchronous reference frame (SRF). Simulations and experimentations have confirmed the validity of the proposed approaches

A primal-dual mimetic finite element scheme for the rotating shallow water equations on polygonal spherical meshes

Copyright © 2015 Elsevier. NOTICE: this is the author’s version of a work that was accepted for publication in Journal of Computational Physics. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Journal of Computational Physics Vol. 290 (2015), DOI: 10.1016/j.jcp.2015.02.045A new numerical method is presented for solving the shallow water equations on a rotating sphere using quasi-uniform polygonal meshes. The method uses special families of finite element function spaces to mimic key mathematical properties of the continuous equations and thereby capture several desirable physical properties related to balance and conservation. The method relies on two novel features. The first is the use of compound finite elements to provide suitable finite element spaces on general polygonal meshes. The second is the use of dual finite element spaces on the dual of the original mesh, along with suitably defined discrete Hodge star operators to map between the primal and dual meshes, enabling the use of a finite volume scheme on the dual mesh to compute potential vorticity fluxes. The resulting method has the same mimetic properties as a finite volume method presented previously, but is more accurate on a number of standard test cases.Natural Environment Research Council under the “GungHo” projec