129 research outputs found
Comparison of different nonlinear solvers for 2D time-implicit stellar hydrodynamics
Time-implicit schemes are attractive since they allow numerical time steps
that are much larger than those permitted by the Courant-Friedrich-Lewy
criterion characterizing time-explicit methods. This advantage comes, however,
with a cost: the solution of a system of nonlinear equations is required at
each time step. In this work, the nonlinear system results from the
discretization of the hydrodynamical equations with the Crank-Nicholson scheme.
We compare the cost of different methods, based on Newton-Raphson iterations,
to solve this nonlinear system, and benchmark their performances against
time-explicit schemes. Since our general scientific objective is to model
stellar interiors, we use as test cases two realistic models for the convective
envelope of a red giant and a young Sun. Focusing on 2D simulations, we show
that the best performances are obtained with the quasi-Newton method proposed
by Broyden. Another important concern is the accuracy of implicit calculations.
Based on the study of an idealized problem, namely the advection of a single
vortex by a uniform flow, we show that there are two aspects: i) the nonlinear
solver has to be accurate enough to resolve the truncation error of the
numerical discretization, and ii) the time step has be small enough to resolve
the advection of eddies. We show that with these two conditions fulfilled, our
implicit methods exhibit similar accuracy to time-explicit schemes, which have
lower values for the time step and higher computational costs. Finally, we
discuss in the conclusion the applicability of these methods to fully implicit
3D calculations.Comment: Accepted for publication in A&
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
Preconditioning the Advection-Diffusion Equation: the Green's Function Approach
We look at the relationship between efficient preconditioners (i.e., good approximations to the discrete inverse operator) and the generalized inverse for the (continuous) advection-diffusion operator -- the Green's function. We find that the continuous Green's function exhibits two important properties -- directionality and rapid downwind decay -- which are preserved by the discrete (grid) Green's functions, if and only if the discretization used produces non-oscillatory solutions. In particular, the downwind decay ensures the locality of the grid Green's functions. Hence, a finite element formulation which produces a good solution will typically use a coefficient matrix with almost lower triangular structure under a "with-the-flow" numbering of the variables. It follows that the block Gauss-Seidel matrix is a first candidate for a preconditioner to use with an iterative solver of Krylov subspace type
- …