Mark4格式器升级为Mark5B采样器的方法及步骤

近年来世界VLBI（Very Large Baseline Interferometry）终端技术发展迅速,主流的记录终端已发展到Mark5B记录系统,Mark5C也正在研制和测试阶段。为了适应未来高速率、海量数据以及数字基带转换器观测的要求,需要将现有的Mark5A升级到Mark5B或Mark5B＋,相应地需要将Mark4格式器升级为MarkB采样器。主要描述了相关设备的原理以及Mark4格式器升级的步骤。     

非均匀层状介质一维波动方程精确解的有限差分算法

平面波的传播问题通常可以归结为一维波动方程的定解问题。在非均匀介质中,即使简单的一维波动方程也需要借助于数值方法获得近似解。3层5点古典差分格式是计算偏微分方程一种常用算法,作为一种显式迭代格式,需要满足稳定性条件 ,其中 为波速, 为空间采样间隔, 为时间采样间隔。当 时, ,古典差分格式达到临界稳定状态。在这种情况下,平面波在 时间内的传播距离恰好等于空间采样间隔,差分格式真实地反映了平面波的传播原理,因而可以得到一维波动方程的精确解。但是,由于在非均匀介质中存在不连续的波阻抗界面,此方法不适于计算非均匀介质的波场。为了将临界稳定情况下的古典差分格式推广应用至非均匀层状介质,提出了一种能够处理波阻抗界面的有限差分格式,并应用傅里叶分析法得到其稳定性条件。模型算例验证了此算法的正确性。     

《甘肃地质学报》征稿简则

《甘肃地质学报》是由甘肃省地矿局、甘肃省地质学会共同主办的国内外公开发行的学术性刊物。以刊登本省及相邻省 (区 )地质行业最新科研成果和论文为主 ,同时对普查、勘探、水文地质、物化探、探矿和测绘、测试、数学地质等方面的新技术、新方法加以报导。欢迎业内人士踊跃投稿。来稿要求　　 ( 1)文稿应有创新 ,立论科学 ,文题明确 ,数据可靠 ,词语准确、规范 ,文句精炼 ,标点正确。文稿主篇以 80 0 0字以内为宜 (包括插图、插表 )。　　 ( 2 )文稿应包括 :题名、作者姓名、工作单位、邮政编码、中图分类号、文献标识码、中、英文摘要、关…     

基于Qt的空间数据存储格式检查软件的设计与实现

面对日益增长的空间数据,网络环境较差的地图、地理信息应用服务系统面临着存储容量有限、地理空间数据格式不开放的挑战。而现有的许多地理空间数据存储依赖于平台,还需要克服跨平台带来的问题,亟需研究设计统一的数据存储标准来解决空间数据格式不一致、不同GIS平台数据转换过程复杂等问题。考虑到不同系统GIS平台、GIS应用程序间数据格式的差异,本研究基于GeoPackage（GPKG）数据容器存储标准,采用Qt 5.2.0+VS2010开发平台,C++编程语言,通过三层架构编程模式研发了GIS数据存储格式的检查软件,以解决数据格式不统一的问题。该软件利用批量数据库查询,快速地实现了核心功能项、可选功能项、注册扩展功能三大类等超过百项功能文件格式的检查,并输出详细的检查结果,为GIS数据格式存储标准化提供支持和帮助,可有效解决不同GIS平台、不同GIS应用程序数据共享时面临的格式不一致问题,简化了数据转换的步骤。     

IAP AGCM-4动力框架的积分方案及模式检验

以灵活性跳点格式和时间分解算法为重点, 介绍了IAP (中国科学院大气物理研究所) AGCM-4的动力框架, 并用R-H (Rossby-Haurwitz) 波和Held-Suarez方案对框架进行了检验。结果表明: 高纬灵活性跳点格式很好地保持了普通跳点格式的各种性质, 未产生计算紊乱现象, 可在一定程度上增大时间步长, 在未采用滤波处理时增大时间步长的效果更为显著; 时间分解算法, 即将平流过程与适应过程分开计算, 二者都采用3次非线性迭代积分方案, 在N=5和N=10 (N为平流过程时间步长与适应过程时间步长的比值) 的情况下, 可分别节省CPU时间10.7%和19.9%; R－H波检验表明, 在积分的前80天, 框架较好地保持了R－H 4波的波形, 总有效能量仅衰减0.1%, 从第80天开始, 风场波形开始变形和破碎, 对应的动能和总能量也开始急剧衰减, 至第365天, 总有效能量衰减约8%, 此时, 风场和高度场亦变得较平直。经分析, 第80天R－H波破碎的主要原因是平流项的计算误差使R－H波失稳, 而其后的高度场和风场的变平直是典型的旋转适应机制; Held－Suarez方案检验也证明该框架是可靠的。     

非均匀网格下的高阶精度中央差分格式:理论推导和理想试验

传统的高阶精度有限差分格式通常是在均匀网格的基础上推导得到的,在非均匀网格的情况下它会出现精度退化的问题。基于泰勒展开方法构造了一种适用于非均匀网格的2阶、4阶和6阶精度中央有限差分方案,利用Burgers方程和一维平流方程对新方案的性能进行测试,着重分析新方案对其误差大小及分布形态的改进效果。数值模拟结果表明:在非均匀网格下,提高差分方案的精度可明显减小数值解误差(降低了70%~88%),特别是当差分精度从2阶提高到4阶的时候。同时,高阶精度方案在梯度变化较大或者网格距较粗区域的模拟结果更有优势,4阶和6阶精度方案在以上区域的误差远小于2阶精度方案。方案可用于提高数值天气预报模式中非均匀分层模式的垂直差分计算精度。      

扩散系数反演及其差分格式研究

空气污染预报属于正问题，而从污染物浓度来求解扩散系数则属于反问题。正问题和反问题有着本质的不同，在解的定义和求解方法上也有很大的区别。从最优控制的角度定义了大气边界层中垂直扩散系数反演问题的解，用伴随模式方法得到目标函数的梯度并求解反问题。研究中发现，反演的结果与模式差分格式的选取有关，与测试源的设置也有直接的关系。在经过多次数值试验后，对于误差的来源进行了理论分析，发现了反演结果与差分格式及测试源之间的联系，得到了满意的反演结果，并为实验测定扩散系数提供了依据。     

二阶Adams-Bashforth格式的垂直敏感性分析

蛙跳(leapfrog)时间差分格式采用Asselin-Robert时间滤波方案去除计算解能够降低原始方程组的时间差分格式的计算精度,采用二阶Adams-Bashforth格式构造的欧拉前差方案可弥补蛙跳格式的不足,即:在不存在计算解的条件下去除滤波的影响,更大程度的保持方程组的计算准确性。本文基于NCAR CAM3.0(Community Atmosphere Model 3.0)完善的软件平台,将原模式的三时间层蛙跳时间差分方案修改为两时间层二阶Adams-Bashforth时间差分格式,对与重力波有关项使用中央差隐式处理,以此构建半隐式大气环流谱模式。利用动力检验的方法探讨模式对垂直分辨率的敏感性,从而寻找模式在较小计算代价下提高计算效果的可能。通过斜压波实验发现,提高垂直分辨率使模式具有更强的斜压波模拟能力,其模拟效果甚至已经与提高水平分辨率的效果相当,可以作为一种弥补模式运算效率不足的可行方案。     

高阶Runge-Kutta-Li算法对二维线性平流方程的计算检验

利用高阶Li空间微分方案(Li, 2005),实现了时间积分为3～6阶Runge-Kutta-Li(RKL)格式的求解算法。二维线性平流方程的试验结果表明:在计算稳定的条件下,各阶算法的计算误差随时间的推移基本上是线性增加的。非转动背景场的平流算例中（高斯型的初值）,高阶RKL算法可以取得较好的计算效果。与3、4、5、6阶RK算法配合的Li空间差分方案有效阶数可以达到5、7、9、10阶。RK 算法的阶数为5（6）阶时,总误差控制在10-7（10-8）以内。随RK阶数增加Li微分的有效阶数有增加趋势,且总误差逐渐减小。定常转速的背景场算例中（偏心的高斯型初值）,当RK阶数为3时,最优空间差分阶数为10;相应的阶数为4、5、6时对应的空间最优阶为16,22,22,总计算误差可以控制在10-15～10-16。随着精度的提高,误差的绝对值减小很迅速,说明算法是非常有效的。对于圆锥型初值（定常转速的背景场）,4、5、6阶RK算法和3阶算法的效果差不多。高阶算法对此类具有导数不连续点的算例,效果不如高斯初始场好,结果不能保持正定,有些地方误差出现下冲和上翘。随着空间差分精度的提高,非正定的解数量和数值减小,误差的绝对值减小,说明了算法在一定程度上是有效的,但并不适合追求极高的算法阶数。这与谱方法中的导数不连续问题有些相似,误差的产生主要源于导数的不连续性,差分类方法仅能获得与导数连续性阶数相当的算法精度。各种算例中,采用恰当的边界条件是必要的,例如旋转背景场算例,比较适合使用无穷远边界条件,否则会出现计算不稳定或无法将计算误差控制到较小的范围内。     

利用算子分裂迎风均衡格式解对流为主溶质运移问题

水污染模拟问题是水流问题与溶质运移问题的耦合问题.各种常见的数值解法在以对流为主溶质运移问题的求解中都会遇到困难,如用有限单元法或有限差分法时,会产生数值弥散与过量这两类误差.引入算子分裂迎风均衡格式法求解对流为主的水污染模拟问题,较好地克服了数值弥散和数值解出现振荡问题,该格式具有良好的稳定性、单调性及守恒性特点.