导热几何形状反演的变分原理及边界条件的确立

本文采用变域变分原理,建立了导热几何形状反演问题的变分原理,同时获得了该问题所需满足的边界条件和附加条件.该变分原理能将未知形状的几何变量及控制方程结合在一个变分泛函中,使得数学描述简洁、紧凑,且几何变量及控制方程的求解能耦合地进行.介绍了运用该变分原理并结合有限元方法进行数值计算的方法. 关键词： 几何形状反演 变分原理 有限元 导热     

变系数瞬态热传导问题边界元格式的特征正交分解降阶方法

提出了一种基于边界元法求解变系数瞬态热传导问题的特征正交分解(POD)降阶方法,重组并推导出变系数瞬态热传导问题适合降阶的边界元离散积分方程,建立了变系数瞬态热传导问题边界元格式的POD降阶模型,并用常数边界条件下建立的瞬态热传导问题的POD降阶模态,对光滑时变边界条件瞬态热传导问题进行降阶分析.首先,对一个变系数瞬态热传导问题,建立其边界域积分方程,并将域积分转换成边界积分;其次,离散并重组积分方程,获得可用于降阶分析的矩阵形式的时间微分方程组;最后,用POD模态矩阵对该时间微分方程组进行降阶处理,建立降阶模型并对其求解.数值算例验证了本文方法的正确性和有效性.研究表明:1)常数边界条件下建立的低阶POD模态矩阵,能够用来准确预测复杂光滑时变边界条件下的温度场结果;2)低阶模型的建立,解决了边界元法中采用时间差分推进技术求解大型时间微分方程组时求解速度慢、算法稳定性差的问题.     

Combining the complex variable reproducing kernel particle method and the finite element method for solving transient heat conduction problems

In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.     

周期孔洞区域中热力耦合问题的双尺度有限元计算

复合材料的研究中经常遇到具有周期孔洞结构的材料,由于区域的小周期性及剧烈振荡性,用传统的有限元计算方法来计算这些材料对应的问题时需要大量的计算机存储空间及计算时间.对这类材料的热力耦合问题给出了一种新型的高阶双尺度渐近解,得到了对应的均匀化常数、均匀化方程及对应的有限元算法.数值算例表明,周期单胞的局部结构对局部应力与应变有较大的影响.算法对数值模拟这类材料的力学行为是高效和可行的. 关键词： 双尺度方法 热力耦合 周期孔洞区域 有限元方法     

The extended auxiliary equation method for the KdV equation with variable coefficients

This paper applies an extended auxiliary equation method to obtain exact solutions of the KdV equation with variable coefficients. As a result, solitary wave solutions, trigonometric function solutions, rational function solutions, Jacobi elliptic doubly periodic wave solutions, and nonsymmetrical kink solution are obtained. It is shown that the extended auxiliary equation method, with the help of a computer symbolic computation system, is reliable and effective in finding exact solutions of variable coefficient nonlinear evolution equations in mathematical physics.     

Fast multipole accelerated boundary element method for the Helmholtz equation in acoustic scattering problems

We apply the fast multipole method (FMM) accelerated boundary element method (BEM) for the three-dimensional (3D) Helmholtz equation, and as a result, large-scale acoustic scattering problems involving 400000 elements are solved efficiently. This is an extension of the fast multipole BEM for two-dimensional (2D) acoustic problems developed by authors recently. Some new improvements are obtained. In this new technique, the improved Burton-Miller formulation is employed to over-come non-uniqueness difficultie...     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势,     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势, 具有良好的工程应用前景.     

The simulation of acoustic radiation from a moving line source with variable speed

The purpose of this work is to simulate and investigate the sound field generated by a moving line source with finite length and variable speed. Expect for the variation of the acoustic pressure at the specific field point, the distribution of the surface pressure along the surface of the line source was also considered. For achieving this purpose, a numerical method which combines the Time Domain Boundary Element Method (TDBEM) and moving sound source theory was developed in the present work. After comparing the results with the constant and the variable speed case, it showed that the effect of the variable speed not only influenced the variation rate of the frequency modulation, i.e., Doppler effect, but also the time about the maximum acoustic pressure being observed. In addition, the simulation results also presented that the difference as to the amplitude variation of the acoustic pressure still existed between the moving case and the stationary case even if the length of the line source is very long.     

Simulating regular wave propagation over a submerged bar by boundary element method model and Boussinesq equation model

Numerical models based on the boundary element method and Boussinesq equation are used to simulate the wave transform over a submerged bar for regular waves.In the boundary-element-method model the linear element is used,and the integrals are computed by analytical formulas.The Boussinesq-equation model is the well-known FUNWAVE from the University of Delaware.We compare the numerical free surface displacements with the laboratory data on both gentle slope and steep slope,and find that both the models simulate the wave transform well.We further compute the agreement indexes between the numerical result and laboratory data,and the results support that the boundary-element-method model has a stable good performance,which is due to the fact that its governing equation has no restriction on nonlinearity and dispersion as compared with Boussinesq equation.     

开口/封闭薄壳体声辐射和散射的统一边界积分方程解法

推导了在二维和三维空间下开口和封闭薄壳体在任意阻抗边界条件下声辐射和散射的统一边界积分方程.相对于以前的求解方法,该方程求解声辐射和散射问题具有相同的影响矩阵,能够同时求解薄壳体气动和振动噪声的辐射和散射现象,以及分析壳体声阻抗对声波传播的影响.推导的方程可以应用于叶轮机械、管道等噪声和消声器消声性能的预测等方面.在此方程基础上,可以进一步考虑运动边界和运动介质对声辐射和散射的影响. 关键词： 薄壳体 声阻抗 积分方程 边界元方法     

用有限电流元法识别HL-2A 装置的等离子体边界

在电流丝模型对HL-2A装置等离子体边界重建取得满意结果[3]、[5]的基础上，本文用有限电流元模型对HL-2A等离子体边界进行了重建研究。计算和实验结果表明，在通常情况下，有限电流元模型比固定电流丝模型的重建误差稍小，前者误差小于3mm，后者误差小于6mm。当部分有限电流元分布在等离子体边界之外时，用有限电流元模型仍然可以成功重建边界。有限电流元的位置在一定范围内变化时，重建的误差都很小。用普通奔腾4 2.4GHz PC机计算一组等离子体放电边界的时间不超过1ms。有限电流元法能准确而快速地识别等离子体的偏滤器位形，这对于HL-2A装置的实时位形控制是基本和重要的。     

Transparent boundary conditions for the Helmholtz equation in some ramified domains with a fractal boundary

The paper addresses a class of boundary value problems in some self-similar ramified domains, with the Laplace or Helmholtz equations. Much stress is placed on transparent boundary conditions which allow the solutions to be computed in subdomains. A self similar finite element method is proposed and tested. It can be used for numerically computing the spectrum of the Laplace operator with Neumann boundary conditions, as well as the eigenmodes. The eigenmodes are normalized by means of a perturbation method and the spectral decomposition of a compactly supported function is carried out. Finally, a numerical method for the wave equation is addressed.     

Finite element approximation of the Fokker-Planck equation for diffuse optical tomography

In diffuse optical tomography, light transport theory is used to describe photon propagation inside turbid medium. A commonly used simplification for the radiative transport equation is the diffusion approximation due to computational feasibility. However, it is known that the diffusion approximation is not valid close to the sources and boundary and in low-scattering regions. Fokker-Planck equation describes light propagation when scattering is forward-peaked. In this article a numerical solution of the Fokker-Planck equation using finite element method is developed. Approach is validated against Monte Carlo simulation and compared with the diffusion approximation. The results show that the Fokker-Planck equation gives equal or better results than the diffusion approximation on the boundary of a homogeneous medium and in turbid medium containing a low-scattering region when scattering is forward-peaked.     

Optimal configuration for a finite high-temperature source heat engine cycle with the complex heat transfer law

The optimal configuration of a heat engine operating between a finite high-temperature source and an infinite low-temperature reservoir is derived by using finite time thermodynamics based on a complex heat transfer law, including Newtonian heat transfer law, linear phenomenological heat transfer law, radiative heat transfer law, Dulong-Petit heat transfer law, generalized convective heat transfer law and generalized radiative heat transfer law, q ∝ (ΔT ⁿ ). In the engine model the only irreversibility of finite rate heat transfer is considered. The optimal relation between the power output and efficiency of the heat engine is also derived by using an equivalent temperature of the hot reservoir. The obtained results include those obtained in recent literature and can provide some theoretical guidance for the designs of practical engines. Supported by the Program for New Century Excellent Talents in University of China (Grant No. 20041006) and the Foundation for the Author of National Excellent Doctoral Dissertation of China (Grant No. 200136)     

一种适用于任意阶空间差分时域有限差分方法的色散介质通用吸收边界条件算法

基于伸展坐标完全匹配层方程和辅助微分方程 方法, 给出了一种在时域有限差分(FDTD) 计算中适用于常见色散介质的通用边界条件算法. 该算法适用于任意阶的FDTD空间差分, 并且由于所采用的D-H方程独立于计算区域, 所以可以被用于截断任意电介质. 数值试验结果表明, 与卷积完全匹配层 算法相比较, 所提出的吸收边界条件算法不仅通用性强、 计算复杂度低、 计算时间短, 并且吸收效果有明显的提高.     

带源参数的热传导反问题的无网格方法

利用无网格有限点法求带有源参数的一维热传导反问题，推导了相应的离散方程.与其他基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不需要划分网格，用配点法离散求解方程，可以直接施加边界条件，不需要在区域内部求积分，减小了计算量.用有限点法求解热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点.最后通过算例验证了该方法的有效性.     

The projective Riccati equation expansion method and variable separation solutions for the nonlinear physical differential equation in physics

Using the projective Riccati equation expansion (PREE) method, new families of variable separation solutions (including solitary wave solutions, periodic wave solutions and rational function solutions) with arbitrary functions for two nonlinear physical models are obtained. Based on one of the variable separation solutions and by choosing appropriate functions, new types of interactions between the multi-valued and single-valued solitons, such as a peakon-like semi-foldon and a peakon, a compacton-like semi-foldon and a compacton, are investigated.     

辅助方程构造带强迫项变系数组合KdV方程的精确解

在辅助方程法的基础上给出第一种椭圆辅助方程和函数变换相结合的一种方法,并借助符号计算系统Mathematica构造了带强迫项变系数组合KdV方程的类Jacobi椭圆函数精确解以及退化后的类孤子解和三角函数解. 关键词： 辅助方程 函数变换 变系数组合KdV方程 精确解