Slip at the surface of a sphere translating perpendicular to a plane wall in micropolar fluid

The Stokes axisymmetrical flow caused by a sphere translating in a micropolar fluid perpendicular to a plane wall at an arbitrary position from the wall is presented using a combined analytical-numerical method. A linear slip, Basset type, boundary condition on the surface of the sphere has been used. To solve the Stokes equations for the fluid velocity field and the microrotation vector, a general solution is constructed from fundamental solutions in both cylindrical, and spherical coordinate systems. Boundary conditions are satisfied first at the plane wall by the Fourier transforms and then on the sphere surface by the collocation method. The drag acting on the sphere is evaluated with good convergence. Numerical results for the hydrodynamic drag force and wall effect with respect to the micropolarity, slip parameters and the separation distance parameter between the sphere and the wall are presented both in tabular and graphical forms. Comparisons are made between the classical fluid and micropolar fluid.      

Stability of rank 2 vector bundles along isomonodromic deformations

We are interested in the stability of holomorphic rank 2 vector bundles of degree 0 over compact Riemann surfaces, which are provided with irreducible meromophic tracefree connections. In the case of a logarithmic connection on the Riemann sphere, such a vector bundle will be trivial up to the isomonodromic deformation associated to a small move of the poles, according to a result of A. Bolibruch. In the general case of meromorphic connections over Riemann surfaces of arbitrary genus, we prove that the vector bundle will be semi-stable, up to a small isomonodromic deformation. More precisely, the vector bundle underlying the universal isomonodromic deformation is generically semi-stable along the deformation, and even maximally stable. For curves of genus g ≥ 2, this result is non-trivial even in the case of non-singular connections. The author was partially supported by ANR SYMPLEXE BLAN06-3-137237.     

Multilevel methods for mixed finite elements in three dimensions

In this paper we consider second order scalar elliptic boundary value problems posed over three–dimensional domains and their discretization by means of mixed Raviart–Thomas finite elements [18]. This leads to saddle point problems featuring a discrete flux vector field as additional unknown. Following Ewing and Wang [26], the proposed solution procedure is based on splitting the flux into divergence free components and a remainder. It leads to a variational problem involving solenoidal Raviart–Thomas vector fields. A fast iterative solution method for this problem is presented. It exploits the representation of divergence free vector fields as s of the –conforming finite element functions introduced by Nédélec [43]. We show that a nodal multilevel splitting of these finite element spaces gives rise to an optimal preconditioner for the solenoidal variational problem: Duality techniques in quotient spaces and modern algebraic multigrid theory [50, 10, 31] are the main tools for the proof. Received November 4, 1996 / Revised version received February 2, 1998     

The three-dimensional problem for an anisotropic plate of arbitrary thickness

By expanding the components of the displacement vector in a certain system of functions of the transverse coordinate, we reduce the solution of the three-dimensional problem of the theory of elasticity of an anisotropic body to a series of two-dimensional problems. To determine the displacements we obtain a system of differential equations of infinite order with two independent variables. We show how to pass from the infinite system to a series of finite systems depending on the form of the external forces. Translated fromTeoreticheskaya i Prikladnaya Mekhanika, No. 28, 1998, pp. 11–19.     

Application of the multigrid approach for solving 3D Navier-Stokes equations on hexahedral grids using the discontinuous Galerkin method

The Galerkin method with discontinuous basis functions is adapted for solving the Euler and Navier-Stokes equations on unstructured hexahedral grids. A hybrid multigrid algorithm involving the finite element and grid stages is used as an iterative solution method. Numerical results of calculating the sphere inviscid flow, viscous flow in a bent pipe, and turbulent flow past a wing are presented. The numerical results and the computational cost are compared with those obtained using the finite volume method.     

Divergence-Free RBFs on Surfaces

This article presents a new tool for fitting a divergence-free vector field tangent to a two-dimensional orientable surface to samples of such a field taken at scattered sites on P. This method, which involves a kernel constructed from radial basis functions, has applications to problems in geophysics, and has the advantage of avoiding problems with poles. Numerical examples testing the method on the sphere are included.     

A boundary parametric approximation to the linearized scalar potential magnetostatic field problem

We consider the linearized scalar potential formulation of the magnetostatic field problem in this paper. Our approach involves a reformulation of the continuous problem as a parametric boundary problem. By the introduction of a spherical interface and the use of spherical harmonics, the infinite boundary conditions can also be satisfied in the parametric framework. That is the field in the exterior of a sphere is expanded in a ‘harmonic series’ of eigenfunctions for the exterior harmonic problem. The approach is essentially a finite element method coupled with a spectral method via a boundary parametric procedure. The reformulated problem is discretized by finite element techniques which leads to a discrete parametric problem which can be solved by well conditioned iteration involving only the solution of decoupled Neumann type elliptic finite element systems and L² projection onto subspaces of spherical harmonics. Error and stability estimates given show exponential convergence in the degree of the spherical harmonics and optimal order convergence with respect to the finite element approximation for the resulting fields in L².     

Moduli spaces of meromorphic connections and quiver varieties

We describe the moduli spaces of meromorphic connections on trivial holomorphic vector bundles over the Riemann sphere with at most one (unramified) irregular singularity and arbitrary number of simple poles as Nakajima's quiver varieties. This result enables us to solve partially the additive irregular Deligne–Simpson problem.     

A pathwise solution for nonlinear parabolic equations with stochastic perturbations

We analyse here a semilinear stochastic partial differential equation of parabolic type where the diffusion vector fields are depending on both the unknown function and its gradient ∂ _xu with respect to the state variable, ∈ ℝⁿ. A local solution is constructed by reducing the original equation to a nonlinear parabolic one without stochastic perturbations and it is based on a finite dimensional Lie algebra generated by the given diffusion vector fields.     

The uniqueness of a solution to the renewal type system of integral equations on the line

We study the uniqueness of a solution to a renewal type system of integral equations z=g+F * z on the line ℝ; here z is the unknown vector function, g is a known vector function, and F is a nonlattice matrix of finite measures on ℝ such that the matrix F(ℝ) is of spectral radius 1 and indecomposable. We show that in a certain class of functions each solution to the corresponding homogeneous system coincides almost everywhere with a right eigenvector of F(ℝ) with eigenvalue 1.     

Convex mean curvature flow with a forcing term in direction of the position vector

A smooth, compact and strictly convex hypersurface evolving in ℝ ⁿ⁺¹ along its mean curvature vector plus a forcing term in the direction of its position vector is studied in this paper. We show that the convexity is preserving as the case of mean curvature flow, and the evolving convex hypersurfaces may shrink to a point in finite time if the forcing term is small, or exist for all time and expand to infinity if it is large enough. The flow can converge to a round sphere if the forcing term satisfies suitable conditions which will be given in the paper. Long-time existence and convergence of normalization of the flow are also investigated.     

Solving the Laplace-Beltrami equation on S2 using spherical triangles

The geometry of the domain −S², causes difficulty in solving the Laplace-Beltrami Equation, for example, in discretization for the differential equation. To overcome this problem, we study a numerical method, which is based on the finite element approximation with a hierarchical refinement of icosahedron for the grid. We construct a geometrically intrinsic base vector field for the Galerkin approximation. In this way, no artificial poles are introduced, and the numerical grids are distributed more evenly. We use radial projection to map the curved triangle onto a flat one, so that existing quadrature schemes can be applied for the numerical integration. The resulting system of linear algebraic equations is solved by using a conjugate gradient method. © 1996 John Wiley & Sons, Inc.     

On splitting schemes in the mixed finite element method

Within the research into some geothermal modes, a 3D heat transfer process was described by a first-order system of differential equations (in terms of “temperature-heat-flow”). This system was solved by an explicit scheme for the mixed finite element spatial approximations based on the Raviart-Thomas degrees of freedom. In this paper, several algorithms based on the splitting technique for the vector heat-flow equation are proposed. Some comparison results of accuracy of the algorithms proposed are presented.     

Pattern recognition with the help of quadratic discriminant functions

The problem of pattern recognition with the help of spherical and elliptic discriminant functions is studied; in so doing the pattern of an object is assumed to be a vector of its characters from a finite-dimensional Euclidean space. Using a conformal mapping of a punctured sphere onto the plane as well as the inversion transformation, a criterion for the error-free recognition of two sets containing a finite number of points of training samples is obtained with the help of spherical discriminant functions. An algorithm for solving approximately a problem of construction of an ellipsoid of the minimal volume containing a given finite set of points is described. Bibliography: 5 titles. Translated fromObchyslyuval'na ta Prykladna Matematyka, No. 80, 1996, pp. 90–105.     

The distribution of the poles of the best approximating rational functions and the analytical properties of the approximated function

The following theorem was proved in our paper in Math. USSR-Sb.9 (1969): Iff(z)∈H ₂ and the poles of the rational functions of best approximation tend to infinity sufficiently quickly, thenf(z) is an entire function. In the present article we weaken the restrictions on the distribution of poles by assuming only that these poles have no finite limit point (Theorem 1). Some generalizations of this result are also given.     

Wavelet Decomposition of Spherical Vector Fields with Respect to Sources

This article is concerned with an approach of modelling the Earth’s magnetic field as measured by satellites in terms of a special system of vector spherical harmonics and in terms of vector kernel functions, called vector scaling functions and wavelets. The main ingredient is the presentation of a system of vector spherical harmonics which separates a given spherical vector field with respect to its sources, i.e., the spherical vector field is separated into a part which is induced by sources inside the sphere, a part which is induced by sources outside the sphere and a part which is induced by sources on the sphere, which are, for example, currents crossing the sphere. Using this special system of vector spherical harmonics vector scaling functions and wavelets are constructed which keep the advantageous property of separating with respect to sources but which also allow a locally reflected modelling of the respective vector field. At the end of the article, the method is tested on real magnetic field data measured by the German geoscientific research satellite CHAMP.     

Spatial discretization of partial differential equations with integrals

We consider the problem of constructing spatial finite-differenceapproximations on an arbitrary fixed grid which preserve anynumber of integrals of the partial differential equation andpreserve some of its symmetries. A basis for the space of suchfinite-difference operators is constructed; most cases of interestinvolve a single such basis element. (The ‘Arakawa’Jacobian is such an element, as are discretizations satisfying‘summation by parts’ identities.) We show how thegrid, its symmetries, and the differential operator interactto affect the complexity of the finite difference.     

Clement-type interpolation on spherical domains--interpolation error estimates and application to a posteriori error estimation

In this paper, a mixed boundary value problem for the Laplace–Beltramioperator is considered for spherical domains in R³, i.e. fordomains on the unit sphere. These domains are parametrized byspherical coordinates (, ). A suitable finite element spaceis introduced. It corresponds to an isotropic triangulationof the underlying domain on the unit sphere. Error estimatesare proven for a Clément-type interpolation operator,where appropriate weighted norms are used. The estimates areapplied to the derivation of a reliable and efficient residualerror estimator for the Laplace–Beltrami operator.     

Semigroups containing proximal linear maps

A linear automorphism of a finite dimensional real vector spaceV is calledproximal if it has a unique eigenvalue—counting multiplicities—of maximal modulus. Goldsheid and Margulis have shown that if a subgroupG of GL(V) contains a proximal element then so does every Zariski dense subsemigroupH ofG, providedV considered as aG-module is strongly irreducible. We here show thatH contains a finite subsetM such that for everyg∈GL(V) at least one of the elements γg, γ∈M, is proximal. We also give extensions and refinements of this result in the following directions: a quantitative version of proximality, reducible representations, several eigenvalues of maximal modulus. Partially supported by NSF grant DMS 9204-720.     

A fixed-center spherical separation algorithm with kernel transformations for classification problems

We consider a special case of the optimal separation, via a sphere, of two discrete point sets in a finite dimensional Euclidean space. In fact we assume that the center of the sphere is fixed. In this case the problem reduces to the minimization of a convex and nonsmooth function of just one variable, which can be solved by means of an “ad hoc” method in O(p log p) time, where p is the dataset size. The approach is suitable for use in connection with kernel transformations of the type adopted in the support vector machine (SVM) approach. Despite of its simplicity the method has provided interesting results on several standard test problems drawn from the binary classification literature. This research has been partially supported by the Italian “Ministero dell’Istruzione, dell’Università e della Ricerca Scientifica”, under PRIN project Numerical Methods for Global Optimization and for some classes of Nonsmooth Optimization Problems (2005017083.002).