A partially mesh-free scheme for representing anisotropic spatial variations along field lines

A common numerical task is to represent functions which are highly spatially anisotropic, and to solve differential equations related to these functions. One way such anisotropy arises is that information transfer along one spatial direction is much faster than in others. In this situation, the derivative of the function is small in the local direction of a vector field B. In order to define a discrete representation, a set of surfaces Mi indexed by an integer i are chosen such that mapping along the field B induces a one-to-one relation between the points on surface Mi to those on Mi+1. For simple cases Mi may be surfaces of constant coordinate value. On each surface Mi, a function description is constructed using basis functions defined on a regular structured mesh. The definition of each basis function is extended from the surface M along the lines of the field B by multiplying it by a smooth compact support function whose argument increases with distance along B. Function values are evaluated by summing contributions associated with each surface Mi. This does not require any special connectivity of the meshes used in the neighbouring surfaces M, which substantially simplifies the meshing problem compared to attempting to find a space filling anisotropic mesh. We explore the numerical properties of the scheme, and show that it can be used to efficiently solve differential equations for certain anisotropic problems.