改进的无奇异局部边界积分方程方法

在局部边界积分方程方法中,当源节点位于分析域的整体边界上时,局部边界积分将出现奇异积分问题,这些奇异积分需要做特别的处理．为此,提出了对域内节点采用局部积分方程,而对边界节点直接采用移动最小二乘近似函数引入边界条件来解决奇异积分问题,这同时也解决了对积分边界进行插值引入近似误差的问题．作为应用和数值实验,对Laplace方程和Helmholtz方程问题进行了分析,取得了很好的数值结果．进而,在Helmholtz方程求解中,采用了含波解信息的修正基函数来代替单项式基函数进行近似．数值结果显示,这样处理是简单高效的,在高波数声传播问题的求解中非常具有前景．     

薄体结构温度场的高阶边界元分析

分析了二维问题边界元法3节点二次单元的几何特征,区分和定义了源点相对高阶单元的Ⅰ型和Ⅱ型接近度.针对二维位势问题高阶边界元中奇异积分核,构造出具有相同Ⅱ型几乎奇异性的近似核函数,在几乎奇异积分单元上分离出积分核中主导的奇异函数部分.原积分核扣除其近似核函数后消除几乎奇异性,成为正则积分核函数,并采用常规Gauss数值方法计算该正则积分;对奇异核函数的积分推导出解析公式,从而建立了一种新的边界元法高阶单元几乎奇异积分半解析算法.应用该算法计算了二维薄体结构温度场算例,计算结果表明高阶单元半解析算法能充分发挥边界元法优势,显著提高计算精度.     

矩形容器中黏性流体的横波谐振:格子Boltzmann浸没边界方法（英文）

将三维格子Boltzmann法(LBM)与浸没边界法(IBM)相结合,研究弹性矩形容器中黏性流体的横波谐振所引起的流动物理特性.提出了一个半微观表达式来计算边界节点处的流体受力.基于薄板弹性变形理论,使用解析变形解法来计算边界所经历的位移.基于该方法的数值模拟结果与固定边界的理论预测结果一致.采用振荡边界模拟展现了与理论预期相符合的流动模式.     

基于Gauss全局径向基函数的近岸浅水变形波高数值计算新方法

应用Gauss全局径向基函数来模拟波浪浅水变形波高变化方程中的未知函数,经实例分析探讨得到了一种可用于求解该方程数值解的新方法,并将其计算结果与常用数值分析方法得到的数值解相互对比印证,证明了基于Gauss全局径向基函数法计算结果的正确性.经验证,Gauss径向基函数法的平均计算误差相比其他方法均要小,表明该方法拥有更高的计算精度.同时,根据Gauss全局径向基函数的逼近结果,得出了浅水变形波高变化微分方程数值解的拟合函数,在实际工程中,可以利用该拟合函数来代替原方程的解析解,研究成果可为求解近岸浅水区域波浪运动提供一种新思路.     

奇异边界法分析含水下障碍物水域中的水波传播问题

研究奇异边界法模拟水波在含水下障碍物水域的传播过程.奇异边界法是一种最近提出的新型边界配点方法,具有无网格和无数值积分、数学简单、编程容易等优点.首先研究了奇异边界法分析典型水波算例的精度及效率,并与边界元法的计算结果进行比较,然后通过数值模拟讨论分析了水下障碍物位置、尺寸及形状等因素对水波传播的影响.发现奇异边界法的计算精度较高,且与边界元法的计算结果吻合较好;数值结果显示水下障碍物的不同高宽比对水波的传播影响明显:障碍物无量纲高度越大对水波的屏障作用越明显;障碍物无量纲宽度增加对水波的屏障作用先增强后变弱.在高宽比一定时,斜率变化对水波的屏障作用不明显;含吸收边界水下障碍物可以得到较低的传递系数和较高的反射系数, 对水波的屏障作用更为明显.     

瞬态热传导的奇异边界法及其MATLAB实现

基于动力学问题时间依赖基本解的奇异边界法是一种无网格边界配点法.该方法引入源点强度因子的概念从而避免了基本解的源点奇异性,具有数学简单、编程容易、精度高等优点.将该方法用于瞬态热传导问题的数值模拟,运用MATLAB实现该问题的数值研究,并创建相应的MATLAB工具箱.针对二维和三维瞬态热传导问题,进行了基于反插值技术和经验公式的奇异边界法MATLAB算例实现.针对支撑圆坯低温瞬态温度场的模拟结果表明,瞬态热传导奇异边界法的MATLAB工具箱具有简单、方便、精确可靠的优点.研究成果有助于发展瞬态热传导的奇异边界法,并为瞬态热传导问题的数值分析和仿真提供了一种简单高效的模拟工具.     

半无限平面裂纹构型横向应力的Green函数

针对各向同性弹性无限大板中半无限裂纹,用解析函数方法求解了裂尖处横向应力的Green函数.加载情况为一任意集中力作用于任意一内点处.用叠加法求解了复势,它给出该平面问题的弹性解.通过渐近分析抽取复势的非奇异部分.基于该非奇异部分,用一种直接方法求解了横向应力的Green函数.进一步,用叠加法得到了一对对称和反对称集中力加载时的Green函数.然后,用得到的Green函数来预测铁电材料双悬臂梁试验中畴变引起的横向应力.用力电联合加载引起的横向应力来判断试验中所观察到的稳定和不稳定裂纹扩展行为.预测结果和试验数据基本吻合.     

板弯曲断裂计算的边界配置解法

采用复变函数理论和边界配置方法,分析计算了Kirchhoff板的弯曲断裂问题.假设了位移及内力的复变函数式,它们能满足一系列的基本方程和支配条件,例如域内的平衡方程、裂纹表面的边界条件、裂纹尖端的应力奇异性质.这样,仅板边界的边界条件需要考虑.它们可用边界配置法和最小二乘法近似满足.对不同边界条件和载荷情形进行了分析计算.数值算例表明,本文方法精度较高,计算量小,是一种有效的半解析、半数值计算方法.     

指数时程差分与有理谱配点法求解奇异摄动Burgers-Huxley问题

带小参数ε的Burgers-Huxley方程是一类非线性、非定常奇异摄动初边值问题,本文用指数时程差分与有理谱配点法求其数值解.对空间方向的边界层,用带sinh变换的有理谱配点法便Chebyshev节点在边界层处加密,只需取较少节点即可达到较高精度;时间方向采用指数时程差分与4阶Runge-Kutta法相结合的格式,并用围线积分计算矩阵甬数的方法克服了求解奇异摄动问题时遇到的的数值不稳定堆题.数值实验表明,本文提出的方法在求解左、右边界层和内部层的奇异摄动Burgers-Huxley问题都有较高的精度.     

非连续变形分析的MLS后处理方法

非连续变形分析(discontinuous deformatrion analysis, DDA)通过引入虚拟节理网格将块体离散成子块体系统进行断裂扩展数值模拟.针对这种方法难以获得精确块体应力分布的问题, 提出一种基于无网格法移动最小二乘(moving least squares, MLS)插值的应力恢复算法.利用DDA计算得到的节点位移, 通过恰当构造MLS形函数及其导数, 推导了块体任意点应力的计算公式.数值算例将基于MLS后处理的结果与解析解及平均值法后处理结果进行比较, 验证了所提出方法的精确性和有效性.     

用时变雷诺方程模型模拟孤立波与半圆型防波堤的相互作用

以时变雷诺方程为控制方程,用k-ε模型封闭该方程,采用体积函数(VOF)方法来跟踪波动自由表面,建立了二维垂向波浪数学模型,并用已有的实验资料进行了验证．随后用该模型模拟了半圆型防波堤与孤立波在淹没、平顶水位、完全露顶且不越浪3种典型工况下的相互作用过程．得到了半圆堤附近的流场、压强场和波面的变形过程．结果表明,在淹没状态下,半圆堤背浪面的底部会产生涡旋;平顶水位时,由于越浪的冲击作用,在半圆堤的背浪面将逐渐形成一对较大的涡旋,而半圆堤背浪面的底部,速度始终相对较小;而在露顶不越浪时,半圆堤的迎浪面会出现波浪的二次爬升的现象．为进一步研究结构物附近的污染物的输移扩散和泥沙运动提供基础．     

Self-adaptive FEM numerical modeling of the mild-slope equation

Based on the linear wave theory, the mild-slope equation (MSE) is a preferred mathematical model to simulate nearshore wave propagation. A numerical model to solve the MSE is developed here on the basis of a self-adaptive finite element model (FEM) combined with an iterative method to determine the wave direction angle to the boundary and thus to improve the treatment of the boundary conditions. The numerical resolution of the waves into ideal domains and multidirectional waves through a breakwater gap shows that the numerical model developed here is effective in representing wave absorption at the absorbing boundaries and can be used to simulate multidirectional wave propagation. Finally, the simulated wave distribution in a real harbor shows that the numerical model can be used for engineering practice.     

An efficient technique for finding the eigenvalues of fourth-order Sturm–Liouville problems

In this paper, we will develop a numerical technique for finding the eigenvalues of fourth-order non-singular Sturm–Liouville problems. We used the variational iteration methods as a basis for this technique. Numerical results and conclusions will be presented. Comparison results with others will be presented.     

Numerical study of the stability of breakwater built on a sloped porous seabed under tsunami loading

Strong earthquake induced huge tsunami has occurred for three times in Pacific ocean in recent ten years; for example, the tsunami triggered by the Sumatra earthquake in 2004, Chile earthquake in 2010 and Tohoku earthquake (Japan) in 2011. Tsunami carrying huge energy always would bring high risks to the population living near to coastline. Breakwater is widely used to dissipate the wave energy, and protect coastline and ports. However, they are vulnerable when being attacked by tsunami wave. At present, the interaction mechanism between tsunami, breakwater and its seabed foundation is not fully understood. In this study, the dynamics and stability of a breakwater under the attacking of tsunami wave is investigated by adopting an integrated model PORO-WSSI 2D, in which the VARANS equation for wave motion, and the Biot’s dynamic equation for soil are used. Based on the numerical results, it is found that offshore breakwater interacts intensively with tsunami wave when it overtopping and overflowing over a breakwater. The impact force on the lateral side of breakwater applied by tsunami wave is huge. The shear failure is likely to occur in the seabed foundation of breakwater. The liquefaction is unlikely to occur due to the fact that there is basically no upward seepage force in seabed foundation in the process of tsunami wave passing through the breakwater.     

Singular boundary method based on time‐dependent fundamental solutions for active noise control

Active noise control is an efficient strategy of noise control. A numerical wave shielding model to inhibit wave propagation, which can be considered as an extension of traditional active noise control, is established using the singular boundary method using time‐dependent fundamental solutions in this study. Two empirical formulas to evaluate the origin intensity factors with Dirichlet and Neumann boundary conditions are derived respectively. In comparison with other similar numerical methods, the method can obtain highly accurate results using very few boundary nodes and small CPU time. These meet the major technical requirements of simulation of active noise control. The subsequent numerical experiments show that the proposed model can shield efficiently from the wave propagation for both inner and exterior problems. By applying the newly derived empirical formulas, the CPU time of the singular boundary method is further reduced significantly, which makes the method a competitive new and efficient meshless method. In addition, the singular boundary method makes active noise control in an online manner via time‐dependent fundamental solutions as its basis functions.     

Determination of a control parameter in a one-dimensional parabolic equation using the method of radial basis functions

In this work, the method of radial basis functions is used for finding the solution of an inverse problem with source control parameter. Because a much wider range of physical phenomena are modelled by nonclassical parabolic initial-boundary value problems, theoretical behavior and numerical approximation of these problems have been active areas of research. The radial basis functions (RBF) method is an efficient mesh free technique for the numerical solution of partial differential equations. The main advantage of numerical methods which use radial basis functions over traditional techniques is the meshless property of these methods. In a meshless method, a set of scattered nodes are used instead of meshing the domain of the problem. The results of numerical experiments are presented and some comparisons are made with several well-known finite difference schemes.     

弹性固体和充满两种互不相溶粘性流体的多孔固体之间界面的反射和折射及其衰减波

在充满两种互不相溶粘性流体的多孔固体中,研究弹性波的传播.用3个数性的势函数描述3个纵波的传播,用1个矢性的势函数单独描述横波的传播.根据这些势函数,在不同的组合相中,定义出质点的位移.可以看出,可能存在3个纵波和1个横波.在一个弹性固体半空间与一个充满两种互不相溶粘性流体的多孔固体半空间之间,研究其界面上入射纵波和横波所引起的反射和折射现象.由于孔隙流体中有粘性,折射到多孔介质中的波,朝垂直界面方向偏离.将入射波引起的反射波和折射波的波幅比,作为非奇异的线性代数方程组计算.进一步通过这些波幅比,计算出各个被离散波在入射波能量中所占的份额.通过一个特殊的数值模型,计算出波幅比和能量比系数随入射角的变化.超过SV波的临界入射角,反射波P将不再出现.越过界面的能量守恒原理得到了验证.绘出了图形并对不同孔隙饱和度以及频率的变化,讨论它们对能量分配的影响.     

Simulating wave transmission in the lee of a breakwater in spectral models due to overtopping

Spectral wave models have experienced constant development and vast improvements over the past decades. They are constantly being extended and refined in order to cover the complex wave transformation processes that take place in the coastal zone. Nevertheless, wave transmission due to overtopping has not been treated similarly yet. In this paper, a methodology to include wave generation due to overtopping in spectral wave models is presented. Incorporation of overtopping aims at better simulating the wave disturbance in the lee side of a system of offshore breakwaters and the induced hydrodynamic processes. So far, the waves generated due to wave overtopping were being neglected. The methodology consists of executing sequential simulations at small time step intervals and whenever wave overtopping occurs in a breakwater, waves are generated and transmitted in the lee side of the structure. This is achieved by modifying the boundary condition at the lee of a coastal structure to account for wave generation due to overtopping. Additionally, the transmitted spectrum source function was modified, to capture the observed transfer of energy in the higher frequencies of the spectrum due to the aforementioned overtopping process. The above methodology was implemented in the open source wave model TOMAWAC and verification with experimental measurements was carried out yielding satisfactory results. Inclusion of wave transmission due to overtopping in spectral wave models is considered to be a valuable asset, especially for the simulation of inshore hydrodynamic processes.