Alexander Bihlo

Symmetry preserving parameterization schemes

Methods for the design of physical parameterization schemes that possess certain invariance properties are discussed. These methods are based on different techniques of group classification and provide means to determine expressions for unclosed terms arising in the course of averaging of nonlinear differential equations. The demand that the averaged equation is invariant with respect to a subalgebra of the maximal Lie invariance algebra of the unaveraged equation leads to a problem of inverse group classification which is solved by the description of differential invariants of the selected subalgebra. Given no prescribed symmetry group, the direct group classification problem is relevant. Within this framework, the algebraic method or direct integration of determining equations for Lie symmetries can be applied. For cumbersome parameterizations, a preliminary group classification can be carried out. The methods presented are exemplified by parameterizing the eddy vorticity flux in the averaged vorticity ...

Complete group classification of a class of nonlinear wave equations

Preliminary group classification became a prominent tool in the symmetry analysis of differential equations due to the paper by Ibragimov, Torrisi, and Valenti [J. Math. Phys. 32, 2988–2995 (1991)10.1063/1.529042]. In this paper the partial preliminary group classification of a class of nonlinear wave equations was carried out via the classification of one-dimensional Lie symmetry extensions related to a fixed finite-dimensional subalgebra of the infinite-dimensional equivalence algebra of the class under consideration. In the present paper we implement the complete group classification of the same class up to both usual and general point equivalence using the algebraic method of group classification. This includes the complete preliminary group classification of the class and finding those Lie symmetry extensions which are not associated with subalgebras of the equivalence algebra. The complete preliminary group classification is based on listing all inequivalent subalgebras of the whole infinite-dimensi...

Enhanced preliminary group classification of a class of generalized diffusion equations

The method of preliminary group classification is rigorously defined, enhanced and related to the theory of group classification of differential equations. Typical weaknesses in papers on this method are discussed and strategies to overcome them are presented. The preliminary group classification of the class of generalized diffusion equations of the form ut=f(x,u)ux2+g(x,u)uxx is carried out. This includes a justification for applying this method to the given class, the simultaneous computation of the equivalence algebra and equivalence (pseudo) group, as well as the classification of inequivalent appropriate subalgebras of the whole infinite-dimensional equivalence algebra. The extensions of the kernel algebra, which are induced by such subalgebras, are exhaustively described. These results improve those recently published in Commun Nonlinear Sci Numer Simul.

Lie symmetries and exact solutions of the barotropic vorticity equation

Lie group methods are used for the study of various issues related to symmetries and exact solutions of the barotropic vorticity equation. The Lie symmetries of the barotropic vorticity equations on the f- and β-planes, as well as on the sphere in rotating and rest reference frames, are determined. A symmetry background for reducing the rotating reference frame to the rest frame is presented. The one- and two-dimensional inequivalent subalgebras of the Lie invariance algebras of both equations are exhaustively classified and then used to compute invariant solutions of the vorticity equations. This provides large classes of exact solutions, which include both Rossby and Rossby–Haurwitz waves as special cases. We also discuss the possibility of partial invariance for the β-plane equation, thereby further extending the family of its exact solutions. This is done in a more systematic and complete way than previously available in literature.

Invariant Discretization Schemes for the Shallow-Water Equations

Invariant discretization schemes are derived for the one- and two-dimensional shallow-water equations with periodic boundary conditions. While originally designed for constructing invariant finite difference schemes, we extend the use of difference invariants to allow construction of invariant finite volume methods as well. It is found that the classical invariant schemes converge to the Lagrangian formulation of the shallow-water equations. These schemes require redistribution of the grid points according to the physical fluid velocity; i.e., the mesh cannot remain fixed in the course of the numerical integration. Invariant Eulerian discretization schemes are proposed for the shallow-water equations in computational coordinates. Instead of using the fluid velocity as the grid velocity, an invariant moving mesh generator is invoked in order to determine the location of the grid points at the subsequent time level. The numerical conservation of energy, mass, and momentum is evaluated for both the invariant...

Rayleigh–Bénard convection as a Nambu-metriplectic problem

The traditional Hamiltonian structure of the equations governing conservative Rayleigh-Benard convection (RBC) is singular, i.e., its Poisson bracket possesses nontrivial Casimir functionals. We show that a special form of one of these Casimirs can be used to extend the bilinear Poisson bracket to a trilinear generalized Nambu bracket. It is further shown that the equations governing dissipative RBC can be written as the superposition of the conservative Nambu bracket with a dissipative symmetric bracket. This leads to a Nambu-metriplectic system, which completes the geometrical picture of RBC.

Lie symmetry analysis and exact solutions of the quasigeostrophic two-layer problem

The quasigeostrophic two-layer model is of superior interest in dynamic meteorology since it is one of the easiest ways to study baroclinic processes in geophysical fluid dynamics. The complete set of point symmetries of the two-layer equations is determined. An optimal set of one- and two-dimensional inequivalent subalgebras of the maximal Lie invariance algebra is constructed. On the basis of these subalgebras, we exhaustively carry out group-invariant reduction and compute various classes of exact solutions. Wherever possible, reference to the physical meaning of the exact solutions is given. In particular, the well-known baroclinic Rossby wave solutions in the two-layer model are rediscovered.

Group classification of linear evolution equations

Abstract The group classification problem for the class of (1+1)-dimensional linear r th order evolution equations is solved for arbitrary values of r > 2 . It is shown that a related maximally gauged class of homogeneous linear evolution equations is uniformly semi-normalized with respect to linear superposition of solutions and hence the complete group classification can be obtained using the algebraic method. We also compute exact solutions for equations from the class under consideration using Lie reduction and its specific generalizations for linear equations.

The Korteweg–de Vries equation and its symmetry-preserving discretization

The Korteweg-de Vries equation is one of the most important nonlinear evolution equations in the mathematical sciences. In this article invariant discretization schemes are constructed for this equation both in the Lagrangian and in the Eulerian form. We also propose invariant schemes that preserve the momentum. Numerical tests are carried out for all invariant discretization schemes and related to standard numerical schemes. We find that the invariant discretization schemes give generally the same level of accuracy as the standard schemes with the added benefit of preserving Galilean transformations which is demonstrated numerically as well.

Parallel stochastic methods for PDE based grid generation

Abstract The efficient generation of meshes is an important step in the numerical solution of various problems in physics and engineering. We are interested in situations where global mesh quality and tight coupling to the physical solution is important. We consider elliptic PDE based mesh generation and present a method for the construction of adaptive meshes in two spatial dimensions using domain decomposition that is suitable for an implementation on parallel computing architectures. The method uses the stochastic representation of the exact solution of a linear mesh generator of Winslow type to find the points of the adaptive mesh along the subdomain interfaces. The meshes over the single subdomains can then be obtained completely independently of each other using the probabilistically computed solutions along the interfaces as boundary conditions for the linear mesh generator. Further to the previously acknowledged performance characteristics, we demonstrate how the stochastic domain decomposition approach is particularly suited to the problem of grid generation — generating quality meshes efficiently. In addition we show further improvements are possible using interpolation of the subdomain interfaces and smoothing of mesh candidates. An optimal placement strategy is introduced to automatically choose the number and placement of points along the interface using the mesh density function. Various examples of meshes constructed using this stochastic–deterministic domain decomposition technique are shown and compared to the respective single domain solutions using a representative mesh quality measure. A brief performance study is included to show the viability of the stochastic domain decomposition approach and to illustrate the effect of algorithmic choices on the solver’s efficiency.