Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

大地电磁二维对称各向异性介质的有限元数值模拟

假定垂直轴为二维对称各向异性介质主轴之一，构造走向与另一主轴方向成任一夹角，用伽勒金（Ｇａｌｅｒｋｉｎ）有限元法和矩形网格，优化地合成总体刚度矩阵，使二维各向异性介质的基本方程形成有限元代数方程组，求出各节点场值，并利用ＭＯＭ法求出辅助场，进而求出张量阻抗等响应函数．同时还对有关文献的计算模型进行了数值模拟检验     

成都粘土裂隙成因研究

根据成都粘土中裂隙特征与地质环境条件的分析，提出成都粘土中途中眩要裂隙类型的成因模式，上升剥蚀引起的垂向卸荷作用产生波状水平裂隙，冲沟切割引起的侧向卸荷作用，产生平行斜列隙裂。     

地震应力触发数值模拟的初步研究

本文基于有限元方法，初步建立了西南地区二维和三维数值模拟有限元模型，并针对川滇地区，重点建立了包括上地壳、下地壳、上地幔等的三维有限元模型。最后，讨论了地震产生的应力变化对整个西南地区的应力调整情况。结果表明：前一地震对后续地震有一定的触发作用。因此，强震活动有限元模型的建立对判定未来强震危险区有一定的作用。     

Resolution Sensitivity and Scaling of Large-Eddy Simulations of the Stable Boundary Layer

Large-eddy simulations (LES) of the continuously turbulent quasi-equilibrium stable boundary layer (SBL) are conducted with grid lengths in the range of 12.5 m to 2 m, in order to explore resolution sensitivity, and determine at what point grid convergence occurs. The structure of the mean potential temperature, winds, and turbulent fluxes varies significantly over this resolution range. The highest resolution simulations show a significant degree of convergence. The dimensionless momentum diffusivity asymptotes to a value of 0.06, corresponding to a limiting flux Richardson number of 0.15.Using the converged simulations, some scaling hypotheses underpinning first-order and second-order closure models are revisited. The effective Richardson number stability functions of the LES are compared with the forms often used in numerical weather prediction (NWP). The mixing implied by the LES is less than that used in NWP. The commonly used similarity profiles for heat and momentum fluxes, and the scalings for dissipation and pressure covariances are compared with the LES. This information could provide guidance for the next generation of SBL parametrization schemes.     

Large-Eddy Simulation of Turbulent Flow Over a Rough Surface

A family of wall models is proposed that exhibits moresatisfactory performance than previousmodels for the large-eddy simulation (LES) of the turbulentboundary layer over a rough surface.The time and horizontally averaged statistics such asmean vertical profiles of windvelocity, Reynolds stress, turbulent intensities, turbulentkinetic energy and alsospectra are compared with wind-tunnel experimental data.The purpose of the present study is to obtain simulatedturbulent flows that are comparable with wind-tunnelmeasurements for use as the wind environment for thenumerical prediction by LES of source dispersion in theneutral atmospheric boundary layer.     

Laboratory and Numerical Studies of the Convective Boundary Layer Capped by a Strong Inversion

The convective boundary layer (CBL) with a wide range of stability is simulated experimentally using a thermally stratified wind tunnel, and numerically by direct numerical simulation (DNS). The turbulence structures and flow characteristics of various CBL flows, capped by a strong temperature inversion and affected by surface shear, are investigated. The various vertical profiles of turbulence statistics similar to those from the observed CBL in the field are successfully simulated in both the wind-tunnel experiment and in DNS. The comparison of the wind-tunnel data and DNS results with those of atmospheric observations and water-tank studies shows the crucial dependence of the turbulence statistics in the upper part of the layer on the strength of the inversion layer, as well as the modification of the CBL turbulence regime by the surface shear.     

Turbulent Pressure Statistics in the Atmospheric Boundary Layer from Large-Eddy Simulation

We use large-eddy simulation (LES) to study the turbulent pressure field in atmospheric boundary layers with free convection, forced convection, and stable stratification. We use the Poisson equation for pressure to represent the pressure field as the sum of mean-shear, turbulence–turbulence, subfilter-scale, Coriolis, and buoyancy contributions. We isolate these contributions and study them separately. We find that in the energy-containing range in the free-convection case the turbulence–turbulence pressure dominates over the entire boundary layer. That part dominates also up to midlayer in the forced-convection case; above that the mean-shear pressure dominates. In the stable case the mean-shear pressure dominates over the entire boundary layer.We find evidence of an inertial subrange in the pressure spectrum in the free and forced-convection cases; it is dominated by the turbulence–turbulence pressure and has a three-dimensional spectral constant of about 4.0. This agrees well with quasi-Gaussian predictions but is a factor of 2 less than recent results from direct numerical simulations at moderate Reynolds numbers. Measurements of the inertial subrange pressure spectral constant at high Reynolds numbers, which might now be possible, would be most useful.     

模拟试验法确定桩基承载力

通过室内中型剪切摩擦试验确定桩与地层间的C、Φ值,依据莫尔—库仑准则确定桩侧摩阻力。桩侧法向应力采用桩土上部荷载引起的作用于桩身的侧向应力,桩端承载力通过室内三轴试验确定,两者之和即为单桩承载力。该方法用于肇源松花江大桥桩基承载力的测试,得出的桩承载力与实测值吻合较好,取得了预期的效果。     

鄂尔多斯盆地奥陶系海相碳酸盐岩生烃性能与中部长庆气田气源成因研究

用激光 -荧光显微镜等有机岩石学方法 ,挑选样品中有明显生、排烃现象的碳酸盐岩 ,再以测定的这类碳酸盐岩的有机碳含量为依据 ,确定本区高成熟碳酸盐气源岩残余有机碳的下限指标用 0 .13%～ 0 .14 %比较合理。根据碳酸盐岩储层孔隙中充填沥青和皮膜状残余沥青的分布和含量 ,论证了陕参 1井等奥陶系碳酸盐岩地层在地质历史中 ,存在若干油气有效储层。采用连续升温 ,步进取样的热模拟实验方法 ,揭示了本区煤与碳酸盐岩干酪根随热演化作用增强 ,烃类气体产物的干燥系数增大 ,δ1 3C1 、δ1 3C2 明显变重的规律 ,并根据样品实验结果推算 ,本区石炭系—二叠系煤最大生气阶段的δ1 3C1 应为 - 2 7‰～ - 2 8‰ ,碳酸盐岩最大生气阶段的δ1 3C1 应为- 31.7‰～ - 33‰。在气源成因判识中 ,根据主生气期阶段对天然气成藏储聚的贡献最大、甲烷碳同位素较重以及长庆气田风化壳气藏天然气δ1 3C1 重值区并不与石炭系—二叠系煤的镜质组反射率高值区匹配的现象等进行气源判识 ,认为长庆奥陶系风化壳气藏 ,具有下古生界海相碳酸盐岩气与上古生界石炭系—二叠系煤成气的混源成因