球谐函数法求解辐射传输方程的假散射和射线效应

假散射和射线效应是辐射传输方程近似解法中出现的特有误差.从辐射传输方程的近似求解过程出发,在定性分析的基础上构造物理模型,通过数值模拟研究球谐函数法(P1和P3近似法)的假散射和射线效应.构造激光平行入射和倾斜入射二维半透明介质的物理模型,通过内部温度场的分布特征研究假散射.构造顶部侧面保持高温而其余侧面保持低温的二维半透明介质方案,通过对比底面边界净热流密度分析射线效应.计算结果表明球谐函数法中同时存在假散射和射线效应,P3近似比P1近似减小了射线效应.同时,球谐函数法的射线效应随光学厚度的增加而减小.

球谐函数法解地球扰动引力并行计算方法

随着导弹射程的增加,传统求解地球扰动引力的球谐函数法模型阶数增大,扰动引力的计算速度明显下降。为提高弹道导弹扰动引力的赋值速度,在现有条件下结合球谐函数法计算模型的特点,在并行计算平台上,应用流水线技术和任务交叉分配的方法,建立了球谐函数法的并行计算模型。通过仿真验证,该方法能有效地提高地球扰动引力的计算速度,得到较高的并行加速比和CPU利用率。

基于球谐函数法的气溶胶辐射强迫分析

本文针对气溶胶的辐射强迫效应开展了数值方法研究。建立了基于球谐函数法的多层半透明介质辐射计算模型。并与射线踪迹-节点分析法的计算结果进行了对比验证。最后,我们使用该模型分析了北京地区的大气气溶胶辐射强迫。结果显示,本文计算出的结果与AERONET的气溶胶光学厚度有很好的关联性。

任意阶球谐函数法求解辐射传输问题的误差分析

针对各向异性散射介质内的红外辐射传输开展数值方法研究。应用球谐函数展开推导了一维吸收、发射、散射性灰介质的辐射传输近似方程组,采用差分方法结合三对角矩阵解法建立了任意阶辐射传输近似方程数值解法。并通过计算算例验证了高阶P_N方法的数值无关性问题。计算结果表明,高阶球谐函数数值解与理论解吻合得很好,同时计算验证了不同展开阶数对数值结果的影响。

超球谐函数方法及其在核三体问题中的应用

本文对超球谐函数方法及超球谐函数方法的绝热近似的数学推导作了较详细的叙述。对这两种方法的收敛问题,通过具体应用,作了比较。

自由段弹道扰动引力计算的球谐函数极点变换

为提高弹道导弹自由段扰动引力赋值速度,提出了自由段扰动引力计算的球谐函数极点变换一般方法。重点针对轨道动量矩与天球交点为极点及主动段关机点为极点两种情况,推导了考虑侧向偏差和地球旋转修正的扰动引力计算的球谐函数极点变换模型。仿真分析表明:从落点偏差来看,前者扰动引力计算相对误差为1.89%,后者为6.07%,满足诸元和制导精度要求。

用球函数方法计算东亚测站的综合海潮负荷参数(英)

用球谐函数方法（SHM）计算了由Schwiderski海潮图和日本海海潮图引起的M2，S2，K1和O1各分潮波的综合海潮负荷参数.计算中采用了Pagiatakis（1990）地球模型的负荷Love数.结果表明SHM具有更高的计算效率和计算精度．为了真正实现±毫米的VLBI地壳运动测定精度和±0.1微伽的绝对重力测量精度，对未来海潮图模型，特别是日本海的海潮图模型的精度要求作了讨论．计算结果还表明不同地球模型的负荷Love数对计算的测站位移之差为0.7毫米，在目前测量精度下尚可以不作考虑.中国近海的海潮M2波对台北，上海，Daejeon,Okinawa,和Kagosima站的径向和水平分量位移影响可分别达±22．4～±1．0毫米和±05-±5.5毫米，表明中国近海高精度的海潮图模型对东亚测站的地壳运动高精度测量是十分重要的．

地壳垂直形变场逼近的球冠谐分析法

基于地壳垂直运动的均衡理论,讨论了球谐函数法用于区域地壳垂直形变场逼近时所存在的问题,提出并推导了球冠谐分析法逼近区域地壳垂直形变场的模型。并利用该模型对两个实际区域进行了试验,取得了良好的效果。

二维参与性介质辐射换热的球谐有限元法

采用球谐函数与有限元相结合的方法(即球谐有限元法),对二维介质内辐射换热问题进行研究。以规则和非规则形状几何体辐射问题为例,使用球谐有限元法可得到较准确的计算结果,且与其它数值计算方法(如球谐有限差分法、离散坐标有限元法等)相比具有以下优势:球谐有限元法方法简单,计算速度快,适用于大型工程应用;球谐有限元法可用于计算各种形状的几何体,计算应用范围广。

全球环流异常特征及其与中国气候变化的关系

采用球谐函数展开方法对1948—2009年NCEP/NCAR(美国国家环境预报中心/美国国家大气研究中心)高度距平场再分析资料进行分析,通过"时间-高度"的二维坐标方式考察环流异常的立体结构特征,发现南、北半球大气环流在1980年前后发生了明显的正负位相(槽脊)突变.与中国冬季温度和夏季降水的相关分析表明:前期春、夏、秋季的北半球环流纬向正异常越强,华北地区的冬季气温越高,反之越低;前一年冬季环流异常与华北—河套地区的夏季降水呈现反相关,春季北半球均匀异常环流与华北及西南夏季降水呈现反相关;大气环流的纬向异常与华北—河套地区的夏季降水呈现正相关,而同期环流异常与华北东部、西南及江南北部的降水呈现负相关.

基于北斗导航系统的我国地基增强格网模型研究

基于全国均匀分布的210个北斗地基增强基准站观测数据,提出基于非差非组合精密单点方法获取空间相关的延迟误差及残差,研究利用反距离加权法、克里金插值法以及球谐函数法建立了中国区域的三种地基增强格网模型。使用单频伪距单点定位方法对中国区域的北斗地基增强格网模型定位精度进行了评估,验证结果表明本文建立的格网模型的定位精度达到米级,可满足大众用户的需求。

SP_3方程的非线性迭代解法研究

SP_3方程是简化球谐函数法三阶展开的中子输运方程,比扩散方程的计算精度更高。本文采用非线性迭代的方法求解SP_3方程,并且提出了带角度离散的耦合修正关系式,解决了SP_3方程在真空边界条件下耦合修正关系复杂和非线性迭代中二阶中子注量率修正数值不稳定的问题。基于该方法编制非线性迭代求解SP_3方程的计算程序NLSP3,对基准题进行了计算。数值结果表明,本文所采用的非线性迭代求解SP_3的方法是正确的,且具有较好的计算精度。

基于平面通量展开的三维Pin by Pin输运SP3计算

基于Pin by Pin输运SP3计算是下一代反应堆物理计算方法中一个重要的研究方向。节块法一般采用高阶通量和中子源展开以适应较大的节块尺寸,因此对于Pin by Pin计算来说效率偏低。由于反应堆堆芯不均匀性更多发生在径向,本文提出一种径向基于二维平面通量展开结合轴向常规节块法的综合方法进行三维Pin by Pin输运SP_(3)计算,其三个方向的节块耦合迭代均采用基于中子净流耦合的格式。通过IAEA 2D/3D基准问题和典型输运基准问题计算,验证该计算方法对于堆芯扩散计算和输运计算的正确性和可行性。