3D裂隙网络随机模拟及其工程应用研究

岩体中存在的大量随机分布的裂隙对岩体的工程性质及稳定性有着极其重要的影响。传统的工程地质研究方法难以客观地描述其空间分布特征。应用概率与统计学的理论 ,通过统计分析得到裂隙的分布规律性 ,采用随机模拟的方法实现与统计分布相适应的裂隙网络是研究随机裂隙宏观特征的有效途径。通过模拟所得的 3D裂隙网络 ,进一步研究了裂隙的连通特征及工程开挖面上的稳定分析与加固措施优化的方法     

景像匹配仿真的一种新方法

分析了影响景像相关匹配精度的各种误差因素，并给出它们的误差仿真模型，该模型可以适应飞行姿态、飞行器特性，天候等条件变化时的匹配概率估计，并可以与图像特征参数分析方法相结合。实验验证了该模型的有效性，且可以快速地给出影像图的适配性估计图，用于飞行器的航迹规划。     

包络地形和重力波拖曳对区域气候模拟效果的影响

用NCAR区域气候模式(RegCM2)的数值试验，研究了包络地形和地形重力波拖曳作用对1991年5、6月份中国区域内气候模拟效果的影响。通过试验结果的对比分析发现，在模式中引入地形重力波拖曳作用可在一定程度上改善模式对气候系统和要素的模拟效果，使模拟结果更符合气候实况。包络地形对模拟结果也有一定的改善，但对于细网格的区域气候模式(RCM)来说，其效果不如重力波拖曳明显。不论是地形重力波拖曳还是包络地形，在模式积分的初期，它们的作用并不明显，随着积分时间的增长，它们对模拟结果的影响程度增大。     

利用GCM模式对全球雷电活动的模拟分析

将修正的参数化方案嵌套入美国马里兰大学海-陆-气研究中心的全球大气环流谱模式(COLAGCM),建立了可用于全球雷电活动模拟与估测的模式系统.通过对模式输出与实况资料的对比,检验其气候模拟性能.虽然模式的分辨率较低但经过较长时间积分,仍可再现实况资料的主要特征.性能良好,计算稳定,且因时间步长较大可大大节省机时.耦合模式系统经长时间的积分后,模拟的多年平均的全球雷电频次时空分布及变化与观测分析较一致.全球范围内每秒钟平均有79次的雷电发生,陆海比为12.2.多发生在北半球夏季且具明显的季节变化特征.     

地理信息系统中仿真接口的设计与实现

讨论了以地理信息系统作为建模与仿真过程的空间信息支持 ,在高层体系结构 (HLA)提供的RTI服务基础上 ,运用多线程技术实现地理信息系统作为仿真成员的HLA接口的设计与实现。实验表明 ,能够满足仿真的需求     

非线性最小二乘参数平差迭代算法

在非线性最小二乘问题现有的3类主要算法－－高斯－牛顿法、阻尼最小二乘法和最小二乘的拟牛顿法的基础上，引入了综合性能更优的非线性规划的SQPM（序列二次规划法）算法，并且为进一步提高SQPM算法迭代的收敛性，对其步长策略进行了改进。改进的SQPM算法成为无需精确计算参数概略值的非线性最小二乘参数平差的实用和有效算法。     

北极海冰的厚度和面积变化对大气环流影响的数值模拟

文中利用中国科学院大气物理研究所设计的两层大气环流模式 ,模拟研究了北极海冰厚度和面积变化对大气环流的影响 ,尤其是对东亚区域气候变化的影响。模式中海冰厚度处理趋于合理分布 ,导致东亚冬、夏季风偏强 ,使冬季西伯利亚高压和冰岛低压的模拟结果更趋合理 ;另一方面 ,海冰厚度变化可以激发出跨越欧亚大陆的行星波传播 ,在低纬度地区 ,该行星波由西太平洋向东太平洋地区传播 ;海冰厚度变化对低纬度地区的对流活动也有影响。冬季北极巴伦支海海冰变化对后期大气环流也有显著的影响。数值模拟结果表明 :冬季巴伦支海海冰偏多 (少 )时 ,春季 (4～ 6月 )北太平洋中部海平面气压升高 (降低 ) ,阿留申低压减弱 (加深 ) ,有利于春季白令海海冰偏少 (多 ) ;而夏季 ,亚洲大陆热低压加深 (减弱 ) ,5 0 0 h Pa西太平洋副热带高压位置偏北 (南 )、强度偏强 (弱 ) ,东亚夏季风易偏强 (弱 )。     

An Evaluation Of Local Similarity At The Top Of The Mixed Layer Based On Large-Eddy Simulations

Local similarity, referred to as type II similarity,in the interfacial, stably-stratified layer at thetop of the atmospheric (or oceanic) mixed layer isdiscussed. Type II scales for scalars are based onthe local values of scalar gradients. Similaritypredictions are derived from the second-orderclosure model of Yamada and Mellor, and also fromsimilarity arguments. The obtainedformulation is verified for active and passive scalarsbased on the large-eddy simulation model.     

An Analysis Of Secondary Circulations And Their Effects Caused By Small-Scale Surface Inhomogeneities Using Large-Eddy Simulation

A parallelized large-eddy simulation model has been used to investigate the effects of two-dimensional, discontinuous, small-scale surface heterogeneities on the turbulence structure of the convective boundary layer.Heterogeneities had a typical size of about the boundary-layer heightz_i. They were produced by a surface sensible heat flux pattern ofchessboard-type and of strong amplitude as typical, e.g., for the marginalice zone. The major objectives of this study were to determinethe effects of such strong amplitude heat flux variations and to specify theinfluence of different speeds and directions of the background wind.Special emphasis has been given to investigate the secondary circulations induced by the heterogeneities by means of three-dimensional phase averages.Compared with earlier studies of continuous inhomogeneities, the same sizeddiscontinuous inhomogeneities in this study show similar but stronger effects.Significant changes compared with uniform surface heating are only observedwhen the scale of the inhomogeneities is increased to z_i. Especially the vertical energy transport is much more vigorous and even the mean emperature profile shows a positive lapse rate within the whole mixed layer. However, the effects are not directly caused by the different shape of the inhomogeneities but can mainly be attributed to the large amplitude of the imposed heat flux,as it is typical for the partially ice covered sea during cold air outbreaks.The structure of the secondary flow is found to be very sensitive to the wavelength and shape of the inhomogeneities as well as to the heatflux amplitude, wind speed and wind direction. The main controlling parameter is the near-surface temperature distribution and the related horizontal pressure gradient perpendicular to the main flow direction. The secondary flow varies from a direct circulation with updraughts mainly above the centre of the heated regions to a more indirect circulation with updraughts beneath the centre and downdraughts above it. For background winds larger than 2.5 m s^–1 a roll-like circulation pattern is observed.From previous findings it has often been stated that moderate backgroundwinds of 5 m s^–1 eliminate all impacts of surface inhomogeneitiesthat could potentially be produced in realistic landscapes. However, this studyshows that the effects caused by increasing the wind speed stronglydepend on the wind direction relative to the orientation of theinhomogeneities. Secondary circulations remain strong, even for abackground wind of 7.5 m s^–1, when the wind direction is orientatedalong one of the two diagonals of the chessboard pattern. On the otherhand, the effects of inhomogeneities are considerably reduced, even undera modest background wind of 2.5 m s^–1, if the wind direction isturned by 45°. Mechanisms for the different flow regimesare discussed.     

Large-Eddy Simulation Of Neutral Turbulent Flow Over Rough Sinusoidal Ridges

The requirements for a credible large-eddy simulation of neutral, turbulent flow over hills with an aerodynamically rough surface are discussed, in order to select a suitable case for simulation. As well as providing adequate resolution within the dynamically important inner region, obtaining a realistic upstream or undisturbed mean velocity profile is also of critical importance. A distributed drag canopy formulation has been introduced to the model to allow it to obtain such a profile while resolving very close to the rough surface. Simulations have then been performed of flow over ridges of varying heights. The results from the steepest case, which is just on the verge of separation, are compared with wind-tunnel observations. It is shown that the large-eddy simulation results are in much better agreement with the experimental data than are the results from a simple first-order mixing-length closure model. An encouraging lack of sensitivity of the simulation results to numerical resolution is also demonstrated.