A new method for Stokes problems with circular boundaries using degenerate kernel and Fourier series

This study is concerned with the Stokes flow of an incompressible fluid of constant density and viscosity with circular boundaries. To fully capture the circular boundary, the boundary densities in the direct and indirect boundary integral equations (BIEs) are expanded in terms of Fourier series. The kernel functions in either the direct BIE or the indirect BIE are expanded to degenerate kernels by using the separation of field and source points. Consequently, the improper integrals are transformed to series sum and are easily calculated. The linear algebraic system can be established by matching the boundary conditions at the collocation points. Then, the unknown Fourier coefficients can be easily determined. Finally, several examples including circular and eccentric domains are presented to demonstrate the validity of the present method. Five gains were obtained: (1) meshless approach; (2) free of boundary‐layer effect; (3) singularity free; (4) exponential convergence; and (5) well‐posed model. Copyright © 2007 John Wiley & Sons, Ltd.     

Analysis of three‐dimensional fracture mechanics problems: A two‐scale approach using coarse‐generalized FEM meshes

This paper presents a generalized finite element method (GFEM) based on the solution of interdependent global (structural) and local (crack)‐scale problems. The local problems focus on the resolution of fine‐scale features of the solution in the vicinity of three‐dimensional cracks, while the global problem addresses the macro‐scale structural behavior. The local solutions are embedded into the solution space for the global problem using the partition of unity method. The local problems are accurately solved using an hp‐GFEM and thus the proposed method does not rely on analytical solutions. The proposed methodology enables accurate modeling of three‐dimensional cracks on meshes with elements that are orders of magnitude larger than the process zone along crack fronts. The boundary conditions for the local problems are provided by the coarse global mesh solution and can be of Dirichlet, Neumann or Cauchy type. The effect of the type of local boundary conditions on the performance of the proposed GFEM is analyzed. Several three‐dimensional fracture mechanics problems aimed at investigating the accuracy of the method and its computational performance, both in terms of problem size and CPU time, are presented. Copyright © 2009 John Wiley & Sons, Ltd.     

Shifted FSAI preconditioners for the efficient parallel solution of non‐linear groundwater flow models

The present paper investigates the performance of a shifted factorized sparse approximate inverse as a parallel preconditioner for the iterative solution to the linear systems arising in the finite element discretization of non‐linear groundwater flow models. The shift strategy is based on an inexpensive preconditioner update exploiting the structure of the coefficient matrix. The proposed algorithm is experimented with in the parallel simulation of a large‐scale real multi‐aquifer system characterized by a stochastic distribution of the hydraulic conductivity. The numerical results show that the shifted factorized sparse approximate inverse algorithm may yield an overall computational gain up to 300% with respect to the non‐shifted scheme with an excellent parallel efficiency. Copyright © 2011 John Wiley & Sons, Ltd.     

Stokes flow around cylinders in a bounded two‐dimensional domain using multipole‐accelerated boundary element methods

The multipole technique has recently received attention in the field of boundary element analysis as a means of reducing the order of data storage and calculation time requirements from O(N²) (iterative solvers) or O(N³) (gaussian elimination) to O(N log N) or O(N), where N is the number of nodes in the discretized system. Such a reduction in the growth of the calculation time and data storage is crucial in applications where N is large, such as when modelling the macroscopic behaviour of suspensions of particles. In such cases, a minimum of 1000 particles is needed to obtain statistically meaningful results, leading to systems with N of the order of 10 000 for the smallest problems. When only boundary velocities are known, the indirect boundary element formulation for Stokes flow results in Fredholm equations of the second kind, which generally produce a well‐posed set of equations when discretized, a necessary requirement for iterative solution methods. The direct boundary element formulation, on the other hand, results in Fredholm equations of the first kind, which, upon discretization, produce ill‐conditioned systems of equations. The model system here is a two‐dimensional wide‐gap couette viscometer, where particles are suspended in the fluid between the cylinders. This is a typical system that is efficiently modelled using boundary element method simulations. The multipolar technique is applied to both direct and indirect formulations. It is found that the indirect approach is sufficiently well‐conditioned to allow the use of fast multipole methods. The direct approach results in severe ill‐conditioning, to a point where application of the multipole method leads to non‐convergence of the solution iteration. Copyright © 1999 John Wiley & Sons, Ltd.     

Bubble‐enhanced smoothed finite element formulation: a variational multi‐scale approach for volume‐constrained problems in two‐dimensional linear elasticity

This paper presents a bubble‐enhanced smoothed finite element formulation for the analysis of volume‐constrained problems in two‐dimensional linear elasticity. The new formulation is derived based on the variational multi‐scale approach in which unequal order displacement‐pressure pairs are used for the mixed finite element approximation and hierarchical bubble function is selected for the fine‐scale displacement approximation. An area‐weighted averaging scheme is employed for the two‐scale smoothed strain calculation under the framework of edge‐based smoothed FEM. The smoothed fine‐scale solution is shown to naturally contain the stress field jump of the smoothed coarse‐scale solution across the boundary of edge‐based smoothing domain and thus provides the possibility to stabilize the global solution for volume‐constrained problems. A global monolithic solution strategy is employed, and the fine‐scale solution is solved without the consideration of approximating the strong form of the fine‐scale equation. Several numerical examples are analyzed to demonstrate the accuracy of the present formulation. Copyright © 2014 John Wiley & Sons, Ltd.     

Insight into the flow‐condition‐based interpolation finite element approach: solution of steady‐state advection–diffusion problems

The flow‐condition‐based interpolation (FCBI) finite element approach is studied in the solution of advection–diffusion problems. Two FCBI procedures are developed and tested with the original FCBI method: in the first scheme, a general solution of the advection–diffusion equation is embedded into the interpolation, and in the second scheme, the link‐cutting bubbles approach is used in the interpolation. In both procedures, as in the original FCBI method, no artificial parameters are included to reach stability for high Péclet number flows. The procedures have been implemented for two‐dimensional analysis and the results of some test problems are presented. These results indicate good stability and accuracy characteristics and the potential of the FCBI solution approach. Copyright © 2005 John Wiley & Sons, Ltd.     

Multiscale Galerkin method using interpolation wavelets for two‐dimensional elliptic problems in general domains

One major hurdle in developing an efficient wavelet‐based numerical method is the difficulty in the treatment of general boundaries bounding two‐ or three‐dimensional domains. The objective of this investigation is to develop an adaptive multiscale wavelet‐based numerical method which can handle general boundary conditions along curved boundaries. The multiscale analysis is achieved in a multi‐resolution setting by employing hat interpolation wavelets in the frame of a fictitious domain method. No penalty term or the Lagrange multiplier need to be used in the present formulation. The validity of the proposed method and the effectiveness of the multiscale adaptive scheme are demonstrated by numerical examples dealing with the Dirichlet and Neumann boundary‐value problems in quadrilateral and quarter circular domains. Copyright © 2003 John Wiley & Sons, Ltd.     

A two‐dimensional flow model for the process simulation of complex shape composite laminates

A numerical flow‐compaction model is developed and implemented in a finite element code to simulate the multiple physical phenomena involved during the autoclave processing of fibre‐reinforced composite laminates. The model is based on the effective stress formulation coupled with a Darcian flow theory. A Galerkin approach is employed to discretize the weak form of the governing equations. The current formulation successfully describes the compaction behaviour of complex shape laminates caused by flow of the resin. A parametric study is performed to investigate the effect of the material properties on the compaction of angle‐shaped composite laminates. It is found that the fibre bed shear modulus significantly affects the compaction behaviour in the corner sections of curved laminates while the resin viscosity and fibre bed permeability affect the compaction rate of the laminate. Copyright © 1999 John Wiley & Sons, Ltd.     

A two‐scale approach for fluid flow in fractured porous media

A two‐scale numerical model is developed for fluid flow in fractured, deforming porous media. At the microscale the flow in the cavity of a fracture is modelled as a viscous fluid. From the micromechanics of the flow in the cavity, coupling equations are derived for the momentum and the mass couplings to the equations for a fluid‐saturated porous medium, which are assumed to hold on the macroscopic scale. The finite element equations are derived for this two‐scale approach and integrated over time. By exploiting the partition‐of‐unity property of the finite element shape functions, the position and direction of the fractures is independent from the underlying discretization. The resulting discrete equations are non‐linear due to the non‐linearity of the coupling terms. A consistent linearization is given for use within a Newton–Raphson iterative procedure. Finally, examples are given to show the versatility and the efficiency of the approach, and show that faults in a deforming porous medium can have a significant effect on the local as well as on the overall flow and deformation patterns. Copyright © 2006 John Wiley & Sons, Ltd.     

A two‐dimensional stochastic algorithm for the solution of the non‐linear Poisson–Boltzmann equation: validation with finite‐difference benchmarks

This paper presents a two‐dimensional floating random walk (FRW) algorithm for the solution of the non‐linear Poisson–Boltzmann (NPB) equation. In the past, the FRW method has not been applied to the solution of the NPB equation which can be attributed to the absence of analytical expressions for volumetric Green's functions. Previous studies using the FRW method have examined only the linearized Poisson–Boltzmann equation. No such linearization is needed for the present approach. Approximate volumetric Green's functions have been derived with the help of perturbation theory, and these expressions have been incorporated within the FRW framework. A unique advantage of this algorithm is that it requires no discretization of either the volume or the surface of the problem domains. Furthermore, each random walk is independent, so that the computational procedure is highly parallelizable. In our previous work, we have presented preliminary calculations for one‐dimensional and quasi‐one‐dimensional benchmark problems. In this paper, we present the detailed formulation of a two‐dimensional algorithm, along with extensive finite‐difference validation on fully two‐dimensional benchmark problems. The solution of the NPB equation has many interesting applications, including the modelling of plasma discharges, semiconductor device modelling and the modelling of biomolecular structures and dynamics. Copyright © 2005 John Wiley & Sons, Ltd.     

Accuracy of Galerkin finite elements for groundwater flow simulations in two and three‐dimensional triangulations

In standard finite element simulations of groundwater flow the correspondence between hydraulic head gradients and groundwater fluxes is represented by the stiffness matrix. In two‐dimensional problems the use of linear triangular elements on Delaunay triangulations guarantees a stiffness matrix of type M. This implies that the local numerical fluxes are physically consistent with Darcy's law. This condition is fundamental to avoid the occurrence of local maxima or minima, and is of crucial importance when the calculated flow field is used in contaminant transport simulations or pathline evaluation. In three spatial dimensions, the linear Galerkin approach on tetrahedra does not lead to M‐matrices even on Delaunay meshes. By interpretation of the Galerkin approach as a subdomain collocation scheme, we develop a new approach (OSC, orthogonal subdomain collocation) that is shown to produce M‐matrices in three‐dimensional Delaunay triangulations. In case of heterogeneous and anisotropic coefficients, extra mesh properties required for M‐stiffness matrices will also be discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

A finite difference scheme for solving a three‐dimensional heat transport equation in a thin film with microscale thickness

Heat transport at the microscale is of vital importance in microtechnology applications. The heat transport equation is different from the traditional heat diffusion equation since a second‐order derivative of temperature with respect to time and a third‐order mixed derivative of temperature with respect to space and time are introduced. In this study, we develop a finite difference scheme with two levels in time for the three‐dimensional heat transport equation. It is shown by the discrete energy method that the scheme is unconditionally stable. The three‐dimensional implicit scheme is then solved by using a preconditioned Richardson iteration, so that only a tridiagonal linear system is solved each iteration. Numerical results show that the solution is accurate. Copyright © 2001 John Wiley & Sons, Ltd.     

Guaranteed energy error bounds for the Poisson equation using a flux‐free approach: Solving the local problems in subdomains

A method to compute guaranteed upper bounds for the energy norm of the exact error in the finite element solution of the Poisson equation is presented. The bounds are guaranteed for any finite element mesh however coarse it may be, not just in the asymptotic regime. The bounds are constructed by employing a subdomain‐based a posteriori error estimate which yields self‐equilibrated residual loads in stars (patches of elements). The proposed approach is an alternative to standard equilibrated residual methods providing sharper bounds. The use of a flux‐free error estimator improves the effectivities of the upper bounds for the energy while retaining the certainty of the bounds. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical method for the solution of non‐linear two‐dimensional inverse heat conduction problem using unstructured meshes

A new method of solving multidimensional heat conduction problems is formulated. The developed space marching method allows to determine quickly and exactly unsteady temperature distributions in the construction elements of irregular geometry. The method which is based on temperature measurements at the outer surface, is especially appropriate for determining transient temperature distribution in thick‐wall pressure components. Two examples are included to demonstrate the capabilities of the new approach. Copyright © 2000 John Wiley & Sons, Ltd.     

An explicit residual‐type error estimator for ‐quadrilateral extended finite element method in two‐dimensional linear elastic fracture mechanics

This paper deals with explicit residual a posteriori error estimation analysis for ‐quadrilateral extended finite element method (XFEM) discretizations applied to the two‐dimensional problem of linear elastic fracture mechanics. The result is twofold. First, to enable estimation procedures with application to XFEM, a specific quasi‐interpolation operator of averaging type is constructed. The main challenge here arises from the different types of enrichments implemented, and hence, to impose the constant‐preserving property of the interpolation operator on an element, we use the idea of an extension operator. An upper bound on the discretization error measured in the energy norm and associated local error indicators are then constructed and analyzed. The second result follows from the error analysis and concerns an alternative choice of branch functions used in XFEM applications. In particular, the branch functions have to be chosen to fulfill the divergence‐free conditions within the crack tip element and traction‐free boundary conditions on the crack faces. Then, the corresponding XFEM solution gains a better accuracy with less degrees of freedom. Finally, numerical examples are provided with comparative results. Copyright © 2012 John Wiley & Sons, Ltd.     

Integration of singular enrichment functions in the generalized/extended finite element method for three‐dimensional problems

A mapping method is developed to integrate weak singularities, which result from enrichment functions in the generalized/extended finite element method. The integration scheme is applicable to 2D and 3D problems including arbitrarily shaped triangles and tetrahedra. Implementation of the proposed scheme in existing codes is straightforward. Numerical examples for 2D and 3D problems demonstrate the accuracy and convergence properties of the technique. Copyright © 2008 John Wiley & Sons, Ltd.     

An algorithm for modelling the interaction of a flexible rod with a two‐dimensional high‐speed flow

We present an algorithm for modelling coupled dynamic interactions of a very thin flexible structure immersed in a high‐speed flow. The modelling approach is based on combining an Eulerian finite volume formulation for the fluid flow and a Lagrangian large‐deformation formulation for the dynamic response of the structure. The coupling between the fluid and the solid response is achieved via an approach based on extrapolation and velocity reconstruction inspired in the Ghost Fluid Method. The algorithm presented does not assume the existence of a region exterior to the fluid domain as it was previously proposed and, thus, enables the consideration of very thin open boundaries and structures where the flow may be relevant on both sides of the interface. We demonstrate the accuracy of the method and its ability to describe disparate flow conditions across a fixed thin rigid interface without pollution of the flow field across the solid interface by comparing with analytical solutions of compressible flows. We also demonstrate the versatility and robustness of the method in a complex fluid–structure interaction problem corresponding to the transient supersonic flow past a highly flexible structure. Copyright © 2005 John Wiley & Sons, Ltd.     

Non‐reflecting artificial boundaries for transient scalar wave propagation in a two‐dimensional infinite homogeneous layer

This paper presents an exact non‐reflecting boundary condition for dealing with transient scalar wave propagation problems in a two‐dimensional infinite homogeneous layer. In order to model the complicated geometry and material properties in the near field, two vertical artificial boundaries are considered in the infinite layer so as to truncate the infinite domain into a finite domain. This treatment requires the appropriate boundary conditions, which are often referred to as the artificial boundary conditions, to be applied on the truncated boundaries. Since the infinite extension direction is different for these two truncated vertical boundaries, namely one extends toward x →∞ and another extends toward x→‐ ∞, the non‐reflecting boundary condition needs to be derived on these two boundaries. Applying the variable separation method to the wave equation results in a reduction in spatial variables by one. The reduced wave equation, which is a time‐dependent partial differential equation with only one spatial variable, can be further changed into a linear first‐order ordinary differential equation by using both the operator splitting method and the modal radiation function concept simultaneously. As a result, the non‐reflecting artificial boundary condition can be obtained by solving the ordinary differential equation whose stability is ensured. Some numerical examples have demonstrated that the non‐reflecting boundary condition is of high accuracy in dealing with scalar wave propagation problems in infinite and semi‐infinite media. Copyright © 2003 John Wiley & Sons, Ltd.