23,638 research outputs found
An adaptive preconditioner for steady incompressible flows
This paper describes an adaptive preconditioner for numerical continuation of
incompressible Navier--Stokes flows. The preconditioner maps the identity (no
preconditioner) to the Stokes preconditioner (preconditioning by Laplacian)
through a continuous parameter and is built on a first order Euler
time-discretization scheme. The preconditioner is tested onto two fluid
configurations: three-dimensional doubly diffusive convection and a reduced
model of shear flows. In the former case, Stokes preconditioning works but a
mixed preconditioner is preferred. In the latter case, the system of equation
is split and solved simultaneously using two different preconditioners, one of
which is parameter dependent. Due to the nature of these applications, this
preconditioner is expected to help a wide range of studies
Gr\"obner Bases and Generation of Difference Schemes for Partial Differential Equations
In this paper we present an algorithmic approach to the generation of fully
conservative difference schemes for linear partial differential equations. The
approach is based on enlargement of the equations in their integral
conservation law form by extra integral relations between unknown functions and
their derivatives, and on discretization of the obtained system. The structure
of the discrete system depends on numerical approximation methods for the
integrals occurring in the enlarged system. As a result of the discretization,
a system of linear polynomial difference equations is derived for the unknown
functions and their partial derivatives. A difference scheme is constructed by
elimination of all the partial derivatives. The elimination can be achieved by
selecting a proper elimination ranking and by computing a Gr\"obner basis of
the linear difference ideal generated by the polynomials in the discrete
system. For these purposes we use the difference form of Janet-like Gr\"obner
bases and their implementation in Maple. As illustration of the described
methods and algorithms, we construct a number of difference schemes for Burgers
and Falkowich-Karman equations and discuss their numerical properties.Comment: Published in SIGMA (Symmetry, Integrability and Geometry: Methods and
Applications) at http://www.emis.de/journals/SIGMA
An adaptive fixed-mesh ALE method for free surface flows
In this work we present a Fixed-Mesh ALE method for the numerical simulation of free surface flows capable of using an adaptive finite element mesh covering a background domain. This mesh is successively refined and unrefined at each time step in order to focus the computational effort on the spatial regions where it is required. Some of the main ingredients of the formulation are the use of an Arbitrary-Lagrangian–Eulerian formulation for computing temporal derivatives, the use of stabilization terms for stabilizing convection, stabilizing the lack of compatibility between velocity and pressure interpolation spaces, and stabilizing the ill-conditioning introduced by the cuts on the background finite element mesh, and the coupling of the algorithm with an adaptive mesh refinement procedure suitable for running on distributed memory environments. Algorithmic steps for the projection between meshes are presented together with the algebraic fractional step approach used for improving the condition number of the linear systems to be solved. The method is tested in several numerical examples. The expected convergence rates both in space and time are observed. Smooth solution fields for both velocity and pressure are obtained (as a result of the contribution of the stabilization terms). Finally, a good agreement between the numerical results and the reference experimental data is obtained.Postprint (published version
- …