4,510 research outputs found
Segregated Runge–Kutta time integration of convection-stabilized mixed finite element schemes for wall-unresolved LES of incompressible flows
In this work, we develop a high-performance numerical framework for the large eddy simulation (LES) of incompressible flows. The spatial discretization of the nonlinear system is carried out using mixed finite element (FE) schemes supplemented with symmetric projection stabilization of the convective term and a penalty term for the divergence constraint. These additional terms introduced at the discrete level have been proved to act as implicit LES models. In order to perform meaningful wall-unresolved simulations, we consider a weak imposition of the boundary conditions using a Nitsche’s-type scheme, where the tangential component penalty term is designed to act as a wall law. Next, segregated Runge–Kutta (SRK) schemes (recently proposed by the authors for laminar flow problems) are applied to the LES simulation of turbulent flows. By the introduction of a penalty term on the trace of the acceleration, these methods exhibit excellent stability properties for both implicit and explicit treatment of the convective terms. SRK schemes are excellent for large-scale simulations, since they reduce the computational cost of the linear system solves by splitting velocity and pressure computations at the time integration level, leading to two uncoupled systems. The pressure system is a Darcy-type problem that can easily be preconditioned using a traditional block-preconditioning scheme that only requires a Poisson solver. At the end, only coercive systems have to be solved, which can be effectively preconditioned by multilevel domain decomposition schemes, which are both optimal and scalable. The framework is applied to the Taylor–Green and turbulent channel flow benchmarks in order to prove the accuracy of the convection-stabilized mixed FEs as LES models and SRK time integrators. The scalability of the preconditioning techniques (in space only) has also been proven for one step of the SRK scheme for the Taylor–Green flow using uniform meshes. Moreover, a turbulent flow around a NACA profile is solved to show the applicability of the proposed algorithms for a realistic problem.Peer ReviewedPostprint (author's final draft
Dual virtual element method for discrete fractures networks
Discrete fracture networks is a key ingredient in the simulation of physical
processes which involve fluid flow in the underground, when the surrounding
rock matrix is considered impervious. In this paper we present two different
models to compute the pressure field and Darcy velocity in the system. The
first allows a normal flow out of a fracture at the intersections, while the
second grants also a tangential flow along the intersections. For the numerical
discretization, we use the mixed virtual finite element method as it is known
to handle grid elements of, almost, any arbitrary shape. The flexibility of the
discretization allows us to loosen the requirements on grid construction, and
thus significantly simplify the flow discretization compared to traditional
discrete fracture network models. A coarsening algorithm, from the algebraic
multigrid literature, is also considered to further speed up the computation.
The performance of the method is validated by numerical experiments
The Grid Dependence of Well Inflow Performance in Reservoir Simulation
Imperial Users onl
Natural preconditioners for saddle point systems
The solution of quadratic or locally quadratic extremum problems subject to linear(ized) constraints gives rise to linear systems in saddle point form. This is true whether in the continuous or discrete setting, so saddle point systems arising from discretization of partial differential equation problems such as those describing electromagnetic problems or incompressible flow lead to equations with this structure as does, for example, the widely used sequential quadratic programming approach to nonlinear optimization.\ud
This article concerns iterative solution methods for these problems and in particular shows how the problem formulation leads to natural preconditioners which guarantee rapid convergence of the relevant iterative methods. These preconditioners are related to the original extremum problem and their effectiveness -- in terms of rapidity of convergence -- is established here via a proof of general bounds on the eigenvalues of the preconditioned saddle point matrix on which iteration convergence depends
Continuous, Semi-discrete, and Fully Discretized Navier-Stokes Equations
The Navier--Stokes equations are commonly used to model and to simulate flow
phenomena. We introduce the basic equations and discuss the standard methods
for the spatial and temporal discretization. We analyse the semi-discrete
equations -- a semi-explicit nonlinear DAE -- in terms of the strangeness index
and quantify the numerical difficulties in the fully discrete schemes, that are
induced by the strangeness of the system. By analyzing the Kronecker index of
the difference-algebraic equations, that represent commonly and successfully
used time stepping schemes for the Navier--Stokes equations, we show that those
time-integration schemes factually remove the strangeness. The theoretical
considerations are backed and illustrated by numerical examples.Comment: 28 pages, 2 figure, code available under DOI: 10.5281/zenodo.998909,
https://doi.org/10.5281/zenodo.99890
- …