Meshless analysis of three-dimensional steady-state heat conduction problems

Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson's equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.     

A second-order particle tracking method

An accurate particle tracking method for a finite difference method model is developed using a constant acceleration method. Being assumed constant temporal and spatial gradients, the new method permits temporal-spatial variability of particle velocity. Test results in a solid rotating flow show that the new method has second-order accuracy. The performance of the new method is compared with that of other methods; the first-order Euler forward method, and the second-order Euler predictorcorrector method. The new method is the most efficient method among the three. It is more accurate and efficient than the other two.     

A two-dimensional numerical wave flume based on SA-MPLS method

A spatially adaptive(SA) two-dimensional(2-D) numerical wave flume is presented based on the quadtree mesh system,in which a new multiple particle level set(MPLS) method is proposed to solve the problem of interface tracking,in which common intersection may be traversed by multiple interfaces.By using the adaptive mesh technique and the MPLS method,mesh resolution is updated automatically with time according to flow characteristics in the modeling process with higher resolution around the free surface and the solid boundary and lower resolution in less important area.The model has good performance in saving computer memory and CPU time and is validated by computational examples of small amplitude wave,second-order Stokes wave and cnoidal wave.Computational results also indicate that standing wave and wave overtopping are also reasonably simulated by the model.     

Numerical analysis of shipping water impact on a deck using a particle method

The objective of this study is to simulate three-dimensional behavior of shipping water and to predict the impact pressure on the deck using a particle method. An experiment carried out in a two-dimensional wave tank with a fixed deck model was numerically analyzed in three dimensions. The fluid motion, the surface elevation, and the impact pressure on the deck were compared between the experiment and the calculation and good agreement was obtained.     

粗颗粒两相流计算方法研究

针对具有颗粒粒级跨度大且粒径大特点的固液两相流系统,如深海采矿系统等,开发新型粗颗粒—均质浆体计算模型,该模型采用欧拉—拉格朗日计算法,大粒径颗粒视为固相粗颗粒,其余视为细颗粒与海水组成液相均质浆体。固—液相间滑移速度v_(slip)作为粗细颗粒划分标准,通过联立颗粒雷诺数Re_p,液相对颗粒绕流阻力系数C_d,颗粒粒径d_p等参数,得到滑移速度v_(slip)与颗粒粒径d_p的相关函数,获得粗细颗粒的划分粒径,进而得到两相流控制方程。利用CFD软件,使用该方法对粗颗粒固—液两相流垂直管提升系统进行仿真计算,并将所获得的仿真结果与已发表成果进行对比,发现数据较为吻合,从一定程度上验证了该计算模型的准确性。     

激光-筛分联合法粒度分析结果的表达及粒度端元的沉积动力学意义

河口海岸沉积动力环境复杂,沉积物粗细组分皆有。对于粗细混杂（同时含有砂砾质和泥质成分）的沉积物样品,单一的粒度分析方法不能覆盖整个粒度分布范围,常采用激光-筛分联合法进行粒度分析。由于两种方法的分析原理不同,经常出现细部的激光法结果与粗部的筛分结果部分重合现象,这会影响整个样品的粒度分布计算。石油天然气行业标准《碎屑岩粒度分析方法》（SY/T 5434-2018）对此问题给出了新的解决方案,但是只有文字描述,没有可操作的计算方法。本研究根据该行业标准提出的方案,推算出激光-筛分联合法粒度分析结果的修正计算公式。以浙江舟山群岛潮间带197个粗细混杂样品为例,采用激光-筛分联合法进行了粒度分析,分别采用直接拼接算法及SY/T 5434-2018算法计算了粒度分布,在此基础上对两套算法所得组分含量和粒度参数进行了比较。结果表明,SY/T 5434-2018算法相较于直接拼接算法,扩大了粗颗粒的体积百分比,减少了细颗粒的体积百分比。两种算法得到的组分含量具有显著的线性相关关系,且砾、砂、粉砂、黏土等组分含量的差异从粗到细依次减少。用两种算法的结果分别计算粒度参数,同样具有显著的线性相关关系,且...     

Transparent boundary condition for simulating nonlinear water waves by a particle method

Open boundaries are important when simulating water waves. In this study, a transparent boundary condition at an open boundary was developed for simulating nonlinear water waves propagating to a distant area using the Moving Particle Semi-implicit method. The novelty of this study is that the technique of wave analysis used in the experiment was introduced into the particle simulation to absorb incident waves; the simulation cost was reduced by employing inflow and outflow regions instead of a long dissipation region. Incident waves in front of the boundary were evaluated using Fourier analysis, and the particles on the transparent boundary were forced to move at the velocity of the analytical solution for Stokes waves in order to absorb the incident waves. The analysis was restricted to periodic waves. Wave propagation was simulated for two wave periods using the developed transparent boundary condition. The results showed that this transparent boundary transmitted the incident waves with small reflection and the simulation cost was lower than that for wave damping by a conventional highly viscous region.     

Shared-Memory parallelization of consistent particle method for violent wave impact problems

A shared-memory parallelization is implemented to the recently developed Consistent Particle. Method (CPM) for violent wave impact problems. The advantages of this relatively new particle method lie in four key aspects: (1) accurate computation of Laplacian and gradient operators based on Taylor series expansion, alleviating spurious pressure fluctuation and being able to model two-phase flows characterized by large density difference, (2) a thermodynamics-based compressible solver for modelling compressible air that eliminates the need of determining artificial sound speed, (3) seamless coupling of the compressible air solver and incompressible water solver, and (4) parallelization of the numerical model based on Open Multi-Processing (OpenMP) and a parallel direct sparse solver (Pardiso) to significantly improve computational efficiency. Strong and weak scaling analyses of the parallelized CPM are conducted, showing an efficiency speedup of 100 times or more depending on the size of simulated problem. To demonstrate the accuracy of the developed numerical model, three numerical examples are studied including the benchmark study of wave impact on seawall, and our experimental studies of violent water sloshing under rotational excitations and sloshing impact with entrapped air pocket. CPM is shown to accurately capture highly deformed breaking waves and violent wave impact pressure including pressure oscillation induced by air cushion effect.     

The complex band structure for armchair graphene nanoribbons

Using a tight binding transfer matrix method, we calculate the complex band structure of armchair graphene nanoribbons. The real part of the complex band structure calculated by the transfer matrix method fits well with the bulk band structure calculated by a Hermitian matrix. The complex band structure gives extra information on carrier's decay behaviour. The imaginary loop connects the conduction and valence band, and can profoundly affect the characteristics of nanoscale electronic device made with graphene nanoribbons. In this work, the complex band structure calculation includes not only the first nearest neighbour interaction, but also the effects of edge bond relaxation and the third nearest neighbour interaction. The band gap is classified into three classes. Due to the edge bond relaxation and the third nearest neighbour interaction term, it opens a band gap for N=3M-1. The band gap is almost unchanged for N=3M+1, but decreased for N=3M. The maximum imaginary wave vector length provides additional information about the electrical characteristics of graphene nanoribbons, and is also classified into three classes.     

An improved particle swarm optimization algorithm with harmony strategy for the location of critical slip surface of slopes

The determination of optimal values for three parameters required in the original particle swarm optimization algorithm is very difficult.It is proposed that two new parameters simulating the harmony search strategy can be adopted instead of the three parameters which are required in the original particle swarm optimization algorithm to update the positions of all the particles.The improved particle swarm optimization is used in the location of the critical slip surface of soil slope,and it is found that the improved particle swarm optimization algorithm is insensitive to the two parameters while the original particle swarm optimization algorithm can be sensitive to its three parameters.     

耦合Galerkin简正波解的抛物方程方法水下声传播计算

研究水平变化的海洋环境下声传播的计算方法.把Galerkin方法的简正波解应用于耦合简正波抛物方程,可同时考虑海水和海底声场计算,对水平变化的海洋环境问题的数值计算表明,在包含海水和海底的声场计算中该方法的计算结果都具有较高的精度.     

拉格朗日描述下三层流体系统中界面内波的三阶渐近解

在流体力学中,描述流体运动有Lagrange方法和Euler方法.Euler方法是通过观测通过空间各固定位置点处流体质点的运动行为来描述流体运动规律,而Lagrange方法是跟踪各个流体质点,通过观测它们在时空运动中所走过的路径来描述流体的运动规律.在数学处理上,Euler方法较Lagrange方法简单,但Lagrange方法可以完全描述流体运动的整个流场的所有特性,而Euler方法却无法描述每个流体质点的运动轨迹.本文,我们研究具有刚性边界的三层流体系统中的界面内波,其中上层流体的密度比下层流体的密度大.通过在界面处引入朗格朗日匹配条件并使用微扰法得到了拉格朗日描述下的界面内波的一阶解、二阶解及三阶解,给出了质量输运速度、波频率、平均水平和质点运动轨迹的解.结果表明对于质量输运速度、波频率、平均水平和质点运动轨迹在界面处会有不连续性,但是我们发现在满足一定的三层流体水深比和密度比条件时这种不连续性将会消失.     

复指数方法降噪技术及其试验研究

实测信号中的噪声,以及模型阶次的不确定性给模态参数的准确识别带来困难。以提高模态参数识别精度为目标,提出基于模型定阶和信噪分离的复指数模态参数识别方法。该方法借助奇异值分解技术确定模型阶次,采用结构低秩逼近方法进行信噪分离。在此基础上,利用复指数法进行模态参数识别。分别选取一维的悬臂梁模型和二维的悬挂板模型进行物理模型实验,结果表明:该方法提高模态参数的识别精度,尤其是阻尼比的识别精度,具有较好的工程应用前景。     

计算结构可靠度的RBF神经网络响应面法

对功能函数不能明确表达的问题进行可靠度分析,常采用响应面法。其中二次多项式响应面法应用较为广泛,采用与此方法相同的思路,提出了RBF神经网络响应面法,并通过算例与常用的BP神经网络响应面法进行了对比分析,该方法在学习速度、迭代次数等方面均优于BP神经网络响应面法。该方法用于大型复杂结构的可靠性分析,可相应提高工作效率和解题质量,具有一定实际应用价值。     

渤海海面太阳辐照强度的观测分析与计算方法研究

对渤海海面(1998-09-24～10-07和1999-04-28～05-11)实测太阳辐射日总量同云量的关系和春、秋季的对比关系及其与插值计算出的晴日太阳辐射日总量的关系进行了分析,并同烟台(福山)测站的结果进行了对比分析.通过分析发现,渤海海面太阳辐射日总量春季远大于秋季,其无云条件下的太阳辐射日总量与插值计算出的晴日太阳辐射日总量的差异秋季大于春季.分析表明,海陆温差引起的海面大气垂直输送与水汽物的产生是造成这种差异的主要原因.采用经验公式对渤海海面太阳辐射日总量进行计算,计算结果与实测结果具有较好的一致性.结果表明所用经验公式及推算出的参数适合于春、秋季渤海海面太阳辐照度的计算.