一种基于特征理论模拟凝聚炸药爆轰的单元中心型Lagrange方法

提出一种数值模拟凝聚炸药爆轰问题的单元中心型Lagrange方法.利用有限体积离散爆轰反应流动方程组,基于双曲型偏微分方程组的特征理论获得离散网格节点的速度与压力,获得的网格节点速度与压力用于更新网格节点位置以及计算网格单元边的数值通量.以这种方式获得的网格节点解是一种"真正多维"的理论解,是一维Godunov格式在二维Riemann问题的推广.有限体积离散得到的爆轰反应流动的半离散系统使用一种显-隐Runge-Kutta格式来离散求解：显式格式处理对流项,隐式格式处理化学反应刚性源项.算例表明,提出的单元中心型Lagrange方法能够较好地模拟凝聚炸药的爆轰反应流动.     

油气藏三维三相组分模型软件研制与应用

提出了一种对压力和饱和度隐式、组分摩尔分数显式的部分隐式组分模型求解方法。该方法所需联立求解的方程组数量少,且兼有全隐式收敛性好的特点。通过与引进软件用相同地质模型、PVT资料对柯克亚X₅¹(3)凝析气藏循环注气进行模拟对比研究,验证了自编的部分隐式组分模型软件计算结果的可靠性,标志着自主版权的组分模型软件已经开始走向成功。     

McMillan Tc公式的解析推导(Ⅲ)

本文把作者在前面两篇文章导出的T_c公式推广成下面形式:T_c=αω_log(ω_log/ω_c)^{(μ*/(λ-μ*))exp{-(1+λ)/(λ-μ*)},并从线性Eliashberg方程出发,导出了计算α的方程组。α一般是λ和μ*的函数。在弱耦合极限下,由上述方程组解得,α=2γ/π,其中lnγ=C=0.5772是Euler常数。这个结果表明了,前面两篇文章得到的T_c公式在弱耦合极限下是正确的。作者进而在Einstein谱和μ*=0情形,用数值计算方法从定α的方程组算出当λ=0.23,0.25,0.38和0.48时,a的数值。结果表明,至少在0.23≤λ≤0.45区间中,α变化很小,近似等于1/1.30。此时,本文的T_c公式实际上就是Allen及Dynes修改后的经验的McMillan T_c公式。 关键词：}     

40Ca20是否稳定?

本文计算表明通常认为高稳定的双满壳层核⁴⁰Ca₂₀其实并不稳定。⁴⁰Ca₂₀将可能通过双K层电子俘获的机制而转化为⁴⁰Ar₁₈,其半衰期约为1.2×10³³年。 关键词：     

用高分辨差分格式解SNS方程

本文从Sweby的显式TVD格式^[2,3]出发,构造发一个推进方向一阶精度、侧向二阶精度的隐式TVD差分格式,并用于解可压缩流定常简化Naviev-Stokes(SNS)方程。对不出现流向分离的平板边界层和激波—边界层相互干扰两个问题的数值计算表明,方法可以有效捕捉流场中的间断。     

在LHC上寻找希格斯粒子

 大型强子对撞机(Large Hadron Collider,简称LHC)2010 年3 月开始运行,对撞能量为7TeV,这是设计能量14 TeV 的一半。瞬时亮度从10²⁷cm^-2s^-1迅速上升,2010 年底达到10³²cm^-2s^-1,2011 年底达到10³³cm^-2s^-1,积分亮度有5/fb, 超过原来的预期。LHC 进入出成果阶段。     

石英基掺Tm3+光纤激光器特性的理论分析

基于速率方程理论建立了石英基掺Tm³⁺光纤激光器³H₆—³H₄抽运方式的理论模型,采用龙格-库塔 (Runge-Kutta) 法对该模型进行了数值分析,得出抽运光与激光功率在掺Tm³⁺光纤中的分布特性,并对Tm³⁺光纤长度、Tm³⁺的掺杂浓度等因素对掺Tm³⁺光纤激光器性能的影响进行 关键词： 3+光纤激光器')" href="#">石英基掺Tm³⁺光纤激光器 龙格-库塔法 抽运效率 掺杂浓度     

考虑量子效应的Zakharov方程组的孤波解

借助Maple程序，利用扩展的双曲函数法和双函数法求解了考虑量子效应的Zakharov方程组，得到了多种孤波解，其中包括亮孤波解、W型孤波解、M型孤波解和奇性孤波解. 关键词： 量子效应 Zakharov方程组 扩展的双曲函数法 孤波解     

Davey-Stewartson方程组的包络周期解和孤立波解

应用Jacobi椭圆函数展开法,求得了Davey-Stewartson方程组的包络周期解和孤立波解． 关键词： Davey-Stewartson方程 Jacobi椭圆函数 包络周期解 孤立波解     

非线性长波方程组和Benjamin方程的新精确孤波解

在齐次平衡法、双曲正切函数法和辅助方程法的基础上引入一个新的辅助方程,并借助符号计算系统Mathematica来构造了非线性长波方程组和Benjamin方程的新精确孤波解, 这种方法也可用于寻找其他非线性发展方程的新的孤波解. 关键词： 新的辅助方程 非线性长波方程组 Benjamin方程 孤波解     

基于JASMIN三维势场快速多极子算法的并行实现

在JASMIN上,基于进程/线程两级并行实现策略,研制三维Laplace核函数FMM(fast multipole method)的解法器模块"JASMIN-3DLapFMM".该解法器已成功应用于三维静电场远场势的并行计算.固定单机问题规模,在上万个处理器核上运行百亿粒子的大规模问题,获得进程级几乎线性的并行可扩展性.固定总的问题规模和1024个进程,4个线程时,获得大约3倍的加速.     

求解Burgers方程的一种交替分段隐式算法

首先提出一个新的求解Burgers方程的差分格式,然后在此差分格式的基础上构造了便于并行计算的交替分段隐格式,并作了线性化稳定性分析.数值结果表明,本方法具有较高的精度,尤适于扩散项系数较小时的计算,且有效避免了数值结果的非物理振荡.     

改进的2PPJ格式

提出了解线性代数方程组的一种方法,称为改进的2PPJ(Modified two parameters parallel Jacobi-type)格式,简记为M2PPJ.讨论了它的收敛性和最优参数的选取.收敛速度比2PPJ格式快一倍以上.给出的数值结果说明本方法的优越性.     

激光等离子体二维相对论电磁粒子模拟

研制了激光等离子体二维相对论电磁粒子模拟程序(2DCIC)。追踪几万甚至百万个模拟粒子在外加激光场和自洽场中运动,自洽地计算电荷和电流密度,求解完全的Maxwel方程,电子的相对论运动方程和离子的牛顿运动方程,辅以灵活的诊断研究波-波,波-粒子相互作用的发生、发展和饱和的细节以及时间演化规律。激光可以正入射,也可以斜入射;等离子体可以是均匀密度,也可以具有密度梯度;为了节约机时,还发展了并行运算。物理模型参数可调,既适用于研究激光聚变等离子体相互作用,也适用于超短脉冲超强激光等离子体相互作用和其它等离子体问题。经过多次试算检验,对等离子体平衡态进行了计算研究,对于超短脉冲超强激光的传播也进行了初步模拟计算。     

Adaptation of homotopy-perturbation method for numeric-analytic solution of system of ODEs

A new reliable algorithm based on an adaptation of the standard Homotopy-Perturbation Method (HPM) is presented. The HPM is treated, for the first time, as an algorithm in a sequence of intervals (i.e., time step) for finding accurate approximate solutions of linear and nonlinear systems of ODEs. Numerical comparisons between the Multistage Homotopy-Perturbation Method (MHPM) and the available exact solution and the classical fourth-order Runge-Kutta (RK4) method reveal that the new technique is a promising tool for linear and nonlinear systems of ODEs.     

A Simple Approach to Study the Isospin Effect in Mass Splitting of Three-Nucleon Systems by Using Hyperspherical Functions

In this work, the binding energy and wavefunctions of three-nucleon systems are obtained by using hyperspherical harmonic approach. We have used a mathematical modification method to obtain the eigenvalues and eigenfunctions of Schrödinger equation for three-nucleon systems in calculation. Next, we have used a simple approach to obtain the difference between binding energy of ³H and ³He where gives us mass splitting of three-nucleon systems. We have compared our results with the other works and experimental values.     

用B-splines技术计算超强平行外磁场中氢分子离子的能级

考虑两核之间的相互作用,采用单中心展开方法和B-splines技术构建氢分子离子的径向和角向波函数,计算超强平行外磁场中氢分子离子的能级E_b(m=0,-1,-2,-3,-4,-5)和两核间的平衡距离2R.计算中,磁场强度从2.35×10⁶T变到2.35×10⁸T.对比绝热变分近似计算结果,计算精度提高的范围在0.00015 Ry~0.0488 Ry之间;对比变分近似计算结果,计算精度提高的范围在0.0029 Ry~0.0257 Ry之间.这些计算结果表明,用二重B-splines基函数展开方法使氢分子离子的波函数的完备性更好,并能得到比绝热变分近似和绝热近似方法更高的计算精度.     

CsI(Tl)晶体的$lt;i$gt;α$lt;/i$gt;/γ 波形数字化实时甄别

基于CsI（Tl）探测器对α/γ 粒子的波形甄别能力,采用电荷比较法设计了一种波形实时甄别系统. 介绍了实时甄别系统的设计原理,利用⁶⁰Co-γ源、²⁴¹Am-α源对实时系统进行了甄别实验,探究了不同参数对甄别效果的影响,并给出了最优甄别效果下的参数设置. 研究表明,设计的数字化实时波形甄别系统体积小,能准确、实时地甄别开α/γ粒子,最佳品质因子大于1.4,事件计数率可达3× 10⁵/s. 关键词： 数字化波形甄别 实时 电荷比较法 CsI(Tl)晶体     

IS 40Ca A STABLE ISOTOPE?

It is shorn that the double magic nucleus ⁴⁰Ca is unstable with respect to the double K-electron capture ⁴⁰Ar. By using the shell model wave functions for both ⁴⁰Ca and ⁴⁰Ar, the lifetime of ⁴⁰Ca is estimated as 1.2×10³³yrs.     

Numerical analysis of motional mode coupling of sympathetically cooled two-ion crystals

A two-ion pair in a linear Paul trap is extensively used in the research of the simplest quantum-logic system; however,there are few quantitative and comprehensive studies on the motional mode coupling of two-ion systems yet. This study proposes a method to investigate the motional mode coupling of sympathetically cooled two-ion crystals by quantifying three-dimensional(3 D) secular spectra of trapped ions using molecular dynamics simulations. The 3 D resonance peaks of the ~(40)Ca~+ – ~(27)Al~+ pair obtained by using this method were in good agreement with the 3 D in-and out-of-phase modes predicted by the mode coupling theory for two ions in equilibrium and the frequency matching errors were lower than 2%.The obtained and predicted amplitudes of these modes were also qualitatively similar. It was observed that the strength of the sympathetic interaction of the ~(40)Ca~+ – ~(27)Al~+ pair was primarily determined by its axial in-phase coupling. In addition,the frequencies and amplitudes of the ion pair's resonance modes(in all dimensions) were sensitive to the relative masses of the ion pair, and a decrease in the mass mismatch enhanced the sympathetic cooling rates. The sympathetic interactions of the ~(40)Ca~+ – ~(27)Al~+ pair were slightly weaker than those of the ~(24)Mg~+– ~(27)Al~+ pair, but significantly stronger than those of ~9Be~+ – ~(27)Al~+ . However, the Doppler cooling limit temperature of ~(40)Ca~+ is comparable to that of ~9Be~+ but lower than approximately half of that of ~(24)Mg~+. Furthermore, laser cooling systems for ~(40)Ca~+ are more reliable than those for ~(24)Mg~+and ~9Be~+ . Therefore, ~(40)Ca~+ is probably the best laser-cooled ion for sympathetic cooling and quantum-logic operations of ~(27)Al~+ and has particularly more notable comprehensive advantages in the development of high reliability, compact, and transportable ~(27)Al~+ optical clocks. This methodology may be extended to multi-ion systems, and it will greatly aid efforts to control the dynamic behaviors of sympathetic cooling as well as the development of low-heating-rate quantum logic clocks.