A scaled boundary node method applied to two-dimensional crack problems

A boundary-type meshless method called the scaled boundary node method(SBNM) is developed to directly evaluate mixed mode stress intensity factors(SIFs) without extra post-processing.The SBNM combines the scaled boundary equations with the moving Kriging(MK) interpolation to retain the dimensionality advantage of the former and the meshless attribute of the latter.As a result,the SBNM requires only a set of scattered nodes on the boundary,and the displacement field is approximated by using the MK interpolation technique,which possesses the δ function property.This makes the developed method efficient and straightforward in imposing the essential boundary conditions,and no special treatment techniques are required.Besides,the SBNM works by weakening the governing differential equations in the circumferential direction and then solving the weakened equations analytically in the radial direction.Therefore,the SBNM permits an accurate representation of the singularities in the radial direction when the scaling center is located at the crack tip.Numerical examples using the SBNM for computing the SIFs are presented.Good agreements with available results in the literature are obtained.     

渗流自由面分析的比例边界有限元法

为确定渗流自由面,利用比例边界有限元法对二维稳定渗流场进行分析.通过比例坐标与直角坐标之间的转换,推导渗流问题的比例边界有限元方程.利用比例边界有限元空间降低一个维度,并在降维方向保持解析的特点,只对求解域边界进行离散.当比例中心选在坝体上游直线边界和坝体下面不可渗透直线边界交点时,只需离散自由面及其下游边界.并给出一种控制溢出点的控制点法.对二维坝体的稳态渗流场自由面问题进行分析并与实验结果比较.结论表明,该方法收敛快、结果较精确、数据准备工作量小.     

Torsional vibration of multi-step non-uniform rods with various concentrated elements

An analytical approach and exact solutions for the torsional vibration of a multi-step non-uniform rod carrying an arbitrary number of concentrated elements such as rigid disks and with classical or non-classical boundary conditions is presented. The exact solutions for the free torsional vibration of non-uniform rods whose variations of cross-section are described by exponential functions and power functions are obtained. Then, the exact solutions for more general cases, non-uniform rods with arbitrary cross-section, are derived for the first time. In order to simplify the analysis for the title problem, the fundamental solutions and recurrence formulas are developed. The advantage of the proposed method is that the resulting frequency equation for torsional vibration of multi-step non-uniform rods with arbitrary number of concentrated elements can be conveniently determined from a homogeneous algebraic equation. As a consequence, the computational time required by the proposed method can be reduced significantly as compared with previously developed analytical procedures. A numerical example shows that the results obtained from the proposed method are in good agreement with those determined from the finite element method (FEM), but the proposed method takes less computational time than FEM, illustrating the present methods are efficient, convenient and accurate.     

Acoustic scattering in dispersions: improvements in the calculation of single particle scattering coefficients

Measurements of ultrasound speed and attenuation can be related to the properties of dispersed systems by applying a scattering model. Rayleigh's method for scattering of sound by a spherical object, and its subsequent developments to include viscous, thermal, and other effects (known as the ECAH model) has been widely adopted. The ECAH method has difficulties, including numerical ill-conditioning, calculation of Bessel functions at large arguments, and inclusion of thermal effects in all cases. The present work develops techniques for improving the ECAH calculations to allow its use in instrumentation. It is shown that thermal terms can be neglected in some boundary equations up to approximately 100 GHz in water, and several simplified solutions result. An analytical solution for the zero-order coefficient is presented, with separate nonthermal and thermal parts, allowing estimation of the thermal contribution. Higher orders have been simplified by estimating the small shear contribution as the inertial limit is approached. The condition of the matrix solutions have been greatly improved by these techniques and by including appropriate scaling factors. A method is presented for calculating the required Bessel functions when the argument is large (high frequency or large particle size). The required number of partial wave orders is also considered.     

An Improved Formulation of Singular Boundary Method

This study proposes a new formulation of singular boundary method (SBM) to solve the 2D potential problems, while retaining its original merits being free of integration and mesh, easy-to-program, accurate and mathematically simple without the requirement of a fictitious boundary as in the method of fundamental solutions (MFS). The key idea of the SBM is to introduce the concept of the origin intensity factor to isolate the singularity of fundamental solution so that the source points can be placed directly on the physical boundary. This paper presents a new approach to derive the analytical solution of the origin intensity factor based on the proposed subtracting and adding-back techniques. And the troublesome sample nodes in the ordinary SBM are avoided and the sample solution is also not necessary for the Neumann boundary condition. Three benchmark problems are tested to demonstrate the feasibility and accuracy of the new formulation through detailed comparisons with the boundary element method (BEM), MFS, regularized meshless method (RMM) and boundary distributed source (BDS) method.     

Efficient enforcement of far-field boundary conditions in the Transformed Field Expansions method

The Method of Transformed Field Expansions (TFE) has been demonstrated to be a robust and highly accurate numerical scheme for simulating solutions of boundary value and free boundary problems from the sciences and engineering. As a Boundary Perturbation Method it builds highly accurate solutions based upon exact solutions in a simple, canonical, geometry and corrects these via Taylor series to fit the actual geometry at hand. The TFE method has significantly enhanced stability properties when compared with other Boundary Perturbation approaches, however, this comes at the cost of requiring a full volumetric discretization as opposed the surface formulation that other methods can realize. In this paper we outline two techniques for ameliorating this shortcoming, first by employing a Legendre Spectral Element Method to implement efficient, graded meshes, and second by utilizing an Artificial Boundary with a Transparent Boundary Condition placed quite close to the boundary of the domain. In this contribution we focus on the specific problem of simulating the Dirichlet–Neumann operator associated to Laplace’s equation on a periodic cell (which arises in the water wave problem). While the details of our results are specific to this problem, the general conclusions are valid for the wider class of problems to which the TFE method can be applied. For each technique we discuss implementation details and display numerical results which support the conclusion that each of these techniques can greatly reduce the computational cost of using the TFE method.     

Burton–Miller-type singular boundary method for acoustic radiation and scattering

This paper proposes the singular boundary method (SBM) in conjunction with Burton and Miller?s formulation for acoustic radiation and scattering. The SBM is a strong-form collocation boundary discretization technique using the singular fundamental solutions, which is mathematically simple, easy-to-program, meshless and introduces the concept of source intensity factors (SIFs) to eliminate the singularities of the fundamental solutions. Therefore, it avoids singular numerical integrals in the boundary element method (BEM) and circumvents the troublesome placement of the fictitious boundary in the method of fundamental solutions (MFS). In the present method, we derive the SIFs of exterior Helmholtz equation by means of the SIFs of exterior Laplace equation owing to the same order of singularities between the Laplace and Helmholtz fundamental solutions. In conjunction with the Burton–Miller formulation, the SBM enhances the quality of the solution, particularly in the vicinity of the corresponding interior eigenfrequencies. Numerical illustrations demonstrate efficiency and accuracy of the present scheme on some benchmark examples under 2D and 3D unbounded domains in comparison with the analytical solutions, the boundary element solutions and Dirichlet-to-Neumann finite element solutions.     

A Chebyshev–Lagrangian method for acoustic analysis of a rectangular cavity with arbitrary impedance walls

A general Chebyshev–Lagrangian method is proposed to obtain the analytical solution for a rectangular acoustic cavity with arbitrary impedance boundary conditions. The originality of the present paper is the successful attempt of applying orthogonal polynomials, such as Chebyshev polynomials of the first kind, to the analysis of a rectangular sound field with general wall impedance. The sound pressure is uniformly expressed as triplicate Chebyshev polynomial series which is independent in each direction. The Chebyshev polynomial series solution is obtained using the Rayleigh–Ritz procedure after considering the influence of boundary impedance on the cavity as the work done by the impedance surfaces in the Lagrangian function. The accuracy and reliability of the proposed method are validated against the analytical solutions and some numerical results available in the literature. Excellent orthogonality and complete properties of the Chebyshev polynomials ensure the rapid convergence, numerical stability, high accuracy of the current solution. The simplicity and low computational cost of the present approach make it preferable to obtain the results of complex models even in the relative high frequency range by choosing enough truncated terms in the sound pressure expression. Numerous cases with various uniform or non-uniform impedance boundary conditions are analyzed numerically and some of the results can be used as benchmark. It is shown that the impedance boundary condition can effectively influence or modify the acoustic characteristics and response of a cavity.     

Numerical simulation of the regularized long wave equation by He's homotopy perturbation method

In this Letter, we present the homotopy perturbation method (shortly HPM) for obtaining the numerical solution of the RLW equation. We obtain the exact and numerical solutions of the Regularized Long Wave (RLW) equation for certain initial condition. The initial approximation can be freely chosen with possible unknown constants which can be determined by imposing the boundary and initial conditions. Comparison of the results with those of other methods have led us to significant consequences. The numerical solutions are compared with the known analytical solutions.     

一种滑移区气体流动的格子Boltzmann曲边界处理新格式

针对滑移区复杂气-固边界存在速度滑移现象,提出了一种基于格子Boltzmann方法的非平衡态外推与有限差分相结合的曲边界处理新格式.该格式具有可考虑实际物理边界与网格线偏移量的优势,较传统half-way DBB(diffusive bounce-back)格式更能准确反映实际边界情况,同时还可获取壁面处气体宏观量及其法向梯度等信息.采用本文所提曲边界处理格式模拟分析了滑移区气体平直/倾斜微通道Poiseuille流、微圆柱绕流和同心微圆柱面旋转Couette流问题.研究结果表明,采用曲边界处理新格式所得结果与理论值以及文献结果符合良好,适用于滑移区气体流动的复杂边界处理,且比half-way DBB格式具有更高的精度,较修正DBB格式具有更好的适应性.     

Two analytic methods applied to radiative transfer in the linear-anisotropic scattering medium with graded index and diffuse gray boundaries

The curved ray-tracing method is extended to radiative transfer in the graded index medium with diffuse gray boundary conditions instead of black boundary conditions and the pseudo-source adding method is extended to the case of the linear-anisotropic scattering medium with graded index from non-scattering medium. Furthermore, the equivalence of the two methods is verified by formulation derivation. As exact analytical solutions, both the methods have high accuracy and fast computational speed. The predicted temperature distributions and dimensionless radiative heat flux at radiative equilibrium are determined by the proposed methods, and the numerical results are compared with the data in references. The results show that the present methods have a good accuracy. Influences of various combinations of refractive index and boundary emissivities on the temperature distributions and dimensionless radiative heat flux are also investigated.     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

Calculation of transient dynamic stress intensity factors at bimaterial interface cracks using a SBFEM-based frequency-domain approach

A frequency-domain approach based on the semi-analytical scaled boundary finite element method (SBFEM) was developed to calculate dynamic stress intensity factors (DSIFs) at bimaterial interface cracks subjected to transient loading. Because the stress solutions of the SBFEM in the frequency domain are analytical in the radial direction, and the complex stress singularity at the bimaterial interface crack tip is explicitly represented in the stress solutions, the mixed-mode DSIFs were calculated directly by definition. The complex frequency-response functions of DSIFs were then used by the fast Fourier transform (FFT) and the inverse FFT to calculate time histories of DSIFs. A benchmark example was modelled. Good results were obtained by modelling the example with a small number of degrees of freedom due to the semi-analytical nature of the SBFEM. Supported by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, Ministry of Education of China (Grant No. J20050924) and the Australian Research Council Discovery Project (Grant No. DP0452681)     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

Three Boundary Meshless Methods for Heat Conduction Analysis in Nonlinear FGMs with Kirchhoff and Laplace Transformation

This paper presents three boundary meshless methods for solving problems of steady-state and transient heat conduction in nonlinear functionally graded materials (FGMs). The three methods are, respectively, the method of fundamental solution (MFS), the boundary knot method (BKM), and the collocation Trefftz method (CTM) in conjunction with Kirchhoff transformation and various variable transformations. In the analysis, Laplace transform technique is employed to handle the time variable in transient heat conduction problem and the Stehfest numerical Laplace inversion is applied to retrieve the corresponding time-dependent solutions. The proposed MFS, BKM and CTM are mathematically simple, easy-to-programming, meshless, highly accurate and integration-free. Three numerical examples of steady state and transient heat conduction in nonlinear FGMs are considered, and the results are compared with those from meshless local boundary integral equation method (LBIEM) and analytical solutions to demonstrate the efficiency of the present schemes.     

Lattice Boltzmann modeling of microchannel flow in slip flow regime

We present the lattice Boltzmann equation (LBE) with multiple relaxation times (MRT) to simulate pressure-driven gaseous flow in a long microchannel. We obtain analytic solutions of the MRT-LBE with various boundary conditions for the incompressible Poiseuille flow with its walls aligned with a lattice axis. The analytical solutions are used to realize the Dirichlet boundary conditions in the LBE. We use the first-order slip boundary conditions at the walls and consistent pressure boundary conditions at both ends of the long microchannel. We validate the LBE results using the compressible Navier–Stokes (NS) equations with a first-order slip velocity, the information-preservation direct simulation Monte Carlo (IP-DSMC) and DSMC methods. As expected, the LBE results agree very well with IP-DSMC and DSMC results in the slip velocity regime, but deviate significantly from IP-DSMC and DSMC results in the transition-flow regime in part due to the inadequacy of the slip velocity model, while still agreeing very well with the slip NS results. Possible extensions of the LBE for transition flows are discussed.     

Solutions of Heat-Like and Wave-Like Equations with Variable Coefficients by Means of the Homotopy Analysis Method

We employ the homotopy analysis method （HAM） to obtain approximate analytical solutions to the heat-like and wave-like equations. The HAM contains the auxiliary parameter h, which provides a convenient way of controlling the convergence region of series solutions. The analysis is accompanied by several linear and nonlinear heat-like and wave-like equations with initial boundary value problems. The results obtained prove that HAM is very effective and simple with less error than the Adomian decomposition method and the variational iteration method.     

Plane wave interpolation in direct collocation boundary element method for radiation and wave scattering: numerical aspects and applications

The classical boundary element formulation for the Helmholtz equation is rehearsed, and its limitations with respect to the number of variables needed to model a wavelength are explained. A new type of interpolation for the potential is then described in which the usual boundary element shape functions are modified by the inclusion of a set of plane waves, propagating in a range of directions. This is termed the plane wave basis boundary element method. The modifications needed to the classical procedures, in terms of integration of the element matrices, and location of collocation points are described. The well-known Singular Value Decomposition solution technique, which is adopted here for the solution of the system matrix equation in its complex form, is briefly outlined. The conditioning of the system matrix is analysed for a simple radiation problem. The corresponding diffraction problem is also analysed and results are compared with analytical and classical boundary element solutions. The CHIEF method is adopted to enhance the quality of the solution, particularly in the vicinity of irregular frequencies. The plane wave basis boundary element method is then applied to two problems: scattering of plane waves by an elliptical cylinder and the multiple circular cylinder plane wave scattering problem. In both cases results are compared with analytical solutions. The results clearly demonstrate that the new method is considerably more efficient than the classical approach. For a given number of degrees of freedom, the frequency for which accurate results can be obtained, using the new technique, can be up to three or four times higher than that of the classical method. This makes the method a powerful new addition to our tools for tackling high-frequency radiation and scattering problems.     

电力线谐波辐射在分层各向异性电离层中的传播特点

电力线谐波辐射特指在电离层或磁层中观测到的来源于地面电力系统输电线的电磁波辐射, 其在电磁场时频功率谱中表现为400 Hz至5 kHz范围内, 频率间隔为50/100 Hz或60/120 Hz 的平行谱线, 已成为近地空间环境的一种人为污染源. 对于该现象的形成机理尚缺乏定量研究. 本文研究了非理想导电大地上方由电偶极子源产生的电磁场在分层各向异性电离层中的传播模型, 提出了一种新的求解方法, 有效解决了编程计算中的数值溢出问题, 并利用已有解析解对所提方法进行了验证. 在此基础上, 利用实际电力线、大地、电离层的相关参数, 研究了偶极子源频率、电离层下边界高度、大地电导率、地磁场方向等对电力线谐波辐射在电离层中的传播的影响. 结果表明, 频率等于地-电离层波导导波模截止频率时透入电离层的电力线谐波辐射强度更大; 谐波电流一定时, 大地电导率小的地区, 电力线谐波辐射的功率更大; 电力线谐波辐射在电离层中沿地磁场方向传播. 本文所得结果有益于阐释电力线谐波辐射现象的形成机理.