Flow of evaporating, gravity-driven thin liquid films over topography

The effect of topography on the free surface and solvent concentration profiles of an evaporating thin film of liquid flowing down an inclined plane is considered. The liquid is assumed to be composed of a resin dissolved in a volatile solvent with the associated solvent concentration equation derived on the basis of the well-mixed approximation. The dynamics of the film is formulated as a lubrication approximation and the effect of a composition-dependent viscosity is included in the model. The resulting time-dependent, nonlinear, coupled set of governing equations is solved using a full approximation storage multigrid method. The approach is first validated against a closed-form analytical solution for the case of a gravity-driven, evaporating thin film flowing down a flat substrate. Analysis of the results for a range of topography shapes reveal that although a full-width, spanwise topography such as a step-up or a step-down does not affect the composition of the film, the same is no longer true for the case of localized topography, such as a peak or a trough, for which clear nonuniformities of the solvent concentration profile can be observed in the wake of the topography

A Fast Parallel Poisson Solver on Irregular Domains Applied to Beam Dynamic Simulations

We discuss the scalable parallel solution of the Poisson equation within a Particle-In-Cell (PIC) code for the simulation of electron beams in particle accelerators of irregular shape. The problem is discretized by Finite Differences. Depending on the treatment of the Dirichlet boundary the resulting system of equations is symmetric or `mildly' nonsymmetric positive definite. In all cases, the system is solved by the preconditioned conjugate gradient algorithm with smoothed aggregation (SA) based algebraic multigrid (AMG) preconditioning. We investigate variants of the implementation of SA-AMG that lead to considerable improvements in the execution times. We demonstrate good scalability of the solver on distributed memory parallel processor with up to 2048 processors. We also compare our SAAMG-PCG solver with an FFT-based solver that is more commonly used for applications in beam dynamics

The evolution of electron overdensities in magnetic fields

When a neutral gas impinges on a stationary magnetized plasma an enhancement in the ionization rate occurs when the neutrals exceed a threshold velocity. This is commonly known as the critical ionization velocity effect. This process has two distinct timescales: an ion–neutral collision time and electron acceleration time. We investigate the energization of an ensemble of electrons by their self-electric field in an applied magnetic field. The evolution of the electrons is simulated under different magnetic field and density conditions. It is found that electrons can be accelerated to speeds capable of electron impact ionization for certain conditions. In the magnetically dominated case the energy distribution of the excited electrons shows that typically 1% of the electron population can exceed the initial electrostatic potential associated with the unbalanced ensemble of electrons

Multigrid optimization for space-time discontinuous Galerkin discretizations of advection dominated flows

The goal of this research is to optimize multigrid methods for higher order accurate space-time discontinuous Galerkin discretizations. The main analysis tool is discrete Fourier analysis of two- and three-level multigrid algorithms. This gives the spectral radius of the error transformation operator which predicts the asymptotic rate of convergence of the multigrid algorithm. In the optimization process we therefore choose to minimize the spectral radius of the error transformation operator. We specifically consider optimizing h-multigrid methods with explicit Runge-Kutta type smoothers for second and third order accurate space-time discontinuous Galerkin finite element discretizations of the 2D advection-diffusion equation. The optimized schemes are compared with current h-multigrid techniques employing Runge-Kutta type smoothers. Also, the efficiency of h-, p- and hp-multigrid methods for solving the Euler equations of gas dynamics with a higher order accurate space-time DG method is investigated

Particle-Particle, Particle-Scaling function (P3S) algorithm for electrostatic problems in free boundary conditions

An algorithm for fast calculation of the Coulombic forces and energies of point particles with free boundary conditions is proposed. Its calculation time scales as N log N for N particles. This novel method has lower crossover point with the full O(N^2) direct summation than the Fast Multipole Method. The forces obtained by our algorithm are analytical derivatives of the energy which guarantees energy conservation during a molecular dynamics simulation. Our algorithm is very simple. An MPI parallelised version of the code can be downloaded under the GNU General Public License from the website of our group.Comment: 19 pages, 11 figures, submitted to: Journal of Chemical Physic

Local Fourier Analysis of the Complex Shifted Laplacian preconditioner for Helmholtz problems

In this paper we solve the Helmholtz equation with multigrid preconditioned Krylov subspace methods. The class of Shifted Laplacian preconditioners are known to significantly speed-up Krylov convergence. However, these preconditioners have a parameter beta, a measure of the complex shift. Due to contradictory requirements for the multigrid and Krylov convergence, the choice of this shift parameter can be a bottleneck in applying the method. In this paper, we propose a wavenumber-dependent minimal complex shift parameter which is predicted by a rigorous k-grid Local Fourier Analysis (LFA) of the multigrid scheme. We claim that, given any (regionally constant) wavenumber, this minimal complex shift parameter provides the reader with a parameter choice that leads to efficient Krylov convergence. Numerical experiments in one and two spatial dimensions validate the theoretical results. It appears that the proposed complex shift is both the minimal requirement for a multigrid V-cycle to converge, as well as being near-optimal in terms of Krylov iteration count.Comment: 20 page

Quantitative phase-field modeling of solidification at high Lewis number

A phase-field model of nonisothermal solidification in dilute binary alloys is used to study the variation of growth velocity, dendrite tip radius, and radius selection parameter as a function of Lewis number at fixed undercooling. By the application of advanced numerical techniques, we have been able to extend the analysis to Lewis numbers of order 10 000, which are realistic for metals. A large variation in the radius selection parameter is found as the Lewis number is increased from 1 to 10 000

Spontaneous deterministic side-branching behavior in phase-field simulations of equiaxed dendritic growth

The accepted view on dendritic side-branching is that side-branches grow as the result of selective amplification of thermal noise and that in the absence of such noise dendrites would grow without the development of side-arms. However, recently there has been renewed speculation about dendrites displaying deterministic side-branching [see, e.g., M. E. Glicksman, Metall. Mater. Trans A 43, 391 (2012)]. Generally, numerical models of dendritic growth, such as phase-field simulation, have tended to display behaviour which is commensurate with the former view, in that simulated dendrites do not develop side-branches unless noise is introduced into the simulation. However, here, we present simulations that show that under certain conditions deterministic side-branching may occur. We use a model formulated in the thin interface limit and a range of advanced numerical techniques to minimise the numerical noise introduced into the solution, including a multigrid solver. Spontaneous side-branching seems to be favoured by high undercoolings and by intermediate values of the capillary anisotropy, with the most branched structures being obtained for an anisotropy strength of 0.03. From an analysis of the tangential thermal gradients on the solid-liquid interface, the mechanism for side-branching appears to have some similarities with the deterministic model proposed by Glicksman

On the indefinite Helmholtz equation: complex stretched absorbing boundary layers, iterative analysis, and preconditioning

This paper studies and analyzes a preconditioned Krylov solver for Helmholtz problems that are formulated with absorbing boundary layers based on complex coordinate stretching. The preconditioner problem is a Helmholtz problem where not only the coordinates in the absorbing layer have an imaginary part, but also the coordinates in the interior region. This results into a preconditioner problem that is invertible with a multigrid cycle. We give a numerical analysis based on the eigenvalues and evaluate the performance with several numerical experiments. The method is an alternative to the complex shifted Laplacian and it gives a comparable performance for the studied model problems

Multigrid elliptic equation solver with adaptive mesh refinement

In this paper we describe in detail the computational algorithm used by our parallel multigrid elliptic equation solver with adaptive mesh refinement. Our code uses truncation error estimates to adaptively refine the grid as part of the solution process. The presentation includes a discussion of the orders of accuracy that we use for prolongation and restriction operators to ensure second order accurate results and to minimize computational work. Code tests are presented that confirm the overall second order accuracy and demonstrate the savings in computational resources provided by adaptive mesh refinement.Comment: 12 pages, 9 figures, Modified in response to reviewer suggestions, added figure, added references. Accepted for publication in J. Comp. Phy