Abstract

We examine the use of natural boundary conditions and conditions of the Sommerfeld type for finite element simulations of convective transport in viscous incompressible flows. We show that natural boundary conditions are superior in the sense that they always provide a correct boundary condition, as opposed to the Sommerfeld‐type conditions, which can lead to a singular formulation and a great loss of accuracy. For the Navier–Stokes equations, the natural boundary conditions must be combined with a simple method to eliminate perturbations on the pressure at the open boundary, which is the source of most errors.

Abstract

The finite point method (FPM) is a gridless numerical procedure based on the combination of weighted least square interpolations on a cloud of points with point collocation for evaluating the approximation integrals. In the paper, details of a procedure for stabilizing the numerical solution for advective-diffusive transport and fluid flow problems using the FPM are given. The method is based on a consistent introduction of the stabilizing terms in the governing differential equations. One example showing the applicability of the FPM is given.

Abstract

We propose a novel approach for flow simulations in fractured porous media based on an equi-dimensional representation of the fractures and a standard continuous finite element (FE) method. We employ an adaptive mesh refinement strategy to automatically adjust an initial regular mesh to any fracture distribution with a desired accuracy. The proposed approach is easily implementable in any FE software, does not involve the coupling of different discretizations, and provides symmetric and positive definite stiffness matrices. We provide a systematic validation of our method and we show that it can be used to simulate fluid flow for fracture networks of realistic complexity.