Anisotropic model for granulated sea ice dynamics

A continuum model describing sea ice as a layer of granulated thick ice, consisting of many rigid, brittle floes, intersected by long and narrow regions of thinner ice, known as leads, is developed. We consider the evolution of mesoscale leads, formed under extension, whose lengths span many floes, so that the surrounding ice is treated as a granular plastic. The leads are sufficiently small with respect to basin scales of sea ice deformation that they may be modelled using a continuum approach. The model includes evolution equations for the orientational distribution of leads, their thickness and width expressed through second-rank tensors and terms requiring closures. The closing assumptions are constructed for the case of negligibly small lead ice thickness and the canonical deformation types of pure and simple shear, pure divergence and pure convergence. We present a new continuum-scale sea ice rheology that depends upon the isotropic, material rheology of sea ice, the orientational distribution of lead properties and the thick ice thickness. A new model of lead and thick ice interaction is presented that successfully describes a number of effects: (i) because of its brittle nature, thick ice does not thin under extension and (ii) the consideration of the thick sea ice as a granular material determines finite lead opening under pure shear, when granular dilation is unimportant.

A continuum anisotropic model of sea-ice dynamics

We develop the essential ingredients of a new, continuum and anisotropic model of sea-ice dynamics designed for eventual use in climate simulation. These ingredients are a constitutive law for sea-ice stress, relating stress to the material properties of sea ice and to internal variables describing the sea-ice state, and equations describing the evolution of these variables. The sea-ice cover is treated as a densely flawed two-dimensional continuum consisting of a uniform field of thick ice that is uniformly permeated with narrow linear regions of thinner ice called leads. Lead orientation, thickness and width distributions are described by second-rank tensor internal variables: the structure, thickness and width tensors, whose dynamics are governed by corresponding evolution equations accounting for processes such as new lead generation and rotation as the ice cover deforms. These evolution equations contain contractions of higher-order tensor expressions that require closures. We develop a sea-ice stress constitutive law that relates sea-ice stress to the structure tensor, thickness tensor and strain rate. For the special case of empty leads (containing no ice), linear closures are adopted and we present calculations for simple shear, convergence and divergence.

The influence of anisotropic F region ion velocity distributions on ionospheric ion outflows into the magnetosphere

The contribution to the field-aligned ionospheric ion momentum equation, due to coupling between pressure anisotropy and the inhomogeneous geomagnetic field, is investigated. We term this contribution the “hydrodynamic mirror force” and investigate its dependence on the ion drift and the resulting deformations of the ion velocity distribution function from an isotropic form. It is shown that this extra upforce increases rapidly with ion drift relative to the neutral gas but is not highly dependent on the ion-neutral collision model employed. An example of a burst of flow observed by EISCAT, thought to be the ionospheric signature of a flux transfer event at the magnetopause, is studied in detail and it is shown that the nonthermal plasma which results is subject to a hydrodynamic mirror force which is roughly 10% of the gravitational downforce. In addition, predictions by the coupled University College London-Sheffield University model of the ionosphere and thermosphere show that the hydrodynamic mirror force in the auroral oval is up to 3% of the gravitational force for Kp of about 3, rising to 10% following a sudden increase in cross-cap potential. The spatial distribution of the upforce shows peaks in the cusp region and in the post-midnight auroral oval, similar to that of observed low-energy heavy ion flows from the ionosphere into the magnetosphere. We suggest the hydrodynamic mirror force may modulate these outflows by controlling the supply of heavy ions to regions of ion acceleration and that future simulations of the effects of Joule heating on ion outflows should make allowance for it.

Modelling the rheology of anisotropic particles adsorbed on a two-dimensional fluid interface

We present a general approach based on nonequilibrium thermodynamics for bridging the gap between a well-defined microscopic model and the macroscopic rheology of particle-stabilised interfaces. Our approach is illustrated by starting with a microscopic model of hard ellipsoids confined to a planar surface, which is intended to simply represent a particle-stabilised fluid–fluid interface. More complex microscopic models can be readily handled using the methods outlined in this paper. From the aforementioned microscopic starting point, we obtain the macroscopic, constitutive equations using a combination of systematic coarse-graining, computer experiments and Hamiltonian dynamics. Exemplary numerical solutions of the constitutive equations are given for a variety of experimentally relevant flow situations to explore the rheological behaviour of our model. In particular, we calculate the shear and dilatational moduli of the interface over a wide range of surface coverages, ranging from the dilute isotropic regime, to the concentrated nematic regime.

The importance of anisotropic Coulomb interaction in LaMnO3

In low-temperature anti-ferromagnetic LaMnO3, strong and localized electronic interactions among Mn 3d electrons prevent a satisfactory description from standard local density and generalized gradient approximations in density functional theory calculations. Here we show that the strong on-site electronic interactions are described well only by using direct and exchange corrections to the intra-orbital Coulomb potential. Only DFT+U calculations with explicit exchange corrections produce a balanced picture of electronic, magnetic and structural observables in agreement with experiment. To understand the reason, a rewriting of the functional form of the +U corrections is presented that leads to a more physical and transparent understanding of the effect of these correction terms. The approach highlights the importance of Hund’s coupling (intra-orbital exchange) in providing anisotropy across the occupation and energy eigenvalues of the Mn d states. This intra-orbital exchange is the key to fully activating the Jahn-Teller distortion, reproducing the experimental band gap and stabilizing the correct magnetic ground state in LaMnO3. The best parameter values for LaMnO3 within the DFT(PBEsol)+U framework are determined to be U = 8 eV and J = 1.9 eV.