色散光学模型中色散积分的解析和数值计算

利用解析和数值计算的方法计算了色散光学模型中色散积分的数值,并与Quesada等人的解析计算作了比较,结果显示两者结果基本一致。     

二维柱几何各向异性输运问题的数值计算方法

本文研究了求解二维输运问题的各向异性三角网间断有限元数值方法。详细讨论了求解过程,给出了解的预估公式,证明了数值解的存在唯一性和数值方法的稳定性。给出了各向异性问题的二维数值计算结果及钨——碳临界装置TCCA^[11]的二维数值计算结果。并与国内外S_N数值计算及实验结果进行比较,取得了令人满意的结果。     

基于SIMPLE求解对流-扩散方程的一种新方法

针对SIMPLE系列算法,通过分别求解控制容积界面和节点两个位置的对流-扩散方程来提高数值计算精度,提出了一种流动与传热数值模拟的新方法。采用新的数值方法对顶盖驱动流与正方腔自然对流进行了数值模拟。数值模拟结果表明,在相同网格划分时,新的数值方法相对迎风格式、乘方格式、QUICK格式SIMLE算法的计算精度高;而在计算精度基本相同时,新方法有较高的计算效率。     

多级轴流压气机全工况特性计算

本文使用三维粘性流动计算软件Fine／Numeca,对某十五级轴流压气机进行了内流流场和全工况特性的数值计算尝试。分析了该压气机在设计工况和非设计工况的性能,同时把整机计算结果和前七级叶片的计算结果进行了比较。计算结果表明,当计算的级数较少时,目前的软件和硬件平台可以比较合理地预测压气机的全工况特性;而当计算的级数较多时,准确的数值模拟仍需要更为准确的多级模型和数值方法。     

中子裂变链统计涨落问题的数值计算方法

研究了弱中子源驱动下,裂变系统中子裂变链统计涨落问题的数值计算方法。进行了数值计算建模、数值方法分析、数值计算检验、一类问题概率分布函数的统计涨落特征量的数值计算示范。特例数值检验表明:只要数值解方程组阶数(截断) N足够大,数值解满足归一(守恒) 律、指数增长律,并与精确解析解一致。对于非定常裂变系统中子裂变链统计涨落问题提出了一维等效模型下数值模拟的方法。     

负离子源引出系统束光学的数值模拟

根据在表面-等离子体型负离子源中,负离子主要由转换电极表面产生这一特点,参考正离子源引出系统数值计算程序,建立了负离子源引出系统的数值计算模型和计算程序;对表面-等离子体型桶式负离子源引出系统束光学的性质进行了数值模拟。对计算结果的检验和分析表明,这个计算模型和程序反映了负离子源引出系统束光学的基本特性。     

流动和传热问题的分区并行计算

针对交错网格下的SIMPLE数值算法实施了分区并行计算方法,在小型局域网下实现了流动和传热问题的并行数值计算.对两个经典的流动和传热问题的数值模拟实验表明,所建立的并行计算环境和分区并行算法能够得到正确的和收敛的数值结果.但与串行计算结果相比,并行计算误差明显大于串行计算误差.对并行算法做出的性能分析表明,所给出的并行算法得到了明显的加速效率.随着计算规模的增大,加速比和并行效率提高更显著.     

计算流管声场的三维有限元数值模型

发展了一种三维有限元数值模型和计算方法来对矩形流管声场进行整体的计算.与以往的二维方法相比,此种数值方法不仅全面反映了矩形流管内声波的传播情况,而且提高了网格精度,从而大大扩展了对铺设有声衬的流管的计算领域.结果表明,该数值模型是有效和准确的,与其它方法和文献的计算结果吻合得非常好.同时,在大大增加计算量的同时,也对程序代码进行了优化工作,提高了计算效率.     

对旋轴流水轮机的非定常数值模拟研究

为了研究采用对旋转轮的新型轴流式水轮机的非定常性能,应用三维非定常数值计算方法,对该新型水轮机进行了数值模拟研究.数值计算中,采用了全三维全流道的数值计算方法,基于DES(Detached Eddy Simulation)方法获得了对旋轴流水轮机的非定常流场和能量性能.计算结果显示,该新型水轮机在一个转动周期的不同时刻,其内部流场存在显著的非定常特性.     

近壁面气泡运动特性的数值计算

采用体积加速度模型确定水下爆炸气泡运动的初始条件,基于MSC.DYTRAN有限元软件开发了定义流场初始条件和边界条件的子程序,对气泡在水平刚性壁面附近的运动特性进行了数值模拟,通过对比发现数值计算结果与实验结果具有较好的一致性,证明了初始条件定义、子程序开发和有限元模型建立的正确性和数值计算的准确性。以此为基础,研究了水深、泡心与刚性底面之间的距离对气泡动态特性和射流速度的影响,通过数值计算得到了一些有规律的曲线。计算模型、方法及结果对相关的工程研究和计算具有一定的参考价值。     

时间分数阶Edwards-Wilkinson方程的标度行为研究

分别采用数值模拟和标度分析方法对1+1维时间分数阶Edwards-Wilkinson方程的标度行为进行研究.利用Caputo分数阶导数数值解求得的生长指数与采用直接标度分析方法得到的结果一致.     

On the rate of convergence of the Legendre spectral collocation method for multi-dimensional nonlinear Volterra-Fredholm integral equations

While the approximate solutions of one-dimensional nonlinear Volterra-Fredholm integral equations with smooth kermels are now well understood,no systematic studies of the numerical solutions of their multi-dimensional counterparts exist.In this paper,we provide an efficient numerical approach for the multi-dimensional nonlinear Volterra-Fredholm integral equations based on the multi-variate Legendre-collocation approach.Spectral collocation methods for multi-dimensional nonlinear integral equations are known to cause major difficulties from a convergence analysis point of view.Consequently,rigorous error estimates are provided in the weighted Sobolev space showing the exponential decay of the numerical errors.The existence and uniqueness of the numerical solution are established.Numerical experiments are provided to support the theoretical convergence analysis.The results indicate that our spectral collocation method is more flexible with better accuracy than the existing ones.     

Numerical study of the one-dimensional quantum compass model

The ground state magnetic phase diagram of the one-dimensional quantum compass model (QCM) is studied using the numerical Lanczos method. A detailed numerical analysis of the low energy excitation spectrum is presented. The energy gap and the spin-spin correlation functions are calculated for finite chains. Two kind of the magnetic long-range orders, the Néel and a type of the stripe-antiferromagnet, in the ground state phase diagram are identified. Based on the numerical analysis, the first and second order quantum phase transitions in the ground state phase diagram are identified.     

EXPRESSIONS FOR PREDICTING FUNDAMENTAL NATURAL FREQUENCIES OF NON-CYLINDRICAL HELICAL SPRINGS

Numerical and analytical studies are performed for the free vibration analysis of non-cylindrical (conical, barrel and hyperboloidal types) helical springs. The stiffness matrix method is used in the numerical analysis. A total of 12 degrees of freedom (six displacements and six rotations) is described for an element. The exact element stiffness matrix and the exact concentrated element inertia matrix are used in the formulation. The rotary inertia, the shear and extensional deformation effects are considered in the analysis. Comparison of the numerical results with the reported results obtained numerically and experimentally gives satisfactory values. After verification of the numerical frequencies, the non-dimensional fundamental frequencies of fixed-fixed non-cylindrical helical springs with circular section are expressed in a simple formula with a maximum absolute relative error of 5% using those numerical values for the constant helix pitch angles (5°, 10°, and 15°). These expressions restricted to the fundamental frequencies are also verified with ANSYS results.     

Verification of numerical simulation method for entropy generation of radiation heat transfer in semitransparent medium

The numerical simulation method of radiative entropy generation in participating media presented by Caldas and Semiao [Entropy generation through radiative transfer in participating media: analysis and numerical computation. JQSRT 2005;96:423-37] is extended to analyze the radiative entropy generation in the enclosures filled with semitransparent media. A discrete ordinates method is used to solve radiative transfer equation and radiative entropy generation. Two different examples are employed to verify the numerical simulation method of radiative entropy generation in the enclosure. Numerical results of dimensionless radiative entropy generation of enclosure are identical to that of entire thermodynamics analysis for the enclosure system. This numerical simulation method can be used in the entropy generation analysis of high-temperature systems such as boilers and furnaces, in which radiation is the dominant mode of heat transfer.     

Hydrothermal convection in moderately thin spherical shells

Hydrothermal convection of pore water with a temperature-dependent viscosity within a permeable, internally heated, moderately thin spherical shell is investigated by both a perturbation analysis and a direct numerical simulation. The analysis and simulation are mainly focused on a thin spherical shell in that convective instabilities are characterized by the spherical harmonic degree l=6 with a 13-fold mathematical degeneracy. Four different three-dimensional analytical solutions of convection are derived by removing the degeneracy through the nonlinear effect. A direct numerical simulation of the nonlinear problem is also carried out, showing satisfactory agreement between the analytical solutions and the numerical simulations.     

串级结构Z扫描特性的分析

对非线性介质串级结构的Z扫描特性进行了分析，分别给出了解析解和数值解. 解析分析主要是基于高斯分解方法（GDM）和分布透镜方法的组合；数值分析采用了Crank-Nicholson有限差分方法和快速汉克尔变换方法(FHT). 同时，也给出了这两种数值分析的使用条件和方法. 将解析结果、数值结果和实验结果进行了比较，结果显示它们具有很好的一致性. 关键词： 串级结构 Z扫描 有限差分方法 快速汉克尔变换     

Harmonic Generation of High-Power Microwave in Plasma Filled Waveguide

The harmonic generation of high power microwave in plasma filled waveguide is studied and the analytical theory is presented, and the numerical calculation and analysis of second harmonic are done in this paper. The theoretical and numerical analysis show that high power microwave can generate harmonic in the plasma filled waveguide.     

A numerical study on the enhancement of heat transfer by pressure distributions in acoustic fields

The enhancement of heat transfer in a cavity was investigated in the absence of and in the presence of acoustic streaming induced by ultrasonic waves. The present study provides the experimental and numerical results of heat transfer in the acoustic fields. The enhancement of heat transfer was experimentally investigated in the presence of acoustic streaming and was compared with the profiles of acoustic pressure calculated by the numerical analysis. A coupled finite element-boundary element method (FE-BEM) was applied for a numerical analysis. Experimental and numerical studies clearly show that pressure variations are closely related to the enhancement of heat transfer in the acoustic fields.     

Numerical Method for the Time Fractional Fokker-Planck Equation

In this paper, a new numerical algorithm for solving the time fractional Fokker-Planck equation is proposed. The analysis of local truncation error and the stability of this method are investigated. Theoretical analysis and numerical experiments show that the proposed method has higher order of accuracy for solving the time fractional Fokker-Planck equation.