Steady filtration problems with seawater intrusion: macro‐hybrid penalized finite element approximations

Macro‐hybrid penalized finite element approximations are studied for steady filtration problems with seawater intrusion. On the basis of nonoverlapping domain decompositions with vertical interfaces, sections of coastal aquifers are decomposed into subsystems with simpler geometries and small scales, interconnected via transmission conditions of pressure and flux continuity. Corresponding local penalized formulations are derived from the global penalized variational formulation of the two‐free boundary flow problem, with continuity transmission conditions modelled variationally in a dual sense. Then, macro‐hybrid finite element approximations are derived for the system, defined on independent subdomain grids. Parallel relaxation penalty‐duality algorithms are proposed from fixed‐point problem characterizations. Numerical experiments exemplify the macro‐hybrid penalized theory, showing a good agreement with previous primal conforming penalized finite element approximations (Comput. Methods Appl. Mech. Engng. 2000; 190 :609–624). Copyright © 2005 John Wiley & Sons, Ltd.     

Pressure‐stabilized maximum‐entropy methods for incompressible Stokes

We present a parameter‐free stable maximum‐entropy method for incompressible Stokes flow. Derived from a least‐biased optimization inspired by information theory, the meshfree maximum‐entropy method appears as an interesting alternative to classical approximation schemes like the finite element method. Especially compared with other meshfree methods, e.g. the moving least‐squares method, it allows for a straightforward imposition of boundary conditions. However, no Eulerian approach has yet been presented for real incompressible flow, encountering the convective and pressure instabilities. In this paper, we exclusively address the pressure instabilities caused by the mixed velocity‐pressure formulation of incompressible Stokes flow. In a preparatory discussion, existing stable and stabilized methods are investigated and evaluated. This is used to develop different approaches towards a stable maximum‐entropy formulation. We show results for two analytical tests, including a presentation of the convergence behavior. As a typical benchmark problem, results are also shown for the leaky lid‐driven cavity. The already presented information‐flux method for convection‐dominated problems in mind, we see this as the last step towards a maximum‐entropy method capable of simulating full incompressible flow problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Time‐domain BEM solution of convection–diffusion‐type MHD equations

The two‐dimensional convection–diffusion‐type equations are solved by using the boundary element method (BEM) based on the time‐dependent fundamental solution. The emphasis is given on the solution of magnetohydrodynamic (MHD) duct flow problems with arbitrary wall conductivity. The boundary and time integrals in the BEM formulation are computed numerically assuming constant variations of the unknowns on both the boundary elements and the time intervals. Then, the solution is advanced to the steady‐state iteratively. Thus, it is possible to use quite large time increments and stability problems are not encountered. The time‐domain BEM solution procedure is tested on some convection–diffusion problems and the MHD duct flow problem with insulated walls to establish the validity of the approach. The numerical results for these sample problems compare very well to analytical results. Then, the BEM formulation of the MHD duct flow problem with arbitrary wall conductivity is obtained for the first time in such a way that the equations are solved together with the coupled boundary conditions. The use of time‐dependent fundamental solution enables us to obtain numerical solutions for this problem for the Hartmann number values up to 300 and for several values of conductivity parameter. Copyright © 2007 John Wiley & Sons, Ltd.     

Geometrical interpretation of the multi‐point flux approximation L‐method

In this paper, we first investigate the influence of different Dirichlet boundary discretizations on the convergence rate of the multi‐point flux approximation (MPFA) L‐method by the numerical comparisons between the MPFA O‐ and L‐method, and show how important it is for this new method to handle Dirichlet boundary conditions in a suitable way. A new Dirichlet boundary strategy is proposed, which in some sense can well recover the superconvergence rate of the normal velocity. In the second part of the work, the MPFA L‐method with homogeneous media is studied. A systematic concept and geometrical interpretations of the L‐method are given and illustrated, which yield more insight into the L‐method. Finally, we apply the MPFA L‐method for two‐phase flow in porous media on different quadrilateral grids and compare its numerical results for the pressure and saturation with the results of the two‐point flux approximation method. Copyright © 2008 John Wiley & Sons, Ltd.     

A stress‐based least‐squares finite‐element model for incompressible Navier–Stokes equations

In this paper we present a stress‐based least‐squares finite‐element formulation for the solution of the Navier–Stokes equations governing flows of viscous incompressible fluids. Stress components are introduced as independent variables to make the system first order. Continuity equation becomes an algebraic equation and is eliminated from the system with suitable modifications. The h and p convergence are verified using the exact solution of Kovasznay flow. Steady flow past a large circular cylinder in a channel is solved to test mass conservation. Transient flow over a backward‐facing step problem is solved on several meshes. Results are compared with that obtained using vorticity‐based first‐order formulation for both benchmark problems. Copyright © 2007 John Wiley & Sons, Ltd.     

The consistency of pressure‐gradient approximations used in multi‐dimensional shock hydrodynamics

This paper presents a study of the consistency properties of the pressure‐gradient approximation used in multi‐dimensional finite‐element shock hydrodynamics codes today. In specific, consideration is given to the so‐called ‘bent‐element blues’ problem associated with the pressure‐gradient approximation when using the Q1Q0 element. On arbitrary grids comprised of distorted elements, the piecewise‐constant representation of the pressure field leads to a low‐order pressure‐gradient approximation at the global (nodal) level. This results in spurious nodal forces that are not aligned with the pressure gradient. There are several side‐effects of this behavior that include (a) incorrectly exciting physical modes in problems that exhibit unstable behavior, e.g. Rayleigh–Taylor problems (both magnetic and hydrodynamic), (b) potentially seeding hourglass modes, and (c) exhibiting non‐stationary behavior for steady‐state problems. A series of commonly used pressure‐gradient approximations are reviewed and evaluated based on linear consistency—the ability of the approximation to annihilate constant terms and exactly reproduce a linear gradient. The deeper theoretical issues associated with the proper selection of function spaces for the finite‐element hydro formulation are not discussed here. There are two gradient approximations that use piecewise‐constant data and deliver a consistent pressure‐gradient approximation on arbitrary grids. The first is the well‐known least‐squares gradient construction, and the second is a corrected gradient approximation that imposes linear consistency at the (global) nodal level. At the time of this writing, the corrected gradient approximation appears to be the most viable candidate for resolving the consistency issues associated with the Q1Q0 element technology. Copyright © 2009 John Wiley & Sons, Ltd.     

FEM‐BEM coupling for the exterior Stokes problem with non‐conforming finite elements and an application to small droplet deformation dynamics

A non‐conforming, discontinuous Galerkin finite element–boundary element coupling procedure is presented for the exterior planar Stokes problem. The novel coupled formulation is developed using that for the conforming case as a guide to the introduction of extra mortar variables used to couple a discontinuous interior finite element solution with a continuous exterior boundary element solution. Convergence results for the new scheme are presented, for a range of different interior penalties, on computational domains discretized with regular structured meshes. To illustrate an application, the excitations required to model two‐phase droplet deformations in an extensional flow, under simple surface tension, with the new scheme are also presented. For a selection of different drop viscocities and exterior flows, with and without a rotational component, the progression to a steady‐state deformation of initially undeformed circular drops is calculated and the results compared with those from both a conforming FEM‐BEM equivalent scheme and from a small perturbation analysis where available. Copyright © 2011 John Wiley & Sons, Ltd.     

Least‐squares meshfree method for incompressible Navier–Stokes problems

A least‐squares meshfree method based on the first‐order velocity–pressure–vorticity formulation for two‐dimensional incompressible Navier–Stokes problem is presented. The convective term is linearized by successive substitution or Newton's method. The discretization of all governing equations is implemented by the least‐squares method. Equal‐order moving least‐squares approximation is employed with Gauss quadrature in the background cells. The boundary conditions are enforced by the penalty method. The matrix‐free element‐by‐element Jacobi preconditioned conjugate method is applied to solve the discretized linear systems. Cavity flow for steady Navier–Stokes problem and the flow over a square obstacle for time‐dependent Navier–Stokes problem are investigated for the presented least‐squares meshfree method. The effects of inaccurate integration on the accuracy of the solution are investigated. Copyright © 2004 John Wiley & Sons, Ltd.     

Positive‐definite q‐families of continuous subcell Darcy‐flux CVD(MPFA) finite‐volume schemes and the mixed finite element method

A new family of locally conservative cell‐centred flux‐continuous schemes is presented for solving the porous media general‐tensor pressure equation. A general geometry‐permeability tensor approximation is introduced that is piecewise constant over the subcells of the control volumes and ensures that the local discrete general tensor is elliptic. A family of control‐volume distributed subcell flux‐continuous schemes are defined in terms of the quadrature parametrization q (Multigrid Methods. Birkhauser: Basel, 1993; Proceedings of the 4th European Conference on the Mathematics of Oil Recovery, Norway, June 1994; Comput. Geosci. 1998; 2 :259–290), where the local position of flux continuity defines the quadrature point and each particular scheme. The subcell tensor approximation ensures that a symmetric positive‐definite (SPD) discretization matrix is obtained for the base member (q=1) of the formulation. The physical‐space schemes are shown to be non‐symmetric for general quadrilateral cells. Conditions for discrete ellipticity of the non‐symmetric schemes are derived with respect to the local symmetric part of the tensor. The relationship with the mixed finite element method is given for both the physical‐space and subcell‐space q‐families of schemes. M‐matrix monotonicity conditions for these schemes are summarized. A numerical convergence study of the schemes shows that while the physical‐space schemes are the most accurate, the subcell tensor approximation reduces solution errors when compared with earlier cell‐wise constant tensor schemes and that subcell tensor approximation using the control‐volume face geometry yields the best SPD scheme results. A particular quadrature point is found to improve numerical convergence of the subcell schemes for the cases tested. Copyright © 2007 John Wiley & Sons, Ltd.     

A spectral/hp least‐squares finite element analysis of the Carreau–Yasuda fluids

A least‐squares finite element model with spectral/hp approximations was developed for steady, two‐dimensional flows of non‐Newtonian fluids obeying the Carreau–Yasuda constitutive model. The finite element model consists of velocity, pressure, and stress fields as independent variables (hence, called a mixed model). Least‐squares models offer an alternative variational setting to the conventional weak‐form Galerkin models for the Navier–Stokes equations, and no compatibility conditions on the approximation spaces used for the velocity, pressure, and stress fields are necessary when the polynomial order (p) used is sufficiently high (say, p > 3, as determined numerically). Also, the use of the spectral/hp elements in conjunction with the least‐squares formulation with high p alleviates various forms of locking, which often appear in low‐order least‐squares finite element models for incompressible viscous fluids, and accurate results can be obtained with exponential convergence. To verify and validate, benchmark problems of Kovasznay flow, backward‐facing step flow, and lid‐driven square cavity flow are used. Then the effect of different parameters of the Carreau–Yasuda constitutive model on the flow characteristics is studied parametrically. Copyright © 2016 John Wiley & Sons, Ltd.     

A stream function–vorticity formulation‐based immersed boundary method and its applications

A new stream function–vorticity formulation‐based immersed boundary method is presented in this paper. Different from the conventional immersed boundary method, the main feature of the present model is to accurately satisfy both governing equations and boundary conditions through velocity correction and vorticity correction procedures. The velocity correction process is performed implicitly based on the requirement that velocity at the immersed boundary interpolated from the corrected velocity field accurately satisfies the nonslip boundary condition. The vorticity correction is made through the stream function formulation rather than the vorticity transport equation. It is evaluated from the firstorder derivatives of velocity correction. Two simple and efficient ways are presented for approximation of velocity‐correction derivatives. One is based on finite difference approximation, while the other is based on derivative expressions of Dirac delta function and velocity correction. It was found that both ways can work very well. The main advantage of the proposed method lies in its simple concept, easy implementation, and robustness in stability. Numerical experiments for both stationary and moving boundary problems were conducted to validate the capability and efficiency of the present method. Good agreements with available data in the literature were achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

Local Heaviside‐weighted LRPIM meshless method and its application to two‐dimensional potential flows

In this paper, the local radial point interpolation meshless method (LRPIM) is used for the analysis of two‐dimensional potential flows, based on a local‐weighted residual method with the Heaviside step function as the weighting function over a local subdomain. Trial functions are constructed using radial basis functions. The present method is a truly meshless method based only on a number of randomly located nodes. Integration over the subdomains requires only a simple integration cell to obtain the solution. No element matrix assembly is required and no special treatment is needed to impose the essential boundary conditions. The novelty of the paper is the use of a local Heaviside weight function in the LRPIM, which does not need local domain integration and integrations only on the boundary of the local domains are needed. Effects of the sizes of local subdomain and interpolation domain on the performance of the present method are investigated. The behavior of shape parameters of multiquadrics has been systematically studied. Two numerical tests in groundwater and fluid flows are presented and compared with closed‐form solutions and finite element method. The results show that the use of a local Heaviside weight function in the LRPIM is highly accurate and possesses no numerical difficulties. Copyright © 2008 John Wiley & Sons, Ltd.     

A family of multi‐point flux approximation schemes for general element types in two and three dimensions with convergence performance

A family of flux‐continuous, locally conservative, control‐volume‐distributed multi‐point flux approximation (CVD‐MPFA) schemes has been developed for solving the general geometry‐permeability tensor pressure equation on structured and unstructured grids. These schemes are applicable to the full‐tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir simulation schemes when applied to full‐tensor flow approximation. The family of flux‐continuous schemes is characterized by a quadrature parameterization. Improved numerical convergence for the family of CVD‐MPFA schemes using the quadrature parameterization has been observed for structured and unstructured grids in two dimensions. The CVD‐MPFA family cell‐vertex formulation is extended to classical general element types in 3‐D including prisms, pyramids, hexahedra and tetrahedra. A numerical convergence study of the CVD‐MPFA schemes on general unstructured grids comprising of triangular elements in 2‐D and prismatic, pyramidal, hexahedral and tetrahedral shape elements in 3‐D is presented. Copyright © 2011 John Wiley & Sons, Ltd.     

The effects of control‐volume distributed multi‐point flux approximation (CVD‐MPFA) on upscaling—A study using the CVD‐MPFA schemes

A family of flux‐continuous, locally conservative, finite‐volume schemes has been developed for solving the general geometry‐permeability tensor (petroleum reservoir‐simulation) pressure equation on structured and unstructured grids and are control‐volume distributed (textit Comput. Geo. 1998; 2 :259–290; Comput. Geo. 2002; 6 :433–452). The schemes are applicable to diagonal and full tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir‐simulation schemes (two‐point flux approximation) when applied to full tensor flow approximation. The family of flux‐continuous schemes is quantified by a quadrature parameterization (Int. J. Numer. Meth. Fluids 2006; 51 :1177–1203). Improved convergence (for two‐ and three‐dimensional formulation) using the quadrature parameterization has been observed for the family of flux‐continuous control‐volume distributed multi‐point flux approximation (CVD‐MPFA) schemes (Ph.D. Thesis, University of Wales, Swansea, U.K., 2007). In this paper family of flux‐continuous (CVD‐MPFA) schemes are used as a part of numerical upscaling procedure for upscaling the fine‐scale grid information (permeability) onto a coarse grid scale. A series of data‐sets (SPE, 2001) are tested where the upscaled permeability tensor is computed on a sequence of grid levels using the same fixed range of quadrature points in each case. The refinement studies presented involve:

• (i) Refinement comparison study: In this study, permeability distribution for cells at each grid level is obtained by upscaling directly from the fine‐scale permeability field as in standard simulation practice.
• (ii) Refinement study with renormalized permeability: In this refinement comparison, the local permeability is upscaled to the next grid level hierarchically, so that permeability values are renormalized to each coarser level. Hence, showing only the effect of increased grid resolution on upscaled permeability, compared with that obtained directly from the fine‐scale solution.
• (iii) Refinement study with invariant permeability distribution: In this study, a classical mathematical convergence test is performed. The same coarse‐scale underlying permeability map is preserved on all grid levels including the fine‐scale reference solution.

The study is carried out for the discretization of the scheme in physical space. The benefit of using specific quadrature points is demonstrated for upscaling in this study and superconvergence is observed. Copyright © 2010 John Wiley & Sons, Ltd.     

A finite element thermally coupled flow formulation for phase‐change problems

A finite element, thermally coupled incompressible flow formulation considering phase‐change effects is presented. This formulation accounts for natural convection, temperature‐dependent material properties and isothermal and non‐isothermal phase‐change models. In this context, the full Navier–Stokes equations are solved using a generalized streamline operator (GSO) technique. The highly non‐linear phase‐change effects are treated with a temperature‐based algorithm, which provides stability and convergence of the numerical solution. The Boussinesq approximation is used in order to consider the temperature‐dependent density variation. Furthermore, the numerical solution of the coupled problem is approached with a staggered incremental‐iterative solution scheme, such that the convergence criteria are written in terms of the residual vectors. Finally, this formulation is used for the solutions of solidification and melting problems validating some numerical results with other existing solutions obtained with different methodologies. Copyright © 2000 John Wiley & Sons, Ltd.     

Smoothed‐profile method for momentum and heat transfer in particulate flows

The smoothed‐profile method for the motion of solid bodies suspended in a fluid phase is investigated when combined with a high‐order spatial discretization. The performance of the combined method is tested for a wide range of flow and geometry parameters as well as for static and for moving particles. Moreover, a sensitivity analysis is conducted with respect to the smoothed‐profile function. The algorithm is extended to include thermal effects in Boussinesq approximation. Several benchmark problems are considered to demonstrate the potential of the technique. The implementation of the energy equation is verified by dedicated tests. All simulations are compared with either theoretical, numerical, or experimental data. The results demonstrate the accuracy and efficiency of the smoothed‐profile method for non‐isothermal problems in combination with a discontinuous finite‐element solver for the fluid flow, which allows for a flexible handling of the grid and the order of spectral approximation in each element. Copyright © 2016 John Wiley & Sons, Ltd.     

Control of hydrothermal waves in a thermocapillary flow using a gradient‐based control strategy

The problem of controlling the hydrothermal waves in a thermocapillary flow is addressed using a gradient‐based control strategy. The state equations are the two‐dimensional unsteady incompressible Navier–Stokes and energy equations under the Boussinesq approximation. The modeled problem is the ‘open boat’ process of crystal growth, the flow which is driven by Marangoni and buoyancy effects. The control is a spatially and temporally varying heat flux boundary condition at the free surface. The control that minimizes the hydrothermal waves is found using a conjugate gradient method, where the gradient of the objective function with respect to the control variables is obtained from solving a set of adjoint equations. The effectiveness of choices of the parameters governing the control algorithm is examined. Almost complete suppression of the hydrothermal waves is obtained for certain choices of the parameters governing the control algorithm. The numerical issues involved with finding the control using the optimizer are discussed, and the features of the resulting control are analyzed with the goal of understanding how it affects the flow.Copyright © 2012 John Wiley & Sons, Ltd.     

Local moving least square‐one‐dimensional integrated radial basis function networks technique for incompressible viscous flows

This paper presents a local moving least square‐one‐dimensional integrated radial basis function networks method for solving incompressible viscous flow problems using stream function‐vorticity formulation. In this method, the partition of unity method is employed as a framework to incorporate the moving least square and one‐dimensional integrated radial basis function networks techniques. The major advantages of the proposed method include the following: (i) a banded sparse system matrix which helps reduce the computational cost; (ii) the Kronecker‐ δ property of the constructed shape function which helps impose the essential boundary condition in an exact manner; and (iii) high accuracy and fast convergence rate owing to the use of integration instead of conventional differentiation to construct the local radial basis function approximations. Several examples including two‐dimensional (2D) Poisson problems, lid‐driven cavity flow and flow past a circular cylinder are considered, and the present results are compared with the exact solutions and numerical results from other methods in the literature to demonstrate the attractiveness of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Three‐dimensional transient free‐surface flow of viscous fluids inside cavities of arbitrary shape

The three‐dimensional transient free‐surface flow inside cavities of arbitrary shape is examined in this study. An adaptive (Lagrangian) boundary‐element approach is proposed for the general three‐dimensional simulation of confined free‐surface flow of viscous incompressible fluids. The method is stable as it includes remeshing capabilities of the deforming free‐surface, and thus can handle large deformations. A simple algorithm is developed for mesh refinement of the deforming free‐surface mesh. Smooth transition between large and small elements is achieved without significant degradation of the aspect ratio of the elements in the mesh. The method is used to determine the flow field and free‐surface evolution inside cubic, rectangular and cylindrical containers. These problems illustrate the transient nature of the flow during the mixing process. Surface tension effects are also explored. Copyright © 2003 John Wiley & Sons, Ltd.