关于正确处理对流－弥散方程和渗流方程中源汇项的若干问题

对流─弥散方程和渗流基本微分方程中的源汇项，分别是决定含水层中任一内点上水质量和溶质质量时变率的重要因素之一。而关于源汇项不适当的乃至是错误的表述，在文献中却不乏其见，其影响之大是显而易见的．本文深入分析了某些表述的不当或错误之处，给出了正确的处理与表述方法。     

高效优化非分裂PML边界二阶标量波方程数值模拟方法

卷积完全匹配层（convolution perfectly matched layer,CPML）吸收边界是一种高效处理波动方程数值模拟中人工边界反射波的方法。本文基于传统的一阶系统CPML吸收边界条件推广并推导了新的二阶系统CPML边界条件（NCPML）。与常规二阶系统CPML边界条件不同,新边界条件推导的核心思想是在复数-频率域中忽略部分衰减因子空间变化特性,避免其在时间域产生复杂卷积算子,然后反变换至时间域得到基于CPML吸收条件的二阶标量波方程,并应用于二阶标量波方程数值模拟。均匀介质模型测试验证了NCPML吸收条件在内存使用上相对于常规二阶系统CPML与SPML（split PML）吸收条件更少。在对人工边界反射的吸收效果上,NCPML稍逊色于常规二阶系统CPML,但二者均相对于SPML优势明显。最后通过层状模型和Marmousi模型测试验证了NCPML的稳定性及其在效率上的优势。     

风暴尺度集合预报最优侧边界条件扰动方法设计:个例分析

考虑到全球预报模式与风暴尺度预报模式在分辨率上的显著差异,在构造风暴尺度集合预报系统的时候需要用一个中间分辨率的中尺度区域模式为风暴尺度模式提供侧边界条件扰动,但如何构造侧边界扰动才能更为有效地提高风暴尺度集合预报系统的预报能力目前仍然未知。本文基于WRF模式,通过一次个例试验设计了风暴尺度集合预报中的3种不同侧边界扰动方案,结果表明:直接通过0.5°水平分辨率全球集合预报扰动插值所得到的侧边界扰动(LBC_DOWN)在预报中可以获得较高的大尺度扰动能量,而在中尺度区域模式(本文中为模式外层)中通过ETKF循环所构造的侧边界条件扰动(LBC_CYCLE)包含较高的中小尺度能量,而将LBC_CYCLE中的中尺度扰动信息与LBC_DOWN中的大尺度扰动信息相混合所得到的混合侧边界扰动(LBC_BLEND)在大尺度能量上更接近于LBC_DOWN,在中小尺度能量上更接近于LBC_CYCLE;LBC_BLEND较前两种方案有着更好的离散度技巧表现;在降水概率预报技巧方面LBC_BLEND与LBC_CYCLE较为接近,且均优于LBC_DOWN。     

一个近岸海流数值计算模式的开边界条件

为了计算的省时性,ＩＡＰ近岸海流数值模式采用了分解算法,即将控制方程分解为３个过程并用不同的时间步长进行积分。章在此基础上对分解后的控制方程作局地线性化,得到了适合不同过程的开边界条件,对于适应过程,控制方程变换为等价的表示沿不同方向传播波动的特征方程。传出计算区域的波动用特征方程来描述,而传入波动则由无反射边界条件来消除。对于演变过程,求出了解析形式的通解。在流出点上,边界条件可以借助解析解由     

基于正交贴体网格的VTI介质地震波数值模拟

为了解决规则网格有限差分法在处理起伏地表自由边界条件时需进行复杂坐标转换和插值运算的问题,以及对起伏地表进行阶梯状离散近似所产生的虚假散射波问题,将计算流体力学中的正交贴体网格生成技术引入到起伏地表下VTI介质的网格剖分中,采用基于同位网格的DRP/opt MacCormack有限差分法模拟起伏地表下VTI介质中的地震波场。网格的正交性使实施自由边界条件时无需做复杂的坐标转换和插值运算。数值算例表明,本文提出的网格剖分法能够有效消除阶梯状网格产生的虚假散射波,从而提高起伏地表下VTI介质地震波数值模拟精度。     

被动桩水平位移荷载施加位置的探讨

被动桩是指一种由于土体水平位移而发生挠曲变形的桩。目前不论是室内试验还是数值计算,常在模型边界对土体施加水平位移荷载（位移边界条件）形成土体的位移场,用来研究被动桩的变形特征,但不同研究者采用的水平位移施加位置（边界位置）与被动桩的距离存在较大差异。结合一个工程案例,采用有限单元法对水平位移加载位置进行了单变量参数分析。研究表明,随着水平位移加载位置与被动桩距离的增加,桩身变形显著减小。这表明了施加位移边界后地基中存在明显的应力扩散。当该距离小于5.5倍桩径时,在模型边界施加同等大小的水平位移所需应力显著增加。同时讨论了水平荷载加载位置与被动桩距离的合理范围,认为应在5～8倍桩径范围之内。     

三角网格谱元法地震波场数值模拟

谱元法结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度,是近年来发展的一种重要的地震波场数值模拟方法.经典的谱元法采用四边形(六面体)网格,利用一维Gauss-Legendre-Lobatto(GLL)积分的张量积得到对角的质量矩阵,以大大提高计算效率,但是四边形(六面体)网格不能够灵活地刻画复杂的几何模型的弯曲界面.为此,在谱元法中引入三角形(四面体)网格到二维(三维)是十分必要的.不同于经典的谱元法,在非结构化网格中不能使用GLL积分的张量积,使得非结构化网格的谱元法的实现存在着诸多的困难.目前,比较流行的三角网格谱元法,通过使用KoornwinderDubiner(KD)正交多项式,并正交化这些KD多项式构建基函数,同时利用重合的插值节点和积分节点以获取对角的质量矩阵;它所使用的积分点为优化的点集——Fekete点,且这些积分点能与四边形网格完全耦合.相比于四边形,三角网格谱元法能显著提高复杂模型的描述能力,对起伏地表模型有很大优势.本文引入高效的最佳匹配层(PML)吸收边界条件,并通过数值试验将三角网格谱元法与经典的谱元法进行对比研究.相比于经典的谱元法,三角网格谱元法显著缺点为较低的计算精度.对于7阶谱元,为了能够精确地模拟面波,三角网格谱元法需要在每个最短的面波波长内至少有11个采样点,然而经典的谱元法仅需4个采样点,并且前者所需的内存量约为后者的5.5倍.     

上边界条件对多年冻土地温场数值模拟结果的影响分析

以玛多地区多年冻土为背景,建立多年冻土地温场的数值计算模型,以不同的方式考虑近 60 a 来的气温变化构成不同的上边界条件,通过模型计算分析不同上边界条件下的不同时期温度场、未来冻土退化特征.结果表明：在上边界条件中采用气象站实测近60 a波动温度值和采用近60 a平均恒定值时,浅层冻土地温差异明显,且越浅层地温与越近时间的上边界条件相关.预测未来100 a冻土地温变化趋势发现,相同升温速率和升温初始温度条件下,上边界采用实测60 a波动温度值对冻土退化过程影响较小;升温初始温度值提高到与趋势线衔接后,冻土退化起始时间从约第45年提前到约第20年;60 a实测温度和升温初始温度值均提高到与其初始温度场上边界条件衔接后,冻土退化起始时间从约第20年提前到约第15年;冻土退化从开始到完全退化经历时间为25 a左右.     

基于GPU并行加速的叠前逆时偏移方法

为了提高复杂地下介质的成像精度和偏移算法的计算效率,提出可高效对地下复杂构造进行准确成像的GPU加速叠前逆时偏移方法.该方法采用双程声波方程进行波场延拓,突破倾角限制,借助于高阶有限差分方法实现叠前逆时偏移成像;利用GPU(Graphic Processing Unit)并行加速技术对波场延拓和成像进行计算,相比于传统算法,其计算效率有较大提高,可以解决叠前逆时偏移算法计算量过大问题;在获取波场信息过程中,也采用随机边界条件,实施以计算换存储策略,解决逆时偏移计算中的海量存储问题.模型测试结果表明,该方法能够高效和高精度地对地下复杂地质体成像.     

浅析发动机冷却系统

一、现代发动机冷却系统的特点 传统冷却系统的作用是可靠地保护发动机,而还应具有改善燃料经济性和降低排放的作用。为此,现代冷却系统要综合考虑下面的因素：发动机内部的摩擦损失;冷却系统水泵的功率;燃烧边界条件,如燃烧室温度、充量密度、充量温度。