用拟合场量的电流确定等离子体边界

在轴对称的等离子体物理实验装置中,如何用一种简易的方法有效地确定放电时的等离子体边界,是一个令人感兴趣的问题。当然,人们可以用求解自由边界平衡的程序来计算,但这仍避免不了要假设等离子体电流分布;另外,由于求解自由边界问题的复杂性,在没有大型计算机的条件下,要在放电后迅速求出等离子体的边界形状是不现实的。本文根据GBH-1装置上的6根等效电流环,以及真空室外边界上的磁场测量值(或ψ值),和外线圈电流的测     

含有大位移动边界的复杂流场的数值模拟

含有运动弹丸的膛口流场是典型的大位移动边界问题,同时弹丸在出膛过程中,流场的结构也会发生变化,增加了流场的复杂程度.在对该流场的数值模拟中,将它分为两个区:弹丸运动区和普通流场区,它们之间用一个特殊的分区边界联系,同时运用网格局部重构技术处理弹丸运动造成的网格变形问题,并将对称轴定义为网格变形边界,使得弹丸在对称轴上运动过程中不会导致对称轴上的网格体积为负.从计算结果可以看出整个膛口波系结构变化过程和弹丸先加速后减速的过程,从而表明该动网格处理方法是成功的.     

静电磁场不规则区域问题的小波插值Galerkin算法

讨论了用小波插值Galerkin方法(WIGM)求解椭圆型偏微分方程,特别是求解区域不规则时的问题.在归纳出WIGM一般形式的基础上,推导出该方法在Sobolev空间范数下的误差界限为C2-m.提出了一种解决不规则区域中静电磁场场分析问题的数值算法,其中选用对称插值尺度函数为基函数,它的对称性及其与平均插值尺度函数的关系可以在一定程度上降低数值求解的计算量.最后通过计算实例说明该算法的有效性.     

在二维叶栅势流计算中边界元法的应用

一、前言 边界元法是近年来兴起的一种新的基于边界积分方程的数值计算方法.Brebbia将其归之为加权剩余法的一个分支,但该法比有限元和有限差分法更具有解析——数值计算特点.有别于区域计算法,边界元法通过引入一个满足场方程的奇异函数作为权函数,将问题的区域计算转化为边界计算.由于所获得的一组边界积分方程仅联系边界上各个     

抛物型方程一类自由边界问题的微分求积区域分裂法

在微分求积法的基础上,结合区域分裂法的优点提出了一种新的数值计算方法——微分求积区域分裂法.数值试验表明,该方法在求解抛物型方程一类初值带有弱奇性的自由边界问题时十分灵活有效.     

基于区域生长理论的波面解包算法研究

在移相干涉技术的相位计算中,周期函数arctan的相位值分布在[-π,π]之间,相位值分布的跃变形成了被压包后的波面。理想无噪声和边界形状规则的压包波面可以直接通过条纹扫描法进行解包。例如当被测物边界形状不规则和存在噪声(灰尘脏点)时,会导致干涉图上出现无效区域,必须采用不依赖于边界形状,能自动绕过无效区域的算法。采用区域生长法对这类干涉图进行了解包。通过编程计算表明,这种解包算法能解决光学检测中的几乎所有的干涉图处理问题。     

搅拌器内部流动的无网格法三维数值模拟

本文将移动粒子半隐式法(MPS)的基本算法由二维扩展至三维。将圆柱坐标系引入到初场粒子的布置中,避免了在笛卡儿坐标系下处理不规则形状(如斜边或曲边)问题时粒子初场布置困难和精确度较低的问题,改善了对计算边界条件表达的精确性。引入移动边界模型,对直叶片搅拌器的内部流动进行了三维数值模拟。还提出了一种新的初始粒子布置简易方法,明显简化粒子初始布置时的复杂程度,提高了对三维复杂几何形状问题的可操作性。     

蓝宝石生长炉热场的全局/局部耦合算法研究

本文对工业用泡生法制备蓝宝石单晶炉热场建立了二维稳态全局及局部数值模型.基于全局模拟结果,比较分析了不同热内边界条件的选取对局部核心区域热场及流场分布的影响.结果表明,基于全局/局部耦合算法的计算方案是可行的,当辐射-对流边界选取第二类边界、辐射-导热边界选取第一类边界、导热-导热边界选取第一类或第二类边界时,局部模拟结果的准确性最高,且局部网格划分与全局相应区域网格划分不一致时,对模拟结果几乎无影响.     

一种基于多矩的有限体积浸入边界法

基于多矩VSIAM3格式及浸入边界法,提出一套在复杂计算区域内求解不可压缩流动的数值格式.不可压N-S方程使用VSIAM3格式进行离散,引入浸入边界法处理复杂、移动边界,使用虚拟网格方法计算动量方程修正项,同时还考虑了对连续方程的修正.使用标准算例对数值模式进行验证.     

无网格方法中的背景积分方案及单颗粒下降过程的数值模拟

在用无网格方法对有边界移动的气固两相流动进行数值模拟时发现,如果采用背景积分网格对求解域进行积分,将会产生压力场的数值振荡,使计算结果失真.详细分析了产生这种数值振荡的原因,提出了解决方案,并计算了一个颗粒在管道内的下降过程.计算结果表明,提出的解决方案能有效地降低压力场的数值振荡,从而使无网格方法能应用于气固两相流动的直接数值模拟.     

抛物方程基于自然边界归化的耦合法

将冯康和余德浩提出的自然边界归化方法^[1~4]应用于求解抛物方程初边值外区域问题,提出一种自然边界元与有限元耦合算法。先将控制方程对时间进行离散化,得到关于时间步长的离散化格式,给出圆外域上的自然积分方程,基于此研究抛物方程无界区域问题的自然边界元与有限元耦合法,最后给出相应的数值例子。     

An approximation for the boundary optimal control problem of a heat equation defined in a variable domain

In this paper, we consider a numerical approximation for the boundary optimal control problem with the control constraint governed by a heat equation defined in a variable domain. For this variable domain problem, the boundary of the domain is moving and the shape of theboundary is defined by a known time-dependent function. By making use of the Galerkin finite element method, we first project the original optimal control problem into a semi-discrete optimal control problem governed by a system of ordinary differential equations. Then, based on the aforementioned semi-discrete problem, we apply the control parameterization method to obtain an optimal parameter selection problem governed by a lumped parameter system, which can be solved as a nonlinear optimization problem by a Sequential Quadratic Programming （SQP） algorithm. The numerical simulation is given to illustrate the effectiveness of our numerical approximation for the variable domain problem with the finite element method and the control parameterization method.     

FGMs稳态热传导分析的重心Lagrange插值配点法

提出一种求解二维功能梯度材料（FGMs）稳态热传导问题的重心Lagrange插值配点法.基于Chebyshev节点构造二维重心Lagrange插值函数及其偏导数,然后基于配点法将其直接代入FGMs热传导问题的控制方程和边界条件,得到系统离散方程.重心Lagrange插值配点法是一种真正的无网格方法,很好地融合了重心Lagrange插值和配点格式的优势,具有高效、稳定、高精度和易于数值实现的优点.采用重心Lagrange插值配点法分别对指数型、二次型和三角型FGMs热传导问题进行数值模拟.结果表明：该方法具有较高的计算效率和计算精度,对材料梯度参数的变化不敏感.可以进一步拓展到FGMs瞬态问题和FGMs的热力耦合分析.     

Element-free Galerkin method for a kind of KdV equation

The present paper deals with the numerical solution of the third-order nonlinear KdV equation using the element-free Galerkin (EFG) method which is based on the moving least-squares approximation. A variational method is used to obtain discrete equations, and the essential boundary conditions are enforced by the penalty method. Compared with numerical methods based on mesh, the EFG method for KdV equations needs only scattered nodes instead of meshing the domain of the problem. It does not require any element connectivity and does not suffer much degradation in accuracy when nodal arrangements are very irregular. The effectiveness of the EFG method for the KdV equation is investigated by two numerical examples in this paper.     

重心Lagrange插值配点法求解二维双曲电报方程

提出一种求解二维双曲电报方程的高精度重心Lagrange插值配点法.采用重心Lagrange插值构造包含时间和空间变量的近似函数.在给定Chebyshev-Gauss-Lobatto节点上,将多变量重心Lagrange插值近似函数代入双曲电报方程及其定解条件,得到离散代数方程组.包含狄里克雷和诺依曼边界条件的数值算例表明,本文方法程序实现方便并具有高精度,可应用于求解高维问题.     

A New Numerical Method for Solving the Acoustic Radiation Problem

A numerical method of solving the problem of acoustic wave radiation in the presence of a rigid scatterer is described. It combines the finite element method and the boundary algebraic equation one. In the proposed method, the exterior domain around the scatterer is discretized, so that there appear an infinite domain with regular discretization and a relatively small layer with irregular mesh. For the infinite regular mesh, the boundary algebraic equation method is used with spurious resonance suppression according to Burton and Miller. In the thin layer with irregular mesh, the finite element method is used. The proposed method is characterized by simple implementation, fair accuracy, and absence of spurious resonances.     

An efficient differential quadrature methodology for free vibration analysis of arbitrary straight-sided quadrilateral thin plates

In this paper, a new differential quadrature (DQ) methodology is employed to study free vibration of irregular quadrilateral straight-sided thin plates. A four-nodded super element is used to map the irregular physical domain into a square domain in the computational domain. Second order transformation schemes with relative ease and less computation are employed to transform the fourth order governing equation of thin plates between the two domains. The only degree of freedom within the domain is the displacement, whereas along the boundaries, the displacement as well as the second order derivative of the displacement with respect to associated normal co-ordinate variable in computational domain are the two degrees of freedom. Implementing the method, the formulation for the DQ method for the free vibration analysis of plates of straight-sided shapes was presented together with the implementation procedure for the different boundary conditions. To demonstrate the accuracy, convergency and stability of the new methodology, detail studies are made on isotropic plates at acute angles with different geometries, boundary and loading conditions including DQ free-edge boundary condition implementations. Accurate results even with fewer degrees of freedom than for those of comparable numerical algorithms were achieved.     

基于平面问题的位移压力混合配点法

引入压力变量,将弹性力学控制方程表达为位移和压力的耦合偏微分方程组,采用重心插值近似未知量,利用重心插值微分矩阵得到平面问题控制方程的矩阵形式离散表达式.采用重心插值离散位移和应力边界条件,采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法求解过约束方程组,得到平面问题位移数值解.数值算例验证了所提方法的有效性和计算精度.     

激光辐照下复合材料热响应的轴对称数值分析

为了模拟激光辐照下碳纤维增强复合材料的瞬态热响应,建立了轴对称计算模型。模型考虑了激光辐照过程中基体热分解、质量迁移、比热容和热导率等物理量的变化情况。采用有限元方法求解控制方程,边界条件包含了激光辐照、对流换热以及辐射换热。在此基础上编写了计算程序,预测了激光辐照下碳纤维增强复合材料的瞬态温度场和基体热分解状况。为了校核模型,开展了激光辐照碳纤维复合材料试验。计算结果与试验数据比较表明,模型预测的复合材料温度-时间曲线与试验结果较好吻合,在较低功率密度激光辐照下复合材料热响应以基体热分解为主,与试验烧蚀形貌观察结果一致。     

Element-free Galerkin （EFG） method for a kind of two-dimensional linear hyperbolic equation

The present paper deals with the numerical solution of a two-dimensional linear hyperbolic equation by using the element-free Galerkin (EFG) method which is based on the moving least-square approximation for the test and trial functions. A variational method is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Compared with numerical methods based on mesh, the EFG method for hyperbolic problems needs only the scattered nodes instead of meshing the domain of the problem. It neither requires any element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. The effectiveness of the EFG method for two-dimensional hyperbolic problems is investigated by two numerical examples in this paper.