多重网格方法求解两类Helmholtz方程

详细给出了多重网格方法的实现过程,借助正定Helmholtz方程及不定Helmholtz方程的求解来探讨多重网格方法的特性。对多重网格V环、W环以及F环三种不同迭代格式的收敛效果进行了对比。通过正定Helmholtz方程的求解,发现多重网格的确有很高的计算效率。对于不定Helmholtz方程,随着波数的增加,利用多重网格方法得到结果不收敛,原因出在细网格光滑和粗网格矫正过程。如何针对此问题对多重网格进行有效改进还有待进一步研究。     

用无网格方法求解声波散射问题

在均匀介质中,对软表面障碍、时间调和声波散射问题归结为Helmholtz方程的Dirichlet外问题.应用无网格方法求解Helmholtz方程的Dirichlet外问题,并给出了一个数值例子,与Nystrom方法进行了对比,表明该方法是较精确的.     

一种针对大波数Helmholtz方程的高性能并行预条件迭代求解算法

针对传统串行迭代法求解大波数Helmholtz方程存在效率低下且受限于单机内存的问题,提出了一种基于消息传递接口(Message Passing Interface,MPI) 的并行预条件迭代法。该算法利用复移位拉普拉斯算子对Helmholtz方程进行预条件处理,联合稳定双共轭梯度法和基于矩阵的多重网格法来求解预条件方程离散后的大规模线性系统,在Linux集群系统上基于 MPI环境实现了求解算法的并行计算,重点解决了多重网格的并行划分、信息传递和多重网格组件的构建问题。数值实验表明,对于大波数问题,提出的算法具有良好的并行加速比,相较于串行算法极大地提高了计算效率。     

Poisson-Boltzmann方程的并行有限元求解及网格自适应生成

说明如何利用并行自适应有限元软件平台PHG 求解生物分子溶液体系的非线性Poisson-Boltzmann方程,并介绍一种解决这类问题的方法,它将网格生成与自适应计算过程结合在一起,可自动产生合适的网格,避免复杂的曲面网格生成步骤.之前的网格生成工作有：（1） TMSmesh生成高斯曲面的三角网格; （2） TransforMesh删除自相交的三角网格; （3） ISO2Mesh提高表面网格质量3个步骤.而基于PHG的自适应加密模块可以在逐次调整网格的同时保持动态负载平衡,高效地得到计算网格用于近似求解非线性Poisson-Boltzmann方程.计算了小球模型和AChE系统,分别从误差指示子下降阶和溶剂化能收敛的角度验证了方法的有效性,并且还将网格生成算法成功地应用于gA离子通道.     

稠密网格上交互指定局部区域的新方法

在造型、绘制以及动画中,在网格上指定局部区域都是一个重要的操作。传统的3D动画软件采用的交互方法在选取稠密网格上的特征局部区域时需要繁琐的用户操作,而现有的其他交互方法又难以做到实时响应。该文提出了一个在网格上指定局部区域的新方法,该方法基于网格上的近似最短路径算法。该算法不对网格进行全局加细,而是利用预处理阶段求出的信息对网格进行局部加细,提高了算法的效率。而且通过网格上近似等距线的抽取,还能指定网格上的拓扑同构于环形的窄带区域。     

结构变形对整流罩分离轨迹的影响

为研究弹性体在稠密大气中的分离问题,基于非结构网格,采用运动网格与局部网格重构相结合的方法求解大位移相对运动的流场,并耦合6自由度刚体运动方程得到整流罩的运动.非定常流动方程使用格心有限体积法进行空间离散,并运用LU-SGS进行求解.应用标准算例验证该方法的准确性,并用于某整流罩飞行轨迹的计算.结果表明结构变形可能会使...     

含曲率的水平集方程在非结构四边形网格上的数值离散方法

在非结构四边形网格上,含曲率的水平集方程采用伽辽金等参有限元方法空间离散,时间离散采用半隐格式．离散形成的线性方程组的系数矩阵是对称的稀疏矩阵,采用共轭梯度法求解．数值算例表明,在笛卡儿网格和随机网格上,含曲率的水平集方程离散格式可达到近似二阶精度．重新初始化方程的离散格式精度可达到近似一阶精度．给出了非结构四边形网格上不光滑界面以曲率收缩的运动过程．在不采用重新初始化的情况下,收缩过程未出现不稳定现象．     

一种改进的基于深度神经网络的偏微分方程求解方法

偏微分方程求解是计算流体力学等科学与工程领域中数值分析的计算核心。由于物理的多尺度特性和对离散网格质量的敏感性,传统的数值求解方法通常包含复杂的人机交互和昂贵的网格剖分开销,限制了其在许多实时模拟和优化设计问题上的应用效率。提出了一种改进的基于深度神经网络的偏微分方程求解方法TaylorPINN。该方法利用深度神经网络的万能逼近定理和泰勒公式的函数拟合能力,实现了无网格的数值求解过程。在Helmholtz、Klein-Gordon和Navier-Stokes方程上的数值实验结果表明,TaylorPINN能够很好地拟合计算域内时空点坐标与待求函数值之间的映射关系,并提供了准确的数值预测结果。与常用的基于物理信息神经网络方法相比,对于不同的数值问题,TaylorPINN将预测精度提升了3~20倍。     

Helmholtz方程的楔形基区域分解法

基于楔形基函数和无网格配点法,提出了一种求解Helmholtz型方程区域分解法。该方法克服了在求解大规模问题时用一般的全域配点法所带来的配置矩阵为非对称满阵,且高度病态的问题。通过数值结果表明,该算法在求解Helmholtz型方程降低系数矩阵条件数的同时,也能够降低误差,并达到满意的收敛效果。     

逼近误差有界的相容性高阶网格生成

文中提出了一种构造逼近误差有界的高质量相容性高阶网格的方法。给定两个定向的、拓扑同构的三角形网格和一组稀疏的对应点，此方法包含两个步骤：(1)生成满足误差有界的相容性高阶网格；(2)在确保逼近误差总是有界的前提下，降低网格的几何复杂度，并在该过程中通过优化控制顶点来降低相容性网格之间的扭曲以及与原始网格之间的几何近似误差。第一步先生成满足误差有界的相容性线性网格，然后升阶为高阶网格。第二步通过迭代地执行基于边长的重新网格化和增加相容性目标边长场，有效地降低了网格几何复杂度。从切空间的角度，推导出了3DBézier三角形之间映射的雅可比矩阵，从而可以有效地优化扭曲能量。通过对扭曲能量和几何近似误差能量的优化，有效地降低了相容性网格之间的扭曲以及相容性网格与原始网格之间的几何近似误差。通过大量实验，证明了此方法对于构造误差有界的高质量相容性高阶网格的有效性和实用性。     

Construction of the matched multiple knot B-spline wavelets on a bounded interval

In order to solve a boundary value problem using Galerkin's method, the selection of basis functions plays a crucial rule. When the solution of a boundary value problem is not enough smooth or the domain is irregular, multiple knot B-spline wavelets (MKBSWs) with locally compact support are appropriate basis functions. However, to have globally continuous basis functions, a matching across the subdomain interfaces is required. In other words, MKBSWs that are non-zero in the interelement boundaries should be matched. In this paper, we present the primal and dual matched multiple knot B-spline scaling and wavelet functions whose main properties of smoothness and biorthogonality are kept.     

B-Spline Curve Approximation by Utilizing Big Bang-Big Crunch Method

The location of knot points and estimation of the number of knots are undoubtedly known as one of the most difficult problems in B-Spline curve approximation. In the literature, different researchers have been seen to use more than one optimization algorithm in order to solve this problem. In this paper, Big Bang-Big Crunch method (BB-BC) which is one of the evolutionary based optimization algorithms was introduced and then the approximation of B-Spline curve knots was conducted by this method. The technique of reverse engineering was implemented for the curve knot approximation. The detection of knot locations and the number of knots were randomly selected in the curve approximation which was performed by using BB-BC method. The experimental results were carried out by utilizing seven different test functions for the curve approximation. The performance of BB-BC algorithm was examined on these functions and their results were compared with the earlier studies performed by the researchers. In comparison with the other studies, it was observed that though the number of the knot in BB-BC algorithm was high, this algorithm approximated the B-Spline curves at the rate of minor error.     

基于MATLAB的双旋翼控制系统的优化设计

对于非线性复杂控制系统，传统的求解优化设计问题的方法难以获得良好的控制效果。为解决这一问题，使非线性复杂控制系统中的PID控制器的参数整定更加简单方便，文中分析了MATLAB优化工具箱求解优化问题与传统的求解方法相比的主要优点，提出了一种运用MATLAB优化工具箱中提供的函数，对双旋翼非线性复杂控制系统进行优化设计的方法。仿真结果证明该方法是简单可行的，而且可使设计者将精力完全集中于优化问题的本身，而不是过多地纠缠在具体的算法实现上，计算简单方便，结果可靠。     

基于Matlab的支持向量机工具箱

介绍了基于MATLAB的支持向量机工具箱,详细说明了工具箱中用于支持向量分类和支持向量回归的函数.并通过两个具体的实例来说明利用SVM工具箱进行分类和回归方面的方法.     

动态区间系统的最优保性能控制--LMI方法

针对动态区间系统和一个给定的二次型性能指标,研究了其保性能控制问题,基于线性矩阵不等式(LMI)提出了最优保性能控制器设计方法,并将相关结果推广到参数不确定系统.利用功能强大的LMI工具,求解非常方便.所给实例表明,该方法用于设计动态区间系统与秩-1型参数不确定系统的最优保性能控制器,非常有效.     

Shape optimization of continua using NURBS as basis functions

The present paper introduces a numerical solution to shape optimization problems of domains in which boundary value problems of partial differential equations are defined. In the present paper, the finite element method using NURBS as basis functions in the Galerkin method is applied to solve the boundary value problems and to solve a reshaping problem generated by the H1 gradient method for shape optimization, which has been developed as a general solution to shape optimization problems. Numerical examples of linear elastic continua illustrate that this solution works as well as using the conventional finite element method.     

离散区间系统的二次稳定性及其稳定裕度分析

研究了动态离散区间系统的鲁棒稳定性问题,提出了系统二次稳定的充分与必要条件,以及相应的稳定裕度的计算方法,并推广到离散参数不确定系统.所有结论以线性矩阵不等式(LMI)的形式给出.利用功能强大的LMI工具,求解非常方便.所给实例表明,该方法用于确定动态离散区间系统的鲁棒稳定性及其稳定裕度,非常有效.     

NURBS-enhanced line integration boundary element method for 2D elasticity problems with body forces

A NURBS-enhanced boundary element method for 2D elasticity problems with body forces is proposed in this paper. The non-uniform rational B-spline (NURBS) basis functions are applied to construct the geometry and the model can be reproduced exactly at all stages since the refinement will not change the shape of the boundary. Both open curves and closed curves are considered. The fields are approximated by the traditional Lagrangian basis functions in parameter space, rather than by the same NURBS basis functions for geometry approximation. The parametric boundary elements and collocation nodes are defined from the knot vector of the curve and the refinement of the NURBS curve is easy. Boundary conditions can be imposed directly since the Lagrangian basis functions have the property of delta function. In addition, most methods for the treatment of singular integrals in traditional boundary element method can be applied in the proposed method. To overcome the difficulty for evaluation of the domain integrals in problems with body forces, a line integration method is further applied in this paper to compute the domain integrals without additional volume discretizations. Numerical examples have shown the accuracy of the proposed method.     

Numerical investigation on convergence of boundary knot method in the analysis of homogeneous Helmholtz, modified Helmholtz, and convection-diffusion problems

This paper concerns a numerical study of convergence properties of the boundary knot method (BKM) applied to the solution of 2D and 3D homogeneous Helmholtz, modified Helmholtz, and convection-diffusion problems. The BKM is a new boundary-type, meshfree radial function basis collocation technique. The method differentiates from the method of fundamental solutions (MFS) in that it does not need the controversial artificial boundary outside physical domain due to the use of non-singular general solutions instead of the singular fundamental solutions. The BKM is also generally applicable to a variety of inhomogeneous problems in conjunction with the dual reciprocity method (DRM). Therefore, when applied to inhomogeneous problems, the error of the DRM confounds the BKM accuracy in approximation of homogeneous solution, while the latter essentially distinguishes the BKM, MFS, and boundary element method. In order to avoid the interference of the DRM, this study focuses on the investigation of the convergence property of the BKM for homogeneous problems. The given numerical experiments reveal rapid convergence, high accuracy and efficiency, mathematical simplicity of the BKM.