Non-nested multi-grid solvers for mixed divergence-free Scott-Vogelius discretizations

Studying high-dimensional Hamiltonian systems with microstructure, it is an important and challenging problem to identify reduced macroscopic models that describe some effective dynamics on large spatial and temporal scales. This paper concerns the question how reasonable macroscopic Lagrangian and Hamiltonian structures can by derived from the microscopic system. In the first part we develop a general approach to this problem by considering non-canonical Hamiltonian structures on the tangent bundle. This approach can be applied to all Hamiltonian lattices (or Hamiltonian PDEs) and involves three building blocks: (i) the embedding of the microscopic system, (ii) an invertible two-scale transformation that encodes the underlying scaling of space and time, (iii) an elementary model reduction that is based on a Principle of Consistent Expansions. In the second part we exemplify the reduction approach and derive various reduced PDE models for the atomic chain. The reduced equations are either related to long wave-length motion or describe the macroscopic modulation of an oscillatory microstructure

Numerical Study of Diffusion of Interacting Particles in a Magnetic Fluid Layer

The paper is devoted to construction of mathematical model and computational algorithm for determining equilibrium shapes of a free surface of a magnetic fluid subject to diffusion of ferromagnetic particles in the fluid. The mathematical model consists of the Maxwell equations to determine the structure of magnetic field, the mass transfer equation for the particles, and the Young-Laplace equations for the fluid-air interface as well. Models of the uniform distribution of the particles, the nonuniform distribution without considering the particle interaction, and the nonuniform distribution subject to the particle interaction are considered. Numerical comparison of the models is carried out on the ferrohydrostatic problem of stability of a magnetic-fluid layer under a uniform magnetic fiel

Robust arbitrary order mixed finite element methods for the incompressible Stokes equations

Standard mixed finite element methods for the incompressible Navier-Stokes equations that relax the divergence constraint are not robust against large irrotational forces in the momentum balance and the velocity error depends on the continuous pressure. This robustness issue can be completely cured by using divergence-free mixed finite elements which deliver pressure-independent velocity error estimates. However, the construction of H1-conforming, divergence-free mixed finite element methods is rather difficult. Instead, we present a novel approach for the construction of arbitrary order mixed finite element methods which deliver pressure-independent velocity errors. The approach does not change the trial functions but replaces discretely divergence-free test functions in some operators of the weak formulation by divergence-free ones. This modification is applied to inf-sup stable conforming and nonconforming mixed finite element methods of arbitrary order in two and three dimensions. Optimal estimates for the incompressible Stokes equations are proved for the H1 and L2 errors of the velocity and the L2 error of the pressure. Moreover, both velocity errors are pressure-independent, demonstrating the improved robustness. Several numerical examples illustrate the results

Robust arbitrary order mixed finite element methods for the incompressible Stokes equations

Standard mixed finite element methods for the incompressible Navier-Stokes equations that relax the divergence constraint are not robust against large irrotational forces in the momentum balance and the velocity error depends on the continuous pressure. This robustness issue can be completely cured by using divergence-free mixed finite elements which deliver pressure-independent velocity error estimates. However, the construction of H1-conforming, divergence-free mixed finite element methods is rather difficult. Instead, we present a novel approach for the construction of arbitrary order mixed finite element methods which deliver pressure-independent velocity errors. The approach does not change the trial functions but replaces discretely divergence-free test functions in some operators of the weak formulation by divergence-free ones. This modification is applied to inf-sup stable conforming and nonconforming mixed finite element methods of arbitrary order in two and three dimensions. Optimal estimates for the incompressible Stokes equations are proved for the H1 and L2 errors of the velocity and the L2 error of the pressure. Moreover, both velocity errors are pressure-independent, demonstrating the improved robustness. Several numerical examples illustrate the results

The Surface Topography of a Magnetic Fluid -- a Quantitative Comparison between Experiment and Numerical Simulation

The normal field instability in magnetic liquids is investigated experimentally by means of a radioscopic technique which allows a precise measurement of the surface topography. The dependence of the topography on the magnetic field is compared to results obtained by numerical simulations via the finite element method. Quantitative agreement has been found for the critical field of the instability, the scaling of the pattern amplitude and the detailed shape of the magnetic spikes. The fundamental Fourier mode approximates the shape to within 10% accuracy for a range of up to 40% of the bifurcation parameter of this subcritical bifurcation. The measured control parameter dependence of the wavenumber differs qualitatively from analytical predictions obtained by minimization of the free energy.Comment: 21 pages, 16 figures; corrected typos, added reference to Kuznetsov and Spector(1976), S.J. Fortune(1995) and Harkins&Jordan (1930). Figures revise

Analysis of a mathematical model related to Czochralski crystal growth

This paper is devoted to the study of a stationary problem consisting of the Boussinesq approximation of the Navier–Stokes equations and two convection–diffusion equations for the temperature and concentration, respectively. The equations are considered in 3D and a velocity–pressure formulation of the Navier–Stokes equations is used. The problem is complicated by nonstandard boundary conditions for velocity on the liquid–gas interface where tangential surface forces proportional to surface gradients of temperature and concentration (Marangoni effect) and zero normal component of the velocity are assumed. The velocity field is coupled through this boundary condition and through the buoyancy term in the Navier–Stokes equations with both the temperature and concentration fields. In this paper a weak formulation of the problem is stated and the existence of a weak solution is proved. For small data, the uniqueness of the solution is established