A fully-coupled discontinuous Galerkin method for two-phase flow in porous media with discontinuous capillary pressure

In this paper we formulate and test numerically a fully-coupled discontinuous Galerkin (DG) method for incompressible two-phase flow with discontinuous capillary pressure. The spatial discretization uses the symmetric interior penalty DG formulation with weighted averages and is based on a wetting-phase potential / capillary potential formulation of the two-phase flow system. After discretizing in time with diagonally implicit Runge-Kutta schemes the resulting systems of nonlinear algebraic equations are solved with Newton's method and the arising systems of linear equations are solved efficiently and in parallel with an algebraic multigrid method. The new scheme is investigated for various test problems from the literature and is also compared to a cell-centered finite volume scheme in terms of accuracy and time to solution. We find that the method is accurate, robust and efficient. In particular no post-processing of the DG velocity field is necessary in contrast to results reported by several authors for decoupled schemes. Moreover, the solver scales well in parallel and three-dimensional problems with up to nearly 100 million degrees of freedom per time step have been computed on 1000 processors

Towards large-scale modelling of fluid flow in fractured porous media

To date, the complexity of fractured porous media still precludes the direct incorporation of small-scale features into field-scale modelling. These features, however, can be instrumental in shaping and triggering coarsening instabilities and other forms of emergent behaviour which need to be considered on the field-scale. Here we develop numerical simulation methods for this purpose and demonstrate their improved performance in single-and two-phase flow simulations with models of fractured porous media. Material discontinuities in fractured porous media strongly influence single-and multi-phase fluid flow. When continuum methods are used to model transport across such interfaces, they smear out jump discontinuities of concentration or saturation. To overcome this drawback, we “explode” hybrid finite-element node-centred finite-volume models along these introducing complementary finite-volumes along the material interfaces. With this embedded discontinuity discretization we develop a transport scheme that realistically represents the dependent variable discontinuities arising at these interfaces. The main advantage of this new scheme is its ability to honour the flow effects that we know that these discontinuities have in physical experiments. We have also developed a new time-stepping control scheme for the transport equation. It allows the user to specify the volume fraction of the model in which he/she is prepared to relax the CFL condition. This scheme is applied in a study of the impact of fracture pattern development on solute transport. These two-dimensional simulations quantify the effect of the fractures on macro-scale dispersion in geomechanically generated fracture geometries, as opposed to stochastically generated ones. Among other insights, the results indicate that fracture density, fracture spacing, and the fracture-matrix flux ratio control anomalous mass transport in such media. We also find that it is crucial to embed discontinuities into large-scale models of heterogeneous porous media

Thermophysical Phenomena in Metal Additive Manufacturing by Selective Laser Melting: Fundamentals, Modeling, Simulation and Experimentation

Among the many additive manufacturing (AM) processes for metallic materials, selective laser melting (SLM) is arguably the most versatile in terms of its potential to realize complex geometries along with tailored microstructure. However, the complexity of the SLM process, and the need for predictive relation of powder and process parameters to the part properties, demands further development of computational and experimental methods. This review addresses the fundamental physical phenomena of SLM, with a special emphasis on the associated thermal behavior. Simulation and experimental methods are discussed according to three primary categories. First, macroscopic approaches aim to answer questions at the component level and consider for example the determination of residual stresses or dimensional distortion effects prevalent in SLM. Second, mesoscopic approaches focus on the detection of defects such as excessive surface roughness, residual porosity or inclusions that occur at the mesoscopic length scale of individual powder particles. Third, microscopic approaches investigate the metallurgical microstructure evolution resulting from the high temperature gradients and extreme heating and cooling rates induced by the SLM process. Consideration of physical phenomena on all of these three length scales is mandatory to establish the understanding needed to realize high part quality in many applications, and to fully exploit the potential of SLM and related metal AM processes

Modelling of Multiphase Fluid flow in Heterogeneous Reservoirs

Computational modeling of multiphase fluid flow in highly heterogenous problems with complicated geometries is a challenging problem for reservoir engineers, with a rich research in establishing best methods and approaches. The novelty in this work is centered around the implementation and comparison of simulation results from two software - the open source ICFESRT and the commercial software ECLIPSE - for a two-phase multiphase problem (oilwater) in both simple and complex geometries. The work involves: (a) implementation and comparison of simulation results from the two software on three different, hypothetical but typical geometries; (b) consideration of a real field case and the associated data analysis, rock characterization, and geostatistics of a real field representative of a highly heterogeneous reservoir; and (c) implementation of both software on the real field case for predictions of oil production at the site, and comparison of the simulation results from the two software. The initial comparison of simulation results for was carried out using three hypothetical (but common) geometries, these being: (a) a quarter five spot with one geological layer; (b) the same geometry as in (a) but with a vertical heterogeneity i.e. 5 different geological layers; (c) and lastly a full 5 spot with 5 different geological layers was implemented. Three different mesh resolutions were applied in both software and comparisons were carried out for mesh-independency. The results showed that in all these three scenarios, good agreement was observed between IC-FERST (coarse mesh) and ECLIPSE (fine mesh) with an average percentage difference at the production well ranging between 2.5% and 10.5% for the oil production and 12% and 26% for the water production. Both the ICFERST and ECLIPSE were subsequently implemented on a real, heterogeneous field – which consisted of 25 producing wells and 8 injections wells. Prior to the software implementation, a data analysis and rock characterization was carried out –Using data from the 33 wells. The logging and core data (a total of 30,000 log readings and 1150 core samples) were utilized and a novel rock characterization technique -Balaha Rock Characterization Code- was implemented to allow for the optimal clustering of rock types within the reservoir, The rock characterization resulted in identifying 7 rock types with their unique porosity-hydraulic permeability relationships. Subsequently, geostatistical methods were implemented – which enabled populating the computational cells of the two software with the corresponding reservoir properties (porosity, hydraulic permeability). To achieve the property population into the unstructured computational domain of the ICFERST software, a newly-developed script was written in Matlab and Python. The rock properties data populated on IC-FERST consist of porosity, permeability, relative permeability, capillary pressure and connate water saturation. A further comparison between the IC-FERST simulation results with the corresponding ECLIPSE simulations was carried out – were all simulations were carried out for a period of 40 years. The percentage differences between the two software simulations were estimated for : (i) ten individual production wells and (ii) the total of all production wells. The results showed that a good agreement exists between the IC-FERST and ECLIPSE simulations, with an average percentage difference for the total oil production of 10.5%, the total water production of 26% and the total water injection of 14%. The results for the ten individual wells showed an average percentage difference of 15.5% ranging from 3 to 29% for the oil production in the late time period. Slightly higher differences were observed when the overall period was considered, due to the large difference at the early time period of the simulation. The results indicated that IC-FERST, when incorporating the necessary rock characterization information – which highlight the heterogeneity of the reservoir – can produce results that can compete with the industry standard ECLIPSE. Additional aspects need to be considered within the current real field IC-FERST simulation, the inclusion of possible fractures and faults, as these were incorporated in the computational domain of ECLIPSE. Additional capabilities also still need to be embedded into IC-FERST, such as the incorporation of the fluid density and viscosity variations with pressure and the consideration of the volume factors, in order to enhance its competitiveness with existing commercial reservoirs simulators such as ECLIPSE

The use of the mesh free methods (radial basis functions) in the modeling of radionuclide migration and moving boundary value problems

Recently, the mesh free methods (radial basis functions-RBFs) have emerged as a novel computing method in the scientific and engineering computing community. The numerical solution of partial differential equations (PDEs) has been usually obtained by finite difference methods (FDM), finite element methods (FEM) and boundary elements methods (BEM). These conventional numerical methods still have some drawbacks. For example, the construction of the mesh in two or more dimensions is a nontrivial problem. Solving PDEs using radial basis function (RBF) collocations is an attractive alternative to these traditional methods because no tedious mesh generation is required. We compare the mesh free method, which uses radial basis functions, with the traditional finite difference scheme and analytical solutions. We will present some examples of using RBFs in geostatistical analysis of radionuclide migration modeling. The advection-dispersion equation will be used in the Eulerian and Lagrangian forms. Stefan's or moving boundary value problems will also be presented. The position of the moving boundary will be simulated by the moving data centers method and level set method

Solution strategies for nonlinear conservation laws

Nonlinear conservation laws form the basis for models for a wide range of physical phenomena. Finding an optimal strategy for solving these problems can be challenging, and a good strategy for one problem may fail spectacularly for others. As different problems have different challenging features, exploiting knowledge about the problem structure is a key factor in achieving an efficient solution strategy. Most strategies found in literature for solving nonlinear problems involve a linearization step, usually using Newton's method, which replaces the original nonlinear problem by an iteration process consisting of a series of linear problems. A large effort is then spent on finding a good strategy for solving these linear problems. This involves choosing suitable preconditioners and linear solvers. This approach is in many cases a good choice and a multitude of different methods have been developed. However, the linearization step to some degree involves a loss of information about the original problem. This is not necessarily critical, but in many cases the structure of the nonlinear problem can be exploited to a larger extent than what is possible when working solely on the linearized problem. This may involve knowledge about dominating physical processes and specifically on whether a process is near equilibrium. By using nonlinear preconditioning techniques developed in recent years, certain attractive features such as automatic localization of computations to parts of the problem domain with the highest degree of nonlinearities arise. In the present work, these methods are further refined to obtain a framework for nonlinear preconditioning that also takes into account equilibrium information. This framework is developed mainly in the context of porous media, but in a general manner, allowing for application to a wide range of problems. A scalability study shows that the method is scalable for challenging two-phase flow problems. It is also demonstrated for nonlinear elasticity problems. Some models arising from nonlinear conservation laws are best solved using completely different strategies than the approach outlined above. One such example can be found in the field of surface gravity waves. For special types of nonlinear waves, such as solitary waves and undular bores, the well-known Korteweg-de Vries (KdV) equation has been shown to be a suitable model. This equation has many interesting properties not typical of nonlinear equations which may be exploited in the solver, and strategies usually reserved to linear problems may be applied. In this work includes a comparative study of two discretization methods with highly different properties for this equation