求解隐式差分方程的并行算法

本文研究了求解隐式差分方程的并行算法,其基本想法是把差分方程组化为若干个子方程组来求解,文中给出了直接法和迭代法,证明了迭代法的收敛性。为了说明分段隐式迭代方法的有效性,文中针对扩散方程的具体例子给出了数值试验计算结果。     

等离子体中Fokker-Planck方程的有限元模拟

利用有限元方法设计了一套相对简单明了的求解Fokker-Planck方程的方案.这个方案不必严格限制计算格点的步长和时间步长,就可以确保分布函数的非负性和粒子数的守恒.通过一维程序模拟,进一步证实了这个方案的可靠性.对于多维问题的分析和一维问题完全一样,所以非常容易将其推广到多维问题.     

粘性流动的局部隐式有限元解法

粘性流动的局部隐式有限元解法依凤鸣，吴遇春（哈尔滨工业大学动力系哈尔滨１５Ｏ００１）关键词粘性流动，局部隐式有限元法，Ｎ．Ｓ．方程１引言有限元法具有对复杂几何边界适应性强，对流场内变量梯度大的区域计算精度高以及易于编制大型通用程序等特点，已经成了叶轮...     

热辐射输运问题的隐式蒙特卡罗方法求解

热辐射与物质相互作用及辐射光子在物质中的传输是惯性约束聚变研究中的重要课题. 介绍了基于隐式蒙特卡罗方法的辐射输运方程,在该方程的积分-微分形式基础上,推导了利于蒙特卡罗方法模拟的等价的积分输运方程;基于积分方程设计数值模拟流程,编写三维蒙特卡罗数值模拟程序;针对热辐射输运典型问题及benchmark问题开展了数值实验,计算结果验证了方法的适应性及程序的正确性. 关键词： 热辐射 惯性约束聚变 输运方程 隐式蒙特卡罗     

应用两步隐式算法求解Navier-Stokes方程

本文提出一种隐式分步算法求解流体的运动方程。研究证明,SIMPLEC算法是这种算法的特例。通过两组层流算例的计算表明,改变加权隐式系数α可以明显加快计算的收敛速度。     

一种求解可压缩Navier-Stokes方程的隐式迭代时间推进方法

针对有壁面边界的可压缩流动问题,提出与基于非等距网格的高精度紧致型差分格式相结合的简化隐式迭代时间推进法,建立求解可压缩Navier-Stokes方程的直接数值模拟方法,提高了计算效率.应用该方法,直接数值模拟两种有壁面边界的二维可压缩流动问题,即可压缩平板边界层流动和可压缩槽道流动.     

基于大涡模拟的移动粒子半隐式算法

移动粒子半隐式方法(MPS)是一种粒子方法,多用于模拟带有自由表面的不可压缩流动。工程实际中的自由表面流动往往是复杂的湍流流动,本文借鉴网格类方法的亚格子应力模型发展了基于Smagorinsky模型的亚粒子应力模型,并将其耦合到MPS方法中,实现了基于大涡模拟的MPS方法并用于研究自由表面湍流问题。为了提高计算的准确性和稳定性,SPS模型中出现的一阶导数项采用最小二乘法拟合得到,SPS项采用显式算法进行计算。使用这一算法模拟了溃坝问题,结果表明,采用亚粒子应力模型的模拟结果与实验的吻合程度明显提高。     

解二维时间相关光子扩散方程的交替方向隐式法

研究了应用交替方向隐式法求解非均匀生物组织中的二维时间相关光子扩散方程,并同解析解以及Monte Carlo模拟的结果作了比较,数值结果表明,该方法具有较好的稳定性和较高的精度.     

一种新的LU型隐式格式及其应用

发展了一种新的ＬＵ型隐式格式。新格式除包含ＬＵ－ＳＧＳ格式全部特色外，尽管采用了相似变换而不是近似处理来精确地构造通量矢量的迎风Ｊａｃｏｂｉａｎ矩阵，仍可避免块对角矩阵求逆．且左端项仍可保持矢量化处理。新格式显示出更快的收敛性和更高的稳定性且没有多余的数值耗散和自由参数，尤其适用于求解三维非定常流动问题。     

隐式格式求解拟压缩性非定常不可压Navier-Stokes方程

采用Rogers发展的双时间步拟压缩方法,数值求解不可压非定常问题.数值通量分别采用三阶精度Roe格式和二阶精度Harten-Yee的TVD格式离散.为了加快收敛,提高求解效率,试验了几种隐式格式(ADI-LU,LGS,LU-SGS).针对经典的低雷诺数(Re=200)圆柱绕流问题,比较了不同隐式方法的计算结果和求解效率,以及两种数值离散格式计算结果的异同.最后采用Roe格式数值求解了两种典型的低速非定常流动问题:绕转动圆柱(ω=1)低雷诺数流动;NACA0015翼型等速拉起数值模拟.     

ISSDE: A Monte Carlo implicit simulation code based on Stratonovich SDE approach of Coulomb collision

A Monte Carlo implicit simulation program, Implicit Stratonovich Stochastic Differential Equations(ISSDE), is developed for solving stochastic differential equations(SDEs) that describe plasmas with Coulomb collision. The basic idea of the program is the stochastic equivalence between the Fokker–Planck equation and the Stratonovich SDEs. The splitting method is used to increase the numerical stability of the algorithm for dynamics of charged particles with Coulomb collision. The cases of Lorentzian plasma, Maxwellian plasma and arbitrary distribution function of background plasma have been considered. The adoption of the implicit midpoint method guarantees exactly the energy conservation for the diffusion term and thus improves the numerical stability compared with conventional Runge–Kutta methods. ISSDE is built with C++ and has standard interfaces and extensible modules. The slowing down processes of electron beams in unmagnetized plasma and relaxation process in magnetized plasma are studied using the ISSDE, which shows its correctness and reliability.     

Alternating Direction Implicit Orthogonal Spline Collocation on Non-Rectangular Regions

The alternating direction implicit (ADI) method is a highly efficient technique for solving multi-dimensional time dependent initial-boundary value problems on rectangles. When the ADI technique is coupled with orthogonal spline collocation (OSC) for discretization in space, we not only obtain the global solution efficiently, but the discretization error with respect to space variables can be of an arbitrarily high order. In [2], we used a Crank Nicolson ADI OSC method for solving general nonlinear parabolic problems with Robin's boundary conditions on rectangular polygons and demonstrated numerically the accuracy in various norms. A natural question that arises is: Does this method have an extension to non-rectangular regions? In this paper, we present a simple idea of how the ADI OSC technique can be extended to some such regions. Our approach depends on the transfer of Dirichlet boundary conditions in the solution of a two-point boundary value problem (TPBVP). We illustrate our idea for the solution of the heat equation on the unit disc using piecewise Hermite cubics.     

一类非线性色散Boussinesq方程的隐式孤立波解

用动力系统分岔方法研究了一类非线性色散Boussinesq方程.在不同的参数条件下,给出了该方程具有隐函数形式的孤立波解的解析表达式.数值模拟进一步验证了所得结果的正确性.     

Density Functional Approach Based on Numerically Obtained Bridge Functional

The Ornstein Zernike equation is solved with the Rogers Young approximation for bulk hard sphere fluidand Lennard-Jones fluid for several state points. Then the resulted bulk fluid radial distribution function combinedwith the test particle method is employed to determine numerically the function relationship of bridge functional as afunction of indirect correlation function. It is found that all of the calculated points from different phase space statepoints for a same type of fluid collapse onto a same smooth curve. Then the numerically obtained curve is used tosubstitute the analytic expression of the bridge functional as a function of indirect correlation function required in themethodology [J. Chem. Phys. 112 (2000) 8079] to deterrnine the density distribution of non-uniform hard spherefluid and Lennard Jones fluid. The good agreement of theoretical predictions with the computer simulation data isobtained. The present numerical procedure incorporates the knowledge of bulk fluid radial distribution function intothe constructing of the density functional approximation and makes the original methodology more accurate and moreflexible for various interaction potential fluid.     

碲基掺铒光纤放大器增益饱和特性的数值模拟

依据均匀加宽四能级结构速率方程组和光功率传输方程组,数值模拟了1520～1620nm波段范围内碲基掺铒光纤放大器(EDTFA)增益谱特性,以及EDTFA增益谱特性与输入信号光功率的关系.模拟结果显示:随着输入信号光功率的增大,EDTFA增益谱由于饱和效应而呈现平坦化特征.模拟结果与报道的实验测量结果达到了很好的一致.     

基于Z型石墨纳米条带的量子点：第一性原理研究

本文采用全自洽的格林函数和电子密度泛函理论方法对Z型石墨纳米条带的电子结构和输运性质进行了理论研究。理论计算结果表明Z型石墨纳米条带可以用来设计量子点。由于分子结的特殊几何结构，一些分立电子态可以被囚禁在中间纳米条带内，表现出明显量子限域效应。我们还发现这些被囚禁电子态的数目和它们的空间分布特征是可以调制的，然而这种基于Z型石墨纳米条带的量子点总是可以实现，而不受掺杂、边界化学修饰以及分子结长度的影响     

第一性原理下铜锰共掺铌酸锂晶体的电子结构和吸收光谱

采用基于密度泛函理论的第一性原理研究了Cu、Mn单掺及共掺LiNbO3晶体的电子结构和光学性质.结果显示,Cu、Mn掺杂LiNbO_3晶体禁带中的杂质能级分别由Cu 3d轨道、Mn 3d轨道贡献;各掺杂体系的带隙均较纯LiNbO_3晶体变窄.共掺晶体中Cu离子形成了较单掺时更浅的能级中心,并在2.87eV处有较强的吸收峰;Mn离子在1.73eV附近的吸收较单掺时减弱且中心略有偏移,在2.24eV处的非光折变峰与Mn~(3+)相关,这对吸收峰的变化被认为与Cu、Mn间电子转移相联系.相对Cu、Fe共掺LiNbO_3晶体,Cu、Mn共掺LiNbO_3晶体可以通过适当提高Cu离子浓度,来改善存储参量中的动态范围和记录灵敏度.由于同一深能级掺杂离子伴以不同浅能级掺离子将呈现出不同的吸收特征并影响存储性能,在共掺离子的配搭选择时对各待选配搭的模拟计算非常必要.     

Computational Simulations of Vocal Fold Vibration: Bernoulli Versus Navier–Stokes

The use of the mechanical energy (ME) equation for fluid flow, an extension of the Bernoulli equation, to predict the aerodynamic loading on a two-dimensional finite element vocal fold model is examined. Three steady, one-dimensional ME flow models, incorporating different methods of flow separation point prediction, were compared. For two models, determination of the flow separation point was based on fixed ratios of the glottal area at separation to the minimum glottal area; for the third model, the separation point determination was based on fluid mechanics boundary layer theory. Results of flow rate, separation point, and intraglottal pressure distribution were compared with those of an unsteady, two-dimensional, finite element Navier-Stokes model. Cases were considered with a rigid glottal profile as well as with a vibrating vocal fold. For small glottal widths, the three ME flow models yielded good predictions of flow rate and intraglottal pressure distribution, but poor predictions of separation location. For larger orifice widths, the ME models were poor predictors of flow rate and intraglottal pressure, but they satisfactorily predicted separation location. For the vibrating vocal fold case, all models resulted in similar predictions of mean intraglottal pressure, maximum orifice area, and vibration frequency, but vastly different predictions of separation location and maximum flow rate.     

Formation of Superheavy Elements: Study Based on Dynamical Approach

Using multi-dimensional Langevin equations for the probability distribution of the distance between the surfaces of two approaching nuclei, we have studied the formation of superheavy elements via calculation of evaporation and fission cross sections of these elements. Evaporation residue cross sections have been calculated for the 1n, 2n, 3n, 4n, and 5n evaporation channels using one and four dimensional Langevin equations for the ⁴⁸Ca+²²⁶Ra, ²³²Th, ²³⁸U, ²³⁷Np, ^{239,240,242,244}Pu, ²⁴³Am, ^245,248Cm, ²⁴⁹Bk, and ²⁴⁹Cf reactions. Our results show that with increasing dimension of Langevin equations the evaporation residue cross section is increased. Also, obtained results based on fourdimensional Langevin are in better agreement with experimental data in comparison with one-dimensional model.     

Eu~(3+)∶Y_2SiO_5晶体中电磁感应透明及群速度减慢研究

利用光与物质相互作用的密度矩阵方程及其数值解研究了Eu3 ∶Y2SiO5晶体中一种Λ型的电磁感应透明(EIT)和群速度减慢(SGV),得到了探测场的透射率随着探测场的失谐与耦合场的拉比频率的变化关系。分析了自旋子能级非均匀线宽及光学跃迁非均匀线宽对电磁感应透明和群速度减慢的影响,这两个线宽的增加不同程度地抑制了电磁感应透明。根据浓度与非均匀线宽的关系发现中心透射率不是随着浓度的增加而单调变化,而是存在一个最佳浓度使电磁感应透明的效果更加明显。分析表明对某一耦合场强,光群速度存在一极小值。