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

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

密度矩阵重正化群的异构并行优化

密度矩阵重正化群方法（DMRG）在求解一维强关联格点模型的基态时可以获得较高的精度,在应用于二维或准二维问题时,要达到类似的精度通常需要较大的计算量与存储空间.本文提出一种新的DMRG异构并行策略,可以同时发挥计算机中央处理器（CPU）和图形处理器（GPU）的计算性能.针对最耗时的哈密顿量对角化部分,实现了数据的分布式存储,并且给出了CPU和GPU之间的负载平衡策略.以费米Hubbard模型为例,测试了异构并行程序在不同DMRG保留状态数下的运行表现,并给出了相应的性能基准.应用于4腿梯子时,观测到了高温超导中常见的电荷密度条纹,此时保留状态数达到10⁴,使用的GPU显存小于12 GB.     

基于等几何分析的边界元法求解Helmholtz问题

将基于一类局部双变量B样条函数的等几何分析方法和Burton-Miller方法相结合,分析三维Helmholtz问题.对于某些从二维参数域映射到三维空间具有奇异点的参数曲面,该方法可以有效地避免奇异点处大量奇异与近奇异积分的计算.数值算例表明该方法具有较好的计算精度和计算效率.复杂问题的分析表明,该方法具有良好的工程应用前景.     

比例边界有限元法在电机电磁场问题中的应用

将比例边界有限元法应用到电机磁场分析计算中,采用加权余量法推导无旋磁场的比例边界有限元方程及其离散格式,给出相邻子域公共边上内部节点磁位求解方法.以具有解析解的电机磁场问题为例进行计算和对比分析,验证离散格式的正确性和比例边界有限元法求解电机磁场的可行性,为该方法在电机磁场问题中的应用打下基础.     

Modelling of water wave interaction with multiple cylinders of arbitrary shape

This paper describes the development of an efficient numerical model, namely scaled boundary finite-element method (SBFEM) for linear waves interaction with cylindrical structures of arbitrary shapes. The two-dimensional Helmholtz equation is firstly weakened in the circumferential direction, so that the governing partial differential equation is transformed to an ordinary matrix differential equation in radial direction, and is solved fully analytically. As a key element, a virtual porous circular cylinder surrounding the cylindrical structures is introduced so that the entire computational domain is partitioned along the virtual cylinder into an unbounded and several bounded sub-domains with common interfaces. The principle innovation is that, the present SBFEM model chooses Hankel function as a base solution for the unbounded sub-domain, while a power series is used for the internal bounded sub-domains. The approach discretises only the common interfaces of the sub-domains with surface finite-elements, and fewer elements are required to obtain very accurate results. Numerical simulations show that the new SBFEM model offers a considerable improvement by far in its numerical performance, as well as in the range of physical phenomena that is capable of simulating. The wave forces and run-ups are presented for a single and multiple cylindrical structures of different cross sectional shapes. Influences of the incident wave parameters and structural configurations on the hydrodynamics are examined.     

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)     

Isogeometric finite element analysis of interior acoustic problems

Isogeometric Analysis (IGA) can bridge the gap between geometrical and numerical modeling. To this end, the same basis functions used in Computer Aided Design are applied to represent geometry and approximate physical field in analysis. In this paper, the IGA is firstly introduced to finite element method (FEM) for interior acoustic problems. The domain is parameterized by Non-Uniform Rational B-Spline (NURBS) in the algorithm, which simplifies the mesh generation greatly and furthermore supplies an exact representation of curved boundaries. In addition, the IGA FEM possesses a distinct feature of flexible order-elevation technique without modifying the geometry. Several numerical examples are presented to validate the accuracy and demonstrate the merits of the IGA FEM in the analysis of interior acoustic problems.     

虚边界元法与强流电子枪的电子轨迹模拟

以位势问题为分析对象,从格林公式出发严格导出了虚边界元法的基本积分方程.并以虚边界元法基本积分方程为出发点,从理论上导出了虚边界元法对虚边界形状位置的基本要求.理论结果表明,虚边界元法本身是一种严格的方法,但在一定条件下,它是一种近似方法.基于虚边界元方法开发了一套适用于强流电子枪的电子轨迹模拟程序,并对强流电子枪进行了计算.对强流电子枪的计算表明,在同等条件下,取定最佳的虚实边界距离与边界单元数,虚边界元法较边界元法的精度要高.     

Alternating Direction Implicit Orthogonal Spline Collocation on Non-Rectangular Regions

The alternating direction implicit (ADI) method is a highly efficient technique for solving multi-dimensional time dependent initial-boundary value problems on rectangles. When the ADI technique is coupled with orthogonal spline collocation (OSC) for discretization in space, we not only obtain the global solution efficiently, but the discretization error with respect to space variables can be of an arbitrarily high order. In [2], we used a Crank Nicolson ADI OSC method for solving general nonlinear parabolic problems with Robin's boundary conditions on rectangular polygons and demonstrated numerically the accuracy in various norms. A natural question that arises is: Does this method have an extension to non-rectangular regions? In this paper, we present a simple idea of how the ADI OSC technique can be extended to some such regions. Our approach depends on the transfer of Dirichlet boundary conditions in the solution of a two-point boundary value problem (TPBVP). We illustrate our idea for the solution of the heat equation on the unit disc using piecewise Hermite cubics.     

High-order compact schemes for incompressible flows: A simple and efficient method with quasi-spectral accuracy

In this paper, a finite difference code for Direct and Large Eddy Simulation (DNS/LES) of incompressible flows is presented. This code is an intermediate tool between fully spectral Navier–Stokes solvers (limited to academic geometry through Fourier or Chebyshev representation) and more versatile codes based on standard numerical schemes (typically only second-order accurate). The interest of high-order schemes is discussed in terms of implementation easiness, computational efficiency and accuracy improvement considered through simplified benchmark problems and practical calculations. The equivalence rules between operations in physical and spectral spaces are efficiently used to solve the Poisson equation introduced by the projection method. It is shown that for the pressure treatment, an accurate Fourier representation can be used for more flexible boundary conditions than periodicity or free-slip. Using the concept of the modified wave number, the incompressibility can be enforced up to the machine accuracy. The benefit offered by this alternative method is found to be very satisfactory, even when a formal second-order error is introduced locally by boundary conditions that are neither periodic nor symmetric. The usefulness of high-order schemes combined with an immersed boundary method (IBM) is also demonstrated despite the second-order accuracy introduced by this wall modelling strategy. In particular, the interest of a partially staggered mesh is exhibited in this specific context. Three-dimensional calculations of transitional and turbulent channel flows emphasize the ability of present high-order schemes to reduce the computational cost for a given accuracy. The main conclusion of this paper is that finite difference schemes with quasi-spectral accuracy can be very efficient for DNS/LES of incompressible flows, while allowing flexibility for the boundary conditions and easiness in the code development. Therefore, this compromise fits particularly well for very high-resolution simulations of turbulent flows with relatively complex geometries without requiring heavy numerical developments.     

Stress analysis of multiple anisotropic elliptical inclusions in composites

A volume integral equation method is introduced for the solution of elastostatic problems in an unbounded isotropic elastic solid containing interacting multiple anisotropic elliptical inclusions subject to uniform remote tension or remote in-plane shear. This method is applied to two-dimensional problems involving long parallel elliptical cylindrical inclusions. A detailed analysis of the stress field at the interface between the matrix and the central elliptical inclusion is carried out for square and hexagonal packing of anisotropic inclusions. The effects of the number of anisotropic inclusions and various inclusion volume fractions on the stress field at the interface between the isotropic matrix and the central elliptical cylindrical inclusion are investigated in detail. The stress field at the interface between the isotropic matrix and the central elliptical inclusion is also compared with that between the isotropic matrix and the central circular inclusion.     

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.     

On Solution of the Integrable Initial Boundary Value Problem for KdV Equation on the Semi-axis

This paper is concerned with the Korteweg-de Vries (KdV) equation on the right semi-axis. The initial boundary value problem with inhomogeneous integrable boundary conditions is studied. We show that, under some conditions on the initial data the problem has a solution and provide the procedure for constructing this solution. The procedure is based on the inverse spectral method and consists of several steps reducing to either solving some linear problems or calculations via some explicit formulas.     

双曲方程外问题的完全无反射人工边界条件

研究了数值求解双曲方程外问题的人工边界方法.在圆形的人工边界上,得到了三类等价的完全无反射的人工边界条件.给出的数值例子验证了方法的有效性.     

Natural Convection in a Concentric Annulus: A Lattice Boltzmann Method Study with Boundary Condition-Enforced Immersed Boundary Method

In this paper, a boundary condition-enforced IBM is introduced into the LBM in order to satisfy the non-slip and temperature boundary conditions, and natural convections in a concentric isothermal annulus between a square outer cylinder and a circular inner cylinder are simulated. The obtained results show that the boundary condition-enforced method gives a better solution for the flow field and the complicated physics of the natural convections in the selected case is correctly captured. The calculated average Nusselt numbers agree well with the previous studies.     

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.     

A solution of a non-homogeneous orthotropic cylindrical shell for axisymmetric plane strain dynamic thermoelastic problems

A solution of a non-homogeneous orthotropic elastic cylindrical shell for axisymmetric plane strain dynamic thermoelastic problems is developed. Firstly, a new dependent variable is introduced to rewrite the governing equation, the boundary conditions as well as the initial conditions. Secondly, a special function is introduced to transform the inhomogeneous boundary conditions to the homogeneous ones. Then by virtue of the orthogonal expansion technique, the equation with respect to the time variable is derived, of which the solution can be obtained. The displacement solution is finally presented, which can degenerate in a rather straightforward way to the solution for a homogeneous orthotropic cylindrical shell and isotropic solid cylinder as well as that for a non-homogeneous isotropic cylindrical shell. Using the present method, integral transform can be avoided. It is fit for a cylindrical shell with arbitrary thickness subjected to arbitrary thermal loads. It is also very convenient to deal with dynamic thermoelastic problems for different boundary conditions. Besides, the numerical calculation involved is very easy to be performed. Several examples are presented.     

The R-function method for the free vibration analysis of thin orthotropic plates of arbitrary shape

Free flexural vibrations of homogeneous, thin, orthotropic plates of an arbitrary shape with mixed boundary conditions are studied using the R-function method. The proposed method is based on the use of the R-function theory and variational methods. In contrast to the widely used methods of the network type (finite differences, finite element, and boundary element methods), in the R-function method all the geometric information given in the boundary value problem statement is represented in an analytical form. This allows one to seek a solution in a form of some formulas called a solution structure. These solution structures contain some indefinite functional components that can be determined by using any variational method. A method of constructing the solution structures satisfying the required mixed boundary conditions for eigenvalue plate bending problems is described. Numerical examples for the vibration analysis of orthotropic plates of complex geometry with mixed boundary conditions for illustrating the aforementioned R-function method and comparison against the other methods are made to demonstrate its merits.     

Exact sum rules for inhomogeneous drums

We derive general expressions for the sum rules of the eigenvalues of drums of arbitrary shape and arbitrary density, obeying different boundary conditions. The formulas that we present are a generalization of the analogous formulas for one dimensional inhomogeneous systems that we have obtained in a previous paper. We also discuss the extension of these formulas to higher dimensions. We show that in the special case of a density depending only on one variable the sum rules of any integer order can be expressed in terms of a single series. As an application of our result we derive exact sum rules for the homogeneous circular annulus with different boundary conditions, for a homogeneous circular sector and for a radially inhomogeneous circular annulus with Dirichlet boundary conditions.     

轴对称回转体瞬态温度场边界元分析

本文从三维瞬态势问题的边界积分方程及基本解出发,推导出轴对称势问题的边界积分方程及基本解,然后离散形成边界元方程,在此基础上对若干瞬态温度场数值例进行了分析计算,结果表明推导出的方法是可行的,并具有较高的精度和稳定性,能应用于工程中复杂的轴对称回转体的瞬态温度场分析。