The purpose of this paper is to derive the anisotropic averaged Euler equations and to study their geometric and analytic properties. These new equations involve the evolution of a mean velocity field and an advected symmetric tensor that captures the fluctuation effects. Besides the derivation of these equations, the new results in the paper are smoothness properties of the equations in material representation, which gives well-posedness of the equations, and the derivation of a corrector to the macroscopic velocity field. The numerical implementation and physical implications of this set of equations will be explored in other publications.

In this paper, discrete analogues of Euler-Poincar {e} and Lie-Poisson reduction theory are developed for systems on finite dimensional Lie groups $G$ with Lagrangians $L:TG \to {\mathbb R}$ that are $G$-invariant. These discrete equations provide ``reduced'' numerical algorithms which manifestly preserve the symplectic structure. The manifold $G \times G$ is used as an approximation of $TG$, and a discrete Langragian ${\mathbb L}:G \times G \to {\mathbb R}$ is construced in such a way that the $G$-invariance property is preserved. Reduction by $G$ results in new ``variational'' principle for the reduced Lagrangian $\ell:G \to {\mathbb R}$, and provides the discrete Euler-Poincar {e} (DEP) equations. Reconstruction of these equations recovers the discrete Euler-Lagrange equations developed in \cite{MPS,WM} which are naturally symplectic-momentum algorithms. Furthermore, the solution of the DEP algorithm immediately leads to a discrete Lie-Poisson (DLP) algorithm. It is shown that when $G=\text{SO} (n)$, the DEP and DLP algorithms for a particular choice of the discrete Lagrangian ${\mathbb L}$ are equivalent to the Moser-Veselov scheme for the generalized rigid body. %As an application, a reduced symplectic integrator for two dimensional %hydrodynamics is constructed using the SU$(n)$ approximation to the volume %preserving diffeomorphism group of ${\mathbb T}^2$.

For a discrete mechanical system on a Lie group $G$ determined by a (reduced) Lagrangian $\ell$ we define a Poisson structure via the pull-back of the Lie-Poisson structure on the dual of the Lie algebra ${\mathfrak g}^*$ by the corresponding Legendre transform. The main result shown in this paper is that this structure coincides with the reduction under the symmetry group $G$ of the canonical discrete Lagrange 2-form $\omega_\mathbb{L}$ on $G \times G$. Its symplectic leaves then become dynamically invariant manifolds for the reduced discrete system. Links between our approach and that of groupoids and algebroids as well as the reduced Hamilton-Jacobi equation are made. The rigid body is discussed as an example.

This paper presents a variational and multisymplectic formulation of both compressible and incompressible models of continuum mechanics on general Riemannian manifolds. A general formalism is developed for non-relativistic first-order multisymplectic field theories with constraints, such as the incompressibility constraint. The results obtained in this paper set the stage for multisymplectic reduction and for the further development of Veselov-type multisymplectic discretizations and numerical algorithms. The latter will be the subject of a companion paper.

Connections on principal bundles play a fundamental role in expressing the equations of motion for mechanical systems with symmetry in an intrinsic fashion. A discrete theory of connections on principal bundles is constructed by introducing the discrete analogue of the Atiyah sequence, with a connection corresponding to the choice of a splitting of the short exact sequence. Equivalent representations of a discrete connection are considered, and an extension of the pair groupoid composition, that takes into account the principal bundle structure, is introduced. Computational issues, such as the order of approximation, are also addressed. Discrete connections provide an intrinsic method for introducing coordinates on the reduced space for discrete mechanics, and provide the necessary discrete geometry to introduce more general discrete symmetry reduction. In addition, discrete analogues of the Levi-Civita connection, and its curvature, are introduced by using the machinery of discrete exterior calculus, and discrete connections.