2,320 research outputs found
On acceleration of Krylov-subspace-based Newton and Arnoldi iterations for incompressible CFD: replacing time steppers and generation of initial guess
We propose two techniques aimed at improving the convergence rate of steady
state and eigenvalue solvers preconditioned by the inverse Stokes operator and
realized via time-stepping. First, we suggest a generalization of the Stokes
operator so that the resulting preconditioner operator depends on several
parameters and whose action preserves zero divergence and boundary conditions.
The parameters can be tuned for each problem to speed up the convergence of a
Krylov-subspace-based linear algebra solver. This operator can be inverted by
the Uzawa-like algorithm, and does not need a time-stepping. Second, we propose
to generate an initial guess of steady flow, leading eigenvalue and eigenvector
using orthogonal projection on a divergence-free basis satisfying all boundary
conditions. The approach, including the two proposed techniques, is illustrated
on the solution of the linear stability problem for laterally heated square and
cubic cavities
Time integration and steady-state continuation for 2d lubrication equations
Lubrication equations allow to describe many structurin processes of thin
liquid films. We develop and apply numerical tools suitable for their analysis
employing a dynamical systems approach. In particular, we present a time
integration algorithm based on exponential propagation and an algorithm for
steady-state continuation. In both algorithms a Cayley transform is employed to
overcome numerical problems resulting from scale separation in space and time.
An adaptive time-step allows to study the dynamics close to hetero- or
homoclinic connections. The developed framework is employed on the one hand to
analyse different phases of the dewetting of a liquid film on a horizontal
homogeneous substrate. On the other hand, we consider the depinning of drops
pinned by a wettability defect. Time-stepping and path-following are used in
both cases to analyse steady-state solutions and their bifurcations as well as
dynamic processes on short and long time-scales. Both examples are treated for
two- and three-dimensional physical settings and prove that the developed
algorithms are reliable and efficient for 1d and 2d lubrication equations,
respectively.Comment: 33 pages, 16 figure
Efficient preconditioning of the method of lines for solving nonlinear two-sided space-fractional diffusion equations
A standard method for the numerical solution of partial differential equations (PDEs) is the method of lines. In this approach the PDE is discretised in space using �finite di�fferences or similar techniques, and the resulting semidiscrete problem in time is integrated using an initial value problem solver. A significant challenge when applying the method of lines to fractional PDEs is that the non-local nature of the fractional derivatives results in a discretised system where each equation involves contributions from many (possibly every) spatial node(s). This has important consequences for the effi�ciency of the numerical solver. First, since the cost of evaluating the discrete equations is high, it is essential to minimise the number of evaluations required to advance the solution in time. Second, since the Jacobian matrix of the system is dense (partially or fully), methods that avoid the need to form and factorise this matrix are preferred. In this paper, we consider a nonlinear two-sided space-fractional di�ffusion equation in one spatial dimension. A key contribution of this paper is to demonstrate how an eff�ective preconditioner is crucial for improving the effi�ciency of the method of lines for solving this equation. In particular, we show how to construct suitable banded approximations to the system Jacobian for preconditioning purposes that permit high orders and large stepsizes to be used in the temporal integration, without requiring dense matrices to be formed. The results of numerical experiments are presented that demonstrate the effectiveness of this approach
- …