158 research outputs found
A constrained pressure-temperature residual (CPTR) method for non-isothermal multiphase flow in porous media
For both isothermal and thermal petroleum reservoir simulation, the
Constrained Pressure Residual (CPR) method is the industry-standard
preconditioner. This method is a two-stage process involving the solution of a
restricted pressure system. While initially designed for the isothermal case,
CPR is also the standard for thermal cases. However, its treatment of the
energy conservation equation does not incorporate heat diffusion, which is
often dominant in thermal cases. In this paper, we present an extension of CPR:
the Constrained Pressure-Temperature Residual (CPTR) method, where a restricted
pressure-temperature system is solved in the first stage. In previous work, we
introduced a block preconditioner with an efficient Schur complement
approximation for a pressure-temperature system. Here, we extend this method
for multiphase flow as the first stage of CPTR. The algorithmic performance of
different two-stage preconditioners is evaluated for reservoir simulation test
cases.Comment: 28 pages, 2 figures. Sources/sinks description in arXiv:1902.0009
Constraint interface preconditioning for topology optimization problems
The discretization of constrained nonlinear optimization problems arising in
the field of topology optimization yields algebraic systems which are
challenging to solve in practice, due to pathological ill-conditioning, strong
nonlinearity and size. In this work we propose a methodology which brings
together existing fast algorithms, namely, interior-point for the optimization
problem and a novel substructuring domain decomposition method for the ensuing
large-scale linear systems. The main contribution is the choice of interface
preconditioner which allows for the acceleration of the domain decomposition
method, leading to performance independent of problem size.Comment: To be published in SIAM J. Sci. Com
The LifeV library: engineering mathematics beyond the proof of concept
LifeV is a library for the finite element (FE) solution of partial
differential equations in one, two, and three dimensions. It is written in C++
and designed to run on diverse parallel architectures, including cloud and high
performance computing facilities. In spite of its academic research nature,
meaning a library for the development and testing of new methods, one
distinguishing feature of LifeV is its use on real world problems and it is
intended to provide a tool for many engineering applications. It has been
actually used in computational hemodynamics, including cardiac mechanics and
fluid-structure interaction problems, in porous media, ice sheets dynamics for
both forward and inverse problems. In this paper we give a short overview of
the features of LifeV and its coding paradigms on simple problems. The main
focus is on the parallel environment which is mainly driven by domain
decomposition methods and based on external libraries such as MPI, the Trilinos
project, HDF5 and ParMetis.
Dedicated to the memory of Fausto Saleri.Comment: Review of the LifeV Finite Element librar
- …