Efficient quadrature for a boundary element method to compute three-dimensional Stokes flow

A collocation-type boundary element method based on bilinear B-splines is used for the numerical solution of the Stokes Dirichlet problem in bounded domains D ? R³. The computation of the influence matrix requires the numerical evaluation of weakly singular integrals on the domain boundary if the usual double-layer potential ansatz is chosen. Here mostly standard methods with disjoint grids for collocation and integration are used. We develop a special integration scheme based on triangular co-ordinates near the singularity and show its efficiency compared with the method mentioned above.     

Two-phase nonlocal integral models with a bi-Helmholtz averaging kernel for nanorods

In this work, the static tensile and free vibration of nanorods are studied via both the strain-driven (StrainD) and stress-driven (StressD) two-phase nonlocal models with a bi-Helmholtz averaging kernel. Merely adjusting the limits of integration, the integral constitutive equation of the Fredholm type is converted to that of the Volterra type and then solved directly via the Laplace transform technique. The unknown constants can be uniquely determined through the standard boundary conditions and two constrained conditions accompanying the Laplace transform process. In the numerical examples, the bi-Helmholtz kernel-based StrainD (or StressD) two-phase model shows consistently softening (or stiffening) effects on both the tension and the free vibration of nanorods with different boundary edges. The effects of the two nonlocal parameters of the bi-Helmholtz kernel-based two-phase nonlocal models are studied and compared with those of the Helmholtz kernel-based models.

     

A short note on the general boundary element method for viscous flows with high Reynolds number

In this paper, the general boundary element method and the parallel computation are employed to solve laminar viscous flows in a driven square cavity, governed by the exact Navier–Stokes equations. Using the solution at Re=0 as the initial approximation, the convergent numerical results for high Reynolds number at Re=7500 are obtained, for the first time, by the boundary element method. This verifies the validity and great potential of the general boundary element method for highly non‐linear problems, which may greatly enlarge application regions of the boundary element method in science and engineering. Copyright © 2003 John Wiley & Sons, Ltd.     

Stream function finite element solution of Stokes flow with penalties

A finite element method for solution of the stream function formulation of Stokes flow is developed. The method involves complete cubic non-conforming (C⁰) triangular Hermite elements. This element fails the patch test. To correct the element and produce a convergent method we employ a penalty method to weakly enforce the desired continuity constraint on the normal derivative across the inter-element boundaries. Successful use of the method is demonstrated to require reduced integration of the inter-element penalty with a 1-point Gauss rule. Error estimates relate the optimal choice of penalty parameter to mesh size and are corroborated by numerical convergence studies. The need for reduced integration is interpreted using rank relations for an associated hybrid method.     

Local discontinuous Galerkin method for a nonlocal viscous conservation laws

The purpose of this article is to investigate high‐order numerical approximations of scalar conservation laws with nonlocal viscous term. The viscous term is given in the form of convolution in space variable. With the help of the characteristic of viscous term, we design a semidiscrete local discontinuous Galerkin (LDG) method to solve the nonlocal model. We prove stability and convergence of semidiscrete LDG method in L² norm. The theoretical analysis reveals that the present numerical scheme is stable with optimal convergence order for the linear case, and it is stable with sub‐optimal convergence order for nonlinear case. To demonstrate the validity and accuracy of our scheme, we test the Burgers equation with two typical nonlocal fractional viscous terms. The numerical results show the convergence order accuracy in space for both linear and nonlinear cases. Some numerical simulations are provided to show the robustness and effectiveness of the present numerical scheme.     

Two notes on nonlocal field theory

In this paper, two fundamental problems completely unsolved in nonlocal field theory are studied. The first is the dependence of nonlocal residuals. By studying this problem, a theorem concerning the relationship between the residuals of nonlocal body force and nonlocal moment of momentum is given and proven. The other problem is how to give the stress boundary conditions in the linear theory of nonlocal elasticity. The stress boundary conditions obtained in this paper can not only answer why the nonlocal stress solution satisfying the boundary conditionst _ji ^(s) n _j ¦O ₂ =p _i (p _i is a constant) on the surface of crack does not exist but also give a model of the molecular cohesive stress on the crack tip.     

A Line-force loading on the surface of a nonlocal elastic half-infinite medium

Transmission of a concentrated force into a half-infinite nonlocal elastic medium is examined. A nonlocal concentrated force boundary condition which is different from that of Nowinski is derived from the nonlocal principle. The nonlocal stress field is obtained analytically. The result eliminates the stress singularity which occurs in the classical stresses and also the nonlocal stresses derived by Nowinski. A theoretical cut-off strength P_c is discussed.     

Structural symmetry and boundary conditions for nonlocal symmetrical problems

The paper deals with the determination of the symmetric model and proper boundary conditions for solving nonlocal elastic symmetric structures. The above concepts, in the context of nonlocal integral elasticity, turn out to be different with respect to the standard ones, classically applied when dealing with local elastic symmetric structures. Indeed, when only a symmetric portion of the structure is analyzed, the nonlocal effects induced by the remaining (cut) portions are lost, this necessitates the consideration of an enlarged symmetric model on which appropriate nonlocal boundary conditions have to be imposed. It has to be pointed out how the width of such an enlarged model depends on the nonlocal material parameters, while the correct unknown nonlocal boundary conditions are here obtained and enforced by an iterative procedure. The accuracy of the proposed approach in solving nonlocal structural symmetric problems is tested with the aid of two numerical examples.     

Finite element solutions for nonhomogeneous nonlocal elastic problems

A finite element procedure for analysing nonhomogeneous nonlocal elastic 2D problems is presented and discussed. The procedure grounds on a variationally consistent approach known, in the relevant literature, as Nonlocal Finite Element Method. The latter is recast making use of a recently theorized phenomenological strain-difference-based nonhomogeneous nonlocal elastic model. The peculiarities of the numerical procedure together with the pertinent nonlocal operators are expounded and discussed. Two simple numerical 2D examples close the paper.     

A pseudospectral approach applicable for time integration of linearized N-S operator that removes pole singularity and physically spurious eigenmodes

This paper addresses a modified singularity removal technique for the eigenvalue or optimal mode problems in pipe flow using a pseudospectral method. The current approach results in the linear stability operator to be devoid of any unstable physically spurious modes, and thus, it provides higher numerical stability during time-based integration. The correctness of the numerical operator is established by calculating the known eigenvalues of pipe Poiseuille flow. Subsequently, the optimal modes are determined with Farrell's approach and compared with the existing literature. The usefulness of this approach is further demonstrated in the time-based numerical integration of the linearized Navier-Stokes operator for the adjoint method–based optimal mode determination. The numerical scheme is implemented with the radial velocity-radial vorticity formulation. Even number of Chebyshev-Lobatto grid points are distributed over the domain r∈[−1,1] omitting the centerline, which also efficiently provides higher resolution near the wall boundary. The boundary conditions are imposed with homogeneous wall boundary conditions, whereas the analytic nature of a proper set of base functions enforces correct centerline conditions. The resulting redundancy introduced in the process is eliminated with the proper usage of parity.     

Quadratic programming algorithm for wall slip and free boundary pressure condition

Wall slip is often observed in a highly sheared fluid film in a solid gap. This makes a difficulty in mathematical analysis for the hydrodynamic effect because fluid velocity at the liquid–solid interfaces is not known a priori. If the gap has a convergent–divergent wedge, a free boundary pressure condition, i.e. Reynolds pressure boundary condition, is usually used in the outlet zone in numerical solution. This paper, based on finite element method and parametric quadratic programming technique, gives a numerical solution technique for a coupled boundary non‐linearity of wall slip and free boundary pressure condition. It is found that the numerical error decreases with the number of elements in a negative power law having an index larger than 2. Our method does not need an iterative process and can simultaneously gives rise to fluid film pressure distribution, wall slip velocity and surface shear stress. Wall slip always decreases the hydrodynamic pressure. Large wall slip even causes a null hydrodynamic pressure in a pure sliding solid gap. Copyright © 2005 John Wiley & Sons, Ltd.     

On steady state computation of turbulent flows using k–ε models approximated by the time splitting method

The time splitting method is frequently used in numerical integration of flow equations with source terms since it allows almost independent programming for the source part. In this paper we will consider the question of convergence to steady state of the time splitting method applied to k–ε turbulence models. This analysis is derived from a properly defined scalar study and is carried out with success for the coupled k–ε equations. It is found that the time splitting method does not allow convergence to steady state for any choice of finite values of the time step. Numerical experiments for some typical turbulent compressible flow problems support the fact that the time splitting method is always nonconvergent, while its nonsplitting counterpart is convergent. Copyright © 2005 John Wiley & Sons, Ltd.     

二维非局部线弹性平面问题的辛分析

将二维非局部线弹性理论引入到Hamilton体系下,基于变分原理推导得出了二维线弹性理论的对偶方程和相应的边界条件.在分析验证对偶方程的准确性的基础上,该套方法被应用于二维弹性平面波问题的求解.将精细积分与扩展的W-W算法相结合在Hamilton体系下建立了求解平面Rayleigh波的数值算法.从推导到计算的保辛性确保了辛体系非局部理论与算法的准确性.通过对不同算例的数值计算,分析和对比了非局部理论方法与传统局部理论方法的差别,并进一步指出了该套算法的适用性和优势所在.     

3D BEM for general anisotropic elasticity

A new method for computing the system matrices G and H of the three-dimensional boundary element method is presented for the elastostatic problems of general anisotropy. Green’s functions, expressed by the line integrals over a semi-circle, have a simple structure that allows analytical integration over the boundary elements. The computation for the system matrices G and H is reduced to numerical integration of regular line integrals over the unit semi-circle. The proposed method has been implemented and tested with numerical examples.     

Experimental and numerical study of developing turbulent flow and heat transfer in convergent/divergent square ducts

 The work reported in this paper is a systematic experimental and numerical study of friction and heat transfer characteristics of divergent/convergent square ducts with an inclination angle of 1^∘ in the two direction at cross section. The ratio of duct length to average hydraulic diameter is 10. For the comparison purpose, measurement and simulation are also conducted for a square duct with constant cross section area, which equals to the average cross section area of the convergent/divergent duct. In the numerical simulation the flow is modeled as being three-dimensional and fully elliptic by using the body-fitted finite volume method and the kɛ turbulence model. The uniform heat flux boundary condition is specified to simulate the electrical heating used in the experiments. The heat transfer performance of the divergent/convergent ducts is compared with the duct with uniform cross section under three constraints (identical mass flow rate, pumping power and pressure drop). The agreement of the experimental and numerical results is quite good except at the duct inlet. Results show that for the three ducts studied there is a weak secondary flow at the cross section, and the circumference distribution of the local heat transfer coefficient is not uniform, with an appreciable reduction in the four corner regions. In addition, the acceleration/deceleration caused by the cross section variation has a profound effect on the turbulent heat transfer: compared with the duct of constant cross section area, the divergent duct generally shows enhanced heat transfer behavior, while the convergent duct has an appreciable reduction in heat transfer performance. Received on 18 September 2000 / Published online: 29 November 2001     

Buckling analysis of functionally graded nanobeams under non-uniform temperature using stress-driven nonlocal elasticity

In this work, the size-dependent buckling of functionally graded(FG)Bernoulli-Euler beams under non-uniform temperature is analyzed based on the stressdriven nonlocal elasticity and nonlocal heat conduction. By utilizing the variational principle of virtual work, the governing equations and the associated standard boundary conditions are systematically extracted, and the thermal effect, equivalent to the induced thermal load, is explicitly assessed by using the nonlocal heat conduction law. The ...     

Variation of the stress intensity factor along the front of a 3-D rectangular crack subjected to mixed-mode load

Summary  The singular integral equation method is applied to the calculation of the stress intensity factor at the front of a rectangular crack subjected to mixed-mode load. The stress field induced by a body force doublet is used as a fundamental solution. The problem is formulated as a system of integral equations with r ⁻³-singularities. In solving the integral equations, unknown functions of body-force densities are approximated by the product of polynomial and fundamental densities. The fundamental densities are chosen to express two-dimensional cracks in an infinite body for the limiting cases of the aspect ratio of the rectangle. The present method yields rapidly converging numerical results and satisfies boundary conditions all over the crack boundary. A smooth distribution of the stress intensity factor along the crack front is presented for various crack shapes and different Poisson's ratio. Received 5 March 2002; accepted for publication 2 July 2002     

A boundary element method for the solution of finite mobility ratio immiscible displacement in a Hele‐Shaw cell

In this paper, the interaction between two immiscible fluids with a finite mobility ratio is investigated numerically within a Hele‐Shaw cell. Fingering instabilities initiated at the interface between a low‐viscosity fluid and a high‐viscosity fluid are analysed at varying capillary numbers and mobility ratios using a finite mobility ratio model. The present work is motivated by the possible development of interfacial instabilities that can occur in porous media during the process of CO₂ sequestration but does not pretend to analyse this complex problem. Instead, we present a detailed study of the analogous problem occurring in a Hele‐Shaw cell, giving indications of possible plume patterns that can develop during the CO₂ injection. The numerical scheme utilises a boundary element method in which the normal velocity at the interface of the two fluids is directly computed through the evaluation of a hypersingular integral. The boundary integral equation is solved using a Neumann convergent series with cubic B‐Spline boundary discretisation, exhibiting sixth‐order spatial convergence. The convergent series allows the long‐term nonlinear dynamics of growing viscous fingers to be explored accurately and efficiently. Simulations in low‐mobility ratio regimes reveal large differences in fingering patterns compared with those predicted by previous high‐mobility ratio models. Most significantly, classical finger shielding between competing fingers is inhibited. Secondary fingers can possess significant velocity, allowing greater interaction with primary fingers compared with high‐mobility ratio flows. Eventually, this interaction can lead to base thinning and the breaking of fingers into separate bubbles. Copyright © 2015 John Wiley & Sons, Ltd.     

A numerical implementation using mid-node collocation for the hypersingular boundary contour method

A variant of the boundary element method, called the boundary contour method (BCM), offers a further reduction in dimensionality. Consequently, boundary contour analysis of two-dimensional problems does not require any numerical integration at all. In another development, a boundary contour implementation of a regularized hypersingular boundary integral equation (HBIE) using quadratic elements and end-node collocation was proposed and the technique is termed the hypersingular boundary contour method (HBCM). As reported in that work, the approach requires highly refined meshes in order to numerically enforce the stress continuity across boundary contour elements. This continuity requirement is very crucial since the regularized HBIE is only valid at collocation points where the stress tensor is continuous, while the computed stress at the endpoints of a boundary contour element, which is a non-conforming element, is generally not. This paper presents a new implementation of the HBCM for which the regularized HBIE is collocated at the mid-node of a boundary contour element. As the computed stress tensor is continuous at these mid-nodes, there is no need for unusually refined meshes. Some numerical tests herein show that, for the same mesh density, the HBCM using mid-node collocation has a comparable accuracy as the BCM.     

Solving high Reynolds‐number viscous flows by the general BEM and domain decomposition method

In this paper, the domain decomposition method (DDM) and the general boundary element method (GBEM) are applied to solve the laminar viscous flow in a driven square cavity, governed by the exact Navier–Stokes equations. The convergent numerical results at high Reynolds number Re = 7500 are obtained. We find that the DDM can considerably improve the efficiency of the GBEM, and that the combination of the domain decomposition techniques and the parallel computation can further greatly improve the efficiency of the GBEM. This verifies the great potential of the GBEM for strongly non‐linear problems in science and engineering. Copyright © 2004 John Wiley & Sons, Ltd.