2,307 research outputs found
Mixed finite elements for numerical weather prediction
We show how two-dimensional mixed finite element methods that satisfy the
conditions of finite element exterior calculus can be used for the horizontal
discretisation of dynamical cores for numerical weather prediction on
pseudo-uniform grids. This family of mixed finite element methods can be
thought of in the numerical weather prediction context as a generalisation of
the popular polygonal C-grid finite difference methods. There are a few major
advantages: the mixed finite element methods do not require an orthogonal grid,
and they allow a degree of flexibility that can be exploited to ensure an
appropriate ratio between the velocity and pressure degrees of freedom so as to
avoid spurious mode branches in the numerical dispersion relation. These
methods preserve several properties of the C-grid method when applied to linear
barotropic wave propagation, namely: a) energy conservation, b) mass
conservation, c) no spurious pressure modes, and d) steady geostrophic modes on
the -plane. We explain how these properties are preserved, and describe two
examples that can be used on pseudo-uniform grids: the recently-developed
modified RT0-Q0 element pair on quadrilaterals and the BDFM1-\pdg element pair
on triangles. All of these mixed finite element methods have an exact 2:1 ratio
of velocity degrees of freedom to pressure degrees of freedom. Finally we
illustrate the properties with some numerical examples.Comment: Revision after referee comment
Numerical wave propagation for the triangular - finite element pair
Inertia-gravity mode and Rossby mode dispersion properties are examined for
discretisations of the linearized rotating shallow-water equations using the
- finite element pair on arbitrary triangulations in planar
geometry. A discrete Helmholtz decomposition of the functions in the velocity
space based on potentials taken from the pressure space is used to provide a
complete description of the numerical wave propagation for the discretised
equations. In the -plane case, this decomposition is used to obtain
decoupled equations for the geostrophic modes, the inertia-gravity modes, and
the inertial oscillations. As has been noticed previously, the geostrophic
modes are steady. The Helmholtz decomposition is used to show that the
resulting inertia-gravity wave equation is third-order accurate in space. In
general the \pdgp finite element pair is second-order accurate, so this leads
to very accurate wave propagation. It is further shown that the only spurious
modes supported by this discretisation are spurious inertial oscillations which
have frequency , and which do not propagate. The Helmholtz decomposition
also allows a simple derivation of the quasi-geostrophic limit of the
discretised - equations in the -plane case, resulting in a
Rossby wave equation which is also third-order accurate.Comment: Revised version prior to final journal submissio
Enforcing the non-negativity constraint and maximum principles for diffusion with decay on general computational grids
In this paper, we consider anisotropic diffusion with decay, and the
diffusivity coefficient to be a second-order symmetric and positive definite
tensor. It is well-known that this particular equation is a second-order
elliptic equation, and satisfies a maximum principle under certain regularity
assumptions. However, the finite element implementation of the classical
Galerkin formulation for both anisotropic and isotropic diffusion with decay
does not respect the maximum principle.
We first show that the numerical accuracy of the classical Galerkin
formulation deteriorates dramatically with increase in the decay coefficient
for isotropic medium and violates the discrete maximum principle. However, in
the case of isotropic medium, the extent of violation decreases with mesh
refinement. We then show that, in the case of anisotropic medium, the classical
Galerkin formulation for anisotropic diffusion with decay violates the discrete
maximum principle even at lower values of decay coefficient and does not vanish
with mesh refinement. We then present a methodology for enforcing maximum
principles under the classical Galerkin formulation for anisotropic diffusion
with decay on general computational grids using optimization techniques.
Representative numerical results (which take into account anisotropy and
heterogeneity) are presented to illustrate the performance of the proposed
formulation
An efficient discontinuous Galerkin method for aeroacoustic propagation
An efficient discontinuous Galerkin formulation is applied to the solution of the linearized Euler equations and the acoustic perturbation equations for the simulation of aeroacoustic propagation in two-dimensional and axisymmetric problems, with triangular and quadrilateral elements. To improve computational efficiency, a new strategy of variable interpolation order is proposed in addition to a quadrature-free approach and parallel implementation. Moreover, an accurate wall boundary condition is formulated on the basis of the solution of the Riemann problem for a reflective wall. Time discretization is based on a low dissipation formulation of a fourth-order, low storage Runge-Kutta scheme. Along the far-field boundaries a perfectly matched layer boundary condition is used. For the far-field computations, the integral formulation of Ffowcs Williams and Hawkings is coupled with the near-field solver. The efficiency and accuracy of the proposed variable order formulation is assessed for realistic geometries, namely sound propagation around a high-lift airfoil and the Munt problem
Spectral methods for CFD
One of the objectives of these notes is to provide a basic introduction to spectral methods with a particular emphasis on applications to computational fluid dynamics. Another objective is to summarize some of the most important developments in spectral methods in the last two years. The fundamentals of spectral methods for simple problems will be covered in depth, and the essential elements of several fluid dynamical applications will be sketched
- …