In this paper, we present a novel semi-analytical boundary handling method for spatially adaptive and divergence-free smoothed particle hydrodynamics (SPH) simulations, including two-way coupling. Our method is consistent under varying particle resolutions and allows for the treatment of boundary features below the particle resolution. We achieve this by first introducing an analytic solution to the interaction of SPH particles with planar boundaries, in 2D and 3D, which we extend to arbitrary boundary geometries using signed distance fields (SDF) to construct locally planar boundaries. Using this boundary-integral-based approach, we can directly evaluate boundary contributions, for any quantity, allowing an easy integration into state of the art simulation methods. Overall, our method improves interactions with small boundary features, readily handles spatially adaptive fluids, preserves particle-boundary interactions across varying resolutions, can directly be implemented in existing SPH methods, and, for non-adaptive simulations, provides a reduction in memory consumption as well as an up to 2× speedup relative to current particle-based boundary handling approaches.

To analyze the waterflooding system of an oil reservoir and predict the effectiveness of geological and technical measures, information is required on the distribution of injection rate between the reacting production wells and the reservoir boundary. The most reliable methods for calculating these characteristics are methods based on hydrodynamic modeling of flow. Modern commercial software implement algorithms for these purposes based on the construction and analysis of streamlines. At the same time, there are no reliable estimates of the accuracy of these algorithms and recommendations for choosing the optimal parameters in the available literature. In this paper, we propose an algorithm for calculating the proportions of the distribution of the total well flow rate between the surrounding wells and the reservoir boundary using streamlines. Streamlines are constructed on the basis of a finite element solution to the flow problem averaged over the formation thickness and determine the boundaries of the streamtubes connecting the corresponding wells. The flow rate through the flow tubes is calculated by numerically integrating the Darcy velocity field of the indicated two-dimensional problem. The algorithm was tested on idealized examples of waterflooding elements of typical well placement schemes, when the exact distribution of the proportions of fluid injected into the formation is known, and on the example of comparison with the solution of the problem of simulating the injection of a tracer into the reservoir. Recommendations for the selection of starting points for tracing streamlines are presented, which allow achieving a minimum level of error in determining the mutual influence of wells in a wide range of the computational grid resolution of the flow model. A more general application of the described method without significant changes is to equip the high resolution flow model along fixed stream tubes with their rate characteristics.

Boundary handling is an important aspect of fluid simulation, and several boundary handling approaches exist in smoothed particle hydrodynamics (SPH), which have individual strengths and weaknesses. However, comparing different boundary handling approaches is challenging as there is no common basis for evaluations, that is, no universal set of experiments with quantitative evaluation across different methods, especially within computer animation where many evaluations rely mainly on visual perception. This article proposes a set of experiments to aid the evaluation of the main categories of fluid‐boundary interactions that are important in computer animation, that is, no motion (resting) fluid, tangential and normal motion of a fluid with respect to the boundary, and a fluid impacting a corner. We propose ten experiments, comprising experimental setup and quantitative evaluation with optional visual inspections, that are arranged in four groups which focus on one of the main category of fluid‐boundary interactions. We use these experiments to evaluate three particle‐based boundary handling methods, that is, pressure mirroring, pressure boundaries, and moving least squares pressure extrapolation, in combination with two incompressible SPH fluid simulation methods, namely IISPH and DFSPH, to establish a quantifiable relation between different combinations of boundary handling with simulation approaches and the main categories of fluid‐boundary interactions. Finally, we summarize all results in a rating table and show how our experiments can be used to determine the promising method for specific requirements regarding a given constellation of fluid‐boundary interaction.

The reasons for the heterogeneity of the reservoir fluid inflow to the cross-sectional contour of a vertical well can be: a) asymmetry of the external pressure field relative to the well axis as a result of the interference of the surrounding wells; b) heterogeneity of the permeability field near the well, which is a consequence of either the heterogeneity of the absolute permeability field of the reservoir or the mobility function of the multiphase mixture of formation fluids. To simulate filtration in a reservoir over a relatively long time interval, the main interest is constant or long-term factors associated with well spacing and the distribution of absolute permeability. In the work, solutions of two model problems were constructed, which allow a quantitative evaluation of the influence of both factors on the degree of inhomogeneity of the inflow to the well and indicate the conditions under which this effect becomes significant. The obtained estimates are intended primarily for computational schemes of streamline and streamtube methods, which require a high degree of solution detailing near wells.