三维高频声波的矩阵压缩边界节点法模拟

本文采用边界节点法(Boundary Knot Method, BKM)求解三维高频声场．由于高频赫姆霍兹方程的解是振荡的，极大影响了数值求解的精确度，需要在计算区域增加离散点，这会增加计算量．同时对于大规模声学问题，依靠边界节点法形成的插值矩阵为满秩，导致计算量过高和存储量过大．所以本文采用矩阵压缩技术(Matrix Compression, MC)，在有效继承边界节点法高精确度的基础上减少计算内存需求和时间，从而提高计算效率．数值实验表明，MC-BKM 求解精度高、收敛速度快、计算时间少，在高频大规模声波问题中应用前景广泛．     

快速求解叶栅流动的边界元方法

本文利用边界元方法，通过ｋｒｉｃｈｈｏｆｆ变换将描述叶栅流动的控制方程转换成线性方程。并将广义ｋ－Ｊ条件与边界积分方程联立求解，避免了非线性项和叶片出气角的迭代计算。完成了一种快速求解任意迥转面叶栅流场的计算程序，实用表明与其它数值方法及实验结果符合较好，具有快带、简明、实用的特点。     

pFFT快速边界元方法模拟三维声散射

研究了用pFFT快速边界元方法模拟声散射问题的关键技术。采用Burton—Miller方程消除了声学边界元方法中外问题解的不唯一现象。为此,文中研究了采用常量元时该方程中超奇异积分的计算方法。最后,通过对平面声波的刚性圆球声散射的数值模拟,验证了建立的声学pFFT快速边界元方法。     

非连续边界元——有限元耦合方法分析

对边界元－有限元耦合方法进行了分析，采用非连续元离散边界积分方程，解决了耦合分析中自由度约束问题，给出了非连续边界元同有限元耦合的具体实施步骤，通过对二维弹性力学和流＝固耦合问题分析，表明了该文方法的有效性。     

基于行消元回代法的多域边界元分析方法

提出了一种用多域边界元技术求解大型工程问题的新算法. 首先, 采用三步变量凝聚技术, 将由内部点、边界点和公共结点表述的每一子域的基本边界元代数方程表述成只有公共结点变量为未知量的代数方程, 然后, 根据公共结点的平衡方程和协调条件组集具有稀疏系数特征的总体系统方程组. 为了有效求解该系统方程组, 首次在边界元法中引进一种能有效求解大型非对称稀疏系数矩阵方程组的行消元回代法(REBSM), 该方法可在方程的每一行组形成时进行消元和回代, 当方程组组集完毕后即可得到方程的解, 不需要最后的回代过程. 因为一些项的重复计算在每一行的处理中合并掉, 因此REBSM要比传统的高斯消元法需要较少的内存, 而且计算速度具有数量级的提高, 可为边界元法求解大型工程问题提供有力的方程求解器.      

边界元技术中的全特解场方法

采用相应基本解来构造全特解场，据而确定各该问题边界元方程系数矩阵的全部元素。解算中不涉及具体插值，不用数值积分，避免了奇异积分的数值处理，而任意点的物理量计算（如应力、声强等等）不依赖于待解的边界未知量，因而算效奇高，精度特好。适于求解各类数学物理和工程技术问题。     

边界元的模态综合技术

1 前言基于质量矩阵的边界元自由振动分析首先由Nardini D.和Brebbia C.A.于1983年给出.以后Ahmad S.和Banerjee P.K.于1986年做了类似的工作.边界元方法由于仅需要处理结构的边界信息,较之区域性方法有明显的优点,可以得到高的计算效率和精度.但由于边界元方法导出的系统方程,其系数矩阵为非对称的满阵,方程阶数较高时,求解费用大.为克服这一困难,可采用有限元中较为成熟的模态综合技术.为了使边界元动力方程解耦,引入了左特征向量和右特征向量,并在此基础上导出了     

边界元技术新探——边界点方法

简要回顾边界元方法,重点介绍边界点方法及其展望,此法的特色是精度和效率相统一,好懂易用,因而将更有利于向工程界推广．     

确定裂纹体等效弹性模量的边界元方法

采用边界元方法计算含有序分布裂纹的裂纹体在压缩载荷作用下的等效弹性模量，利用一种能适当考虑裂纹有间相作用的自洽理论，建立了相应的迭代格式，通过算例研究了裂纹方向，裂纹面间摩擦系数对裂纹体等效弹性模量的影响。     

A Symmetric Galerkin BEM for Harmonic Problems and Multiconnected Bodies

A steady state harmonic solution of Navier equations of motion is presented using the symmetric Galerkin Boundary Element Method, which leads to a symmetric system of equations. For the integration, a two-step regularisation process is proposed to deal with the singularities of the involved kernels. It is also shown that the standard application of the symmetric Galerkin methodology for multiply connected bodies leads to the deterioration of the conditioning of the system for low frequencies. Two examples are presented comparing the solutions obtained with the proposed formulation with those obtained with the standard collocation method. A third example is presented to illustrate the behaviour of a simple multi-connected body.     

Applications of wavelet Galerkin FEM to bending of plate structure

A kind of calculating method for high order differentials expanded by the wavelet scaling functions and the integral of their product used in Galerkin FEM is proposed, so that we can use the wavelet Galerkin FEM to solve boundary-value differential equations with orders higher than two. To combine this method with the Generalized Gaussian integral method in wavelet theory, we can find that this method has many merits in solving mechanical problems, such as the bending of plates and those with variable thickness. The numerical results show that this method is accurate. This project is supported by the National Natural Science Foundation of China (No. 19772014) and the National Outstanding Young Scientist Foundation of China (No. 19725207).     

Divergence‐free discontinuous Galerkin schemes for the Stokes equations and the MAC scheme

We show that recently studied discontinuous Galerkin discretizations in their lowest order version are very similar to the marker and cell (MAC) finite difference scheme. Indeed, applying a slight modification, the exact MAC scheme can be recovered. Therefore, the analysis applied to the DG methods applies to the MAC scheme as well and the DG methods provide a natural generalization of the MAC scheme to higher order and irregular meshes. Copyright © 2007 John Wiley & Sons, Ltd.     

Modified extended BDF scheme for the discontinuous Galerkin solution of unsteady compressible flows

In this paper, a high‐order DG method coupled with a modified extended backward differentiation formulae (MEBDF) time integration scheme is proposed for the solution of unsteady compressible flows. The objective is to assess the performance and the potential of the temporal scheme and to investigate its advantages with respect to the second‐order BDF. Furthermore, a strategy to adapt the time step and the order of the temporal scheme based on the local truncation error is considered. The proposed DG‐MEBDF method has been evaluated for three unsteady test cases: (i) the convection of an inviscid isentropic vortex; (ii) the laminar flow around a cylinder; and (iii) the subsonic turbulent flow through a turbine cascade. Copyright © 2014 John Wiley & Sons, Ltd.     

A discontinuous Galerkin scheme for the numerical solution of flow problems with discontinuities

In this paper, a new discontinuous Galerkin finite element method for the numerical solution of flow problems with discontinuities is presented. The method is based on the limitation in every cell of the difference between the extrema values and the mean value of the numerical solution. The algorithm and technical details for the implementation of the method are presented in one‐and two‐dimensional problems. Numerical experiments for classical test problems are solved on unstructured triangulations to demonstrate the performance of the proposed method. Copyright © 2011 John Wiley & Sons, Ltd.     

一种新的摄动粘弹性边界元法

本文在拉普拉斯变换空间中,运用摄动理论和“弹性-粘弹性”对应性原理,提出了一种新的粘弹性问题边界元求解法。文中结合广义变分原理,导出了摄动粘弹性边界元方程,并给出了摄动粘弹性问题的基本解。最后,详述了一阶摄动的求解过程,并给出了算例。     

基于边界元方法的气动弹性降阶模型及仿真

传统气动弹性的时域计算耗费了大量时间,为了提高计算效率,本文发展了基于边界元方法的降阶模型技术。首先基于边界元方法建立非定常流场的求解模型,结合特征值分析技术建立了非定常气动力的低阶模型;然后,利用边界元方法建立了气动网格和结构网格之间的信息转换矩阵;最后将非定常气动力降阶模型和结构动力学方程联合,建立了气动弹性系统的低阶状态空间模型。将所发展的降阶模型方法应用于NACA0012翼型的非定常气动力求解中,结果表明降阶模型可以在保证原系统计算精度的同时提高了计算效率;将降阶模型技术应用到三维机翼的气动弹性响应计算中,在系统阶数仅为12阶的情况下可以得到与原系统一致的极限环响应,说明降阶模型技术在求解气动弹性问题中的巨大优势。     

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.     

Revisit of dilation-based shock capturing for discontinuous Galerkin methods

The idea of using velocity dilation for shock capturing is revisited in this paper, combined with the discontinuous Galerkin method. The value of artificial viscosity is determined using direct dilation instead of its higher order derivatives to reduce cost and degree of difficulty in computing derivatives. Alternative methods for estimating the element size of large aspect ratio and smooth artificial viscosity are proposed to further improve robustness and accuracy of the model. Several benchmark tests are conducted, ranging from subsonic to hypersonic flows involving strong shocks. Instead of adjusting empirical parameters to achieve optimum results for each case, all tests use a constant parameter for the model with reasonable success, indicating excellent robustness of the method. The model is only limited to third-order accuracy for smooth flows. This limitation may be relaxed by using a switch or a wall function. Overall, the model is a good candidate for compressible flows with potentials of further improvement.     

A posteriori error of a discontinuous Galerkin scheme for compressible miscible displacement problems with molecular diffusion and dispersion

A kind of compressible miscible displacement problems including molecular diffusion and dispersion in porous media is studied. A symmetric interior penalty discontinuous Galerkin method is applied to the coupled system of flow and transport. Explicit a posteriori error estimates are obtained based on the duality argument, approximation properties and residual notations. The resulting error bound is shown to be reliable and worthwhile for guiding dynamic mesh adaptivity. Copyright © 2009 John Wiley & Sons, Ltd.