数值模式中长波辐射参数化方案的对比试验

应用一个简化的长波辐射参数方案，与郭晓岚－钱永甫模式原有长波辐射参数化方案进行了对比，分析了它们对１００ｈＰａ，３００ｈＰａ，５００ｈＰａ，７００ｈＰａ及地面的高度场、温度场、流场及气压场模拟结果的影响，发现不仅两者结果相当一致，而且与实测结果也符合很好。     

陆面特征非均一作用参数化及其对区域气候影响的数值模拟试验

本文首先改进了陆地下垫面特征非均一性的次网格尺度参数化方法，然后利用三维地气耦合的区域气候模式，设计一系列值试验，研究了下垫面特征改变对区域气候环境变化的影响，主要分析了陆地表面特征变化对我国苏南附近地区夏季温度化的影响，结果表明采用地下垫面特征非均一作用的次网格尺度参数化方法对于改进数值模拟结果的质量有一定效果。     

地形对大地电磁测深(MTS)资料的影响

地形影响是山区ＭＴＳ工作的主要问题之一，对此做了理论分析并用均匀半空间表面二维地形模型进行了数值模拟。结果表明，地形影响不仅与ＭＴＳ资料的种类、极化模式、场源电磁波周期及测点位置有关，而且还依赖于不同的地形形态和尺度。此外，还根据悬崖模型计算，给出了关于地形影响范围的估     

水汽方程积分计算格式的数值试验

该文把Easter(1993)的水汽方程积分方案应用到球坐标的有限区模式（LAFS）中，并作了两组比较性试验。Easter的方案是把水汽的通量方程和连续方程联立求解。在计算中不会出现水汽负值。而且稳定性好。本文所做的两组试验表明。Easter的方案比LAFS的原方案所产生的计算性扩散小。因而降水预报更接近实况。但CPU时间增加了16%。     

中国西北部“4.5”沙尘暴过程中尺度低压的数值模拟

利用改进型ＰＳＵ／ＮＣＡＲ中尺度数值模式（ＭＭ４标准版），取模式水平格距４０ｋｍ，４６×６１网格，垂直方向ａ取１５层，即从地面到模式顶（１００ｈＰａ），σ＝０．０、０．１、０．２、０．３、０．４、０．５、０．６、０．７、０．７８、０．８４、０．８９、０．９３、０．９６、０．９８、０．９９、１．００，采用ＮＣＡＲ的３０＇×３０＇地形资料，以常规观测资料作为初始场，较好地模拟了此次沙尘暴过程的海平面气压的演变和分布，特别是张掖、柴达本盆地以及敦煌附近的三个中低压。同时，模拟了张掖中低压与蒙古冷高压之间的甘肃河西沙尘暴东大风。敏感性试验表明，沙尘暴中低压的形成发展主要是受于物理过程制约；沙尘暴中尺度系统的研究与暴雨中尺度系统的研究是有区别的，积云对流参数化并不是特别重要，在设计研究沙尘暴的数值模式中，应当合理地处理其他的热力、动力过程及大气外强迫源的作用。模式水平格距、地形真实程度对模拟中低压的位置、中心强度有重要贡献；下垫面变化中低压强度有一定影响。张掖热低压的形成发展主要是在有利的环境形势下，特殊地势起了重要作用，表现为直接动力强迫和间接热力强迫。     

斜向坡变形破坏机制的数值模拟研究

本文讨论了斜向坡的旋转式变形和破坏机制，用有限元法计算了此种斜坡中的剩余下滑力的分布方式，并以此来解释在滑动过程中发生一定旋转的原因；论文还运用离散元法，模拟了此种滑坡的滑动过程。     

The mean and turbulence structure simulation of the monsoon trough boundary layer using a one-dimensional model withe-l ande-ε closures

An attempt has been made here to study the sensitivity of the mean and the turbulence structure of the monsoon trough boundary layer to the choice of the constants in the dissipation equation for two stations Delhi and Calcutta, using one-dimensional atmospheric boundary layer model withe-ε turbulence closure. An analytical discussion of the problems associated with the constants of the dissipation equation is presented. It is shown here that the choice of the constants in the dissipation equation is quite crucial and the turbulence structure is very sensitive to these constants. The modification of the dissipation equation adopted by earlier studies, that is, approximating the Tke generation (due to shear and buoyancy production) in theε-equation by max (shear production, shear + buoyancy production), can be avoided by a suitable choice of the constants suggested here. The observed turbulence structure is better simulated with these constants. The turbulence structure simulation with the constants recommended by Aupoixet al (1989) (which are interactive in time) for the monsoon region is shown to be qualitatively similar to the simulation obtained with the constants suggested here, thus implying that no universal constants exist to regulate dissipation rate. Simulations of the mean structure show little sensitivity to the type of the closure parameterization betweene-l ande-ε closures. However the turbulence structure simulation withe-ε. closure is far better compared to thee-l model simulations. The model simulations of temperature profiles compare quite well with the observations whenever the boundary layer is well mixed (neutral) or unstable. However the models are not able to simulate the nocturnal boundary layer (stable) temperature profiles. Moisture profiles are simulated reasonably better. With one-dimensional models, capturing observed wind variations is not up to the mark.     

Application of an optimum design technique for determining the coefficient of consolidation by using piezocone test data

For normally consolidated clay, several researchers have developed a number of theoretical time factors to determine the coefficient of consolidation from piezocone test results. However, depending on assumptions and analytical techniques, it could vary considerably, even for a specific degree of consolidation. In this paper a method is proposed to determine a consistent coefficient of consolidation by applying the concept of an optimum design technique over all ranges of the degree of consolidation. Initial excess pore pressure distribution is assumed to be capable of being obtained by the successive spherical cavity expansion theory. The dissipation of pore pressure is simulated by means of a two-dimensional linear-uncoupled axi-symmetric consolidation analysis. The minimization of differences between measured and predicted excess pore pressure was carried out by the BFGS unconstrained optimum design algorithm with a one-dimensional golden section search technique. By analyzing numerical examples and in-situ test results, it was found that the adopted optimum design technique gives consistent and convergent results.     

导水裂隙带高度预测途径探讨

导水裂隙带高度是设计防水煤（岩）柱尺寸的主要技术参数，是煤矿防治水工作的重要内容。详细论述了其研究现状及发展趋势，并提出了履岩体地质环境，履岩体力学环境和数值模拟技术三位一体的综合预测评价方法。     

Optimal weighting in the finite difference solution of the convection-dispersion equation

Oscillation and numerical dispersion limit the reliability of numerical solutions of the convection-dispersion equation when finite difference methods are used. To eliminate oscillation and reduce the numerical dispersion, an optimal upstream weighting with finite differences is proposed. The optimal values of upstream weighting coefficients numerically obtained are a function of the mesh Peclet number used. The accuracy of the proposed numerical method is tested against two classical problems for which analytical solutions exist. The comparison of the numerical results obtained with different numerical schemes and those obtained by the analytical solutions demonstrates the possibility of a real gain in precision using the proposed optimal weighting method. This gain in precision is verified by interpreting a tracer experiment performed in a laboratory column.